[논문 리뷰] Simulating Protein Conformations through Global Optimization

Antonio Mucherino, Onur Şeref외 1명|ArXiv.org|2008. 11. 19.

Protein Structure and Dynamics참고 문헌 23인용 수 3

한 줄 요약

이 논문은 단순한 구조적 제약 조건—밀도, 자체 교차 방지, 알파-헤르름 형성—기반으로 전역 최적화 문제로 공식화된 기하 모델을 제안한다. 원숭이 탐색 메타휴리스틱을 사용하여 높은 품질의 구조를 효율적으로 생성하며, 1v66 및 2cro와 같은 실제 단백질과 비교해 낮은 RMSD(최저 4.97Å)를 기록함으로써, 전체 물리 기반 힘장이 필요 없이도 단순한 기하 규칙만으로 생물학적으로 타당한 구조를 도출할 수 있음을 보여준다.

ABSTRACT

Many researches have been working on the protein folding problem from more than half century. Protein folding is indeed one of the major unsolved problems in science. In this work, we discuss a model for the simulation of protein conformations. This simple model is based on the idea of imposing few geometric requirements on chains of atoms representing the backbone of a protein conformation. The model leads to the formulation of a global optimization problem, whose solutions correspond to conformations satisfying the desired requirements. The global optimization problem is solved by the recently proposed Monkey Search algorithm. The simplicity of the optimization problem and the effectiveness of the used meta-heuristic search allowed the simulation of a large set of high-quality conformations. We show that, even though only few geometric requirements are imposed, some of the simulated conformation results to be similar (in terms of RMSD) to conformations real proteins actually have in nature.

연구 동기 및 목표

• 복잡한 물리 기반 힘장에 의존하지 않고도 단순하고 기하 기반의 모델을 개발하여 단백질 구조를 시뮬레이션하는 것.
• 구조 시뮬레이션을 핵심 기하 조건에 의해 제약된 전역 최적화 문제로 공식화하는 것.
• 이러한 최소 모델이 천연 단백질의 구조와 유사한 구조를 생성할 수 있는지 평가하는 것.
• 원숭이 탐색 알고리즘이 단백질의 구조적 다양성 공간을 탐색하는 데 있어 효율성과 효과성을 입증하는 것.
• 기하 제약 조건이 전통적인 에너지 기반 모델의 타당한 대안이 될 수 있는지 탐색하는 것.

제안 방법

• 모델은 세 가지 핵심 기하 조건을 적용한다: 밀도, 자기 교차 방지, 알파-헤르름 세그먼트 형성.
• 각 조건은 페널티 함수로 변환되며, 총 목적 함수는 이들 함수의 가중치 합으로 최소화된다.
• 전역 최적화 문제는 복잡하고 다모드적인 검색 공간에서 효율적인 것으로 알려진 원숭이 탐색 메타휴리스틱을 사용하여 해결한다.
• 기하 제약 조건의 가중치를 변화시켜 다양한 유효한 구조를 생성함으로써 구조적 다양성 공간을 탐색한다.
• 고정된 길이와 이차 구조 조성비를 가진 단백질 사슬을 시뮬레이션하며, 주로 전부 알파-헤르름 단백질 유사 접힘에 집중한다.
• 실제 단백질과의 구조 유사도는 Cα 원자 정렬 후 루트 평균 제곱 오차(RMSD)를 사용하여 정량화한다.

실험 결과

연구 질문

• RQ1 소수의 구조적 제약 조건에 기반한 최소 기하 모델이 자연에서 관찰되는 단백질 구조와 유사한 구조를 생성할 수 있는가?
• RQ2 이러한 기하 제약 조건 하에서 원숭이 탐색 알고리즘이 단백질의 구조적 다양성 공간을 얼마나 효율적으로 탐색할 수 있는가?
• RQ3 기하 규칙만을 사용할 경우, 시뮬레이션된 구조와 실제 단백질 구조 간의 RMSD 값은 어떻게 비교되는가?
• RQ4 세부 에너지 함수가 없더라도 기하 최적화를 통해 생물학적으로 관련성이 있는 단백질 구조를 도출할 수 있는가?
• RQ5 기하 제약 조건의 수를 늘일수록 시뮬레이션된 구조가 천연 유사 접힘으로 수렴하는가?

주요 결과

• 모델은 단순한 기하 제약 조건만을 사용하여 합리적인 시간 내에 고품질의 단백질 구조를 다수 생성하였다.
• 한 시뮬레이션 구조는 실제 단백질 1v66와 RMSD 4.97Å를 기록하여 강력한 구조 유사성을 보였다.
• 다른 구조는 2cro 단백질과 RMSD 5.36Å를 기록하여 결과의 생물학적 관련성을 추가로 확인하였다.
• 모델의 단순성에도 불구하고, 시뮬레이션된 구조는 밀도가 높고 자기 교차가 없으며 안정적인 알파-헤르름 세그먼트를 포함하고 있었다.
• 결과는 기하 성질만으로도 세부 에너지 함수 없이도 단백질 유사 접힘의 형성이 이끌릴 수 있음을 시사한다.
• 원숭이 탐색 알고리즘이 전역 최적화 문제를 효과적으로 해결하여 구조적 다양성 공간의 효율적 탐색이 가능함을 입증하였다.