[논문 리뷰] Use of Polynomial Interpolation to Reduce Bias and Uncertainty of Free Energy Estimates via Thermodynamic Integration

Conrad Shyu, F. Marty Ytreberg|arXiv (Cornell University)|2008. 09. 05.

Phase Equilibria and Thermodynamics참고 문헌 24인용 수 2

한 줄 요약

이 논문은 열역학적 통합을 통한 자유에너지 차이 추정에서 편향과 불확실성을 줄이기 위해 등간격이 아닌 라미다 값과 다항식 보간을 사용하는 방법을 제안한다. 이 방법은 소프트코어 포텐셜이나 추가 시뮬레이션 없이도 용해 자유에너지 계산의 정확도와 정밀도를 향상시키며, 해석적으로 풀 수 있는 시스템과 소분자에서 검증되었다.

ABSTRACT

This report details an approach to improve the accuracy and precision of free energy difference estimates using thermodynamic integration data (slope of the free energy with respect to the switching variable lambda) and its application to calculating solvation free energy. The central idea is to utilize polynomial fitting schemes to approximate the thermodynamic integration data to improve the accuracy and precision of the free energy difference estimates. In this report we introduce polynomial and spline interpolation techniques. Two systems with analytically solvable relative free energies are used to test the accuracy and precision of the interpolation approach (Shyu and Ytreberg, J Comput Chem 30: 2297-2304, 2009). We also use both interpolation and extrapolation methods to determine a small molecule salvation free energy. Our simulations show that, using such polynomial techniques and non-equidistant lambda values, the solvation free energy can be estimated with high accuracy without using soft-core scaling and separate simulations for Lennard-Jones and partial charges. The results from our study suggest these polynomial techniques, especially with use of non-equidistant lambda values, improve the accuracy and precision for dF estimates without demanding additional simulations. To allow researchers to immediately utilize these methods, free software and documentation is provided via this http URL

연구 동기 및 목표

• 열역학적 통합 데이터로부터 자유에너지 차이 추정에서의 편향과 높은 불확실성 문제를 해결한다.
• 소프트코어 스케일링이나 별도의 시뮬레이션에 의존하지 않고도 용해 자유에너지 계산의 정확도와 정밀도를 향상시키는 방법을 개발한다.
• 해석적으로 풀 수 있는 시스템에서 다항식 및 스퍼플린 보간 기법의 성능을 평가한다.
• 비등간격 라미다 간격이 추정 품질 향상에 기여하는지를 입증한다.
• 연구 공동체가 즉시 활용할 수 있도록 오픈소스 소프트웨어를 제공한다.

제안 방법

• 자유에너지가 스위칭 파라미터 라미다의 함수임을 가정하고 열역학적 통합 데이터에 다항식 피팅을 적용한다.
• 자유에너지 곡률이 높은 영역을 더 잘 포착하기 위해 등간격이 아닌 라미다 값을 사용하여 보간 정확도를 향상시킨다.
• 시뮬레이션 데이터로부터 자유에너지 차이를 추정하기 위해 보간 및 외삽 기법을 모두 구현한다.
• 상대 자유에너지가 해석적으로 구할 수 있는 두 시스템에서 방법을 검증하여 편향과 정밀도를 평가한다.
• 표준 방법과의 비교를 통해 정확도 향상과 불확실성 감소 정도를 정량화한다.
• 재현성과 공동체의 도입을 지원하기 위해 무료 소프트웨어와 문서를 배포한다.

실험 결과

연구 질문

• RQ1 열역학적 통합에서 유도된 자유에너지 차이 추정에서 다항식 보간은 편향과 불확실성을 줄일 수 있는가?
• RQ2 등간격 간격 대비 비등간격 라미다 값 사용이 자유에너지 추정의 정확도와 정밀도에 어떤 영향을 미치는가?
• RQ3 이 방법은 레너드존스 및 쿨롱 상호작용 항에 대해 소프트코어 포텐셜과 별도의 시뮬레이션을 완전히 제거할 수 있는가?
• RQ4 이 보간 기법은 해석적으로 알려진 자유에너지 차이를 가진 시스템에서 얼마나 잘 작동하는가?
• RQ5 이 방법은 실생활 용해 자유에너지 계산, 예를 들어 소분자에 대해 신뢰성 있게 적용될 수 있는가?

주요 결과

• 등간격이 아닌 라미다 값과 다항식 보간을 사용하면 자유에너지 차이 추정에서 편향과 불확실성이 크게 감소한다.
• 이 방법은 레너드존스 및 부분 전하 항에 대해 소프트코어 스케일링이나 별도의 시뮬레이션 없이도 높은 정확도의 용해 자유에너지 추정을 달성한다.
• 해석적으로 알려진 자유에너지 차이가 있는 시험 시스템에서 표준 방법보다 이 기법이 뛰어난 성능을 보인다.
• 비등간격 라미다 간격은 특히 자유에너지 변화가 급격한 영역에서 보간의 정밀도를 향상시킨다.
• 최적화된 라미다 간격을 사용하는 단일 시뮬레이션 세트만으로도 정밀한 용해 자유에너지 추정이 가능하다.
• 저자들은 자유롭게 이용 가능한 소프트웨어를 제공하여 즉각적인 구현과 보다 넓은 도입을 가능하게 하였다.