[論文レビュー] A new computational method for a model of C. elegans biomechanics: Insights into elasticity and locomotion performance

Netta Cohen, Thomas Ranner|arXiv (Cornell University)|Feb 16, 2017

Genetics, Aging, and Longevity in Model Organisms参考文献 40被引用数 14

ひとこと要約

本稿では、線形変形不能性と局所的長さ制約を用いて、C. elegans の体を非伸縮性な弾性ビームとしてモデル化することで、多様な波状歩行と体形状を、さまざまな環境下で正確にシミュレートするための新しい有限要素法を提示する。主な貢献は、安定で半陽的数値スキームであり、効果的な運動を実現するための有効なヤング率の範囲（10–1000 kPa）を予測している。これは、抵抗性媒体を通過する際の速度調節が、筋力の変化によらず主要な制御機構であることを裏付けている。

ABSTRACT

An organism's ability to move freely is a fundamental behaviour in the animal kingdom. To understand animal locomotion requires a characterisation of the material properties, as well as the biomechanics and physiology. We present a biomechanical model of C. elegans locomotion together with a novel finite element method. We formulate our model as a nonlinear initial-boundary value problem which allows the study of the dynamics of arbitrary body shapes, undulation gaits and the link between the animal's material properties and its performance across a range of environments. Our model replicates behaviours across a wide range of environments. It makes strong predictions on the viable range of the worm's Young's modulus and suggests that animals can control speed via the known mechanism of gait modulation that is observed across different media.

研究の動機と目的

• ニュートン流体および粘弾性ゲルを含む、多様な物理的環境下での C. elegans の運動をシミュレートするための堅牢な計算フレームワークの構築。
• 任意の体形状と歩行ダイナミクスを捉える非線形初期境界値問題を用いた波状運動のバイオメカニクスのモデル化。
• 抵抗性環境下での運動性能と材料特性を関連付けることで、虫の体の有効な弾性率（ヤング率）の推定。
• 筋力の変化によらず、歩行の変更がどのように異なる媒体での速度最適化に寄与するかを解明。
• 神経・機械的統合応用を想定した、軟体非セグメンテッド運動体のシミュレーションに適した数値的安定性と効率性を備えた手法の提供。

提案手法

• 非線形初期境界値問題として定式化されたモデルは、質量保存則に基づく局所的長さ制約を用いて非伸縮性を組み込む。
• 中線曲線の区分線形近似と、区分定数の線張力による局所的非伸縮性を実現する、新規のパラメトリック有限要素法を考案。
• 半陽的時間ステッピング戦略により、各時間ステップで線形方程式系が得られ、数値的安定性と効率性が向上。
• 曲率を新たな独立変数として導入することで、連続的区分線形有限要素の使用が可能となり、収束性が向上。
• バーレットら（2010）のインスパイアを受けた人工接線運動を用いることで、メッシュの歪みを回避し、大変形時におけるロバスト性を確保。
• 内部圧力力の自然な適用が可能であり、局所的に保存されない距離を持つ系（例：ペルストラチックな波状運動体）にも適応可能。

実験結果

リサーチクエスチョン

• RQ1 C. elegans がアガー様および流体環境で効果的な前進力発生を実現するためのヤング率の有効範囲は何か？
• RQ2 筋力の変化によらず、歩行の変更が異なる抵抗性媒体下での運動速度に与える影響は何か？
• RQ3 実際のバイオメカニカル制約のもとで、ねじれやオメガターンを含む多様な体形状を、計算モデルが正確にシミュレートできるか？
• RQ4 局所的非伸縮性と曲率に基づく有限要素法が、軟体波状運動体のシミュレーションにおける数値的安定性をどのように向上させるか？
• RQ5 神経制御を明示的に組み込まず、僅かに筋肉活性化関数を指定するだけで、運動性能をどれほど正確に予測できるか？

主な発見

• 本モデルは、C. elegans が抵抗性環境下で推進力を発生させられる有効なヤング率の範囲を10–1000 kPaと予測している。下限は推進力発生能力、上限は破壊を避ける能力に由来する。
• 蠕虫は、媒体の変化に応じて歩行を変更することで、速度を最適化できることを確認した。これは、歩行の変更が運動性能の主要な制御機構であることを裏付けている。
• 本手法は、従来のシミュレータが困難とする、ねじれやループ状の体形状を含む複雑な形状のシミュレーションにおいても安定した結果を得られた。
• モデルの予測は、実験的観察および先行の定性的データと整合しており、スニットマンら（2010a）およびファン・エンら（2010）による初期の弾性率推定値を精緻化している。
• 本手法は、ニュートン流体および粘弾性環境の両方で正確な運動シミュレーションが可能であり、今後の内部粘性や遊泳への移行の拡張にも応用可能である。
• 本定式化は、将来的に神経筋制御システムとの統合が可能であり、複雑なナビゲーションおよび環境適応のクローズドループシミュレーションを可能にする。