광공포 반응에서 편모 위상 동기화의 시간‑주파수 추적

초록

본 연구는 고강도 빛 자극(광쇼크) 후 초미세조류 Chlamydomonas reinhardtii 의 두 편모가 어떻게 재동기화되는지를 연속 웨이브렛 변환 기반의 복소 위상 동기화 지수(PSI)를 이용해 시간‑주파수 영역에서 정량화한다. 분석 결과는 40 Hz의 전형적인 가슴스트로크, 80 Hz까지 상승하는 역방향 유동, 그리고 재동기화 단계라는 세 가지 뚜렷한 주파수 대역과 위상 변화를 밝혀낸다. 또한, 기본 비트와 그 고조파가 광쇼크 전후에 지속적으로 존재함을 보여주어, 세포 간 변이와 무관하게 스펙트럼 모드의 균형 조절이 환경 스트레스에 대한 강인한 적응 메커니즘임을 시사한다.

상세 분석

이 논문은 미세조류 C. reinhardtii 의 편모 운동을 ‘비정상적이고 비선형적인’ 생물학적 진동자로 정의하고, 기존의 고정 주파수 기반 동기화 검출 방법이 이러한 전이 현상을 포착하기에 한계가 있음을 지적한다. 이를 극복하기 위해 저자들은 연속 웨이브렛 변환(CWT)을 이용해 신호의 순간 주파수와 위상을 동시에 추출하고, 복소 형태의 위상 동기화 지수(PSI, ΥΨ(t,f))를 도입한다. Ψ는 각 편모에서 얻은 웨이브렛 스펙트럼 W₁(t,f), W₂(t,f) 의 정규화된 곱으로 정의되며, 시간‑주파수 영역에서 위상 일치 정도를 0~1 사이의 값으로 나타낸다. 기존의 복소 코히어런스와 차별화되는 점은 진폭 정보를 사전에 정규화함으로써 순수 위상 관계만을 강조한다는 점이다.

실험적으로는 13개의 단일 세포를 마이크로파이펫에 고정하고, 1000 fps로 25 s 길이의 영상들을 촬영하였다. 흐름장 측정을 위해 250 nm 폴리스티렌 비드를 이용한 고속 입자 영상측정(GPV)을 적용했으며, 세 개의 ROI(좌·우 Y=±1, 전방 Y=0)에서 x·y 축 속도 신호를 추출했다. 광쇼크는 470 nm LED로 50 ms, 약 10 W m⁻²의 펄스를 3‑5 s 간격으로 가했다.

시뮬레이션 부분에서는 두 위상 진동자를 아델러 방정식 형태로 모델링하고, Gaussian 백색 잡음(디퓨전 계수 D)과 다양한 결합 강도 B 를 적용해 동기화·비동기화 전이를 재현한다. CWT는 Morlet 웨이브렛(중심 주파수 f* ≈ 1 Hz)으로 수행했으며, 주파수 해상도와 시간 해상도 사이의 트레이드오프를 최소화하도록 Q‑factor가 일정하도록 설계하였다.

시간‑주파수 PSI 분석 결과는 다음과 같다. (i) 광자극 전, 두 편모는 40 Hz 대역에서 높은 |ΥΨ| (≈0.9)와 위상 차 0 rad를 보이며 가슴스트로크 형태로 동기화된다. (ii) 광쇼크 직후, Ψ는 60‑80 Hz 대역으로 급격히 이동하고, |ΥΨ|는 일시적으로 감소하지만 80 Hz에서 새로운 위상 고정(≈π rad)과 함께 강한 반위상 동기화가 나타난다. 이는 고주파수 ‘역방향 유동’ 모드가 지배함을 의미한다. (iii) 자극 후 200‑300 ms 사이에 Ψ는 다시 40 Hz 대역으로 회귀하면서 |ΥΨ|가 서서히 회복되고, 위상 차는 0 rad로 복귀한다. 이 과정에서 1차 고조파(80 Hz)와 2차 고조파(120 Hz)의 에너지 비율이 변하는데, 특히 역방향 단계에서 1차 고조파가 주된 에너지원을 차지한다는 점이 눈에 띈다.

이러한 결과는 편모가 단순히 두 개의 이산적인 비트(가슴스트로크 vs 역방향) 사이를 전환하는 것이 아니라, 여러 주파수 모드와 그 고조파가 동시 존재하며, 외부 자극에 따라 그 상대적 가중치가 조절된다는 새로운 관점을 제공한다. 또한, PSI의 시간‑주파수 구현은 비정상적인 생물학적 진동자들의 순간적 결합 강도와 위상 관계를 정량화하는 강력한 도구로, 향후 신경전달, 심장 박동, 섬모 운동 등 다양한 분야에 적용 가능성을 시사한다.