Tech Blog

未知反応を触媒する酵素探索

はじめに

事業開発部の礒崎です。弊社では未知反応を担う酵素探索を行っています。既知酵素反応との反応類似度および既知反応を担う酵素配列との配列相同性を基にして、標的の未知反応を担う酵素配列候補を出力します。本ブログでは、その具体例としてエイズ治療薬候補の1つIslatravirという化合物の合成反応を実際に担う酵素配列を予測した結果をご紹介します。

酵素探索に使用した材料

Huffman et al., 2019のデータを使わせていただきました。この論文ではIslatravirの新規合成経路を設計し、経路中の反応それぞれを触媒する酵素を見つけ、その実証実験を行っていました。Islatravirの合成反応を図1に示しました。化合物6 → 化合物7 or 8 → 化合物5 → 化合物4 → 化合物3 + 化合物2 → Islatravirという順番に合成されます。この合成経路の各反応を担う酵素弊社の酵素探索技術で予測し、論文で使用された酵素と比較しました。

図1. Islatravirの合成経路。Huffman et al., 2019 より抜粋。

結果

まずは、合成経路の5反応それぞれと類似する反応を探索しました。

1. 類似反応探索

出発物質である6→7 (または8→5)酸化反応の類似反応

水酸基をアルデヒド基に酸化する反応が類似反応として多く抽出されました。標的反応との反応類似度が高い類似反応の一部とその順位を図2に示しました。この反応は既知の代謝反応にはないため複数類似反応が抽出されました。

図2. 6→7(8→5)の赤い丸で囲った酸化反応の類似反応探索結果の一部。

出発物質6→8(または7→5)リン酸化反応の類似反応

水酸基をリン酸化する反応が類似反応として多く抽出されました。標的反応との反応類似度が高い類似反応の一部とその順位を図3に示しました。この反応は上記反応同様、既知の代謝反応にはないため複数候補の類似反応が抽出されました。

図3. 6→ 8(7→5)の赤い丸で囲ったリン酸化反応の類似反応探索結果の一部。

中間体5→4 のリボース合成反応

アセトアルデヒドを添加することで閉環してデオキシリボースとなる反応が抽出されました(図4)。論文中でこの反応は既知の代謝反応を模倣しているため、アルキニル基以外全く同じ反応が類似反応として得られました。

図4. 5→4のリボース合成反応の類似反応探索結果。

中間体4→3のリン酸基転移反応

リン酸基をヒドロキシアルキル基から水酸基へ転移する反応が抽出されました(図5)。上記リボース合成反応と同様にこの反応も既知の代謝反応を模倣しているため、アルキニル基以外全く同じ反応が類似反応として得られました。

図5. 4→3のリン酸基転移反応の類似反応探索結果。

中間体3→Islatravirのヌクレオシド合成反応

プリンをデオキシリボースに付加する反応が抽出されました(図6)。上記リン酸基転移反応と同様にこの反応も既知の代謝反応を模倣しているため、アルキニル基およびフッ素以外全く同じ反応が類似反応として得られました。

図6. 3→Islatravirのヌクレオシド合成反応の類似反応探索結果。

2. 類似反応該当酵素探索

結果1で5反応それぞれの類似反応を抽出しました。この類似反応を担う酵素配列taxonごとに抽出しました。絞り込んだ配列の中に論文中で使用された酵素が含まれたかを5反応それぞれで調べました。さらに、5反応それぞれにおいて、全taxon由来の酵素配列から系統的位置を用いたスクリーニングにより配列を絞り込みました。

類似反応のTaxonごとの酵素配列抽出

結果1の類似反応を担う酵素を探索し、Escherichia属、バクテリア、全taxon由来の3段階で抽出しました。抽出した配列数を以下の表1に示しました。今回の論文で使用している酵素と一致するものが含まれていたか調べました。5反応中4反応で、弊社酵素探索技術で抽出した酵素配列に論文中で使用された酵素が含まれていました。

TechBlog2-T1
表1. 類似反応を触媒する酵素のtaxonごとに抽出した配列数、および論文中で使用された酵素の有無。論文中で使用された酵素のUniprotおよびUniParc IDは配列相同性から推定。

系統的位置によるスクリーニング

全taxonより抽出した上記の類似反応酵素配列から、さらに系統的位置により配列を絞り込みます。全配列をクラスタリング・系統樹を生成、系統的にまとまった配列群から1配列ずつ選定しました。この選定の際、その配列の種間での保存度の高いものを優先して選定します(表2、図7)。

TechBlog2-T2
表2 . 全taxonより抽出した類似反応酵素を系統的位置によりスクリーニングした結果。
TechBlog2-7
図7. 3→Islatravir反応の類似反応酵素を系統的位置によって選出した過程。赤のエッジが選出したcentroid配列。さらにこのcentroid配列に属す配列中で最も保存度の高い配列を抽出する。

6→7(8→5)の酸化反応についての考察

上記結果にあるように、論文中でこの反応に使用された酵素は今回見つかりませんでした。その要因としては標的反応6→7(8→5)と論文中で使用した酵素反応がそれほど似ていないことにあると考えます(図8)。今回の類似反応探索で見つけた反応も論文中で使用されたO2を利用したOxidoreductaseであり、標的反応を触媒する可能性があると考えます。図9の反応を触媒する酵素として論文中ではUniParc ID: UPI0001E112C2を使用していると推定されます。これはRHEA_24161を触媒することが確認済みの配列のUniRef50のメンバーに含まれています。UPI0001E112C2自体はこの反応を触媒するかcurationされていません。

図8. 類似反応探索に含まれなかった論文中の酵素反応。論文中の酵素反応は標的反応にあまり類似していない。

終わりに

本ブログでは、未知反応を担う酵素の探索を実演しました。材料として論文で触媒する酵素を見つけたIslatravirの新規合成経路を用いました。この経路の各反応を担う酵素を類似反応酵素探索技術を用いて見つけることができるか試しました。5つの未知反応に対して、実際に類似反応をそれぞれ複数抽出しました。この際、任意のtaxonごとに配列を抽出し、それぞれの配列数を示しました。4つの反応では実際に論文で使っていた酵素を含む複数の配列を抽出することができました。そして、系統的位置を加味したスクリーニングで、全taxonから抽出した類似反応酵素を絞り込みました。通常のスクリーニングではさらに、酵素自体の性質(細胞局在など3次元構造など他の指標も使って候補配列を絞り込めます。

謝辞

今回の類似反応酵素探索には以下の論文のデータを利用させていただきました。

Huffman et al., (2019) Design of an in vitro biocatalytic cascade for the manufacture of islatravir. Science.


Page top