진화 경로 예측 가능성: 단백질 적합도 지형의 부드러움과 경로 발산

초록

오프‑격자 단백질 접힘 모델을 이용해 만든 적합도 지형을 전면 분석하였다. 적합도를 오작동에 대한 내성으로 정의하고, ‘평균 경로 발산(mean path divergence)’이라는 새로운 지표로 진화 경로의 결정성을 정량화했다. 전역적인 지형 거칠기(roughness) 지표와 평균 경로 발산은 강하게 상관했으며, 모델 및 실험 지형은 무작위 지형보다 훨씬 부드럽고, 소음이 가미된 가산(additive) 지형에 근접했다. 특히, 이러한 부드러운 지형은 피크(최적점)의 수가 현저히 적어, 단백질 진화가 물리적 접힘 제약에 의해 크게 제한된다는 결론을 제시한다.

상세 분석

이 연구는 단백질 진화의 예측 가능성을 정량적으로 평가하기 위해 두 가지 핵심 방법론을 도입한다. 첫째, 적합도를 ‘오작동된 복제 수’를 최소화하는 능력, 즉 미접힘에 대한 내성으로 정의한 오프‑격자 단백질 접힘 모델을 사용하였다. 이 모델은 아미노산 간 상호작용을 물리적으로 구현하면서도, 기존 실험 데이터가 보여준 진화 속도 분포와 일치하는 보편적 특성을 재현한다. 둘째, ‘평균 경로 발산(mean path divergence)’이라는 새로운 지표를 제시하였다. 이는 동일한 시작점과 종착점을 갖는 모든 단조(monotonic) 경로 쌍에 대해, 각 경로 간의 평균 해밍 거리의 가중 평균으로 정의된다. 경로가 서로 크게 다를수록 발산값이 커지고, 이는 진화 과정이 덜 결정적임을 의미한다.

지형의 전반적인 거칠기를 평가하기 위해 세 가지 보조 지표를 도입했다. (1) 지역 거칠기(local roughness)는 각 점과 이웃 점들 사이의 평균 제곱 차이로, 지형의 ‘가파름’과 ‘변동성’을 동시에 포착한다. (2) 피크 비율(peak fraction)은 주변에 더 높은 적합도가 없는 점들의 비율로, 순수 가산 지형에서는 1/N에 불과하지만 실제 모델에서는 현저히 낮았다. (3) 트리 컴포넌트(tree component)까지의 평균 거리로, 이는 단조 경로가 형성하는 트리 구조의 깊이를 측정한다.

실험적으로는 기존의 효소 진화 실험(β‑락타메이스, 세퀴터페인 신테이스 등)과 DNA‑단백질 결합 마이크로어레이 데이터를 활용한 실제 적합도 지형과 비교하였다. 모든 경우에서 모델 지형은 무작위 NK 지형보다 부드럽고, 가산 지형에 소량의 잡음(noise)을 추가한 형태와 거의 일치했다. 특히, 모델 지형은 동일한 거칠기 수준의 잡음이 가미된 가산 지형보다 피크가 현저히 적었다. 이는 단백질 접힘 물리학이 ‘서브옵티멀 피크’를 억제하고, 전체 적합도 지형을 전반적으로 평탄하게 만든다는 강력한 증거이다.

또한, 평균 경로 발산과 전역 거칠기 지표 사이의 강한 양의 상관관계를 발견했다. 부드러운 지형에서는 경로 간 차이가 작아 예측 가능성이 높으며, 거칠기가 큰 지형에서는 경로가 흩어져 예측이 어려워진다. 이는 기존 연구에서 제시된 ‘단조 경로의 수’만으로는 진화 경로의 예측 가능성을 충분히 설명하지 못한다는 점을 보완한다.

결과적으로, 단백질 진화는 물리적 접힘 제약에 의해 자연스럽게 ‘부드러운’ 적합도 지형을 형성하고, 이는 진화 경로의 제한성을 강화한다는 결론에 도달한다. 이러한 부드러움과 피크 결핍은 단백질이 기능을 유지하면서도 변이를 견디는 ‘robustness’를 확보하도록 진화했음을 시사한다.