黃芪活性成分結合EphA2的分子對接分析

韓 雪，金家莉，周蜜兒，武體林，萬子豪，李志強，卞冰芝 （1.蚌埠醫學院：a.檢驗醫學院免疫學教研室；b.感染與免疫安徽省重點實驗室；c.慢性疾病免疫學基礎與臨床安徽省重點實驗室；d.檢驗醫學院；e.臨床醫學院，安徽 蚌埠 33030；.安徽醫科大學附屬六安醫院重癥醫學科，安徽 六安 37005）

在藥物研發領域，尋找有效的藥物分子并將其與目標蛋白結合是一個極具挑戰性的任務。分子對接技術是一種常用的計算機輔助藥物設計方法[1-3]，它可以預測分子之間的相互作用，為藥物研發提供有力支持。黃芪作為中藥材，具有廣泛的藥理活性，其主要活性成分黃芪總苷、黃芪總黃酮和黃芪總糖等化合物已被證明具有多種生物活性，具有很高的應用前景[4]。本研究旨在通過分析黃芪主要活性成分與目標蛋白EphA2 的結合模式，探討藥物分子與目標蛋白之間相互作用的機制，為藥物研發提供新的思路和方法。

1 材料和方法

1.1 分子準備

黃芪主要活性成分一般為三大類，分別為黃芪總糖、黃芪總苷和黃芪總黃酮。黃芪總糖主要成分為葡萄糖、阿拉伯糖、半乳糖、甘露糖、果糖、木糖、葡糖醛酸和半乳糖醛酸等，這些單體主要通過α-型糖苷鍵進行連接，β-型糖苷鍵較少。黃芪總苷主要成分為黃芪皂苷甲，是一種羊毛酯醇形的四環三萜皂苷。黃芪總黃酮主要成分為黃芪苷和毛蕊異黃酮苷，兩者的含量在黃芪總黃酮中較高，可達到30%~40%左右。此外，黃芪總黃酮中還含有其他黃酮苷類化合物，如木犀草苷、異黃酮苷、甘草苷等，它們的含量也可能在5%~20%之間。本實驗選取了阿拉伯糖、黃芪皂苷甲和毛蕊異黃酮苷作為黃芪的主要活性成分進行后續研究。

1.2 蛋白質準備

從PDB 數據庫獲得蛋白結構[5]，在gHgL 蛋白與EphA2 蛋白結合的結構7CZF 中提取出EphA2 蛋白結構用于分子對接。再使用pymol和autodocktools軟件處理蛋白質分子。從兩個相同的蛋白質單體EphA2 結合形成的二聚體分子中提取出一條EphA2蛋白，移除結構中的所有水分子，并對蛋白進行加氫處理。

1.3 分子對接方法

使用Autodock4 和Autodock vina 兩種不同的分子對接方法，Autodock4（版本號4.2.6，美國）和Autodock vina（版本號1.2.3，美國）軟件均為蚌埠醫學院購買。在對接過程中將搜索參數設置為60 ?×60 ?×60 ? 的網格框，中心位于蛋白的活性位點。網格間距設置為0.375 ?，這是Autodock4 中的默認值。關于蛋白的活性位點，實驗分析了gHgL 與EphA2 蛋白的結合結構，重點關注兩者結合的區域，確定參與結合的重要的EphA2 蛋白位點部分。此外，還允許蛋白側鏈的靈活性參與對接過程，這是準確預測小分子結合的重要特征。

2 結果

2.1 分子對接結果

2.1.1 EphA2與阿拉伯糖對接結果

圖1所示是使用Autodock4對接得到的最穩定的結構，結合能為-4.02 kcal/mol。對接的主要位點為阿拉伯糖的羥基與EphA2 蛋白的殘基SER74、VAL72、MET73、GLN77、LYS50，一共形成了6 條氫鍵。結合中起主要貢獻的是范德華作用、氫鍵以及脫溶能，三者一共貢獻了-4.81 kcal/mol 的能量，靜電能的貢獻較小，僅為-0.04 kcal/mol。圖2所示是使用Autodock vina進行對接得到的最穩定的結構，結合能為-3.71 kcal/mol。阿拉伯糖一共與EphA2 形成了1條氫鍵，為阿拉伯糖的羥基氧與氨基酸的氫原子成鍵。表1 所示的是使用Autodock4 對接得到的10 個構象及其結合能、氫鍵數目以及氫鍵連接的殘基，涉及的主要殘基是SER74、ASP76、ASP78、GLY51、ASN71、ASP53、ASN71、THR45、VAL102 等。表2 是使用Autodock vina 進行對接得到的10個構象及其及核能、氫鍵數目以及氫鍵連接的殘基，涉及的殘基有GLN77、GLU117、VAL102、ASN164、THR45、TYR67、ASP53、ARG82、ASP129、ASN60、ASP61、LEU128、ARG168、PHE134 等。綜合兩種對接方法，阿拉伯糖與EphA2 結合的殘基主要有GLN77、VAL102、THR45、TYR67、ASP53、ASN164。

表1 使用Autodock4對接方法所得的所有構象結合能以及結合的殘基Tab.1 All conformational binding energies and bound residues obtained using the Autodock4 docking method

表2 使用Autodock vina 對接方法所得的所有構象結合能以及結合的殘基Tab.2 All conformational binding energies and bound residues obtained using the Autodock vina docking method

圖1 EphA2與阿拉伯糖Autodock4對接結果Fig.1 Docking results of EphA2 and arabinose by Autodock4

圖2 EphA2與阿拉伯糖Autodock vina對接結果Fig.2 Docking results of EphA2 and arabinose by Autodock vina

2.1.2 EphA2與黃芪皂苷甲對接結果

圖3所示是使用Autodock4對接得到的最穩定的結構，結合能為-5.06 kcal/mol。對接的主要位點為黃芪皂苷甲羥基與EphA2 蛋白的殘基ASP53、ASN71、VAL72、MET73、ARG103，一共形成了7 條氫鍵。結合中起主要貢獻的是范德華作用、氫鍵以及脫溶能。圖4 所示是使用Autodock vina 進行對接得到的最穩定的結構，結合能為-3.85 kcal/mol。黃芪皂苷甲一共與EphA2形成了1條氫鍵，為黃芪皂苷甲的羥基氧與氨基酸的氫原子成鍵。表3 所示的是使用Autodock4對接得到的10個構象及其結合能、氫鍵數目以及氫鍵連接的殘基，涉及的主要殘基是ASP53、ASN71、ARG103、MET73、VAL72、LYS162、THR101等。表4是使用Autodock vina進行對接得到的10 個構象及其及核能、氫鍵數目以及氫鍵連接的殘基，涉及的主要殘基有ARG137、GLY49、ASP76、LYS50、THR45、GLN56等。綜合兩種對接方法，阿拉伯糖與EphA2 結合的殘基主要有ASN71、ARG103、THR101。

表3 使用Autodock4 對接方法所得的所有構象結合能以及結合的殘基Tab.3 All conformational binding energies and bound residues obtained using the Autodock4 docking method

表4 使用Autodock vina 對接方法所得的所有構象結合能以及結合的殘基Tab.4 All conformational binding energies and bound residues obtained using the Autodock vina docking method

圖3 EphA2與黃芪皂苷甲Autodock4對接結果Fig.3 Docking results of EphA2 and Astragalus saponin A Autodock4

圖4 EphA2與黃芪皂苷甲Autodock vina對接結果Fig.4 Docking results of EphA2 and Astragalus saponin A Autodock vina

2.1.3 EphA2與毛蕊異黃酮苷對接結果

圖5所示是使用Autodock4對接得到的最穩定的結構，結合能為-5.87 kcal/mol。對接的主要位點為毛蕊異黃酮苷羥基與EphA2 蛋白的殘基LYS50、ARG103、MET73、GLN77、VAL72，一共形成了5 條氫鍵，結合中起主要貢獻的是范德華作用、氫鍵以及脫溶能。圖6 所示是使用Autodock vina 進行對接得到的最穩定的結構，結合能為-4.16 kcal/mol，毛蕊異黃酮苷一共與EphA2形成了3條氫鍵，為毛蕊異黃酮苷的羥基氧與氨基酸的氫原子成鍵。表5 所示的是使用Autodock4對接得到的10個構象及其結合能、氫鍵數目以及氫鍵連接的殘基，涉及的主要殘基是LYS50、ARG103、MET73、GLN77、VAL72、VAL102、LYS162、SER68、ASN106 等。表6 是使用Autodock vina 進行對接得到的10 個構象及其及核能、氫鍵數目以及氫鍵連接的殘基，涉及的殘基有CYS70、ASN71、GLY51、PHE134、ASN106、TYR48、GLY49、SER169、VAL165、ARG103、PHE108 等。 綜合兩個對接方法，阿拉伯糖與EphA2 結合的殘基主要有ARG103、ASN106。

表5 使用Autodock4 對接方法所得的所有構象結合能以及結合的殘基Tab.5 All conformational binding energies and bound residues obtained using the Autodock4 docking method

表6 使用Autodock vina 對接方法所得的所有構象結合能以及結合的殘基Tab.6 All conformational binding energies and binding residues obtained using the Autodock vina docking method

圖5 EphA2與毛蕊異黃酮苷Autodock4對接結果Fig.5 Docking results of EphA2 and mullein isoflavonoe glycosides by Autodock4

圖6 EphA2與毛蕊異黃酮苷Autodock vina對接結果Fig.6 Docking results of EphA2 and mullein isoflavonoe glycosides by Autodock vina

2.2 gL與EphA2結合模式分析

據相關文獻報道，gHgL 通過結合EphA2 的不同位點而介導KSHV、EBV 等γ-皰疹病毒感染宿主細胞[5-7]。因此，本研究進一步分析gHgL 與EphA2 的結合模式。從PDB 數據庫下載gHgL 與EphA2 蛋白結合的結構，去除水分子、加氫之后載入chimera 軟件。分析顯示與EphA2 結合的主要是gL 蛋白，再使用氫鍵工具分析gL蛋白鏈與EphA2蛋白鏈之間的相互作用，得到表7，其中A 鏈為EphA2 鏈，C 鏈為gL 鏈。從表7 可以看出，gL 鏈與EphA2 蛋白鏈之間存在多個氫鍵作用，其中包括D..A 距離和D-H..A 距離兩個參數。D..A 距離表示給體原子（Donor atom）與受體原子（Acceptor atom）之間的距離，D-H..A 距離表示給體原子上的氫原子與受體原子之間的距離。在這些氫鍵中，一些給體原子是來自gL 蛋白鏈上的酰胺氮或羥基氧原子，而一些受體原子則是來自EphA2 蛋白鏈上的羰基碳或羥基氧原子。這些給體和受體原子之間形成了穩定的氫鍵作用，促進了gL 鏈與EphA2蛋白鏈之間的結合。

3 討論

黃芪總苷、黃芪總黃酮和黃芪總糖等化合物是黃芪的主要活性成分，它們通過不同的方式與目標蛋白結合，發揮著不同的生物活性。本研究主要探討黃芪主要活性成分阿拉伯糖、黃芪皂苷甲和毛蕊異黃酮苷與目標蛋白EphA2的結合模式，并通過分子對接技術對其進行預測和驗證。分子對接是一種計算機輔助的藥物設計方法，它可以預測小分子與目標蛋白之間的相互作用，并為藥物研發提供有力支持。Autodock4 是一種廣泛使用的分子對接軟件[8-9]，它基于Lamarckian 遺傳算法和能量評估函數來進行分子對接[10]。這種方法可以考慮分子的柔性和蛋白質的柔性，適用于大規模的分子對接研究。其優點在于能夠處理大量的分子和復雜的蛋白質結構，具有很高的準確性和穩定性。Autodock vina是一種基于梯度優化算法的分子對接軟件[11-12]，它可以快速地預測分子之間的相互作用。相對于Autodock4，Autodock vina 的優勢在于速度更快、準確性更高。Autodock vina可以搜索更大的對接空間，同時也考慮了分子的柔性，因此適用于大規模的分子對接研究。在進行分子對接研究時，可以根據實際需求和研究對象的不同，選擇適合的分子對接方法，以獲得更加準確可靠的結果。本研究同時使用Autodock4 和Autodock vina這兩種不同分子對接方法，能夠更全面地探索分子間的相互作用，提高分子對接的準確性和可靠性。

本研究使用Autodock4 和Autodock vina 最終得到了多個結合構象。通過分析不同構象的結合能、氫鍵數目以及氫鍵連接的殘基，對黃芪主要活性成分與EphA2 蛋白的結合進行了研究。同時我們也分析了gHgL 與EphA2 蛋白的結合模式，確定了參與結合的重要的EphA2 蛋白位點部分。通過比較黃芪主要活性成分與gHgL 蛋白對EphA2 的結合模式，發現它們都是通過形成多個氫鍵作用來實現與目標蛋白的結合。在藥物設計中，氫鍵作用通常是一種重要的相互作用方式，是一個需要深入研究的領域。通過分析氫鍵的距離和角度等參數進一步優化藥物結構，可以增強其與目標蛋白結合的親和力和特異性[13-14]。

綜上所述，本研究為黃芪主要活性成分對gHgL與EphA2 蛋白結合的干預提供了一定程度上的解釋，為靶向EphA2蛋白的干預藥物設計提供了有價值的參考，并為進一步優化黃芪主要活性成分的結構、增強其與目標蛋白結合的親和力和特異性提供了重要的參考。