1、 氨基酸在冠醚手性固定相上定量结构 对映体保留相关性的研究作者：施介华 肖科科 吕园园【摘要】 运用量子化学中的 HartreeFock 程序(631G 基组)方法计算质子化 氨基酸溶质的分子结构参数，借助于多元线性回归法建立了 氨基酸对映体在冠醚手性固定相上的色谱保留与其分子结构参数之间的定量结构 对映异构体保留(QSERR)模型。结果表明， 氨基酸光学异构体的容量因子对数(logk)与质子化 氨基酸溶质的分子结构描述参数之间具有较好的线性相关性。在 QSERR模型中，溶质结构描述参数 EHOMO， DIP， ElcE， Ang 和 logP具有较为明确的物理意义，这些参数反映了固定相与溶质
2、分子之间的静电作用、 作用力、色散力、立体位阻和疏水作用，能较好地解释 氨基酸对映体在联萘冠醚 CSP上保留机理。建立的 QSERR模型具有较好的稳定性和预测能力。 【关键词】 定量结构 对映异构体保留关系， 氨基酸， 冠醚手性固定相， 色谱保留机理， HartreeFock 程序1 引 言运用高效液相色谱(HPLC) 分离对映异构体已成为现代合成化学、生物医药学及农业化学等领域中常用的分析方法。大部分色谱手性分离很大程度上依赖于分析测试者的经验。建立具有明确物理含义的定量结构 对映体保留关系(QSERR)能有效地预测色谱保留值，可以为选择最佳分离条件提供理论依据，还有助于理解色谱保留及手性识
3、别的机理14。近年来, 随着计算机技术和量子化学方法的不断发展, 量化计算已越来越受到人们的关注57。在 QSRR/ QSERR研究中，独自应用量化参数或将量化参数与其它常规参数联合已经越来越多地被使用。文献8，9应用量化参数建立了 QSERR模型，研究了 氨基膦酸酯类化合物在两种不同的衍生化环糊精固定相上手性识别机理; Suzuki等10通过 2D和 3D分子描述参数建立了一系列芳基醇类化合物在 4种刷型固定相上的QSERR模型，指出影响容量因子(k)和分离因子()的分子描述参数不同;施介华等6,11运用量化参数分别建立烃类和酯类化合物的色谱保留关系的 QSRR模型。冠醚作为一类色谱手性固定
4、相，常被用来识别具有 NH2官能团的对映体1216，其独特的冠环结构对客体分子的选择性识别起着重要作用。但有关手性冠醚类固定相上 氨基酸的色谱保留和对映异构体结构间的定量构效关系的报道鲜见。本研究运用量子化学方法中的HartreeFock 程序(631G 基组)计算了 20个质子化 氨基酸溶质分子的结构描述参数，通过多元线性回归建立了 氨基酸类化合物在键合联萘冠醚手性固定相(CSP)上的定量结构与对映体保留的 QSERR模型，探讨了 氨基酸对映体在键合联萘冠醚固定相上的色谱保留机理。2 实验部分2.1 数据来源本研究所用 氨基酸化合物在新型联萘冠醚手性柱上的色谱保留数据取自文献15(见表 1)
5、。2.2 溶质分子描述参数的计算及 QSERR模型建立由于 氨基酸的 NH2基团在酸性流动相中将发生质子化，形成强极性的铵离子(NH+4)，实际上色谱分离过程中与冠醚类手性固定相发生相互作用的是质子化的氨基酸16(见图 1b)。所有质子化表 1 氨基酸在联萘冠醚手性固定相上保留值15Table 1 Retention value of amino acids on the crown chiral stationary phase15No.氨基酸 Amino acidlogk1logk2No.氨基酸 Amino acidlogk1logk21 丙氨酸 Alanine -0.23660.3385
6、11甲硫氨酸 Methionine-0.11350.53912精氨酸 Arginine-0.2366-0.045812苯基丙氨酸 Phenylalanine-0.10240.37663天冬酰胺 Asparagine-0.5528-0.337213苯甘氨酸 Phenylglycine0.15490.59884天冬氨酸 Aspartic acid-0.37680.008614脯氨酸 Proline-0.6383-0.63835半胱胺酸 Cysteine-0.4202-0.013215丝氨酸 Serine-0.5686-0.30107谷氨酸 Glutanic acid-0.22910.541616苏
7、氨酸 Threonine-0.6198-0.19386谷氨酰胺 Glutamine-0.44370.245517 色氨酸 Tryptophan0.14610.60538组胺酸 Histidine-0.4815-0.309818酪胺酸 Tyrosine-0.20760.28109异亮胺酸 Isoleucine-0.5086-0.236619缬氨酸 Valine-0.5528-0.284010亮氨酸 Leucine-0.19380.4997204 羟基苯甘氨酸4Hydroxyphenylglycine-0.17390.6730k1: L 氨基酸容量因子(Capacity factor of LAm
8、ino acid)， k2: D 氨基酸容量因子(Capacity factor of DAmino acid)。的 氨基酸分子结构首先用ChenDraw作出平面图，然后用 Chem3D中的 MOPAC 8.00对分子几何结构进行初步优化后，再用 Gaussian 03软件包中的 HartreeFock 程序(631G 基组)进行结构优化，振动分析计算的结果中无虚频，证明优化得到的分子结构对应于能量极小点。获取质子化氨基酸溶质的量子化学参数：DIP(分子总体偶极距)、 TE(分子总能量)、 ELUMO(最低分子空轨道的能量)、 EHOMO(最高分子占有轨道的能量)、 Ang(原子 NC*C之间
9、平面角)、 MR(分子折射率)、 logP(疏水性参数)和 ElcE等分子描述参数在 Chem3D中计算得到。所有量化参数的计算和分子模拟在 PD2.80PC机上运行完成。3 结果与讨论3.1 QSERR 模型的建立由于分子“手性”的差别，手性分子与手性固定相相互作用的过程中存在着一定差异，即手性分子的结构参数与色谱参数之间存在一定相关性这一假设是完全合理的6。在 QSERR模型建立过程中，最大的困难是如何找到能够区分不同对映体的手性分子结构参数，因为对映体分子诸多性质参数(如分子尺寸、形状、电性参数和疏水性参数等)有时很接近，甚至完全相同。与许多手性分子一样， 氨基酸对映体分子结构参数大部分
10、相同。本研究通过 Gaussian 03和 Chem3D软件计算，提取了质子化后的 L 型和 D 型 氨基酸的分子结构描述参数(见表 2和表 3)。为保证所选取的结构参数与容量因子对数(logk)之间有较好的相关性，尽量选择与容量因子相关性较大，彼此之间相关性较小的参数。本研究首先对各参数进行相关性分析，发现 MR与 ElcE和 TE之间有明显相关性，因此分别选取 MR， ElcE和 TE与其它结构描述参数进行多元线性回归，得到 氨基酸对映体在键合联萘冠醚 CSP上的色谱保留参数的 QSERR模型为：logk1=-13.400-0.402EHOMO-0.047DIP+0.115Ang-5.90
11、10-5ElcE(n=20, r=0.959, SD=0.067)(1)表 2 质子化 L 基酸分子参数Table 2 Molecular descriptors of Lamino acid cationsNo.DIPElcEEHOMOELUMOTElogPMRAng 13.8840-4882.0-1.42461.8078-1340.9-1.4091.9489108.4722.1710-12247.4-1.32411.5128-2447.1-2.7254.2690105.3133.0044-8496.3-1.08701.4144-2009.9-2.9232.8171107.3042.4240
12、-8797.0-1.38020.9798-2110.0-2.2812.6015106.8252.0069-5938.5-1.60210.7928-1534.8-1.3062.7552106.2273.7806-12972.2-1.49831.1399-2266.0-1.8253.0653105.4065.1110-9873.5-0.76881.2623-2165.3-2.4673.2809106.6284.0373-10590.2-0.70971.1597-2164.1-2.2353.6785106.3592.4120-12247.4-1.32411.5128-2447.1-2.7254.26
13、90105.31104.2044-8616.1-1.25371.6917-1808.1-0.2933.3403105.32113.1407-11670.7-1.46750.6253-1846.4-0.8603.9776105.20122.8048-11314.1-1.45410.5860-2163.30.0664.4601105.48133.8687-9734.0-1.89341.0413-2007.8-0.3903.9963105.15142.1267-6106.1-0.70761.5282-1624.3-0.7332.7878106.30152.2410-6260.0-1.10271.56
14、18-1661.3-2.5742.1020107.88162.2684-7661.3-1.07131.5946-1816.9-2.3162.5658107.80172.0277-16153.7-1.56800.2930-2658.50.0125.5776105.32182.1460-11972.2-1.43460.5580-2483.80.0125.5776105.47192.7089-7398.3-1.35201.7728-1652.2-0.7492.8765104.00204.8453-11317.9-1.67890.8915-2328.3-0.7794.1494104.95氨基酸编号同表
15、 1(The number of aminoacids are the same as in Table 1)。 DIP: Dipole moment (debye); ElcE: Electronic energy (eV); EHOMO: Highest occupied molecular orbital (eV); ELUMO: Lowest unoccupied molecular orbital(eV); TE: Total energy(eV); logP: Hydrophobic parameters; MR: Molar refractivity; Ang: Plane an
16、gle between NC*C。表 3 质子化 D 氨基酸分子参数式中, n为线性回归样本的个数， r 为拟合方程的相关系数，SD 为标准偏差。由式(1)和(2)可知，所建立的 氨基酸对映体在键合联萘冠醚 CSP上的色谱保留参数的QSERR模型具有较好的多元线性相关性。同时，为了识别异常值，对参与回归样本的保留值进行去一法检验，未出现异常值。这表明此 QSERR模型具有较好的稳定性。3.2 手性识别机理的探讨根据液相色谱分离原理，溶质在色谱柱上的保留行为主要由溶质与固定相以及流动相之间的分子作用力所决定的。所建立的 氨基酸对映体在键合联萘冠醚 CSP上的色谱保留参数的 QSERR模型中，每个
17、分子描述参数均表征色谱过程中溶质同固定相或流动相的某种相互作用。在液相色谱中 ELUMO和 EHOMO表示溶质和固定相之间电荷的传递反应，也就是在溶质和固定相之间存在氢键作用力或 相互作用17。logP与溶质疏水性相关，体现了溶质分子与固定相之间的疏水作用。Ang表示对映体色谱保留过程中对分子构象的要求9。DIP 可以表示溶质分子和固定相之间的静电相互作用能力的大小2。MR 反映了溶质分子参与分子色散力的能力18，而 ElcE和 TE与 MR具有较好的相关性，因此ElcE和 TE包含一定色散力的信息,也可以表示溶质分子参与色散力的能力。尽管对分子间力(包括氢键力、范德华力、空间位阻、疏水作用等
18、)的研究己比较深入、全面，但是如何协同这些作用力而增加或减弱非对映体络合物的稳定性却并不很清楚。从所建立的 QSERR模型可知，在联萘冠醚 CSP上，质子化氨基酸的保留行为主要取决于固定相与溶质分子之间的静电作用、 作用力、色散力、立体位阻和疏水作用。3.3 保留值的预测所建立的 QSERR模型可预测 氨基酸对映体在联萘冠醚手 CSP上的色谱保留值，即通过计算 氨基酸对映体的结构描述参数来预测其对映体的色谱保留值。结果表明，在这种手性冠醚固定相上， 氨基酸对映体的容量因子对数(logk)的预测值与实验值之间具有较好的线性相关性(图 2)，其线性回归方程分别为：logk1(pred)=0.920
19、1logk1(exp)+0.0069(n=20, r=0.9572, SD=0.064)(3)logk2(pred)=0.9199 logk2(exp)+0.0161(n=20, r=0.9566, SD=0.064)(4)两者回归方程的斜率接近于 1，基本通过原点，这表明所建立的QSERR模型用于预测 氨基酸对映体在这两种手性冠醚固定相上的保留行为是切实可行的。图 2 预测值和测定值Fig.2 Plot of experimental vs. predicted logk values3.4 小结本研究建立了 氨基酸对映体在键合联萘冠醚 CSP上的色谱保留参数的 QSERR模型，各结构描述参
20、数具有明确的物理含义。能较好地解释 氨基酸对映体在联萘冠醚 CSP上保留机理。所建模型在预测其保留值时具有较好的准确性和稳定性。【参考文献】1 Jiang Feng(江 峰), Shi JieHua( 施介华), Yan Wei(严 巍), Li Bing(李 兵). Chinese J. Anal. Chem.(分析化学), 2004, 32(11): 155015552 Gao RuYu( 高如瑜), Wang HeFang( 王荷芳). Chinese Journal of Organic Chemistry(有机化学)， 1999, 19(2): 1351403 Lipkowitz K
21、 B. J. Chromatogr. A, 1994, 666: 4935034 Lipkowitz K B. J. Chromatogr. A, 1995, 694: 15375 Ouyang Liang(欧阳 亮), He Gu(何 谷), Guo Li(郭 丽). Acta Chim. Sinica(化学学报), 2006, 64:137913846 Shi JieHua( 施介华), Jiang Feng(江 峰), Yan Wei(严 巍), Wang JiKang (王纪康), Wang GuiLin( 王桂林). Chinese J. Anal. Chem.(分析化学), 200
22、5, 33(2): 1811867 Zeng XiaoLan( 曾小兰), Liu HongXia( 柳红霞), Liu HongYan( 刘红艳). Acta Chim. Sinica(化学学报), 2007, 65: 17978 Huang JunMin( 黄君珉), Chen Hui(陈 慧), Gao RuYu( 高如瑜), Wang QinSun( 王琴孙). Chem. J. Chinese Universities(高等学校化学学报), 2001, 22(11): 183818429 Chen Hui(陈 慧 ), L XianYu( 吕宪禹), Huang JunMin( 黄君
23、珉), Gao RuYu( 高如瑜), Wang QinSun( 王琴孙). Chem. J. Chinese Universities(高等学校化学学报), 2000, 21(4): 56256510 Suzuki T, Timofei S, Iuoras B E, Uray G, Verdino P, Fabian Walter M F. J. Chromatogr. A, 2001, 922(12): 132311 Shi JieHua( 施介华), Jiang Feng(江 峰), Wang GuiLin( 王桂林), Yan Wei(严 巍). J. Zhejiang Univers
24、ity of Technology(浙江工业大学学报), 2004, 32(4): 41842212 Chen JuanJuan( 陈娟娟), Zhang DaTong( 张大同), Shen BaoChun( 沈报春), Zhang XueJun( 张雪君), Xu BeiJia( 徐贝佳), Xu XiuZhu( 徐秀珠). Chinese J. Anal. Chem.(分析化学), 2006, 34(11): 1535154013 Machida Y, Nishi H, Nakamura K, Nakai H, Sato T. J. Chromatogr. A, 1998, 805(12
25、): 859214 Hyun M H, Jin J S, Koo H J, Lee W. J. Chromatogr. A, 1999, 837(12): 758215 Hyun M H, Han S C, Lipshutz B H, Shin Y J, Welch C J. J. Chromatogr. A, 2001, 910(2): 35936516 Sousa L R, Sogah G. D Y, Hoffman D H, Cram D J. J. Am. Chem. Soc., 1978, 100(14): 4569457617 Matlin S A， Tiritan M E, Crawford A J, Cass Q B, Boyd D R. Chirality, 1994, 6： 13518 Wang Peng(王 鹏). Method of study quantitative structureactivityrelationships( 定量构效关系及研究方法), Harbin(哈尔滨)： Harbin Ind. Tech. Press(哈尔滨工业大学出版社), 2004: 118
