DNA base trimers: empirical and quantum chemical ab initio calculations versus experiment in vacuo.

A complete scan of the potential-energy surfaces for selected DNA base trimers has been performed by a molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The resulting most stable/populated structures were then reoptimized at a correlated ab initio level by employing resolution of the identity, Møller-Plesset second-order perturbation theory (RI-MP2).

On the origin of red and blue shifts of X-H and C-H stretching vibrations in formic acid (formate ion) and proton donor complexes.

Complexes between formic acid or formate anion and various proton donors (HF, H(2)O, NH(3), and CH(4)) are studied by the MP2 and B3LYP methods with the 6-311++G(3df,3pd) basis set. Formation of a complex is characterized by electron-density transfer from electron donor to ligands. This transfer is much larger with the formate anion, for which it exceeds 0.

Theoretical study of photoacidity of HCN: the effect of complexation with water.

The character of the hydrogen bonding and the excited state proton transfer (ESPT) in the model system HCN...

Na+, Mg2+, and Zn2+ binding to all tautomers of adenine, cytosine, and thymine and the eight most stable keto/enol tautomers of guanine: a correlated ab initio quantum chemical study.

Interactions of adenine, cytosine, guanine, and thymine with Na(+), Mg(2+), and Zn(2+) cations were studied using an approximate resolution of identity correlated second-order MP2 (RI-MP2) method with the TZVPP ([5s3p2d1f/3s2p1d]) basis set. All existing tautomers of adenine, cytosine, and thymine and the eight most stable keto/enol tautomers of guanine were considered. Cations bind mostly in a bidentate manner, and stabilization energies of these complexes are larger than those in the case when cations bind in a unidentate manner.

Hydration-induced changes of structure and vibrational frequencies of methylphosphocholine studied as a model of biomembrane lipids.

The chemical characteristics of the polar parts of phospholipids as the main components of biological membranes were investigated by using infrared (IR) spectroscopy and theoretical calculations with water as a probe molecule. The logical key molecule used in this study is methylphosphocholine (MePC) as it is not only a representative model for a polar lipid headgroup but itself has biological significance. Isolated MePC forms a compact (folded) structure which is essentially stabilized by two intramolecular C-H.

Potential energy surfaces of an adenine-thymine base pair and its methylated analogue in the presence of one and two water molecules: molecular mechanics and correlated ab initio study.

Potential energy surfaces of monohydrated and dihydrated adenine-thymine and 9-methyladenine-1-methylthymine base pairs were examined by the molecular dynamics/quenching technique using the Cornell et al. force field (J. Am.

Vibrational spectroscopy of the G...C base pair: experiment, harmonic and anharmonic calculations, and the nature of the anharmonic couplings.

The results of harmonic and anharmonic frequency calculations on a guanine-cytosine complex with an enolic structure (a tautomeric form with cytosine in the enol form and with a hydrogen at the 7-position on guanine) are presented and compared to gas-phase IR-UV double resonance spectral data. Harmonic frequencies were obtained at the RI-MP2/cc-pVDZ, RI-MP2/TZVPP, and semiempirical PM3 levels of electronic structure theory. Anharmonic frequencies were obtained by the CC-VSCF method with improved PM3 potential surfaces; the improved PM3 potential surfaces are obtained from standard PM3 theory by coordinate scaling such that the improved PM3 harmonic frequencies are the same as those computed at the RI-MP2/cc-pVDZ level.

Stabilization energies of the hydrogen-bonded and stacked structures of nucleic acid base pairs in the crystal geometries of CG, AT, and AC DNA steps and in the NMR geometry of the 5'-d(GCGAAGC)-3' hairpin: Complete basis set calculations at the MP2 and CCSD(T) levels.

Stabilization energies of the H-bonded and stacked structures of a DNA base pair were studied in the crystal structures of adenine-thymine, cytosine-guanine, and adenine-cytosine steps as well as in the 5'-d(GCGAAGC)-3' hairpin (utilizing the NMR geometry). Stabilization energies were determined as the sum of the complete basis set (CBS) limit of MP2 stabilization energies and the Delta E(CCSD(T)) - Delta E(MP2) correction term evaluated with the 6-31G*(0.25) basis set.

Theoretical study of the ground and excited states of 7-methyl guanine and 9-methyl guanine: comparison with experiment.

The keto-enol tautomerization of 7-methyl-guanine and 9-methyl-guanine in the excited state was investigated using the time-dependent DFT (TDDFT) method. For both species, the potential energy surfaces of the ground state and two lowest singlet excited states (due to pi-->pi* and n-->pi* transitions) have been investigated and their features discussed in terms of consequences on the excited state dynamics. The findings suggest that, for both species, the state due to the n-->pi* transition, suspected to be an intermediate in the excited state deactivation, exhibits two minima with the second minimum characterized by an elongated N1-H distance.

Interaction of carboranes with biomolecules: formation of dihydrogen bonds.

Noncovalent interactions of the polyhedral carborane 1-carba-closo-dodecaborane (CB(11)H(12))(-) with building blocks of biomolecules, modelled by glycine (GLY), serine (SER), phenylalanine (PHE), glutamic acid (GLU), lysine (LYS) and arginine (ARG), were investigated in vacuo by molecular dynamics simulations with the UFF empirical potential. Selected structures were further studied by accurate ab initio quantum chemical procedures. Interactions with a peptide bond (GLY-SER dipeptide) and a nucleic acid building block (guanine) were also considered.

Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs.

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ).

Interaction energies for the purine inhibitor roscovitine with cyclin-dependent kinase 2: correlated ab initio quantum-chemical, DFT and empirical calculations.

The interaction between roscovitine and cyclin-dependent kinase 2 (cdk2) was investigated by performing correlated ab initio quantum-chemical calculations. The whole protein was fragmented into smaller systems consisting of one or a few amino acids, and the interaction energies of these fragments with roscovitine were determined by using the MP2 method with the extended aug-cc-pVDZ basis set. For selected complexes, the complete basis set limit MP2 interaction energies, as well as the coupled-cluster corrections with inclusion of single, double and noninteractive triples contributions [CCSD(T)], were also evaluated.

Structure of isolated tryptophyl-glycine dipeptide and tryptophyl-glycyl-glycine tripeptide: ab initio SCC-DFTB-D molecular dynamics simulations and high-level correlated ab initio quantum chemical calculations.

The tryptophyl-glycine (Trp-Gly) and tryptophyl-glycyl-glycine (Trp-Gly-Gly) peptides have been studied by means of molecular dynamic simulations combined with high-level correlated ab initio quantum chemical and statistical thermodynamic calculations. The lowest energy conformers were localized in the free energy surface. The structures of the different Trp-Gly and Trp-Gly-Gly conformers coexisting in the gas phase have been for the first time reported and their scaled theoretical IR spectra unambiguously assigned and compared with previous gas-phase experimental results.

Theoretical investigations into the blue-shifting hydrogen bond in benzene complexes.

The benzene...

How the stabilization of INK4 tumor suppressor 3D structure evaluated by quantum chemical and molecular mechanics calculations corresponds well with experimental results: interplay of association enthalpy, entropy, and solvation effects.

The folding free energy of the INK4c tumor suppressor core, consisting of 10 helices, was determined as the sum of gas-phase interaction enthalpy, gas-phase interaction entropy, and dehydration and hydration free energy. The interaction energy and the hydration free energy were determined using the nonempirical density functional theory (DFT) method, augmented by a dispersion-energy correction term, the semiempirical density-functional tight-binding method covering the dispersion energy, and the density functional theory/conductor-like screening model (DFT/COSMO) procedure, whereas the interaction entropy was calculated with the empirical Cornell et al. force field.

Exceptional thermodynamic stability of DNA duplexes modified by nonpolar base analogues is due to increased stacking interactions and favorable solvation: Correlated ab initio calculations and molecular dynamics simulations.

The geometries of DNA hexamer (5'-GGAACC-3') and DNA 13-mer (5'-GCGTACACATGCG-3') have been determined by molecular dynamics (MD) simulations using an empirical force field. The central canonical base pair was replaced by a pair of nonpolar base analogues, 2,2'-bipyridyl and 3-methylisocarbostyril. The stabilization energy of the model system (model A) consisting of a central base pair (base-analogue pair) and two neighboring base pairs was determined by the RI-MP2 method using an extended aug-cc-pVDZ basis set.

Potential energy surfaces of the microhydrated guanine...cytosine base pair and its methylated analogue.

A complete scan of the potential and free-energy surfaces of monohydrated and dihydrated guanine...

Nature of base stacking: reference quantum-chemical stacking energies in ten unique B-DNA base-pair steps.

Base-stacking energies in ten unique B-DNA base-pair steps and some other arrangements were evaluated by the second-order Møller-Plesset (MP2) method, complete basis set (CBS) extrapolation, and correction for triple (T) electron-correlation contributions. The CBS(T) calculations were compared with decade-old MP2/6-31G*(0.25) reference data and AMBER force field.

Binding of cationic and neutral phenanthridine intercalators to a DNA oligomer is controlled by dispersion energy: quantum chemical calculations and molecular mechanics simulations.

Correlated ab initio as well as semiempirical quantum chemical calculations and molecular dynamics simulations were used to study the intercalation of cationic ethidium, cationic 5-ethyl-6-phenylphenanthridinium and uncharged 3,8-diamino-6-phenylphenanthridine to DNA. The stabilization energy of the cationic intercalators is considerably larger than that of the uncharged one. The dominant energy contribution with all intercalators is represented by dispersion energy.

IR-UV double resonance spectroscopy of guanine-H2O clusters.

We present the IR-UV double resonance spectrum of guanine monohydrate in the region 3100 cm(-1) to 3800 cm(-1) along with the energies and frequencies of these structures calculated at the non-empirical correlated ab initio RI-MP2/cc-pVDZ level. We assign the structures of guanine-water clusters by comparing the experimental spectra with the ab initio calculations and with the IR spectra of the bare guanine monomer. We find two clusters with guanine in the enol-amino tautomeric form and one structure with guanine in the keto-amino form.

Structure and IR spectrum of phenylalanyl-glycyl-glycine tripetide in the gas-phase: IR/UV experiments, ab initio quantum chemical calculations, and molecular dynamic simulations.

We investigated the potential-energy surface (PES) of the phenylalanyl-glycyl-glycine tripeptide in the gas phase by means of IR/UV double-resonance spectroscopy, and quantum chemical and statistical thermodynamic calculations. Experimentally, we observed four conformational structures and we recorded their IR spectra in the spectral region of 3000-4000 cm(-1). Computationally, we investigated the PES by a combination of molecular dynamics/quenching procedures with high-level correlated ab initio calculations.

Are the hydrogen bonds of RNA (AU) stronger than those of DNA (AT)? A quantum mechanics study.

The intrinsic stability of Watson-Crick d(AT) and r(AU) hydrogen bonds was analyzed by employing a variety of quantum-mechanical techniques, such as energy calculations, determination of reactivity indexes, and analysis of electron density topology. The analyses were performed not only for equilibrium gas-phase geometries, but also on hundreds of conformations derived from molecular dynamics (MD) and database analysis. None of our results support the idea that r(AU) hydrogen bonds are intrinsically more stable than those of d(AT).

On geometries of stacked and H-bonded nucleic acid base pairs determined at various DFT, MP2, and CCSD(T) levels up to the CCSD(T)/complete basis set limit level.

The geometries and interaction energies of stacked and hydrogen-bonded uracil dimers and a stacked adeninecdots, three dots, centeredthymine pair were studied by means of high-level quantum chemical calculations. Specifically, standard as well as counterpoise-corrected optimizations were performed at second-order Moller-Plesset (MP2) and coupled cluster level of theory with single, double, and perturbative triple excitations [CCSD(T)] levels with various basis sets up to the complete basis set limit. The results can be summarized as follows: (i) standard geometry optimization with small basis set (e.

On differences between hydrogen bonding and improper blue-shifting hydrogen bonding.

Twenty two hydrogen-bonded and improper blue-shifting hydrogen-bonded complexes were studied by means of the HF, MP2 and B3LYP methods using the 6-31G(d,p) and 6--311 ++G(d,p) basis sets. In contrast to the standard H bonding, the origin of the improper blue-shifting H bonding is still not fully understood. Contrary to a frequently presented idea, the electric field of the proton acceptor cannot solely explain the different behavior of the H-bonded and improper blue-shifting H-bonded complexes.

Correlated ab initio study of nucleic acid bases and their tautomers in the gas phase, in a microhydrated environment and in aqueous solution. Part 4. Uracil and thymine.

Altogether 13 keto and enol tautomers of uracil and 13 keto and enol tautomers of thymine were studied theoretically in the gas phase, in a microhydrated environment (1 and 2 water molecules) and in a water environment. Bulk water was described using the thermodynamic integration method, Conductor-like polarizable continuum model (C-PCM, COSMO) and hybrid model (C-PCM + 1-2 explicit water molecules). The structures of various tautomers were determined at the RI-MP2 level using the TZVPP basis set while relative energies were determined at the CCSD(T) level.