차분 에너지 표면 분석을 이용한 폴딩 지형 최적 재구성

초록

본 논문은 여러 개의 조화적 바이어스 포텐셜 하에서 수집된 단일 분자 실험·시뮬레이션 데이터를 이용해 자유에너지 표면의 기울기(dE/dx)를 직접 추정하는 차분 에너지 표면 분석(DESA) 방법을 제시한다. DESA는 최대우도 추정에 기반하며, 기존의 가중 히스토그램 분석(WHAM)과 동일한 정확도를 제공한다. 또한 두 가지 자체 일관성 검증 기준을 도입해 반응좌표가 자유에너지의 단일값 함수를 이루는지를 판단한다. DNA 헤어핀의 말단‑말단 거리 측정 실험에 적용한 결과, DESA가 WHAM과 일치하는 에너지 지형을 복원함을 확인했으며, 반응좌표로서 거리 변수가 적절함을 검증하였다.

상세 요약

DESA는 전통적인 WHAM이 전체 자유에너지 E(x)를 직접 복원하는 방식과 달리, 에너지 표면의 미분값인 dE/dx를 개별 바이어스 실험에서 얻은 확률 분포의 로그 기울기로부터 추정한다. 구체적으로, 각 실험 i는 조화적 바이어스 V_i(x)=½k_i(x−x_i^0)^2를 적용받으며, 관측된 확률 p_i(x)는 Boltzmann 분포 p_i(x)∝exp