Pathway variability, coat stiffening and mechanical adaptation during clathrin-mediated endocytosis

기술 요약: 클라트린 매개 엔도사이토시스에서의 경로 가변성, 코트 경직도 및 기계적 적응

문제 제기
클라트린 매개 엔도사이토시스(CME)는 진핵세포가 수용체와 영양분을 내부화하는 주요 메커니즘이다. 이 과정에서 클라트린 트리스켈리아(triskelia)는 다각형 격자를 형성하여 막을 리모델링하고, 최종적으로 코티드 피트(coated pit)를 거쳐 소낭(vesicle)으로 변형시킨다. 그러나 실험적 관찰에 따르면 상당한 경로 가변성이 나타난다. 즉, 조립체는 평평한 플라크(plaque) 상태로 지속되거나, 부분적인 함입 단계에서 중단되거나, 혹은 성공적으로 소낭으로 닫히기도 한다. 이러한 서로 다른 운명을 결정하는 결정 요인은 아직 해결되지 않은 상태이다. 기존 모델들은 코트의 경직도와 선호되는 곡률을 고정된 재료 파라미터로 취급하는 경우가 많으며, 이러한 특성이 조립 과정 자체로부터 어떻게 동적으로 발생하는지는 설명하지 못한다. 또한, 미시적인 격자 조립, 위상학적 결함(오각형/칠각형), 그리고 거시적인 막 곡률 생성 사이의 관계도 불분명하다.

방법론
저자들은 이산적인 클라트린 에이전트를 적응형 연속체 막에 결합하는 하이브리드 키네틱 몬테카를로(KMC) 시뮬레이션 프레임워크를 개발하였다.

• 이산적 구성 요소: 개별 클라트린 트리스켈리아는 연속적인 막 표면 위의 이동 가능한 노드로 모델링된다. 각 트리스켈리온은 선호하는 다리 길이($L_0 = 18.4$ nm), 평면 내 각도($\phi_0$), 그리고 퍼커 각도($\chi_0$)로부터의 편차에 대한 비용을 갖는 탄성 해밀토니안(Hamiltonian)을 가진다. 노드는 확산, 결합/해리하며, 확률적이고 에너지 의존적인 비율로 결합을 형성하거나 해제한다.
• 연속체 구성 요소: 막은 면적($A$)과 평균 곡률($H$)로 정의되는 구형 캡(spherical-cap) 기하학으로 기술된다. 막의 곡률은 확률적인 조립 이벤트 사이에서 코트의 기계적 상태에 따라 단열적(adiabatically)으로 적응한다.
• 핵심 혁신: 이전 모델들과 달리, 이 프레임워크는 유효 코트 굽힘 강성($\kappa_{cl}^{eff}$)과 선호되는 코트 곡률($H_{min}^{eff}$)이 성장 과정 중에 자연스럽게 발현되도록 허용한다.
• 검증: 본 모델은 코트 기하 구조의 고속 원자 힘 현미경(AFM) 재구성 및 HS-AFM 절단 연구를 모사한 in silico 나노해부(nanodissection) 실험을 포함한 실험 데이터와 검증되었다.

주요 기여 및 결과

1. 동적 경직화 및 곡률 기억:
   고정된 곡률을 가진 막 위에서 성장하는 코트를 시뮬레이션한 결과, 코트가 성숙함에 따라 유효 굽힘 강성이 거의 두 자릿수(약 $20 k_B T$에서 $\sim 10^3 k_B T$까지) 증가함을 밝혔다. 이러한 경직화는 기하학적 전이에서 기인한다. 평평하거나 느슨하게 연결된 격자에서는 곡률 변화가 부드러운 각도 재배열을 통해 수용된다. 반면, 연결된 곡선 격자에서는 곡률 변화가 잔류 면적 변형(residual area strain)을 유발하며, 이로 인해 시스템은 훨씬 더 단단한 신축 모드를 통해 저항하게 된다.
   동시에, 코트는 "곡률 기억"을 발달시킨다. 유효 에너지 최소화 곡률($H_{min}^{eff}$)은 개별 트리스켈리아의 미시적 선호도에서 벗어나, 재료가 조립되었던 곡률 쪽으로 이동한다. 이는 각 삽입된 트리스켈리온이 국소적 참조 기하 구조를 격자에 기록하는 방식에 따라, 성장 곡률들의 가중 평균으로 모델링된다.

2. 두 단계의 순차적 운명 게이트:
   본 모델은 이산적 연결성과 연속체 역학 사이의 상호작용을 바탕으로 조립의 운명(평평함, 중단됨, 또는 닫힘)을 결정하는 두 개의 순차적 게이트를 식별한다.

• 게 Gate 1 (연결성): 코트는 일관된 굽힘 응력을 전달할 수 있을 만큼 충분한 격자 연결성을 확보해야 한다. 연결성이 낮은 코트는 흐물흐물한 상태로 남아 평평한 플라크로 성장한다. 충분한 연결성이 확보되면 굽힘-신축 전이가 트리거되어, 경직화와 곡률 생성을 유도한다.
• Gate 2 (에너지 지형): 경직화가 발생하면, 운명은 시스템의 이동하는 에너지 최소점과 폐쇄 경계($\phi = 1$) 사이의 경쟁에 의해 결정된다. 이 최소점의 위치는 막 장력, 중합 에너지, 그리고 성장 이력(곡률 기억)에 따라 달라진다. 시스템이 폐쇄 경계에 도달하기 전에 최소점이 경계 너ars를 넘어가면 코트가 닫히고, 시스템이 최소점을 접근 가능한 영역 내에서 포착하면 피트가 멈춘다(중단됨).

3. 위상학적 결함의 역할:
   시뮬레이션은 위상학적 결함(오각형 및 칠각형)이 곡률 생성의 주요 트리거가 아님을 보여준다. 곡률은 눈에 보이는 비육각형 결함이 카운트되기 전에도 시작될 수 있다. 대신, 결함은 충분히 연결되고 경직된 코트에 의해 생성된 곡률를 뒤따르며 안정화하며, 곡률이 이미 확립된 후에야 닫힌 케이지를 위한 위상학적 요구 사항을 충족시킨다.

4. 비유클리드 탄성 이론:
   저자들은 관찰된 성장 법칙을 예측하는 비유클리드 탄성 기반의 해석적 이론을 도출하였다. 이 이론은 유효 강성과 선호되는 곡률이 기하학적 가중치($w \propto A^2 H^2$)와 이력 의존적 곡률 변수의 함수라고 상정한다. 이 이론은 피팅 파라미터 없이 다양한 성장 곡률과 막 장력으로부터의 시뮬레이션 데이터를 보편적인 경향성으로 성공적으로 통합(collapse)시켰다.

5. *In Silico 절단을 통한 검증:*
   예응력(prestress)과 기억 개념을 테스트하기 위해, 저자들은 연결된 코트를 절단하는 것(나노해부 모사)을 시뮬레이션하였다. 결과는 절단이 유효 굽 bending 강도를 감소시키고 유효 에너지 최소화 곡률을 증가시킨다는 것을 보여주었다. 이는 코트가 연결된 격자 내에 성장 이력을 저장하고 있음을 확인시켜 준다. 즉, 초기 성장 단계의 저곡률 재료를 제거하면 남은 코트의 이완된 상태가 더 높은 곡률 쪽으로 이동한다.

의의
본 논문은 각 운명에 대해 별도의 메커니즘을 도입하지 않고도 CME 결과의 가변성을 설명할 수 있는 통합적인 물리적 프레임워크를 제공한다고 주장한다. 코트의 강성과 선호되는 곡률이 조립 과정 자체에 의해 생성되는, 이력에 의존하는 창발적 특성임을 입증함으로써, 저자들은 이산적인 분자 조립과 연속체 막 역학 사이의 간극을 메웠다.

이 연구는 CME 모델링의 입력과 출력을 재정의한다. 즉, 재료 파라미터를 미리 규정하는 대신, 환경(예: 막 장력)과 조립 이력이 어떻게 유효한 기계적 상태를 공동으로 결정하는지를 보여준다. 이는 경로 가변성이 단순한 오류가 아니라, 코트의 성숙이 자신이 존재하는 에너지 지형을 재형성하는 기계적 피드백 루프에 의해 제어되는 자기 조립 과정의 본질적인 특징임을 시사한다. 본 연구의 결과는 평평한 코트, 정지된 코트, 닫힌 코트를 단일한 물리 법칙 아래에서 화해시키며, 클라트린 코트의 "운명"이 그것이 성장한 기계적 환경의 반영임을 강조한다.