[論文レビュー] Quantitative modeling and data analysis of SELEX experiments

Marko Djordjević, Anirvan M. Sengupta|arXiv (Cornell University)|Nov 30, 2005

Advanced biosensing and bioanalysis techniques参考文献 39被引用数 10

ひとこと要約

本稿では、転写因子-DNA 相互作用のパラメータ推定の正確性を向上させるために、生物理学的根拠に基づいた定量的モデルを SELEX 実験に提案する。化学ポテンシャルをラウンド間で固定するように SELEX プロトコルを変更し、有限温度エネルギー行列法を適用することで、従来の情報理論的および経験的手法と比較して、誤検出と誤検出の両方を顕著に低減した。CFT/NFI データを用いた実証により、結合部位の予測性能が向上したことが示された。

ABSTRACT

SELEX (Systematic Evolution of Ligands by Exponential Enrichment) is an experimental procedure that allows extracting, from an initially random pool of DNA, those oligomers with high affinity for a given DNA-binding protein. We address what is a suitable experimental and computational procedure to infer parameters of transcription factor-DNA interaction from SELEX experiments. To answer this, we use a biophysical model of transcription factor-DNA interactions to quantitatively model SELEX. We show that a standard procedure is unsuitable for obtaining accurate interaction parameters. However, we theoretically show that a modified experiment in which chemical potential is fixed through different rounds of the experiment allows robust generation of an appropriate data set. Based on our quantitative model, we propose a novel bioinformatic method of data analysis for such modified experiment and apply it to extract the interaction parameters for a mammalian transcription factor CTF/NFI. From a practical point of view, our method results in a significantly improved false positive/false negative trade-off, as compared to both the standard information theory based method and a widely used empirically formulated procedure.

研究の動機と目的

• 転写因子結合部位を予測するための従来の位置加重行列（PWMs）の感受性および特異性の低さを是正すること。
• SELEX 実験手順の標準的手法に内在する欠陥を同定し、それらが転写因子-DNA 相互作用のパラメータ推定を損なう要因となることの解明。
• 生物理学的相互作用パラメータの堅牢な推定を可能にする、改良された SELEX プロトコルおよび新しい計算手法の開発。
• 提案手法から導出されたエネルギー行列が、情報理論的および経験的 PWM と比較して、結合部位予測性能において優れていることを実証すること。

提案手法

• 結合自由エネルギーと統計力学を組み込んだ、転写因子-DNA 結合の生物理学的モデルを用いて、SELEX 実験を定量的にシミュレートする。
• 自己無撞着方程式（式 47）を解いて得られる有限温度エネルギー行列（FTEM）を導入し、結合の飽和および熱力学的一致性を反映する。
• 化学ポテンシャルをラウンド間で一定に保つことで、パラメータ推定に適したデータセットを生成する、改良された SELEX プロトコルを提案する。
• 熱力学的制約下で最大尤度に類似したアプローチを用いて、実験的 SELEX データにモデルを適合させることで、エネルギー行列のパラメータを推定する。
• 二量体転写因子（例：CTF/NFI）に対して、結合部位行列に回文対称性を強制することで、対称性制約を組み込む。
• 検出誤りトレードオフ（DET）曲線を用いて、得られたエネルギー行列の性能を評価し、情報理論的および MatInspector 行列と比較する。

実験結果

リサーチクエスチョン

• RQ1 標準的な SELEX プrotocol は、転写因子-DNA 相互作用のパラメータ推定に適したデータを信頼性高く生成できるか？
• RQ2 結合の飽和における熱力学的一致性が、SELEX データにおける従来の情報理論的加重行列構築法の信頼性にどのように影響するか？
• RQ3 SELEX データからの堅牢かつ再現可能なパラメータ推定を保証するためには、どのような実験的修正が必要か？
• RQ4 提案手法の有限温度エネルギー行列法は、従来の手法と比較して誤検出と誤検出のトレードオフにおいてどのように優れているか？
• RQ5 個々の行列要素の差が、結合部位予測性能の差にどの程度寄与するか？

主な発見

• 標準的な SELEX 手順では、ラウンド間で化学ポテンシャルが一貫しないため、エネルギー行列の正確な推定に適したデータが得られない。
• 化学ポテンシャルを固定した改良された SELEX プロトコルにより、堅牢なデータ生成と信頼性の高いパラメータ推定が可能になる。
• 有限温度エネルギー行列法は、情報理論的および経験的 PWM と比較して、誤検出/誤検出トレードオフを顕著に低減する。
• 位置 2G において、有限温度エネルギー行列は情報理論的行列と約 150% の差を示し、位置 4C および 5A では約 100% の差を示しており、顕著なパラメータ差が確認された。
• 個々の要素の差はやや小さいものの、累積的効果により予測性能の向上が顕著に現れ、DET 曲線の比較で裏付けられた。
• Roulet らのデータから導出した CTF/NFI のエネルギー行列は、保存された位置（2,3,4）で MatInspector と強い一致を示し、保存度が低い位置（1,5,6）ではより大きな差を示しており、機能的影響と整合的である。