陈斌+王石慧+李蓉 樊安

摘要利用前沿色谱法，通过Cu2+、Ni2+和Co2+与螯合配体IDA 键合的准确度（R2>0.98）与精密度（RSD <5%）实验，研究了前沿色谱法同时测定络合稳定常数KML与总键合位点数Λ0值的可行性。为了进一步证明前沿色谱法的普适性，以Cu2+、Ni2+和Co2+为代表，在3种键合缓冲体系（NaAcHAc、NaPB、TrisHCl）中，考察了金属离子在3种氨羧类螯合配体（IDA、Asp、Glu）上络合稳定常数KML的变化规律。结果表明，螯合配体对金属离子键合强度遵循IDA>Asp>Glu；金属离子对螯合配体键合强度遵循Cu2+>Ni2+>Co2+；3种键合缓冲体系中，NaAcHAc键合效果最好。按照M06/6311++G （d， p） 方法对螯合配体与金属离子的结合能（ΔE）与吉布斯自由能（ΔG）进行相关的量子计算。根据ΔE 与ΔG 的大小，从理论上推测出螯合配体与金属离子的键合规律，此规律与上述实验结果基本相符。本研究为金属离子与螯合配体间键合参数的求取提供了有效的方法和手段，从而为今后提高IMAC柱的稳定性，解决固定金属亲和柱在应用过程中尤其是竞争洗脱过程中金属离子的流失问题奠定了良好的基础。

关键词前沿色谱法； 量子计算； 氨羧类螯合配体； 络合稳定常数； 总键合位点数

1引 言

自1975 年Porath 等[1]提出固定金屬离子亲和色谱（Immobilized metal affinity chromatography，IMAC）的概念以来，该技术在蛋白质的分离纯化中得到了广泛应用[2，3]。该方法利用螯合物中固定化的金属离子与蛋白质间的亲和作用，实现目标的分离和纯化[4]。相对于传统的蛋白质分离技术，IMAC 对目标蛋白具有高度的选择性。此外，与其它亲和色谱相比，蛋白质负载量高，更容易实现工业化[5]。然而，金属离子的泄漏是IMAC 柱在实际应用过程中的致命缺陷[6，7]。该泄漏一方面与外部洗脱条件有关，另一方面与IMAC 柱自身即金属离子与螯合配体的键合强度有关，后者起主导作用。键合强度越高，键合稳定常数KML 越大，则IMAC 柱越稳定，在蛋白质竞争洗脱时金属离子越不容易泄漏。因此，建立一种快速准确测定键合强度的方法尤为重要。根据测得的键合参数可有效制备出高强度IMAC 柱，从而缓解甚至解决IMAC 柱上金属离子的流失问题。

键合稳定常数KML和总键合位点数Λ0直接体现金属离子在螯合配体上的键合特性。KML和Λ0的传统检测方法存在诸多不足[8～13]。前沿色谱法（Frontal chromatography，FC）能够在动态条件下同时测定溶质分子在配体上的KML和Λ0值，广泛应用于分子生物、药物机理和药物先导化合物的筛选研究[14～16]。在本课题组前期工作中，提出用FC法进行IMAC柱上KML和Λ0值的测定，结合均匀设计实验，制备出高稳定性能的GluCu柱[17]。

IMAC 柱的稳定性主要同金属离子与螯合配体的键合强弱有关，而键合强弱又与金属离子、络合配基及键合溶液相关。因此，采用FC 法测定IMAC 柱上键合参数是否具有普适性，是本研究的工作重点。根据IMAC 固定金属离子的选择原则（金属离子既要具有适宜的空轨道数目与螯合剂结合，还要有剩余的轨道亲合待分离的蛋白质），本研究以常见的过渡金属离子Cu2+、Ni2+和Co2+为代表，分别考察了在NaAcHAc、NaPB、TrisHCl 的键合溶液体系下，金属离子与IDA、Glu、Asp 的键合参数KML和Λ0。通过分析不同条件下键合参数的变化规律，以及与量子计算结果的比较，论证了方法的可行性。

2实验部分

2.1仪器与试剂

AKTA purifier10 生物色谱制备仪（Amersham Biosciences公司）；124PP 装填机（日本Chemico公司）；纯水仪（西安优普仪器设备公司）。 硅胶（粒径7 μm，孔径30 nm，兰州物理化学研究所）；γ缩水甘油醚氧丙基三甲氧基硅烷（γGLDP，分析纯，辽宁盖县化工研究所）；亚氨基二乙酸（IDA，分析纯，上海国药集团化学试剂有限公司）；谷氨酸（Glu）、天冬氨酸（Asp）均为分析纯；实验用水为超纯水（18.25 MΩ·cm）。

2.2固定相的制备

按照文献[6]分别制得IDA硅胶填料、Asp硅胶填料和Glu硅胶填料。采用匀浆法，将上述不同螯合配体硅胶填料分别装填至规格为50 mm × 4.6 mm的不锈钢柱内，在40 MPa 的压力下装填12 min， 依次得到IDA硅胶柱、Asp硅胶柱和Glu硅胶柱。

2.3KML和Λ0的测定

2.3.1空白突破曲线的测定用5.00 mmol/L NaNO225 mmol/L NaAcHAc 缓冲溶液（pH 4.0）以0.5 mL/min 的流速依次通过未键合金属离子的不同螯合配体硅胶柱，直至在356 nm 处的紫外可见吸收值达到一个平台，即得到空白突破曲线。

2.3.2金属离子突破曲线的测定用0.30 mmol/L CuSO425 mmol/L NaAcHAc 缓冲溶液（pH 4.0）以0.5 mL/min 的流速通过IDA硅胶柱，直至在700 nm 处的紫外可见吸收值达到一个平台，即得到一条金属离子突破曲线，该突破曲线为金属离子在螯合柱上达到饱和吸附量时所对应的平台曲线。依次用0.50， 0.60， 0.80 和1.00 mmol/L 等一系列浓度的CuSO425 mmol/L NaAcHAc（pH=4.0） 重复上述步骤，得到一系列不同浓度Cu2+在IDA硅胶柱上的突破曲线。每测完一条突破曲线，用0.05 mmol/L EDTANa2溶液将键合在色谱柱上的金属Cu2+洗脱。

不同浓度金属Cu2+的其它键合缓冲溶液（NaPB、TrisHCl）分别在Asp硅胶和Glu硅胶柱，以及不同浓度金属Ni2+，Co2+的3种键合缓冲溶液（NaAcHAc、NaPB、TrisHCl）分别在IDA硅胶柱，Asp硅胶柱和Glu硅胶柱上的突破曲线均按上述步骤获得。endprint

2.3.3表观摩尔吸附量的计算用Peakfit（V4.12） 软件分别对空白突破曲线和金属离子突破曲线进行一阶微分处理，得到曲线拐点处所对应的时间，即为空白突破时间to和金属离子突破时间tm。由式（1）计算达到突破时间时的表观摩尔吸附量[18]。

Qe=V×（tm-t0）×[M]60×1000（1）

其中，Qe 为达到平衡时的表观摩尔吸附量（mol）；v 为流速（0.5 mL/min）；tm 为不同浓度金属离子的突破时间（min）；to 为NaNO2 的空白突破时间（min）；[M] 为金属离子的浓度（mol/L）。

2.3.4KML值与Λ0值的求取金属离子在螯合色谱柱上的吸附符合单位点的Langmuir 吸附等温模型[19]：

Qe=KMLΛ0[M]1+KML[M]（2）

将（2）式取倒数可得：

1Qe=1KMLΛ0[M]+1Λ0（3）

其中，KML为金属离子在络合配基上的络合稳定常数（L/mol）； Λ0为金属离子在络合配基上的总键合位点数（mol）。由公式（3）可知1/Qe和1/[M] 呈线性关系，根据斜率和截距可求出KML和Λ0。

3结果与讨论

3.1KML和Λ0的测定

以金属Cu2+为例，在25 mmol/L NaAcHAc键合缓冲溶液（pH=4.0）体系下，按照2.3.2节的方法测得不同浓度Cu2+在IDA 螯合配体上的前沿突破曲线，得到图1的结果。

根据2.3.3节的方法得到不同浓度Cu2+的tm和Qe。由图1和表1可知，随着Cu2+浓度的升高，tm 减小，Qe 增大。表明当Cu2+浓度升高时，有利于Cu2+在IDA硅胶柱上的吸附，短时间内金属Cu2+在IDA硅胶柱上迅速达到吸附平衡，且吸附量逐渐增大。

根据公式（3）将1/Qe 和1/[M] 进行线性拟合，结果如图2所示。 Cu2+浓度在0.30～1.00 mmol/L 的范围内，1/Qe与1/[M] 呈良好的线性相关性，线性相关系数R2=0.985。这充分说明了在NaAcHAc 的键合缓冲体系中，金属Cu2+在IDA螯合配体上的吸附特性符合单位点的Langmuir 吸附等温模型。根据1/Qe与1/[M] 的线性方程可以得到在25 mmol/L NaAcHAc键合缓冲体系（pH=4.0）的条件下，金属Cu2+在IDA 螯合柱上的络合稳定常数KML=1.12×104 L/mol， 键合位点数 Λ0=9.31×10

Symbolm@@ 6。

3.2FC方法的准确度与精密度实验

为了考察本法的准确性与可靠性，在25 mmol/L NaAcHAc（pH=4.0）键合缓冲溶液的体系下，对不同的金属离子，在IDA螯合柱上进行吸附特性的研究，结果如表2所示。

由表2可知，3种金属离子（Cu2+、Ni2+和Co2+）浓度在0.30～1.00 mmol/L范围内，1/Qe和1/[M] 的线性关系良好，线性相关系数R2>0.98；吸附特性参数KML和Λ0的RSD分别低于2.7% 和3.9%。说明方法的准确度高，重现性好。这是因为IDA是一种含有3个配位原子的氨羧类络合剂，可以与Cu2+、 Ni2+和Co2+以1∶1的络合比进行络合，金属离子在其上的吸附符合单位点的Langmuir 吸附等温模型。此外，Cu2+、Ni2+和Co2+与IDA 的络合稳定常数遵循如下顺序KML（CuIDA） > KML（NiIDA） >KML（CoIDA），此结果与理论推测相一致[20～22]。

由此可见，利用FC 法可准确有效地求得表征金属离子在氨羧类螯合剂上键和特性的参数KML和Λ0， 并且通过评价参数KML直观地反映出不同的金属离子与螯合剂结合能力的强弱，这也为分析金属离子与螯合剂的键合强度以及间接地判断IMAC 柱的稳定性提供了依据。

3.3金属离子在氨羧类螯合配体上的键合特性

为了进一步论证FC 法的通用性， 将研究范围扩大到Glu，Asp螯合配体以及NaPB 和TrisHCl 键合缓冲体系中，依次考察了Cu2+、 Ni2+和Co2+的键合特性，结果如表3所示。

从表3可知，在相同的键合缓冲体系下：（1）3种金属离子对不同螯合配体柱的络合稳定常数遵循：KML（CuIDA） > KML（CuAsp） >KML（CuGlu）；KML（NiIDA）>KML（NiAsp）>KML（NiGlu）；KML（CoIDA） >KML（CoAsp） >KML（CoGlu），均表现为螯合配体对金属离子键合强度：IDA>Asp>Glu。这一规律与本课题组先前用AAS法测得的结果相符[23，24]；（2）3种金属离子对相同螯合配体柱的络合稳定常数遵循：KML（CuIDA） > KML（NiIDA） >KML（CoIDA）；KML（CuAsp） > KML（NiAsp） >KML（CoAsp）；KML（CuGlu） > KML（NiGlu） >KML（CoGlu），均表现为金属离子对螯合配体键合强度：Cu2+>Ni2+>Co2+。这是因为Cu2+、Ni2+和Co2+的半径分别为69、72和74 pm，当金属离子的电荷数相等时，随着金属离子半径的依次增大高，金属离子与络合剂间的配位作用會逐步减弱[25]。 此外，3种键合缓冲体系中，NaAcHAc 的键合效果最好，这与文献[25]的研究结论一致。

为了进一步论证实验结果的可靠性，用Gaussian 09 软件对上述金属离子与螯合配体进行量子模拟计算：采用B3LYP 和631G（d） 基组进行几何构型优化，在M06/6311++G（d， p）水平计算螯合配体与金属离子的结合能（ΔE）与吉普斯自由能（ΔG），结果见表4。ΔE 在一定程度上体现了金属离子与螯合配体构型的稳定性：ΔE 越高， 说明这种构型在配位过程中越容易产生且能稳定存在；ΔG 反映了反应正向进行的趋势， 其为负值说明反应能自发进行，且数值的绝对值越大，螯合物构型越稳定，说明金属离子与螯合配体的键合力越大。endprint

由表4 可知，ΔE>0， ΔG<0 ， 说明这9种螯合物都能稳定存在。对于同一种螯合配体而言，金属离子与其的键合能力均表现为Cu2+>Ni2+>Co2+。这是因为当金属离子所帶电荷相同时，螯合物的稳定性与金属离子的半径有关[24～26]：金属离子的半径越小，螯合物的ΔE 和ΔG （绝对值）越大，则结构更加稳定。因此量子计算的结果与实验结果一致。

金属离子与螯合配体是按照1∶1 的络合比形成稳定的螯合物[27]，因此不同配体的空间构型将对金属离子与其螯合的稳定性产生一定的影响。由螯合配体的结构可知，对于提供同样配位原子的螯合体，

主链上含有配位原子的螯合体能够与金属离子形成更加稳定的结构；此外，螯合配体的链长也会影响其与金属离子的键合能力。计算结果表明，对于同一种金属Cu2+而言，IDA与其的键合能力最强；与Asp 与Glu相比，由于IDA 主链上含有提供配位的N原子，使得金属Cu2+与其配位后形成的络合物更加稳定[27]。而Asp与Glu 结构相似，主链上含有更多的CH2， 使得这两种螯合配体与金属离子相结合时，空间位阻较大，导致形成的螯合物的稳定性较弱[28，29]。而对于Ni2+和Co2+而言，计算结果与实验结果有偏差，这可能是因为理论的计算是在水溶液进行的，忽略了键合缓冲溶液的影响。

4结 论

本研究在FC 法的基础上，利用Langmuir 单位点的等温吸附模型，考察了在不同缓冲体系条件下，不同金属离子与不同螯合配体间的键合稳定参数之间的关系。通过数据分析与比较，论证了FC 法在IMAC 制备过程中金属离子对螯合配体键合参数测定的普适性。本研究为金属离子与螯合配体间键合参数方便快速的求取提供了有效的方法和手段，从而为今后制备高强度IMAC 柱，解决IMAC 柱在使用过程中金属离子的流失问题奠定了良好的基础。

References

1Porath J， Carlsson J， Olsson I， Belfrage G. Nature， 1975， 258（5536）： 598-599

2WANG Yong. J. Graduates. Sun YATSEN. Univ： Nat. Sci. Med. Ed.， 2012， 33（2）： 15-21

王 勇. 中山大学研究生学刊： 自然科学与医学版， 2012， 33（2）： 15-21

3Lyu Y， Bao X L， Liu H， Ren J H， Guo S T. Food Chem.， 2013， 141（3）： 1645-1650

4 YU ShiLin. Affinity Chromatography and Application. Beijing： Chemical Industry Press， 2008： 235-241

于世林. 亲和色谱方法及应用. 北京： 化学工业出版社， 2008： 235-241

5Papaneophytou C P， Kontopidis G. Protein Expr. Purif.， 2014， 94（2）： 22-32

6Li R， Wang Y， Chen G L， Shi M， Wang X G， Chen B， Zheng J B. Sep. Sci. Technol.， 2010， 46（2）： 309-314

7GabercPorekar V， Menart V. Chem. Eng. Technol.， 2005， 28（11）： 1306-1314

8WANG Yan， LI Rong， CHEN GuoLiang， WANG XiaoGang. Chemistry， 2008， 11： 838-845

王 燕， 李 蓉， 陈国亮， 王小刚. 化学通报， 2008， 11： 838-845

9Harrington J M， Jones S B， Hancock R D. Inorg. Chim. Acta， 2005， 358（15）： 4473-4480

10CHEN Ying， ZHANG ManLi， ZHANG Wei， HE HongBo， GUAN LianZhu， YAN Li. Journal of Liaoning Technical University （Natura Science Edition）， 2008， 27（3）： 478-480

陈 盈， 张满利， 张 威， 何洪波， 关连珠， 颜 丽. 辽宁工程技术大学学报（自然科学版）， 2008， 27（3）： 478-480

11Byrne L A， Hynes M J， Connolly C D， Murphy R A. J. Inorg. Biochem.， 2011， 105（12）： 1656-1661

12Liu X， Wang H， Liang A， Li Y， Gai H， Lin B. J. Chromatogr. A， 2012， 1270（24）： 340-343

13LIU YingHong， SHA Ou， WANG GuiFen. J. Huaihai. Inst. Technol（Nat.Sci.Ed）.， 2005， 14： 56-58

刘英红， 沙 鸥， 王桂芬. 淮海工学院学报（自然科学版）， 2005， 14： 56-58

14Hage D S， Schiel J E. Encyclopedia Chromatogr.， 2009： 184-191

15Bian L J， Qian L， Xu J. J. Chromatogr. A， 2015， 1401： 42-51endprint

16Calleri E， Temporini C， Massolini G. J. Pharmaceut. Biomed.， 2011， 54（5）： 911-925

17Chen B， Li R， Li S H， Chen X L， Yang K D， Chen G L， Ma X X. J. Sep. Sci.， 2016， 39（3）： 518-524

18Bian L J， Geng X D. Chin. J. Chrom.， 1994， 12（5）： 317-320

19Iwata H， Saito K， Furusaki S， Sugo T， Okamoto J. Biotechnol. Prog.， 1991， 7（5）： 412-418

20Sengupta S， Sengupta A K. Environ. Sci. Technol.， 1993， 27（10）： 2133-2140

21Foresti M L， Nyholm L. J. Electroanal. Chem.， 1989， 269（1）： 2144-2140

22LIN XiuLi， ZHU ChenFu， ZOU ShiFu. Chinese J. Anal. Chem.， 1996， 24（2）： 175-179

林秀麗，主沉浮，邹时复. 分析化学， 1996， 24（2）： 175-179

23Li R， Wang Y， Chen G L， Wang X G， Zheng J B. Chromatographia， 2009， 70（56）： 731-737

24Li R， Ju M Y， Chen B， Sun Q Y， Chen G L， Shi M， Wang X G， Zheng J B. Korean Chem. Soc.， 2011， 32（2）： 590-594

25Sharma S， Agarwal G P. Sep. Sci. Technol.， 2002， 37（15）： 3491-3511

26ZHANG JiXin， SU KeMan， WANG Yan， HU Ping. Analytical Chemistry （6th ed）. Beijing： Higher Education Press， 2009： 102-126

张济新， 苏克曼， 王 燕， 胡 坪. 分析化学（第六版）. 北京： 高等教育出版社， 2009： 102-126

27Shircliff A D， Wilson K R， CannonSmith D J， Jones D G， Zhang Z， Chen Z Q， Yin G C， Prior T J， Hubin T J. Inorg. Chem. Commun.， 2015， 59： 71-75

28Sharma R P， Saini A， Venugopalan P， Jezierska J， Ferretti V. Inorg. Chem. Commun.， 2012， 20（6）： 209-213

29Sóvágó， Vrnagy K. Met. Ions Life Sci.， 2013， 11（11）： 275-302

Study on Binding Characteristics of Aminocarboxyl

Chelating Ligands for Metal Ions

CHEN Bin， WANG ShiHui， LI Rong*， FAN An

（School of Chemical Engineering， Northwest University， Xi'an 710069， China）

AbstractThe feasibility of frontal chromatography for determining the complexation stability constant KML and total mole of binding site Λo was demonstrated by the accuracy and precision binding experiments between metal ions （Cu2+， Ni2+ and Co2+） and chelating ligand （IDA）， in which R2>0.98 and RSD <5%. To further prove the universality of the frontal chromatography， the changing rules of KML values between Cu2+， Ni2+， Co2+ and IDA， Asp， LGlu were examined under NaAcHAc， NaPB and TrisHCl buffer systems. The results showed that the binding strength of chelating ligands for metal ions followed IDA>Asp>Glu； binding strength of metal ions for chelate ligands followed Cu2+>Ni2+>Co2+； and the binding effect with NaAcHAc buffer was the best. In aqueous solution， quantum computing of binding energy （ΔE） and gibbs free energy （ΔG） between chelating ligand and metal ion was performed at the M06/6311++G （d， p） level. According to ΔE and ΔG， the binding rules between chelating ligand and metal ion were predicted theoretically. These rules were basically in agreement with above experimental results. The present work provided effective method for studying on binding characteristics of metal ions for aminocarboxyl chelating ligands， thus exhibited a good foundation for improving the stability of immobilized metal affinity chromatographic column and solving the leakage of metal ions from the column in the process of competitive elution.

KeywordsFrontal chromatography； Quantum computing； Aminocarboxyl chelating ligand；Complexation stability constant； Total mole of binding site

（Received 15 March 2017； accepted 22 June 2017）

This work was supported by the National Natural Science Foundation of China （No.21376191） and the Natural Science Foundation of Shaanxi Province， China （No. 2011JM2013）， the Service Local Special Project of Shaanxi Province Education Department （No.14JF027） and the Industrialization Cultivation Item of Shaanxi Province Educational Department （No.2013jc23）endprint