당기는 힘에 따라 변하는 FnIII10의 기계적 전개 경로

초록

본 연구는 전산적 전산모델을 이용해 섬유소단백질 FnIII10 도메인의 기계적 전개를 상수 힘과 상수 속도 조건에서 조사하였다. 두 단계 전개와 세 가지 상호 배타적인 중간 상태가 관찰되었으며, 전개 경로는 인장 강도에 따라 크게 달라진다. 낮은 힘·속도에서는 세 중간 상태가 모두 나타나지만, 높은 힘·속도에서는 N‑말단과 C‑말단 베타 스트랜드가 동시에 탈구된 중간 상태가 우세한다. 확장된 Jarzynski 등식으로 자유에너지 지형을 추정한 결과, 네이티브, 완전 전개, 그리고 중간 상태에 해당하는 세 개의 최소점이 확인되었다.

상세 분석

이 논문은 전통적인 원자 수준의 명시적 용매 시뮬레이션 대신, 암시적 물 모델과 단순화된 에너지 함수(E_loc, E_ev, E_hb, E_hp)를 적용한 전산적 접근법을 채택하였다. 이러한 모델은 계산 효율성을 크게 향상시켜 수천 개의 전개 궤적을 생성할 수 있게 하였으며, 이는 다중 전개 경로를 통계적으로 분석하는 데 필수적이다. 상수 힘 시뮬레이션에서는 50 pN부터 192 pN까지 여섯 가지 힘을 적용했으며, 상수 속도 시뮬레이션에서는 0.03 fm/MC step부터 1.0 fm/MC step까지 네 가지 속도를 사용하였다. MC 업데이트는 백본 각도 변화를 위한 편향 가우시안 스텝(BGS), 측쇄 각도에 대한 메트로폴리스 스텝, 그리고 전체 사슬의 강체 회전을 포함한다.

전개 궤적 분석은 β‑스트랜드 간 수소 결합 존재 여부를 기준으로 중간 상태를 정의하였다. β‑스트랜드 쌍이 30 % 이상의 원래 수소 결합을 유지하면 해당 쌍이 형성된 것으로 간주하고, 이를 통해 A‑G 7개의 스트랜드 중 어느 두 개가 탈구되었는지를 판별한다. 상수 힘 조건에서는 엔드‑투‑엔드 거리(L) 히스토그램을 이용해 메타안정 상태를 식별했으며, 피크의 위치는 중간 상태의 평균 거리(L_I)로, 피크 면적은 이론적으로 상태의 평균 수명을 나타낸다. 상수 속도 조건에서는 힘‑시간 곡선의 급격한 감소를 ‘파열 사건’으로 정의하고, 파열 전 최대 힘을 해당 상태의 안정성 지표(F_I)로 사용하였다.

Jarzynski 평등을 확장한 EJE 방법을 통해, 속도 구동 궤적에서 측정된 작업(W_t)과 거리(L)를 이용해 비변형 자유에너지 G_0(L)를 재구성하였다. 가중 히스토그램 기법을 적용해 G_0(L) 곡선을 얻었으며, 이는 네이티브 상태(L≈5 nm), 중간 상태(L≈12–16 nm), 완전 전개 상태(L≈30 nm)에 각각 최소점을 보인다. 특히, 높은 힘·속도에서 관찰된 ‘N‑C 말단 탈구’ 중간 상태는 자유에너지 장벽이 낮아 전개가 빠르게 진행되는 메커니즘을 제공한다.

결과적으로, 전개 경로는 인장 강도에 민감하게 반응하며, 중간 상태의 상대적 발생 빈도는 힘·속도에 따라 크게 변한다. 이는 실험적 AFM 측정에서 보고된 다중 단계 전개와 일치하며, FnIII10의 기계적 안정성이 세포 외 기계 신호 전달에 어떻게 기여할 수 있는지를 설명한다. 또한, 암시적 물 모델과 EJE 기반 자유에너지 재구성이 복잡한 단백질 전개의 정량적 이해에 유용함을 입증한다.