(a)	(b)
	图2 PM3优化得到的a-环糊精(a)和b-环糊精(b)的结构

        1995年，Bodor小组采用AM1方法对a-和b-环糊精的结构进行了详细的研究。结果表明，传统认为的环糊精具有C_n的对称性(对于a-、b-和g-环糊精，n分别为6，7和8)并不正确，晶体结构和分子力学优化给出的环糊精也并不处在其最优构象。相比之下，AM1优化得到的环糊精能量最低，因而被认为与真实构象最为接近。偶极分析指出，不同构象下的环糊精的偶极矩显著不同，而且随着不断优化，环糊精的偶极矩将逐渐地减小。比较AM1优化结果和晶体结构数据，发现优化没有显著地影响到环糊精的各个键长和键角，但是，明显地改变了环糊精的整体形状。在此基础上，Bodor等人又研究了2,6-二甲基-b-环糊精、2,6-二羟丙基-b-环糊精和2,3,6-三羟丙基-b-环糊精，并发现烷基化和羟丙基化并不会给环糊精带来大的空间位阻。
    在最近的研究中，我们采用AM1和PM3半经验方法研究了a-和b-环糊精的结构。结果发现AM1和PM3给出的环糊精在能量上都比晶体结构更加稳定。晶体中的结晶水被认为是导致这种差别的重要原因。因此，在研究自由状态下的环糊精与客体分子的包络时，使用优化后的结构还是很有必要的。进一步比较发现AM1优化的环糊精在整体形状上过于偏离圆形，使得部分用于稳定环糊精构象的分子内氢键受到了破坏。相比之下，PM3优化的环糊精很好地保持了这些氢键(图2)。这是由于PM3在处理氢键上较AM1更加可靠所致。所以，在目前的计算机条件下，PM3是能够用于环糊精结构研究的较好的量子化学方法。
    对于环糊精的理化性质，我们使用了两种不同级别(HF/STO-3G和HF/3-21G*)的基组，分别对晶体和PM3优化的a-和b-环糊精进行了单点计算。计算结果列入表1。从表1可以看出，PM3优化的环糊精的确比其晶体结构更加稳定。晶体中的环糊精具有很大的偶极矩，这可能来自于结晶水的影响。相比之下，PM3优化得到的环糊精偶极矩较低。前线轨道分析表明，环糊精在紫外-可见光谱区域中没有吸收，这与实验观测一致。值得一提的是，以上数据是国际上在环糊精领域中最先的从头算(Ab Initio)结果。
2 环糊精分子识别的驱动力
    主客体包络物形成的驱动力以及包络物的几何结构是超分子化学研究中的两个重要问题。其中，驱动力问题虽经多年研究，目前仍无定论。大多数的观点认为范德华作用、疏水作用以及静电作用等是环糊精分子识别的主要驱动力。但是，这些观点仍然停留在定性阶段^[13]。

        Chujo小组根据文献中的晶体数据，采用CNDO/2方法，计算了a-环糊精与对硝基苯酚、对羟基苯甲酸和苯甲酸等形成的包络物的能量和偶极矩。结果表明，在包络物中环糊精分子的偶极与客体分子的偶极几乎完全反平行。基于这一点，该小组认为，静电作用是环糊精分子识别的主要驱动力。但是，随后的研究发现并不是所有的环糊精包络物都采用这种偶极反平行的取向。苯甲醛就是其中一例^[15]。因此，以上结论不够全面。
    Starikov等人研究了b-环糊精和1,4-二羟基-2-丁炔之间分子识别的驱动力。分析表明，环糊精腔壁和客体分子之间可能存在有C-H · · · O和C-H · · · p作用。对于这些作用的大小，Starikov等人对参与作用的片段使用从头算方法进行了量子化学研究。结果显示，这些作用的能量大致落在0.7至1.1kcal/mol这一范围，远远小于常见氢键的能量，但是大于普通范德华作用的能量。
    最近，我们采用PM3和HF/3-21G*的方法，计算了a-环糊精与一系列等电子客体分子之间的分子识别。通过对计算结果进行密立根(Mulliken)电荷分析发现，在主客体包络物中环糊精和客体分子都带上了可观的电荷。进一步的分子轨道分析发现，环糊精与客体分子在形成超分子体系的过程中前线轨道发生了显著的重叠。对于中性苯的衍生物，环糊精总是作为电子的受体而成为路易斯(Lewis)酸，而客体分子总是作为电子的给体而成为路易斯碱。
    基于这一结果，我们发现由Morokuma提出的能量分解理论能够很好地说明环糊精分子识别的驱动力^[16]。根据该理论，在非共价超分子体系的形成过程中，各组分化合物的电子将不再属该化合物所有。因此，一个超分子的形成可以看成是四种作用的共同结果：
    静电作用(Electrostatic Interaction)。这种作用的本质是库仑相互作用。它产生于各组分分子自身的电荷分布，包括永久电荷、永久偶极以及永久多极等。这种作用既可以是相互吸引的，也可以是相互排斥的。通常较大的永久电荷和偶极将强化这种作用。
    极化作用(Polarization Interaction)。这种作用产生于各组分分子电荷分布的相互变形。它总是有利于超分子体系的形成。一般较大的体积和极化率将使这种作用增强。
    交换能(Exchange Energy)。它产生于各组分分子的满轨道之间的重叠， 其本质是泡利(Pauli)排斥。因此，这种作用总是不利的。
    电荷转移能(Charge-transfer Energy)。它产生于超分子体系中某个组分分子的满轨道与另一组分分子的空轨道之间的重叠。这种作用总是朝着能量降低的方向进行。它的大小与作用轨道之间的能隙密切联系。通常，这种作用中最强的项来自于某一组分的最高占据轨道(HOMO)和另一组分的最低空轨道(LUMO)。
根据以上分析可以看出，在超分子的形成过程中各组分分子之间尽管没有共价键的生成， 但是它们仍然是通过各自的分子轨道与分子轨道之间的相互作用而聚为一体的。除了能量上的差别外，这些分子轨道与分子轨道之间的相互作用与传统原子轨道相互作用并无二致。它们都要满足同样的量子力学原理。
    值得一提的是，在以上的四种作用中，前三种在环糊精分子识别的领域中已经被充分地认识和理解。然而，电荷转移作用作为一种非经典的量子化学效应，一直没有被人们看作是环糊精分子识别的驱动力之一。在我们的最新研究中，发现这种作用的贡献是不可忽视的。
3 环糊精包络物的结构
    量子化学一直是研究分子结构的有力手段。但是，由于超分子体系过于庞大，量子化学在超分子结构方面的研究还不广泛。
    Chujo等人曾经用CNDO/2方法研究了a-环糊精与对硝基苯酚和对羟基苯甲酸之间的分子识别。结果表明，客体分子的酚羟基处于环糊精腔内伯羟基端的取向在能量上不如硝基或羧基处于腔内伯羟基端的取向有利。其中对于a-环糊精与对硝基苯酚的包络物，两种取向之间的能量差别达到了141 kJ/mol，远远大于Mastui用分子力学方法计算得到的结果。基于这一点，Chujo等人认为偶极-偶极作用是决定环糊精包络物结构的主要因素，而相比之下范德华作用并不产生重要影响。这一研究虽然十分有趣，但是它也有不足之处。其一，这项研究采用的是文献中的晶体结构，而晶体中的晶胞排列及结晶水往往对包络物的结构产生不可忽视的影响。计算中仅仅考虑到环糊精和客体分子之间的作用是不够全面的。其二，对于客体分子的酚羟基处于环糊精腔内伯羟基端的研究简单地假设其苯环位置与酚羟基处于环糊精仲羟基端结构中苯环的位置一致。由于环糊精的包络物中主体分子的结构往往随客体分子结构的改变而改变，所以这种假设是不可靠的。总之，在环糊精包络物结构的研究中，简单地使用晶体结构而不真正地使用量子化学方法对包络物结构进行优化，得到的关于影响包络物稳定性的因素是不严格的。
    在随后的研究中，Chujo小组^[17]采用了“双层”模型(即人为地假设环糊精腔内外的介电常数不同)^[18]，使用CNDO/2方法模拟了环糊精与间硝基苯酚、苯甲酸和对羟基苯甲酸在水溶液中的分子识别。计算结果与¹³C NMR的实验结果一致。
    1996年，Botsi等人使用AM1方法计算了b-环糊精与1，7-二氧螺[5, 5]十一烷及壬醛之间的分子识别。研究中对包络物的结构进行了优化。结果表明，包络并不明显地影响到主体和客体分子的键长和键角。包络物也可以采用多种构象。对于长链的客体分子，它的碳链在环糊精腔内还可以适当地折叠。
    Huang等人使用AM1方法对环糊精的分子识别进行了一系列研究，发现g-环糊精与C₆₀之间按照1:1的比例形成的包络物是不稳定的。相反，两个g-环糊精分子仲羟基端面对面地靠近能够稳定地“笼”住一个C₆₀分子，形成2:1的包络物。
    对于a-和b-环糊精与苯酚以及苯甲酸之间的分子识别，Huang等人考虑了两种包络结构^[19]。在其中的一种结构中，苯环上的取代基被包络在环糊精腔内；而在另一种结构中，苯环上的取代基则被留在环糊精的腔外。AM1计算的结果表明：除了b-环糊精与苯酚的包络物以外，其余的包络物都采用的是苯环上的取代基被包络在环糊精腔内的结构。这与实验观察一致。分析显示，包络物的稳定与否取决于包络前后主客体分子的结构是否发生明显的改变。在进一步对b-环糊精与取代苯甲酸的包络进行计算后，Huang等人又发现环糊精包络物采取的构型与包络物中能够形成的分子内以及分子间氢键数目有关。
    Balabai等人使用AM1方法研究了b-环糊精与9-羟基-3-异吩噁唑酮之间的分子识别。研究中使用刚性的分子结构，采用模拟退火(Simulated Annealing)算法得到了最优的包络物构型。发现b-环糊精只部分包络了9-羟基-3-异吩噁唑酮分子。这一结构与实验观测能够吻合。

	(a)		(b)
	图3 PM3描述的a-环糊精与苯乙酮之间的分子识别。(a) 乙酰基处于环糊精的伯羟基端； (b) 乙酰基处于环糊精的仲羟基端。

        最近，我们研究了a-环糊精与苯乙酮之间的分子识别^[20]，用PM3方法模拟了客体分子一步步插入环糊精腔内的全过程，并在模拟的每一步中对包络物的结构进行了完全非限制性的优化。计算结果表明，对于该包络物乙酰基被包络在环糊精腔内的结构在能量上不如乙酰基被留在环糊精腔外的结构稳定(图3)。这一结果与Chujo等人认为的偶极相反的包络结构不同，但是与最新的实验观察结果相一致。分析表明，在乙酰基被包络在环糊精腔内的取向中，乙酰基上的氢原子与环糊精的内壁过于接近，导致了十分强烈的范德华排斥作用。这就使得这种取向在能量上不利。比较包络前后的几何参数，发现在包络物中环糊精和客体分子的结构都发生了显著的改变。这可以看成是“诱导契合”概念的一个实例。它也表明传统计算中假设环糊精和客体分子的结构不受包络的影响这一观点是不合适的。
4 结论和展望
    环糊精作为一种多功能分子，在许多领域中受到了越来越多的关注。环糊精化学的发展也将为人们理解酶与底物，受体与配体的相互作用提供帮助。量子化学作为一种大有前途的研究手段，将会在环糊精化学中得到日益广泛的应用。
    在目前计算机资源的条件下，PM3半经验方法是用于环糊精化学研究的良好的量子化学工具。它提供的结果的确能够帮助我们理解环糊精分子识别中的一些深层次问题。当然，从头算的方法在环糊精化学中的应用还有待于大力发展。尤其是在目前的条件下，在从头算的基础上，使用小的模型化合物来模拟环糊精分子识别的某些细节还是很有价值的。
    值得指出的是，迄今为止在环糊精领域的量化研究仍然被限制在真空气态的环境中(in vacuum)。因此，这些计算无法反映出疏水作用对环糊精分子识别的贡献。 鉴于疏水作用在超分子识别中的重要地位，发展出一套能够用于水溶液环境中超分子体系的量子化学计算方法必将具有重大的价值。
5 参考文献
[1] Saenger, W., Angew. Chem. Int. Ed. Engl., 1980, 19: 344.
[2] Szejtli, J., Chem. Rev., 1998, 98: 1743.
[3] Breslow, R., Dong, S. D. Chem. Rev., 1998, 98: 1997.
[4] Takahashi, K., Chem. Rev., 1998, 98: 2013.
[5] Hedges, A. R., Chem. Rev., 1998, 98, 2035.
[6] Uekama, K., Hirayama, F., Irie, T., Chem. Rev., 1998, 98: 2045.
[7] Connors, K. A., Chem. Rev., 1997, 98: 1325.
[8] Rekharsky, M. V., Inoue, Y., Chem. Rev., 1998, 98: 1875.
[9] Schneider, H.-J., Hacket, F., Rudiger, V., et al, Chem. Rev., 1998, 98: 1755.
[10]  Harata, K. Chem. Rev., 1998, 98: 1803.
[11] Guo, Q.-X., Liu, L., Cai, W.-S., et al, Chem. Phys. Lett., 1998, 290: 514.
[12] Liu, L., Quo, Q.-X., J. Chem. Inf. Comput. Sci., 1999, 39: 133.
[13] Liu, L., Guo, Q.-X., J. Phys. Chem. B, 1999, 103: 3461.
[14] Lipkowitz, K. B., Chem. Rev., 1998, 98: 1829.
[15] Rao, V. P., Turro, N. J., Tetrahedron Lett., 1989, 35: 4641.
[16] Morokuma, K., Acc. Chem. Res., 1977, 10: 294.
[17] Sakurai, M., Hoshi, H., Inoue, Y., et al, Chem. Phys. Lett., 1989, 163: 217.
[18] Inoue, Y., Hoshi, H., Sakurai, M., et al, R. J. Am. Chem. Soc., 1985, 107: 2319.
[19] Huang, M.-J., Watts, J. D., Bodor, N., Int. J. Quantum Chem., 1997, 64: 711.
[20] Li, X.-S., Liu, L., Guo, Q.-X., et al, Chem. Phys. Lett., 1999, 307, 117.

---

刘磊 男，22岁，安徽滁州人，学士学位，现美国哥伦比亚大学化学系研究生。**联系人
国家自然科学基金资助。1999-06-20收稿，1999-10-08修回。