Likelihood-Based Heterogeneity Inference Reveals Non-Stationary Effects in Biohybrid Cell-Cargo Transport

기술 요약: 가능도 기반 이질성 추론을 통한 바이오하이브리드 세포-화물 수송의 비정상성 효과 규명

문제 정의
유전적으로 동일한 미생물로 구성된 경우라 할지라도, 생물학적 활성 물질(active matter) 시스템은 고유한 개체 간 운동성 변동성을 나타낸다. 수동적 물체(화물)가 이러한 이질적인 세포 집단에 의해 구동될 때, 화물의 역학은 이러한 변동성을 물려받게 되며, 이는 종종 비가우시안(non-Gaussian) 변위 통계로 이어져 집단 평균량의 해석을 복잡하게 만든다. 이전 연구들은 앙상블 평균 지표를 사용하여 이러한 시스템을 특징지어 왔으나, 궤적 데이터로부터 기저의 변동성(이질성)을 직접 추론하는 것은 특히 데이터가 희소하거나 시스템의 특성이 시간에 따라 변할 수 있는 경우(비정상성) 매우 어렵다. 각 궤적에 대해 먼저 매개변수를 추정하고 그 다음 분포를 피팅하는 표준적인 2단계 추정 방식은 정보를 손실하며, 특히 짧은 궤적에 대해 견고한 불확실성 경계를 제공하지 못한다.

방법론
저자들은 Dictyostelium discoideum (D. discoideum) 세포가 수동적인 폴리스티렌 비드(직경 61 µm)를 구동하는 바이오하이브리드 시스템을 조사한다. 실험 설정은 4시간 동안의 타임랩스 이미징을 포함하며, 이를 통해 177개의 비드 궤적이 생성되었다.

본 연구는 비드의 이질성을 추정하기 위해 전체 가능도 기반 추론(full likelihood-based inference) 접근법을 사용한다. 방법론은 다음과 같이 진행된다:

1. 모델 정식화: 평균 제곱 변위(MSD) 및 변위 자기상관함수(DACF)의 통계적 분석을 바탕으로, 2분 이상의 시간 척도에서 비드의 역학은 궤적 특이적 노이즈 강도 $\sigma^2_n$을 갖는 표준 브라운 운동으로 모델링된다.
2. 계층적 구조: 집단 내의 노이즈 강도 $\sigma^2_n$은 이질성 매개변수 $\theta = (\alpha, \beta)$에 의해 특징지어지는 감마 분포 $\Gamma(\alpha, \beta)$를 따른다고 가정한다.
3. 추론 체계: 각 궤적을 하나의 $\sigma^2_n$에 대한 최대 가능도 추정치(MLE)로 축소시킨 후 분포를 피팅하는 "2단계" 접근 방식과 달리, 이 방법은 모든 가능한 노이즈 강도에 대해 적분한다. $\theta$에 대한 로그 가능도는 전체 데이터셋의 로그 가능도를 합산함으로써 구축되며, 여기서 각 궤적의 기여도는 감마 분포에 의해 가중치가 부여된 궤적 가능도의 적분이다 (식 9).
4. 불확실성 정량화: 이 방법은 최대 가능도 추정치에서의 로그 가능도 함수의 헤시안(Hessian) 행렬을 활용하여 추론된 매개변수에 대한 자연스러운 불확실성 경계를 제공한다.
5. 비정상성 분석: 시간 의존성을 테스트하기 위해, 저자들은 4시간의 실험 전체에 걸쳐 슬라이딩 시간 창(16분 지속 시간)에 전체 가능도 접근법을 적용하며, 이질성을 평가하기 위해 짧고 희소한 궤적 조각(2분 간격)을 분석한다.

주요 결과

• 확산 역학의 검증: 통계적 분석은 2분 이상의 시간 척도에서 비드의 움직임이 확산적(Fickian)임을 확인해주지만, 변위 분포는 단일 가우시안 과정보다는 이질성을 나타내는 지수 꼬리(exponential tails)를 보인다.
• 전체 가능도의 우수성: 전체 가능도 접근법은 이질성 분포(노이즈 강도의 감마 분포)를 성공적으로 추론하며, 2단계 접근법보다 더 조밀하고 신뢰할 수 있는 불확실성 경계를 제공한다. 2단계 추정치는 전체 가능도 방법의 불확실성 경계 밖에 위치하며, 이는 후자에서의 정보 손실을 강조한다.
• 비정상성의 발견: 슬라이딩 시간 창 분석은 시스템이 정상 상태가 아님을 보여준다.
  • 시간적 드리프트: 실험 첫 2시간 동안, 추론된 이질성 분포의 평균과 분산이 크게 감소한다.
  • 준정상성(Quasi-Stationarity): 약 2시간 후, 시스템은 이질성 매개변수가 안정화되는 준정상 상태에 도달한다.
• 물리적 해석: 관찰된 "냉각 현상"(평균 확산 계수 및 분산의 감소)은 비드 주변의 비특이적 부착으로 인한 세포 축적에 기인한다. 이러한 축적은 힘의 상쇄를 유발하며, 입체적 반발력으로 인해 비드당 최대 세포 수를 제한함으로써 유효 확산과 세포 밀도 변동에 의한 변동성을 모두 감소시킨다. 세포의 정족수 인식(quorum-sensing) 행동 또한 전반적인 속도 저하에 기여할 수 있다.

의의 및 주장
논문은 가능도 기반 접근법이 짧은 궤적을 포함하는 "정보가 부족한" 상황에서 복잡한 활성 물질 시스템을 특징짓는 견고한 프레임워크를 제공한다고 주장한다. 중간 매개변수 추정을 건너뜀으로써, 이 방법은 정보 병목 현상을 피하고 자연스러운 불확실성 추정치를 제공한다.

본 연구의 주요 의의는 바이오하이브리드 수송 시스템에서 이질성의 시간 의존적 성질을 밝혀낸 데 있다. 저자들은 이질성을 정량화할 때 정상성을 가정하는 것이 왜곡된 결과를 초래할 수 있음을 입증한다. 이 방법을 통해 정상 역학의 창(window)을 식별할 수 있는 능력은, 이질성 자체를 정량화하는 것이 최종 목표가 아닌 경우에도 복합 시스템을 특징짓는 데 유용한 도구임을 보여준다. 연구는 개별 하위 시스템(세포와 비드)은 평형 상태에 도달할 수 있을지라도, 복합 시스템은 정상 상태에 도달하기 위해 훨씬 더 많은 시간이 필요하다는 동적 특징을 가지고 있다고 결론짓는다. 이는 표준 앙상블 평균 분석으로는 놓칠 수 있는 부분이다.