气相脯氨酸Ca（Ⅱ）配合物的旋光异构及水分子（簇）的催化机制

气相脯氨酸Ca（Ⅱ）配合物的旋光异构及水分子（簇）的催化机制 Optical Isomerism of Proline Ca（Ⅱ） Complex in Gas Phase and Catalytic Mechanism of Water Molecules（Clusters） 1 first-author 黄煜纯，女，实验师，硕士，现从事计算机应用化学研究。E-mail：huangyc@ccsfu.edu.cn 黄煜纯，女，实验师，硕士，现从事计算机应用化学研究。E-mail：huangyc@ccsfu.edu.cn huangyc@ccsfu.edu.cn 2 2 2 1 2 corresp E-mail： gaof317@163.com E-mail： gaof317@163.com gaof317@163.com 2 corresp wangzc188@163.com wangzc188@163.com wangzc188@163.com 长春师范大学 网络中心，吉林 长春 130000 长春师范大学 网络中心，吉林 长春 130000 Network Center， Changchun Normal University， Changchun 130000， Jilin， China Network Center， Changchun Normal University， Changchun 130000， Jilin， China 白城师范学院 物理学院，吉林 白城 137000 白城师范学院 物理学院，吉林 白城 137000 College of Physics， Baicheng Normal University， Baicheng 137000， Jilin， China College of Physics， Baicheng Normal University， Baicheng 137000， Jilin， China 采用密度泛函理论的M06-2X和MN15方法，在6-31++G（d，p）和6-311++G（2df，pd）基组水平下，研究了气相 S -型两性脯氨酸（ S -Pro_1）和中性脯氨酸（ S -Pro_2）与Ca 2+ 配合物的旋光异构及水分子（簇）的催化作用。反应通道研究发现： S -Pro_1·Ca 2+ 的旋光异构有2个通道a和b，a通道是 α -氢以氧为桥迁移，b通道是 α -氢从 α -碳迁移至氧后，再从氮向 α -碳迁移； S -Pro_2·Ca 2+ 的旋光异构也有2个通道a和b，a通道是螯合环打开后向 S -Pro_1·Ca 2+ 异构，b通道是羧基内质子迁移后， α -氢再以羰基氧为桥迁移，向 R -Pro_2·Ca 2+ 异构。势能面计算表明： S -Pro_1·Ca 2+ 旋光异构的a通道最具优势，反应活化能为279.4 kJ/mol，水分子（簇）的催化作用使反应活化能降至149.1 kJ/mol； S -Pro_2·Ca 2+ 旋光异构的a通道也最具优势，反应活化能为217.5 kJ/mol。结果表明，气相脯氨酸钙配合物很难发生旋光异构，水分子存在时脯氨酸钙能痕量地消旋。 We have studied the optical isomerism of S -amphoteric proline ( S -Pro_1) Ca 2+ complexes and neutral proline ( S -Pro_2) Ca 2+ complexes in gas phase and the catalytic action of water molecules(clusters) using the methods of M06-2X and MN15 based on the density functional theory at 6-31++G(d,p) and 6-311++G(2df,pd) basis set level. It is identified that the optical isomerism of S -Pro_1·Ca 2+ occurs in two pathways a and b. α -H is transferred using the O atom in pathway a. In pathway b, α -H is transferred to O atom from α -C and then it is transferred to α -C from the N atom. The optical isomerism of S -Pro_2·Ca 2+ also occurs in two pathways a and b. In pathway a, the chelate ring is broken and S -Pro_2·Ca 2+ isomerizes to S -Pro_1·Ca 2+ . In pathway b, α -H is transferred using carbonyl O atom as a bridge after the transition of proton on the carboxyl group, and then it isomerizes to R -Pro_2·Ca 2+ . The calculations of potential energy surface show that the optical isomerism of S -Pro_1·Ca 2+ is the most advantageous one in the pathway a, and the activation energy of the reaction is 279.4 kJ/mol. The activation energy of the reaction is reduced to 149.1 kJ/mol under the catalytic action of water molecules (clusters); the optical isomerism of S -Pro_2·Ca 2+ is the most advantageous in the pathway a, and the activation energy of the reaction is 217.5 kJ/mol. The results show that the optical isomerism of S -Pro·Ca 2+ complexes are difficult to occur in gas phase, and S -Pro·Ca 2+ can only racemize at an extreme trace when water molecules are there. 脯氨酸钙 水分子 密度泛函理论 旋光异构 amino acid metal complex water molecules density functional theory optical isomerism 0　 引 言 1 氨基酸金属配合物（amino acid metal complex）是蛋白质、酶等生物大分子的活性中心，参与生命体的代谢及生化过程，在生命过程中扮演着重要的角色，是生命科学、医学和农业等领域的研究热点 ［ 1 1 ］ 。Marino等 ［ 2 2 ］ 和Lamsabhi等 ［ 3 3 ］ 研究了Cu 2+ 、Ca 2+ 分别与丙氨酸和丝氨酸的配合物的结构稳定性；Nidhi等 ［ 4 4 ］ 研究了丙氨酸与Ca 2+ 、Cu 2+ 和Zn 2+ 的配合物由于质子转移诱导而发生的构型变化（非手性变化）。手性氨基酸是生物体内的重要配体，其手性特征使其金属离子配合物也具有手性。手性药物的药理作用是药物分子通过与体内的大分子之间严格的手性识别和匹配实现的，通常手性分子只有一种对映体具有所要求的药理活性，而另一种对映体没有药理作用甚至有严重的副作用 ［ 5 5 ， 6 6 ］ 。“反应停”惨剧 ［ 7 7 ］ 给予了科研工作者一个重要的警示：手性药物使用前必须明确其旋光异构机制，优构体必须能较好地保持其手性特征。因此，研究手性氨基酸金属配合物的旋光异构，对于指导人们科学地选择使用氨基酸金属盐进行氨基酸与金属离子同补具有重要意义。基于此，王佐成课题组已采用密度泛函理论，研究了孤立条件下 α -Ala（丙氨酸，结构最简单的手性氨基酸）与生命体必需的Cu 2+［8］ 、Na +［9］ 、Zn 2+［10］ 、Mn 2+［11］ 、K +［12］ 、Fe 2+［13］ 、Mg 2+ ［ 14 14 ］ 和Ca 2+［15］ 等金属离子配合物的旋光异构。结果表明，这些配合物消旋能垒分别为123.9、253.3、197.2、247.1、254.7、233.5、283.9和297.5 kJ/mol，除 α -Ala·Cu 2+ 可以缓慢消旋，其他均不能消旋。文献［ 14 14 ， 15 15 ］的密度泛函理论研究还表明，水分子（簇）对 α -Ala与Ca 2+ 和Mg 2+ 配合物的手性转变有很好的催化作用，消旋能垒分别降至159.3和157.9 kJ/mol，表明水汽环境下 α -Ala与Ca 2+ 、Mg 2+ 的配合物均不能消旋，故两种配合物均易于保存。 脯氨酸（proline，Pro）是构成蛋白质的主要成分，按其构型分为 S -Pro和 R -Pro，按旋光作用分为 L -Pro（左旋体）和 D -Pro（右旋体）。在生命体内， L -Pro具有生物活性，可合成胶原蛋白，促进软骨修复，具有治疗骨关节炎、减轻持续性软组织疼痛和慢性背痛的作用 ［ 16 16 ］ 。以往研究 ［ 17 17 ， 18 18 ］ 表明，单体脯氨酸在生命体内可以缓慢消旋，单独用于生命体会产生一定的副作用。钙是生命体必需的宏量元素，在生命体内以2价态存在，对维持骨骼的构成、心脏的规律跳动以及协调酸碱平衡等均有重要作用 ［ 19 19 ~ 21 21 19~21 ］ 。食物中的钙被人体吸收过程为：食物分解出的钙离子在小肠处与氨基酸及短肽链螯合，形成手性氨基酸钙，再以配合物的形式被整体吸收 ［ 22 22 ］ 。氨基酸螯合钙作为补钙产品更容易被吸收，在临床、食品及畜牧业中已经被应用 ［ 23 23 ］ 。因此，考察脯氨酸螯合钙用于生命体同补脯氨酸和钙是否安全具有重要意义。 对孤立体系的研究，可获得体系本质的物理和化学性质，可以为深入研究该体系在复杂环境下的行为奠定基础 ［ 24 24 ， 25 25 ］ 。基于此，考虑到大气中的水分子 ［ 26 26 ］ 多以水分子（簇）形式存在 ［ 27 27 ］ ，并与其他分子形成氢键网络 ［ 28 28 ］ ，结合文献［ 8 8 ~ 15 15 8~15 ］的经验，研究了 S -Pro·Ca 2+ 配合物在气相及水分子（簇）作用下的手性旋光异构。本文详细讨论了水分子（簇）的作用下 S -Pro·Ca 2+ 在优势通道上完整的旋光异构过程及其他反应通道重要的质子迁移过程，以探讨水分子（簇）对质子迁移反应的催化作用。 1　 研究与计算方法 1 采用基于第一原理的色散校正泛函M06-2X方法 ［ 29 29 ］ ，结合基组6-31++G（d， p）计算单重态（该体系的最稳定态）反应势能面上的驻点结构（自由能热校正温度为298.15 K）；采用NBO（natural bond orbital）自然键轨道理论计算NPA（natural population analysis）电荷；对搜索的过渡态 ［ 30 30 ］ 计算了内禀反应坐标（IRC） ［ 31 31 ］ ，证明过渡态可以信赖。采用具有更高精度的泛函方法MN15 ［ 32 32 ］ ，结合更高级别基组6-311++G（2df， pd）计算驻点的单点能。总自由能为自由能热校正和单点能之和，势能面上势能曲线的零势点为该曲线上驻点相对的势能零点。 S -Pro构型1与Ca 2+ 的配合物记作 S -P_1； S -P_1在反应通道b的第4个 R -型过渡态用 R -T4_1 b 表示，第4个 R -型中间体用 R -I4_1 b 表示；在反应通道a和b共用的第1个没有手性的过渡态用T1_1 a（b） 表示， S -P_2在a1通道异构得到的 S -P_1用 S -P_1_2 a1 表示；6个H 2 O与 S -P_1中的Ca 2+ 配位，同时2个H 2 O与 S -P_1氢键作用的配合物络合物用 S -P_1←6H 2 O·（H 2 O） 2 表示，其他体系表示法相似。文中计算由Gaussian16程序 ［ 33 33 ］ 完成。 2　 结果与讨论 1 Pro与Ca 2+ 配位有两性和中性两种配合物。计算表明两性Pro与Ca 2+ 配位相对稳定，优化的两对手性对映体见 图1 图1 。下文对气相 S -P_1和 S -P_2的旋光异构以及水分子（簇）的催化作用分别进行讨论。 2.1　 S -P_1的旋光异构 2 S -P_1旋光异构有2个通道a和b，反应历程和驻点构象见 图2 图2 ，反应历程的势能面见 图3 图3 。 第1基元反应：反应物 S -P_1经过渡态T1_1 a（b） ， α -氢13H从4C迁移至16O，异构成中间体I1_1 a（b） 。从 S -P_1到T1_1 a（b） ，13H—4C从0.109 1 nm拉伸至0.151 6 nm，断裂，二面角5N—4C—3C—14C从118.5°增至133.3°，C—H键的断裂和骨架的形变使T1_1 a（b） 产生了279.4 kJ/mol的内禀能垒。 第2基元反应：I1_1 a（b） 经过渡态T2_1 a（b） ，五元环形变异构，得到第2中间体I2_1 a（b） 。从I1_1 a（b） 到T2_1 a（b） ，二面角2C—1C—5N—4C从-33.3°变为-21.9°，二面角3C—4C—14C—15O从-178.3°变为-170.5°，二面角形变需要能量很少，T2_1 a（b） 产生的内禀能垒只有13.1 kJ/mol。接下来分别介绍a和b通道专属的反应历程。 2.1.1　a通道专属反应 3 第3基元反应：I2_1 a（b） 经和T1_1 a（b） 镜像对称的过渡态T3_1 a ，异构成两性的 R -型产物 R -P_1 a ， S -P_1实现旋光异构。结构分析表明， R -P_1 a 是 S -P_1的手性对映体。从I2_1 a（b） 到T3_1 a ，13H—16O从0.096 8 nm拉伸到0.124 0 nm，T3_1 a 产生了181.6 kJ/mol的内禀能垒。该能垒小于T2_1 a（b） 产生的内禀能垒，原因是此基元O—H键拉伸幅度小于第2基元C—H键拉伸幅度，且从I2_1 a（b） 到T3_1 a ， α -碳从sp 2 杂化向sp 3 杂化过渡，体系释放能量。 从a通道上 S -P_1的驻点构象和势能面可以看出，T2_1 a（b） 是对称点。 2.1.2　b通道专属反应 3 第3基元反应：I2_1 a（b） 经过渡态T3_1 b ，17H在纸面内侧从N向 α -碳迁移，得到 R -I3_1 b ， S -P_1在b通道实现旋光异构。从I2_1 a（b） 到T3_1 b ，5N—17H键从0.102 5 nm拉伸至0.139 3 nm ，断裂，T3_1 b 产生了224.0 kJ/mol的内禀能垒，与孤立的Pro构型反转在此基元的内禀能垒111.2 kJ/mol ［ 17 17 ］ 相比升高显著，Ca 2+ 表现了极强的负催化作用。这是由于Ca 2+ 与O配位使I2_1 a（b） 中17H的电荷从0.398 e（1 e=1.602×10 -19 C）大幅增加至0.435 e，5N的电荷从0.763 e增加至0.785 e，4C的电荷量从-0.256 e急剧下降至-0.093 e，从而N对17H的库仑引力大幅增加，4C对17H的库仑引力大幅减小。 第4基元反应： R -I3 b _1经14C—4C右视顺时针内旋转的过渡态 R -T4_1 b ，异构到两性产物 R -P_1 b 。从 R -I3 b _1到 R -T4_1 b ，二面角15O—14C—4C—5N从20.2°增加至67.8°，化学键的内旋转需要的能量很小， R -T4_1 b 产生的内禀能垒仅10.4 kJ/mol。结构分析表明， R -P_1 b 是 S -P_1的手性对映体。 从 图3 图3 可以看出， S -P_1在a通道的旋光异构具有优势，在b通道的异构处于劣势。优势反应通道的总包能垒（气相反应的活化能为总包能垒）为279.4 kJ/mol，来自过渡态T1_1 a（b） ，该能垒远高于极限反应能垒167.0 kJ/mol ［ 34 34 ］ ，说明孤立的 S -P_1不能消旋；劣势通道的总包能垒为321.8 kJ/mol，来自过渡态T3_1 b 。 2.2　 S -P_2的旋光异构 2 S -P_2的旋光异构在a和b两个通道实现，反应历程见 图4 图4 ，反应历程的自由能势能面见 图5 图5 。 2.2.1　a通道 3 第1基元反应： S -P_2经16O—14C内旋转的过渡态 S -T1m_2 a 或 S -T1n_2 a ，羧羟基氢17H从羧基内侧转到外侧，异构成中间体 S -I1_2 a 。从 S -P_2到 S -T1m_2 a ，二面角17H—16O—14C—15O从1.1°增加至98.1°， S -T1m_2 a 产生了53.7 kJ/mol的内禀能垒，相似地， S -T1n_2 a 产生了55.8 kJ/mol的内禀能垒。 第2基元反应： S -I1_2 a 经18Ca—15O—14C做翻转运动的过渡态 S -T2_2 a ，异构成 S -I2_2 a 。从 S -I1_2 a 到 S -T2_2 a ，Ca—N配位键断裂，Ca离子从羧基外转到羧基内， S -T2_2 a 产生了107.9 kJ/mol的内禀能垒。接着 S -I2_2 a 在两个分路径a1和a2上异构。 2.2.1.1 a1路径专属反应 4 第3基元反应： S -I2_2 a 经过渡态 S -T3_2 a1 ，14C—4C右视逆时针内旋转，异构成 S -P_1_2 a1 。此基元与2.1节b通道的第4基元相似， S -T3_2 a1 产生了10.1 kJ/mol的内禀能垒，不再赘述。结构分析表明， S -P_1_2 a1 与 S -P_1完全相同， S -P_1_2 a1 的异构与 S -P_1的异构完全相同，见2.1节。 2.2.1.2 a2路径专属的反应 4 第3基元反应： S -I2_2 a 经氨基在纸面内侧外翻转的过渡态 S -T3_2 a2 ，异构成 S -I3_2 a2 。从 S -I2_2 a 到 S -T3_2 a2 ，二面角6H—5N—1C—4C从-130.6°变为164.6°，二面角异构所需能量很少， S -T3_2 a2 产生的内禀能垒为7.0 kJ/mol。 第4基元反应： S -I3_2 a2 经13H从4C向5N迁移的过渡态T4_2 a2 ，异构成I4_2 a2 。从 S -I3_2 a2 到T4_2 a2 ，13H—4C从0.110 2 nm拉伸至0.122 5 nm，断裂，二面角5N—4C—3C—14C从123.9°增加至165.8°，T4_2 a2 产生了123.8 kJ/mol的内禀能垒。与孤立的Pro_2此基元反应的内禀能垒246.0 kJ/mol ［ 17 17 ］ 相比有较大幅度下降，表明Ca 2+ 的配位产生了极好的催化作用。原因是Ca 2+ 与O配位导致13H—4C从0.100 1 nm拉伸至0.110 2 nm， α -氢的红外伸缩振动频率从3 029.1 cm -1 降至2 890.3 cm -1 ，C—H键被活化。I4_2 a2 与I2_1 a（b） 构象相同，其与I2_1 a（b） 的异构相同。 2.2.2　b通道 3 第1基元反应： S -P_2经过渡态 S -T1_2 b ，17H在羧基内迁移，异构成 S -I1_2 b 。从 S -P_2到 S -T1_2 b ，17H—16O从0.097 5 nm被拉伸了0.037 5 nm， S -T1_2 b 产生了183.9 kJ/mol的内禀能垒。 第2基元反应： S -I1_2 b 经质子在4C和16O间迁移的过渡态T2_2 b ， α -氢从4C向16O迁移，异构成I2_2 b 。从 S -I1_2 b 到T2_2 b ，4C—113H键从0.109 4 nm拉伸至0.152 1 nm，断裂，二面角5N—4C—3C—14C从120.4°变为137.3°，化学键断裂及骨架形变需较大能量，产生了287.2 kJ/mol的内禀能垒。 第3基元反应：I2_2 b 经16O—14H内旋转的过渡态T3_2 b ，13H绕16O—14H仰视逆时针旋转到纸面内侧，异构成I3_2 b 。从I2_2 b 到T3_2 b ，化学键内旋转所需能量不多，T3_2 b 产生的内禀能垒为71.8 kJ/mol。 第4基元反应：I3_2 b 经过渡态T4_2 b ，13H在纸面内侧从16O迁移至4C，异构成 R -I4_2 b ， S -P_2实现旋光异构。从I3_2 b 到T4_2 b ，16O—13H键从0.096 9 nm拉伸至0.135 5 nm，断裂，T4_2 b 产生的内禀能垒为301.1 kJ/mol。此基元反应的内禀能垒较高，原因之一是I3_2 b 分子中的13H距离4C较远，原因之二是16O—13H键很稳定（键长较短，红外伸缩振动频率很大，是3 796.6 cm -1 ）。 第5基元反应： R -I4_2 b 经过渡态 R -T5_2 b ，实现了17H从15O向16O的迁移，异构成产物 R -P_2 b 。从 R -I4_2 b 到 R -T5_2 b ，17H—15O从0.097 6 nm拉伸至0.126 4 nm，导致 R -T5_2 b 产生了107.7 kJ/mol的内禀能垒。此基元反应比第1基元反应的内禀能垒小很多，原因之一是 R -I4_2 b 分子中的15O—17H比 S -P_2的16O—17H键的活化程度高（前者比后者键长），原因之二是此基元O—H键的拉伸幅度较小，原因之三是此基元的质子迁移是逆偶极矩方向，体系的电场力会助力该过程的质子迁移。 结合 图5 图5 和 图3 图3 可以看出， S -P_2在a1通道向 S -P_1异构后，接 S -P_1的a通道的旋光异构具有优势，总包能垒为279.4-61.9=217.5 kJ/mol（来自过渡态T1_1 a（b） ）。该能垒高于化学反应的极限能垒167.0 kJ/mol ［ 34 34 ］ 很多，说明通常情况下气相 S -P_2的消旋能垒难以逾越。从 图5 图5 还可以看出， S -P_1比 S -P_2稳定， S -P_2向 S -P_1异构的总包能垒为131.2 kJ/mol，来自羧基钙基团旋转的过渡态 S -T3_2 a1 。该能垒高于温和反应能垒80.0 kJ/mol ［ 34 34 ］ 很多，距离极限反应能垒167.0 kJ/mol ［ 34 34 ］ 不大，故 S -P_2向 S -P_1的异构也十分缓慢。从 图5 图5 还可以看出， S -P_1向 S -P_2异构的总包能垒为131.2+61.9=193.1 kJ/mol，该能垒高于极限反应能垒167.0 kJ/mol ［ 34 34 ］ ，表明通常情况下 S -P_1向 S -P_2的异构难以进行。结合 图5 图5 和 图3 图3 还可以看出， S -P_2在a2通道向 S -P_1的异构为亚优势通道， S -P_2在a2通道向 S -P_1异构的活化能为250.5 kJ/mol，来自过渡态T4_2 a2 。 2.3　水分子（簇）对旋光异构的影响 2 S -P_2在a1通道向 S -P_1异构的所有基元反应均为非质子迁移反应，而对于非质子迁移反应，水分子簇的作用很小，可以不予考虑。篇幅所限，本文仅对 S -P_1在a通道旋光异构的整个过程、 S -P_1在b通道旋光异构的决速步骤、 S -P_2在a2通道向 S -P_1异构的决速步骤进行讨论。 2.3.1　水分子（簇）作用下 S -P_1在a通道的旋光异构 3 水分子（簇）与底物配合物作用，配位键强于氢键，所以水分子优先与Ca“配位饱和”，再与底物氢键作用。Ca 2+ 与O最多是八配位，底物的Ca 2+ 已与羰基O二配位，故再与钙配位的水分子数取最大值6。计算表明，质子在 α -碳和O或N之间迁移时，二聚水和三聚水的催化作用基本相同，均远大于单个水分子做质子转移媒介的情形，因此只讨论2个水分子簇作为H迁移媒介的情况。对于第2基元杂环异构，水分子（簇）的氢键作用对能垒影响很小，只讨论水分子与Ca 2+ 配位的情况。水分子（簇）作用下 S -P_1在a通道旋光异构历程及过渡态见 图6 图6 ，反应的自由能势能面见 图7 图7 。 第1基元反应： S -P_1与6个水分子配位，二聚水再与其 α -碳和羰基O氢键作用并与和Ca 2+ 配位的水分子形成稳定的氢键网络，得到较稳定的反应物络合物 S -P_1 a（b） ←6H 2 O·（H 2 O） 2 m（m表示与底物氢键作用的二聚水在纸面外侧，相似之处不再说明）。 S -P_1 a（b） ←6H 2 O·（H 2 O） 2 m经过渡态T1_1 a（b） ←6H 2 O·（H 2 O） 2 m，实现了氢质子从 α -碳向羰基O的迁移，异构成I1_1 a（b） ←6H 2 O·（H 2 O） 2 m。 从 S -P_1 a（b） ←6H 2 O·（H 2 O） 2 m到T1_1 a（b） ←6H 2 O·（H 2 O） 2 m过程，4C—13H键从0.109 0 nm拉伸至0.133 3 nm，断裂、36O—38H和39O—41H分别从0.097 5和0.097 3 nm拉伸断裂；骨架二面角5N—4C—3C—14C从116.5°增加至121.6°。3个化学键的拉伸断裂及骨架的形变需要一定的能量，T1_1 a（b） ←6H 2 O·（H 2 O） 2 m产生的内禀能垒为149.1 kJ/mol，远低于T1_1 a（b） 产生的内禀能垒279.4 kJ/mol，水分子（簇）对此基元反应起到了极好的催化作用。其原因有三：一是T1_1 a（b） ←6H 2 O·（H 2 O） 2 m的氢键角4C—13H—36O、36O—38H—39O和39O—41H—46O分别为164.7°、158.2°和156.2°，接近平角，3个氢键均较强；二是水分子（簇）作用下 S -P_1 a（b） ←6H 2 O·（H 2 O） 2 m到T1_1 a（b） ←6H 2 O·（H 2 O） 2 m过程骨架5N—4C—3C—14C形变显著减小；三是过渡态的八元环结构共面程度较好，过渡态的水分子形成了比较稳定的氢键网络。T1_1 a（b） ←6H 2 O·（H 2 O） 2 m靠近产物中间体，是晚期过渡态。 第2基元反应：I1 a（b） ←6H 2 O经过渡态T2_1 a（b） ←6H 2 O，五元杂环异构，二面角5N—1C—2C—3C从40.3°变为-40.3°，构象异构成I2_1 a（b） ←6H 2 O。从I1 a（b） ←6H 2 O到T2_1 a（b） ←6H 2 O过程，5N—4C—3C—14C从40.3°增加至-0.5°，二面角异构所需能量不多，T2_1 a（b） ←6H 2 O产生的内禀能垒为12.6 kJ/mol，与T2_1 a（b） 产生的内禀能垒13.1 kJ/mol相差无几，可以认为没改变。这说明水分子与Ca 2+ 配位对 R -基杂环的形变能垒影响微乎其微。 第3基元反应：I2_1 a（b） ←6H 2 O的4C、13H、20H和32H与二聚水氢键作用，形成络合物I2_1 a（b） ←6H 2 O·（H 2 O） 2 n（n表示（H 2 O） 2 在纸面内侧与底物氢键作用）经3质子协同迁移的过渡态T3_1 a ←6H 2 O·（H 2 O） 2 n，实现了质子在纸面内侧从O向 α -碳净迁移，异构成产物复合物 R -P_1 a ←6H 2 O·（H 2 O） 2 n。至此， S -P_1 a（b） ←6H 2 O·（H 2 O） 2 m实现了手性对映体转变。 从I2_1 a（b） ←6H 2 O·（H 2 O） 2 n到反应的过渡态T3_1 a（b） ←6H 2 O·（H 2 O） 2 n，16O—13H、41O—40H和37O—39H键长分别从0.097 6、0.097 2和0.097 9 nm拉伸至0.100 8、0.100 9 和0.128 5 nm。两个很小的化学键拉伸及1个化学键小幅拉伸断裂，T3_1 a ←6H 2 O·（H 2 O） 2 n产生的能垒仅35.8 kJ/mol，远低于T3_1 a 产生的能垒181.6 kJ/mol，二聚水起了极好的催化作用。原因有三：一是T3_1 a（b） ←6H 2 O·（H 2 O） 2 n的水分子间存在稳定的氢键网络；二是八元环存在3个较强的氢键；三是从I2_1 a（b） ←6H 2 O·（H 2 O） 2 n到T3_1 a（b） ←6H 2 O·（H 2 O） 2 n过程，4C从sp 2 杂化向sp 3 杂化过渡（从空配向满配），体系释放能量。T3_1 a（b） ←6H 2 O·（H 2 O） 2 n靠近中间体反应物，是早期过渡态。 从 图7 图7 可以看出，水分子（簇）作用下 S -P_1在a通道旋光异构的自由能垒为149.1 kJ/mol，距反应的极限能垒167.0 kJ/mol ［ 34 34 ］ 较近，说明水分子簇作用下 S -P_1的旋光异构只能极其缓慢地进行。 2.3.2　水分子（簇）作用下 S -P_1在b通道旋光异构的决速步骤 3 S -P_1在b通道旋光异构的决速步骤是I2_1 b 经过渡态T3_1 b 异构成 R -I3_1 b 。水分子（簇）存在时，I2_1 b 与6个水分子配位，在其 α -碳和氨基氢17H与纸面内侧二聚水氢键作用同时与和钙配位的水分子形成稳定的氢键网络，得到稳定的中间体反应物复合物I2_1 b ←6H 2 O·（H 2 O） 2 n。I2_1 b ←6H 2 O·（H 2 O） 2 n经过渡态T3_1 b ←6H 2 O·（H 2 O） 2 n，实现了质子从氨基N向 α -碳的净迁移，异构成中间体产物复合物 R -I3_1 b ←6H 2 O·（H 2 O） 2 n。 该基元的反应历程见 图8 图8 ，自由能势能面见 图9 图9 。 从I2_1 b ←6H 2 O·（H 2 O） 2 n到T3_1 b ←6H 2 O·（H 2 O） 2 n，5N—17H断裂，37O—39H和40O—41H键长从0.097 5 nm分别拉伸至0.136 4和0.112 7 nm，需要一些能量。但由于4C从sp 2 杂化向sp 3 杂化过渡，体系释放能量，并且T3_1 b ←6H 2 O·（H 2 O） 2 n的七元环存在较强的氢键，构象相对稳定。因此，过渡态T3_1 b ←6H 2 O·（H 2 O） 2 n产生的能垒不是很高，是107.3 kJ/mol。I2_1 b 在势能面的位置是97.8 kJ/mol，因此在水分子（簇）作用下， S -P_1在b通道旋光异构的活化能（总包能垒）变为107.3 + 97.8=205.1 kJ/mol。该能垒高于水分子（簇）作用下 S -P_1在a通道旋光异构的自由能垒149.1 kJ/mol，因此，水分子（簇）作用下 S -P_1旋光异构的优势通道仍然是a通道。T3_1 b ←6H 2 O·（H 2 O） 2 n靠近中间体反应物，是早期过渡态。 2.3.3　水分子（簇）作用下 S -P_2在a2通道向 S -P_1异构的决速步骤 3 S -P_2在a2通道向 S -P_1异构的决速步骤是 S -I3_2 a2 经T4_2 a2 向I4_2 a2 异构。水分子（簇）存在时， S -I3_2 a2 与6个水分子配位，其 α -氢和氨基N与纸面外侧2个水分子簇氢键作用并与和钙配位的水分子形成稳定的氢键网络，得到中间体反应物复合物 S -I3_2 a2 ←6H 2 O·（H 2 O） 2 m。 S -I3_2 a2 ←6H 2 O·（H 2 O） 2 m经3个质子协同迁移的过渡态T4_2 a2 ←6H 2 O·（H 2 O） 2 m，质子从 α -碳向N净迁移，异构成产物中间体复合物I4_2 a2 ←6H 2 O·（H 2 O） 2 m。该基元的反应历程见 图10 图10 ，自由能势能面见 图11 图11 。 从 S -I3_2 a2 ←6H 2 O·（H 2 O） 2 m到T4_2 a2 ←6H 2 O·（H 2 O） 2 m过程，4C—13H、40O—41H和37O—38H键长从0.109 4、0.097 7和0.100 1 nm分别拉伸至0.153 1、0.107 5和0.100 3 nm，需要一些能量；骨架二面角5N—4C—3C—14C从121.4°增加至147.5°；水分子间的氢键网络被破坏。这些变化使T4_2 a2 ←6H 2 O·（H 2 O） 2 m产生了151.5 kJ/mol的内禀能垒，比T4_2 a2 产生的内禀能垒123.8 kJ/mol高22.3%，水分子（簇）表现出较显著的负催化作用。这有三个方面的原因：一是水分子（簇）的作用使4C—13H键的键长变小（从0.110 2 nm变为0.109 4 nm），共价键作用增强，变得不易断裂；二是使4C—13H的拉升幅度显著增加（从0.012 3 nm增加至0.043 7 nm）；三是从 S -I3_2 a2 ←6H 2 O·（H 2 O） 2 m到T4_2 a2 ←6H 2 O·（H 2 O） 2 m过程破坏了水分子间的氢键网络。因为 S -I3_2 a2 在势能面的位置是126.7 kJ/mol，因此在水分子（簇）作用下， S -P_2在a2通道向 S -P_1异构的活化能（总包能垒）变为126.7+151.5=278.2 kJ/mol， S -P_2在a2通道向 S -P_1异构的劣势进一步增加。T4_2 a2 ←6H 2 O·（H 2 O） 2 m靠近中间体反应物，是早期过渡态。 3　 结 论 1 采用密度泛函理论的M06-2X和MN15方法，在6-311+G（d，p）和6-311++G（2df，pd）基组水平，研究了两性和中性的Pro与Ca 2+ 配合物的旋光异构及水分子（簇）的催化作用，得到如下结论： 1） S -Pro_1·Ca 2+ 旋光异构优势通道的总包能垒为279.4 kJ/mol（来自质子从 α -碳向羰基O迁移的过渡态）； S -Pro_2·Ca 2+ 先越过131.2 kJ/mol的能垒（来自羧基钙基团旋转的过渡态）向 S -Pro_1·Ca 2+ 异构，然后再实现旋光异构的通道具有优势； 2） 水分子（簇）的催化使 S -Pro_1·Ca 2+ 旋光异构反应的总包能垒降至149.1 kJ/mol；3） 水分子（簇）对质子迁移反应有正催化作用或负催化作用。结果表明， S -Pro_2·Ca 2+ 向 S -Pro_1·Ca 2+ 的异构很缓慢，气相脯氨酸钙的手性难以改变，水汽相下脯氨酸钙只能痕量地消旋，因此脯氨酸钙在气固相储存具有安全性。Pro分子中的 R -基为五元杂环，其金属配合物手性对映体异构的研究远比 α -Ala的对映异构研究复杂，因此本研究可望为水液相下 S -Pro·Ca 2+ 及更复杂的氨基酸金属配合物对映异构研究提供借鉴。 参考文献 董璐 . 氨基酸金属配合物的合成与表征 ［D］. 北京 ： 北京理工大学 ， 2015 . DOI: 10.21698/simi2015 http://dx.doi.org/10.21698/simi2015 . 氨基酸金属配合物的合成与表征 DONG L . The Synthesis and Characterization of Coordination Complexes with Amino Acid Ligands ［D］. Beijing ： Beijing Institute of Technology ， 2015 （Ch） . DOI: 10.21698/simi2015 http://dx.doi.org/10.21698/simi2015 . The Synthesis and Characterization of Coordination Complexes with Amino Acid Ligands MARINO T ， RUSSO N ， TOSCANO M . Interaction of Cu + and Cu 2+ ions with α -alanine. A density functional study ［J］. Journal of Mass Spectrometry ， 2002 ， 37 （ 8 ）： 786 - 791 . DOI： 10.1002/jms.331 http://dx.doi.org/10.1002/jms.331 . Interaction of Cu+ and Cu2+ ions with α-alanine. A density functional study LAMSABHI A M ， MÓ O ， YÁÑEZ M . Serine-Ca 2+ versus serine-Cu 2+ complexes—A theoretical perspective ［J］. Canadian Journal of Chemistry ， 2010 ， 88 ： 759 - 768 . DOI： 10.1139/V10-038 http://dx.doi.org/10.1139/V10-038 . Serine-Ca2+ versus serine-Cu2+ complexes—A theoretical perspective VYAS N ， OJHA A K . Investigation on transition states of ［alanine + M 2+ ］ （M = Ca， Cu， and Zn） complexes： A quantum chemical study ［J］. International Journal of Quantum Chemistry ， 2012 ， 112 （ 5 ）： 1526 - 1536 . DOI： 10.1002/qua.23141 http://dx.doi.org/10.1002/qua.23141 . Investigation on transition states of ［alanine + M2+］ （M = Ca， Cu， and Zn） complexes： A quantum chemical study CHOI J R ， YONG K W ， CHOI J Y ， et al . Progress in molecularly imprinted polymers for biomedical applications ［J］. Combinatorial Chemistry & High Throughput Screening ， 22 （ 2 ）： 78 - 88 . DOI： 10.2174/ 1386207322666190325115526 http://dx.doi.org/10.2174/1386207322666190325115526 . Progress in molecularly imprinted polymers for biomedical applications WARD T J ， WARD K D . Chiral separations： Fundamental review 2010 ［J］. Analytical Chemistry ， 2010 ， 82 （ 12 ）： 4712 - 4722 . DOI： 10.1021/ac1010926 http://dx.doi.org/10.1021/ac1010926 . Chiral separations： Fundamental review 2010 齐晓霞 ． 药害事故防范与救济制度研究 ［D］. 上海 ： 复旦大学 ， 2011 . DOI: 10.15385/yb.miracle.2011 http://dx.doi.org/10.15385/yb.miracle.2011 . 药害事故防范与救济制度研究 QI X X . The Study on the Drug - Induced Accident and Insurance Relief ［D］. Shanghai ： Fudan University ， 2011 （Ch） . DOI: 10.15385/yb.miracle.2011 http://dx.doi.org/10.15385/yb.miracle.2011 . The Study on the Drug-Induced Accident and Insurance Relief 李冰 ， 徐锐英 ， 姜春旭 ， 等 . 气相环境下Cu 2+ 催化 α -丙氨酸手性对映体转变的机理 ［J］. 武汉大学学报（理学版） ， 2019 ， 65 （ 6 ）： 621 - 628 . DOI： 10.14188/j.1671-8836.2019.06.013 http://dx.doi.org/10.14188/j.1671-8836.2019.06.013 . 气相环境下Cu2+催化α-丙氨酸手性对映体转变的机理 LI B ， XU R Y ， JIANG C X ， et al . Mechanism of chiral enantiomer transition of α -Ala catalyzed by Cu 2+ in gas phase environment ［J］. Journal of Wuhan University（Science Edition） ， 2019 ， 65 （ 6 ）： 621 - 628 . DOI： 10.14188/j.1671-8836.2019.06.013（Ch http://dx.doi.org/10.14188/j.1671-8836.2019.06.013（Ch ）. Mechanism of chiral enantiomer transition of α-Ala catalyzed by Cu2+ in gas phase environment 李冰 ， 陈静思 ， 徐锐英 ， 等 . 气相环境下Na + 催化丙氨酸分子手性转变的机理 ［J］. 复旦学报（自然科学版） ， 2020 ， 59 （ 4 ）： 495 - 504 . DOI： 10.15943/j.cnki.fdxb-jns.2020.04.015 http://dx.doi.org/10.15943/j.cnki.fdxb-jns.2020.04.015 . 气相环境下Na+催化丙氨酸分子手性转变的机理 LI B ， CHEN J S ， XU R Y ， et al . Mechanism of chiral transition of alanine catalyzed by Na + in gas phase ［J］. Journal of Fudan University （Natural Sciences） ， 2020 ， 59 （ 4 ） ： 495 - 504 . DOI： 10.15943/j.cnki.fdxb-jns.2020.04.015（Ch http://dx.doi.org/10.15943/j.cnki.fdxb-jns.2020.04.015（Ch ）. Mechanism of chiral transition of alanine catalyzed by Na+ in gas phase 车立新 ， 张雪娇 ， 刘芳 ， 等 . 气相 α -丙氨酸二价锌配合物手性转变的密度泛函理论研究 ［J］. 云南大学学报（自然科学版） ， 2020 ， 42 （ 5 ）： 966 - 976 . DOI： 10.7540/j.ynu. 20200038 http://dx.doi.org/10.7540/j.ynu.20200038 . 气相α-丙氨酸二价锌配合物手性转变的密度泛函理论研究 CHE L X ， ZHANG X J ， LIU F ， et al . Chiral transition mechanism of Zn（Ⅱ） and α -Ala complexes in gas phase： A density functional theory study ［J］. Journal of Yunnan University （Natural Sciences Edition） ， 2020 ， 42 （ 5 ）： 966 - 976 . DOI： 10.7540/j.ynu.20200038（Ch http://dx.doi.org/10.7540/j.ynu.20200038（Ch ）. Chiral transition mechanism of Zn（Ⅱ） and α-Ala complexes in gas phase： A density functional theory study 张雪娇 ， 李冰 ， 姜春旭 ， 等 . 气相环境下丙氨酸Mn（Ⅱ）配合物的手性转变机理 ［J］. 武汉大学学报（理学版） ， 2020 ， 66 （ 4 ）： 345 - 355 . DOI： 10.14188/j.1671-8836.2020.0024 http://dx.doi.org/10.14188/j.1671-8836.2020.0024 . 气相环境下丙氨酸Mn（Ⅱ）配合物的手性转变机理 ZHANG X J ， LI B ， JIANG C X ， et al . Chiral transition mechanism of Mn（Ⅱ） and Ala complexes in gas phase environment ［J］. Journal of Wuhan University（Science Edition） ， 2020 ， 66 （ 4 ）： 345 - 355 . DOI： 10.14188/j.1671-8836.2020.0024（Ch http://dx.doi.org/10.14188/j.1671-8836.2020.0024（Ch ）. Chiral transition mechanism of Mn（Ⅱ） and Ala complexes in gas phase environment 刘芳 ， 马宏源 ， 潘宇 ， 等 . 气相环境下K + 催化丙氨酸分子手性对映体转变的机理 ［J］. 厦门大学学报（自然科学版） ， 2020 ， 59 （ 6 ）： 912 - 920 . DOI： 10.6043/j.issn.0438-0479.201910017 http://dx.doi.org/10.6043/j.issn.0438-0479.201910017 . 气相环境下K+催化丙氨酸分子手性对映体转变的机理 LIU F ， MA H Y ， PAN Y ， et al . Mechanism of alanine chiral enantiomer transition catalyzed by K + in gas phase ［J］. Journal of Xiamen University （Natural Sciences） ， 2020 ， 59 （ 6 ）： 912 - 920 . DOI： 10.6043/j.issn. 0438-0479. 201910017（Ch http://dx.doi.org/10.6043/j.issn.0438-0479.201910017（Ch ）. Mechanism of alanine chiral enantiomer transition catalyzed by K+ in gas phase 姜春旭 ， 陈凤清 ， 张雪娇 ， 等 . 基于密度泛函理论的气相下 α -Ala·Fe 2+ 配合物手性转变机理 ［J］. 复旦学报（自然科学版） ， 2020 ， 59 （ 5 ）： 690 - 701 . DOI： 10.15943/j.cnki.fdxb-jns.2020.05.014 http://dx.doi.org/10.15943/j.cnki.fdxb-jns.2020.05.014 . 基于密度泛函理论的气相下α-Ala·Fe2+配合物手性转变机理 JIANG C X ， CHEN F Q ， ZHANG X J ， et al . Chiral transition mechanism of Fe 2+ and α -Ala complexes in gas phase based on density functional theory ［J］. Journal of Fudan University （Natural Sciences） ， 2020 ， 59 （ 5 ）： 690 - 701 . 690- 701 . DOI： 10.15943/j. cnki.fdxb-jns.2020.05.014 （Ch http://dx.doi.org/10.15943/j.cnki.fdxb-jns.2020.05.014（Ch ）. Chiral transition mechanism of Fe2+ and α-Ala complexes in gas phase based on density functional theory 张雪娇 ， 刘芳 ， 吴梓昊 ， 等 . 丙氨酸Mg 2+ 配合物的手性转变机理水分子（簇）的作用及水溶剂效应 ［J］. 中山大学学报（自然科学版） ， 2020 ， 59 （ 5 ）： 156 - 168 . DOI： 10.13471/j.cnki.acta.snus.2020.02.02.2020c003 http://dx.doi.org/10.13471/j.cnki.acta.snus.2020.02.02.2020c003 . 丙氨酸Mg2+配合物的手性转变机理水分子（簇）的作用及水溶剂效应 ZHANG X J ， LIU F ， WU Z H ， et al . Chiral transition mechanism of Mg 2+ and Ala complexes and effect of water molecules （clusters） and water solvents ［J］. Acta Scientiarum Naturalium Universitatis Sunyatseni ， 2020 ， 59 （ 5 ）： 156 - 168 . DOI： 10.13471/j.cnki.acta.snus.2020.02.02.2020c003（Ch http://dx.doi.org/10.13471/j.cnki.acta.snus.2020.02.02.2020c003（Ch ）. Chiral transition mechanism of Mg2+and Ala complexes and effect of water molecules （clusters） and water solvents 徐锐英 ， 刘芳 ， 马宏源 ， 等 . 气相环境下丙氨酸Ca 2+ 配合物的手性转变机理及水分子的催化作用 ［J］. 浙江大学学报（理学版） ， 2020 ， 47 （ 5 ）： 630 - 641 . DOI： 10.3785/j.issn.1008-9497.2020.05.015 http://dx.doi.org/10.3785/j.issn.1008-9497.2020.05.015 . 气相环境下丙氨酸Ca2+配合物的手性转变机理及水分子的催化作用 XU R Y ， LIU F ， MA H Y ， et al . Chiral transition mechanism of Ala and Ca 2+ complexes in gas phase and catalysis of water molecules ［J］. Journal of Zhejiang University （Science Edition） ， 2020 ， 47 （ 5 ）： 630 - 641 . DOI： 10.3785/j.issn.1008-9497.2020.05.015（Ch http://dx.doi.org/10.3785/j.issn.1008-9497.2020.05.015（Ch ）. Chiral transition mechanism of Ala and Ca2+ complexes in gas phase and catalysis of water molecules 漆剑 . L -脯氨酸和 L -缬氨酸构型转换的研究 ［D］. 南昌 ： 南昌大学 ， 2006 . DOI: 10.15385/yb.miracle.2006 http://dx.doi.org/10.15385/yb.miracle.2006 . L-脯氨酸和L-缬氨酸构型转换的研究 QI J . Study on Configuration Transformation of L - Proline and L - Valine ［D］. Nanchang ： Nanchang university ， 2006 （Ch） . DOI: 10.15385/yb.miracle.2006 http://dx.doi.org/10.15385/yb.miracle.2006 . Study on Configuration Transformation of L-Proline and L-Valine 王佐成 ， 喻小继 ， 闫红彦 ， 等 . 脯氨酸分子的手性转变机理及水的作用 ［J］. 武汉大学学报（理学版） ， 2017 ， 63 （ 4 ）： 317 - 324 . DOI： 10.14188/j.1671-8836.2017.04.006 http://dx.doi.org/10.14188/j.1671-8836.2017.04.006 . 脯氨酸分子的手性转变机理及水的作用 WANG Z C ， YU X J ， YAN H Y ， et al . Chiral enantiomers transformation of Proline molecule and effect of water ［J］. Journal of Wuhan University（Science Edition） ， 2017 ， 63 （ 4 ）： 317 - 324 . DOI： 10.14188/j. 1671-8836. 2017.04.006 （Ch http://dx.doi.org/10.14188/j.1671-8836.2017.04.006（Ch ）. Chiral enantiomers transformation of Proline molecule and effect of water 刘军 ， 潘宇 ， 李加欢 ， 等 . 水液相环境下两性脯氨酸分子旋光异构的理论研究 ［J］. 武汉大学学报（理学版） ， 2020 ， 66 （ 6 ）： 572-2-582 . DOI： 10.14188/j.1671-8836.2020.0167 http://dx.doi.org/10.14188/j.1671-8836.2020.0167 . 水液相环境下两性脯氨酸分子旋光异构的理论研究 LIU J ， PAN Y ， LI J H ， et al . Theoretical study on optical isomerization of amphoteric Proline molecule in aqueous-liquid environment ［J］. Journal of Wuhan University（Science Edition） ， 2020 ， 66 （ 6 ）： 572-2-582 . DOI： 10.14188/j.1671-8836.2020.0167（Ch http://dx.doi.org/10.14188/j.1671-8836.2020.0167（Ch ）. Theoretical study on optical isomerization of amphoteric Proline molecule in aqueous-liquid environment HAMASE K ， INOUE T ， MORIKAWA A ， et al . Determination of free D -proline and D -leucine in the brains of mutant mice lacking D -amino acid oxidase activity .［J］. Anal Biochem ， 2001 ， 298 （ 2 ）： 253 - 258 . DOI： 10.1006/abio.2001.5382 http://dx.doi.org/10.1006/abio.2001.5382 . Determination of free D-proline and D-leucine in the brains of mutant mice lacking D-amino acid oxidase activity 李铭 ， 邹颖 ， 郭寒 ， 等 . 钙离子在内质网应激中的作用 ［J］. 动物医学进展 ， 2018 ， 39 （ 9 ）： 112 - 116 . DOI： 10.16437/j.cnki.1007-5038.2018.09.021 http://dx.doi.org/10.16437/j.cnki.1007-5038.2018.09.021 . 钙离子在内质网应激中的作用 LI M ， ZOU Y ， GUO H ， et al . The role of calcium ions in endoplasmic reticulum stress ［J］. Progress in Veterinary Medicine ， 2018 ， 39 （ 9 ）： 112 - 116 . DOI： 10.16437/j.cnki.1007-5038.2018. 09.021（Ch http://dx.doi.org/10.16437/j.cnki.1007-5038.2018.09.021（Ch ）. The role of calcium ions in endoplasmic reticulum stress 陈瑞仪 ， 谭剑斌 ， 周轶琳 ， 等 . 氨基酸螯合钙增加骨密度作用的实验研究 ［J］. 氨基酸和生物资源 ， 2011 ， 33 （ 1 ）： 60 - 62 . DOI： 10.14188/j.ajsh.2011.01.012 http://dx.doi.org/10.14188/j.ajsh.2011.01.012 . 氨基酸螯合钙增加骨密度作用的实验研究 CHEN R Y ， TAN J B ， ZHOU Y L ， et al . Experimental study on raising bone density by calcium Amino acid chelate compound in rats ［J］. Amino Acids ＆ Biotic Resources ， 2011 ， 33 （ 1 ）： 60 - 62 . DOI： 10.14188/ j.ajsh.2011.01.012（Ch http://dx.doi.org/10.14188/j.ajsh.2011.01.012（Ch ）. Experimental study on raising bone density by calcium Amino acid chelate compound in rats 梁媛 ， 王昕 ， 赵婷婷 ， 等 . 食品中钙吸收的机理与钙体外转化新方法探讨 ［J］. 食品研究与开发 ， 2016 ， 37 （ 4 ）： 195 - 199 . DOI： 10.3969/j.issn.1005-6521.2016.04.048 http://dx.doi.org/10.3969/j.issn.1005-6521.2016.04.048 . 食品中钙吸收的机理与钙体外转化新方法探讨 LIANG Y ， WANG X ， ZHAO T T ， et al . Mechanism and method of calcium absorption in food and discussion of transforming calcium in vitro ［J］. Food Research And Development ， 2016 ， 37 （ 4 ）： 195 - 199 . DOI： 10.3969/j.issn.1005-6521.2016.04.048（Ch http://dx.doi.org/10.3969/j.issn.1005-6521.2016.04.048（Ch ）. Mechanism and method of calcium absorption in food and discussion of transforming calcium in vitro 薛荣涛 ， 李翠芹 ， 何腊平 . 复合氨基酸螯合钙的研究进展 ［J］. 食品工业科技 ， 2014 ， 35 （ 21 ）： 390 - 394 . DOI: 10.13386/j.issn1002-0306.2014.21.076 http://dx.doi.org/10.13386/j.issn1002-0306.2014.21.076 . 复合氨基酸螯合钙的研究进展 XUE R T ， LI C Q ， HE L P . Research progress in compound amino acid chelated calcium ［J］. Science and Technology of Food Industry ， 2014 ， 35 （ 21 ）： 390 - 394 （Ch） . DOI: 10.13386/j.issn1002-0306.2014.21.076 http://dx.doi.org/10.13386/j.issn1002-0306.2014.21.076 . Research progress in compound amino acid chelated calcium STRITTMATTER E F ， LEMOFF A S ， WILLIAMS E R . Structure of cationized glycine， Gly·M 2+ （M=Be， Mg， Ca， Sr， Ba） in the gas phase：Intrinsic effect of cation size on zwitterion stability ［J］. The Journal of Physical Chemistry A ， 2000 ， 104 （ 43 ）： 9793 - 9796 . DOI： 10.1021/jp002970e http://dx.doi.org/10.1021/jp002970e . Structure of cationized glycine， Gly·M2+（M=Be， Mg， Ca， Sr， Ba） in the gas phase：Intrinsic effect of cation size on zwitterion stability REDDY A S ， SASTRY N G . Cation［M=H + ， Li + ， Na + ， K + ， Ca 2+ ， Mg 2+ ， NH +4 ］ inter-actions with the aromatic motifs of naturally oc-curring amino acids： A theoretical study ［J］. The Journal of Chemical Physics A ， 2005 ， 109 （ 39 ）： 8893 - 8903 . DOI： 10.1021/jp0525179 http://dx.doi.org/10.1021/jp0525179 . Cation［M=H+， Li+， Na+， K+， Ca2+， Mg2+， NH+4］ inter-actions with the aromatic motifs of naturally oc-curring amino acids： A theoretical study TRENBERTH K E ， SMITH L . The mass of the atmosphere： A constraint on global analyses ［J］. Journal of Climate ， 2005 ， 18 ： 864 - 875 . DOI: 10.1175/jcli-3299.1 http://dx.doi.org/10.1175/jcli-3299.1 . The mass of the atmosphere： A constraint on global analyses SENNIKOV P G ， IGNATOV S K ， SCHREMS O . Complexes and clusters of water relevant to atmospheric chemistry： H 2 O complexes with oxidants ［J］. ChemPhysChem ， 2005 ， 6 ： 392 - 412 . DOI: 10.1002/cphc.200400405 http://dx.doi.org/10.1002/cphc.200400405 . Complexes and clusters of water relevant to atmospheric chemistry： H2O complexes with oxidants ALOISIO S ， FRANCISCO J S . Radical-water complexes in Earth’s atmosphere ［J］. Accounts of Chemical Research ， 2000 ， 33 ： 825 - 830 . DOI: 10.1021/ar000097u http://dx.doi.org/10.1021/ar000097u . Radical-water complexes in Earth’s atmosphere WANG Y ， VERMA P ， ZHANG L J ， et al . M06-SX screened-exchange density functional for chemistry and solid-state physics ［J］. Proceedings of the National Academy of Sciences of the United States of America ， 2020 ， 117 （ 5 ）： 2294 - 2301 . DOI： 10.1073/pnas.1913699117 http://dx.doi.org/10.1073/pnas.1913699117 . M06-SX screened-exchange density functional for chemistry and solid-state physics GARRETT B C ， TRUHLAR D G . Criterion of minimum state density in the transition state theory of bimolecular reactions ［J］. The Journal of Chemical Physics ， 1979 ， 70 （ 4 ）： 1593 - 1598 . DOI： 10.1063/1. 437698 http://dx.doi.org/10.1063/1.437698 . Criterion of minimum state density in the transition state theory of bimolecular reactions HRATCHIAN H P ， SCHLEGEL H B . Using Hessian updating toincrease the efficiency of a hessian based predictor-corrector reaction path following method ［J］. Journal of Chemical Theory and Computation ， 2005 ， 1 （ 1 ）： 61 - 69 . DOI： 10.1021/ct0499783 http://dx.doi.org/10.1021/ct0499783 . Using Hessian updating toincrease the efficiency of a hessian based predictor-corrector reaction path following method YU H S ， HE X ， LI S L ， et al . MN15： A kohn-sham global-hybrid exchange-correlation density functional with broad accuracy for multi-reference and single-reference systems and noncovalent interactions ［J］. Chemical Science ， 2016 ， 7 （ 8 ）： 5032 - 5051 . DOI： 10.1039/c6sc00705h http://dx.doi.org/10.1039/c6sc00705h .. MN15： A kohn-sham global-hybrid exchange-correlation density functional with broad accuracy for multi-reference and single-reference systems and noncovalent interactions FRISCH M J ， TRUCKS G W ， SCHLEGEL H B ， et al . Gaussian 16 Revision C .01［CP］. Pittsburgh： Gaussian ， Inc， 2019 . Gaussian 16 Revision C GORB L ， LESZCZYNSKI J . Intramolecular proton transfer in mono- and dihydrated tautomers of guanine：  An ab initio post Hartree-Fock study ［J］. Journal of the American Chemical Society ， 1998 ， 120 （ 20 ）： 5024 - 5032 . DOI： 10.1021/ja972017w http://dx.doi.org/10.1021/ja972017w . Intramolecular proton transfer in mono- and dihydrated tautomers of guanine：  An ab initio post Hartree-Fock study 10.14188/j.1671-8836.2020.0287 O641 A 1671-8836（2021）05-0478-11 HUANG Yuchun，SUN Yongxin，ZHANG Xuejiao， et al .Optical Isomerism of Proline Ca（Ⅱ） Complex in Gas Phase and Catalytic Mechanism of Water Molecules（Clusters） ［J］. J Wuhan Univ （ Nat Sci Ed ），2021，67（5）：478-488. DOI： 10.14188/j.1671-8836.2020.0287（Ch http://dx.doi.org/10.14188/j.1671-8836.2020.0287(Ch ）. HUANG Yuchun，SUN Yongxin，ZHANG Xuejiao，et al.Optical Isomerism of Proline Ca（Ⅱ） Complex in Gas Phase and Catalytic Mechanism of Water Molecules（Clusters） ［J］.J Wuhan Univ（Nat Sci Ed），2021，67（5）：478-488. DOI：10.14188/j.1671-8836.2020.0287（Ch）. 黄煜纯,孙永欣,张雪娇, 等.气相脯氨酸Ca（Ⅱ）配合物的旋光异构及水分子（簇）的催化机制[J].武汉大学学报（理学版）,2021,67(5):478-488. DOI:10.14188/j.1671-8836.2020.0287 黄煜纯,孙永欣,张雪娇, 等.气相脯氨酸Ca（Ⅱ）配合物的旋光异构及水分子（簇）的催化机制[J].武汉大学学报（理学版）,2021,67(5):478-488. DOI:10.14188/j.1671-8836.2020.0287 吉林省教育厅科学研究规划项目（JJKH20200002KJ）；吉林省自然科学基金（20130101131JC） 2020-11-26 2021-10-24 2021-10-12 武汉大学学报（理学版） 5 67