劉海波 崔 巍 徐 峻 彭 勇,* 周家駒 肖培根

(1中國醫(yī)學(xué)科學(xué)院北京協(xié)和醫(yī)學(xué)院藥用植物研究所,北京 100193; 2中國科學(xué)院研究生院化學(xué)與化學(xué)工程學(xué)院,北京 100049; 3中山大學(xué)藥學(xué)院,廣州 510275;4中國科學(xué)院過程工程研究所,北京 100190)

中藥活性成分對血栓素A2受體抑制作用的分子模擬

劉海波1崔 巍2徐 峻3,*彭 勇1,*周家駒4肖培根1

(1中國醫(yī)學(xué)科學(xué)院北京協(xié)和醫(yī)學(xué)院藥用植物研究所,北京 100193;2中國科學(xué)院研究生院化學(xué)與化學(xué)工程學(xué)院,北京 100049;3中山大學(xué)藥學(xué)院,廣州 510275;4中國科學(xué)院過程工程研究所,北京 100190)

中藥中黃酮類化合物和白藜蘆醇等活性成分對血栓素A2受體具有抑制作用,但具體機(jī)理不詳.本研究通過同源模建方法,以墨魚視紫紅質(zhì)蛋白為模板,構(gòu)建血栓素A2受體的蛋白質(zhì)結(jié)構(gòu)模型.并使用分子對接方法研究中藥活性成分白藜蘆醇和芹菜苷元與血栓素A2受體模型的作用方式,據(jù)此建立藥效團(tuán)模型,篩選其他潛在的血栓素A2受體抑制劑.結(jié)果表明:白藜蘆醇等中藥活性成分能與血栓素A2受體活性口袋中的殘基發(fā)生氫鍵作用,結(jié)合方式與血栓素相似.血栓素與Ser201、Leu198、Arg295和Thr298發(fā)生氫鍵作用,白藜蘆醇等活性成分與Ser201、Leu198和Arg295發(fā)生氫鍵作用.建立的藥效團(tuán)模型由7個藥效元素以及排斥性空間元素組成,經(jīng)測試對高活性的血栓素A2受體抑制劑有比較好的選擇性.使用該藥效團(tuán)模型對中藥天然產(chǎn)物數(shù)據(jù)庫進(jìn)行篩選,命中了一批可能具有血栓素A2受體抑制作用的活性化合物.其中一些已經(jīng)報道有抑制血小板凝聚活性.本研究表明血栓素A2受體可能是活血化瘀類中藥的一個潛在的靶點.

中藥; 藥物靶標(biāo); 血栓素受體; 活血

血栓素A2(thromboxane A2,TXA2)是花生四烯酸代謝通路的產(chǎn)物之一,由前列腺素H2(prostaglandin H2,PGH2)在TXA2合成酶(thromboxane A2 synthase, TBXAS)的作用下發(fā)生異構(gòu)化生成,該反應(yīng)主要發(fā)生在血小板中.TXA2結(jié)構(gòu)不穩(wěn)定,很快水解成為血栓素B2(thromboxane B2,TXB2)(圖1),半衰期只有30 s[1].TXA2是與血小板聚集相關(guān)的重要內(nèi)源性激素,通過作用于血栓素A2受體(thromboxane A2 receptors,TP)促進(jìn)血小板凝集.PGH2也有類似活性,而水解產(chǎn)物TXB2沒有該作用[2].

自從上個世紀(jì)70年代被發(fā)現(xiàn)以來,TXA2一直是心血管藥物開發(fā)的重要研究對象[3].前期研究重點放在對TBXAS抑制劑的開發(fā)上,后來發(fā)現(xiàn)由于TXA2前體PGH2也能促進(jìn)血小板聚集,因此單純抑制TBXAS活性并不能降低血栓形成.90年代以后,研究重點轉(zhuǎn)向了TXA2/PGH2共同受體抑制劑,以及TBXAS和TP雙抑制劑的開發(fā)上[4-5].已經(jīng)證實TP抑制劑可以有效降低TXA2和PGH2的血小板聚集和血管收縮作用[6].

TP屬于G蛋白偶聯(lián)受體超家族(G-protein-coupled receptor,GPCR).GPCR在化學(xué)藥物的研發(fā)中占有重要地位,GPCR為靶標(biāo)的藥物占有很高的市場份額[7].該家族一級序列相似性不高,但具有高度的結(jié)構(gòu)相似性,每個受體都是單亞基膜蛋白,由7個α螺旋跨膜區(qū)(transmembrane,TM),3個外袢(extracellular loop,eLP),3個內(nèi)袢(cytoplasmic loop,cLP),1個細(xì)胞外N末端和1個膜內(nèi)C末端構(gòu)成[8].

TP的蛋白質(zhì)結(jié)構(gòu)研究有較長的歷史[9],但到目前為止還沒能通過實驗方法獲得蛋白質(zhì)晶體.Ruan等[10-11]使用核磁共振方法得到了eLP2和eLP3在溶液中的構(gòu)象,并據(jù)此通過分子模擬軟件,構(gòu)建了TP與小分子結(jié)合的活性口袋[12].TXA2與TP作用方式的研究一般使用與TXA2結(jié)構(gòu)相似的不易水解的化合物[9,13],通過突變、光親合標(biāo)記和特異位點抗體等方法尋找與配體發(fā)生作用的關(guān)鍵殘基[14-15].

中藥中的黃酮類化合物[16]和白藜蘆醇[17-19]被報道對TP有抑制作用,但是分子作用機(jī)制仍不清楚.由于這兩類化合物有共同結(jié)構(gòu)特征,因此很可能與TP有相似的作用方式.本研究擬通過同源建模方法構(gòu)建TP的蛋白質(zhì)結(jié)構(gòu),使用分子對接方法研究TXA2與TP的作用機(jī)制,之后進(jìn)一步研究黃酮類化合物和白藜蘆醇對TP抑制作用的機(jī)理,據(jù)此建立藥效團(tuán)模型,在中藥天然化合物中尋找其他可能的TP抑制劑.

1 數(shù)據(jù)和方法

1.1 TP結(jié)構(gòu)同源建模

1.1.1 軟件和算法

同源建模工作在靶分子信息學(xué)平臺(target informatics platform,TIP)上完成.該平臺是由Eidogen-Sertanty公司開發(fā)的結(jié)構(gòu)信息知識庫[20],包括了來自SWISS-PROT和IPI(international protein index)的人、小鼠、大鼠的序列數(shù)據(jù),以及PDB(protein data bank)數(shù)據(jù)庫中晶體結(jié)構(gòu)數(shù)據(jù).目前共有超過75000個人類蛋白結(jié)構(gòu)(包括已經(jīng)發(fā)表的晶體模型和使用STRUCTFAST建立的模型)和超過125000的活性位點信息.

圖1 血栓素A2(TXA2)代謝反應(yīng)Fig.1 Metabolic reactions of thromboxane A2(TXA2)TBXAS:thromboxane A2 synthase

STRUCTFAST(structure realization utilizing cogent tips from aligned structural templates)是TIP特有的用于蛋白質(zhì)同源建模的一種算法.該算法利用一種新的動力學(xué)程序引擎,將同族蛋白結(jié)構(gòu)特征信息結(jié)合到比對過程中,能夠篩選蛋白質(zhì)序列間的弱的相似性.有研究證明,STRUCTFAST和專家同源建模的結(jié)果相當(dāng)或者更好[21].由于TP所屬的GPCR家族具有序列一致性低,結(jié)構(gòu)相似性高的特點,因此更適合使用STRUCTFAST進(jìn)行建模.

1.1.2 操作及模板選擇

TIP平臺每次收錄新的蛋白質(zhì)序列后,會計算該序列與目前所有的蛋白質(zhì)晶體結(jié)構(gòu)的序列一致性,尋找最合適的模板進(jìn)行同源建模.登錄TIP平臺,使用TP受體序列P21731-2(SWISS-PROT編號)建立同源模型TP model,使用的模板為墨魚視紫紅質(zhì)蛋白(PDB編號:2z73)[22-23]的A鏈,序列一致性為15%,該模板也屬于GPCR家族.

1.1.3 動力學(xué)優(yōu)化

應(yīng)用AMBER03力場參數(shù)處理蛋白質(zhì)分子,并應(yīng)用gaff力場處理配體分子,配體分子電荷使用Gaussian 03[24]計算,應(yīng)用靜電勢收斂算法,在分子力學(xué)優(yōu)化的結(jié)構(gòu)上進(jìn)行單點計算以確定各原子上的電荷.分子動力學(xué)模擬計算以AMBER 10軟件[25]運行,應(yīng)用tleap程序添加缺失的氫原子,體系以削角正八面體的TIP3P水周期盒子包裹,根據(jù)體系靜電勢分布,在正電勢最集中區(qū)域加入7個Cl-離子使體系總電荷平衡.

首先對體系進(jìn)行能量優(yōu)化,能量優(yōu)化分三步進(jìn)行:第一步約束蛋白質(zhì)和配體分子中的重原子,優(yōu)化全部溶劑分子和溶質(zhì)分子中的氫原子;第二步約束蛋白質(zhì)骨架重原子,優(yōu)化全部側(cè)鏈原子和氫;第三步自由優(yōu)化全部原子.優(yōu)化結(jié)束后,開始進(jìn)行分子動力學(xué)模擬,首先在60 ps時間中分7步應(yīng)用溫度耦合算法將體系從0 K升溫至310 K,隨后在310 K進(jìn)行穩(wěn)態(tài)模擬,應(yīng)用shake算法穩(wěn)定含氫化學(xué)鍵,應(yīng)用PME算法處理長程靜電相互作用,動力學(xué)模擬步長選為2 fs,每隔0.2 ps輸出一次軌跡構(gòu)象和體系能量等信息.最終運行了6 ns模擬,從軌跡中選取勢能穩(wěn)定階段的平均構(gòu)象作為分子動力學(xué)優(yōu)化結(jié)果.

1.2 TXA2的分子對接及結(jié)合模式分析

用molecular operating environment(MOE)[26]進(jìn)行配體-受體對接研究,步驟如下:(1)將由TIP建立的模型導(dǎo)入MOE,根據(jù)TP結(jié)構(gòu)的相關(guān)實驗研究結(jié)果進(jìn)行相應(yīng)調(diào)整,使用AMBER99力場進(jìn)行能量優(yōu)化;(2)導(dǎo)入TXA2分子,使用MMFF94x力場進(jìn)行構(gòu)象能量優(yōu)化;(3)以模板蛋白中的配體分子視黃醛(retinal)為小分子模板,使用 MOE中的 flexible alignment模塊將TXA2與之進(jìn)行疊合,產(chǎn)生一個疊合構(gòu)象庫,從中根據(jù)經(jīng)驗選擇最合理的幾個構(gòu)象作為TXA2在TP中的初始構(gòu)象;(4)選取TXA2分子和受體活性口袋中0.8 nm范圍內(nèi)的氨基酸殘基,其余部分固定,使用MMFF94x力場進(jìn)行能量優(yōu)化; (5)使用potential energy模塊計算結(jié)合自由能,使用ligand interactions分析活性小分子與活性口袋中殘基間的相互作用;(6)取疊合構(gòu)象庫中TXA2的其他構(gòu)象,重復(fù)步驟4、5,將結(jié)果與實驗研究進(jìn)行比較,最終確定最合理的結(jié)合模式.

1.3 黃酮和白藜蘆醇對TP抑制的分子對接研究

基本步驟與1.2相同,黃酮類化合物選擇活性最高的芹菜苷元(apigenin)作為研究對象[16],芹菜苷元和白藜蘆醇的初始位置和構(gòu)象參考1.2節(jié)得到的模型中TXA2的結(jié)合方式.

1.4 TP抑制劑藥效團(tuán)模型的建立及數(shù)據(jù)庫篩選

1.4.1 ATCMD數(shù)據(jù)庫及3D構(gòu)象庫的建立

ATCMD(annotated traditional Chinese medicine database)是根據(jù)《中藥原植物化學(xué)成分集》[27]中的數(shù)據(jù)建立的小分子庫.該書收集整理了中外科學(xué)家發(fā)表的5507篇(部)文獻(xiàn)中的中藥成分化合物23033個,涉及藥用植物6760種,其中7819種化合物有藥理實驗數(shù)據(jù).該數(shù)據(jù)庫的化合物結(jié)構(gòu)非常精準(zhǔn)可靠,藥理數(shù)據(jù)也比較豐富,因此本研究采用該數(shù)據(jù)庫作為藥效團(tuán)篩選對象.首先將數(shù)據(jù)庫中全部化學(xué)結(jié)構(gòu)導(dǎo)出到sdf文件,之后使用德國Mol-Net公司Corina軟件[28]生成低能構(gòu)象庫,每個分子除了生成最低能量構(gòu)象外,還在最低能量構(gòu)象以上84 kJ·mol-1范圍內(nèi)保留最多20個低能構(gòu)象.原數(shù)據(jù)庫23033個分子共生成255879個構(gòu)象結(jié)構(gòu),構(gòu)象庫使用MOE database管理.

1.4.2 建立抑制劑藥效團(tuán)模型

使用MOE的Ph4 query editor模塊建立藥效團(tuán)模型.由于MOE在藥效團(tuán)建立過程中不考慮小分子的生物活性,側(cè)重于尋找小分子間三維結(jié)構(gòu)的相似性,所建立的藥效團(tuán)模型比較類似于目前流行的藥效團(tuán)軟件CATALYST[29]建立的common features hypotheses.這種算法比較適合于在配體小分子數(shù)量比較少或者缺少充足的活性數(shù)據(jù)的情況下建立藥效團(tuán)模型.MOE在建立藥效團(tuán)的過程中能加入受體口袋的空間因素,因此所建立的藥效團(tuán)模型可以快速篩選到體積適合,并且具有相似三維構(gòu)象的分子,但不能對生物活性進(jìn)行預(yù)測.

具體步驟如下:(1)從1.3節(jié)建立的模型中分離出芹菜苷元和白藜蘆醇分子,保持構(gòu)象疊合在一起; (2)在 Ph4 query editor中設(shè)定 PCHD(polaritycharge-hydrophobicity-direction)藥效元素組合,使用Consen-sus tool枚舉出所有可能的藥效團(tuán)元素;(3)結(jié)合1.3節(jié)中小分子與TP相互作用模式及相關(guān)實驗研究結(jié)果,選擇藥效團(tuán)元素組合方式,并對每個元素的體積等條件進(jìn)行調(diào)整;(4)選擇TP活性口袋中與配體小分子距離0.6 nm范圍內(nèi)的氨基酸殘基,據(jù)此建立排除體積(excluded volume)模型.

為了考查建立藥效團(tuán)模型的選擇效率,建立了一個測試數(shù)據(jù)集.根據(jù)文獻(xiàn)報道[16],收集了20種黃酮類化合物,對TP的抑制率在0%-64%之間.使用與1.4.1節(jié)中相同的方法建立測試分子3D結(jié)構(gòu)數(shù)據(jù)庫,使用建立的模型對數(shù)據(jù)庫進(jìn)行藥效團(tuán)篩選,根據(jù)命中分子的情況調(diào)整藥效元素直到確定最好的藥效團(tuán)模型.

1.4.3 3D分子庫的藥效團(tuán)篩選

在MOE database的Ph4 database search模塊中使用1.4.2節(jié)中生成的藥效團(tuán)模型對ATCMD構(gòu)象庫進(jìn)行篩選,對命中的天然產(chǎn)物的化學(xué)結(jié)構(gòu)類型和藥理活性進(jìn)行分析整理.

2 結(jié)果與討論

2.1 TP同源建模結(jié)果

Ramachandran圖是蛋白質(zhì)主鏈的ψ-φ角度分布圖,是蛋白質(zhì)構(gòu)象合理性的一種直觀的表述.如果高能角出現(xiàn)頻率過高,則說明該構(gòu)象不合理或者與實際構(gòu)象差別比較大.MOE的protein geometry模塊能對蛋白質(zhì)構(gòu)象生成四種Ramachandran圖[30],分別對應(yīng)于一般氨基酸,甘氨酸,脯氨酸和前脯氨酸(pre-proline),多數(shù)高能角出現(xiàn)在一般氨基酸圖中.對經(jīng)過能量優(yōu)化的同源建模模型TP model及其模板2z73A作Ramachandran圖(圖2).可以看到,模板蛋白質(zhì)只有一個高能角,TP model在一般氨基酸圖中有9個高能角,只占總數(shù)的2.7%,主要分布在內(nèi)袢和外袢區(qū),7個過膜區(qū)沒有高能角,可以認(rèn)為由TIP建立的模型基本上是合理的.

由TIP模型為基礎(chǔ)進(jìn)行動力學(xué)模擬,結(jié)果顯示體系的勢能在500 ps內(nèi)達(dá)到穩(wěn)態(tài)狀態(tài).均方根偏差(RMSD)1000 ps后也基本達(dá)到平衡狀態(tài)(圖3).對最終得到的平均能量構(gòu)象進(jìn)行穩(wěn)定性分析,發(fā)現(xiàn)高能角的數(shù)量由9個降低到3個(圖2),說明模型已經(jīng)得到了優(yōu)化,非常穩(wěn)定,可以進(jìn)行下面的分子對接研究.

2.2 TXA2對接結(jié)果

TXA2與TP的對接結(jié)果顯示(圖4),TXA2與TP的結(jié)合部位位于第三外袢以下,7個過膜螺旋之間的部位,分子的羧基端與Arg295形成兩個氫鍵,與Thr298形成一個氫鍵,羥基端與Ser201和Leu198形成一個氫鍵.結(jié)合自由能也比較大,達(dá)到-633.14kJ·mol-1,TXA2與TP間的靜電作用為主要作用力(表1).

圖2 同源建模模型和模板的蛋白質(zhì)構(gòu)象以及Ramachandran圖Fig.2 Conformation and Ramachandran plots of homology model and templateTIP:target information platform,MD:molecular dynamics

文獻(xiàn)中TXA2與TP結(jié)合模式的實驗研究結(jié)果主要有:(1)TXA2類似物和受體特異性拮抗作用發(fā)現(xiàn),TP的配體結(jié)合部位為TM5的Ser201與TM7的Arg295之間的一個疏水性口袋[31],從圖4a可以看到對接模型與該實驗結(jié)果吻合;(2)TXA2及其類似物的羧基與Arg295,羥基與Ser201發(fā)生氫鍵作用,而分子中央的疏水性部分位于受體的疏水性口袋中[31],從圖4b看,模型與實驗吻合;(3)Turek等[14]用光親合標(biāo)記法和特異位點抗體法研究發(fā)現(xiàn),第三外袢Cys183到Asp193是直接參與TP和配體結(jié)合的氨基酸序列,從圖4c看,模型中TXA2正位于eLP3的Cys183到Asp193的部分之下,與該段肽鏈有緊密的作用,與實驗吻合;(4)Ruan等[13]使用PGH2相似分子進(jìn)行對接研究發(fā)現(xiàn),TP和配體結(jié)合時,配體的兩條臂充分打開,羧基一端盡量伸展,羥基一端在羥基以外的碳鏈向回折,形成一種長方形的形式與TP結(jié)合,從圖4d看,模型中的TXA2分子的構(gòu)象與實驗基本吻合.

表1 分子對接的結(jié)合自由能Table 1 Binding free energy obtained from docking results

圖5 芹菜苷元和白藜蘆醇與TP受體間的對接結(jié)果Fig.5 Docking results for the apigenin,resveratrol and TP(a)hydrogen bond interaction;(b)docking result

綜上所述,本研究得到的對接模型與上面四個實驗結(jié)果都保持一致,因此認(rèn)為這種結(jié)合方式比較符合實際情況.

2.3 芹菜苷元和白藜蘆醇的分子對接結(jié)果

結(jié)果顯示,芹菜苷元和白藜蘆醇以相同的方式結(jié)合到 TP活性口袋里,與 Ser201、Leu198和Arg295發(fā)生氫鍵作用(圖5).這種作用方式與TXA2相似,不同之處在于沒有與Thr298產(chǎn)生氫鍵的作用.從結(jié)合自由能看,白藜蘆醇的結(jié)合強(qiáng)度比較大,芹菜苷元較弱(表1).

圖6 芹菜苷元和白藜蘆醇的藥效團(tuán)模型Fig.6 Pharmacophore model for the apigenin and resveratrolF1,F2:hydrogen bond acceptor;F3:hydrogen bond donor&acceptor; F4:hydrogen bond acceptor 2;F5:hydrogen bond donor 2; F6:hydrogen bond donor 2/acceptor 2;F7:hydrophobes/aromatic

圖7 藥效團(tuán)模型性能測試結(jié)果Fig.7 Test results for the pharmacophore modelNHDC:neohesperidin dihydrochalcone;the hit compounds have been marked with asterisk.

GPCR激動劑和抑制劑有不同的結(jié)合方式,前者能夠引起跨膜螺旋區(qū)發(fā)生重排,從而改變細(xì)胞內(nèi)C端的構(gòu)象,暴露出偶聯(lián)G蛋白識別位點,最終與G蛋白形成復(fù)合物.而抑制劑不能引起相應(yīng)的構(gòu)象變化,因此使GPCR的活性受到抑制[8].分子對接研究的結(jié)果顯示,TXA2與TP作用于Ser201、Leu98、 Arg295和Thr298,而芹菜苷元和白藜蘆醇等抑制劑與Ser201、Thr298、Arg295均有氫鍵作用,但沒有與Thr298形成氫鍵.這種結(jié)合方式的差異可能導(dǎo)致了TP構(gòu)象的不同變化.因此,TXA2對TP有激活作用,芹菜苷元和白藜蘆醇等抑制劑不能激活TP,而且對TXA2產(chǎn)生了競爭性抑制作用.這可能是抑制作用的分子機(jī)理.

2.4 藥效團(tuán)模型及篩選結(jié)果

2.4.1 芹菜苷元和白藜蘆醇的藥效團(tuán)模型

如圖6所示,建立的藥效團(tuán)模型由7個藥效元素構(gòu)成,F1、F2、F3為氫鍵受體或者給體,F4、F5、F6為外部氫鍵受體或者給體,分別對應(yīng)于芹菜苷元和白藜蘆醇兩端與受體發(fā)生氫鍵結(jié)合的位置.其中, F2-F5限制比較嚴(yán)格,半徑在0.03-0.05 nm之間. F1和F6比較寬松,半徑為0.15 nm.中間疏水部分由一個疏水/芳香元素F7構(gòu)成,半徑均為0.6 nm.此外還根據(jù)受體活性口袋建立了排斥性空間元素,每個空間元素的半徑為0.13 nm.

2.4.2 藥效團(tuán)篩選及結(jié)果分析

使用藥效團(tuán)模型對測試數(shù)據(jù)庫中20個化合物進(jìn)行篩選,共命中10個化合物(圖7).檢中的化合物多數(shù)集中在抑制率20%以上的活性區(qū),抑制率小于5%低活性化合物沒有發(fā)生誤檢的情況,說明該藥效團(tuán)模型對天然產(chǎn)物中的TP抑制劑有一定的篩選能力,可以滿足初步篩選的需要.但是仍有一些問題存在,例如對黃酮骨架的依賴性比較強(qiáng),因此對于genistein等異黃酮結(jié)構(gòu)的化合物預(yù)測能力較差.使用藥效團(tuán)模型對ATCMD-3D構(gòu)象庫進(jìn)行篩選,共命中261種化合物,其中8種有報道具有抑制血小板聚集活性.其他命中的分子是否有TP抑制作用,還有待進(jìn)一步研究.

3 結(jié) 論

中藥藥效的物質(zhì)基礎(chǔ)研究是中藥現(xiàn)代化的重要組成部分,中藥作用的靶點分析又是中藥物質(zhì)基礎(chǔ)的重點研究方向.本文通過同源建模方法,建立了血栓素A2受體蛋白的結(jié)構(gòu)模型,模擬了血栓素與受體的作用方式,并與實驗結(jié)果比對驗證了該模型的可靠性.據(jù)此,進(jìn)一步分析芹菜苷元和白藜蘆醇對TP的抑制機(jī)理,并構(gòu)建了抑制劑的藥效團(tuán)模型,在ATCMD構(gòu)象數(shù)據(jù)庫中進(jìn)行篩選,命中了一批可能具有TP抑制活性的天然產(chǎn)物.研究結(jié)果顯示,TP很可能是中藥活血作用的重要靶標(biāo).這個結(jié)果對基于中藥的先導(dǎo)化合物發(fā)現(xiàn)和新藥開發(fā)有一定的啟發(fā)作用,更深入的研究將在今后逐步展開.

1 http://www.genome.jp/kegg/

2 Chen,X.Cardiovascular pharmacology.3rd ed.Beijing:People′s Medical Publishing House,2002:129 [陳 修.心血管藥理學(xué).第三版.北京:人民衛(wèi)生出版社,2002:129]

3 Hamberg,M.;Svensson,J.;Samuelsson,B.Proc.Natl.Acad.Sci. U.S.A.,1975,72:2994

4 The RAPT Investigators.Circulation,1994,89:588

5 Squadrito,F.;Ioculano,M.;Altavilla,D.;Zingarelli,B.;Canale,P.; Campo,G.M.;Saitta,A.;Oriti,S.;Spignoli,G.;Caputi,A.P. Pharmacology,1993,47:167

6 Albuquerque,M.G.;Hopfinger,A.J.;Barreiro,E.J.;de Alencastro,R.B.J.Chem.Inf.Comput.Sci.,1998,38:925

7 Saunders,J.Bioorganic&Medicinal Chemistry Letters,2005,15:3653

8 Devlin,T.M.Textbook of biochemistry with clinical correlations. 6th ed.Trans.Wang,H.Y.Beijing:Science Press,2008:601 [Devlin,T.M.生物化學(xué)——基礎(chǔ)理論與臨床.王紅陽,譯.北京:科學(xué)出版社,2008:601]

9 Hirata,M.;Hayashi,Y.;Ushikubi,F.;Yokota,Y.;Kageyama,R.; Nakanishi,S.;Narumiya,S.Nature,1991,349:617

10 Ruan,K.H.;So,S.P.;Wu,J.;Li,D.;Huang,A.;Kung,J. Biochemistry,2001,40:275

11 Wu,J.;So,S.P.;Ruan,K.H.Arch.Biochem.Biophys.,2003,414: 287

12 Ruan,K.H.;Wu,J.;So,S.P.;Jenkins,L.A.;Ruan,C.H.Eur.J. Biochem.,2004,271:3006

13 Ruan,K.H.;Wijaya,C.;Cervantes,V.;Wu,J.Arch.Biochem. Biophys.,2008,477:396

14 Turek,J.W.;Halmos,T.;Sullivan,N.L.;Antonakis,K.;Le Breton,G.C.J.Biol.Chem.,2002,277:16791

15 Funk,C.;Furci,L.;Moran,N.;Fitzgerald,G.Mol.Pharmacol., 1993,44:934

16 Navarro-Nunez,L.;Castillo,J.;Lozano,M.L.;Martinez,C.; Benavente-Garcia,O.;Vicente,V.;Rivera,J.J.Agric.Food Chem., 2009,57:1589

17 Gehm,B.D.;McAndrews,J.M.;Chien,P.Y.;Jameson,J.L.Proc. Natl.Acad.Sci.U.S.A.,1997,94:14138

18 Olas,B.;Wachowicz,B.Platelets,2005,16:251

19 Yang,Y.M.;Chen,J.Z.;Wang,X.X.;Wang,S.J.;Hu,H.;Wang, H.Q.Eur.J.Pharmacol.,2008,583:148

20 Target informatics platform.San Diego:Eidogen-Sertanty,2009

21 Debe,D.A.;Danzer,J.F.;Goddard,W.A.;Poleksic,A.Proteins, 2006,64:960

22 http://www.rcsb.org/pdb/home/home.do

23 Murakami,M.;Kouyama,T.Nature,2008,453:363

24 Frisch,M.J.;Trucks,G.W.;Schlegel,H.B.;et al.Gaussian 03. Revision C.02.Wallingford,CT:Gaussian Inc.,2004

25 Case,D.A.;Cheatham,T.E.;Darden,T.;Gohlke,H.;Luo,R.; Merz,K.M.;Onufriev,A.;Simmerling,C.;Wang,B.;Woods,R.J. J.Comput.Chem.,2005,26:1668

26 Molecular operating environment.Version 2007.09.Sherbrooke Street West,Suite 910,Montreal,H3A 2R7,Canada:Chemical Computing Group Inc.,2007

27 Zhou,J.J.;Xie,G.R.;Yan,X.J.Chemical components of source plants in traditional Chinese medicine.Beijing:Science Press,2009 [周家駒,謝桂榮,嚴(yán)新建.中藥原植物化學(xué)成分集.北京:科學(xué)出版社,2009]

28 Corina 3.4.Erlangen,Germany:Molecular Networks GmbH,2008

29 Catalyst software package.Release 4.7.2003

30 Ho,B.K.;Brasseur,R.BMC Struct.Biol.,2005,5:14

31 Yamamoto,Y.;Kamiya,K.;Terao,S.J.Med.Chem.,1993,36: 820

Molecular Simulations of the Inhibition of Active Components in Traditional Chinese Medicine on the Thromboxane A2 Receptor

LIU Hai-Bo1CUI Wei2XU Jun3,*PENG Yong1,*ZHOU Jia-Ju4XIAO Pei-Gen1
(1Institute of Medicinal Plant Development,Peking Union Medical College,Chinese Academy of Medical Sciences, Beijing 100193,P.R.China;2College of Chemistry and Chemical Engineering,Graduate University of Chinese Academy of Sciences,Beijing 100049,P.R.China;3School of Pharmaceutical Sciences,Sun Yat-sen University,Guangzhou 510275,P.R.China;4Institute of Process Engineering,Chinese Academy of Sciences,Beijing 100190,P.R.China)

Some natural products from traditional Chinese medicine(TCM)such as flavones and resveratrol have been reported to have thromboxane A2 receptor(TP)inhibiting activity.We investigated a possible inhibition mechanism by a homology model of TP,which was built based on the crystal structure of squid rhodopsin.After that,docking methods were used to investigate the binding modes of resveratrol and apigenin in the active pocket of TP.Furthermore,a threedimensional pharmacophore model was generated for screening other potential natural TP inhibitors.The results indicate that resveratrol and apigenin bind to the active site of TP similar to the way that thromboxane A2 binds to Ser201,Leu198,Arg295,and Thr298.The former three key residues can form hydrogen bonds with the inhibitors.The pharmacophore model consisted of seven features and a set of volume spheres,which has been proven to be efficient in identifying compounds with high TP inhibition activity.In this way,a set of potential TP inhibitors were screened from a natural product database.Some of them were reported to have platelet aggregation inhibiting activities.This research indicates that TP could be an important target of TCM drugs with blood circulation activation effects.

Traditional Chinese medicine;Drug target;Thromboxane A2 receptor; Blood circulation activating

O641

Received:February 10,2010;Revised:May 5,2010;Published on Web:July 7,2010.

*Corresponding authors.XU Jun,Email:junxu@biochemomes.com.PENG Yong,Email:ypeng@implad.ac.cn;Tel:+86-10-51296683,+86-10-62894462.

The project was supported by the National Natural Science Foundation of China(30530860).

國家自然科學(xué)基金重點項目(30530860)資助

?Editorial office of Acta Physico-Chimica Sinica