[論文レビュー] Numerical Techniques for Applications of Analytical Theories to Sequence-Dependent Phase Separations of Intrinsically Disordered Proteins

Yi‐Hsuan Lin, Jonas Wessén+3 more|arXiv (Cornell University)|Jan 6, 2022

RNA Research and Splicing被引用数 5

ひとこと要約

本稿では、自己無秩序タンパク質（IDP）における配列依存性液体-液体相分離（LLPS）に解析理論を適用する実用的な計算フレームワークを提示する。Flory-Huggins理論およびランダム位相近似（rG-RPA）理論を用いて、スピンodalおよびビノーダル相境界を数値的に計算する手法を詳述し、複素ランジュバンダイナミクスを用いた場理論的シミュレーション（FTS）を可能にするとともに、粗粒度分子動力学（CG-MD）のためのスラブプロトコルを提示し、抽象的な理論と明示的シミュレーションの間をつなぐものとなる。

ABSTRACT

Biomolecular condensates, physically underpinned to a significant extent by liquid-liquid phase separation (LLPS), are now widely recognized by numerous experimental studies to be of fundamental biological, biomedical, and biophysical importance. In the face of experimental discoveries, analytical formulations emerged as a powerful yet tractable tool in recent theoretical investigations of the role of LLPS in the assembly and dissociation of these condensates. The pertinent LLPS often involves, though not exclusively, intrinsically disordered proteins engaging in multivalent interactions that are governed by their amino acid sequences. For researchers interested in applying these theoretical methods, here we provide a practical guide to a set of computational techniques devised for extracting sequence-dependent LLPS properties from analytical formulations. The numerical procedures covered include those for the determinination of spinodal and binodal phase boundaries from a general free energy function with examples based on the random phase approximation in polymer theory, construction of tie lines for multiple-component LLPS, and field-theoretic simulation of multiple-chain heteropolymeric systems using complex Langevin dynamics. Since a more accurate physical picture often requires comparing analytical theory against explicit-chain model predictions, a commonly utilized methodology for coarse-grained molecular dynamics simulations of sequence-specific LLPS is also briefly outlined.

研究の動機と目的

• 研究者に、解析理論から得られる配列依存性LLPS特性を実装可能な数値的手法で抽出するための技術を提供すること。
• 抽象的な理論的定式化と生物学的条件下の現実的シミュレーションとの間のギャップを埋めること。
• アミノ酸配列に基づいて、多価IDP系のスピンodalおよびビノーダル相境界を正確に予測できること。
• 『スラブ』法を用いた、配列特異的LLPSのための粗粒度分子動力学シミュレーションの体系的プロトコルを提供すること。
• 簡略化モデル系において、解析的予測（例：Flory-Huggins、rG-RPA）と明示的鎖シミュレーションの結果を比較し、理論モデルの妥当性を検証すること。

提案手法

• Flory-Huggins（FH）およびランダム位相近似（rG-RPA）理論から導出された自由エネルギー関数を用いて、スピンodalおよびビノーダル相境界を数値的に決定する。
• スピンodal境界を特定するための曲率の変化点（f''(ρ) = 0）の計算と、ビノーダル共存濃度を特定するための共通接線法の適用。
• 配列特異的相互作用を持つヘテロポリマー系をモデル化するため、複素ランジュバルンダイナミクスを用いた場理論的シミュレーション（FTS）の適用。
• HOOMD-blueにおける『スラブ』プロトコルの実装：Packmolによる初期構造生成、エネルギー最適化、NPT圧縮、z軸方向へのボックス拡張を含む。
• 多成分系における共存濃度を特定するためのタイライン構築法の適用。
• 簡略化モデル系における相挙動の比較的分析を通じて、解析理論と明示的シミュレーションを統合すること。

実験結果

リサーチクエスチョン

• RQ1 一般の自由エネルギー関数から、配列依存性LLPSにおけるスピンodalおよびビノーダル相境界をどのように数値的に計算できるか？
• RQ2 Flory-Huggins理論およびrG-RPA理論は、特定の配列パターンを持つIDPの相分離挙動を予測するために果たす役割は何か？
• RQ3 複素ランジュバルンダイナミクスを用いた場理論的シミュレーションは、平均場理論の近似を避けて、配列特異的LLPSをどのように正確に捉えられるか？
• RQ4 『スラブ』法を用いた粗粒度明示的鎖モデルにおけるバイオ分子LLPSをシミュレートする最適なプロトコルは何か？
• RQ5 モデル系において、解析的理論の予測と明示的鎖シミュレーションの結果を定量的に比較すると、どのような相違・一致が得られるか？

主な発見

• スピンodal境界は、自由エネルギー曲線f(ρ)の曲率変化点（f''(ρ) = 0）を特定することで、ビノーダル境界は共通接線法を用いて、異なる温度における相境界を成功裏に構築した。
• N=3、χ=1.7のFlory-Hugginsモデルでは、f''(ρ) = 0となるρ値においてスピンダル領域が特定され、共通接線法によりビノーダル共存濃度ρ^αおよびρ^βが決定された。
• 複素ランジュバルンダイナミクスを用いた場理論的シミュレーションにより、ヘテロポリマー系における配列依存性LLPSが高精度にシミュレートされ、平均場理論を逸脱する効果が捉えられた。
• HOOMD-blueにおける『スラブ』プロトコルは、密度プロファイルおよび異なる鎖種の空間的相分離の観察により、相分離を示すドロイプを有する平衡系を成功裏に生成した。
• sv1およびsv28鎖のシミュレーションでは、凝縮相において両鎖が大きく相分離しており、配列特異的分離が確認された。
• 双極子溶媒分子を含むsv30鎖のシミュレーションでは、ドロイプ内でのポリマー鎖と溶媒双極子の共局在が観察され、溶媒が相分離に参加していることが示された。