From populations to absolute binding affinities in molecular simulations: exact volumetric terms and practical estimators

기술 요약: 분자 시뮬레이션에서의 인구수로부터 절대 결합 친화도 도출

문제 제기
분자 시뮬레이션으로부터 정확한 절대 결합 자유 에너지($\Delta F^0$) 및 평형 결합 상수($K$)를 계산하는 것은 계산 생물물리학 및 약물 설계 분야의 핵심적인 과제로 남아 있다. 상대적 결합 자유 에너지에 대한 방법론은 잘 확립되어 있는 반면, 절대적 계산은 비결합 상태(unbound state)의 부피 기여를 엄밀하게 다루어야 한다. 기존의 접근 방식들은 이 부피 항을 부가적인 보정치로 취급하거나, '결합된(bound)' 상태와 '비결립된(unbound)' 상태의 정의에 대한 암묵적인 가정을 바탕으로 한 방법론 특유의 유도 방식에 의존하곤 한다. 결과적으로, 실무자들은 어떤 부피가 표준 상태 참조값에 들어가는지, 결합 상태의 경계는 어디인지, 그리고 유한한 샘플링에 적합한 통계적 추정량은 무엇인지와 같은 미해결 문제에 직면하게 된다. 이러한 모호함은 복잡한 열역학적 상태의 전체 통계적 가중치를 포착하지 못하는 단순화된 "단일 빈(single-bin)" 추정량을 사용할 때 체계적인 오차를 초래할 수 있다.

방법론
저자들은 정준 앙상블 분배 함수(canonical partition function)로부터 시작하여, 결합 상수 $K$를 결합 상태($P_b$)와 비결합 상태($P_u$)의 상대적 통계적 인구수의 직접적인 관계로 표현하는 일반적인 통계 역학적 프레임워크를 도출한다. 핵심 관계식은 다음과 같다:
$$K = \frac{P_b}{P_u/V}$$
여기서 $V$는 시스템의 부피이다. 이 공식은 비결합 상태의 부피 기여를 상태 간의 자유 에너지 차이로부터 명시적으로 분리한다.

본 프레임워크는 생체 분자 시뮬레이션에서 흔히 나타나는 세 가지 기하학적 대칭성을 다룬다:

1. 구형 대칭(Spherical symmetry): 정의된 포켓이 없는 용액 내 단백질에 리간드가 결합하는 경우.
2. 원형 대칭(Circular symmetry): 지질 막 내에서의 이량체화(dimerization).
3. 원통형 대칭(Cylindrical symmetry): 매립된 포켓이나 막 삽입 단백질에 리간드가 결합하는 경우로, 반응 좌표가 특정 축을 따른 변위인 경우.

저자들은 구속되지 않은(unrestrained) 시뮬레이션과 구속된(restrained) 시뮬레이션 모두를 위한 실용적인 추정량을 도출한다. (수렴을 가속화하기 위해 포텐셜로 리간드를 부피 $V^R$ 내에 가두는) 구속된 시뮬레이션의 경우, 비결합 상태가 구속 영역 내에서 균질하다는 물리적으로 완만한 가정하에 정확한 결과를 회복하기 위한 해석적 보정치를 제공한다. 결정적으로, 이 유도는 다음을 구분한다:

• 확장 영역 추정량(Extended-region estimators): 정의된 비결합 구간에 대해 확률 밀도를 적분함.
• 하이브리드 추정량(Hybrid estimators): 결합 상태를 적분하면서 비결합 상태에는 단일 빈 참조를 사용함.
• 단일 빈 추정량(Single-bin estimators): 두 상태 모두에 단일 빈을 사용함 (구조화된 결합 상태의 경우 신뢰도가 낮은 것으로 식별됨).

프레임워크는 유한한 샘플링으로 인한 노이즈가 있는 히스토그램을 처리하기 위해 커널 밀도 추정(kernel density estimation)을 포함하며, 차원적 일관성을 보장하기 위해 부피 항에 대한 빈 너비($\Delta x$)에 대한 명시적인 보정을 제공한다.

주요 기여

1. 정확한 유도: 본 논문은 모든 항, 특히 비결합 상태의 부피 항의 기원을 명확히 드러내는, $K$에 대한 단일하고 자기 완결적인 유도 과정을 제공한다. 이는 표준 상태 참조 부피가 오직 비결합 상태에만 접근 가능한 부피임을 명확히 한다.
2. 명시적 상태 정의: 이 연구는 결합 상태의 정의를 숨겨진 가정에서 통제 가능한 입력값으로 변환한다. 결합 상수가 결합된 상태의 통계적 가중치의 합(포즈 또는 결합 모드)에 포함되는 어떤 열역학적 하위 배스인(sub-basin)들을 포함하느냐에 따라 달라짐을 입증한다.
3. 구속에 대한 보정: 저자들은 기하학적 구속을 사용하는 시뮬레이션에 대한 폐쇄형 해석적 보정치를 유도하여, 비결합 상태가 균질하다는 조건하에 조정 가능한 파라미터 없이 정확한 결합 친화도를 회복하는 방법을 보여준다.
4. 체계적 오차 식별: 본 연구는 흔히 사용되는 단일 빈 추정량이 이론적으로 올바른 확장 영역 또는 하이브리드 추정량과 크게 달라질 수 있음을 식별하며, 특히 결합 상태가 구조적 복잡성(다중 포즈 또는 하위 배스)을 가질 때 그러하다.

결과
프레임워크는 세 가지 단계적으로 복잡해지는 시스템에 적용되었다:

• Lennard-Jones 이량체: 단순한 테스트 케이스로서, 부피 보정이 제한된 부피에 대한 원시 자유 에너지 차이의 의존성을 정확히 설명하여 서로 다른 구속 반경에서도 일관된 $\Delta F^0$를 산출함을 확인하였다.
• Cucurbit[7]uril/1-Adamantanol: 벤치마크 호스트-게스트 시스템이다. 결과에 따르면 확장 영역 및 하이브리드 추정량은 서로 일치하는 결합 친화도($\approx -17.6$ kcal/mol)를 나타냈다. 그러나 단일 빈 추정량은 약 1 kcal/mol 정도 편차를 보였는데, 이는 결합 상태(hydroxyl 그룹의 두 가지 뚜렷한 방향 포함)의 전체 통계적 가중치를 포착하는 데 단일 빈이 실패했기 때문이다.
• Galactonate–DgoT 리간드-단백질 복합체: 구조화된 결합 상태를 가진 수송체 시스템이다. 확장 영역 및 하이브리드 추정량은 일치하였으나($\Delta F^0 \approx -1.1$ kcal/mol), 단일 빈 추정량은 약 1 kcal/mol의 편차를 보였다. 또한, 본 연구는 '부분 접근(partial-access)' 설정(구속이 이차 하위 배스에 영향을 줄 수 있는 상황)에서 결합 상태 경계의 배치가 친화도를 약 $-0.6$ kcal/mol만큼 변화시킬 수 있음을 보여주었다.

의의
본 논문은 자신의 유도가 절대 결합 친화도를 계산하기 위한 "신뢰할 수 있는 토대"와 "물리적으로 직관적인 근거"를 제공한다고 주장한다. 모든 근사치와 상태 정의를 명시함으로써, 본 프레임워크를 통해 실무자들은 다음과 같은 작업을 수행할 수 있다:

• 이전에 다뤄지지 않은 설정에 대한 새로운 추정량 구축.
• 공통된 참조점을 기준으로 기존 결과의 타당성 검증.
• 시뮬레이션 정의를 실험적 해상도(예: 벌크 평형 측정과 부위별 NMR 데이터의 구분)와 일치시킴.

저자들은 올바른 추정량과 잘못된 추정량 사이에서 관찰된 약 1 kcal/mol의 차이가 실험과의 의미 있는 비교를 위해 "결정적(decisive)"이라고 강조하는데, 이는 절대 결합 자유 에너지 계산에서 통상적으로 주장되는 정확도를 상회하는 수치이기 때문이다. 본 연구는 새로운 실험 프로토콜을 제안하는 것이 아니라, 계산적 예측이 모델링하고자 하는 물리적 상태의 정의와 일치하도록 보장하는 엄격한 이론적 경로를 제공한다.