[논문 리뷰] Established DFT methods calculation of conjugation disturbed in the presence of torsional hyperconjugation

John N. Sharley|arXiv (Cornell University)|2015. 12. 14.

Advanced Chemical Physics Studies참고 문헌 66인용 수 4

한 줄 요약

이 논문은 고전적인 DFT 방법이 비결정적 토크르성 허공공명을 올바르게 처리하지 못함으로써 카보닐 기능에서 시그마/파이 오비탈 분리가 붕괴되는 중요한 실패를 밝혀낸다. NBO 분석을 통해, 단백질 계산에서 흔히 사용되는 작은 기저함수를 사용한 DFT는 적절한 전자 구조를 유지하지 못함을 보여주며, 이는 β-시트와 같은 2차 단백질 구조에서 수소 결합과 공간적 상호작용을 잘못 예측하게 한다.

ABSTRACT

Accurate treatment of amide resonance is important in electronic structure calculation of protein, for Resonance-Assisted Hydrogen Bonding in the hydrogen bonded chains of backbone amides of protein secondary structures such as beta sheets and alpha helices is determined by amide resonance. Variation in amide resonance is the means by which the hydrogen bonding in these chains is cooperative. Amide carbonyl orbitals are revealed by Natural Bond Orbital, NBO, analysis to substantially maintain sigma/pi separation in the presence of torsional hyperconjugative interactions with wavefunction methods but not with established Density Functional Theory, DFT, methods. This DFT error is most pronounced with small basis sets such as are used with DFT for proteins to reduce the basis function count. This error disturbs calculation of a range of amide donor-acceptor and steric interactions. This finding has important implications for the selection of electronic structure methods and basis sets for protein calculations. For example, great caution is needed in interpreting the results of applying established DFT methods to proteins containing any beta sheets. We recommend that every protein DFT calculation be accompanied by NBO assessment of maintenance of amide carbonyl sigma/pi separation and absence of carbonyl bond bending. Further, we propose that these metrics be standard benchmarks of electronic structure methods and basis sets.

연구 동기 및 목표

• 고전적인 DFT 방법이 단백질 2차 구조에서 아미드 공명을 모델링하는 데 정확도를 갖는지 조사하기.
• 비결정적 토크르성 허공공명이 DFT 계산에서 전자 구조를 어떻게 왜곡시키는지, 특히 작은 기저함수를 사용할 경우를 중심으로 밝히기.
• DFT가 β-시트와 α-헤리스에 있는 수소 결합과 공간적 상호작용을 예측하는 데 신뢰성 있는지 평가하기.
• NBO 기반 진단 지표를 단백질 모델링에서 전자 구조 방법을 검증하는 데 필수적인 기준으로 제안하기.
• 공명과 오비탈 혼합을 잘못 처리함으로써 DFT 기반 단백질 시뮬레이션에서 오해의 소지가 있는 결론이 도출될 위험을 다루기.

제안 방법

• 아미드 카보닐 기능에서 시그마/파이 오비탈 분리를 평가하기 위해 자연 결합 오비탈(NBO) 분석을 수행하였다.
• 정확히 오비탈 분리를 유지하는 파동함수 방법의 결과와 기존 DFT 함수의 결과를 비교하였다.
• 단백질 DFT 계산에서 일반적으로 사용되는 작은 기저함수를 사용하여 실제 계산 조건을 시뮬레이션하였다.
• 다양한 이황각에서 아미드 N–H와 C=O 오비탈 간의 비결정적 허공공명 상호작용을 분석하였다.
• DFT 오류가 결합 휘어짐과 기여자-수용자 상호작용과 같은 핵심 전자 구조적 성질에 미치는 영향을 평가하였다.
• 특히 카보닐 오비탈 분리와 결합 휘어짐을 포함한 NBO 지표를 단백질 생물분자 응용 분야에서 DFT 방법의 표준 검증 절차로 포함시키는 것을 제안하였다.

실험 결과

연구 질문

• RQ1 비결정적 허공공명이 존재할 때, 고전적인 DFT 방법이 아미드 카보닐에서 시그마와 파이 오비탈 분리를 얼마나 잘 유지하지 못하는가?
• RQ2 작은 기저함수의 사용이 단백질에서 아미드 공명을 모델링할 때 DFT 오류를 얼마나 악화시키는가?
• RQ3 잘못된 오비탈 혼합이 β-시트와 α-헤리스에서 수소 결합의 공조성 예측에 어떤 영향을 미치는가?
• RQ4 NBO 분석은 DFT에 의해 유도된 아미드 전자 구조 오류를 신뢰성 있게 감지할 수 있는가?
• RQ5 NBO 진단 지표는 단백질 시뮬레이션을 위한 전자 구조 방법 검증에 얼마나 표준화되어야 하는가?

주요 결과

• 비결정적 토크르성 허공공명이 존재할 경우, 고전적인 DFT 방법은 아미드 카보닐에서 시그마와 파이 오비탈 간의 분리를 제대로 유지하지 못한다.
• 이 오류는 계산 비용 제약으로 인해 단백질 DFT 계산에서 일반적으로 사용되는 작은 기저함수를 사용할 경우 가장 심각해진다.
• 파동함수 방법은 정확히 오비탈 분리를 유지하므로, 오류는 DFT 함수와 기저함수의 한계에 기인한 것이며 특정한 DFT 기반 오류임을 시사한다.
• DFT 오류는 시그마와 파이 성분이 인위적으로 혼합되게 하여 아미드 기능의 전자 구조를 왜곡시킨다.
• DFT 계산에서는 인위적으로 카보닐 결합 휘어짐이 유도되며, 이는 수소 결합과 공간적 상호작용 예측의 정확도에 영향을 미친다.
• NBO 분석은 DFT로 계산된 단백질 시스템에서 이러한 전자 구조 오류를 감지하는 신뢰할 수 있는 진단 도구로 밝혀졌다.