단백질 변이 표에 숨겨진 열역학 정보

초록

본 논문은 1992년 제시된 20×20 대체 행렬 BLOSUM62와 1982년 KD, 2007년 MZ 두 가지 친수성·소수성 척도를 결합하여 변이율을 분석한다. KD는 1차 상전이, MZ는 2차 상전이 특성을 갖으며, 변이율은 2×10 이진 성질을 4×5 하위 그룹(뮤톤)으로 정량화한다. MZ 스케일이 KD보다 일치도가 높다.

상세 분석

이 연구는 기존의 단백질 서열 비교에 널리 사용되는 BLOSUM62 행렬을 두 개의 물리화학적 친수성·소수성 척도와 결합함으로써 변이 패턴에 내재된 열역학적 정보를 추출한다. KD 척도는 1952년 Kyte‑Doolittle가 제안한 20개 아미노산의 평균 전이 자유에너지(ΔG) 기반으로, 1차 상전이와 유사한 급격한 변화를 나타낸다. 반면 MZ 척도는 Moret‑Zebende가 2007년에 제시한, 아미노산의 표면 접근성 및 구조적 상관성을 반영한 2차 상전이 특성을 지닌 스케일이다. 두 스케일을 BLOSUM62와 정규화하여 각 아미노산 쌍의 변이 확률을 재구성하면, 변이율이 단순히 개별 친소성(친수성·소수성) 값의 차이만으로 설명되지 않음을 확인한다. 대신, 20개의 아미노산을 2×10 이진 특성(친소성 vs 비친소성)으로 분류하고, 이를 다시 4×5 형태의 사분면(쿼터너리) 그룹으로 재조합한다. 이러한 그룹을 ‘뮤톤(Mu5)’이라 명명하고, 각 뮤톤 내에서 변이가 집중되는 경향을 발견하였다. 특히 MZ 기반 뮤톤은 KD 기반 뮤톤에 비해 BLOSUM62와의 상관계수가 현저히 높으며, 변이 패턴이 더 정교하게 재현된다. 이는 MZ 스케일이 아미노산 간 상호작용의 미세한 구조적 변화를 포착함을 의미한다. 또한, 5×4 형태의 5차원 그룹(퀸터리)으로 재구성할 경우, 변이율과의 일치도가 급격히 감소한다는 점에서, 단순히 더 세분화된 그룹화가 오히려 열역학적 의미를 손실한다는 결론을 도출한다. 이러한 결과는 단백질 진화와 설계에서 변이율을 예측할 때, 1차 상전이 기반 KD보다 2차 상전이 기반 MZ를 활용하는 것이 더 신뢰할 만함을 시사한다.