Accurate description of van der Waals complexes by density functional theory including empirical corrections

An empirical method to account for van der Waals interactions in practical calculations with the density functional theory (termed DFT-D) is tested for a wide variety of molecular complexes. As in previous schemes, the dispersive energy is described by damped interatomic potentials of the form C6R(-6). The use of pure, gradient-corrected density functionals (BLYP and PBE), together with the resolution-of-the-identity (RI) approximation for the Coulomb operator, allows very efficient computations for large systems. Opposed to previous work, extended AO basis sets of polarized TZV or QZV quality are employed, which reduces the basis set superposition error to a negligible extend. By using a global scaling factor for the atomic C6 coefficients, the functional dependence of the results could be strongly reduced. The "double counting" of correlation effects for strongly bound complexes is found to be insignificant if steep damping functions are employed. The method is applied to a total of 29 complexes of atoms and small molecules (Ne, CH4, NH3, H2O, CH3F, N2, F2, formic acid, ethene, and ethine) with each other and with benzene, to benzene, naphthalene, pyrene, and coronene dimers, the naphthalene trimer, coronene. H2O and four H-bonded and stacked DNA base pairs (AT and GC). In almost all cases, very good agreement with reliable theoretical or experimental results for binding energies and intermolecular distances is obtained. For stacked aromatic systems and the important base pairs, the DFT-D-BLYP model seems to be even superior to standard MP2 treatments that systematically overbind. The good results obtained suggest the approach as a practical tool to describe the properties of many important van der Waals systems in chemistry. Furthermore, the DFT-D data may either be used to calibrate much simpler (e.g., force-field) potentials or the optimized structures can be used as input for more accurate ab initio calculations of the interaction energies.     

Conformational analysis of [Met5]-enkephalin: Solvation and ionization considerations

[Met5]-Enkephalin has the sequence Tyr-Gly-Gly-Phe-Met. Only the extended conformation of the peptide has been observed by X-ray crystallography. Nuclear magnetic resonance spectroscopy supports the presence of a turn at Gly 3 and Phe 4 in dimethyl sulfoxide. In this study, the peptide conformational states and thermodynamic properties are understood in terms of ionization state and solvent environment. In the calculation, final conformations obtained from multiple independent Monte Carlo simulated annealing conformational searches are starting points for molecular dynamics simulations. In an aqueous environment given by the use of solvation free energy and the zwitterionic state, dominant structural motifs computed are G-P Type II bend, G-G Type II bend, and G-G Type I bend motifs, in order of increasing free energy. In the calculation of the peptide with neutral N- and C-termini and solvation free energy, the extended conformer dominates (by at least a factor of 2.5), and the conformation of another low free energy conformer superimposes well on the pharmacophoric groups of morphine. Neutralization of charge and solvation induce and stabilize the extended conformation, respectively. A mechanism of inter-conversion between the extended conformer and three bent conformers is supported by /-scatter plots, and by the conformer relative free energies. An estimate of the entropy change of receptor unbinding is 8.3 cal K-1 mol-1, which gives rise to a -2.5 kcal/mol entropy contribution to the free energy of unbinding at 25 °C. The conformational analysis methodology described here should be useful in studies on short peptides and flexible protein surface loops that have important biological implications.     

Ionization of five aqueous fluorobenzoic acids: Conductance and thermodynamics

The three monofluorobenzoic acids together with 2,4-difluoro and 2,6-difluorobenzoic acids in aqueous solution are the subject of precision conductance measurements. The experimental data are analyzed to give ionization constants and limiting conductances at temperatures from 0 to 100°C. Walden products for the acid anions are derived from the limiting conductances while the ionization consatants are fitted by statistical methods to the function pK _a (m)=A+B/T+ C logT+DT. Only the 2,6- acid requires the fourth term of the function to fit the data to a precision of better than 0.03%. Mathematical analysis of the pK function gives the standard changes in enthalpy, entropy, and heat capacity. All the acids studied are more acidic than the parent, benzoic acid, as well as more acidic than the isoelectronic methylbenzoic acids. In general the increased acidity is tied to decreases in enthalpy while entropy changes on ionization differn little from those found for the parent acid.     

Net information measures for modified Yukawa and Hulthén potentials

The dimensional analyses of the position and momentum variances‐based quantum mechanical Heisenberg uncertainty measure, as well as the entropic information measures given by the Shannon information entropy sum and the product of Fisher information measures are carried out for two widely used nonrelativistic isotropic exponential‐cosine screened Coulomb potentials generated by multiplying the superpositions of (i) Yukawa‐like, ?Z(e^?μr/r), and (ii) Hulthén‐like, ?Zμ(1/(e^μr ? 1)), potentials by cos(bμr) followed by addition of the term a/r², where a and b ≥ 0, μ are the screening parameters and Z, in case of atoms, denotes the nuclear charge. Under the spherical symmetry, all the information measures considered are shown to be independent of the scaling of the set [μ, Z] at a fixed value of μ/Z, a, and b and the other parameters defining the superpositions of the potentials. Numerical results are presented, which support the validity of the scaling properties. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Quantitative Correlations for the Estimation of Thermodynamic and Molecular Properties of Pyridines and 2,2′-Bipyridyines

Standard heats of formation, entropies, ionization potentials, and molecular dipole moments of a series of pyridines have been calculated by MNDO, AM1 and PM3 methods. Linear relationship have been established which permit a priori estimation of thermodynamic and molecular characteristics of pyridines. Correlation have been found between the values of pK_a for 2,2′-bipyridines for aqueous solutions and their gas phase proton affinities.__________Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 3, 391–402, March, 2005.     

Time-resolved fluorescence analysis of the mobile flavin cofactor in p -hydroxybenzoate hydroxylase

Conformational heterogeneity of the FAD cofactor in p-hydroxybenzoate hydroxylase (PHBH) was investigated with time-resolved polarized flavin fluorescence. For binary enzyme/substrate (analogue) complexes of wild-type PHBH and Tyr222 mutants, crystallographic studies have revealed two distinct flavin conformations; the ‘in’ conformation with the isoalloxazine ring located in the active site, and the ‘out’ conformation with the isoalloxazine ring disposed towards the protein surface. Fluorescence-lifetime analysis of these complexes revealed similar lifetime distributions for the ‘in’ and ‘out’ conformations. The reason for this is twofold. First, the active site of PHBH contains various potential fluorescence-quenching sites close to the flavin. Fluorescence analysis of uncomplexed PHBH Y222V and Y222A showed that Tyr222 is responsible for picosecond fluorescence quenching free enzyme. In addition, other potential quenching sites, including a tryptophan and two tyrosines involved in substrate binding, are located nearby. Since the shortest distance between these quenching sites and the isoalloxazine ring differs only little on average, these aromatic residues are likely to contribute to fluorescence quenching. Second, the effect of flavin conformation on the fluorescence lifetime distribution is blurred by binding of the aromatic substrates: saturation with aromatic substrates induces highly efficient fluorescence quenching. The flavin conformation is therefore only reflected in the small relative contributions of the longer lifetimes.     

Comparative study on the inclusion behaviour of cyclodextrin derivatives with venoruton and rutin by thin layer chromatography

The interaction of rutin and venoruton (troxerutin), with alpha-, beta- and gamma-cyclodextrin (CD), hydroxypropyl-beta-cyclodextrin (HP-beta-CD) and methyl-beta-cyclodextrin (M-beta-CD) was investigated by reversed-phase thin layer chromatography on polyamide plates. A mobile phase consisted of NH(4)OH; NH(4)Cl buffer solution containing various CD concentrations (pH = 9.7, 20 degrees C) was used as mobile phase. The equilibrium constants (K(f)) and the retention factor (R(f)) were determined and used to study the inclusion process. The in fluence of CDs on the solubility of rutin and venoruton was characterized by R(M) values and the increasing hydrophilicity of drugs. The results show that the inclusion capacity of cyclodextrins follows the order HP-beta-CD > M-beta-CD > beta-CD > gamma-CD, and rutin is more easily included by the studied cyclodextrins than venoruton. In addition, the thermodynamic parameters (Delta H, Delta S) for the formation of complexes were obtained from the van't Hoff equation, displaying the enthalpy-entropy compensation effect.     

Comparison of thermodynamic properties of coarse-grained and atomic-level simulation models.

Thermodynamic data are often used to calibrate or test amomic-level (AL) force fields for molecular dynamics (MD) simulations. In contrast, the majority of coarse-grained (CG) force fields do not rely extensively on thermodynamic quantities. Recently, a CG force field for lipids, hydrocarbons, ions, and water, in which approximately four non-hydrogen atoms are mapped onto one interaction site, has been proposed and applied to study various aspects of lipid systems. To date, no extensive investigation of its capability to describe salvation thermodynamics has been undertaken. In the present study, a detailed picture of vaporization, solvation, and phase-partitioning thermodynamics for liquid hydrocarbons and water was obtained at CG and AL resolutions, in order to compare the two types or models and evaluate their ability to describe thermodynamic properties in the temperature range between 263 and 343 K. Both CG and AL models capture the experimental dependence of the thermodynamic properties on the temperature, albeit a systematically weaker dependence is found for the CG model. Moreover, deviations are found for solvation thermodynamics and for the corresponding enthalpy-entropy compensation for the CG model. Particularly water/oil repulsion seems to be overestimated. However, the results suggest that the thermodynamic properties considered should be reproducible by a CG model provided it is reparametrized on the basis of these liquid-phase properties.     

The third law of thermodynamics and the degeneracy of the ground state for lattice systems

The third law of thermodynamics, in the sense that the entropy per unit volume goes to zero as the temperature goes to zero, is investigated within the framework of statistical mechanics for quantum and classical lattice models. We present two main results: (i) For all models the question of whether the third law is satisfied can be decided completely in terms of ground-state degeneracies alone, provided these are computed for all possible boundary conditions. In principle, there is no need to investigate possible entropy contributions from low-lying excited states, (ii) The third law is shown to hold for ferromagnetic models by an analysis of the ground states.Dedicated to Pierre Résibois. Work supported in part by NSF grant PHY-7825390 A01.     

van der Waals corrections to density functional theory calculations: Methane,ethane, ethylene,benzene, formaldehyde,ammonia, water,PBE, and CPMD

Parameters are developed for a practical application of the empirical van der Waals (vdW) correction infrastructure available in the CPMD density functional theory (DFT) code. The binding energy, geometry, and potential energy surface (PES) are examined for methane, ethane, ethylene, formaldehyde, ammonia, three benzene dimer geometries, and three benzene–water geometries. The vdW corrected results compare favorably with MP2 and CCSD(T) calculations near the complete basis set limits, and with experimental results where they are available.