1、收稿日期：基金项目：吉林省自然科学基金（），白城医学高等专科学校科研平台重点课题（），白城医学高等专科学校高层次人才团队建设基金（），教育部国家级大学生创新基金（）和吉林省白城市科研规划课题（）资助项目作者简介：柳国洪（），女，吉林白城人，教授，主要从事手性药物化学研究：通信作者：潘 宇（），男，吉林大安人，副教授，主要从事计算机应用化学研究：赵丽红（），女，吉林大安人，教授，主要从事手性药物化学研究：柳国洪，彭国强，张栩宾，等 水液相下赖氨酸钙（）配合物旋光异构的 研究 江西师范大学学报（自然科学版），（）：，（）（），（）：文章编号：（）水液相下赖氨酸钙（）配合物旋光异构的 研究柳国洪，彭
2、国强，张栩宾，郝成欣，刘 芳，王佐成，潘 宇，赵丽红（白城医学高等专科学校基础医学院，吉林 白城；白城师范学院理论计算中心，吉林 白城）摘要：该文采用 的 和 方法结合极化介质的 模型，研究了在水液相下两性赖氨酸分子 价钙配合物（）的旋光异构 反应通道研究发现：（）的旋光异构可在 从两性异构为中性后，以氨基 作桥和羰基 作桥迁移的 个通道上实现 势能面计算结果表明：，以氨基 作桥的反应通道最具优势，在隐性水溶剂效应下决速步能垒为 ，从手性 向氨基 迁移的过渡态在显性水溶剂效应下的能垒降至 左右 研究结果表明：在水液相下手性（）的消旋过程十分缓慢，其用于生命体同补 和（）具有较好的安全性关键词：
3、赖氨酸；价钙；密度泛函理论；旋光异构；自由能垒中图分类号：文献标志码：引言 是生命体必需的手性氨基酸 根据构型不同分为 和，根据旋光性分为左旋体（）和右旋体（）在人体内有生物活性，能促进人体发育、提高免疫力、增强胃液分泌和骨髓造血机能等，无活性 钙是生命体重要的元素，在体内是 价态，其对骨骼形成、保证心脏规律跳动有重要作用 传统的补钙产品有醋酸钙、碳酸钙和葡萄糖酸钙，临床试验表明服用这些补钙产品或会导致肾结石，或会产生胃肠道不良反应在食物中的钙被人体吸收是通过 在小肠中与氨基酸螯合后，氨基酸钙被人体整体吸收进行的 赖氨酸钙用于补钙更易被生命体吸收，同时还可补充 手性 的金属配合物也有手性，手性
4、药的作用是药物分子通过与体内大分子严格地手性匹配实现的 手性不同的药物分子在活性、毒性和代谢过程方面差别很大，通常一个有效，另一个有负作用 如沙利度胺可止吐和镇静，而 沙利度胺则致畸，沙利度胺的消旋导致了“沙利度胺事件”因为优构体能否向劣构体转化以及转化速率均会严重影响手性药的疗效，所以手性分子旋光异构机理的研究对在临床上安全使用手性药具有重要意义 基于手性分子消旋反应的重要性，人们对 的旋光异构进行了广泛研究 文献采用密度泛函理论研究了气相 的旋光异构以及水分子、纳米孔道（分子筛及碳纳米管）和水溶剂的催化作用，结果表明：气相 不能消旋，纳米孔道有限域催化作用，水分子（簇）及水溶剂有较好的催化
5、作用水溶剂可能改变反应的机理，生命体是富水环境，但是关于在水液相下（）旋光异构的研究鲜见文献报道 为揭示（）用于生第 卷 第 期 江西师范大学学报（自然科学版）年 月 （）命体同补 和 是否安全，基于文献的研究经验，本文对水液相下赖氨酸钙（）配合物旋光异构的机理进行了研究 计算方法本文采用密度泛函理论的对主族元素分子间弱作用有较好表现的 方法，结合自洽反应场理论的 方法，采用 （，）基组优化反应势能面的驻点，对过渡态进行（内禀反应坐标）计算，考察过渡态的可靠性 为计算出高水平的势能面及配位键能，采用密度泛函理论具有较高精度的 方法，在高角动量弥散基组 （，）水平上计算单点能，然后计算相关体系的
6、自由能变和焓变 采用在 程序中的自然键轨道理论计算物种的电荷，采用 程序计算成键和成环临界点（和）的电荷密度（和）以及电荷密度的拉普拉斯值（），分析关键物种的成键特征 的 配合物记作（），其在 通道上的第 个无手性过渡态记作，中间体记作，其在 通道上旋光异构的最终产物记作（）；（）与 氢键作用的体系记作（）；个 与 的 配位，同时（）与 作用的体系记作（）；其他体系表示方法与此相似 有关计算工作采用 程序进行 结果与讨论文献研究表明在水液相下两性氨基酸的构象稳定 优化的两性 与 配合物的手性对映体如图 所示 （）（）（）图 在水液相下（）的手性对映体构象 在水液相下（）的旋光异构是在极性水溶剂
7、和水分子（簇）的共同作用下进行的 在水溶剂中的水分子有单分子和水簇，与底物弱作用没参与反应的（）（为构成水簇的水分子个数，）对反应能垒影响较小，与金属配位对非 迁移能垒影响也很小 为节省资源，本文对非 迁移过程不考虑（）的作用 先讨论在隐性水溶剂效应（水溶剂的极性作用）下（）的异构，然后讨论在显性水溶剂效应（水溶剂的极性与（）共同作用）下（）旋光异构决速步的反应过程，对重要的 迁移反应进行详细讨论，对其他基元反应做一般讨论 在隐性水溶剂效应下（）的旋光异构（）的旋光异构可在两性 向中性异构后 向 迁移或向 迁移的 通道和 通道上实现，异构过程如图 和图 所示，势能面如图 所示 通道 第 基元为
8、（）经质子在和 之间迁移的过渡态 沿过渡态矢量的正方向从 迁移到，构型异构成 电荷计算结果表明：从（）经 到，的电荷量从 变为，再变为，该基元反应是质子迁移过程（电荷计算结果表明本文的 迁移均为质子迁移，下同）从（）到，键键长从 拉伸至 ，从 减小至 ，其 从 变为 ，共价键减弱；的 从 增加到 ，其 从 变为，过渡态的 成环，存在大 键作用；共价键键长从 小幅拉伸至 ；配位键键长从 小幅拉伸至 个共价键的拉伸断裂以及 个共价键和配位键的小幅拉伸需要一些能量，由于过渡态成环较稳定，所以产生的内禀能垒（）并不高第 基元为 经 内旋转过渡态 羧羟基 绕 旋转，从羧基外侧转到内侧，构象异构成 在从
9、到 的过程中，键俯视逆时针旋转 ，化学键内旋转所需能量不多，产生了 的内禀能垒第 基元为 经 在 和 间迁移的过渡态 沿过渡态矢量的正向从 迁移到氨基，构型异构成中间体 从 到，键键长从 拉伸至 ，从第 期柳国洪，等：水液相下赖氨酸钙（）配合物旋光异构的 研究 减小至 ，其 从 变为，共价键断裂；键键长从 拉伸至 ，共 价 键 减 弱；骨 架 二 面 角从 增加至 ，键右视逆时针旋转 ，产生的内禀能垒为 第基元为经氢迁移的过渡态 在纸面里沿过渡态矢量的负向从 迁移到，异构成，（）旋光异构 从 到，键键长从 拉伸至 ，从 减小至 ，其 从 变为 ，共价键减弱；键键长从 拉伸至 ，键键长从 拉伸至
10、 ，共价键减弱，产生了 的内禀能垒 （）后视图 正视图 （）图 在隐性水溶剂效应下（）在 通道上的旋光异构过程江西师范大学学报（自然科学版）年 （）后视图 （）（）正视图 正视图 （同）下接 图 在隐性水溶剂效应下（）在 通道上的旋光异构过程图 （）旋光异构反应的自由能势能面 第 基元为 经氨基氢 和 在纸面内外翻转的过渡态 和 从纸面内侧翻转到外侧（后视图），构象异构成 在从 到 的过程中二面角 从 增加至 ，非骨架形变所需能量较少，产生的内禀能垒为 第 基元为 经 内旋转过渡态 羧羟基 绕 俯视逆时针旋转，从羧基内侧转到外侧，构象异构成 在从 到的过程中 键俯视逆时针旋转 ，第 期柳国洪，
11、等：水液相下赖氨酸钙（）配合物旋光异构的 研究产生了 的内禀能垒 该能垒比相同的 键逆时针旋转过渡态 产生的内禀能垒（）小很多，其原因是从 到 的过程中 的运动方向与体系偶极矩方向相反，体系的电场力对 绕 的俯视逆时针旋转做正功第 基元为 经质子在 和 之间迁移的过渡态 沿过渡态矢量的负向从 迁移到，构型异构成 从 到，键键长从 拉伸至 ，及计算结果均表明在此过程中共价键有少许减弱 的 的 为 ，过渡态的 成五元环，较稳定 个共价键的少许拉伸，且过渡态较稳定，再加上 从 向 迁移过程是体系的电场力做正功（此时的迁移方向与体系的偶极矩矢量相反），产生的内禀能垒只有 第 基元为 经螯合环开合（成键
12、与断裂）的过渡态 螯合环闭合，成键，构象异构成（），得到了（）旋光异构的稳定产物 在从 到 的过程中，只是键角 小幅改变，所需能量很少 产生的电子能垒只有（见势能面），吉布斯自由能经校正后变为负势垒（，没在势能面上标识），该基元反应可视为无势垒过程 其原因是过渡态有一个虚频，不计入热力学贡献，的自由能校正量比 小 ，加上该基元电子能量对应的势垒只有，导致 比 的自由能低 因此，该基元的过渡态 只具有理论意义，向（）的异构是无势垒过程 通道 第 基元为（）经 在 和羰基 间迁移的过渡态 从 迁移到，异构成 从（）到，键键长从 拉伸至 ，从 减小至 ，其 从 变为 ，共价键大幅减弱；的从 增加至
13、，其 从 变为，（）成环；键键长从 拉伸至 ，共价键减弱；键键长从 拉伸至 ，配位键减弱；二面角微增 ，键右视顺时针旋转 ，产生了 的内禀能垒 接下来 通过 路径和 路径向（）异构）路径 第 基元为 经质子在羰基 和 间迁移的过渡态 沿过渡态矢量的正向从迁移到，异构至产物，在 通道 路径上实现旋光异构 结构分析表明 同于，记作（）从 到，键键长从 拉伸至 ，和 计算结果均表明共价键减弱；的 和计算结果均表明 成环；键键长从 拉伸至 ，共价键减弱；键键长从 拉伸至 ，配位键减弱，产生了 的内禀能垒）路径为第 基元 经质子迁移的过渡态 在纸面外（后视图）从 迁移到，异构成，（）旋光异构 从 到，键
14、键长从 拉伸至 ，计算结果表明 减小，共价键明显减弱；键键长从 拉伸至 ，键键长从 拉伸至 ，个共价键少许减弱；键键长从 拉伸至 ，配位键减弱，产生了 的内禀能垒第 基元为（）经氨基氢 和 左右翻转的过渡态 和 从左侧翻转到右侧（后视图），构象异构成 从（）到，非骨架二面角 从 变为 ，所需能量较少，产生的内禀能垒为 第 基元为 经 旋转的过渡态 键右视逆时针旋转 ，构象异构成，产生的能垒为 结构分析表明 同于，其接下来异构得到产物（），记作（）（）从图 可以看出：在隐性水溶剂效应下（）在 通道上异构具有优势，决速步骤是第 基元，能垒为 该能垒比极限反应能垒（）还高很多，这说明在隐性水溶剂效应
15、下（）不能消旋 通道比 通道的决速步能垒低，其原因主要有 个：）比（）的 键活化程度更高（前者键长长且红外振动频率低（分别为 和 ）；）通道的反应物比 通道从反应物到过渡态 键的键长拉伸幅度更小正负反应能垒在 以下时反应物和产物可以共存，因此从图 还可以看出：在 通道上的产物是、和（）共存，只是（）的分布更高江西师范大学学报（自然科学版）年 在显性水溶剂效应下（）旋光异构的决速步骤结合图 和图 可以看出：在隐性水溶剂效应下（）在 通道和 通道上旋光异构的决速步基元反应过程分别是 和（）在显性水溶剂效应下，水分子（簇）除了与 和（）的钙原子配位外，同时还与 和（）有氢键作用 当单个水分子或水簇与
16、钙离子配位以及 聚水和 聚水做 迁移媒介时，质子在 和羰基 以及 和氨基 间迁移的能垒相差很小，为使问题简便并节省篇幅，只考虑单个水分子与钙离子配位，做 迁移媒介的水簇只考虑 聚水情况（）最高配位数是，配位键键能计算结果表明：（）和 体系的 与（）的配位键能（严格地说还含有分子间的氢键能）在配位水分子个数 时基本相同，均明显大于在 时的配位键能，因此物种（）和 在 时的分布最高（前者 为 ，后者 为 ）（见表）它们在 聚水作用下的质子转移反应历程如图 所示，反应势能面如图 所示 为节省篇幅，对 个配位水的情况做详细讨论，对满配的情况（和 配位、（）和 配位）做一般讨论，对其他无配位水及存在、和 个配位水的情况只给出反应势能面表 和 与水分子的配位键能物种焓校正 高水平单点能 与底物 配位的（）个 与底物 的配位键能（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）第 期柳国洪，等：水液相下赖氨酸钙（）配合物旋光异构的 研究 （）左视图 （）（）（）左视图 （）（）通道（）（）俯视图 （）（）（）（）俯视图 （）（）通道图 在显性水溶剂效
