馮 惠，秦正龍，李 靖，馮長君

(1. 徐州工程學院材料與化學工程學院，中國 徐州 221018；2.江蘇師范大學化學與材料科學學院，中國 徐州 221116)

因拓撲異構酶既是抗菌作用靶標，也是抗腫瘤藥物的重要靶點，由此可通過結構修飾將以喹啉-4-酮-3-羧酸為優(yōu)勢藥效團的氟喹諾酮藥物的抗菌活性轉化為抗腫瘤活性[1,2]。值得關注的是DNA拓撲異構酶Ⅱ(Topo II)可能是一些雜環(huán)縮氨基硫脲類化合物的抗腫瘤作用靶點，而該靶點也是抗腫瘤藥物研究的重要靶點[3]。胡國強等[4,5]發(fā)現(xiàn)氟喹諾酮中C-3羧基并非是抗腫瘤活性所必要的藥效團，被其等排體如酰腙、唑雜環(huán)等替代可提高其抗腫瘤活性。因此，謝玉鎖等[6]利用藥效團拼合原理，用均三唑雜環(huán)作為抗菌培氟沙星(Pefloxacin)C-3羧基的等排體，并引入縮氨基硫脲，成功合成24個培氟沙星均三唑硫醚衍生物。他們采用MTT法測定上述化合物對小鼠白血病細胞(L1210)的體外抗增殖活性(IC50)，其半數(shù)抑制濃度明顯低于前體藥培氟沙星。馮惠等[7]采用半經(jīng)驗計算法 MOPAC-AM1研究15種氟喹諾酮均三唑稠雜環(huán)類衍生物( FQTHs)體外抗腫瘤活性，獲得良好的擬合性。趙曉輝等[8]釆用 Sybyl-X 2.0軟件中的比較分子力場分析(CoMFA)、比較分子相似性指數(shù)分析(CoMSIA)方法研究12種氟喹諾酮(FQs)同系物光降解半衰期(t1/2)，預測能力令人滿意。目前尚未見到有關培氟沙星(Pefloxacin)同系物的定量結構-活性相關(Quantitative Structure-Activity Relationship, QSAR)研究[9-12]的報道。

徐光憲曾稱定量結構-活性相關為21世紀化學的4個“世紀難題”之一，因它涉及化學的根本問題。QSAR是基于數(shù)理及統(tǒng)計方法研究分子結構與性能之間定量關系的學科，已在預測化合物理化性質與生物活性方面成為化學、醫(yī)藥、環(huán)境等學科的熱門領域。目前，QSAR主要分為二維(2D)和三維(3D)兩類[13]。后者由于引入分子結構的三維構象信息，通常能夠獲得質量更佳的QSAR模型。在3D-QSAR方法中最為常用的是1988年Crameretal提出的比較分子力場分析(comparative molecular field analysis，CoMFA)[14-21]。本文基于CoMFA方法建立培氟沙星均三唑硫醚衍生物對小鼠白血病細胞(L1210)體外抗增殖活性(pIC)[6]的3D-QSAR模型，以揭示影響pIC的分子結構因素，探討抑制小鼠白血病細胞(L1210)作用的微觀機理。

1 數(shù)據(jù)來源

24種培氟沙星均三唑硫醚衍生物的分子基本結構[6]見圖1，其中取代基R和X見表1。在表1中的“IC50”是上述化合物對小鼠白血病細胞(L1210)的半數(shù)抑制濃度，其單位為μmol·L-1。

根據(jù)化學平衡理論，化學反應自由能變(ΔG)與體系平衡濃度為對數(shù)關系。故令：

pIC= -lg[IC50/(mol·L-1)]。

(1)

它們的pIC數(shù)據(jù)見表1。

2 建模方法

培氟沙星均三唑硫醚衍生物對小鼠白血病細胞體外抗增殖活性[6]的3D-QSAR模型，采用Tripos公司Sybyl-x2.1.1藥物分子模擬軟件予以建立。

2.1 化合物分子藥效構象的確定及分子疊合

使用Sketch molecule模塊構建24個培氟沙星均三唑硫醚衍生物分子[6]和新設計3個分子的初始三維結構；對上述分子進行分子力學能量優(yōu)化的程序為Minimize模塊，以獲得它們能量最低構象，即藥效構象(亦稱生物活性構象)，用于分子疊合。

圖2 訓練集(a)、測試集(b)的疊合圖 Fig. 2 The plots of all the aligned molecules in training set(a) and test set(b)

隨機抽取第3，9，16，22及23號分子以及新設計的3個分子(25～27號分子)為測試集(Test set)，余下20個分子作為訓練集(Training set)。在Align database中采用原子契合的公共骨架疊合方式，其公共骨架為樣本中生物活性最強的分子(即模板分子)所提供。該疊合方式能夠盡量使分子間相互重疊時取向相同，疊合誤差最小。Align database模塊以對L1210細胞體外抗增殖活性最強的22號分子提供疊合的公共骨架，Align database模塊給出的訓練集、測試集疊合圖見圖2，圖中顯示各個分子間密切疊合。

表1 培氟沙星均三唑硫醚衍生物分子結構與抗增殖活性(pIC)

2.2 CoMFA模型的建立

采用CoMFA計算訓練集中每個疊合分子周圍所有網(wǎng)格點上的立體場(Steric，St)及靜電場(Electrostatic，El)的場能值。以對小鼠白血病細胞(L1210)體外抗增殖活性pIC為因變量，場能值為自變量，用偏最小二乘法(Partial least squares，PLS)進行建模。一是以交叉驗證系數(shù)Rcv2衡量模型的預測能力，要求模型Rcv2>0.5，其可信度優(yōu)于95%，統(tǒng)計意義十分顯著[22]；二是進行回歸分析，即以最佳主成分數(shù)N建立非交叉QSAR模型，獲得相關系數(shù)R2、統(tǒng)計方差比F、估計標準誤差SD；最后通過View CoMFA模塊呈現(xiàn)三維等勢圖，以直觀反映立體場和靜電場對培氟沙星均三唑硫醚衍生物pIC的貢獻。

3 結果與討論

3.1 pIC的CoMFA模型

培氟沙星均三唑硫醚衍生物對小鼠白血病細胞(L1210)體外抗增殖活性pIC的3D-QSAR模型見表2。

表2 訓練集pIC的CoMFA模型

模型質量檢驗如下：

(1)交叉驗證部分：Rcv2>0.5，表明模型的或然性很低，具有良好的預測能力與穩(wěn)健性。

(2)非交叉驗證部分：①判定系數(shù)檢驗：具有良好擬合性的QSAR模型，其R2≥0.80[23]。R2=0.926>0.8，顯示良好的相關性。由于判定系數(shù)R2又稱為削減誤差比例，說明該模型包含了影響pIC的92.6%以上因素，未被揭示的因素所占份額不足7.4%，因此可用于pIC的估算與預測。②顯著性檢驗：在95%顯著水平下，訓練集CoMFA模型的F臨界值為：F0.05(4,15) =3.06。所建模型的F是其15倍之多，表明模型規(guī)定的因果關系顯著，有95%以上的可信度。

(3)外部預測能力：利用訓練集的3D-QSAR模型對測試集中分子的pIC進行預測，以檢驗其外部預測能力。由表1可見，對pIC預測的平均絕對預測誤差僅為0.139，與相應實驗值非常吻合，顯示很好的外部預測能力。

3.2 CoMFA等勢圖

CoMFA模型的三維等勢圖分為立體作用分布圖和靜電作用分布圖。

圖3 CoMFA模型的立體場(a)與靜電場(b)等勢圖Fig. 3 CoMFA contour maps (a) steric field; (b) electrostatic field

以化合物22為模板分子的pIC的三維立體場等勢圖見圖3(a)(彩圖見封三)。在縮氨基硫脲部位有2個綠色區(qū)域，該區(qū)域引入體積較大的取代基有助于提高化合物的抗增殖活性。因此，1～12號化合物沒有縮氨基硫脲部分，其抗增殖活性都弱于相應13～24化合物。另外，在苯環(huán)的3,4-位之間(偏4-位)有個黃色區(qū)域，該區(qū)域引入體積小的基團會使化合物的生物活性增強。例如6和18號化合物的3,4-位基團為—OCH2—，7和19號化合物的3,4-位為MeO—，它們的體積均比2～5及14～17號化合物的取代基大，因此，化合物6，7，18及19的抗增殖活性較弱。

以化合物22為模板分子的pIC的三維靜電場等勢圖見圖3(b)。在縮氨基硫脲部位有個靠近氨基中氫原子較大的藍色區(qū)域，在此區(qū)域引入正電性基團利于生物活性增強。后12個化合物具有縮氨基硫脲，因此，它們的抗增殖活性都強于(相應)前12個化合物。在苯環(huán)的4-位有個紅色區(qū)域(忽略較遠處的一個小紅色區(qū)域)，此區(qū)域引入負電性取代基有助于增強化合物的抗增殖活性。例如4-位上氟原子是最強的負電性基團，故10號在前12個化合物中、22號在后12個化合物均有最強的抗增殖活性。尚需指出，12號和24號化合物中苯環(huán)的4-位為硝基(—NO2)，雖然—NO2中存在2個電負性僅次于氟原子的氧原子，但其抗增殖活性均為前12和后12中最弱(4.246與4.629)。因這2個負電性氧原子均部分處于紅色區(qū)域，而正電性的氮原子(其電負性低于氧原子)卻處于紅色區(qū)域，導致抗增殖活性急劇下降。

表2給出立體場和靜電場對pIC的貢獻分別為72.8%和27.2%，顯示取代基的立體作用2.6倍于靜電作用。立體場通常含有疏水作用和空間位阻，庫侖力、氫鍵及配位歸于靜電場。因此，影響pIC主要取決于取代基的疏水作用和空間位阻，其次是與靶標形成庫侖力、氫鍵及配位作用。

3.3 分子設計

由3.2分析可知，在苯環(huán)的4-位上引入小體積的負電性基團，能夠增加培氟沙星均三唑硫醚衍生物對小鼠白血病細胞(L1210)體外抗增殖活性pIC。據(jù)此設計3個化合物(見表1中第25～27號化合物)，CoMFA模型給出的pIC預測值為5.224～5.332，均優(yōu)于22號分子預測值(5.130)。

4 結論

(1)基于CoMFA建立培氟沙星均三唑硫醚衍生物對小鼠白血病細胞(L1210)體外抗增殖活性的3D-QSAR模型：Rcv2=0.659，R2=0.926，F(xiàn)=47.21，顯示優(yōu)良的魯棒性、擬合性和預測能力。

(2)CoMFA模型的立體場和靜電場對pIC的貢獻分別為72.8%和27.2%，顯示立體作用近3倍于靜電作用。

(3) CoMFA模型三維等勢圖顯示在苯環(huán)的4-位上引入小體積的負電性基團，使培氟沙星均三唑硫醚衍生物對小鼠白血病細胞(L1210)的體外抗增殖活性增強。據(jù)此設計的3個化合物，顯示較大的pIC預測值，有待生物醫(yī)學實驗確認。