네트워크 동역학에 큰 영향을 미치는 단백질은 더 천천히 진화하지만 필수성은 높지 않다

초록

이 연구는 단백질의 기능적 중요성을 정량화하기 위해 ‘동역학적 영향력(dynamical influence)’이라는 새로운 지표를 도입한다. 생화학적 모델을 이용해 각 반응 속도 상수가 네트워크 동역학에 미치는 영향을 계산하고, 해당 단백질이 관여한 모든 반응의 평균값을 그 단백질의 동역학적 영향력으로 정의한다. 12개의 생화학 시스템과 척추동물 유전체 데이터를 분석한 결과, 동역학적 영향력이 큰 단백질일수록 진화 속도가 느리다는 부정적 상관관계가 발견되었다. 그러나 이러한 단백질이 필수 유전자일 확률은 높지 않아, 필수성만으로는 기능적 중요성을 충분히 설명할 수 없음을 시사한다.

상세 분석

본 논문은 단백질 진화 속도와 기능적 중요성 사이의 관계를 보다 정밀하게 탐구하기 위해 기존의 ‘필수성(essentiality)’ 개념을 넘어서는 ‘동역학적 영향력(dynamical influence)’이라는 새로운 메트릭을 제시한다. 동역학적 영향력은 두 단계로 측정된다. 첫째, 각 생화학 반응에 대한 상세한 ODE 기반 모델을 구축하고, 파라미터(반응 속도 상수)마다 민감도 분석을 수행한다. 여기서 민감도는 특정 파라미터가 변했을 때 전체 네트워크의 시간‑진화(예: 농도 변화, 신호 전달)의 변동에 미치는 기여도를 정량화한다. 둘째, 특정 단백질이 촉매하거나 조절하는 모든 반응의 민감도 값을 평균하여 해당 단백질의 동역학적 영향력을 산출한다. 이 접근법은 단백질이 ‘존재’만으로도 필수적일 수 있지만, 그 효율성(속도 상수)의 미세한 변동이 실제 피트니스에 미치는 영향을 구분한다는 점에서 혁신적이다.

연구자는 12개의 서로 다른 생화학 네트워크(예: MAPK 경로, 포도당 대사, 세포주기 조절 등)를 선택하고, 각 네트워크에 대해 기존 문헌에서 검증된 모델 파라미터를 활용했다. 이후 vertebrate 종들의 orthologous 유전자 서열을 수집해 dN/dS 비율을 계산함으로써 각 단백질의 진화 속도를 추정하였다. 통계적 분석(피어슨 상관, 다중 회귀) 결과, 동역학적 영향력과 dN/dS 사이에 유의한 음의 상관관계(r ≈ –0.45, p < 0.01)가 나타났으며, 이는 높은 영향력을 가진 단백질일수록 아미노산 치환이 억제되는 경향을 의미한다.

흥미롭게도, 동역학적 영향력과 ‘필수성’(knock‑out 실험에서 생존 여부) 사이에는 뚜렷한 상관이 없었다. 즉, 필수적인 유전자는 반드시 높은 동역학적 영향력을 갖지 않으며, 반대로 높은 영향력을 가진 단백질도 비필수일 수 있다. 이는 ‘필수성’이라는 이분법적 지표가 단백질 기능의 미세 조정 수준을 포착하지 못한다는 점을 강조한다.

또한, 저자들은 동역학적 영향력이 단백질 발현 수준, 조직 특이성, 유전자 길이, 네트워크 내 연결도(노드 차수)와 같은 기존 알려진 요인들과 독립적인 영향을 미친다는 것을 다변량 회귀 분석을 통해 검증하였다. 이는 동역학적 영향력이 진화 속도에 미치는 효과가 다른 공변량에 의해 매개되지 않음을 시사한다.

이러한 결과는 진화생물학에서 ‘기능적 제약’이 어떻게 작용하는지를 재해석하게 만든다. 전통적으로는 ‘핵심 경로’에 속한 단백질이 느리게 진화한다는 가설이 있었지만, 본 연구는 ‘정량적 효율성’이 중요한 경우에만 강한 제약이 작용한다는 새로운 관점을 제공한다. 따라서, 단백질의 기능적 중요성을 평가할 때는 단순히 존재 여부가 아니라, 그 단백질이 네트워크 동역학에 미치는 정량적 기여를 고려해야 함을 강조한다.