Strittmatter E F, Williams E R
Department of Chemistry, University of California, Berkeley, California 94720 Received: January 5, 2000; In Final Form: April 11, 2000.
J Phys Chem A. 2000 Jun 29;104(25):6069-76. doi: 10.1021/jp000038y.
The gas-phase structures and energetics of both protonated arginine dimer and protonated bradykinin were investigated using a combination of molecular mechanics with conformational searching to identify candidate low-energy structures, and density functional theory for subsequent minimization and energy calculations. For protonated arginine dimer, a good correlation (R = 0.88) was obtained between the molecular mechanics and EDF1 6-31+G* energies, indicating that mechanics with MMFF is suitable for finding low-energy conformers. For this ion, the salt-bridge or ion-zwitterion form was found to be 5.7 and 7.2 kcal/mol more stable than the simple protonated or ion-molecule form at the EDF1 6-31++G** and B3LYP 6-311++G** levels. For bradykinin, the correlation between the molecular mechanics and DFT energies was poor (R = 0.28), indicating that many low-energy structures are likely passed over in the mechanics conformational searching. This result suggests that structures of this larger peptide ion obtained using mechanics calculations alone are not necessarily reliable. The lowest energy structure of the salt-bridge form of bradykinin is 10.6 kcal/mol lower in energy (EDF1) than the lowest energy simple protonated form at the 6-311G* level. Similarly, the average energy of all salt-bridge structures investigated is 13.6 kcal/mol lower than the average of all the protonated forms investigated. To the extent that a sufficient number of structures are investigated, these results provide some additional support for the salt-bridge form of bradykinin in the gas phase.
采用分子力学与构象搜索相结合的方法来确定候选低能结构,并运用密度泛函理论进行后续的结构优化和能量计算,对质子化精氨酸二聚体和质子化缓激肽的气相结构及能量进行了研究。对于质子化精氨酸二聚体,分子力学与EDF1 6 - 31 + G能量之间获得了良好的相关性(R = 0.88),表明采用MMFF的分子力学方法适用于寻找低能构象。对于该离子,在EDF1 6 - 31++G和B3LYP 6 - 311++G水平下,盐桥或离子两性离子形式比简单的质子化或离子 - 分子形式分别稳定5.7和7.2 kcal/mol。对于缓激肽,分子力学与DFT能量之间的相关性较差(R = 0.28),这表明在分子力学构象搜索中可能遗漏了许多低能结构。该结果表明,仅使用分子力学计算获得的这种较大肽离子的结构不一定可靠。在6 - 311G水平下,缓激肽盐桥形式的最低能量结构比最低能量的简单质子化形式能量低10.6 kcal/mol(EDF1)。同样,所研究的所有盐桥结构的平均能量比所研究的所有质子化形式的平均能量低13.6 kcal/mol。在研究了足够数量结构的程度上,这些结果为气相中缓激肽的盐桥形式提供了一些额外的支持。