Atomic Density Distributions in Proteins: Structural and Functional Implications

기술 요약: 단백질의 원자 밀도 분포: 구조적 및 기능적 함의

문제 제기
원자 패킹(Atomic packing)은 단백질 구조를 특징짓는 근본적인 지표로, 안정성, 진화율 및 기능적 역할에 영향을 미친다. 단백질 내부가 일반적으로 표면보다 더 조밀하게 패킹되어 있으며, 패킹 결함(공극, cavities)이 안정성을 저해할 수 있다는 점은 이미 확립되어 있다. 그러나 원자 밀도의 분포는 모든 단백질에서 균일하게 나타나지는 않는다. 기존 연구들은 주로 평균적인 통계적 특성이나 특정 구조적 지표(예: 접촉 수)에 의존하여 패킹을 분석해 왔다. 하지만 방대한 단백지 데이터셋을 통해 개별 밀도 프로파일을 직접 비교함으로써 체계적인 구조적 및 기능적 패턴을 식별하고자 하는 포괄적인 분석은 부족한 실정이었다. 본 연구는 평균값이 아닌, 대규모의 비중복 단백질 구조 집합에서 유도된 개별 밀도 분포를 직접 비교함으로써 단백질 밀도를 특징짓고자 하는 필요성을 다룬다.

방법론
저자들은 Protein Data Bank(PDB)에서 얻은 21,255개의 비중복 단백질 구조(80,000개 이상의 원자를 가진 구조를 제외하기 위해 필터링된 초기 22,478개의 집합에서 유래) 데이터셋을 분석하였다. 선택 기준에는 해상도 컷오프 2.2 Å, R-factor 컷오프 0.25, 그리고 서열 동일성 컷오프 ≤50%가 포함되었다.

핵심 방법론은 다음 단계들을 포함한다:

1. 시스템 준비: 모든 계산에는 생물학적 조립 단위(Biological assembly units)가 사용되었다. 누락된 수소 원자는 OpenBabel을 사용하여 추가되었으며, in vitro 조건을 모사하기 위해 사전 평형화된 물 분자들로 용매화(solvation)를 진행하되, 매몰된 공극(buried cavities) 내부에 물이 추가되는 것은 방지하였다.
2. 밀도 계산: 각 구조의 원자에 대해, 해당 원자를 중심으로 하는 구(sphere) 내의 국소 원자 밀도(Da/ų)를 계산하였다. 구 내에 포함된 모든 원자의 질량을 합산한 후 구의 부피로 나누었다. 이 과정은 단백질의 모든 원자에 대해 재귀적으로 수행되었다.
3. 파라미터 테스트: 임의의 반경 선택에 따른 편향을 완화하기 위해 세 가지 구 반경(5 Å, 6 Å, 7 Å)을 테스트하였다. 7 Å를 초과하는 반경은 과도한 평균화로 인해 필요한 해상도가 부족한 것으로 나타났다.
4. 데이터 처리: 각 단백질의 밀도 값은 100개의 빈(bin)을 가진 히스토그램으로 변환되었다. 이 분포들은 전체 원자 수에 의해 정규화되어 빈도 분포를 생성하였다.
5. 통계 분석:
   • 주성분 분석 (PCA): 데이터의 분포를 시각화하고 이상치를 식별하는 데 사용되었다.
   • 계층적 군집화 (Hierarchical Clustering): 모든 밀도 분포 쌍 사이의 유클리드 거리를 계산하여 거리 행렬을 생성하였다. 이 행렬들은 단백질들을 유사한 밀도 프로파일로 그룹화하기 위해 계층적 군집화("complete" 연결법 사용)에 적용되었다.
   • 상관관계 분석: 밀도 클러스터와 구조적 특징(크기, 2차 구조, B-factor, 물의 풍부도) 및 기능적 분류(EC 패밀리) 간의 관계를 조사하였다.

주요 결과

• 유의미한 변이성: 누적 밀도 분포를 통해 대부분의 단백질이 평균 주변에 군집을 이루지만, 통계적으로 유의미한 수의 구조들이 이상치로서 편차(일부 경우 >3σ)를 보임을 확인하였다.
• 분포의 군집화: 6 Å 반경 데이터(정보 함량의 최적 균형점으로 선택됨)의 계층적 군집화 결과 12개의 뚜렷한 클러스터가 식별되었다. 가장 뚜렷한 5개 클러스터(작은 그룹 제외)에 대한 추가 분석 결과, 밀도의 체계적인 변화가 드러났다:
  • 저밀도 클러스터: 구아니디노부티레이스(guanidinobutyrase)와 같은 단백질을 포함하였다.
  • 고밀도 클러스터: 사이토크롬(cytochromes) 및 코일드 코일(coiled-coil) 구조를 포함하였다.
• 구조적 상관관계:
  • 단백질 크기: 단백질 크기(잔기 수)와 중앙값 원자 밀도 사이에 음의 상관관계가 관찰되었다. 큰 단백질은 더 느슨하게 패킹되는 경향이 있으며 밀도 값의 범위가 좁은 반면, 작은 단백질은 더 큰 변동성을 보인다.
  • 안정성 지표: 높은 패킹 밀도는 정렬된 결정학적 물 분자의 풍부도와 양의 상관관계를 보였고, 원자 B-factor(온도 인자)와는 음의 상관관계를 보였다. 이는 조밀하게 패킹된 구조가 더 단단하고 안정적임을 시사한다.
  • 2차 구조: 밀도가 높은 클러스터는 α-헬릭스의 비율이 높고 β-스트랜드의 비율이 낮았다. 특히, "코일드 코일" 모티프는 기계적 안정성을 위한 "knobs-into-holes" 패킹 요구 조건과 일치하게 가장 밀도가 높은 클러스터에서 유의미하게 과표집되었다.
• 기능적 상관관계:
  • 사이토크롬: 이 전자 전달 단백질들은 가장 밀도가 높은 클러스터에서 유의미하게 과표집되었으며, 이는 조밀한 패킹이 효율적인 전자 전달을 촉진한다는 가설을 뒷받침한다.
  • 효소 패밀리: 가수분해효소(Hydrolases)는 밀도가 높은 클러스터를 선호한 반데, 전이효소(transferases)와 연결효소(ligases)는 더 느슨한 클러스터를 선호했다. 저자들은 이것이 "안정성-활성" 트레이드오프 및 세포 환경에 존재하는 물(substrate로서의 물)을 다루는 메커니즘적 차이와 관련이 있을 수 있다고 제안한다.

의의 및 주장
저자들은 개별 밀도 프로파일을 비교하는 방식이 단백질의 구조적 및 기능적 패턴을 정량적으로 특징짓는 강력한 방법임을 입증한다고 주장한다. 본 연구는 패킹, 안정성 및 단백질 크기 사이의 관계에 관한 기존의 발견들을 검증하고 확장한다.

주요 주장은 다음과 같다:

1. 방법론적 유용성: 계산된 원자 밀도 분포는 단백질 구조의 패킹 수준을 분류하고 독특한 구조적 특성을 가진 이상치를 식별할 수 있는 품질 지표 역할을 한다.
2. 구조-기능 연결: 특정 구조적 모티프(코일드 코일)와 기능적 패밀리(사이토크롬, 가수분해효소)는 뚜렷한 패킹 선호도를 보인다. 조밀한 패킹은 기계적 안정성과 전자 전달을 위한 기능적 요구사항처럼 보이는 반면, 느슨한 패킹은 특정 효소 메커니즘(예: 전이효소)과 연관될 수 있다.
3. 진화적 관점: 관찰된 원자 밀도의 변이성은 특정 밀도에 대한 능동적인 현재의 선택 결과라기보다, 분자 진화 과정의 잔재로 해석된다. 저자들은 안정성과 새로운 효소 기능 사이의 트레이드오프가 서로 다른 단백질 패밀리에 걸쳐 관찰되는 패킹 밀도의 분포를 유도할 수 있다고 제안한다.

본 논문은 원자 밀도에 체계적인 패턴이 존재하지만, 그것이 구조적 또는 기능적 특성과 고정된 일대일 관계로 대응하지는 않으며, 이는 진화적 역사와 기능적 제약 사이의 복잡한 상호작용을 반영한다고 결론짓는다.