[논문 리뷰] Rejection-free kinetic Monte Carlo simulation of multivalent biomolecular interactions

Jin Yang, William S. Hlavacek|arXiv (Cornell University)|2008. 12. 25.

Gene Regulatory Network Analysis참고 문헌 35인용 수 3

한 줄 요약

이 논문은 다가액 생분자 상호작용을 시뮬레이션하기 위한 거부 없는 운동 몬테카를로 방법을 제시한다. 이 방법은 반응 확률을 결합체 상태에 기반해 직접 샘플링함으로써 기존의 비효율적인 시행착오적 거부 단계를 피한다. 이는 집합체의 수에 비례하는 스케일링을 통해 기존의 거부 기반 접근보다 더 뛰어난 효율성과 확장성을 보이며, 특히 광범위한 집합체 형성이 일어나는 시스템에서 뛰어난 성능을 발휘한다.

ABSTRACT

The system-level dynamics of multivalent biomolecular interactions can be simulated using a rule-based kinetic Monte Carlo method in which a rejection sampling strategy is used to generate reaction events. This method becomes inefficient when simulating aggregation processes with large biomolecular complexes. Here, we present a rejection-free method for determining the kinetics of multivalent biomolecular interactions, and we apply the method to simulate simple models for ligand-receptor interactions. Simulation results show that performance of the rejection-free method is equal to or better than that of the rejection method over wide parameter ranges, and the rejection-free method is more efficient for simulating systems in which aggregation is extensive. The rejection-free method reported here should be useful for simulating a variety of systems in which multisite molecular interactions yield large molecular aggregates.

연구 동기 및 목표

• 다가액 리간드-수용체 상호작용의 운동 몬테카를로 시뮬레이션에서 거부 샘플링의 비효율성을 해결하기 위해, 특히 큰 집합체가 지배하는 솔-젤 영역에서의 문제를 다루는 것.
• 모든 시도하는 위치를 거부하지 않고도 정확한 반응 확률을 계산함으로써 통계적으로 정확한 시뮬레이션 방법을 개발하는 것.
• 모든 개별 분자 위치가 아닌 집합체만 추적함으로써 광범위한 집합체 형성 시스템에서 계산 비용을 줄이는 것.
• 연속체 해법과의 비교를 통해 방법의 타당성을 검증하고, 다양한 운동 파rameter에서 거부 기반 방법과의 성능을 비교하는 것.
• 다중위치, 다가액 상호작용 및 큰 분자 집합체를 포함하는 복잡한 생화학적 시스템의 효율적이고 확장 가능한 시뮬레이션을 가능하게 하는 것.

제안 방법

• 규칙 목록에서 반응 규칙를 샘플링하고, 각 위치의 상태와 집합체 내 연결성에 기반해 그 위치가 반응에 참여할 정확한 확률을 계산한다.
• 무작위로 위치를 샘플링하고 불일치하는 경우를 거부하는 대신, 알고리즘은 반응 확률 분포에 따라 직접적으로 위치를 선택함으로써 거부 이벤트를 완전히 제거한다.
• 동적으로 갱신되는 집합체 목록을 유지하고, 각 집합체 내의 위치 수와 연결성을 기반으로 반응 속도를 계산함으로써 전체 위치에서의 검색 공간을 활성 집합체로만 제한한다.
• 모든 집합체의 누적 반응 속도를 기반으로 반응 이벤트를 선택함으로써 통계적 정확성과 세부 균형을 보장한다.
• 전체 네트워크의 전수 조사 없이도 다음 반응 이벤트를 효율적으로 식별하기 위해 우선순위 큐 또는 유사한 자료구조를 사용한다.
• 비국소적 적용 조건(예: 사이클릭 복합체 내에서의 집합체 내 결합 금지)을 정확히 계산하기 위해 연결성 정보를 추적한다.

실험 결과

연구 질문

• RQ1 큰 집합체가 지배하는 시스템에서 다가액 생분자 상호작용의 운동 몬테카를로 시뮬레이션을 어떻게 더 효율적으로 만들 수 있는가?
• RQ2 거부 없는 접근 방식이 전통적인 거부 샘플링에서 발생하는 높은 거부 비율로 인한 성능 저하 문제를 해결할 수 있는가?
• RQ3 거부 없는 방법의 스케일링 행동은 집합체 수에 비해 분자 위치 수에 비해 어떻게 달라지는가?
• RQ4 다양한 이탈률과 결합 가용성에 따라 거부 없는 방법과 거부 기반 방법의 성능은 어떻게 비교되는가?
• RQ5 거부 없는 방법이 저차원의 다가액 결합 모델에 대해 연속체 해법을 정확하게 재현할 수 있는가?

주요 결과

• 거부 없는 방법은 광범위한 운동 파rameter 범위에서 기존의 거부 기반 방법과 동일하거나 더 뛰어난 성능을 발휘하며, 특히 광범위한 집합체 형성이 일어나는 시스템에서 뛰어난 성능을 보인다.
• 이 방법의 계산 비용은 총 분자 위치 수가 아닌 집합체 수에 비례하므로, 솔-젤 영역에서 상당한 효율성 향상을 이룬다.
• 평형 상태에서 평균 집합체 수 ⟨NA⟩는 수용체 수 NR=5000보다 4분의 1 이하였으며, 이는 집합체 목록이 작고 관리 가능함을 시사한다.
• 거부 비율 θ가 1에 가까워지는 시스템, 즉 집합체 형성이 높은 경향이 있는 시스템에서는 거부 없는 방법이 거부 방법을 능가한다.
• 반응 속도와 비국소적 연결성 조건(예: 집합체 내 결합 금지)을 정확히 고려함으로써 통계적 정확성이 유지된다.
• 큰 집합체에 대한 반복적인 그래프 탐색을 피함으로써 표준 방법에서의 주요 성능 저하 요소를 줄일 수 있어 알고리즘의 효율성이 더욱 향상된다.