On the nature of DNA-duplex stability.

The unwinding free energy of 128 DNA octamers was correlated with the sum of interaction energies among DNA bases and their solvation energies. The former energies were determined by using the recently developed density functional theory procedure augmented by London dispersion energy (RI-DFT-D) that provides accurate hydrogen-bonding and stacking energies highly comparable with CCSD(T)/complete basis set limit benchmark data. Efficient tight-binding DFT covering dispersion energy was also used and yielded satisfactory results.

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

Quantification and rationalization of the higher affinity of sodium over potassium to protein surfaces.

For a series of different proteins, including a structural protein, enzyme, inhibitor, protein marker, and a charge-transfer system, we have quantified the higher affinity of Na+ over K+ to the protein surface by means of molecular dynamics simulations and conductivity measurements. Both approaches show that sodium binds at least twice as strongly to the protein surface than potassium does with this effect being present in all proteins under study. Different parts of the protein exterior are responsible to a varying degree for the higher surface affinity of sodium, with the charged carboxylic groups of aspartate and glutamate playing the most important role.

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.

Unified molecular picture of the surfaces of aqueous acid, base, and salt solutions.

The molecular structure of the interfacial regions of aqueous electrolytes is poorly understood, despite its crucial importance in many biological, technological, and atmospheric processes. A long-term controversy pertains between the standard picture of an ion-free surface layer and the strongly ion specific behavior indicating in many cases significant propensities of simple inorganic ions for the interface. Here, we present a unified and consistent view of the structure of the air/solution interface of aqueous electrolytes containing monovalent inorganic ions.

Interior and interfacial aqueous solvation of benzene dicarboxylate dianions and their methylated analogues: A combined molecular dynamics and photoelectron spectroscopy study.

Aqueous solvation of benzene dicarboxylate dianions (BCD(2-)) was studied by means of photoelectron spectroscopy and molecular dynamics simulations. Photoelectron spectra of hydrated o- and p-BCD(2-) with up to 25 water molecules were obtained. An even-odd effect was observed for the p-BCD(2-) system as a result of the alternate solvation of the two negative charges.

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.

Specific ion effects at protein surfaces: a molecular dynamics study of bovine pancreatic trypsin inhibitor and horseradish peroxidase in selected salt solutions.

The distribution of sodium, choline, sulfate, and chloride ions around two proteins, horseradish peroxidase (HRP) and bovine pancreatic trypsin inhibitor (BPTI), is investigated by means of molecular dynamics simulations with the aim to elucidate ion adsorption at the protein surface. Although the two proteins under investigation are very different from each other, the ion distributions around them are remarkably similar. Sulfate is always strongly attached to the proteins, choline shows a significant, but unspecific, propensity for the protein surfaces, and sodium ions have a weak surface affinity, while chloride has virtually no preference for the protein surface.

Aqueous ionic and complementary zwitterionic soluble surfactants: molecular dynamics simulations and sum frequency generation spectroscopy of the surfaces.

An aqueous ionic surfactant, 1-dodecyl-4-(dimethylamino)pyridinium (DMP) bromide, and the corresponding zwitterion 2-[4-(dimethylamino)pyridinio]dodecanoate (DPN) were explored by means of molecular dynamics (MD) simulations and, for the ionic system, by infrared-visible sum frequency generation (IR-vis SFG). The molecular structure of the interfacial layer was investigated for the ionic and zwitterionic systems as a function of surfactant concentration, both in water and in salt (KF or KBr) solutions, by MD simulations in a slab geometry. The buildup of the surface monolayer and a sublayer was monitored, and density and orientational profiles of the surfactants were evaluated.

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.

Different mechanisms of CDK5 and CDK2 activation as revealed by CDK5/p25 and CDK2/cyclin A dynamics.

A detailed analysis is presented of the dynamics of human CDK5 in complexes with the protein activator p25 and the purine-like inhibitor roscovitine. These and other findings related to the activation of CDK5 are critically reviewed from a molecular perspective. In addition, the results obtained on the behavior of CDK5 are compared with data on CDK2 to assess the differences and similarities between the two kinases in terms of (i) roscovitine binding, (ii) regulatory subunit association, (iii) conformational changes in the T-loop following CDK/regulatory subunit complex formation, and (iv) specificity in CDK/regulatory subunit recognition.

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.

Brine rejection from freezing salt solutions: a molecular dynamics study.

The atmospherically and technologically very important process of brine rejection from freezing salt solutions is investigated with atomic resolution using molecular dynamics simulations. The present calculations allow us to follow the motion of each water molecule and salt ion and to propose a microscopic mechanism of brine rejection, in which a fluctuation (reduction) of the ion density in the vicinity of the ice front is followed by the growth of a new ice layer. The presence of salt slows down the freezing process, which leads to the formation of an almost neat ice next to a disordered brine layer.

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.