婁永峰向明禮楊明理

（1．四川大學華西醫(yī)院生物治療國家重點實驗室，成都，610041；2．四川大學化工學院，成都，610065）

科研與開發(fā)

手性咪唑環(huán)芳識別氨基酸甲酯的作用模式研究

婁永峰1，2向明禮1楊明理1

（1．四川大學華西醫(yī)院生物治療國家重點實驗室，成都，610041；2．四川大學化工學院，成都，610065）

采用Gaussian 03在B3LYP／6－31G（d）水平上分別優(yōu)化手性咪唑環(huán)芳受體和氨基酸甲酯配體，并在同樣條件下計算NBO電荷。用AutoDock Vina將優(yōu)化后的受體、配體進行對接，獲得相應復合物的最穩(wěn)定構象。以復合物親和能和配體側鏈電子電荷指數，建立了表達復合物穩(wěn)定常數的定量關系模型。從單體的NBO電荷分布、復合物最穩(wěn)定構象等方面分析探討，發(fā)現(xiàn)靜電吸引、受體配體間的空間匹配和氫鍵作用等，是影響手性咪唑環(huán)芳選擇性識別氨基酸甲酯的主要影響因素。

咪唑環(huán)芳 分子力學 密度泛函理論 分子識別 氨基酸甲酯

1 前言

環(huán)芳是一類具有特殊結構的有機化合物。1949年，C．J．Brown和A．C．Farthing發(fā)現(xiàn)了一種很獨特的化合物，其結構為由兩個次乙基橋從對位連接起來的面對面的苯環(huán)。該化合物被直觀地命名為“二對二甲苯”［1］。1951年，Donald J．Cram和H．Steinberg為研究剛性分子中不飽和中心間的電子相互作用，合成得到了一系列包括C．J．Brown和A．C．Farthing之前發(fā)現(xiàn)的化合物。他們將這些“面對面”的苯環(huán)環(huán)狀大分子取名為“對位環(huán)芳（paracyclophane）”［2］。隨著研究的深入，“環(huán)芳”的內涵和外延不斷地得到一步步拓展［3］，并逐漸形成為一門倍受關注的環(huán)芳化學［4］，Donald J．Cram也因為對環(huán)芳研究的獨特貢獻而獲得1987年的諾貝爾化學獎［5］。

幾十年來，有關環(huán)芳的研究一直受到持續(xù)關注。環(huán)芳已廣泛應用于人工酶模擬［5－6］、手性識別［7－9］、超分子［10］、分子催化［5，11］、抗菌［12］等領域。咪唑環(huán)芳更為特殊。咪唑是生物體系的重要基團［8］。含咪唑環(huán)的一些化合物具有抗炎活性［13］，并對細胞色素P450酶有抑制活性［14］。再加上咪唑衍生物難于合成等諸多因素［15－17］，使得咪唑環(huán)芳的研究，更引人矚目。

但是迄今，環(huán)芳研究工作的注意力幾乎都集中在合成品種不同功能各異的環(huán)芳以及發(fā)展新的合成方法［4，6，18－19］方面，人們對于環(huán)芳識別氨基酸甲酯的分子作用機制的認識［9，20］，未見深入研究的相關文獻。而手性咪唑環(huán)芳受體對D構型氨基酸甲酯配體，表現(xiàn)出不同尋常的識別能力［17］。其作用機制值得深入探討。

2 研究方法

2．1 單體結構和復合物構象的獲取

目前還沒有該手性咪唑環(huán)芳受體［17］的晶體結構，也沒有該受體識別丙氨酸甲酯（LA）、苯丙氨酸甲酯（LF）、亮氨酸甲酯（LL）色氨酸甲酯（LT）、和纈氨酸甲酯（LV）后形成復合物的晶體結構數據。而沒有環(huán)芳受體與配體間形成的復合物的可靠三維結構數據，要深入探究其間的分子識別機制，幾乎無從下手。

對于受體和配體的單體結構，我們采用Chem－Bio3D構建，用Gaussian 03［21］在B3LYP／6－31G（d）理論水平上進行幾何構型優(yōu)化得到。優(yōu)化后的手性咪唑環(huán)芳受體的三維結構圖，如圖1（a）所示。

用AutoDock Vina［22］將此前優(yōu)化過的單體進行對接。為了充分考慮分子識別時的空間互補效應，將配體作全柔性處理。AutoDock在評價分子間的相互作用時，特別考慮了氫鍵和疏水相互作用。經“迭代局域搜索全局優(yōu)化”［22］算法計算，得到最穩(wěn)定的復合物的構象。

2．2 識別模式和復合物穩(wěn)定常數的定量表達

在選擇性的分子絡合中，受體、配體的結合位點間的互補性和復合物中配體的空間位阻是關鍵因素，而靜電作用、范德法力、氫鍵等也是形成分子復合物的推動力［5－6，23－25］。有鑒于此，對手性咪唑環(huán)芳受體和氨基酸甲酯配體，在構型優(yōu)化的基礎上，在同樣的理論水平上進行了NBO電荷計算。

對Vina對接時所計算的親和能的構成［22］進行分析可知，其親和能的計算考慮了空間效應、排斥效應、氫鍵、疏水作用，以及分子的柔性。為此，以Vina所計算的親和能為核心參數，建立起表達復合物穩(wěn)定常數的定量關系模型，以便進一步深入探討該手性咪唑環(huán)芳與配體分子間的識別機制。

式（1）中，K為復合物穩(wěn)定常數計算值；ΔE為Vina計算的分子間相互作用的親和能；ECI為氨基酸側鏈的電子電荷指數［26］，α、β為待定參數，γ為截距。

3 結果與討論

3．1 手性咪唑環(huán)芳識別氨基酸甲酯的作用模式

NBO電荷計算表明，手性咪唑環(huán)芳的負電荷，大都集中在環(huán)芳的幾個N、O原子上。在咪唑環(huán)兩端起橋接作用的4個C原子，也帶有相當數量的負電荷。而在氨基酸甲酯上，氨基N上的負電荷最多，達－0．70左右，其次是羰基O和酯O，所帶負電荷約為－0．50；而氨基H最具極性，所帶正電荷在0．30左右。手性咪唑環(huán)芳受體對氨基酸甲酯的識別，靜電相互作用不可忽視。

從實驗數據（見表1）可知，手性咪唑環(huán)芳受體A與五個配體相互作用形成的復合物，其穩(wěn)定程度可分為三檔；復合物A－LT最穩(wěn)定；其次是A－LF。余下的三種穩(wěn)定性最差。

圖1 經優(yōu)化后的手性咪唑環(huán)芳（A）三維棒狀結構圖（a），受體與丙氨酸甲酯（b）受體與色氨酸甲酯（c）相互作用的最低能量構象圖（為便于觀察，隱去了配體中的非極性氫原子；紅色:氧原子；藍色:氮原子；灰色:碳原子；白色:氫原子）

將最穩(wěn)定的復合物A－LT（圖1c）和穩(wěn)定性最差一檔的復合物A－LA（圖1b）的構象進行對比分析，可以發(fā)現(xiàn)該手性咪唑環(huán)芳識別氨基酸甲酯的作用模式特征。首先是靜電作用。色氨酸甲酯帶負電的羰基O對著環(huán)芳上一個咪唑環(huán)的H；而其帶負電的酯甲基C，指向環(huán)芳上一個吡啶的H。同時，兩個極性的氨基H分別指向環(huán)芳上的兩個羰基O。在復合物A－LA中，也有類似的電荷匹配現(xiàn)象。其次是空間匹配。在最穩(wěn)定的A－LT復合物中，色氨酸甲酯的C端嵌在由吡啶環(huán)為墻、咪唑環(huán)為底的凹槽中，而其苯并吡咯側鏈則“枕”在鏈接另一個吡啶環(huán)和咪唑環(huán)的僑鏈上。猶如人坐在沙發(fā)。而在A－LA復合物中，由于受體的空間大，配體的體積小，丙氨酸甲酯在其中猶如小腳穿大鞋般，匹配性差。另外還有氫鍵作用。在A－LT中存在氫鍵作用。色氨酸甲酯的一個氨基H與受體上的一個羰基O僅距2．08?，且O－H－N間呈現(xiàn)接近180°的角度。在A－LA中，丙氨酸甲酯的氨基 H 與受體上的羰基 O相距達2．40?，且O－H－N夾角接近150°。因此可以認為，該手性咪唑環(huán)芳對氨基酸甲酯的選擇性識別，是靜電吸引、空間匹配和氫鍵等因素在起重要作用。

3．2 復合物穩(wěn)定常數的定量表達

簡單地分析一下表１中的數據，地注意到，Vina測定的復合物間的親和能ΔE與實驗測定的復合物的穩(wěn)定常數間有較好的相關性。同時我們認為，靜電相互作用可能對復合物穩(wěn)定性有很大的貢獻。因此，將Vina計算的復合物的親和能以及氨基酸側鏈的電子電荷指數作為主要貢獻項（具體數據詳見表1），對式（1）進行最小二乘擬合，得到表達復合物穩(wěn)定常數的定量表達式:

表1 與各配體相應的復合物穩(wěn)定常數、親和能，以及氨基酸側鏈電子電荷指數

根據式（2），考察各因素對復合物穩(wěn)定度的貢獻可以看到，A－LT的親和能最大，色氨酸甲酯的側鏈電子電荷指數最大，該復合物也最穩(wěn)定；A－LA、ALL和A－LV的親和能以及他們的配體側鏈電子電荷指數，都屬最小一檔，它們的穩(wěn)定度也在最差的一檔。

圖2表達的復合物穩(wěn)定常數計算值與實驗值間的關系，表明以Vina計算的親和能為核心參數，外加氨基酸側鏈的電子電荷指數，可以用于定量表達該手性咪唑環(huán)芳對一系列氨基酸甲酯的選擇性識別。

4 結論

（1）靜電吸引、空間互補和氫鍵作用等，是影響手性咪唑環(huán)芳對氨基酸甲酯選擇性識別的主要因素。為實現(xiàn)環(huán)芳對氨基酸甲酯的高度選擇性識別，應該至少從上述三個方面考慮，設計相應的環(huán)芳。

（2）在缺少準確的三維結構數據的情況下，用量子力學方法優(yōu)化單體結構，應用AutoDock Vina，將受體、配體對接，可以獲得復合物的最穩(wěn)定構象。

（3）根據親和能等數據，可以建立表達復合物穩(wěn)定常數的定量關系模型。這可能將有助于環(huán)芳在分子識別、人工酶模擬等方面的認識和應用。

圖2 復合物穩(wěn)定常數實驗值與計算值間的關系曲線

［1］Brown，C．J．；Farthing，A．C．，Preparation and Structure of Di－p－Xylylene．Nature 1949，164，915－916．

［2］Cram，D．J．；Steinberg，H．，Macro Rings．I．Preparation and Spectra of the Paracyclophanes．J．Am．Chem．Soc 1951，73，（12），5691－5704．

［3］Cram，D．J．；Cram，J．M．，Cyclophane chemistry:bent and battered benzene rings．Acc．Chem．Res 1971，4，（6），204－213．［4］Boekelheide，V．，Syntheses and properties of the［2n］cyclophanes．Top．Curr．Chem 1983，113，87－143．

［5］Cram，D．J．，The design of molecular hosts，guests，and their complexes．J．Incl．Phenom．1988，6，（4），397－413．

［6］Cram，D．J．；Cram，J．M．，Design of complexes between synthetic hosts and organic guests．Acc．Chem．Res 1978，11，（1），8－14．

［7］Newcomb，M．；Toner，J．L．；Helgeson，R．C．；et al．，Hostguest complexation．20．Chiral recognition in transport as a molecular basis for a catalytic resolving machine．J．Am．Chem．Soc 1979，101，（17），4941－4947．

［8］Yuan，Y．；Gao，G．；Jiang，Z．；et al．，Synthesis and selective anion recognition of imidazolium cyclophanes．Tetrahedron 2002，58，（44），8993－8999．

［9］Xu，X．－W．；Wang，X．－L．；Wu，A．－M．；et al．，Synthesis and molecular recognition of novel multiimidazole cyclophanes．J．Heterocyclic Chem．2009，46，（6），1137－1141．

［10］Camacho，D．H．；Salo，E．V．；Ziller，J．W．；et al．，Cyclophane－based highly active late－transition－metal catalysts for ethylene polymerization．Angew Chem Int Ed 2004，43，1．

［11］Poyatos，M．；Mata，J．；Peris，E．，Complexes with Poly（N－h(huán)eterocyclic carbene）Ligands:Structural Features and Catalytic Ap－plications．Chem．Rev．2009．

［12］Melaiye，A．；Sun，Z．；Hindi，K．；et al．，Silver（I）－midazole Cyclophane gem－Diol Complexes Encapsulated by Electrospun Tecophilic Nanofibers:Formation of Nanosilver Particles and Antimicrobial Activity．J．Am．Chem．Soc 2005，127，（7），2285－2291．

［13］Lombardino，J．G．；Wiseman，E．H．，Preparation and antiinflammatory activity of some nonacidic trisubstituted imidazoles．J．Med．Chem．1974，17，（11），1182－1188．

［14］Verras，A．；Kuntz，I．D．；Ortiz de Montellano，P．R．，Computer－assisted design of selective imidazole inhibitors for cytochrome p450 enzymes．J Med Chem 2004，47，（14），3572－9．

［15］Yan，J．M．；Jiang，Z．L．；Liu，W．Q．；et al．，Synthesis of Imidazole Cyclophane as Ligand of Binuclear Metalloprotein Models．Chinese Chem．Lett．1997，8，393－396．

［16］Garrison，J．C．；Simons，R．S．；Talley，J．M．；et al．，Synthesis and structural characterization of an imidazolium－linked cyclophane and the silver complex of an N－h(huán)eterocyclic carbene－linked cyclophane．Organometallics 2001，20，（7），1276－1278．

［17］You，J．；Yu，X．；Xiang，Q．；Lan，et al．，Novel chiral imidazole cyclophane receptors:synthesis and enantioselective recognition for amino acid derivatives．Chem．Commun．2001，2001，（18），1816－1817．

［18］Alcalde，E．；Alemany，M．；Gisbert，M．，A Direct Synthetic Approach to Novel Quadrupolar［14］Azolophanes．Tetrahedron 1996，52，（48）．［19］Hindi，K．M．；Siciliano，T．J．；Durmus，S．；et al．J．Med．Chem．2008，51，（6），1577－1583．

［20］Du，C．－P．；You，J．－S．；Yu，X．－Q．；et al．，Homochiral molecular tweezers as hosts for the highly enantioselective recognition of amino acid derivatives．Tetrahedron:Asymmetry 2003，14，3651－3656．

［21］Frisch，M．J．；Trucks，G．W．；Schlegel，H．B．；et al，Gaussian 03，Revision C．02；Gaussian，Inc．:Wallingford CT，2004．

［22］Trott，O．；Olson，A．J．，AutoDock Vina:Improving the Speed and Accuracy of Docking with a New Scoring Function，Efficient Optimization，and Multithreading．J Comput Chem 2009，31，（2），455－461．

［23］Alcalde，E．；Mesquida，N．；Pérez－García，L．；Hydrogen bonded driven anion binding by dicationic［14］imidazoliophanes．Chem．Commun．1999，（3），295－296．

［24］Murakami，Y．；Kikuchi，J．－i．；Hisaeda，et al．，Artificial Enzymes．Chem．Rev．1996，96，721－758．

［25］Langford，S．；Stoddart，J．，Self－assembly in chemical systems．Pure Appl．Chem．1996，68，1255－1260．

［26］Collantes，E．R．；III，W．J．D．，Amino Acid Side Chain Descriptors for Quantitative Structure－Activity Relationship Studies of Peptide Analogues．J．Med．Chem．1995，38，2705－2713．

Abstract:The structures of a chiral imidazole cyclophane acceptor and a series of amino acid methyl ester ligands were optimized by employing Gaussian 03 at the B3LYP with 6－31G（d）basis set．Then NBO charge distribution of the acceptor and ligands was calculated at the same level．The ligands were docked to the acceptor separately by using AutoDock Vina that resulted in the most stable conformations of complexes and affinity energy of the complexes were got consequently．A quantitative relationship between the experimentally determined stable constant and the calculated value including affinity energy and electronic charge index of side chain of amino acids was built．It was pointed out by discussion that electrostatic attraction，steric complementarity and hydrogen bond interaction are the main factors which mainly influence the recognition of chiral imidazole cyclophane for amino acid methyl esters．

Key words:imidazole cyclophane；molecular mechanics；density functional theory；molecular recognition；amino acid methyl ester

Exploration on the Molecular Recognition of a Chiral Imidazole Cyclophane for a Series of Amino Acid Methyl Esters

Lou Yongfeng1，2，Xiang Mingli1，Yang Mingli1
（1．National Key Laboratory of Biotherapy，West China Hospital，West China Medical School，Sichuan University，Chengdu 610041；2．Chemical Engineering college，Sichuan University，Chengdu 610065）