단백질 양과 전사 영향력 재평가: 대규모 데이터의 함정과 새로운 교정

초록

대규모 단백질 정량 연구가 실제 단백질 복제 수와 전사 기여도를 크게 낮게 추정한 것을, 61개의 하우스키핑 단백질을 기준으로 재보정함으로써 교정하였다. 교정 후 중위수 단백질은 세포당 약 17만 분자이며, mRNA 수준이 단백질 변동을 설명하는 비율은 최소 56%에서 81%까지 상승한다. 번역 효율의 실제 변동은 기존 추정보다 훨씬 작으며, 전사가 단백질 발현을 주도한다는 결론을 제시한다.

상세 분석

Schwanhausser 등(2011)의 전사·번역·단백질 분해 모델은 전체 단백질 정량을 질량분석 기반 ‘스펙트럼 카운트’로 추정했으며, 그 결과 평균 단백질 복제 수가 8 000~16 000 분자, mRNA와 단백질 간 상관계수가 R²≈0.4 수준이라고 보고했다. 그러나 이 접근법에는 두 가지 근본적인 오류가 존재한다. 첫째, 질량분석 감도와 펩타이드 검출 효율의 차이로 인해 저발현 단백질이 과소평가되고, 고발현 단백질은 포화 현상으로 실제 양보다 낮게 측정된다. 둘째, 정량 표준이 충분히 검증되지 않아 전체 데이터셋에 스케일링 오류가 전파된다.

본 연구는 61개의 하우스키핑 단백질(핵심 대사 효소, 구조 단백질 등)을 독립적인 정량법(예: SILAC, AQUA 펩타이드)으로 정확히 측정하고, 이를 기준으로 Schwanhausser 데이터 전체를 선형 보정하였다. 보정 결과, 전체 단백질 풀의 중위수 복제 수는 기존 추정치의 약 10배에 해당하는 1.7 × 10⁵ 분자로 증가하였다. 이와 동시에 mRNA와 보정된 단백질 양 사이의 상관계수는 R²≈0.56으로 상승했으며, 이는 전사 단계가 단백질 수준 변동을 절반 이상 설명한다는 강력한 증거다.

또한, 번역 효율을 직접 측정한 리보솜 프로파일링 데이터와 Schwanhausser가 추정한 번역 효율을 비교하였다. 두 데이터 간의 상관은 R²=0.13에 불과했으며, 실제 번역 효율 변동은 전체 변동의 약 12%에 불과했다. 이를 바탕으로 두 번째 독립적인 추정법을 적용하면, mRNA 수준이 단백질 변동을 설명하는 비율은 약 81%에 달한다.

오차 전파 분석에서는 mRNA 측정 오차(≈15%)와 단백질 측정 오차(≈30%)를 고려했으며, 아직 정량되지 않은 시스템적 편향(예: 샘플 손실, 펩타이드 선택성)으로 인해 실제 전사 기여도는 이보다 더 클 가능성이 있다. 최종적으로, 전사, RNA 분해, 번역, 단백질 분해 네 단계가 각각 차지하는 변동 기여도를 추정했을 때, 전사는 5060% 수준으로 가장 큰 비중을 차지하고, 번역은 1015% 수준으로 크게 축소된다.

이러한 결과는 기존 대규모 ‘단백질-전사’ 관계 연구가 과소평가한 전사의 역할을 재조명하고, 시스템 수준에서 정량 정확성을 확보하기 위한 교정 전략의 필요성을 강조한다. 특히, 절대적인 단백질 복제 수를 정확히 알 수 있는 표준화된 정량법이 없을 경우, 상대적 비교에만 의존하는 대규모 데이터는 생물학적 해석에 큰 오차를 초래할 수 있다. 향후 연구에서는 더 많은 하우스키핑 표준을 포함한 다중 교정 모델과, 리보솜 프로파일링 및 단백질 반감기 측정을 통합한 통합 모델을 구축함으로써 전사·번역·분해 네트워크의 정밀한 정량화를 목표로 해야 한다.