Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains
Peter I Nagy, Paul W Erhardt
Journal of Physical Chemistry. A 2008 May 8, 112 (18): 4342-54
18373368
Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.
Full Text Links
Find Full Text Links for this Article
You are not logged in. Sign Up or Log In to join the discussion.