YAP/TAZ 기계감지 메커니즘을 밝히는 새로운 컴퓨테이셔널 모델

초록

본 논문은 ECM 강도와 같은 외부 기계적 신호를 세포 내 생화학적 신호로 변환하는 YAP/TAZ 경로를 정량화한 수학 모델을 제시한다. 모델은 FAK‑RhoA‑스트레스 섬유‑LIMK‑LATS‑YAP/TAZ 순환을 포함하며, 연성 기질에서도 FAK 활성화가 YAP/TAZ 활성을 회복시킬 수 있음을 예측한다. 또한, 분자 농도 변화가 강도‑반응 곡선에 미치는 영향을 분석하고, Hippo 경로와의 협동 효과를 LIMK와 LATS의 상호작용으로 설명한다.

상세 분석

이 연구는 YAP/TAZ 메커노센싱을 이해하기 위해 기존의 Hippo 신호전달 모델에 기계적 입력을 정량적으로 연결한 새로운 수학적 프레임워크를 구축하였다. 모델은 크게 세 부분으로 나뉜다. 첫 번째는 ECM 강도와 리간드 밀도에 의해 활성화되는 FAK( focal adhesion kinase )의 동역학이며, 이는 RhoA‑GTPase를 통해 스트레스 섬유 형성을 촉진한다. 두 번째는 RhoA‑활성화된 ROCK와 mDia가 주도하는 액틴 섬유 조직화 과정으로, 여기서 LIMK은 cofilin을 인산화해 섬유의 안정성을 높인다. 세 번째는 Hippo 핵심 인자인 LATS1/2와 YAP/TAZ의 인산화‑탈인산화 균형이다. 모델은 미분 방정식 형태로 각 단백질의 활성화/비활성화 속도와 상호작용 계수를 정의하고, 파라미터는 문헌값과 실험적 보정값을 혼합해 설정하였다.

주요 발견은 다음과 같다. (1) 연성(soft) 기질에서도 FAK를 과발현하거나 활성화하면 RhoA‑ROCK‑LIMK 경로가 강화되어 스트레스 섬유가 충분히 형성되고, 결과적으로 LATS 억제가 감소해 YAP/TAZ 핵 이동이 촉진된다. 이는 “기계적 구원(mechanical rescue)” 현상으로, 기존 실험에서 관찰된 FAK‑mediated YAP 활성화와 일치한다. (2) 모델은 FAK, RhoA, LIMK, LATS 각각의 초기 농도 변화를 시뮬레이션함으로써, 강도‑반응 곡선이 단일 임계점이 아니라 여러 구간으로 나뉘는 다중 스위치 특성을 보인다는 것을 보여준다. 예를 들어, LATS 농도가 높을 경우 강도에 대한 민감도가 크게 감소하지만, LIMK이 과발현되면 다시 민감도가 회복된다. (3) Hippo 경로와 기계적 경로의 협동 효과는 LIMK이 LATS와 직접적인 물리적 결합을 형성하거나, 간접적으로 ROCK‑mediated phosphorylation을 통해 LATS 활성을 억제한다는 가설을 통해 설명된다. 이 메커니즘은 기존에 “불명확한 시너지”로 남아 있던 현상을 수학적으로 정량화한다.

민감도 분석 결과, 모델 출력인 YAP/TAZ 핵 농도는 ECM 강도, FAK 활성화율, RhoA‑GTP 비율에 가장 높은 민감도를 보였으며, 특히 LIMK‑cofilin 경로는 스트레스 섬유의 동역학을 매개해 간접적으로 YAP/TAZ에 영향을 미친다. 파라미터 스페이스 탐색을 통해, 특정 조합(예: 높은 FAK·낮은 LATS·중간 LIMK)에서는 강도에 거의 무관하게 YAP/TAZ가 지속적으로 활성화되는 “강도 독립적” 상태가 도출되었다. 이는 암세포가 경질 기질에서도 성장 신호를 유지하는 메커니즘을 설명하는 데 활용될 수 있다.

전반적으로 이 모델은 실험적 데이터와의 정량적 일치를 통해 YAP/TAZ 메커노센싱의 복합성을 해소하고, 향후 약물 타깃 탐색이나 조직공학 설계에 있어 중요한 예측 도구로 활용될 가능성을 제시한다.