Describing weak interactions of biomolecules with dispersion-corrected density functional theory

Interaction energies of the biomolecules in the JSCH-2005 database are calculated with density functional theory using the exchange-correlation functional BLYP augmented with dispersion-corrected atom-centered potentials (DCACPs). The results are in excellent agreement with extrapolated CCSD(T) complete basis set limit references with unsigned mean errors of less than 1.6 kcal mol(-1). Geometry optimisations all reach stable configurations that are close to the MP2-optimised reference geometries.     

Supramolecular binding thermodynamics by dispersion-corrected density functional theory

The equilibrium association free enthalpies ΔG(a) for typical supramolecular complexes in solution are calculated by ab initio quantum chemical methods. Ten neutral and three positively charged complexes with experimental ΔG(a) values in the range 0 to -21?kcal?mol(-1) (on average -6?kcal?mol(-1) ) are investigated. The theoretical approach employs a (nondynamic) single-structure model, but computes the various energy terms accurately without any special empirical adjustments. Dispersion corrected density functional theory (DFT-D3) with extended basis sets (triple-ζ and quadruple-ζ quality) is used to determine structures and gas-phase interaction energies (ΔE), the COSMO-RS continuum solvation model (based on DFT data) provides solvation free enthalpies and the remaining ro-vibrational enthalpic/entropic contributions are obtained from harmonic frequency calculations. Low-lying vibrational modes are treated by a free-rotor approximation. The accurate account of London dispersion interactions is mandatory with contributions in the range -5 to -60?kcal?mol(-1) (up to 200?% of ΔE). Inclusion of three-body dispersion effects improves the results considerably. A semilocal (TPSS) and a hybrid density functional (PW6B95) have been tested. Although the ΔG(a) values result as a sum of individually large terms with opposite sign (ΔE vs. solvation and entropy change), the approach provides unprecedented accuracy for ΔG(a) values with errors of only 2?kcal?mol(-1) on average. Relative affinities for different guests inside the same host are always obtained correctly. The procedure is suggested as a predictive tool in supramolecular chemistry and can be applied routinely to semirigid systems with 300-400 atoms. The various contributions to binding and enthalpy-entropy compensations are discussed.     

Performance of dispersion-corrected density functional theory for the interactions in ionic liquids

Potential energy curves for the dissociation of cation-anion associates representing the building units of ionic liquids have been computed with dispersion corrected DFT methods. Non-local van der Waals density functionals (DFT-NL) for the first time as well as an atom pair-wise correction method (DFT-D3) have been tested. Reference data have been computed at the extrapolated MP2/CBS and estimated CCSD(T)/CBS levels of theory. The investigated systems are combined from two cations (1-butyl-3-methylimidazolium and tributyl(methyl)posphonium) and three anions (chloride, dicyanamide, acetate). We find substantial stabilization from London dispersion energy near equilibrium of 5-7 kcal mol(-1) (about 5-6% of the interaction energy). Equilibrium distances are shortened by 0.03-0.09 ? and fundamental (inter-fragment) vibrational frequencies (which are in the range 140-180 cm(-1)) are increased by typically 10 cm(-1) when dispersion corrections are made. The dispersion-corrected hybrid functional potentials are in general in excellent agreement with the corresponding CCSD(T) reference data (typical deviations of about 1-2%). The DFT-D3 method performs unexpectedly well presumably because of cancellation of errors between the dispersion coefficients of the cations and anions. Due to self-interaction error, semi-local density functionals exhibit severe SCF convergence problems, and provide artificial charge-transfer and inaccurate interaction energies for larger inter-fragment distances. Although these problems may be alleviated in condensed phase simulations by effective Coulomb screening, only dispersion-corrected hybrid functionals with larger amounts of Fock-exchange can in general be recommended for such ionic systems.     

Validation of dispersion-corrected density functional theory approaches for ionic liquid systems

The performance of several general gradient approximation, meta general gradient approximation, and hybrid functionals is tested against M?ller-Plesset perturbation theory second-order for ionic liquid systems. Additionally, two dispersion-corrected approaches (addition of van der Waals forces by a 1/r(6) term and employing a dispersion-corrected atom-center dispersion pseudopotential) were studied. For the 1-butyl-3-methylimidazolium cation neglecting dispersion results in different trends for structural stabilities. The two applied correction schemes for density functional theory improve the results tremendously. Investigating several 1-butyl-3-methylimidazolium dicianamide ion pairs shows a mean absolute deviation from M?ller-Plesset perturbation theory of 35.7 kJ/mol for Hartree-Fock and up to 33.2 kJ/mol for the density functional theory methods. The dispersion-corrected methods reduce the mean absolute deviation to less than 10 kJ/mol. Comparing adducts of the 1-ethyl-3-methylimidazolium dicianamide ion pair with Diels-Alder educts (cyclopentadiene and methylacrylate) shows similar energetic differences as for the ion pairs. Furthermore large deviations in geometries for the intermolecular distances were found for the Hartree-Fock approach (mean absolute deviation: 190 pm) and density functional theory (mean absolute deviation up to 178 pm) while for the dispersion-corrected methods the mean absolute deviation is less than 50 pm.     

Searching of potential energy curves for the benzene dimer using dispersion-corrected density functional theory

The present work aims to establish the utility of dispersion-corrected density functional theory for potential energy curves of the benzene dimer, a problem that has received significant attention for a long time. The interaction energies of parallel-stacked, T-shaped and parallel-displaced benzene dimer configurations have been evaluated using both dispersion- and normal gradient-corrected Perdew-Burke-Ernzerhof functionals along with Dunning's augmented correlation-consistent polarized valence triple-zeta (aug-cc-pVTZ) basis functions and compared with explicit correlation methods. The potential energy curves for the parallel-stacked and parallel-displaced benzene dimers are in excellent agreement with highly accurate coupled cluster (CCSD(T)) results, while for the T-shaped benzene dimer the dispersion-corrected results show a distinct deviation, being closer in that case to the MP2 level of results. The overestimation of interaction energy in the T-shaped dimer may be attributed to the presence of a permanent dipole moment in this configuration and indicates a structural dependence of the dispersion-corrected density functional method.     

低阶煤中氢键作用的色散校正密度泛函理论研究

本研究以苯酚…苯酚、苯酚…苯、苯酚…二苯醚、苯酚…喹啉和苯甲酸…苯甲酸为对象，采用色散校正的密度泛函理论分别研究褐煤中自缔合OH、OH-π、OH-醚O、OH-N和COOH-COOH之间形成的氢键。此外，还研究了氢键供体中取代基(CH₃-、CH₃O-、OH-、NH₂-、COOH-和NO₂-)对氢键的影响。对上述复合物进行了几何优化，并计算了能量、Mulliken电荷分布及振动频率。从优化的结构中可以看出上述复合物之间都存在氢键，所有复合物中O-H键键长都比苯酚中自由羟基的长，这表明这些复合物之间存在相互作用。其中，羧酸…羧酸复合物中O-H键的键长最长。此外，通过Mulliken电荷分布可看出上述复合物之间存在电荷转移。基于振动频率分析，所有的O-H键伸缩振动都发生了红移，尤其是羧酸…羧酸和苯酚…喹啉复合物，这可为煤中羟基振动的红外光谱分析提供依据。根据键能不同氢键强度按以下顺序依次递减：COOH-COOH>OH-N > 自缔合OH≈OH-醚O > OH-π，这与振动频率的分析结果一致。此外，不同取代基对氢键作用的影响不同。     

Effects of solvent on weak halogen bonds: Density functional theory calculations

In recent years, many applications of solution‐phase halogen bonding in anion recognition, catalysis, and pseudorotaxane formation have been reported. Moreover, a number of thermodynamic data of halogen bonding interactions in organic solution are now available. To obtain detailed information of the influence of the surrounding medium on weak halogen bonds, a series of dimeric complexes of halobenzene (PhX) with three electron donors (H₂O, HCHO, and NH₃) were investigated by means of DFT/PBE calculations in this work. The PCM implicit solvation approach was utilized to include the effects of three solvents (cyclohexane, chloroform, and water) as representatives for a wide range of dielectric constant. In some cases, halogen‐bond distances are shown to shorten in solution, accompanied by concomitant elongation of the C? X bonds. For the remaining systems, the intermolecular distances tend to increase or remain almost unchanged under solvent effects. In general, the solvent has a slight destabilizing effect on weak halogen bonds; the strength order of halogen bonds observed in vacuum remains unchanged in liquid phases. Particularly, the interaction strength attenuates in the order I > Br > Cl in solution, consistent with the experimental measurements of weak halogen bond door abilities. The similarities between halogen and hydrogen bonding in solution were also elucidated. The results presented herein would be very useful in future applications of halogen bonding in molecular recognition and medicinal chemistry. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Calculated vibrational spectra of weak hydrogen bonds

Potential energy and dipole moment surfaces for the H-bond SH... S in the dimeric methanethiol have been calculated by the SCF-MO-LCGO method, and the vibrational spectrum — transition frequencies and IR absorption intensities at 20 °K — computed. This spectrum is compared with that of the monomeric species and with experimental results. The resulting dimerization energy is 1.4 kcal/mole.     

CO2 adsorption on TiO2(101) anatase: a dispersion-corrected density functional theory study

Adsorption, diffusion, and dissociation of CO(2) on the anatase (101) surface were investigated using dispersion-corrected density functional theory. On the oxidized surface several different local minima were identified of which the most stable corresponds to a CO(2) molecule adsorbed at a five-fold coordinated Ti site in a tilted configuration. Surface diffusion is characterized by relatively small activation barriers. Preferential diffusion takes place along Ti rows and involves a cartwheel type of motion. The presence of a bridging oxygen defect or a surface interstitial Ti atom allows creation of several new strong binding configurations the most stable of which have bent CO(2) structures with simultaneous bonding to two surface Ti atoms. Subsurface oxygen vacancy or interstitial Ti defects are found to enhance the bonding of CO(2) molecules to the surface. CO(2) dissociation from these defect sites is calculated to be exothermic with barriers less than 21 kcal/mol. The use of such defects for catalytic activation of CO(2) on anatase (101) surface would require a mechanism for their regeneration.     

The influence of solvation on short strong hydrogen bonds: a density functional theory study of the Asp-His interaction in subtilisins

The effect of a structural water molecule on the electronic nature of the His64-Asp32 hydrogen bond in subtilisins is examined by DFT calculations; the structural water is found to favor a short strong hydrogen bond in the catalytic triad in sharp contrast to some current beliefs.     

The performance and relationship among range-separated schemes for density functional theory

The performance and relationship among different range-separated (RS) hybrid functional schemes are examined using the Coulomb-attenuating method (CAM) with different values for the fractions of exact Hartree-Fock (HF) exchange (α), long-range HF (β), and a range-separation parameter (μ), where the cases of α + β = 1 and α + β = 0 were designated as CA and CA0, respectively. Attenuated PBE exchange-correlation functionals with α = 0.20 and μ = 0.20 (CA-PBE) and α = 0.25 and μ = 0.11 (CA0-PBE) are closely related to the LRC-ωPBEh and HSE functionals, respectively. Time-dependent density functional theory calculations were carried out for a number of classes of molecules with varying degrees of charge-transfer (CT) character to provide an assessment of the accuracy of excitation energies from the CA functionals and a number of other functionals with different exchange hole models. Functionals that provided reasonable estimates for local and short-range CT transitions were found to give large errors for long-range CT excitations. In contrast, functionals that afforded accurate long-range CT excitation energies significantly overestimated energies for short-range CT and local transitions. The effects of exchange hole models and parameters developed for RS functionals for CT excitations were analyzed in detail. The comparative analysis across compound classes provides a useful benchmark for CT excitations.     

The role of weak hydrogen bonds in chiral recognition

Chiral recognition has been studied in neutral or ionic weakly bound complexes isolated in the gas phase by combining laser spectroscopy and quantum chemical calculations. Neutral complexes of the two enantiomers of lactic ester derivatives with chiral chromophores have been formed in a supersonic expansion. Their structure has been elucidated by means of IR-UV double resonance spectroscopy in the 3 μm region. In both systems described here, the main interaction ensuring the cohesion of the complex is a strong hydrogen bond between the chromophore and methyl-lactate. However, an additional hydrogen bond of much weaker strength plays a discriminative role between the two enantiomers. For example, the 1:1 heterochiral complex between R-(+)-2-naphthyl-ethanol and S-(+) methyl-lactate is observed, in contrast with the 1:1 homochiral complex which lacks this additional hydrogen bond. On the other hand, the same kind of insertion structures is formed for the complex between S-(±)-cis-1-amino-indan-2-ol and the two enantiomers of methyl-lactate, but an additional addition complex is formed for R-methyl-lactate only. This selectivity rests on the formation of a weak CHπ interaction which is not possible for the other enantiomer. The protonated dimers of Cinchona alkaloids, namely quinine, quinidine, cinchonine and cinchonidine, have been isolated in an ion trap and studied by IRMPD spectroscopy in the region of the ν(OH) and ν(NH) stretch modes. The protonation site is located on the alkaloid nitrogen which acts as a strong hydrogen bond donor in all the dimers studied. While the nature of the intermolecular hydrogen bond is similar in the homochiral and heterochiral complexes, the heterochiral complex displays an additional weak CHO hydrogen bond located on its neutral part, which results in slightly different spectroscopic fingerprints in the ν(OH) stretch region. This first spectroscopic evidence of chiral recognition in protonated dimers opens the way to the study of the complexes of Cinchona alkaloids involved in enantioselective catalysis. These examples show how secondary hydrogen bonds controlled by stereochemical factors govern molecular recognition processes.     

A density functional theory study of hydrogen adsorption in MOF-5

Ab initio molecular dynamics in the generalized gradient approximation to density functional theory and ground-state relaxations are used to study the interaction between molecular hydrogen and the metal-organic framework with formula unit Zn4O(O2C-C6H4-CO2)3. Five symmetrically unique adsorption sites are identified, and calculations indicate that the sites with the strongest interaction with hydrogen are located near the Zn4O clusters. Twenty total adsorption sites are found around each Zn4O cluster, but after 16 of these are populated, the interaction energy at the remaining four sites falls off significantly. The adsorption of hydrogen on the pore walls creates an attractive potential well for hydrogen in the center of the pore. The effect of the framework on the physical structure and electronic structure of the organic linker is calculated, suggesting ways by which the interaction between the framework and hydrogen could be modified.     

Performance of density functional theory methods to describe intramolecular hydrogen shifts

The performance of three exchange and correlation density functionals, LDA, BLYP and B3LYP, with four basis sets is tested in three intramolecular hydrogen shift reactions. The best reaction and activation energies come from the hybrid functional B3LYP with triple-ζ basis sets, when they are compared with high-level post-Hartree-Fock results from the literature. For a fixed molecular geometry, the electrophilic Fukui function is computed from a finite difference approximation. Fukui function shows a small dependence with both the exchange and correlation functional and the basis set. Evolution of the Fukui function along the reaction path describes important changes in the basic sites of the corresponding molecules. These results are in agreement with the chemical behavior of those species.     

Comparison of random-walk density functional theory to simulation for bead-spring homopolymer melts

Density profiles for a homopolymer melt near a surface are calculated using a random-walk polymeric density functional theory, and compared to results from molecular dynamics simulations. All interactions are of a Lennard-Jones form, for both monomer-monomer interactions and surface-monomer interactions, rather than the hard core interactions which have been most investigated in the literature. For repulsive systems, the theory somewhat overpredicts the density oscillations near a surface. Nevertheless, near quantitative agreement with simulation can be obtained with an empirical scaling of the direct correlation function. Use of the random phase approximation to treat attractive interactions between polymer chains gives reasonable agreement with simulation of dense liquids near neutral and attractive surfaces.     

Comparative study of weak interactions in molecular crystals: H-H bonds vs hydrogen bonds

The crystal structures of tetraphenylphosphonium squarate, bianthrone, and bis(benzophenone)azine are shown to contain a variety of C-H(delta+)...(delta+)H-C interactions, as well as a variety of C-H...O and C-H...C(pi) interactions. Each of these molecules possesses interactions that can possibly be characterized as either H-H bonds or weak hydrogen bonds based on the first four criteria proposed by Koch and Popelier. These interactions have been completely characterized topologically after the multipole refinement of the structures. It appears that weak interactions of the form C-H(delta+)...(delta+)H-C possess certain correlations between the various properties of the electron density at the bond critical points. The coexistence of the three types of interactions makes it possible to establish similarities and differences in the correlations of these weak interactions. This all leads to a better understanding of H-H interactions and how they fit into the hierarchy of weak interactions.