システム免疫学は、 システム生物学 の一分野であり、数学的 アプローチと計算 手法を用いて、免疫システム の細胞・分子ネットワーク 内の相互作用を解析する。[ 1 ] 免疫システムは、その構成要素と機能に関して「還元主義的」アプローチを用いて徹底的に分析されてきたが、免疫システムは多数の構成要素間の相互作用に大きく依存しているため、個々の構成要素の特性を研究するだけでは、その全体的な機能を容易に予測することはできない。システム免疫 学 は 、生体内 実験よりもコンピューター内 実験 に重点を置いている。
実験免疫学および臨床免疫学 における近年の研究は、自然 免疫系と獲得免疫系 の両方のダイナミクスを 議論する数理モデル の開発につながっています。[ 2 ] これらの数理モデル の多くは、生体内で は検証できないプロセスをコンピューター上で 検証するために用いられています。これらのプロセスには、 T細胞 の活性化、がん免疫相互作用 、様々な免疫細胞 (T細胞 、B細胞 、好中球 など)の移動 と死、そして 臨床試験 を実施することなく特定のワクチン や薬剤に対する 免疫系が どのように反応するかなどが含まれます。[ 3 ]
免疫細胞のモデリング技術 免疫学における数学モデルの使用方法を説明するスキーム。 免疫学における モデリング 技術には、定量的 アプローチと定性的な アプローチがあり、それぞれに長所と短所がある。定量的モデルは 、特定の反応速度パラメータと、特定の 時点 または濃度 におけるシステムの挙動を予測する。短所は、少数の反応にしか適用できず、一部の反応速度パラメータ に関する事前知識が必要となることである。一方、定性モデルは より多くの反応を考慮できるが、その代わりにシステムの反応速度 に関する詳細は少なくなる。唯一の共通点は、どちらのアプローチも、成分数が大幅に増加すると簡潔さを失い、役に立たなくなることである。[ 4 ]
常微分方程式モデル 常微分方程式(ODE)は 、生物システム のダイナミクスを 記述するために使用されます。ODEは、 ミクロ 、メソスコピック、マクロの スケールで連続変数 を調べるために使用されます。方程式は 、タンパク質の濃度 、転写因子、 細胞 タイプの数など、観測された変数 の時間的変化 を表します。これらは通常、免疫シナプス 、微生物 認識、細胞移動を モデル化する ために使用されます。過去10年間、これらのモデルは、TCR のアゴニストリガンドに対する感受性と CD4 およびCD8 コレセプター の役割を研究するために使用されてきました。 これらの方程式の運動速度は、相互作用する種の 結合 速度と解離速度 で表されます。これらのモデルは、ネットワーク 内の相互作用する各分子 の濃度 と定常状態 を表すことができます。ODE モデルは 線形方程式 と非線形 方程式で定義されますが、非線形 方程式はコンピューター(in silico )でのシミュレーション や分析 が容易なため、より頻繁に使用されます。このモデルの限界は 、このモデルを適用するには、ネットワーク ごとに各分子の運動学を知る必要があることである。 [ 5 ]
ODEモデルは、 抗原が B細胞受容体 に結合する仕組みを調べるために使用されました。このモデルは1122の方程式と6つのシグナル伝達タンパク質 で表現されるため、非常に複雑でした。この研究に使用されたソフトウェアツールはBioNetGenでした。[ 6 ] モデルの 結果はin vivo実験に基づいてい ます 。[ 7 ]
エプスタイン・バーウイルス(EBV)は、若年層における 単核球症 の発生率が高いことを説明する3つの仮説 を検証するため、12の方程式を用いて数学的にモデル化され ました。数値シミュレーションを行った結果、最初の2つの仮説のみがモデルによって支持されました。[ 8 ]
偏微分方程式モデル 偏微分方程式 (PDE)モデルは ODE モデルの拡張版であり、時間 と空間の 両方で各変数の 時間発展を 記述します。PDEは 、病原体の経路の感知と認識における 連続変数を モデル化する ためにミクロレベル で使用されます。また、タンパク質が どのように相互作用し、免疫シナプス 内でその動きがどこに向けられているかを記述するための生理学的モデリング[ 9 ] にも適用されます。これらの導関数は 時間だけでなく 空間 についても計算されるため偏微分です。細胞分裂 の年齢などの生理学的変数を空間変数の代わりに使用できる場合があります。細胞の空間分布 を考慮したPDEモデルと ODE モデルを比較すると、PDE の方が 計算 負荷が高いです。空間ダイナミクスは、3 次元区画内での細胞の動きを 記述するため、細胞シグナリング の重要な側面です。T細胞は3 次元の リンパ節 内を移動しますが、TCR は 細胞膜 の表面にあるため 2 次元区画内を移動します。[ 10 ] タンパク質 の空間分布は 、特にT細胞 刺激時に免疫シナプスが形成される際に重要であるため、このモデルは T細胞が 弱いアゴニストペプチド によって活性化される研究に使用された。[ 11 ]
粒子ベース確率モデル 粒子ベースの確率モデルは、 ODE モデルのダイナミクス に基づいて得られます。このモデルが他のモデルと異なるのは、モデルのコンポーネントを、以前のモデルのように連続的ではなく、離散変数 として扱う点です。これらは、免疫特異的伝達経路と免疫細胞 -ガン相互作用において、それぞれ ミクロレベル とメソスコピックレベルの粒子 を調べます。モデルのダイナミクスは マルコフ 過程によって決定され、この場合、システムの各可能な状態の確率を 時間とともに 微分方程式 の形式で表します。この方程式は解析的に解くのが難しいため、コンピューター上でのシミュレーションは運動論的 モンテカルロ 法として実行されます。シミュレーション は通常、化学反応速度定数から導出される反応定数を使用して、反応が起こるかどうかを予測するGillespie アルゴリズムを使用して実行されます。 確率的シミュレーション は計算量が多いため、モデルの サイズと範囲は制限されます。
確率論的シミュレーションは、 T細胞 において重要なシグナル 伝達分子であるRasタンパク質 が活性型と不活性型を取り得ることを示しました。このシミュレーションは、刺激を受けたリンパ球 集団が活性型と不活性型のサブポピュレーションを持つことを明らかにしました。[ 12 ]
共受容体は T細胞活性化 の初期段階で重要な役割を果たしており、確率論的シミュレーションは 相互作用を説明するとともにリンパ節 内の遊走細胞 をモデル化するために使用された。[ 13 ]
このモデルはリンパ系 におけるT細胞の 増殖を調べるために使用された。[ 14 ]
エージェントベースモデル 糖尿病I型におけるCD8+ T細胞とβ細胞の相互作用の要約 エージェントベースモデリング(ABM) は、観測対象となるシステムの構成要素を離散 エージェントとして扱い、個々の分子 または細胞を表現する モデリング の一種です。このシステムでは、構成要素(エージェント)は他のエージェントや環境と相互作用することができます。ABM は マルチスケール レベルでイベントを観測する可能性を秘めており、他の分野でも普及が進んでいます。例えば、糖尿病I 型におけるCD8+ T細胞 とβ細胞 の相互作用のモデリング[ 15 ] や、白血球 のローリングと活性化のモデリング [ 16 ] などに利用されています。
ブールモデル ロジックモデルは、 細胞 のライフサイクル 、免疫シナプス 、病原体認識、ウイルスの侵入などを ミクロレベル およびメソスコピックレベルでモデル化するために使用されます。ODEモデルとは異なり、 ロジスティックモデルでは、相互作用する種の 動態 や濃度 に関する詳細は必要とされません。各生化学種は ネットワーク 内のノード として表され、通常は2つの、例えばON/OFF、高/低、活性/非活性など、有限数 の離散状態を持つことができます。通常、2つの状態のみを持つロジックモデル はブールモデルと見なされます。分子が OFF状態にある場合、それは分子が システムに変化をもたらすほど十分に高いレベルで存在しないことを意味し、濃度が ゼロであることを意味するわけではありません。したがって、ON状態にある場合、分子は反応を開始するのに十分な量に達しています。この手法はカウフマンによって初めて導入されました。このモデルの限界は、システムの定性的な近似値しか提供できず、同時発生しているイベントを完全にモデル化できないことです。[ 17 ]
この手法は、体液性免疫系における親和性成熟や高変異[ 18 ] 、病的リウマチ因子に対する耐性[ 19 ] など、免疫系 における特殊な経路の探索に用いられてきた。このモデルをサポートするシミュレーションツールとしては、DDlab [ 20 ] 、 Cell-Devs [ 21 ] 、 IMMSIM-C などがある。IMMSIM -C はコンピュータプログラミングの知識を必要としないため、他のツールよりも頻繁に利用されている。このプラットフォームは公開ウェブアプリケーションとして利用可能であり、プリンストン大学、ジェノバ大学など、様々な大学の学部レベルの免疫学コースで利用されている。[ 22 ]
システム免疫学において、ステートチャート を用いたモデリングには、これまでRhapsodyのみが使用されてきました。Rhapsodyは、 ステートチャートを実行可能な Java およびC++ コード に変換できます。
この手法はインフルエンザウイルス感染 モデル の構築にも用いられました。結果の一部は以前の研究論文とは一致せず、ブールネットワークは活性化 マクロファージ の量が若いマウスと老齢マウスの両方で増加したことを示しましたが、一方で減少を示唆する結果もありました。[ 23 ]
SBML (システム生物学マークアップ言語)は、 常微分方程式 を用いたモデルのみを扱うことを目的としていましたが、最近ブールモデル も適用できるようにアップグレードされました。ほぼすべてのモデリングツールはSBMLと互換性があります。 ブールモデル を用いたモデリング 用のソフトウェアパッケージとしては、BoolNet [ 24 ] 、 GINsim [ 25 ] 、Cell Collective [ 26 ] などがあります。
微分方程式 を用いてシステムをモデル化するには、コンピュータツールがモデル構築 、キャリブレーション 、検証 、解析 、シミュレーション 、可視化 といった様々なタスクを実行する必要があります。これらの基準を満たす単一のソフトウェアツールは存在しないため、複数のツールを組み合わせる必要があります。[ 27 ]
ジンシム GINsim [ 28 ] は、離散変数 に基づいて遺伝子ネットワーク を生成し、シミュレートするコンピュータツールです。制御グラフと論理パラメータに基づいて、GINsim [ 29 ] はシステムの時間的進化 を計算し、状態がノード で、遷移が 矢印で 表された状態遷移グラフ (STG)として返されます。これは、 TCR およびTLR5パスウェイが活性化されたときに T 細胞が どのように反応する かを調べるために使用されました。これらのプロセスは、個別にも組み合わせても観察されました。最初に、TCR およびTLR5パスウェイの両方の分子マップとロジックモデルが構築され、次にマージされました。分子マップは、文献のデータや KEGG [ 31 ]や Reactome [ 32 ] などのさまざまなデータベースに基づいてCellDesigner [ 30 ] で作成されました。論理モデル は GINsim [ 33 ] によって生成され、各コンポーネントは 0 または 1 の値を持ち、変更されると追加の値を持ちます。次に、このネットワーク において論理ノードと呼ばれる各構成要素に論理ルール が適用されます。これらの構成要素を統合することで、最終的なモデルは128個のノードで構成されます。モデリングの結果は実験結果と一致し、TLR5は CD4+ T細胞 の共刺激 受容体であることが実証されました。[ 34 ]
ブールネット Boolnet [ 35 ] はブールネットワークの再構築、分析、可視化のためのツール を含むRパッケージ です。[ 36 ]
セルコレクティブ Cell Collective [ 37 ] は、科学者が数式やコーディング をすることなく、生物学的モデルを構築、分析、シミュレーションすることを可能にする科学プラットフォームです。Cell Collectiveには、個々の実体(タンパク質 、遺伝子 、細胞 など)に関する知識を動的モデル へと拡張する知識ベースコンポーネントが組み込まれています。データは定性的 ですが、相互作用する種間の動的な関係を考慮に入れています。モデルはリアルタイムでシミュレーションされ、すべてウェブ上で行われます。[ 38 ]
バイオネットジェン BioNetGen(BNG)は、遺伝子制御 、細胞シグナル伝達 、代謝 といった複雑なシステムのルールベースモデリングに用いられるオープンソースソフトウェアパッケージです。このソフトウェアは、 グラフ を用いて様々な分子 とその機能ドメイン を表現し、それらの相互作用を説明するルールを用いています。免疫学の分野では、TLR-4カスケードの細胞内シグナル伝達経路のモデル化に用いられました。[ 39 ]
DSAIRM DSAIRM(免疫応答モデリングへの動的システムアプローチ)は、コーディングの事前知識なしに感染 と免疫 応答のダイナミクスを研究するために設計されたRパッケージ です。[ 40 ]
その他の有用なアプリケーションと学習環境としては、Gepasi、[ 41 ] [ 42 ] Copasi、[ 43 ] BioUML 、[ 44 ] Simbiology (MATLAB) [ 45 ] Bio-SPICE [ 46 ]などがある。
会議 合成免疫学 およびシステム免疫学に関する第1回会議が、CSFとETHチューリッヒの共催でアスコナで開催されました。[ 47 ] 2019年5月初旬に開催され、様々な科学分野から50名以上の研究者が参加しました。行われた全ての発表の中で、最も優れた発表は、TCRエピトープのスクリーニングのためのプラットフォームを発明したゴビンダ・シャルマ博士によるものでした。
ニューヨークのコールド・スプリング・ハーバー研究所(CSHL)[ 48 ] は、2019年3月に、免疫系を深く研究する実験生物学、計算生物学、数理生物学者間の意見交換を目的とした会議を主催しました。会議の議題は、モデリングと制御ネットワーク、合成生物学とシステム生物学の将来、そして免疫受容体でした。[ 49 ]
さらに読む
参考文献 ^ Villani, Alexandra - Chloè; Sarkizova; Hacohen (2018年4月26日). 「システム免疫学:免疫システムのルールを学ぶ」 . Annual Review of Immunology . 36 : 813–842 . doi : 10.1146/annurev-immunol-042617-053035 . PMC 6597491. PMID 29677477 . ^ Eftimie, Raluca; Gillard, Joseph J.; Cantrell, Doreen A. (2016年10月6日). 「免疫学の数理モデル:最新技術と将来の研究方向」 . Bulletin of Mathematical Biology . 78 (10): 2091– 2134. doi : 10.1007/s11538-016-0214-9 . PMC 5069344. PMID 27714570 . ^ Schubert, Charlotte (2011年5月4日). 「システム免疫学:複雑性の解明」 . Nature . 473 (7345): 113–114 . doi : 10.1038/nj7345-113a . PMID 21548192 . ^ Thakar, Juilee; Poss, Mary; Albert, Réka; Long, Gráinne H; Zhang, Ranran (2010年8月20日). 「 免疫応答の動的モデル:理想的な詳細レベルとは?」 . 理論生物学および医療モデリング . 7 35. doi : 10.1186/1742-4682-7-35 . PMC 2933642. PMID 20727155 . ^ Kim, Peter S; Levy, Doron; Lee, Peter P (2009年1月1日). 自己制御ネットワークとしての免疫システムのモデリングとシミュレーション . Methods in Enzymology. Vol. 467. pp. 79– 109. doi : 10.1016/S0076-6879(09)67004-X . ISBN 978-0-12-375023-5 . PMID 19897090 .^ "BioNetGen" . BioNetGen . ^ Barua, Dipak; Hlavacek, William S; Lipniacki, Tomasz (2012年6月18日). 「B細胞抗原受容 体 シグナル伝達における初期イベントの計算モデル:LynとFynの役割の解析」 . The Journal of Immunology . 189 (2): 646– 658. doi : 10.4049/jimmunol.1102003 . PMC 3392547. PMID 22711887 . ^ Huynh; Adler (2012年9月29日). 「エプスタイン・バーウイルス関連感染性単核球症の年齢依存性の数理モデル化」. 数理医学生物学 . 29 (3): 245– 261. doi : 10.1093/imammb/dqr007 . PMID 21700566 . ^ 生物医学工学の理論と実践/生理学的モデリングとシミュレーション 。 ^ クローチェ、ロベルタ;グロンデル、リエンク・ヴァン。アメロンゲン、ハーバート・ヴァン。ストックム、イヴォ・ヴァン。バグショー、クライブ R. (2018 年 9 月 6 日)ダス、ジャヤジット。ジャヤプラカシュ、チリヤム (編)。 システム免疫学: 科学者のための Modelin メソッドの紹介 (第 1 版)。 CRCプレス。 ISBN 978-1-4987-1740-3 。^ Čemerski, Sašo; Das, Jayajit; Locasale, Jason; Arnold, Phoebe; Guirisato, Emanuele; Markiewicz, Mary A.; Fremont, Daved; Allen, Paul M.; Chakraborty, Arup K.; Shaw, Andrey S. (2007年3月26日). 「T細胞抗原 の 刺激効力は免疫シナプスの形成に影響される」 . Immunity . 26 (3): 345– 355. doi : 10.1016/j.immuni.2007.01.013 . PMC 2763191. PMID 17346997 . ^ バグショー、クライヴ R.;バリック、ダグラス E.ノードランド、トーマス M. (2019)。リュブチェンコ、ユーリ・L.クローチェ、ロベルタ。ヴァン・グロンデール、リーンク。ファン・アメロンゲン、ハーバート;ファン・ストックム、イヴォ。ナドー、ジェイ L. (編)。 システム免疫学: 科学者のためのモデリング方法の紹介 。 CRCプレス。 ISBN 978-1-4987-1740-3 。^ バグショー、クライヴ R.;バリック、ダグラス E.ノードランド、トーマス M. (2019)。リュブチェンコ、ユーリ・L.クローチェ、ロベルタ。ヴァン・グロンデール、リーンク。ファン・アメロンゲン、ハーバート;ファン・ストックム、イヴォ。ナドー、ジェイ L. (編)。 システム免疫学: 科学者のためのモデリング方法の紹介 。 CRCプレス。 ISBN 978-1-4987-1740-3 。^ Boianelli, Allessandro; Pettini, Ellena; Prota, Gennaro; Medaglini, Donata; Vicino, Antonio (2015年8月24日). 「一次免疫応答におけるCD4+ T細胞の 増殖 および播種ネットワークの確率モデル」 . PLOS ONE . 10 (8) e0135787. Bibcode : 2015PLoSO..1035787B . doi : 10.1371/journal.pone.0135787 . PMC 4547705. PMID 26301680 . ^ Cagdas Ozturk, Mustafa; Xu, Qian; Cinar, Ali (2018年1月10日). 「1型糖尿病におけるCD8+ T細胞とベータ細胞 の 相互作用に関するエージェントベースモデリング」 . PLOS ONE . 13 (1) e0190349. Bibcode : 2018PLoSO..1390349O . doi : 10.1371/journal.pone.0190349 . PMC 5761894. PMID 29320541 . ^ Tang, Jonathan; Hunt, Anthony C. (2010年2月19日). 「白血球のローリング、活性化、接着を可能にする交戦規則の特定」 . PLOS Computational Biology . 6 (2) e1000681. Bibcode : 2010PLSCB...6E0681T . doi : 10.1371/journal.pcbi.1000681 . PMC 2824748. PMID 20174606 . ^ Kauffman, SA (1969年3月). 「ランダムに構築された遺伝的ネットにおける代謝安定性とエピジェネシス」. Journal of Theoretical Biology . 22 (3): 437– 467. Bibcode : 1969JThBi..22..437K . doi : 10.1016/0022-5193(69)90015-0 . PMID 5803332 . ^ Celada, Franco; Seiden, Philip E. (1996年6月). 「体液性免疫応答のシミュレーションにおける親和性成熟と超変異」. European Journal of Immunology . 26 (6): 1350–8 . doi : 10.1002/eji.1830260626 . PMID 8647216 . S2CID 23503933 . ^ Stewart, JJ; Agosto, H; Litwin, S; Welsh, JD; Shlomchik, M; Weigert, M; Seiden, PE (1997年8月15日). 「リウマトイド因子パラドックスの解決策:病的なリウマトイド因子は、天然のリウマトイド因子との競合によって寛容化される可能性がある」 . Journal of Immunology . 159 (4): 1728–38 . doi : 10.4049 / jimmunol.159.4.1728 . PMID 9257834. S2CID 2443677 . ^ 「DDlab」 。 ^ 「Cell-Devs」 。 ^ Castiglione, Filippo; Celada, Franco (2015年4月7日). 免疫システムとモデリングシミュレーション . CRC Press. ISBN 978-1-4665-9749-5 。^ Keef, Ericka; Zhang, Li Ang; Clermont, Gilles; Swingon, David; Urbano, Alisa; Ermentrout, G. Bard; Matuszewski, Michael; R. Toapanta, Franklin; Ross, Ted M.; Parker, Robert S. (2019年11月). 「インフルエンザウイルス感染の離散的動的モデリングは 、 免疫における年齢依存的な違いを示唆する」 . Journal of Virology . 91 (23) e00395-17. doi : 10.1128/JVI.00395-17 . PMC 5686742. PMID 28904202 . ^ 「ブールネット」 。 ^ 「GINsim」 。 ^ 「セルコレクティブ」 。 ^ Narang, Vipin; Decraene, James; Wong, Shek-Yoon; Aiswarya, Bindu; Wasem, Andrew; Leong, Shiang; Gouaillard, Alexandre (2012年4月13日). 「システム免疫学:モデリング形式論、応用、シミュレーションツールの概説」. 免疫 学的研究 . 53 ( 1–3 ): 251– 265. doi : 10.1007/s12026-012-8305-7 . PMID 22528121. S2CID 18190414 . ^ 「GINsim」 。 ^ 「GINsim」 。 ^ 「セルデザイナー」 。 ^ 「京都遺伝子ゲノム事典」 。 ^ 「リアクトーム」 。 ^ 「GINsim」 。 ^ Rodriguez-Jorge, Otoniel; Kempis-Calanis, Linda A.; Abou-Jaude, Wassim; Gutierrez-Reyna, Darely Y.; 2019年4月16日 (2019). 「CD4+ T細胞活性化におけるT細胞受容体とToll様受容体5シグナル伝達の協調」 ( PDF) . Science Signaling . 12 (577) eaar3641. doi : 10.1126/scisignal.aar3641 . PMID 30992399. S2CID 119514242 . {{cite journal }}: CS1 maint: 数値名: 著者リスト (リンク )^ Müssel, C.; Hopfensitz, M.; Kestler, HA. (2010). 「BoolNet - ブールネットワークの生成、再構築、解析のためのRパッケージ」 . バイオインフォマティクス . 26 (10): 1378– 1380. doi : 10.1093/bioinformatics/btq124 . PMID 20378558 . ^ 「Cranパッケージ - Boolnet」 。2023年2月27日。 ^ 「セルコレクティブ」 。 ^ Helikar, T; Kowal, B; Rogers, JA (2013年2月20日). 「細胞シミュレータプラットフォーム:The Cell Collective」 . Clinical Pharmacology & Therapeutics . 93 (5): 393– 395. doi : 10.1038/clpt.2013.41 . PMC 5242230. PMID 23549147 . ^ ハリス、レナード A.ホッグ、ジャスティン S.タピア、ホセ・フアン。セカール、ジョンAP;グプタ、サンジャナ。コルスンスキー、イリヤ。アローラ、アルシ。バルア、ディパック。シーハン、ロバート P.フェーダー、ジェームス R. (2016 年 11 月 1 日) 「BioNetGen 2.2: ルールベースのモデリングの進歩」 。 バイオインフォマティクス 。 32 (21): 3366–3368 。 土井 : 10.1093/bioinformatics/btw469 。 PMC 5079481 。 PMID 27402907 。 ^ Handel, A. (2020年1月). 「 免疫学者がシミュレーションモデリングを学ぶためのソフトウェアパッケージ」 . BMC Immunol . 21 (1) 1. doi : 10.1186/s12865-019-0321-0 . PMC 6941246. PMID 31898481 . ^ 「ゲパシ」 。 ^ 「Gepasiチュートリアル」 。 ^ 「コパシ」 。 ^ 「BioUML」 。 ^ 「シンバイオロジー」 。 ^ 「バイオSPICE」 。 ^ 「ETHチューリッヒ」 。 ^ 「CSHL」 。 ^ 「会議」 。 ^ Christophe, Christophe; Germain, Ronald; Mathis, Diane (2006年4月). 「『システム免疫学』のためのプラドワイアー」 「免疫学的レビュー . 210 (1). Wiley-Blackwell: 229– 234. doi : 10.1111/j.0105-2896.2006.00374.x . PMID 16623774. S2CID 27355198 .^ Singh, Shailza編(2020年4月25日). Systems and Synthetic Immunology (第1版). Springer Singapore. ISBN 978-981-15-3349-5 。^ Maly, Ivan V. (2009). システム生物学 . Humana Press. ISBN 978-1-934115-64-0 。^ カッツェ、マイケル G. 編(2013年)。 システム生物学 。スプリンガー。 ISBN 978-3-642-33098-8 。^ Houston, MT; Gutierrez, JB (2019年7月13日). 「FRiNDモデル:筋ジストロフィー病態におけるマクロファージ可塑 性 を表現する数理モデル」 . Bull Math Biol . 81 (10): 3976– 3997. doi : 10.1007/s11538-019-00635-8 . PMC 6764940. PMID 31302876 . ^ Mitha, F.; Lucas, TA; Feng, F.; Kepler, TB; Chan, C (2008). 「マルチスケールシステム免疫学プロジェクト:細胞ベースの免疫学的シミュレーションのためのソフトウェア」 . Source Code for Biology and Medicine . 3 6. doi : 10.1186 /1751-0473-3-6 . PMC 2426691. PMID 18442405 . ^ Faeder, JR; Blinov, ML; Hlavacek, WS. (2009). 「BioNetGenを用いた生化学システムのルールベースモデリング」. Systems Biology . Methods Mol Biol. Vol. 500. pp. 113– 167. doi : 10.1007/978-1-59745-525-1_5 . ISBN 978-1-934115-64-0 . PMID 19399430 .