Interaction between cell membranes and protein inclusions in the large-deformation regime

技術要約：大変形領域における細胞膜とタンパク質包含物の相互作用

問題提起
生物学的膜は動的な表面であり、その形態と機能はタンパク質包含物（PI）によって決定的に支配されている。PIと脂質膜の相互作用については、平坦な構成の周りでの摂動展開を用いた小変形（SD）領域において広範にモデル化されてきたが、多くの生物学的に重要なプロセスは、強い膜変形を伴う。これには、PIが高速度で移動して顕著な陥入を引き起こすシナリオや、複数のPIが近接して配置され、線形理論では記述できない急峻な勾配が生じるシナリオが含まれる。本研究は、SD近似の破綻と複雑な弾性挙動の出現に焦点を当て、大変形（LD）領域における膜とタンパク質の相互作用に関する理解の空白を埋めるものである。

手法
著者らは、有限要素（FE）数値シミュレーションと近似的な解析的導出を組み合わせたハイブリッドアプローチを採用している。

• 数値フレームワーク： 本研究では、ヘルフィヒヒ自由エネルギーから導出された完全な非線形形状方程式を解くために、有限要素ソフトウェアIRENEを利用している。膜は、モンジュ・ゲージ（$z(x)$）を用いた3次元に埋め込まれた2次元多様体としてモデル化されている。支配方程式には、平均曲率（$H$）、ガウス曲率（$K$）、曲げ剛性（$\kappa$）、および表面張力（$\sigma$）の項が含まれる。
• 境界条件： シミュレーションでは、半径$R$における固定された膜境界、タンパク質界面（無限に硬いものとしてモデル化）における固定された接触角（$\alpha$）、および膜流のノンスリップ条件を含む特定の制約を課している。
• 解析的手法： 物理的な洞察を得るために、著者らは平均曲率がゼロ（$H=0$）であると仮定することで、LD領域に対する近似的な解析解を導出している。これはカタノイド（懸垂面）のプロファイルをもたらし、接触角と変位の関数としての膜の高さと力の明示的な表現を可能にする。
• 流動ダイナミクス： 膜流が存在する場合、本研究では質量保存、運動量バランス（対流項と粘性応力を含む）、および形状方程式を記述する結合偏微分方程式（PDE）を解いている。ただし、計算領域内でのストークスのパラドックスを避けるため、外部流体は無視している。

主要な貢献と結果

1. 線形理論の限界：
   FE解を線形化方程式（勾配が小さい場合に有効）と比較することにより、著者らは、接触角が $\tan \alpha \approx 0.4$ を超える場合、または垂直変位が大きい場合に、線形近似が成立しなくなることを示している。これらのLD領域においては、正しい膜の幾何学構造を捉えるために、完全な非線形形状方程式を解く必要がある。

2. 近似的解析解（平均曲率ゼロ）：
   著者らは、平均曲率がゼロであるという仮定に基づく解析解を提案しており、これはカタノイドのプロファイルを与える。これは一般的なケースに対する近似ではあるが、LD領域の定性的な特徴をうまく捉えている。これは、固定された接触角を持つタンパク質の自発的な垂直変位（$h$）が、領域サイズに対して対数的にスケールすること（$h \sim \ln R$）を予測している。

3. 力と変位の関係：
   本研究では、膜がPIに及ぼす力を計算している。重要な知見は、この力がタンパク質の垂直変位に対して非単調な挙動を示すことである。力は極大値まで増加した後、減少する。この特徴は特定のタンパク質の幾何学形状に依存しないように見える。この非単調性は、チュビュール（小管）形成や退縮に関連する動的な不安定性の発生を示唆している。力の大きさは1–10 pNの範囲と推定されており、これは分子モーターによって及ぼされる力と一致している。

4. 2つのPI間の膜媒介相互作用：
   2つの相互作用するPIの系において、膜媒介ポテンシャルエネルギーは、SD領域で見られる単純なべき乗則の減衰に従わない。代わりに、ポテンシャルは劣べき乗則（sub-power-law）の減衰を示し、タンパク質間距離に対してべき関数よりも緩やかに減少する。これは、複雑な弾性媒体の性質を反映している。さらに、相互作用ポテンシャルはタンパクスの向きの奇関数である：

   • 同一の向き（同じ「電荷」）を持つPI同士は反発する。
   • 反対の向き（反対の「電荷」）を持つPI同士は引き合う。
     これは、LD領域におけるPIの向きと電気的な電荷との間の類似性を裏付けている。系は反対称な構成をエネルギー的に好むため、タンパク質が交互に並ぶ反強磁性的パターンを生じる可能性がある。

5. 流動誘起変形と特性速度：
   膜流が存在する場合、著者らは特性速度スケール $v^* = \kappa / (L\eta)$ を特定しており、これはScriven-Love数が1となることに対応する。この速度は、以下の2つの異なる領域を分ける：

   • $v < v^*$：流動が膜の形状に与える影響は無視でき、曲げ剛性が支配的である。
   • $v > v^*$：流動誘起の変形が現れ、その特性波長は $\lambda \sim \kappa/(v\eta)$ となる。
     本研究は、この閾値を超えると、粘性力と曲げ力の競合により顕著な形状変化が生じることを示唆している。

意義と主張
本論文は、摂動的なSD理論の限界を超え、生物学的に関連のあるシナリオにおける膜-タンパク質系の定量的予測を提供することを目的としている。結果は以下の事項に直接的な示唆を与える：

• タンパク質のソーティングとクラスタリング： 非単調な力と、向きに基づく特定の引力・斥力の規則は、タンパク質がどのように膜上でクラスターやパターンを形成するかというメカニズムを示唆している。
• 膜輸送： 流動誘起変形の臨界速度の特定は、タンパク質の運動が能動的に膜の形状変化を駆動し、小胞形成や輸送に影響を与える可能性を示している。
• 生物学的例： 著者らは、巨大単層小胞（GUV）上のバクテリオロドプシン（BR）を用いて、自身の知見の妥当性を示している。BR分子の横方向拡散速度は、流動駆動型の変形を誘起するのに十分な大きさであり、その波長はタンパク質間の間隔と同程度であることから、流動媒介型の相互作用がBRのパターニングに大きく影響する可能性があることを示唆している。

総じて、本研究は、膜の弾性、タンパク質の幾何学、および流体力学的な流れの間の複雑で非線形な相互作用を分析するための枠組みを確立しており、動的な細胞環境における膜力学のより現実的な記述を提供している。