Biomedical subjects

Pavel Hobza

Publications and source records attributed to Pavel Hobza.

At least 19 recordsLinked to original sources

PubMed2007

Density functional theory augmented with an empirical dispersion term. Interaction energies and geometries of 80 noncovalent complexes compared with ab initio quantum mechanics calculations.

Standard density functional theory (DFT) is augmented with a damped empirical dispersion term. The damping function is optimized on a small, well balanced set of 22 van der Waals (vdW) complexes and verified on a validation set of 58 vdW complexes. Both sets contain biologically relevant molecules such as nucleic acid bases. Results are in remarkable agreement with reference high-level wave function data based on the CCSD(T) method. The geometries obtained by full gradient optimization are in very good agreement with the best available theoretical reference. In terms of the standard deviation and average errors, results including the empirical dispersion term are clearly superior to all pure density functionals investigated-B-LYP, B3-LYP, PBE, TPSS, TPSSh, and BH-LYP-and even surpass the MP2/cc-pVTZ method. The combination of empirical dispersion with the TPSS functional performs remarkably well. The most critical part of the empirical dispersion approach is the damping function. The damping parameters should be optimized for each density functional/basis set combination separately. To keep the method simple, we optimized mainly a single factor, s(R), scaling globally the vdW radii. For good results, a basis set of at least triple-zeta quality is required and diffuse functions are recommended, since the basis set superposition error seriously deteriorates the results. On average, the dispersion contribution to the interaction energy missing in the DFT functionals examined here is about 15 and 100% for the hydrogen-bonded and stacked complexes considered, respectively.

Hydrogen Bonding↗

PubMed2007

Hydrogen-bonded nucleic acid base pairs containing unusual base tautomers: complete basis set calculations at the MP2 and CCSD(T) levels.

The total interaction energies of altogether 15 hydrogen-bonded nucleic acid base pairs containing unusual base tautomers were calculated. The geometry properties of all selected adenine-thymine and guanine-cytosine hydrogen-bonded base pairs enable their incorporation into DNA. Unusual base pairing patterns were compared with Watson-Crick H-bonded structures of the adenine-thymine and guanine-cytosine pairs. The complete basis set (CBS) limit of the MP2 interaction energy and the CCSD(T) correction term, determined as the difference between the CCSD(T) and MP2 interaction energies, was evaluated. Extrapolation to the MP2 CBS limit was done using the aug-cc-pVDZ and aug-cc-pVTZ results, and the CCSD(T) correction term was determined with the 6-31G*(0.25) basis set. Final interaction energies were corrected while taking into account both tautomeric penalization determined at the CBS level and solvation/desolvation free energies. The situation for the adenine-thymine pairs is straightforward, and tautomeric pairs are significantly less stable than the Watson-Crick pair consisting of the canonical forms. In the case of the guanine-cytosine pair, the Watson-Crick structure made by canonical forms is again the most stable. The other two structures are, however, energetically rather similar (by 5 and 6 kcal/mol), which provides a very small but non-negligible chance of detecting these structures in the DNA double helix (1:5000). Due to the fact that DNA bases and base pairs incorporated into DNA are solvated less favorably than in isolated systems, this probability represents the very upper limit. The results clearly show how precisely the canonical building blocks of DNA molecules were chosen and how well their stability is maintained.

Chemical Phenomena↗

PubMed2007

How the shape of the NH2 group depends on the substituent effect and H-bond formation in derivatives of aniline.

The geometry and electronic structure of the amino group in aniline and its derivatives are very sensitive to both intramolecular interactions such as substituent effects and intermolecular ones such as H-bonding. An analysis of experimental geometries retrieved from the CSD base and computational modeling of aniline and its derivatives and their H-bonded complexes by use of B3LYP/6-311+G** and MP2/aug-cc-pVDZ showed that the degree of pyramidalization of the amino group depends on H-bonding, which exists in two forms, (i) NH...B (base) and (ii) N...HB (Brønsted acid), both of which affect the shape of the NH2 group. The effect may be significantly enhanced by a substituent through resonance interaction from electron-attracting substituents. The NH...B interactions lead to a substantial planarization of the group, whereas N...HB interactions do not. The natural bond orbital analysis allowed the authors to show that the changes in occupancy of the "lone pair" orbital and in geometry parameters describing pyramidalization of the group depend on the substituent constants.

Journal Article↗

PubMed2007

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.1 e. Electron-density transfer from electron lone pairs of the electron donor is directed into sigma* antibonding orbitals of X--H bonds of the electron acceptor and leads to elongation of the bond and a red shift of the X--H stretching frequency (standard H-bonding). However, pronounced electron-density transfer from electron lone pairs of the electron donor also leads to reorganization of the electron density in the electron donor, which results in changes in geometry and vibrational frequency. These changes are largest for the C--H bonds of formic acid and formate anion, which do not participate in H-bonding. The resulting blue shift of this stretching frequency is substantial and amounts to almost 35 and 170 cm(-1), respectively.

Journal Article↗

PubMed2007

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). A systematic study of these trimers at such a complete level of electronic structure theory is presented for the first time. We show that prior experimental and theoretical interpretations were incorrect in assuming that the most stable structures of the methylated trimers corresponded to planar systems characterized by cyclic intermolecular hydrogen bonding. We found that stacked structures of two bases with the third base in a T-shape arrangement are the global minima in all of the methylated systems: they are more stable than the cyclic planar structures by about 10 kcal mol(-1). The different behaviors of nonmethylated and methylated trimers is also discussed. The high-level geometries and interaction energies computed for the trimers serve also as a reference for the testing of recently developed density functional theory (DFT) functionals with respect to their ability to correctly describe the balance between the electrostatic and dispersion contributions that bind these trimers together. The recently reported M052X functional with a polarized triple-zeta basis set predicts 11 uracil trimer interaction energies with a root-mean-square error of 2.3 kcal mol(-1) relative to highly correlated ab initio theoretical calculations.

Chemistry, Physical↗

PubMed2007

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. The latter method can be used for extended systems. The solvation energy was determined by using a C-PCM continuum solvent at HF level calculations. Various models were adopted to correlate theoretical energies with experimental unwinding free energies. Unless all energy components (hydrogen-bonding, intra- and interstrand-stacking, and solvation energies) were included and weighted individually, no satisfactory correlation resulted. The most advanced model yielded very close correlation (RMSE=0.32 kcal mol(-1)) fully comparable with the entirely empirical correlation introduced in the original paper. Analysis of the theoretical results shows the importance of inter- and intramolecular stacking energies, and especially the latter term plays a key role in determining DNA-duplex stabilization.

Algorithms↗

PubMed2006

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...H(2)O is investigated. The PES of the two lowest excited states of the H(2)O...HCN complex was calculated using the CASPT2 method. The nonadiabatic coupling of the two states of the (pi-->pi*) and (pi-->sigma*) character is responsible for the excited state proton/hydrogen transfer. Compared to the ground state, the barrier for this process is significantly smaller. An increased number of water molecules in the complex with cyclic hydrogen-bonded network causes a large blue shift of the state of the (pi-->sigma*) character. The question of the dissociation of the complex in its excited state is also addressed.

Computer Simulation↗

PubMed2006

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. The cation...Y (Y equal to N or O) distances for divalent metals are shorter than those for Na(+) and for Zn(2+) are mostly shorter than the Mg(2+)...Y distance. The intermolecular distances between the cation and the base for complexes containing adenine and cytosine are systematically shorter than those for complexes containing guanine and thymine. Only for cytosine the canonical keto/amino tautomer structure with ions represents the global minimum. For guanine, the metalated canonical form is again the most stable, but its stabilization energy is within less than 5% of the stabilization energies of the two other rare tautomers, which indicates that the canonical form and these two rare tautomers could coexist. The canonical structures of adenine and thymine in the presence of ions are considerably less stable (by more than 10%) than the complexes of the rare tautomers. It can be concluded that the interaction of Na(+), Mg(2+), and Zn(2+) cations with cytosine in the gas phase will not induce the change of the canonical form to the rare tautomeric form. In the case of isolated guanine, the equilibrium of the canonical form with rare tautomers can be found. For isolated adenine and thymine the presence of rare tautomers is highly probable.

Adenine↗

PubMed2006

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. The energies of interaction between roscovitine and small fragments and between roscovitine and substantial sections of protein (722 atoms) were also computed by using density-functional tight-binding methods covering dispersion energy (DFTB-D) and the Cornell empirical potential. Total stabilisation energy originates predominantly from dispersion energy and methods that do not account for the dispersion energy cannot, therefore, be recommended for the study of protein-inhibitor interactions. The Cornell empirical potential describes reasonably well the interaction between roscovitine and protein; therefore, this method can be applied in future thermodynamic calculations. A limited number of amino acid residues contribute significantly to the binding of roscovitine and cdk2, whereas a rather large number of amino acids make a negligible contribution.

Amino Acid Sequence↗

PubMed2006

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. The RESP and NPA charges of carboranes and small model systems are compared and their use is discussed. The dominant interaction between carboranes and biomolecules is the formation of unconventional proton-hydride hydrogen bonds (dihydrogen bonds) characterized by a short distance between hydrogen atoms (as close as 1.8 A) and an average strength in the range of 4.2-5.8 kcal mol(-1). The total stabilization energy of complexes investigated is rather large, and the largest value (approximately 15 kcal mol(-1)) was found for the carborane complexes with ARG and the GLY-SER dipeptide. These interactions are ubiquitous under geometrical constraints influencing the strength of the interaction. The carborane forms dihydrogen bonds with biomolecules preferably with the hydrogen atoms of its lower hemisphere (i.e. the part of the cage opposite to the carbon atom). These two geometrical factors can be used to explain the specificity of inhibition of HIV protease by carboranes.

Biophysics↗

PubMed2006

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. This structure, intermediate between enol and keto tautomers, might play a role in the excited state relaxation. The existence of this second well, however, is observed in both 7- and 9-methyl-guanine, which suggests that it cannot account alone for the different photophysical behavior of these species.

Computer Simulation↗

PubMed2006

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. The geometry of the model was averaged from structures determined by MD simulations. The role of the solvent was covered by the COSMO continuum solvent model and calculations were performed for a larger model system (model B) which also contained a sugar-phosphate backbone. The total stabilization energies of the unperturbed system and the system perturbed by a base-analogue pair (model A) were comparable to the stability of both duplexes experimentally determined. This is due to large stacking interaction energy of the base-analogue self-pair which compensates for the missing hydrogen-bonding energy of the replaced adenine...thymine base pair. The selectivity of the base-analogue pair was reproduced (model B) when their desolvation energy was included with the interaction energy of both strands determined by the approximate SCC-DFTB-D method.

Base Pairing↗

PubMed2006

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. Common geometrical features have been systematically observed for the sequence Trp, Trp-Gly, and Trp-Gly-Gly. In addition, the peptide backbone of Trp-Gly-Gly has been compared with that of the previously studied Phe-Gly-Gly (Reha, D. et. al. Chem. Eur. J. 2005, 11, 6803). From the observed systematic structural behavior between these peptide analogues, it is expected that the gas-phase conformers of other similar aromatic small peptides would present equivalent geometries. The DFT methodology failed to describe the potential energy surface of the studied peptides since the London dispersion energy (not covered in DFT) plays a significant role in the stabilization of most stable conformers.

Oligopeptides↗

PubMed2006

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. The new calculations show modest increases in stacking stabilization compared to the MP2/6-31G*(0.25) data and surprisingly large sequence-dependent variation of stacking energies. The absolute force-field values are in better agreement with the new reference data, while relative discrepancies between quantum-chemical (QM) and force-field values increase modestly. Nevertheless, the force field provides good qualitative description of stacking, and there is no need to introduce additional pair-additive electrostatic terms, such as distributed multipoles or out-of-plane charges. There is a rather surprising difference of about 0.1 A between the vertical separation of base pairs predicted by quantum chemistry and derived from crystal structures. Evaluations of different local arrangements of the 5'-CG-3' step indicate a sensitivity of the relative stacking energies to the level of calculation. Thus, describing quantitative relations between local DNA geometrical variations and stacking may be more complicated than usually assumed. The reference calculations are complemented by continuum-solvent assessment of solvent-screening effects.

Base Pairing↗

PubMed2006

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

The benzene...X complexes (X=benzene, antracene, ovalene) were optimised at the MP2/6-31G** level with the C2v symmetry of the complex and planarity of the proton acceptor being preserved. The resulting stabilisation energies amount to 1.2, 2.3 and 2.9 kcal mol(-1), and the C-H bond of the proton donor is contracted by 0.0035, 0.0052 and 0.0055 A, respectively. The contraction is connected with a blue-shift of the C-H stretch vibration frequency. A two-dimensional anharmonic vibration treatment based on a MP2/6-31G** potential energy surface yields the following blue shifts for the complexes studied: 28, 42 and 43 cm(-1). The dominant attraction in the complexes is London dispersion, while the attractive contribution from electrostatic quadrupole-quadrupole interactions is considerably smaller.

Journal Article↗

PubMed2006

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. Alternatively, all contributions were evaluated consistently using empirical methods. All the values of the interaction energy of helix pairs are stabilizing, and the dominant stabilizing terms stem from the London dispersion energy and, in the case of charged systems, the electrostatic energy. The stabilization energy of the core, determined as the difference of the energy of the core and 10 separate helices, amounts to approximately 450 kcal/mol. Systematically, the difference in the hydration free energy of a helix pair and its separate components is smaller in magnitude than the interaction energy, and it is negative for some pairs while positive for others. The average total free energy of a core formation amounts to -29.6 kcal/mol (yielded by scaled quantum-chemical methods) and +13.9 kcal/mol (resulting from empirical methods). These values are considerably smaller than their single components, which are dominated by the interaction energy. The computationally predicted interval encloses the experimental value of the folding free energy (-2.8 kcal/mol).

Amino Acid Motifs↗

PubMed2006

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). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.

Amino Acids↗

PubMed2006

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...cytosine and 9-methylguanine...1-methylcytosine base pairs was realized by the molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The most stable and populated structures localized were further fully reoptimized at the correlated ab initio level employing the resolution of identity Møller-Plesset method with a large basis set. A systematic study of microhydration of these systems using a high-level correlated ab initio approach is presented for the first time. The different behavior of guanine...cytosine and adenine...thymine complexes is also discussed. These studies of nucleic acid base pairs are important for finding binding sites of water molecules around bases and for better understanding of the influence of the solvent on the stability of the structure of DNA.

Base Pairing↗