단백질 구조의 포장 단위 탐색: H‑폼 기반 보존 접촉 분석

초록

본 논문은 서로 다른 2차 구조 요소(헬릭스 또는 베타‑시트) 사이에서 Cα 원자가 8.5 Å 이내로 접촉하는 경우, 각 접촉 중심을 기준으로 폭 5개의 짧은 백본 조각을 추출한 “H‑폼”을 정의한다. H‑폼의 거리·각도 등 3차원 기하학적 특성을 정량화하고, 구조 정렬 도구를 이용해 여러 단백질 패밀리와 슈퍼패밀리에서 보존되는 H‑폼을 탐색한다. 보존된 H‑폼은 서열·구조적 일치를 동시에 보이며, 이러한 포장 단위가 폴딩 코어와 구조 예측에 유용한 제약 조건이 될 수 있음을 제시한다.

상세 분석

논문은 먼저 단백질 접힘 과정에서 소수의 보존된 접촉이 핵심적인 역할을 한다는 기존 연구를 출발점으로 삼는다. 이를 구체화하기 위해 저자들은 “H‑폼”이라는 새로운 구조 단위를 제안한다. H‑폼은 두 개의 서로 다른 2차 구조 요소(헬릭스 혹은 베타‑시트)에서 각각 중심 Cα 원자를 선택하고, 그 중심을 기준으로 앞뒤 두 잔기의 총 5잔기 길이의 백본 조각을 추출한다. 두 중심 Cα 원자 사이 거리가 8.5 Å 이하이면 접촉으로 간주한다.

기하학적 특성은 네 개의 연속 Cα 좌표를 이용해 축 방향을 추정하고, 축‑벡터 u와 중심 간 벡터 rba 사이의 각도 θa, θb, 그리고 두 축 사이의 비틀림 각 τab를 계산한다. 이렇게 정의된 (d, θa, τab, θb) 네 파라미터는 회전·이동에 불변이며, H‑폼을 정량적으로 비교할 수 있는 기준이 된다. 또한 서열적 거리 ℓ(두 중심 잔기의 서열 인덱스 차)과 필요시 용매 접근성 등 부가 정보를 함께 고려한다.

구조 정렬을 수행할 때는 기존의 전역 정렬 도구가 H‑폼 중심의 국소 정보를 놓치기 쉬운 점을 보완하고자, 저자들은 “zoom‑in” 기법을 적용한 맞춤형 정렬 파이프라인을 개발했다. 초기에는 두 H‑폼을 겹치게 하는 변환을 구하고, 일정 거리 임계값 이하의 추가 대응 잔기를 확대해가며 반복적으로 변환을 업데이트한다. 이 과정을 통해 H‑폼 쌍이 전체 구조 정렬에 포함되는지를 판단한다.

실험에서는 작은 단백질 CI2(2ci2)를 대상으로 7개의 H‑폼을 도출하고, SCOP 슈퍼패밀리 d.122.1 등에서 유사한 H‑폼을 탐색했다. CI2와 구조적으로 유사한 1vbw, 1mit 등과의 비교에서, 일부 H‑폼은 서열 유사도(BLOSUM62 점수 ≥0)와 기하학적 차이(Δd ≤1.5 Å, Δθ ≤0.6°, Δτ ≤0.8°)를 동시에 만족하며 정렬에 일관되게 포함되었다. 특히 A16‑I20, L49‑V47 등은 보존된 접촉으로서 폴딩 핵심에 기여함을 확인했다. 슈퍼패밀리 수준에서는 80여 쌍의 유사 H‑폼이 발견됐으며, 이 중 54쌍은 구조 정렬과 일치했다. 이러한 결과는 H‑폼이 단백질 패밀리 간에 공유되는 “포장 단위”로 작동할 가능성을 시사한다.

마지막으로 저자들은 H‑폼 기반 포장 단위가 기존의 로컬 구조 모티프(예: β‑핸들, 헬릭스‑헬릭스 접촉)와 차별화되는 점을 강조한다. H‑폼은 짧은 백본 조각이지만, 두 2차 구조 요소 사이의 정확한 입체적 배치를 포착하므로, 구조 예측 시 거리 제약이나 에너지 함수에 직접 삽입하기에 적합하다. 또한 보존된 서열 신호가 강해 진화적 분석에도 활용 가능하다. 전체적으로 본 연구는 단백질 구조의 핵심적인 포장 메커니즘을 정량화하고, 이를 통해 향후 ab‑initio 구조 예측 및 폴딩 코어 식별에 새로운 도구를 제공한다.