Accurate Reaction Energies in Proteins Obtained by Combining QM/MM and Large QM Calculations

详细信息    查看全文

• 作者：LiHong Hu ; P盲r S枚derhjelm ; Ulf Ryde
• 刊名：Journal of Chemical Theory and Computation
• 出版年：2013
• 出版时间：January 8, 2013
• 年：2013
• 卷：9
• 期：1
• 页码：640-649
• 全文大小：352K
• 年卷期：v.9,no.1(January 8, 2013)
• ISSN：1549-9626

文摘

We here suggest and test a new method to obtain stable energies in proteins for charge-neutral reactions by running large quantum mechanical (QM) calculations on structures obtained by combined QM and molecular mechanics (QM/MM) geometry optimization on several snapshots from molecular dynamics simulations. As a test case, we use a proton transfer between a metal-bound cysteine residue and a second-sphere histidine residue in the active site of [Ni,Fe] hydrogenase, which has been shown to be very sensitive to the surroundings. We include in the QM calculations all residues within 4.5 脜 of the active site, two capped residues on each side of the active-site residues, and all charged groups that are buried inside the protein, which for this enzyme includes three iron鈥搒ulfur clusters, in total, 930 atoms. These calculations are performed at the BP86/def2-SV(P) level, but the energies are then extrapolated to the B3LYP/def2-TZVP level with a smaller QM system, and zero-point energy, entropy, and thermal effects are added. We test three approaches to model the remaining atoms of the protein solvent, viz., by standard QM/MM approaches using either mechanical or electrostatic embedding or by using a continuum solvation model for the large QM systems. Quite encouragingly, the three approaches give the same results within 14 kJ/mol, and variations in the size of the QM system do not change the energies by more than 8 kJ/mol, provided that the QM/MM junctions are not moved closer to the QM system. The statistical precision for the average over 10 snapshots is 1鈥? kJ/mol.