[논문 리뷰] Bridge functional for the molecular density functional theory with consistent pressure and surface tension and its importance for solvation in water

Cédric Gageat, Luc Belloni외 2명|arXiv (Cornell University)|2017. 09. 28.

Protein Structure and Dynamics참고 문헌 6인용 수 3

한 줄 요약

이 논문은 분자 밀도 기반 함수 이론(MDFT)를 위한 새로운 군집화된 다이어그램 기반 기능을 소개하며, 이는 액체-기화 공존, 균질 압력, 표면 장력 등을 정확히 재현함으로써 열역학적 일致성을 확보한다. 이 방법은 실험치와 1 kcal/mol 이내의 용해 자유 에너지를 달성하고, 분석 시간이 분 단위로 끝나는 빠른 속도로 정확한 용매 밀도 프로파일을 제공한다. 이는 초과 연결 체인 근사법(hypernetted-chain approximation)을 크게 뛰어넘는 성능를 보이며, 수치적 효율성은 유지한다.

ABSTRACT

We address the problem of predicting the solvation free energy and equilibrium solvent density profile in fews minutes from the molecular density functional theory beyond the usual hypernetted-chain approximation. We introduce a bridge functional of a coarse-grained, weighted solvent density. In few minutes at most, for solutes of sizes ranging from small compounds to large proteins, we produce (i) an estimation of the free energy of solvation within 1 kcal/mol of the experimental data for the hydrophobic solutes presented here, and (ii) the solvent distribution around the solute. Contrary to previous propositions, this bridge functional is thermodynamically consistent in that it produces the correct liquid-vapor coexistence and the experimental surface tension. We show this consistency to be of crucial importance for water at room temperature and pressure. This bridge functional is designed to be simple, local, and thus numerically efficient. Finally, we illustrate this new level of molecular theory of solutions with the study of the hydration shell of a protein.

연구 동기 및 목표

• 액체-기화 공존과 물에서 표면 장력을 정확히 반영하는 열역학적으로 일관된 MDFT 브리지 기능을 개발한다.
• 기본적인 초과 연결 체인(HNC) 근사법을 뛰어넘어 MDFT에서 용해 자유 에너지 및 용매 밀도 프로파일 예측의 정확도를 향상시킨다.
• 소분자에서 단백질에 이르기까지 대규모 생고분자에 적용 가능한 빠른 계산 속도(분 단위)를 확보하면서도 높은 정확도를 유지한다.
• 푸리에 공간에서 국소적 기능으로 설계하여 푸리에 변환(FFT)을 활용해 수치적 효율성을 확보한다.
• 수소 결합 화합물과 단백질 4m7g을 포함한 다양한 용질에 대해 실험 데이터 및 기준 분자 동역학 시뮬레이션과의 비교를 통해 방법을 검증한다.

제안 방법

• 비국소적 상관관계를 수치적으로 효율적으로 기술하기 위해 군집화된 가중 평균 용매 밀도를 기반으로 브리지 기능을 구성한다.
• 푸리에 공간에서 국소적 기능으로 설계되어 빠른 푸리에 변환(FFT)을 통해 신속하게 평가할 수 있다.
• 초과, 이상, 외부 기여를 포함한 총 자유 에너지 기능을 최소화하는 방식으로 MDFT 프레임워크에 통합한다.
• 변분 접근법을 사용하여 분자 밀도 ρ(r,ω)에 대해 자유 에너지 기능의 합을 최소화함으로써 용매 밀도 프로파일을 결정한다.
• 정확한 액체-기화 공존과 표면 장력을 재현하도록 캘리브레이션하여 열역학적 일관성을 확보한다.
• 용해 자유 에너지 및 방사 분포 함수에 대해 기준 분자 동역학 시뮬레이션과 실험 데이터와의 비교를 통해 방법을 검증한다.

실험 결과

연구 질문

• RQ1 MDFT에서 브리지 기능을 설계하여 열역학적 일관성을 확보할 수 있는가? 특히 물에서 정확한 액체-기화 공존과 표면 장력을 반영할 수 있는가?
• RQ2 제안된 군집화된 브리지 기능은 기존 HNC 근사법에 비해 용해 자유 에너지 예측 정확도를 크게 향상시키는가?
• RQ3 분 단위 내에 대규모 생고분자(예: 단백질)에 대해서도 높은 정확도로 용매 밀도 프로파일과 용해 자유 에너지를 예측할 수 있는가?
• RQ4 새로운 브리지 기능은 수용성 껍질의 구조를 예측하는 데 있어 기준 분자 동역학 시뮬레이션과 비교해 어떤가?
• RQ5 큰 수소 결합 화합물에 대해 기존의 MDFT-HNC 접근법에서 흔히 나타나는 비정상적인 피크(prepeaks)를 브리지 기능이 얼마나 잘 제거하는가?

주요 결과

• 제안된 브리지 기능은 용해 자유 에너지 예측의 평균 절대 오차(MUE)를 HNC의 23.61 kJ/mol에서 3.88 kJ/mol으로 감소시켜, 수소 결합 화합물에 대해 실험치와 1 kcal/mol 이내로 일치시킨다.
• 이전의 브리지 기능이 실패한 바 있는 물의 표면 장력과 액체-기화 공존 조건을 정확히 재현하여, 정확한 수소 결합 모델링에 필수적인 열역학적 일관성을 확보한다.
• 메탄, neon, 아르곤, 크립톤, 크세논에 대해 용해 자유 에너지 예측이 실험치와 1 kcal/mol 이내이며, MUE는 23.61에서 3.88 kJ/mol으로 감소한다.
• 단백질 4m7g의 수용성 껍질은 신규 브리지 기능을 사용한 MDFT로 20분 내에 계산되었으며, 3일간의 분자 동역학 시뮬레이션과 3× 밀도 기준 등면에서 강한 일치를 보였다.
• 기존의 MDFT-HNC 접근법에서 흔히 나타나는 큰 경직된 구형 용질에 대한 방사 분포 함수의 비정상적인 피크를 브리지 기능이 제거한다.
• 모든 원자 기반 분자 동역학에 비해 계산 비용을 4개 지수 감소시켜 고속이고 정밀한 용해 분석을 가능하게 한다.