scholarly journals Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments

Author(s):  
Ondrej Gutten ◽  
Petr Jurečka ◽  
Zahra Aliakbar Tehrani ◽  
Miloš Buděšínský ◽  
Jan Řezáč ◽  
...  

Computational “error bars” for modelling cyclic dinucleotides – NMR experiment vs. quantum mechanics and molecular dynamics.

Author(s):  
Walker M. Jones ◽  
Aaron G. Davis ◽  
R. Hunter Wilson ◽  
Katherine L. Elliott ◽  
Isaiah Sumner

We present classical molecular dynamics (MD), Born-Oppenheimer molecular dynamics (BOMD), and hybrid quantum mechanics/molecular mechanics (QM/MM) data. MD was performed using the GPU accelerated pmemd module of the AMBER14MD package. BOMD was performed using CP2K version 2.6. The reaction rates in BOMD were accelerated using the Metadynamics method. QM/MM was performed using ONIOM in the Gaussian09 suite of programs. Relevant input files for BOMD and QM/MM are available.


2018 ◽  
Vol 2018 ◽  
pp. 1-14 ◽  
Author(s):  
Eric Paquet ◽  
Herna L. Viktor

Ab initio molecular dynamics is an irreplaceable technique for the realistic simulation of complex molecular systems and processes from first principles. This paper proposes a comprehensive and self-contained review of ab initio molecular dynamics from a computational perspective and from first principles. Quantum mechanics is presented from a molecular dynamics perspective. Various approximations and formulations are proposed, including the Ehrenfest, Born–Oppenheimer, and Hartree–Fock molecular dynamics. Subsequently, the Kohn–Sham formulation of molecular dynamics is introduced as well as the afferent concept of density functional. As a result, Car–Parrinello molecular dynamics is discussed, together with its extension to isothermal and isobaric processes. Car–Parrinello molecular dynamics is then reformulated in terms of path integrals. Finally, some implementation issues are analysed, namely, the pseudopotential, the orbital functional basis, and hybrid molecular dynamics.


2011 ◽  
Vol 10 (03) ◽  
pp. 359-370 ◽  
Author(s):  
JUAN PANG ◽  
HU YANG ◽  
JING MA ◽  
RONGSHI CHENG

Poly(N-alkylacrylamide) is a group of thermo-sensitive polymers that include poly (N-isopropylacrylamide), poly(N-n-propylacrylamide), poly(N-isopropylmethacryl-amide), and so on. The polymers exhibit different levels of lower critical solution temperatures (LCST) in aqueous solutions. In this article, their monomers and oligomers with 10 repeating units are selected, respectively, to demonstrate the cause of different LCST levels of the polymers in aqueous solutions using molecular dynamics simulations and quantum mechanics calculations. The monomers have functional groups of different steric volume that greatly affect the conformational transition of chains and LCST levels of the polymers. A branched chain of N-propyl group in N-isopropylacrylamide and an additional methyl group at α-carbon in N-isopropylmethacrylamide both increase the steric effect, making it more difficult for monomers to draw closer and resulting in higher LCST levels of the polymers. In addition, the simulated results from their corresponding oligomers exhibit the similar trend to those from the monomers.


2013 ◽  
Vol 9 ◽  
pp. 118-134 ◽  
Author(s):  
Jutta Erika Helga Köhler ◽  
Nicole Grczelschak-Mick

Four highly ordered hydrogen-bonded models of β-cyclodextrin (β-CD) and its inclusion complex with benzene were investigated by three different theoretical methods: classical quantum mechanics (QM) on AM1 and on the BP/TZVP-DISP3 level of approximation, and thirdly by classical molecular dynamics simulations (MD) at different temperatures (120 K and 273 to 300 K). The hydrogen bonds at the larger O2/O3 rim of empty β-CDs prefer the right-hand orientation, e.g., O3-H…O2-H in the same glucose unit and bifurcated towards …O4 and O3 of the next glucose unit on the right side. On AM1 level the complex energy was −2.75 kcal mol−1 when the benzene molecule was located parallel inside the β-CD cavity and −2.46 kcal mol−1 when it was positioned vertically. The AM1 HOMO/LUMO gap of the empty β-CD with about 12 eV is lowered to about 10 eV in the complex, in agreement with data from the literature. AM1 IR spectra displayed a splitting of the O–H frequencies of cyclodextrin upon complex formation. At the BP/TZVP-DISP3 level the parallel and vertical positions from the starting structures converged to a structure where benzene assumes a more oblique position (−20.16 kcal mol−1 and −20.22 kcal mol−1, resp.) as was reported in the literature. The character of the COSMO-RS σ-surface of β-CD was much more hydrophobic on its O6 rim than on its O2/O3 side when all hydrogen bonds were arranged in a concerted mode. This static QM picture of the β-CD/benzene complex at 0 K was extended by MD simulations. At 120 K benzene was mobile but always stayed inside the cavity of β-CD. The trajectories at 273, 280, 290 and 300 K certainly no longer displayed the highly ordered hydrogen bonds of β-CD and benzene occupied many different positions inside the cavity, before it left the β-CD finally at its O2/O3 side.


2012 ◽  
Vol 34 (9) ◽  
pp. 750-756 ◽  
Author(s):  
Marcin Nowosielski ◽  
Marcin Hoffmann ◽  
Aneta Kuron ◽  
Malgorzata Korycka-Machala ◽  
Jaroslaw Dziadek

Sign in / Sign up

Export Citation Format

Share Document