[論文レビュー] Simulation of mitochondrial metabolism using multi-agents system

Charles Lales, Nicolas Parisey+2 more|ArXiv.org|Jan 25, 2009

Mitochondrial Function and Pathology参考文献 6被引用数 5

ひとこと要約

本論文は、三次元空間的・動的リアリズムを備えたミトコンドリア代謝をシミュレートするためのマルチエージェントシステム（MAS）モデルを提案する。3次元エージェントベースのシミュレーションにより、分子間相互作用、酵素の局在、内膜における電子伝達を捉える。モデルはQサイクルやコenzymeQの酸化還元ダイナミクスといった複雑なプロセスを効果的に再現し、ミトコンドリア疾患や抗酸化機構に関する新たな知見を提供する。

ABSTRACT

Metabolic pathways describe chains of enzymatic reactions. Their modelling is a key point to understand living systems. An enzymatic reaction is an interaction between one or several metabolites (substrates) and an enzyme (simple protein or enzymatic complex build of several subunits). In our Mitochondria in Silico Project, MitoScop, we study the metabolism of the mitochondria, an intra-cellular organelle. Many ordinary differential equation models are available in the literature. They well fit experimental results on flux values inside the metabolic pathways, but many parameters are di$\pm$cult to transcribe with such models: localization of enzymes, rules about the reactions scheduler, etc Moreover, a model of a significant part of mitochondrial metabolism could become very complex and contain more than 50 equations. In this context, the multi-agents systems appear as an alternative to model the metabolic pathways. Firstly, we have looked after membrane design. The mitochondria is a particular case because the inner mitochondrial space, ie matricial space, is delimited by two membranes: the inner and the outer one. In addition to matricial enzymes, other enzymes are located inside the membranes or in the inter-membrane space. Analysis of mitochondrial metabolism must take into account this kind of architecture.

研究の動機と目的

• ミトコンドリア代謝における空間的組織、一時的相互作用、低分子数の動的特性を捉えることが難しい通常の常微分方程式（ODE）モデルの限界を克服すること。
• 二重膜構造を反映した空間的に明示的なシミュレーションフレームワークを構築し、内膜および膜間空間における酵素の局在を考慮すること。
• Qサイクルのような複雑な生化学的プロセスをシミュレートすること。これは、反応部位の不均一性とコエンザイムQの複数の酸化状態に起因し、従来のODEでは取り扱いにくい。
• 遺伝子変異や酵素の欠損がミトコンドリア機能に与える影響を、特にミトコンドリア疾患の文脈で調査すること。
• 病理状態下における還元型コエンザイムQの挙動と安定性をモデル化することで、抗酸化療法の設計を支援すること。

提案手法

• エージェントは、位置、速度、方向を有する三次元連続空間内の分子（代謝物、酵素、CoQ）を表し、空間的に分解能の高い相互作用を可能にする。
• 分子動力学はニュートン力学に基づいてシミュレートされ、力とモーメントは分子間距離、相互作用点、物理的性質（質量、慣性モーメント、摩擦係数）に基づいて計算される。
• 離散的な時間ステップを用いて、力とモーメントの統合により位置と方向を更新する：$\Delta\overrightarrow{X}_{i}(t+\delta t) \approx \overrightarrow{V}_{i}(t)\delta t + \frac{1}{2m}\overrightarrow{F}_{i}(t)\delta t^2$ および $\Delta\overrightarrow{\Theta}_{i}(t+\delta t) \approx \overrightarrow{W}_{i}(t)\delta t + \frac{1}{2I}\overrightarrow{T}_{i}(t)\delta t^2$。
• グリッドベースの近隣構造を用いてエージェントが周囲のエージェントを局所的に認識し、空間的相互作用の複雑さを簡略化しながらも三次元的リアリズムを維持する。
• リン脂質からの膜自己集合やタンパク質複合体の空間的配置といった、顕在的性質をモデルに組み込み、内膜の不均一性を現実的に再現する。
• 反応スケジューリングは近接性と結合ルールに基づき、温度揺らぎは分子種に依存するランダムトルク成分（$\beta_{rot}$）でモデル化される。

実験結果

リサーチクエスチョン

• RQ1 マルチエージェントシステムは、特に内膜における電子伝達鎖の空間的・動的複雑性を効果的にシミュレートできるか？
• RQ2 コエンザイムQの複数の酸化状態と複雑な空間的配置に依存するQサイクル機構を、エージェントベースの相互作用によって正確にモデル化できるか？
• RQ3 MASフレームワークは、酵素の欠損や変異がコエンザイムQの分布および酸化還元状態ダイナミクスに与える影響をどの程度正確に捉えることができるか？
• RQ4 リン脂質の自己集合とタンパク質複合体の統合から膜構造がどのように顕在的に形成されるかをモデルはシミュレートできるか？
• RQ5 還元型CoQの抗酸化挙動をシミュレートすることで、ミトコンドリア疾患の治療標的の同定にモデルはどのように寄与できるか？

主な発見

• マルチエージェントシステムは、ミトコンドリア酵素および代謝物の空間的ダイナミクスを効果的にシミュレートでき、内膜および膜間空間における複合体の局在を再現した。
• 異なる酸化状態のCoQを表すエージェントとComplex IIIの結合部位との間のルールベースの相互作用により、CoQを介した複雑な電子伝達プロセス、特にQサイクルを捉えた。
• 呼吸鎖活性や変異の変動に応じて、酸化型、セミキニン型、還元型の複数のCoQ種の分布を自然に再現した。
• 還元型CoQがその抗酸化形として安定化されることを示し、CoQをベースとする治療化合物の有効性を予測する根拠を提供した。
• エージェントベースのアプローチにより、遺伝子変異や基質供給不足に起因する低濃度および一時的相互作用といった、ODEモデルではうまく捉えきれない分子的相互作用をシミュレート可能であった。
• リン脂質エージェントとタンパク質複合体からの顕在的膜形成は、内膜の構造的複雑性を反映し、現実的な空間的組織を支持した。