Modeling entangled dynamics: comparison between stochastic single-chain and multichain models

To test the effectiveness of stochastic single-chain models in describing the dynamics of entangled polymers, we systematically compare one such model; the slip-spring model; to a multichain model solved using stochastic molecular dynamics(MD) simulations (the Kremer-Grest model). The comparison involves investigating if the single-chain model can adequately describe both a microscopic dynamical and a macroscopic rheological quantity for a range of chain lengths. Choosing a particular chain length in the slip-spring model, the parameter values that best reproduce the mean-square displacement of a group of monomers is determined by fitting toMDdata. Using the same set of parameters we then test if the predictions of the mean-square displacements for other chain lengths agree with the MD calculations. We followed this by a comparison of the time dependent stress relaxation moduli obtained from the two models for a range of chain lengths. After identifying a limitation of the original slip-spring model in describing the static structure of the polymer chain as seen in MD, we remedy this by introducing a pairwise repulsive potential between the monomers in the chains. Poor agreement of the mean-square monomer displacements at short times can be rectified by the use of generalized Langevin equations for the dynamics and resulted in significantly improved agreement.

Temperature dependence of protein dynamics simulated with three different water models

The effect of variation of the water model on the temperature dependence of protein and hydration water dynamics is examined by performing molecular dynamics simulations of myoglobin with the TIP3P, TIP4P, and TIP5P water models and the CHARMM protein force field at temperatures between 20 and 300 K. The atomic mean-square displacements, solvent reorientational relaxation times, pair angular correlations between surface water molecules, and time-averaged structures of the protein are all found to be similar, and the protein dynamical transition is described almost indistinguishably for the three water potentials. The results provide evidence that for some purposes changing the water model in protein simulations without a loss of accuracy may be possible.

The molecular structure of propylene sulphide (methylthiirane) by gas-phase electron diffraction and theoretical calculations: a molecule used in MOCVD

Gas-phase electron-diffraction (GED) data together with results from ab initio molecular orbital calculations have been used to determine the structure of propylene sulphide. Values found for the main structural parameters for the molecule are consistent with those obtained from microwave studies and are compared here with those found for similar sulphur containing rings of general formula S(CH2)n (n = 2–5). A high ring strain enthalpy was calculated for propylene sulphide which is consistent with the small C–S–C angle (48.2(6)degrees) and the relatively long C–S bond lengths (ra = 1.831(2) Å). This is thought to account for the ease of ring opening in propylene sulphide observed in MOCVD reactions and the ready polymerisation of the molecule.

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Segmental dynamics in entangled linear polymer melts

We present molecular dynamics (MD) and slip-springs model simulations of the chain segmental dynamics in entangled linear polymer melts. The time-dependent behavior of the segmental orientation autocorrelation functions and mean-square segmental displacements are analyzed for both flexible and semiflexible chains, with particular attention paid to the scaling relations among these dynamic quantities. Effective combination of the two simulation methods at different coarse-graining levels allows us to explore the chain dynamics for chain lengths ranging from Z ≈ 2 to 90 entanglements. For a given chain length of Z ≈ 15, the time scales accessed span for more than 10 decades, covering all of the interesting relaxation regimes. The obtained time dependence of the monomer mean square displacements, g1(t), is in good agreement with the tube theory predictions. Results on the first- and second-order segmental orientation autocorrelation functions, C1(t) and C2(t), demonstrate a clear power law relationship of C2(t) C1(t)m with m = 3, 2, and 1 in the initial, free Rouse, and entangled (constrained Rouse) regimes, respectively. The return-to-origin hypothesis, which leads to inverse proportionality between the segmental orientation autocorrelation functions and g1(t) in the entangled regime, is convincingly verified by the simulation result of C1(t) g1(t)−1 t–1/4 in the constrained Rouse regime, where for well-entangled chains both C1(t) and g1(t) are rather insensitive to the constraint release effects. However, the second-order correlation function, C2(t), shows much stronger sensitivity to the constraint release effects and experiences a protracted crossover from the free Rouse to entangled regime. This crossover region extends for at least one decade in time longer than that of C1(t). The predicted time scaling behavior of C2(t) t–1/4 is observed in slip-springs simulations only at chain length of 90 entanglements, whereas shorter chains show higher scaling exponents. The reported simulation work can be applied to understand the observations of the NMR experiments.

GSHMC: an efficient method for molecular simulation

The hybrid Monte Carlo (HMC) method is a popular and rigorous method for sampling from a canonical ensemble. The HMC method is based on classical molecular dynamics simulations combined with a Metropolis acceptance criterion and a momentum resampling step. While the HMC method completely resamples the momentum after each Monte Carlo step, the generalized hybrid Monte Carlo (GHMC) method can be implemented with a partial momentum refreshment step. This property seems desirable for keeping some of the dynamic information throughout the sampling process similar to stochastic Langevin and Brownian dynamics simulations. It is, however, ultimate to the success of the GHMC method that the rejection rate in the molecular dynamics part is kept at a minimum. Otherwise an undesirable Zitterbewegung in the Monte Carlo samples is observed. In this paper, we describe a method to achieve very low rejection rates by using a modified energy, which is preserved to high-order along molecular dynamics trajectories. The modified energy is based on backward error results for symplectic time-stepping methods. The proposed generalized shadow hybrid Monte Carlo (GSHMC) method is applicable to NVT as well as NPT ensemble simulations.

Molecular insights into aggregates made of amphiphilic Fmoc-tetrapeptides

Pre-assembled aggregates made of Fmoc-conjugated RGDS and GRDS peptides, where Fmoc refers to fluorenylmethoxycarbonyl, have been investigated using atomistic molecular dynamics simulations. The structural characteristics of twelve different models involving two sheets packed with the Fmoc-aligned or with the charged side groups oriented face-to-face, each one containing seven explicit peptide molecules arranged in parallel or antiparallel, have been evaluated for each Fmoc-tetrapeptide. Structural criteria have been used to select the preferred assembly for each Fmoc-tetrapeptide. The two peptides have been found to prefer b-sheet assemblies with a parallel configuration under simulated low concentration conditions. Furthermore, the assembly is dominated by the interactions among Fmoc units. The overall results provide a complete atomistic view of the interactions between Fmoc-peptide molecules comprised within the same sheet or in different sheets that was not achieved experimentally.

Repulsion-induced surface-migration, by ballistics and bounce

The motion of adsorbate molecules across surfaces is fundamental to self-assembly, material growth, and heterogeneous catalysis. Recent Scanning Tunneling Microscopy studies have demonstrated the electron-induced long-range surface-migration of ethylene, benzene, and related molecules, moving tens of Angstroms across Si(100). We present a model of the previously unexplained long-range recoil of chemisorbed ethylene across the surface of silicon. The molecular dynamics reveal two key elements for directed long-range migration: first ‘ballistic’ motion that causes the molecule to leave the ab initio slab of the surface traveling 3–8 Å above it out of range of its roughness, and thereafter skipping-stone ‘bounces’ that transport it further to the observed long distances. Using a previously tested Impulsive Two-State model, we predict comparable long-range recoil of atomic chlorine following electron-induced dissociation of chlorophenyl chemisorbed at Cu(110)

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

Classical simulation of deposition of thiophene oligomers on TiO(2)-anatase: Relevance of long-range electrostatic interactions

We performed classical molecular dynamics simulations of the vapor-deposition of alpha-T4 oligomers on the TiO(2)-anatase (101) surface, comparing different sets of charges associated with the atoms of the model. The potential energy surfaces for alpha-T4 and TiO(2) were described by re-parametrizations of the Universal force field with charges given by the charge equilibration (QEq) scheme, or with fixed charges obtained by an ab initio method using the Hirshfeld partition. The two sets of charges lead to completely different results for the interface formation, and for the characteristics of the organic film, with a clearly defined alpha-T4 contact layer in the QEq case, and a more homogeneous molecular distribution when using Hirshfeld charges. The main reason for the discrepancy was found to be the incorrect charge assignment given by QEq to the sulfur and alpha-carbon atoms in thiophenes, and highlight the relevance of long-range interactions in the organization of molecular films. (C) 2009 Elsevier B.V. All rights reserved.

Theoretical investigation of Hf and Zr defects in c-Ge

The introduction of high-permittivity gate dielectric materials into complementary metal oxide semiconductor technology has reopened the interest in Ge as a channel material mainly due to its high hole mobility. Since HfO(2) and ZrO(2) are two of the most promising dielectric candidates, it is important to investigate if Hf and Zr may diffuse into the Ge channel. Therefore, using ab initio density functional theory calculations, we have studied substitutional and interstitial Hf and Zr impurities in c-Ge, looking for neutral defects. We find that (i) substitutional Zr and Hf defects are energetically more favorable than interstitial defects; (ii) under oxygen-rich conditions, neither Zr nor Hf migration towards the channel is likely to occur; (iii) either under Hf- or Zr-rich conditions it is very likely, particularly for Zr, that defects will be incorporated in the channel.

Adsorption Behavior of 5-Fluorouracil on Au(111): An In Situ STM Study

The biological effects of chemical substitution of DNA bases triggered several investigations of their physicochemical properties This paper studies the adsorption behavior of a halogenated uracil, 5-fluorouracil (5FU). at the electrochemical interface of Au(111) and sulfuric acid solution. Upon modulation of the electric field across the interface, four distinct phases could be inferred by means of cyclic voltammetry (CV) At negative potentials relative to the SCE electrode, limited by the threshold of hydrogen evolution, no molecular species could be detected by scanning tunneling microscopy (STM) at the reconstructed Au(111)-(23 x root 3) surface, indicating that any physisorbed molecules are randomly distributed Incursion into more positive potentials increases the surface population but doer not form any two-dimensional (2D) physisorbed ordered structure Instead, we observed metastable structures that are only detectable. on surfaces with high defect density At sufficiently high positive potentials. limited by gold oxidation, the molecules are chemisorbed in a (3 x 2 root 3) ordered structure. with the aromatic ring perpendicular to the surface We report the densest chemisorbed monolayer for pyrimidine-derivative molecules (area per molecule 0 14 +/- 0 04 nm(2)). A comparison of the adsorption behavior of uracil derivatives has been made based on recent results of chemical substitution and solvent effects. We propose that pi-stacking is enhanced when halogens are incorporated in the uracil structure, in a similar fashion to what is observed in then crystal structure

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.

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.

Defeitos pontuais nos compostos intermetálicos ZrNi e Zr/sub 2/Ni estudados por dinâmica molecular

Empregamos a técnica de Dinâmica Molecular para estudar propriedades de defeitos pontuais nos compostos intermetálicos ZrNi e Zr2Ni. Descrevemos as configurações estáveis de defeitos e mecanismos de migração, assim como as energias envolvidas. Os potenciais interatômicos foram derivados do Embedded Atom Model. No intuito de levar em conta a variação de estequiometria causada pela presença de alguns tipos de defeitos em intermetálicos, apresentamos um método numérico que fornece a energia efetiva de formação de defeitos e aplicamos o método ao ZrNi e Zr2Ni. Os resultados mostraram que vacâncias são mais estáveis na sub—rede do Ni, com energia de formação de 0,83 e-0,61 eV em ZrNi e Zr2Ni, respectivamente. Vacâncias de Zr são instáveis em ambos compostos; elas decaem espontaneamente em pares anti—sítio e vacância de Ni. Configurações e energias de formação de intersticiais também foram calculadas e mostraram comportamentos similares. Em ZrNi, a migração de vacâncias ocorre preferencialmente nas direções [025] e [100], com as respectivas energias de migração 0,67 e 0,73 eV, e é um processo essencialmente bidimensional no plano (001). Em Zr,Ni, a migração de vacâncias é unidimensional, ocorrendo na direção [001], com energia de migração de 0,67 eV. Em ambos compostos a presença de defeitos de anti—sítio de Ni diminui a energia de migração da vacância de Ni em até 3 vezes e facilita a movimentação em três dimensões. Mecanismos de anel não são energeticamente eficientes em comparação com saltos diretos. As configurações estáveis de intersticiais em ambos compostos consistem em um átomo de Ni sobre o plano (001) entre dois vizinhos de Zr fora do plano. Intersticiais de Zr são instáveis e tendem a deslocar um átomo de Ni, ocupando seu sítio. Energias de deslocamento foram estudadas através de simulações de irradiação de ambos compostos. Durante o processo de colisão binária, um potencial universal ZBL foi usado para colisões a curta distância. Para distâncias intermediárias usamos um potencial de união arbitrário. Zr mostrou—se mais difícil de ser arrancado de seu sítio do que Ni. Encontramos valores de energia de deslocamento no intervalo de aproximadamente 29 eV até 546 eV. Alguns resultados experimentais são mostrados e apresentam boa concordância com os cálculos.