Cell-cell adhesion cannot sustain extended follower streams in a minimal non-local model of leader-follower migration
牛津大學模型證實,在最高 100 微米交互範圍下最多僅能帶動 35% 細胞,單純短距黏附無法支撐長程遷徙。
- 一維偏微分方程模擬顯示,短距細胞黏附與共吸引力無法單獨維持長距離隨從細胞遷徙。
- 參數掃描呈現兩極化:過強的內部吸引會阻滯群體,過強的領導者吸引會拉斷細胞隊伍。
- 在 100 微米的交互極限下,模型最高只能帶動 35% 細胞,凸顯現有理論框架的結構性極限。
在交互距離被推至 100 微米的情況下,非局部對流擴散模型最多只能帶動 35% 的隨從細胞(follower cells)隨行。牛津大學數學生物學團隊透過一維偏微分方程模擬證實,單純依賴細胞間的短距離黏附力,根本無法在無增殖狀態下維持體內常見的長距離集體遷徙流。這項物理與數學層面的驗證,直接挑戰了生物學界解釋細胞移動的經典假設。
體內細胞長程遷徙在非局部對流擴散模型的挑戰
群體運動在生物界無處不在,而在多細胞系統中,由少數領導細胞引導大量隨從細胞前進的動態,是胚胎發育、傷口癒合與癌症侵襲的關鍵基礎。體內的領導細胞通常會對化學或機械梯度產生反應並主導方向,隨從細胞則依賴與鄰近細胞的互動來協調運動。在缺乏細胞增殖或長程化學感測的前提下,學界普遍將細胞間黏附視為維持長列隨從細胞流不致斷裂的核心機制。
為了以數學語言描述這種黏附交互作用,非局部對流擴散方程式(non-local advection-diffusion equations,考慮空間中一段距離內物質交互影響的偏微分方程系統)成為了標準的分析框架。這類方程式假定細胞的遷徙方向取決於有限範圍內感測到的細胞密度梯度。然而,過去運用此框架的模型大多只能產生範圍受限的靜止細胞簇,或是規模無法突破交互作用長度限制的移動脈衝。牛津大學團隊為此建立了一個極簡的遷徙模型,試圖驗證單靠局部吸引與黏附,是否足以支撐遠超細胞交互距離的長篇遷徙隊伍。
建立2000微米一維網格的偏微分方程系統
研究團隊建構了一個總長 2000 微米(μm) 的一維網格模型,將領導者與隨從者視為連續變化的密度分佈,並以耦合的偏微分方程(PDEs)來推演時間演化。方程式中規定,領導細胞以固定的每分鐘 1 微米速度向正 $x$ 方向前進,藉此模擬微環境訊號帶來的定向反應。隨從細胞的動態則由兩項參數決定:代表無方向性擴散的係數(設定為 $D_f = 10$ μm²/min),以及對領導者與其他隨從者的短距離非局部吸引力。
吸引力的強弱由交互作用強度 $\mu_{ff}$(隨從者彼此間)與 $\mu_{fl}$(隨從者對領導者)來參數化,並搭配一個關鍵的空間交互長度尺度 $\xi$。這個長度尺度界定了細胞能感測並產生互動的典型距離,在模擬實驗中,$\xi$ 被分別設定為 10、50 與 100 微米,大約相當於 1 到 10 個細胞的真實長度。模型同時納入了體積填充約束,以防止總體細胞密度超過最大值 $\kappa$ 時產生過度擁擠的情形。為了量化遷徙效率,團隊會精確計算模擬結束後,緊隨在領導細胞後方的連續隨從細胞質量($M_f$)佔初始總質量($M_0$)的比例。
參數掃描揭示細胞交互吸引力的兩種極端動態
針對不同 $\xi$ 值的系統性參數掃描顯示,模型在運算過程中穩定呈現出兩種截然不同的動態分岔。在第一種極端狀態中,當隨從者之間的吸引力 $\mu_{ff}$ 足夠強大時,領導細胞向前的拉力無法克服隨從細胞群內部的凝聚力。整個隨從細胞群會維持在初始位置並緊密抱團,導致能夠跟隨領導者前進的細胞質量比例直接歸零。
切換到第二種極端狀態,如果隨從者對領導者的吸引力 $\mu_{fl}$ 佔據主導優勢,起初位於領導細胞感測範圍內的隨從細胞確實會開始跟隨移動。但由於此時群體內部的 $\mu_{ff}$ 相對弱勢,這些已經起步的前端細胞無法把向前的拉力往後傳遞給位於大部隊尾端的人。這導致龐大的後方陣列被遺留在原地,僅有最前端的一小批細胞成功被帶走。這兩種動態機制清楚劃分了熱力圖上的零遷徙區與部分遷徙區,證明在現有連續體框架下,單純調整吸引力比例必定會面臨顧此失彼的物理限制。
100微米交互尺度極限下測得35%最高運輸率
運算數據進一步表明,被成功帶走的隨從細胞比例,雖然會隨著領導者吸引力的增加而單調上升,但很快就會觸及理論上的天花板。當交互長度 $\xi = 10$ 微米時,只有不到 5% 的隨從細胞能跟隨領導者移動。把 $\xi$ 提升到 50 微米,運輸比例的極限也只停滯在 20% 上下。即便將交互長度放寬到生理上幾乎達極限的 100 微米,被運輸的隨從細胞質量最大值也僅能達到 35%。
這組數據意味著對於給定的交互長度尺度,不管再怎麼無上限地增強細胞黏附力,隨從隊伍的長度都存在一個不可逾越的飽和點。若要在模型中進一步提高細胞運輸質量,唯一的方法是繼續加大 $\xi$ 值,但這會讓細胞間的通訊距離遠遠超出生物體內真實的短程黏附極限。模擬結果從根本上確認,在所有符合生物學合理性的參數範圍內,多數隨從細胞都無法單靠非局部對流擴散方程式完成長程遷徙。
數學連續體模型整合各向異性交互作用的必要
本次研究明確指出了標準短程非局部黏附數學框架的結構性限制:局部共吸引力與黏附力能形成小規模的緊密細胞團塊,卻不足以衍生出活體內動輒橫跨長距離的連貫遷徙流。這個運輸極限並非單純源自運算參數設定不當,而是該數學方程式家族固有的本質特徵。即使有極少數的孤立參數點能讓兩種吸引力勉強達到完美平衡,這種脆弱的狀態在該模型類別中也完全缺乏結構穩健性。
解開這項限制並不代表黏附作用在真實細胞遷徙中失去重要性,而是現有的連續體模型在描述機制的過程中,可能遺漏了決定性的動態變數。未來要精確透過數學運算模擬大規模集體細胞運動,可能必須引入各向異性交互作用(anisotropic interactions)、更長距離的生化信號傳導途徑,或是領導者對隨從者行為的反向應答。唯有將模型與外部基質的耦合效應一併考量,才有機會在維持質量守恆的鐵律下,重現綿延不絕的細胞遷徙網絡。
標準非局部對流擴散方程無法僅靠細胞黏附力重現長程集體遷徙,未來的連續體模型必須整合各向異性交互或長距信號傳導機制。