Pi and H2 PO4- are kinetically equivalent as fatigue factors in skeletal muscle

技術的要約：骨格筋における疲労因子として、PiとH2PO4-は速度論的に等価である

問題提起
筋疲労は、運動中に起こる可逆的なパワーまたは力の低下と定義され、多くの場合、実際の出力が利用可能なパワーと一致したときに終了する。無機リン酸（Pi）およびプロトン（H+）は、主要な疲労関連代謝物として広く考えられているが、特定のメカニズム的な「真の」疲労因子をめぐっては議論が存在する。一部のエビデンスは、二プロトン化された形態の無機リン酸（Pi-、H2PO4-と表記）こそが、スキン化線維において力を抑制する特定の因子であることを示唆している。なぜなら、その濃度は全Pi量とpHの両方に依存するためである。これに対し、バイオエナジェティック・モデルでは、全Piがしばしば主要な駆動因子として扱われる。本研究で取り組む中心的な問題は、Pi-に基づく疲労メカニズムが、全Piに基づくメカニズムによって生成される動的なシステム特性を速度論的に複製できるのか、あるいはその区別が単にメカニズム上の違いに過ぎず、速度論的な違いではないのかという点である。

方法論
本研究では、コンピュータによる骨格筋バイオエナジェティック・システムの動的モデルを用いて、in silico（コンピュータ内）実験を行う。モデルには以下が含まれる：

• バイオエナジェティック構成要素： ATP使用量（アクトミオシンATPase、Ca2+-ATPase、安静時使用量）、酸化リン酸化（OXPHOS）、NADH供給（TCAサイクル、脂肪酸酸化）、解糖系、クレアチンキナーゼ、およびプロトンフラックス。
• 刺激メカニズム： 「各ステップ活性化（each-step-activation: ESA）」メカニズム。ここでは、安静から運動への移行期において、OXPHOS、NADH供給、および解糖系がATP使用量と並行して直接活性化される。
• 疲労メカニズム： 以下の2つのバージョンのモデルを比較した：
  1. PiDT（オリジナル）： Pi-ダブルスレッショルド（二重閾値）メカニズムに基づく。追加のATP使用（V̇O2スローコンポーネントを駆動する）は、全Piが臨界閾値（$P_{icrit} = 18$ mM）を超えたときに開始される。運動は、Piがピーク閾値（$P_{ipeak} = 25$ mM）に達したときに終了する。
  2. Pi-DT（修正版）： オリジナルと同一であるが、疲労因子としての全Piを二プロトン化リン酸（Pi-）に置き換えたもの。閾値は平衡関係（$P_i^- = P_i / (1 + 10^{pH - pKa})$、ただし $pKa = 6.8$）に基づいて調整された。その結果、$P_{icrit}^- = 8$ mM および $P_{ipeak}^- = 13$ mM となった。

シミュレーションは、中強度、高強度、および過酷（非常に高強度）の強度領域における定常パワー運動（CPE）に対して実施された。モデルは、V̇O2、代謝物（PCr、Pi、ADP、pH）のタイムコース、および限界時間（time to exhaustion）を生成した。また、クリティカルパワー（CP）およびW'（CPを超える仕事容量）を決定するために、パワー・デュレーション関係もシミュレートされた。

主な貢献と結果
本研究は、Pi-DTメカニズムが、コンピュータモデルの文脈においてPiDTメカニズムと速度論的に等価であることを示している。主な知見は以下の通りである：

• 動的キネティクス： Pi-DTバージョンは、安静から運動への移行期におけるV̇O2および代謝物のタイムコースを、PiDTバージョンと実質的に同一のものとして生成した。両モデルは、それぞれの臨界閾値を超えた際に現れるV̇O2スローコンポーネントを含め、中強度、高強度、および過酷な運動強度領域を正常に再現した。
• 運動の終了： 過酷な運動において、両モデルは、それぞれの疲労因子がピーク閾値に達したことによる疲労による運動終了を予測した。運動強度の増加に伴い、両バージョンにおいて限界時間は減少した。
• パワー・デュレーション関係： 両メカニズムは双曲線的なパワー・デュレーション関係を生み出した。パラメータはほぼ同一であった：
  • PiDT： 臨界ATP使用量（$AUT_{crit}$）= 75、W' = 53 任意単位。
  • Pi-DT： $AUT_{crit}$ = 73、W' = 53 任意単位。
• 軽微な相違： 運動終了時のADP、pH、および正確な限界時間の値において、無視できる程度の差異のみが観察された。著者はこれを、根本的な速度論的差異ではなく、選択された特定の閾値の値に起因するものとしている。

意義と主張
本論文は、このバイオエナジェティック・モデルの枠組み内において、PiとPi-は疲労因子として速度論的に等価であるが、メカニズム的には依然として異なるものであると結論付けている。

• 「スーパー代謝物」の概念： モデルは、速度論的な出力のみに基づいてPiとPi-の作用を区別できないため、著者はPi、Pi-、およびH+を、単一の疲労関連「スーパー代謝物」として一括して扱うことができると提案している。それらの構成成分は、pH依存的な平衡により、運動移行期において一意に関連している。
• モデリングの限界： 本研究は、モデルはこれらの代謝物の影響をシミュレートすることはできるが、アクトミオシン阻害やカルシウム沈殿といった疲労プロセスに対する特定の形態（Pi vs. Pi- vs. H+）のメカニズム的な寄与を解決することはできないと主張している。
• 実験の必要性： 著者は、特に（代謝物を独立して操作できる）in vitro条件と（それらが並行して変化する）in vivo条件を区別する実験研究が、真のメカニズム的疲労因子を特定するために必要であると強調している。本論文は、メカニズムの論争を解決することを目的としたものではなく、速度論的なモデリングの観点からは、主要な変数としてPiとPi-のどちらを選択しても、システム力学は区別不能であることを確立することを目指している。

本研究は、Pi-ダブルスレッショルドの枠組みが、トリガーが全Piであってもその二プロトン化形態であっても、筋疲労の複雑な動的特性（スローコンポーネント、強度領域、パワー・デュレーション曲線）を説明する上で有用であることを裏付けている。