A polarized charge model for RNA: enhanced stability of G- quadruplexes and improved folding free energy landscapes

地球上の生命は、リボ核酸、すなわちRNAに書かれた複雑な分子言語に依存しています。DNAはしばしば遺伝のマスター設計図として称賛されますが、RNAは遺伝的指示をタンパク質へと翻訳し、無数の生物学的プロセスを制御する「能動的な作業者」として機能します。RNAがどのようにこれらのタスクを遂行するのかを理解するために、科学者たちは、RNAがどのように特定の三次元形状へと折り畳まれるのかを見極めなければなりません。分子の形状はその機能を決定します。それは、鍵が錠前にはまるためには正しい輪郭を持っていなければならないのと似ています。数十年にわたり、研究者たちは強力なコンピュータシミュレーションを用いて、これらの分子が動き、折り畳まれる様子を観察し、実験室で実際に構築することなくその挙動を予測しようと試みてきました。しかし、これらのデジタルモデルはRNAの真の性質を捉えることに苦戦しており、しばしば実際の実験で見られるものよりも緩すぎたり、不安定だったりする形状を生み出してしまいます。

この問題の核心は、シミュレーションが電気をどのように扱うかという点にあります。RNA分子は高度に帯電しており、水分子が絶えず電気的な力を遮蔽し、形を変え続ける水環境の中に存在しています。標準的なコンピュータモデルは、伝統的に、すべての原子の電荷を固定された不変の値として扱っており、周囲の水や他の原子がその電荷をわずかに変化させる可能性を無視してきました。この簡略化により、シミュレーションでは異なる部分同士がどの程度引き付け合い、あるいは反発し合うかの計算を誤ることになります。その結果、デジタル上の分子はバラバラになったり、自然界が作り出すような緻密で複雑な構造を形成できなかったりすることがあるのです。

最近の研究において、済寧大学と深セン先端技術大学の研究チームは、この問題を解決する新しい方法を提案しました。彼らは、RNA原子の電荷を固定された数値としてではなく、分子の特定の形状や環境に応じて変化する値として計算する方法を開発しました。彼らはこのアプローチを「分極電荷モデル」と呼んでいます。この新手法を2つの異なるタイプのRNA構造に適用することで、研究チームは、従来のモデルよりもはるかに効果的にこれらの分子を結合させられることを実証しました。この結果は、変化する電気的な力を考慮に入れることが、RNAの正確な「デジタルツイン」を作成するために不可欠であることを示唆しています。

研究チームはまず、Gクアドルプレックス（G-quadruplex）として知られる構造でモデルのテストを行いました。これらは、グアニンと呼ばれる構成要素に富むRNA配列によって形成される、複雑な4本鎖の形状です。Gクアドルプレックスは、グアニンの平らな環の積み重ねによって安定化されており、極めて重要な点として、構造を保持するために中心のチャネルにカリウムやナトリウムなどの金属イオンが位置している必要があります。従来のモデルを用いた標準的なシミュレーションでは、これらの中心イオンは周囲の水へと漂い去ってしまうことが多く、その結果、構造全体が解けてしまいます。研究者がこの新しい分極電荷モデルを適用したところ、イオンはしっかりと所定の位置に留まりました。モデルは、中心のイオンと周囲のRNA環との間の強い電気的引力を再現することに成功し、シロレーションの間、構造を安定させ続けました。この安定性は、異なる種類のイオンや様々な標準的なシミュレーション設定においても観察され、この新手法が、従来の方法で計算されていた相互作用における根本的な欠陥を修正していることを示しています。

次に、チームはより単純ですが同様に重要な構造である、RNAヘアピンに注目しました。これは、単一のRNA鎖が自身に折り返してステム（茎）とループ（環）を形成する、一般的な折り畳み構造です。モデルをテストするため、研究者は、緩やかなRNA鎖がこのヘアピン形状へと折り畳まれるプロセスをシミュレートしました。彼らは、コンピュータが多くの異なる形状を探索し、それぞれの形状に到達するために必要なエネルギーを計算できる特殊な手法を用いました。標準的なモデルを用いたシミュレーションでは、エネルギー地形は乱雑で混乱しており、折り畳まれたヘアピンが最も安定した状態であることを示すことに失敗することがよくありました。対照的に、研究者がこの分極電荷モデルを使用したところ、シミュレーションは一貫して、折り畳まれたヘアピンが最もエネルギー的に有利な形状であることを示しました。モデルは、仮想のRNA鎖を正しいコンパクトな形態へと導くことに成功し、実世界の実験で観察されるものと一致しました。

この研究の重要な洞察は、原子の電気的電荷は静的な特性ではなく、その周囲の環境に依存するということです。折り畳まれたヘアピンにおいては、水素結合に関与する原子は、緩やかな未折り畳みの鎖の状態よりも、わずかに電荷が強くなります。新しいモデルはこの変化を捉えますが、古いモデルは、原子が緻密な結び目の内部に埋もれていようと、水にさらされていようと、すべての原子が同じ電荷を持っているものとして扱います。この感度を取り入れることで、研究者たちは、シミュレーションソフトウェアの基礎となる数学的規則を変更することなく、分子の正しい安定性を再現できることを見出しました。これは、この手法が標準的なツールを使用している他の科学者たちによって容易に採用できることを意味しており、分野への実用的なアップグレードを提供しています。

また、この研究は一つの限界についても明らかにしました。新しい電荷の値は、その電荷を計算するために使用された特定の形状に依存するという点です。もし研究者が、緩やかな未折り畳みの鎖に基づいて電荷を計算し、その後に折り畳みプロセスをシミュレートしようとした場合、結果の精度は低下しました。モデルは、電荷が最終的な安定構造に似た形状から導き出された場合に最もよく機能します。このことは、最も正確な結果を得るためには、科学者が分子の振る舞いをシミュレートする前に、既知の安定した形状を用いて電荷パラメータを生成すべきであることを示唆しています。

結局のところ、この研究はRNAの微視的な世界へのより明確な窓を提供しています。電気的な力を計算する方法を洗練させることで、チームは、現実に対してはるかに忠実なシミュレーションを作成することが可能であることを示しました。Gクアドルプレックスを安定させ、ヘアピンの折り畳みを正しく予測できる能力は、健康と疾患におけるRNAの機能を理解しようとする科学者たちを、より真実に近づけます。この手法はあらゆるシナリオに対する完璧な解決策ではありませんが、原子レベルで生命の動的かつ複雑な機構をモデル化しようとする探求における、重要な前進となっています。