[논문 리뷰] Understanding hierarchical protein evolution from first principles

Nikolay V. Dokholyan, Eugene I. Shakhnovich|arXiv (Cornell University)|2001. 04. 24.

Protein Structure and Dynamics인용 수 6

한 줄 요약

이 논문은 단백질의 천연 상태 안정성에 기반한 동적 Z-스코어 모델을 제안하여 계층적인 단백질 진화를 설명한다. 열역학적 선택 조건 하에서의 시퀀스 진화 시뮬레이션을 통해 관측된 보존 패턴을 재현하고, 진화적 시간 스케일의 분리가 드러남: 단기적으로는 보존된 핵심 잔기들을 가진 동원가족이 형성되고, 장기적으로는 낮은 시퀀스 유사도이지만 유사한 접힘 구조를 가진 유사가족이 나타나며, 이는 생물학적 자료와 일치한다.

ABSTRACT

We propose a model that explains the hierarchical organization of proteins in fold families. The model, which is based on the evolutionary selection of proteins by their native state stability, reproduces patterns of amino acids conserved across protein families. Due to its dynamic nature, the model sheds light on the evolutionary time scales. By studying the relaxation of the correlation function between consecutive mutations at a given position in proteins, we observe separation of the evolutionary time scales: at the short time intervals families of proteins with similar sequences and structures are formed, while at long time intervals the families of structurally similar proteins that have low sequence similarity are formed. We discuss the evolutionary implications of our model. We provide a ``profile'' solution to our model and find agreement between predicted patterns of conserved amino acids and those actually observed in nature.

연구 동기 및 목표

• 자연에서 관측된 단백질 접힘 가족의 계층적 조직을 설명하기 위해.
• 진화적 시간 스케일 전반에 걸쳐 시퀀스 보존과 구조 유사성이 어떻게 공존하는지 이해하기 위해.
• 단백질 안정성을 주요 선택 압력으로 삼아 제1원리에서 진화 과정을 모델링하기 위해.
• 실제 단백질 가족에서 관측된 아미노산 보존 패턴을 재현하기 위해.
• 단백질 시퀀스 분화에서 진화적 시간 스케일의 분리 현상을 조사하기 위해.

제안 방법

• Z-스코어 모델은 천연 상태를 안정화하기 위해 Z-스코어를 최소화하며, 접힘 안정성에 대한 진화적 선택을 시뮬레이션한다.
• 몬테카를로 시뮬레이션을 사용하여 시퀀스 공간을 탐색하고, 안정성 향상에 기여하는 돌연변이만 수용한다.
• 모델를 해석적으로 해결하고 각 위치에서의 잔기 보존을 예측하기 위해 프로파일 근사법을 유도한다.
• 예측된 보존과 관측된 보존을 비교하기 위해 정렬된 단백질 가족에서 시퀀스 엔트로피를 계산한다.
• 에너지 표면에서의 구조적 데코이를 모델링하기 위해 동중합체 시뮬레이션에서 접촉 빈도를 추정한다.
• 모델는 시간 진화를 포함하여 연속된 돌연변이 간 상관 함수 분석을 통해 시간 스케일을 추론할 수 있다.

실험 결과

연구 질문

• RQ1 천연 상태 안정성에 대한 진화적 선택이 단백질 접힘 가족의 계층적 조직을 어떻게 유도하는가?
• RQ2 단백질 시퀀스 진화에서 진화적 시간 스케일의 분리는 무엇에 의해 발생하는가?
• RQ3 고도로 유사한 구조를 가진 동원가족 내에서 왜 오직 3–4%의 아미노산만 보존되는가?
• RQ4 안정성에 기반한 동적 모델이 실제 단백질 가족에서 관측된 아미노산 보존 패턴을 재현할 수 있는가?
• RQ5 장기적 시간 스케일에서 상관관계가 있는 돌연변이가 어떻게 시퀀스가 다를 뿐만 아니라 구조적으로도 유사한 유사체를 유도하는가?

주요 결과

• 모델은 단백질 가족 내에서 관측된 아미노산 보존 패턴을 성공적으로 재현하였으며, 예측된 보존 프로파일과 실제 프로파일 간에 양호한 일치를 보였다.
• 명확한 진화적 시간 스케일의 분리가 관찰됨: 단기적 진화는 보존된 핵심 잔기를 가진 동원가족을 형성하고, 장기적 진화는 시퀀스가 다를지라도 구조적으로 유사한 유사가족을 생성한다.
• 모델은 동원가족 내에서 오직 3–4%의 잔기만 보존된다고 예측하며, 생물학적 관측과 일치한다.
• 모델의 프로파일 해법은 실제 단백질 가족의 엔트로피를 정확히 캡처하여 그 예측 능력을 검증한다.
• Z-스코어 최소화 메커니즘은 안정적이고 빠르게 접히는 단백질을 효과적으로 선택하며, 진화적 압력과 일치한다.
• 모델은 功能적으로 중요한 위치에서의 돌연변이가 흔치 않으며, 보존적인 성격을 띠어 천연 구조와 기능을 유지함을 설명한다.