大孔吸附树脂对大黄素吸附性能的研究

期刊菜单

大孔吸附树脂对大黄素吸附性能的研究
Study on the Adsorption Performance of Macroporous Adsorption Resin for Emodin

DOI: 10.12677/hjcet.2024.146050, PDF, HTML, XML,
作者: 王盼盼, 吴慧敏, 黄永腾：湖北科技学院医学部药学院，湖北咸宁；钱文斌^*：湖北科技学院医学部基础医学院，湖北咸宁
关键词: 大孔吸附树脂；大黄素；吸附行为；吸附动力学；吸附等温线；Macroporous Adsorption Resins； Emodin； Adsorption Behavior； Adsorption Kinetics； Adsorption Isotherms

摘要: 通过研究大孔树脂对大黄素的吸附行为，为工业化分离大黄素提供一定的理论指导。采用静态吸附实验，研究不同大孔树脂对大黄素的吸附能力；选用较好吸附树脂考察其对大黄素的吸附热力学及动力学特性。H103树脂吸附能力优于其他4种大孔树脂，平衡吸附量为78 mg/g；通过分子对接吸附材料H103树脂主要结构与大黄素分子的芳香环间形成2个弱作用力π-π键，键长分别为2.9 Å、2.64 Å；Freundlich模型、Dubinin Radushkevich模型、Temkin模型能较好地模拟其等温吸附过程；吸附焓变Δ H > 0、吸附自由能Δ G < 0、吸附熵变Δ S > 0；吸附过程为自发的物理反应过程；大孔树脂对大黄素的吸附在130 min后达到平衡吸附过程符合伪一级动力学方程。H103大孔树脂可用于吸附大黄素；热力学和动力学研究结果为分离大黄素提供理论指导和科学依据。

Abstract: Studying the adsorption behavior of macroporous resin on emodin provides some theoretical guidance for the industrial separation of emodin. The static adsorption experiment was used to study the adsorption capacity of different macroporous resin, and to investigate the thermodynamic and kinetic characteristics of emodin. The adsorption capacity of H103 resin is better than the other four macroporous resin, and the equilibrium adsorption capacity is 78 mg/g; Two weak force π-π bonds are formed between the main structure of H103 resin and the aromatic ring of emodin molecule, and the bond lengths are 2.9 Å and 2.64 Å respectively; Freundlich model, Dubinin Radushkevich model and Temkin model can better simulate the isothermal adsorption process; the adsorption enthalpy is Δ H > 0, the adsorption free energy is Δ G < 0, and the adsorption entropy changes Δ S > 0; the adsorption process is a spontaneous physical reaction process; the adsorption of macroporous resin reaches the equilibrium adsorption process after 130 min conforms to the pseudo first-level kinetic equation. H103 macrohole resin can be used to absorb emodin; thermodynamic and kinetic study results provide theoretical guidance and scientific basis for separating emodin.

文章引用：王盼盼, 钱文斌, 吴慧敏, 黄永腾. 大孔吸附树脂对大黄素吸附性能的研究[J]. 化学工程与技术, 2024, 14(6): 460-472. https://doi.org/10.12677/hjcet.2024.146050

1. 引言

大黄素羟基蒽醌类化合物，其化学名为1,3,8-三羟基-6-甲基蒽醌，结构见图1。大黄素被用作传统药物已有2000多年的历史，其在800多种中药制剂中普遍存在，其在虎杖科、豆科和鼠李科三个植物科中含量较丰富[1]。大黄素具有多种药理作用包括其作为抗癌剂、抗炎、抗病毒、和抗氧化的作用，而在过去的COVID-19 [2]中利用分子对接揭示了大黄素治疗的靶点和药理学机制。通过研究的不断深入大黄素成为中药制剂和膳食补充剂中的宝贵成分[3]。

Figure 1. 1,3,8-trihydroxy-6-methylanthraquinone

图1. 1,3,8-三羟基-6-甲基蒽醌

大黄素的纯化方法有机溶剂萃取法、柱层析法、聚酰胺法、葡聚糖凝胶法、离子交换树脂法等，这些方法存在毒性、对操作人员和设备要求高，得到的纯品大黄素产量较低[4]，和目前市场随着需求的增大在加工和处理残渣或清洗提取设备过程中会产生大量高浓度、高色度废水，未经处理直接排放将对生态环境带来极大危害[5]。而大孔吸附树脂(MARs)是一种功能高分子材料，与其他吸附剂相比，具有结构多样性、成本低、耐酸碱性好、多孔、比表面积高、环境友好、使用寿命长等优点，但目前尚无吸附大黄素理论方面的研究，限制了大黄素分离纯化研究的发展。本实验通过研究SD300、D101、X-5、H103、SV300大孔吸附树脂对大黄素的吸附容量，筛选出吸附效果最好大孔吸附树脂，然后通过分子对接预测大孔树脂主要结构与大黄素分子间的作用力，使用Gaussian软件通过量子化学理论计算大黄素分子的大小和比表面积测试(BET)实验计算树脂孔径的大小，从理论上预测大孔树脂是否符合吸附要求。随后进行吸附热力学与动力学考察，做出大黄素分子在树脂表面的吸附等温线，结合Langmuir与Freundlich等模型和热力学参数判断大黄素分子在大孔吸附树脂表面的吸附方式，和通过伪一级动力学方程、伪二级动力学方程和颗粒内扩散方程对吸附过程进行拟合，分析H103大孔树脂的吸附机理及大黄素分子在颗粒内的扩散行为，为后续大黄素分离纯化工艺研究提供参考，为资源的合理利用提供科学依据。

2. 实验部分

2.1. 材料、试剂与仪器

材料：大黄素(九州康源生物科技，质量分数 ≥ 95%)；大孔树脂SD300、D101、X-5、H103、SV300，浙江争光实业股份有限公司，杭州慧赢仪器科技有限公司。

试剂：乙醇、盐酸、氢氧化钠等均购自国药集团化学试剂有限公司，所用试剂若无特殊说明均为分析纯。

仪器：UV-2100可见光分光光度计，尤尼可(上海)仪器有限公司；分析天平、电极，瑞士梅特勒公司；超声波清洗器，深圳市方奥微电子有限公司；多功能恒温水浴振荡器，上海一恒仪器有限公司。

2.2. 方法

2.2.1. 大黄素分子尺寸的计算

使用Gaussian软件通过量子化学理论计算，依据密度泛函理论(DFT)中的B3LYP [6] [7]方法，并且考虑了D3 (BJ)色散校正，并且考虑了溶剂化效应对结构的影响，对优化好后的结果进行频率计算，确保每个结构的频率为正，优化后的大黄素分子构型见图2所示。通过Multiwfn 3.8程序计算大黄素分子的长宽高数值(见表1)。

Figure 2. The optimized molecular structure of emodin based on DFT

图2. 基于DFT优化后的大黄素分子结构

Table 1. Molecular structures and physicochemical properties of the adsorbates

表1. 吸附物的分子结构和理化性质

吸附物	分子式	分子量(g/mol)	摩尔体积(cm³/g)	分子尺寸(Å)
Emodin	C₁₅H₁₀O₅	270.24	188.541	长：1.37 宽：0.85 高：0.41

2.2.2. 标准曲线的绘制

参考文献[8] [9]的方法作适当修改，精密称取大黄素标准品4.5 mg置于100 mL容量瓶中，加乙醇溶解并摇匀、定容，配成0.045 g/L的母液。依次配置成质量浓度为0.005、0.009、0.015、0.030、0.045 g/L的溶液，分别移取1 mL不同浓度的标准品溶液到10 mL的容量瓶中，加入配置的5 mol/L的氢氧化钠溶液1 mL，去离子水定容至刻度并摇匀放置5 min。5 mol/L的氢氧化钠溶液1 mL，去离子水定容至刻度的空白溶液为对照在530 nm处用紫外分光光度计测定并记录其吸光度值，以大黄素的含量(mg/L)为横坐标，吸光度值(A)为纵坐标，绘制标准曲线，得到线性回归方程如图3，A = 0.0034 × C + 0.014，R² = 0.998；表明大黄素含量在0.001~0.045 g/L范围内线性关系良好。

Figure 3. Standard curve of the emodin solution

图3. 大黄素溶液的标准曲线

2.2.3. 树脂预处理

取适量大孔吸附树脂SD300、D101、H103、X-5、SV300，用无水乙醇浸泡24 h，去除上层悬浮颗粒，用蒸馏水洗至无醇味，再用1.5 mol/L盐酸洗、蒸馏水洗、1.5 mol/L氢氧化钠洗、蒸馏水洗至中性具体树脂类型及其参数见表2。

Table 2. Resin types and their parameters

表2. 树脂类型及其参数

类型		极性	比表面积 (m²/g)	平均孔径 (nm)	粒径范围 (mm)
大孔吸附树脂	SD300	非极性	900~1200	3~3.5	0.3~1.25
	D101	极性	500~550	9~10	0.2~0.84
	X-5	非极性	500~600	29~30	0.2~0.84
	H103	非极性	900-1100	8~9	0.3~1.25
	SV300	非极性	800-1000	2.8~3.5	0.31~1.25

2.2.4. 吸附物溶液的制备

大黄素理化性质为几乎不溶于水溶于乙醇与碱溶液，配置浓度为8 g/L吸附物大黄素溶液150 mL，分别用水，5%、25%、75%的乙醇并调pH为9.0、10.0、11.0、12.0作为吸附物溶液置于具塞三角瓶中，称取已经预处理并抽干的H103大孔树脂9 g，温度293.15 K振荡速度150 r/min，恒温振荡吸附8 h，使之充分吸附，取上清液按“2.2”中相同的操作计算大黄素含量，代入公式(1)计算平衡吸附量。

$q_{e} = \frac{C_{0} - C_{e}}{W} \times V$ (1)

其中q_e为平衡吸附量，mg/g；V为实验中加入溶液的体积，mL；W为干树脂的重量，g。

2.2.5. 大孔树脂筛选

分别精密称取前面处理好的5种大孔树脂SD300、D101、H103、X-5、SV300各9 g，置于250 mL具塞三角瓶中，加入准确配置的浓度为8 g/L的大黄素溶液150 mL (用体积分数25%的乙醇溶液配制，并调节pH为11.0)，温度293.15 K下振荡速度设置为150 r/min的条件下恒温振荡器中避光振摇8 h，使之充分吸附，取上清液按“2.2”中相同的操作计算大黄素含量，代入公式(1)计算平衡吸附量。

2.2.6. 静态吸附热力学实验

称取已经预处理并抽干的H103大孔树脂9 g，配置溶液浓度1、4、8、10 g/L共4个浓度150 mL置于具塞三角瓶中(用体积分数25%的乙醇溶液配制，并调节pH为11.0)放入恒温振荡器中，分别在温度283.15、293.15、303.15 K下以150 r/min振荡吸附8 h，振荡吸附结束后取上清液加水稀释后定容，代入标准曲线测定其吸附前后溶液大黄素浓度，按公式(1)计算平衡吸附量q_e。

2.2.7. 吸附动力学实验

称取已经预处理并抽干的H103大孔树脂9 g，将溶液浓度8 g/L大黄素溶液150 mL置于具塞三角瓶中(用体积分数25%的乙醇溶液配制，并调节pH为11.0)，温度293.15 K、振荡速度150 r/min，恒温振荡吸附8 h，于不同时间间隔(5、10、16、24、42、110、230、320、350 min)振摇暂停，取上清液加水稀释后定容，按上述同法于530 nm处测定其吸光度值，按公式(2)计算t时间的吸附量(q_t)的关系。

$q_{t} = \frac{C_{0} - C_{t}}{W} \times V$ (2)

其中q_t为t时间大孔树脂的吸附量，mg/g；C_t为t时间溶液中大黄素的浓度，g/L。

3. 结果与分析

3.1. 吸附物溶液的选择与树脂的筛选

在四种不同吸附物溶液与pH条件下大黄素的吸附量结果如图4，在水与5%的乙醇溶液中大黄素的吸附量较差。随着乙醇浓度与pH值的增加吸附量也随之增加，但在50%乙醇与25%乙醇并调节pH为11.0与12.0的溶液中大黄素的吸附量变化不大，分别为78.4、79.1、79.4、79.9 mg/g，从节能与保护环境方面考虑溶液使用25%乙醇并调节pH为11.0此条件下进行后续实验。不同型号树脂对大黄素溶液静态平衡吸附量的影响如图5，结果表明大孔树脂H103的平均吸附量最大，符合树脂和大黄素的性质。

H103属于非极性聚苯乙烯型大孔树脂其吸附原理主要靠物理吸附分子间作用力的影响其特点吸附物质不发生任何化学反应，吸附过程进行得较快。将大黄素分子作为受体，H103树脂主要结构聚苯乙烯运用Chemdraw画出化合物的二维结构，随后导入到Chemdraw 3D中利用MM2模块进行能量最小化并保存为mol2文件作为配体，然后用Mgtools 1.5.6通过去水、加氢、计算电荷、合并非极性氢等处理后通过Autodock vina 1.1.2软件进行分子对接，利用pymol进行可视化如图6，结果表明结合能为−4.0 kcal/mol，其结合能小于0表示大黄素分子能够自发性结合到吸附材料中，随后通过分析三维相互作用，吸附材料H103树脂主要结构与大黄素分子的芳香环间形成2个弱作用力π-π键，键长分别为2.9、2.64 Å。

Figure 4. Content of macroflavin in different adsorbent solutions

图4. 不同吸附物溶液下大黄素的含量

Figure 5. Equilibrium adsorption amount of different macroporous resins

图5. 不同大孔树脂的平衡吸附量

Figure 6. Molecular docking was visualized by using pymol

图6. 分子对接利用pymol进行可视化

H103树脂通过分子间弱相互作用与大黄素分子内的疏水部分进行吸附最适于吸附极性小的分子，而大黄素分子存在极性小的问题[10]而H103型树脂比表面积较大，有利于大黄素分子灵活通过孔径被吸附。通过比表面积测试(BET)实验进行验证，根据各参数得出H103树脂在77 K的氮气吸附/解吸附等温线和代入BJH模型分析孔径分布图分别如图7、图8所示，其等温线属于IUPAC分类中的IV型[11]，吸附过程为介孔壁上初发生的单层到多层吸附是工业常出现的吸附过程。孔径在2~50 nm范围的孔称为介孔、将孔径大于50 nm的孔称为大孔，BJH模型计算得H103树脂中主要存在98%、1.99%介孔和大孔属于介孔大孔吸附剂见表3，计算得平均孔径为13.8 nm与表2所示树脂参数相接近。根据2.1所述方法采用Gaussian软件计算可得，大黄素分子的长为1.37 Å，通常认为当吸附剂的孔径是吸附质分子尺寸的2~6倍时，比较有利于吸附[12]。本文吸附剂孔径是大黄素分子直径的100倍，因此H103树脂可适合于吸附大黄素，因而表现出最佳的吸附性能。

Table 3. Aore types and volumes of BET experiments

表3. BET实验孔径类型与体积

孔径类型	微孔	介孔	大孔	总孔径
孔径体积(cm²/g)	0.004	1.262	0.009~0.026	1.275~1.292

Figure 7. N₂ adsorption and desorption of H 103 resin

图7. H103树脂N₂吸附和解吸等温线

Figure 8. The pore size distribution

图8. 孔径分布

3.2. 吸附热力学

3.2.1. 等温吸附模型

图9表明随着温度的升高，q_e先增大后减小，在实验过程中适当增加温度更有利于吸附，温度在303.15 K与293.15 K时的吸附量较为接近，从节能上考虑选择吸附温度设置为293.15 K最为适宜。

Figure 9. Adsorption isotherm line

图9. 吸附等温线

为了进一步研究，将利用不同吸附等温线数学模型进行拟合分析。本研究采用Langmuir、Freundlich、Dubinin-Radushkevich和Temkin等温模型[13]-[15]对吸附实验所得数据进行线性拟合，4种模型方程见公式(3)~(6)。

Langmuir吸附等温方程：

$\frac{C_{e}}{q_{e}} = \frac{C_{e}}{q_{m}} + \frac{1}{K_{L} q_{m}}$ (3)

其中C_e为振荡吸附后溶液中的大黄素浓度，g/L；q_e为吸附平衡后树脂的平衡吸附量，mg/g；K_L为树脂与大黄素间的结合常数，L/g；q_m为H103树脂饱和吸附量，mg/g。

Freundlich吸附等温方程：

$\ln q_{e} = n \ln C_{e} + \ln K_{F}$ (4)

其中K_F为平衡吸附系数，表示吸附量的大小，1/n为特征常数，表明吸附剂表面的不均匀性和相对强度的大小。

Dubinin-Radushkevich吸附等温方程：

$\ln q_{e} = \ln q_{m} - K_{D R} ε^{2}$ (5)

其中K_DR为与吸附能量有关的常数，mol²/J²；ε为Polanyi潜能， $ε = R T \ln (1 + C_{e}^{- 1})$ ；R为气体常数，8.314 J/(mol∙K)；T为温度，K。

Temkin吸附等温方程：

$q_{e} = \frac{R T}{b_{T}} \ln A_{T} + \frac{R T}{b_{T}} \ln C_{e}$ (6)

其中A_T、b_T分别是热力学常数、Temkin方程常数。

由下表4可知在Freundlich方程中R²值均大于0.99，的参数K_F与n是从 $\ln q_{e}$ 对 $\ln C_{e}$ 所得的线性回归方程中计算得到的，K_F是与吸附剂的平衡吸附容量有关的常数，其值与实验中大黄素的吸附容量数据具有相同的趋势，拟合参数K_F值随温度的提高而下降，表明该过程为放热反应，不同温度下其拟合参数1/n值均小于1，说明该吸附过程为优惠吸附[16]。

由表4可知，在Dubinin Radushkevich方程中R²值均大于0.95说明大黄素在H103树脂上的吸附数据可以利用Dubinin Radushkevich等温线模型较好地进行拟合，通过常数K_DR计算出每摩尔吸附质吸附的平均自由能 $E = 1 / \sqrt{2 K_{D R}}$ (283.15、293.15、303.15 K，E分别为1.13、1.31、1.60 kJ/mol)，自由能均小于8 kJ/mol [17]，说明大黄素在H103大孔上的吸附属于简单的物理吸附。

由表4可知，在Temkin方程[18]中R²值均大于0.96说明大黄素在H103树脂上的吸附数据可以利用Temkin等温线模型较好地进行拟合[19]。吸附热(b_T)的值较大，表明利用H103树脂吸附大黄素是可行的，大黄素H103树脂对大黄素的亲和力也较好。

综上Freundlich方程、在Dubinin Radushkevich方程、Temkin方程的相关系数R²值均大于0.95，而Freundlich方程的R²值均大于0.99能相对较好地拟合吸附过程，由模型可知H103大孔树脂对大黄素的亲和力较好且属于简单的物理吸附，其过程为放热反应属于优惠吸附。

Table 4. Fitting parameters of H 103 resin

表4. H103树脂对大黄素吸附等温的拟合参数

T (K)	Langmuir模型			Freundlich模型
T (K)	K_L (L/g)	q_m (mg/g)	R²	K_F	1/n	R²
283.15	0.81	68.11	0.93	44.42	0.31	0.99
293.15	0.94	78.17	0.96	38.99	0.37	0.98
303.15	1.03	83.73	0.94	24.52	0.42	0.99
T (K)	Dubinin Radushkevich模型			Temkin模型
T (K)	K_DR (10⁻⁷ mol²/J²)	q_m (mg/g)	R²	b_T	A_T	R²
283.15	3.91	54.59	0.95	116.67	2.38	0.98
293.15	2.99	73.67	0.98	118.77	6.29	0.96
303.15	1.84	70.81	0.95	187.81	41.26	0.95

3.2.2. 吸附热力学参数

热力学理论有三个重要参数，它们分别是标准吉布斯自由能变(ΔG₀)、标准焓变(ΔH⁰)、标准熵变(ΔS⁰)，它们对于揭示吸附过程中的能量变化有重要指导意义，公式(7)~(9)计算得出ΔG⁰、ΔH⁰、ΔS⁰。

$Δ G^{0} = - R T \cdot \ln K_{0}$ (7)

$\ln K_{0} = - \frac{Δ H^{0}}{R T} + \frac{Δ S^{0}}{R}$ (8)

$Δ S^{0} = \frac{Δ H^{0} - Δ G^{0}}{T}$ (9)

其中K₀为热力学平衡常数；通过ln(q_e/C_e)对q_e作图所得截距即为K₀值；T为绝对温度，K；R为气体常数。以lnK₀对1/T作图，并进行线性拟合，再根据公式(9)计算出吸附焓变ΔH⁰。

由表5可知吉布斯自由能ΔG⁰均为负值，说明该吸附过程是自发的，且容易发生[20]。研究表明∆H⁰的绝对值小于43 kJ/ mol时，吸附过程为物理吸附，说明该吸附可能受物理机制控制，而不是化学过程，这与H103型树脂吸附动力工作原理一致。熵变值∆S⁰为161.79 (J/(mol∙K))结果为正值，表明大黄素与H103树脂的亲和力较好，同时大黄素在吸附过程中的混乱程度是增大的[21] [22]。

Table 5. Thermodynamic parameters of emodin adsorption by H 103 resin

表5. H103树脂对大黄素吸附的热力学参数

T (K)	ΔG⁰ (kJ/mol)	∆H⁰ (kJ/mol)	∆S⁰ (J/(mol∙K)
283.15	−21.44	24.35	161.79
298.15	−24.03
303.15	−24.61

3.3. 吸附动力学

将数据进行分析，通过拟合常用伪一阶动力学模型和伪二阶动力学模型，Weber-Morris颗粒内扩散模型公式(10)~(12)，阐述大黄素在H103树脂上的吸附机理，并研究哪一种模型更适合描述吸附效果，模型方程如下在开始时刻t = 0时有q_t = 0以及任意时刻t = t时都有q_t = q_t的边界条件下进行整合方程式并化简得：

伪一级动力学模型的方程为：

$\frac{d q}{d t} = k_{1} (q_{e} - q_{t})$ 整合后 $\ln (q_{e} - q_{t}) = \ln q_{e} - k_{1} t$ ，通过两边取对数得

$q_{t} = q_{e} (1 - e^{- k t})$ (10)

伪二级动力学模型的方程为：

$\frac{d q}{d t} = k_{2} {(q_{e} - q_{t})}^{2}$ 整合后 $\frac{t}{q_{t}} = \frac{t}{q_{e}} + \frac{1}{k_{2} q 2_{e}}$ ，通过公式变形将t/q_t作为y，t作为x得

$y = \frac{1}{q_{e}} x + \frac{1}{k_{2} q 2_{e}}$ (11)

Weber-Morris颗粒内扩散模型：

$q_{t} = k_{3} t^{\frac{1}{2}} + C$ (12)

其中k₁伪一级动力学模型的速率常数，min⁻¹；k₂伪二级动力学模型的速率常数，(g/(mg∙min))；K₃为Weber-Morris颗粒内扩散模型速率常数，mg/(g∙min^1/2)；t为吸附时间，min；C为Weber-Morris颗粒内扩散模型常数。

图10~12为吸附动力学拟合图，H103树脂对大黄素吸附变化主要有三个阶段。第一阶段(0~50 min)吸附速率较大；第二阶段(50~110 min)该阶段吸附量仍在增加但吸附速率较小；第三阶段(110~350 min)树脂吸附量达到了平衡，实验得到H103树脂对大黄素的平衡吸附量为65.84 mg/g。由表6中各动力学模型参数可知，伪一阶动力学模型实际测得值相当接近，相关系数R² > 0.99说明吸附过程属于物理吸附与吸附热力学得出的结果一致[23]，通过Weber-Morris颗粒内扩散模型图中可以看出方程均不通过原点其吸附速率应该是受薄膜扩散和颗粒内扩散共同作用。

Figure 10. Pseudo-first-order kinetic model

图10. 伪一级动力学模型

Figure 11. Pseudo-second-order kinetic model

图11. 伪二级动力学模型

Figure 12. The Weber-Morris intraparticle diffusion model

图12. Weber-Morris颗粒内扩散模型

Table 6. Dynamic parameters of emodin adsorption by H103 resin

表6. H103树脂对大黄素吸附动力学参数

Kinetic model	拟合方程	R²	动力学参数
Pseudo-first-order	$y = - 0.04 x + 4.15$	0.99	q_e = 63.91
			K₁ = 0.04
Pseudo-second-order	$y = 0.49 x + 0.014$	0.97	q_e = 71.42
			K₂ × 10³ = 0.4
Weber-Morris 颗粒内扩散模型	$y = 9.61 t^{\frac{1}{2}} - 9.74$ (薄膜内扩散)	0.95	K₃ = 9.61
			C = −9.74
	$y = 1.17 t^{\frac{1}{2}} + 44.78$ (颗粒内扩散)	0.88	K₃ = 1.17
			C = 44.78

4. 结论

(1) 通过对比大孔吸附树脂SD300、D101、H103、X-5、SV300的吸附量，其中H103树脂对大黄素的吸附量最高，平衡吸附量为78 mg/g。通过分子对接吸附材料H103树脂主要结构与大黄素分子的芳香环间存在2个π-π键，键长分别为2.9、2.64 Å。计算得到H103树脂孔径是大黄素分子直径的100倍，证明在理论H103树脂对大黄素分子吸附性能较好。

(2) 吸附热力学研究表明，Freundlich方程能更好地描述H103树脂对大黄素的吸附过程，其吸附为优惠吸附，热力学参数ΔG < 0、ΔH > 0、ΔS > 0。结果表明H103树脂对大黄素的吸附过程为自发进行、易发生、亲和力较好的过程，具有明显的物理吸附特征。

(3) 吸附动力学研究结果表明，大黄素在大孔吸附树脂H103上吸附平衡时间为130 min，吸附过程符合伪一级动力学方程R²值大于0.99，吸附速率由颗粒内扩散和膜扩散步骤共同控制。

NOTES

^*通讯作者。

参考文献

[1]	Wang, Y., Zhao, M., Li, B. and Geng, X. (2024) Advances in the Mechanism of Emodin-Induced Hepatotoxicity. Heliyon, 10, e33631. https://doi.org/10.1016/j.heliyon.2024.e33631
[2]	Du, H., Zhu, J., Chen, J., Zhou, H., Yang, J. and Wan, H. (2021) Revealing the Therapeutic Targets and Molecular Mechanisms of Emodin-Treated Coronavirus Disease 2019 via a Systematic Study of Network Pharmacology. Aging, 13, 14571-14589. https://doi.org/10.18632/aging.203098
[3]	Mitra, S., Anjum, J., Muni, M., Das, R., Rauf, A., Islam, F., et al. (2022) Exploring the Journey of Emodin as a Potential Neuroprotective Agent: Novel Therapeutic Insights with Molecular Mechanism of Action. Biomedicine & Pharmacotherapy, 149, Article 112877. https://doi.org/10.1016/j.biopha.2022.112877
[4]	于彦民, 王真真. 大黄素来源的研究进展综述[J]. 科技视界, 2020(16): 199-201.
[5]	Huang, C., Huang, Z., Hu, Y., Li, Z., Wu, Y. and Gao, J. (2021) Treatment of Wastewater Generated from Traditional Chinese Medicine Processing and Utilization: Recent Advances and Future Outlook. Journal of Cleaner Production, 291, Article 125927. https://doi.org/10.1016/j.jclepro.2021.125927
[6]	Tirado-Rives, J. and Jorgensen, W.L. (2008) Performance of B3LYP Density Functional Methods for a Large Set of Organic Molecules. Journal of Chemical Theory and Computation, 4, 297-306. https://doi.org/10.1021/ct700248k
[7]	Masnabadi, N., Thalji, M.R., Alhasan, H.S., Mahmoodi, Z., Soldatov, A.V. and Ali, G.A.M. (2022) Structural, Electronic, Reactivity, and Conformational Features of 2,5,5-Trimethyl-1,3,2-Diheterophosphinane-2-Sulfide, and Its Derivatives: DFT, MEP, and NBO Calculations. Molecules, 27, Article 4011. https://doi.org/10.3390/molecules27134011
[8]	孙小梅, 马聪, 吴士筠, 等. 大黄素的分光光度测定及其应用[J]. 中南民族学院学报(自然科学版), 2001, 20(4): 13-16.
[9]	周昊, 杨瑾, 张玉红. 望江南茎中大黄素提取工艺优化、纯化及鉴定[J]. 东北林业大学学报, 2023, 51(8): 65-70+76.
[10]	李东辉, 王临艳, 吴红伟, 等. 大黄素药理作用研究新进展[J]. 中华中医药学刊, 2021, 39(9): 100-104.
[11]	Mozaffari Majd, M., Kordzadeh-Kermani, V., Ghalandari, V., Askari, A. and Sillanpää, M. (2022) Adsorption Isotherm Models: A Comprehensive and Systematic Review (2010-2020). Science of the Total Environment, 812, Article 151334. https://doi.org/10.1016/j.scitotenv.2021.151334
[12]	Huang, J., Wang, X. and Huang, K. (2009) Adsorption of P-Nitroaniline by Phenolic Hydroxyl Groups Modified Hyper-Cross-Linked Polymeric Adsorbent and XAD-4: A Comparative Study. Chemical Engineering Journal, 155, 722-727. https://doi.org/10.1016/j.cej.2009.09.012
[13]	杨志钧, 蒯翠, 吴锴, 等. 达托霉素在大孔吸附树脂上的吸附性能研究[J]. 化学工程, 2020, 48(4): 11-15.
[14]	Wang, J. and Guo, X. (2020) Adsorption Isotherm Models: Classification, Physical Meaning, Application and Solving Method. Chemosphere, 258, Article 127279. https://doi.org/10.1016/j.chemosphere.2020.127279
[15]	Yang, R., Meng, D., Song, Y., Li, J., Zhang, Y., Hu, X., et al. (2012) Simultaneous Decoloration and Deproteinization of Crude Polysaccharide from Pumpkin Residues by Cross-Linked Polystyrene Macroporous Resin. Journal of Agricultural and Food Chemistry, 60, 8450-8456. https://doi.org/10.1021/jf3031315
[16]	龚剑, 宁方红, 徐环昕, 等. 大孔树脂对纽莫康定B0的吸附行为研究[J]. 离子交换与吸附, 2023, 39(5): 431-439.
[17]	叶佳丰, 马晶菁, 张雨瑞, 等.油茶蒲总黄酮大孔吸附树脂纯化工艺优化及其体外抗氧化、抗炎活性[J/OL]. 食品工业科技, 1-20. https://doi.org/10.13386/j.issn1002-0306.2024040151, 2024-11-13.
[18]	Wang, J. and Guo, X. (2020) Adsorption Isotherm Models: Classification, Physical Meaning, Application and Solving Method. Chemosphere, 258, Article 127279. https://doi.org/10.1016/j.chemosphere.2020.127279
[19]	刘毅, 刘永峰, 裴栋, 等. 大孔吸附树脂对驼血多肽的吸附行为研究[J]. 离子交换与吸附, 2016, 32(5): 385-393.
[20]	Chen, Y. and Zhang, D. (2014) Adsorption Kinetics, Isotherm and Thermodynamics Studies of Flavones from Vaccinium Bracteatum Thunb Leaves on NKA-2 Resin. Chemical Engineering Journal, 254, 579-585. https://doi.org/10.1016/j.cej.2014.05.120
[21]	Gökmen, V. and Serpen, A. (2002) Equilibrium and Kinetic Studies on the Adsorption of Dark Colored Compounds from Apple Juice Using Adsorbent Resin. Journal of Food Engineering, 53, 221-227. https://doi.org/10.1016/s0260-8774(01)00160-1
[22]	Rodrigues, L.A., Maschio, L.J., da Silva, R.E. and da Silva, M.L.C.P. (2010) Adsorption of Cr(vi) from Aqueous Solution by Hydrous Zirconium Oxide. Journal of Hazardous Materials, 173, 630-636. https://doi.org/10.1016/j.jhazmat.2009.08.131
[23]	Lanjekar, K.J. and Rathod, V.K. (2023) Recovery and Separation of Glycyrrhizic Acid from Natural Deep Eutectic Solvent (NADES) Extract by Macroporous Resin: Adsorption Kinetics and Isotherm Studies. Preparative Biochemistry & Biotechnology, 54, 39-48. https://doi.org/10.1080/10826068.2023.2204485

为你推荐

友情链接