Zhang Yan, Wang Chang-Sheng
Department of Chemistry, Liaoning Normal University, Dalian 116029, People's Republic of China.
J Comput Chem. 2009 Jun;30(8):1251-60. doi: 10.1002/jcc.21141.
Computation of accurate intramolecular hydrogen-bonding energies for peptides is of great importance in understanding the conformational stabilities of peptides and developing a more accurate force field for proteins. We have proposed a method to determine the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides. In this article, the method is further applied to evaluate the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in peptides. The optimal structures of the intramolecular 10-membered ring N-H...O=C hydrogen bonds in glycine and alanine tripetide molecules are obtained at the MP2 level with 6-31G(d), 6-31G(d,p), and 6-31+G(d,p) basis sets. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are then evaluated based on our method at the MP2/6-311++G(3df,2p) level with basis set superposition error correction. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are calculated to be in the range of -6.84 to -7.66, -4.44 to -4.98, and -6.95 to -7.88 kcal/mol. The method is also applied to estimate the individual intermolecular hydrogen-bonding energies in the dimers of amino-acetaldehyde, 2-amino-acetamide, formamide, and oxalamide, each dimer having two identical intermolecular hydrogen bonds. According to our method, the individual intermolecular hydrogen-bonding energies in the four dimers are calculated to be -1.77, -1.67, -6.35, and -4.82 kcal/mol at the MP2/6-311++G(d,p) level, which are in good agreement with the values of -1.84, -1.72, -6.23, and -4.93 kcal/mol predicted by the supermolecular method.
计算肽的精确分子内氢键能对于理解肽的构象稳定性以及开发更精确的蛋白质力场至关重要。我们提出了一种确定甘氨酸和丙氨酸肽中分子内七元环N-H...O=C氢键能的方法。在本文中,该方法进一步应用于评估肽中分子内十元环N-H...O=C氢键能。在MP2水平下,使用6-31G(d)、6-31G(d,p)和6-31+G(d,p)基组获得了甘氨酸和丙氨酸三肽分子中分子内十元环N-H...O=C氢键的最佳结构。然后基于我们的方法,在MP2/6-311++G(3df,2p)水平下并进行基组叠加误差校正,评估分子内十元环N-H...O=C氢键能。计算得到分子内十元环N-H...O=C氢键能在-6.84至-7.66、-4.44至-4.98和-6.95至-7.88 kcal/mol范围内。该方法还应用于估计氨基乙醛、2-氨基乙酰胺、甲酰胺和草酰胺二聚体中各个分子间氢键能,每个二聚体具有两个相同的分子间氢键。根据我们的方法,在MP2/6-311++G(d,p)水平下,四个二聚体中各个分子间氢键能计算为-1.77、-1.67、-6.35和-4.82 kcal/mol,与超分子方法预测的-1.84、-1.72、-6.2和-4.93 kcal/mol值吻合良好。