[논문 리뷰] Combination of genetic crossover and replica-exchange method for conformational search of protein systems

Yoshitake Sakae, Tomoyuki Hiroyasu외 3명|arXiv (Cornell University)|2015. 05. 21.

Protein Structure and Dynamics참고 문헌 20인용 수 4

한 줄 요약

이 논문은 단백질 접힘 시뮬레이션에서 샘플링 효율성을 향상시키기 위해 유전적 교배와 복제체 교환 분자 동역학(REMD)을 조합한 하이브리드 구형 샘플링 방법을 제안한다. 단백질 구조의 전역적 교환을 가능하게 하고 온도 복제체 교환과 통합함으로써, 20개 아미노산으로 이루어진 Trp-cage 마이크로단백질의 구형 영역을 효율적으로 탐색하며, 실험 데이터와 근접한 구조를 도출하고 전통적인 MD 및 REMD에 비해 국소 최소값에 갇히는 것을 피하는 데 성공한다.

ABSTRACT

We combined the genetic crossover, which is one of the operations of genetic algorithm, and replica-exchange method in parallel molecular dynamics simulations. The genetic crossover and replica-exchange method can search the global conformational space by exchanging the corresponding parts between a pair of conformations of a protein. In this study, we applied this method to an $α$-helical protein, Trp-cage mini protein, which has 20 amino-acid residues. The conformations obtained from the simulations are in good agreement with the experimental results.

연구 동기 및 목표

• 단백질의 구형 탐색 과정에서 국소 에너지 최소값에서 벗어나지 못하는 전통적인 분자 동역학 및 몬테카를로 시뮬레이션의 한계를 극복하기 위해.
• 유전적 교배를 통한 전역적 구조 업데이트를 복제체 교환 방법과 융합하여 생체분자 시뮬레이션의 샘플링 효율성을 향상시키기 위해.
• 기준 시스템인 Trp-cage 마이크로단백질에 대해 제안된 방법의 정확성과 수렴성을 평가하기 위해.
• 이 방법이 더 큰 단백질에 적용 가능한 확장성과 병렬 처리 가능성 여부를 입증하기 위해.

제안 방법

• 고정된 온도에서 M개의 초기 구조를 동시에 진화시키는 병렬 분자 동역학(PMD/GAc)을 사용하며, 이 과정에서 유전적 교배를 적용한다.
• 두 개의 구조에서 무작위로 n개의 연속된 잔류물을 선택하고, 이들의 이황각을 교환한 후 에너지 최소화 및 배경 구조 고정 조건을 통해 평형화한다.
• 에너지 차이에 기반하여 부모와 자식 구조 중 어느 쪽을 선택할지 결정하기 위해 메트로폴리스 기준을 사용한다: w(p→c) = min(1, exp{−β[E_c − E_p]}).
• 이 방법은 복제체 교환 분자 동역학(REMD/GAc)로 확장되며, 서로 다른 온도에서 실행되는 여러 복제체가 구성 구조를 교환함으로써 구형 샘플링을 향상시킨다.
• 각 복제체에서 별도로 교배 작동을 수행하며, 교배 사건의 수용 비율을 모니터링하여 효율성을 평가한다.
• 시뮬레이션 프로토콜은 선택 이전에 새로운 구조를 안정화시키기 위해 교환 영역에 대해 배경 이황각 고정 조건을 포함한 짧은 평형 단계를 포함한다.

실험 결과

연구 질문

• RQ1 유전적 교배를 복제체 교환 분자 동역학와 융합함으로써 단백질 접힘 시뮬레이션에서의 구형 샘플링이 향상되는가?
• RQ2 PMD/GAc 및 REMD/GAc 방법은 전통적인 MD 및 REMD에 비해 샘플링 효율성과 천연 상태에의 수렴성에서 어떻게 비교되는가?
• RQ3 유전적 교배 작동이 국소 에너지 최소값에 갇히는 정도를 어느 정도 감소시키는가?
• RQ4 교배 작동이 RMSD 값의 분포와 터널링 사건의 수에 어떤 영향을 미치는가?
• RQ5 이 방법은 Trp-cage 마이크로단백질에 대해 실험적 구조 데이터를 얼마나 잘 재현하는가?

주요 결과

• PMD/GAc에서의 평균 RMSD(5.50 Å)는 전통적인 MD(7.06 Å)보다 유의미하게 낮아 천연 구조 근처의 샘플링 향상을 시사한다.
• REMD/GAc 시뮬레이션은 PDB 구조 기준으로 RMSD가 1.93 Å로 전통적인 REMD(2.03 Å)보다 略적으로 우수했으며, 천연 유사 구조 기준으로는 1.27 Å로 더 높은 정확도를 보였다.
• REMD/GAc 시뮬레이션은 RMSD 분포에서 6.68 Å에 단일 주요 피크를 보였고, 전통적인 REMD는 3.26 Å에 보조 피크를 보이며 국소 최소값에 갇히는 현상이 줄어든 것으로 나타났다.
• 터널링 사건 수는 REMD/GAc에서 107건, 전통적인 REMD에서는 54건으로, 교배 작동 덕분에 더 효율적인 구형 샘플링이 이루어졌음을 보여준다.
• 교배 작동의 평균 수용 비율은 0.26로, 시스템과 교배 길이에 따라 중간 수준이지만 효과적인 교환 빈도를 나타낸다.
• 이 방법은 실험적 구조에 가까운 구조를 성공적으로 생성하여, 작은 단백질에서의 구형 탐색에 효과적임을 검증하였다.