Anti-hepatic injury effect of enzymatic hydrolysate from soft tissue of Pinctada martensii
-
摘要:
为探究马氏珠母贝肉酶解产物 (Enzymatic hydrolysate from Pinctada martensii, EP)对酒精性肝损伤 (Alcoholic liver damage, ALD)的保护作用,该研究将EP超滤分级为截留分子量>10 kD (EP-Ⅰ)、3~10 kD (EP-Ⅱ)和<3 kD (EP-Ⅲ) 3个组分,检测其体外抗肝损伤活性及对ALD小鼠肝保护作用的影响。体外试验结果显示,EP-Ⅲ可显著激活体外乙醇脱氢酶 (ADH)活性 (P<0.01),3个超滤组分均具有一定的体外抗氧化能力且ep-ⅲ>EP-Ⅱ>EP-Ⅰ;动物试验结果显示,与模型对照组相比,各超滤组分均能够显著降低小鼠血清中谷丙转氨酶 (ALT)和谷草转氨酶 (AST)活力、小鼠肝脏指数及肝脏中丙二醛 (MDA)和甘油三酯 (TG)含量,同时显著增强小鼠肝脏中超氧化物歧化酶 (SOD)、乙醇脱氢酶 (ADH)和乙醛脱氢酶 (ALDH)活力,提高肝脏中的谷胱甘肽 (GSH)含量。综上,马氏珠母贝肉酶解超滤组分对急性ALD具有一定的辅助保护作用,其中EP-Ⅲ的保护作用效果最佳,其机制可能与加快机体乙醇代谢和减缓乙醇对机体造成的氧化损伤相关。
-
关键词:
- 马氏珠母贝肉酶解产物 /
- 酒精性肝损伤 /
- 体外抗肝损伤活性 /
- 乙醇代谢
Abstract:In order to analyze the anti-hepatic injury effect of enzymatic hydrolysate from Pinctada martensii (EP), we seperated EP by membrane ultrafiltration into three molecular size fractions [MW> 10 kD (EP-Ⅰ), MW=3–10 kD (EP-Ⅱ) and MW < 3 kD (EP-Ⅲ)], and then measured the effects of three ultrafiltration fractions on the anti-ALD in vitro and liver-protection on mice. The results of experiments in vitro show that EP-III could activate alcohol dehydrogenase (ADH) significantly (P<0.01), and="" three="" fractions="" demonstrated="" different="" antioxidant="" capacities="">EP-Ⅱ>EP-Ⅰ). The results of mice experiment show that compared with the alcohol model group, the activities of ALT and AST in serum, liver index, the levels of MDA and TG in liver decreased in each ultrafiltration fraction group significantly, while the activities of SOD, ADH and ALDH and the levels of GSH in liver increased significantly. Therefore, the ultrafiltration fractions of enzymatic hydrolysate from soft tissue of P. martensii showed good protective effect on alcoholic liver damage, and the effect of EP-Ⅲ with low molecular mass was the best. Its mechanism may be related to accelerating ethanol metabolism and slowing down the oxidative damage caused by ethanol.
-
随着水产养殖业的迅速发展和养殖规模的不断扩大,精准预测和监测水质参数对于保障养殖效益[1]和生物健康变得至关重要。然而,现今水产养殖业面临着许多挑战,维持良好的水质环境便是其中之一。在养殖系统中,溶解氧[2]是一个非常重要的水质参数,溶解氧不足可能导致鱼类窒息、生长受限甚至死亡。传统的溶解氧测量方法需要采集水样在实验室进行分析,这种方法耗时久、成本高,且无法提供实时监测数据。随着养殖方式向深远海[3]养殖工船发展,对养殖水质的要求日趋严格,溶解氧的精准预测变得越来越重要。近年来,大数据和机器学习技术[4]为溶解氧预测提供了新的解决方案,可以利用历史监测数据和水质参数建立预测模型,准确预测溶解氧浓度。这些模型可以基于气象、水质、养殖方式和其他相关因素进行训练,能及时响应环境变化和养殖条件调整,以维持良好的水质环境,从而促进鱼类健康生长。
精准预测溶解氧对水产养殖意义重大。近年来,国内外学者对水质预测进行了广泛研究,并对提高预测精度的方法作了改进。沈时宇和陈明[5]采用一种Prophet时序模型的方法来提高水质预测精度,但该方法存在模型训练欠拟合缺陷。Ren等[6]选择基于深度信念网络的溶解氧预测模型,成功实现了溶解氧含量预测。Cao等[7]提出了基于K均值聚类和GRU (Gate Recurrent Unit) 神经网络的池塘养殖溶解氧预测模型,该方法虽然提高了预测精度,但GRU神经网络仍无法完全解决梯度消失问题,且随着数据量和模型规模的增加,效果不理想。此外,研究者还探索了其他方法。宦娟等[8]提出了一种基于集合经验模态分解、游程检测法重构、适宜的单项预测算法建模和BP神经网络非线性叠加的组合预测模型,但在高频分量波动较大的情况下,该算法参数难以确定,预测效果不佳。林彬彬等[9]提出基于麻雀搜索算法 (Sparrow Search Algorithm, SSA)[10-11]和长短期记忆神经网络 (Long Short-Term Memory, LSTM) [12]的黄鳝 (Monopterus albus) 池塘溶解氧浓度预测模型,该模型具有良好的准确性和鲁棒性,但仍存在局部最优解问题,导致模型参数的准确性不高,进而降低了预测的准确性。
当考虑到模型参数的后期优化时,通常会使用算法代替人工参数设置。为解决这一问题,本研究采用了两种启发式算法—麻雀搜索算法和海洋捕食算法 (Marine Predators Algorithm, MPA)[13],以及两种群体智能算法—灰狼优化算法 (Grey Wolf Optimizer, GWO) [14]和鲸鱼算法 (Whale Optimization Algorithm, WOA)[15]。这些算法在全局搜索和多模态问题的处理上均表现出色,并且在不同领域有成功的应用。基于这些算法的特点及其在应用领域的成功案例[16-19],本文比较了这些算法在LSTM模型中进行参数优化的效果,并通过优化模型参数实现既定目标。
为了提高溶解氧参数预测的准确性,本文提出了一种增强型麻雀搜索算法 (Enhance Sparrow Search Algorithm, ESSA) 用于优化LSTM模型参数,更好地调整模型的结构和初始状态,以满足溶解氧参数的预测需求,从而提高预测结果的精确度。
1. 材料与方法
1.1 数据来源
本文用于验证模型的数据来源于某养殖工船,其相关信息为:全长249.9 m,型宽14 m,型深21.5 m,排水量130 000 t,航速10.6 kn,养殖水体90 000 m3,单舱水体6 000 m3,定员3人,年产量370 t,总投资4.5亿元。该工船在福建宁德海域养殖大黄鱼 (Larimichthys crocea),传感器布置于鱼舱内深度14 m处,采集数据开始时间为2023-04-11 00:05:00,结束时间为2023-04-24 00:00:00。采集数据的步长为5 min。一共有3 745组数据。选取MPS-400多参数传感器进行数据采集,该传感器具有自清洗、多参数测量、信号输出协议为Modbus特点,安装传感器的方式为投入式。采集pH、温度 (℃)、盐度 (‰)、溶解氧 (mg·L−1) 4类参数。本文中水质预测目标为溶氧,水质特征变量为pH、温度、盐度。表1为养殖水质参数溶解氧的原始数据,图1、图2分别为养殖工船场景和原始数据溶解氧的变化曲线。
表 1 养殖水质原始数据表Table 1 Raw data sheet for aquaculture water quality时间
TimepH 盐度
Salinity/‰温度
Temperature/℃溶解氧
DO/(mg·L−1)2023-04-11 00:05:00 7.86 36.099 998 20.114 88 9.119 555 2023-04-11 01:05:00 7.84 35.900 002 20.127 55 9.047 465 2023-04-11 02:05:00 7.84 36.200 001 20.123 21 8.998 263 ··· ··· ··· ··· ··· 2023-04-23 22:05:00 7.69 34.200 001 22.060 74 6.736 285 2023-04-23 23:05:00 7.67 34.200 001 22.063 86 6.419 692 2023-04-24 00:00:00 7.67 34.299 999 22.083 02 6.358 287 1.2 研究方法
1.2.1 数据处理
1) 数据预处理。采集的数据,可能会因数据传输异常等导致丢失,需要对数据作预处理,对缺失的数值需要通过线性插值进行填补,公式如下:
$$ {x}_{k}{\text{=}}{x}_{f}{\text{+}}({x}_{b}-{x}_{f})(k-f)/(b-k) $$ (1) 式中:
$ {x}_{k} $ 表示需要填补的值;$ {x}_{f} $ 表示前一个非缺失数据点的值;$ {x}_{b} $ 表示后一个非缺失数据点的值;f 和$ b $ 分别表示前后非缺失数据点的索引。数据归一化[20]是将不同特征的数据按照一定规则进行缩放,以消除特征之间的量纲差异,使得不同特征具有相同的尺度。本文使用Min-Max归一化的具体表达式如下:
$$ \begin{array}{c} {x}_{{\rm{nor}}}{\text{=}}\dfrac{{x}_{i}{\text{−}}{x}_{{\rm{min}}}}{{x}_{{\rm{max}}}{\text{−}}{x}_{{\rm{min}}}} \end{array} $$ (2) 式中:
$ {x}_{\mathrm{m}\mathrm{a}\mathrm{x}} $ 表示样本数据最大值;$ {x}_{\mathrm{m}\mathrm{i}\mathrm{n}} $ 表示样本数据最小值;$ {x}_{i} $ 表示样本数据中的第$ i $ 个数据;$ {x}_{\mathrm{n}\mathrm{o}\mathrm{r}} $ 表示归一化后的数据。在经过归一化的训练和测试后,为了得到真实的预测结果,需要对预测值进行反归一化,将其重新映射回原始数据的范围。2) 数据集的划分。首先将采集到的数据划分为训练集和测试集,本实验是将前12 d的数据作为训练集,最后1 d的数据作为测试集。使用训练集来训练模型,使用测试集评估模型在未见过的数据上的性能。
1.2.2 长短时记忆网络
为了有效捕捉具有长期依赖关系的序列数据,如连续时间上测得的溶解氧参数,传统的神经网络模型往往无法提供令人满意的结果。为解决该问题,长短时记忆网络被引入并广泛应用。LSTM是一种特殊的循环神经网络 (Recurrent Neural Network, RNN)[21],专门设计用于处理序列数据。LSTM与其他神经网络模型之间的区别在于独特的内部单元结构 (图3)。
图3中,f、
$ i $ 、$ {o} $ 分别表示遗忘门、输入门、输出门;$ \sigma $ 表示全连接层。$ {x}_{t} $ 是t时刻的输入,$ {c}_{t-1} $ 和$ {c}_{t} $ 分别表示上一时刻和这一时刻的输出,$ {\tilde{c}}_{t} $ 为t时刻候选记忆细胞的信息,$ {h}_{t-1}{\mathrm{、}h}_{t} $ 分别表示为上一时刻和这一时刻隐藏状态的信息,tanh为激活函数。3个门的更新计算公式为:
$$ f\left(x\right){\text{=}}\sigma (Wx+b) , \sigma \left(x\right){\text{=}}\frac{1}{1{\text{+}}{\mathrm{e}}^{{\text{−}}x}} $$ (3) 式中:
$ W$ 和 b 分别表示网络的权重矩阵和偏置向量;$ \sigma \left(x\right) $ 是Sigmoid激励函数,为平滑的阶梯函数。它可以将任何值转化为0~1,此函数的输出介于0~1,代表具体有多少信息能够流过Sigmoid层,0表示不能通过,1表示能全部通过。1.2.3 ESSA-LSTM预测模型
水质参数溶解氧受多种因素的影响,由于数据采集于养殖工船,其空间特征是一定的,因此本文仅考虑水质中的pH、温度、盐度对溶解氧的影响。为了更好地预测溶解氧,引入长短时记忆网络,并针对该网络模型中存在的参数问题引入增强型麻雀搜索算法进行改进。在基本的LSTM模型中,权重和偏置是通过矩阵乘法和加法来更新的,而在本文使用的LSTM模型中,参数的更新是通过训练过程中的反向传播和优化算法来实现。增强型麻雀搜索算法优化LSTM模型结构如图4所示,模型流程如图5所示。
在LSTM模型中,隐藏层神经元个数是影响模型性能关键性的超参数。本文在LSTM模型的训练模型中引入了Adam (Adaptive Moment Estimation) 梯度下降[22]算法,为此,该算法的初始化学习速率也成为影响模型性能的关键。通过采用增强型麻雀搜索算法,来确定隐藏层的神经元个数和初始化学习速率。具体步骤如下:
1) 初始化参数,将数据进行归一化处理,并且将数据分为训练集和测试集。
2) 确定LSTM模型的拓扑结构。
3) 通过增强型麻雀搜索算法对LSTM模型参数进行优化。
4) 利用优化后的参数建立LSTM模型和训练神经网络模型。
5) 将训练集和测试集数据来对建立好的模型进行实验和验证,并通过模型评价指标来评价模型的预测效果。
1.2.4 模型评价标准
通过平均绝对误差 (Mean Absolute Error, MAE)、平均绝对百分比误差 (Mean Absolute Percentage Error, MAPE)、均方根绝对误差 (Root Mean Square Error, RMSE)和决定系数 (Coefficient of Determination,又称R2值) 对各个模型的预测效果进行评判其表达式如式 (4)、 (5)、 (6) 和 (7) 所示:
$$ \mathrm{M}\mathrm{A}\mathrm{E}{\text{=}}\frac{\sum _{i{\text{=}}1}^{n}\left|{y}_{i}{\text{−}}{x}_{i}\right|}{n} $$ (4) $$ \mathrm{M}\mathrm{A}\mathrm{P}\mathrm{E}{\text{=}}\frac{1}{n}\sum _{i{\text{=}}1}^{n}\left|\frac{{y}_{i}{\text{−}}{x}_{i}}{{y}_{i}}\right| \times 100\mathrm{\%}$$ (5) $$ \mathrm{R}\mathrm{M}\mathrm{S}\mathrm{E}{\text{=}}\sqrt{\frac{\sum _{i{\text{=}}1}^{n}{\left({y}_{i}{\text{−}}{x}_{i}\right)}^{2}}{n}} $$ (6) $$ {R}^{2}{\text{=}}1-\frac{\sum _{i{\text{=}}1}^{n}{({y}_{i}{\text{−}}{x}_{i})}^{2}}{\sum _{i{\text{=}}1}^{n}{({y}_{i}{\text{−}}\overline{y})}^{2}} $$ (7) 式中:
$ {y}_{{i}} $ 表示实际值;$ {x}_{{i}} $ 表示预测值;$ \overline{y} $ 表示实际值的平均值;$ n $ 表示样本个数。当MAE、MAPE和RMSE越小时,所建模型的预测效果越好;$ {R}^{2} $ 值越接近1,表明模型对数据的拟合程度越高。1.2.5 基本麻雀搜索算法
在麻雀[23-26]捕食过程中,分为侦察者和跟踪者两个群体。侦察者负责在整个区域进行觅食和引导跟踪者觅食的方向,而跟踪者负责追随侦察者进行觅食。当危险来临时,麻雀种群会做出反捕食行为。
在每次迭代的过程中,侦察者的位置更新如下:
$$ {X}_{i,j}^{t+1}{\text{=}}\left\{\begin{array}{l}{\mathrm{X}}_{i,j}^{t}\cdot {\rm{exp}}\left(\frac{{\text{−}}i}{\alpha \cdot {{\rm{iter}}}_{{\rm{max}}}}\right)\quad {\rm{if}}\;\;{R}_{2} < {\rm{ST}}\\ {\mathrm{X}}_{i,j}^{t}{\text{+}}Q\cdot L\quad \quad \quad \quad {\rm{if}}\;\;{R}_{2}\geqslant {\rm{ST}}\end{array}\right. $$ (8) 式中:
$ t $ 表示迭代的次数;$ {\mathrm{X}}_{i,j}^{t} $ 表示迭代 t 次时第$ i $ 只麻雀的第$ j $ 维数的值;$ {\mathrm{X}}_{i,j}^{t+1} $ 表示迭代$ t+1 $ 次时第$ i $ 只麻雀的第$ j $ 维数的值;$ {{\rm{iter}}}_{\mathrm{m}\mathrm{a}\mathrm{x}} $ 表示迭代次数最多的一个常数;$ \alpha $ 为介于0~1的一个随机数;$ {R}_{2} $ 表示报警值,范围在$ 0 $ ~1;ST表示安全阈值,范围在0.5~1;Q是一个服从正态分布的随机数;L表示1$ \times d $ 的矩阵,其中,d矩阵的元素全为1。跟踪者的位置更新如下:
$$ {X}_{i,j}^{t+1}{\text{=}}\left\{\begin{array}{l}Q\cdot {\rm{exp}}\left(\frac{{X}_{{\rm{worst}}}^{t}{\text{−}}{X}_{i,j}^{t}}{{i}^{2}}\right)\qquad \quad \quad {\rm{if}}\;\;i > n/2\\ {X}_{P}^{t{\text{+}}1}{\text{+}}\mid {X}_{i,j}^{t}{\text{−}}{X}_{P}^{t{\text{+}}1}\mid {A}^{\text{+}}L\quad \quad {\rm{otherwise}}\end{array}\right. $$ (9) 式中:
$ Q $ 是一个服从正态分布的随机数;$ {X}_{p} $ 表示侦察者处于一个最佳的位置;$ {X}_{P}^{t+1} $ 是在迭代$ t+1 $ 时点上的全局最优解;而$ {X}_{{\rm{worst}}}^{t} $ 表示当前处于全局最差的位置;A是一个1$ \times d $ 的矩阵,里面的元素随机分配为1或−1,并且满足$ {A}^{+}{\text{=}}{A}^{T}(A{A}^{T}{)}^{-1} $ ;$ n $ 代表麻雀的数量,当$ i > n/2 $ 表示适应值比较差的第i个跟踪者觅食能力低,捕获的食物比较少;$ \mathrm{o}\mathrm{t}\mathrm{h}\mathrm{e}\mathrm{r}\mathrm{w}\mathrm{i}\mathrm{s}\mathrm{e} $ 代表其他情况。当危险来临时,麻雀会做出预警行为,其预警者位置更新如下:
$$ {X}_{i,j}^{t+1}{\text{=}}\left\{\begin{array}{l}{X}_{\mathrm{b}\mathrm{e}\mathrm{s}\mathrm{t}}^{t}{\text{+}}\beta \mid {X}_{i,j}^{t}{\text{−}}{X}_{\mathrm{b}\mathrm{e}\mathrm{s}\mathrm{t}}^{t}\mid \quad {\rm{if}}\;\; {f}_{i} > {f}_{g}\\ {X}_{i,j}^{t}{\text{+}}K\cdot \left(\dfrac{\left|{X}_{i,j}^{t}{\text{−}}{X}_{{\rm{wotst}}}^{t}\right|}{({f}_{i}{\text{−}}{f}_{w}){\text{+}}\epsilon }\right)\quad {\rm{if}}\;\;{f}_{i}\;{\text{=}}\;{f}_{g}\end{array}\right. $$ (10) 式中:
$ \beta $ 表示步长控制参数,是一个均值为0、方差为1的正态随机分布;$ K\mathrm{是}\mathrm{一}\mathrm{个} $ 随机数,其范围是 −1~1,表示麻雀移动方向的步长控制系数;$ {f}_{i} $ 表示当前麻雀的适应值;$ {f}_{g} $ 则是当前全局的最佳适应值;$ {f}_{w} $ 表示当前全局的最差适应值;$ \epsilon $ 是一个极小的常数,避免分母出现0的情况。1.2.6 增强型麻雀搜索算法
1.2.6.1 Circle混沌映射初始化种群
混沌映射[27]是指在非线性系统中产生混沌行为的一种数学映射,常见的混沌映射方法有Logistic混沌映射、Henon混沌映射、Tent混沌映射和Circle混沌映射,不同的混沌映射有其自身的特点,由于Circle混沌[23]算法易于实现,且具有周期性和可控性,其表达式如下:
$$ {x}_{i+1}{\text{=}}{\rm{mod}}\left[{x}_{i}{\text{+}}0.2{\text{−}}\frac{0.5}{2\pi }\mathrm{s}\mathrm{i}\mathrm{n}\left(2\pi {x}_{i}\right),1\right] $$ (11) 式中:
$ {\rm{mod}} $ 代表取余;$ {x}_{i} $ 表示当前迭代的值;$ {x}_{i+1} $ 表示下一个步骤迭代的值。基本的SSA算法是采用随机生成的方式来初始化种群,这就导致了分布不均匀的问题,为此,可以使用Circle混沌映射的方式来初始化种群。1.2.6.2 正弦余弦算法对侦察者位置优化
麻雀搜索算法中侦察者位置更新主要依赖于个体之间的位置调整和信息共享,这也就导致了侦察者无法探索整个搜索空间,难以跳出局部最优解。针对该问题,引入了正弦余弦算法[28]对麻雀搜索算法中的侦察者位置进行更新,在该算法中,加入了非线性正弦学习因子如公式 (12) 所示,优化后的侦察者位置如公式 (13) 所示:
$$ \omega {\text{=}}{\omega }_{\mathrm{m}\mathrm{i}\mathrm{n}}{\text{+}}({\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}}{\text{−}}{\omega }_{\mathrm{m}\mathrm{i}\mathrm{n}})\cdot \mathrm{s}\mathrm{i}\mathrm{n}({t}\pi /{{\rm{iter}}}_{\mathrm{m}\mathrm{a}\mathrm{x}}) $$ (12) $$ \begin{array}{c} {X}_{i,j}^{t+1}{\text{=}}\\ \left\{\begin{array}{ll}(1{\text{−}}\omega )\cdot {X}_{i,j}^{t}{\text{+}}\omega \cdot \mathrm{s}\mathrm{i}\mathrm{n}\left({r}_{1}\right)\cdot \mid {r}_{2}\cdot {X}_{{\rm{best}}}{\text{−}}{X}_{i,j}^{t}\mid ,& {R}_{2} < {\rm{ST}}\\ (1{\text{−}}\omega )\cdot {X}_{i,j}^{t}{\text{+}}\omega \cdot \mathrm{c}\mathrm{o}\mathrm{s}\left({r}_{1}\right)\cdot \mid {r}_{2}\cdot {X}_{{\rm{best}}}{\text{−}}{X}_{i,j}^{t}\mid , & {R}_{2}\geqslant {\rm{ST}}\end{array} \right. \end{array}$$ (13) 式中:
$ {\omega }_{\mathrm{m}\mathrm{i}\mathrm{n}}$ 和${\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}} $ 表示惯性权重的最小和最大值;$ \omega $ 表示适应权重值;$ {{\rm{iter}}}_{\mathrm{m}\mathrm{a}\mathrm{x}} $ 表示迭代次数最多的一个常数;$ {r}_{1} $ 为 0~2$ \pi $ 内的随机数;$ {r}_{2} $ 是 0~2 内的随机数;$ {X}_{\mathrm{b}\mathrm{e}\mathrm{s}\mathrm{t}}\mathrm{表}\mathrm{示}\mathrm{群}\mathrm{体}\mathrm{中}\mathrm{最}\mathrm{佳}\mathrm{的}\mathrm{位}\mathrm{置} $ ;$ {\mathrm{X}}_{i,j}^{t} $ 表示迭代 t 次时第 i 个麻雀的第 j 维数的值。1.2.6.3 Levy飞行策略对跟踪者位置优化
由于在麻雀搜索算法中,跟踪者的位置更新主要依赖于当前最佳位置和全局最佳位置,使得算法在搜索空间中的多样性不足,从而限制了算法的全局搜索能力。并且当全局最佳位置附近存在较多的最优解时,算法容易陷入局部最优解而无法跳出。为此,本文提出通过引入Levy策略[29-30]对跟踪者的位置进行更新,Levy[31]步长向量生成公式如下:
$$ S{\text{=}}\frac{u}{{\mid v\mid }^{\frac{1}{\beta }}} $$ (14) $$ \sigma {\text{=}}{\left\{\frac{\mathrm{\Gamma }(1{\text{+}}\beta )\mathrm{s}\mathrm{i}\mathrm{n}\left(\dfrac{\pi \beta }{2}\right)}{\beta \mathrm{\Gamma }\left(\dfrac{1{\text{+}}\beta }{2}\right){2}^{\frac{\beta {\text{−}}1}{2}}}\right\}}^{\frac{1}{\beta }} $$ (15) 式中:S为飞行路径;
$ \mathrm{\Gamma } $ 为 gamma函数;$ \beta $ 取值为1.5;$ \sigma $ 代表正态分布中的方差;u 和 v 是符合正态分布的随机数,$ u~N(0,{\sigma }^{2}) $ ,$ v~N\left(\mathrm{0,1}\right) $ ,则跟踪者位置更新公式如下:$$ {X}_{i,j}^{t{\text{+}}1}{\text{=}}\left\{\begin{array}{l}Q\cdot {\rm{exp}}\left(\dfrac{{X}_{{\rm{worst}}}^{t}-{X}_{i,j}^{t}}{{i}^{2}}\right)\qquad {\rm{if}}\;\;i > n/2\\ {X}_{P}^{t{\text{+}}1}{\text{+}}{X}_{P}^{t{\text{+}}1}\cdot (0.01\cdot S) \quad {\rm{otherwise}}\end{array}\right. $$ (16) 增强型麻雀搜索算法具体流程如下:
1) 设置参数:设置最大的迭代次数,麻雀的种群数量,每个维度上的变量最小值,每个维度上的变量最大值。
2) 初始化麻雀种群位置、优化全局最优个体和最优适应度,计算每个个体的适应度并且对其进行排序。
3) 在每次迭代的过程中,利用公式 (13) 更新侦察者的位置,实现对全局的搜索,并且进行边界检查。
4) 使用公式 (16) 更新跟踪者位置,使得跟踪者不断靠近侦察者位置。
5) 如果预警者发现危险,则通过公式远离危险位置。
6) 在麻雀搜索食物之后,对麻雀个体的最优适应度值进行排序。
7) 更新全局最优解和最佳适应度。
8) 进入循环迭代,转至步骤3)。
9) 判断是否符合中止条件,若符合,则输出全局最优解和最佳适应度。若不符合,回到步骤 3) 重新开始计算。
2. 实验与分析
2.1 基准函数优化性能测试
为了验证增强型麻雀搜索算法在求解优化过程中的可行性和优越性,将ESSA与SSA、GWO、MPA、WOA算法在5种基准函数上进行测试对比。
使用Matlab 2019b软件对5种不同的算法进行对比实验,为了避免实验出现的偶然性造成误差,并增加实验的准确性,选取5种基准测试函数分别独立运行30次,取平均值和标准差作为实验的评价标准。实验采用的基准测试函数如表2所示,F1和F2为高维单峰函数,F3和F4为高维多峰函数,F5为低维多峰函数。种群大小为40个,最大迭代次数为100次。实验结果如表3所示,优化测试函数收敛曲线对比如图6所示。
表 2 基准测试函数Table 2 Process of benchmarking基准测试函数
Benchmark function维度
Dimensionality搜索空间
Search space最优值
Optimal value$ {F}_{1}\left(x\right){\text{=}}{\sum }_{i{\text{=}}1}^{n}\mid {x}_{i}\mid {\text{+}}{\prod }_{i{\text{=}}1}^{n}\mid {x}_{i}\mid $ 30 $ [{\text{−}}\mathrm{10,\;10}] $ 0 $ {F}_{2}\left(x\right){\text{=}}{\sum }_{i{\text{=}}1}^{n}i{x}_{i}^{4}{\text{+}}\mathrm{r}\mathrm{a}\mathrm{n}\mathrm{d}\mathrm{o}\mathrm{m}\left[\mathrm{0,\;1}\right) $ 30 $ [{\text{−}}\mathrm{1.28,\;1.28}] $ 0 $ {F}_{3}\left(x\right){\text{=}}{\sum }_{i{\text{=}}1}^{n}{\text{−}}{x}_{i}\mathrm{s}\mathrm{i}\mathrm{n}\left(\sqrt{\mid {x}_{i}\mid }\right) $ 30 $ [{\text{−}}\mathrm{500,\;500}] $ $ {\text{−}}418.9829\times n $ ${F}_{4}\left(x\right){\text{=}}\dfrac{1}{4\;000}{\sum }_{i{\text{=}}1}^{n}{x}_{i}^{2}{\text{−}}{\prod }_{i{\text{=}}1}^{n}\mathrm{c}\mathrm{o}\mathrm{s}\dfrac{ {x}_{i} }{\sqrt{i} }{\text{+}}1$ 30 $ [{\text{−}}\mathrm{600,\;600}] $ 0 $ {F}_{5}\left(x\right){\text{=}}-\sum _{i{\text{=}}1}^{4}{c}_{i}\mathrm{e}\mathrm{x}\mathrm{p}\left[-\sum _{j{\text{=}}1}^{6}{a}_{ij}({x}_{i}-{p}_{ij}{)}^{2}\right] $ 30 $ \left[\mathrm{0,\;1}\right] $ $ -3.32 $ 表 3 基准函数优化结果对比Table 3 Comparison of benchmark function optimization results算法种类
Type of algorithm$ {F}_{1} $ $ {F}_{2} $ $ {F}_{3} $ $ {F}_{4} $ $ {F}_{5} $ ESSA 最优值 Ave $ 2.351\,3{\times 10}^{{\text{−}}16} $ $ 4.306\,9{\times 10}^{{\text{−}}4} $ $ {\text{−}}1.030\,1{\times 10}^{4} $ 0 −3.266 5 标准差 Std $ 1.242\,5{\times 10}^{{\text{−}}15} $ $ 5.081\,4{\times 10}^{{\text{−}}4} $ $ 1.842\,8{\times 10}^{3} $ 0 0.059 3 SSA 最优值 Ave 0.057 3 0.003 0 $ {\text{−}}3.682\,2{\times 10}^{3} $ $ 4.899\,3{\times 10}^{{\text{−}}5} $ −3.138 5 标准差 Std 0.063 1 0.002 7 359.466 6 $ 9.669\,2{\times 10}^{{\text{−}}5} $ 0.081 4 GWO 最优值 Ave 0.025 5 0.015 9 ${\text{−} }5.646\,0{\times 10}^{ 3}$ 0.082 2 −3.238 5 标准差 Std 0.013 1 0.007 8 879.281 3 0.059 1 −3.322 0 MPA 最优值 Ave 0.091 3 0.006 8 $ {\text{−}}7.235\,6{\times 10}^{3} $ 0.290 3 0.075 0 标准差 Std 0.031 1 0.003 8 685.562 5 0.128 3 1.750 6$ {\times 10}^{{\text{−}}4} $ WOA 最优值 Ave $ 7.826\,5{\times 10}^{{\text{−}}9} $ 0.019 4 $ {\text{−}}8.801\,8{\times 10}^{3} $ 0.086 9 −3.104 0 标准差 Std $ 1.722\,1{\times 10}^{{\text{−}}8} $ 0.017 9 $ 1.354\,2{\times 10}^{3} $ 0.262 4 0.340 7 根据表2中各种算法的最优值 (Ave) 和标准差 (Std) ,可以看出ESSA算法尽管在低维多峰函数上的性能差异和其他算法相差无几,但是在高维单峰、高维多峰函数上的性能均优于其他算法。由图6也可以明显看出,ESSA算法在收敛性和优化结果方面比其他4种算法的表现更佳。针对F1和F2高维单峰函数的收敛情况,不论迭代次数增加到何种程度,使用ESSA算法获得的最优值始终优于其他算法。此外,在经过100次迭代后,ESSA算法获得的最优值数量级明显高于其他算法。对于F3和F4高维多峰函数而言,ESSA算法的最优值优于WOA和MPA算法,且最优值几乎是后两种算法的2倍。这表明ESSA算法在处理高维多峰函数时优势明显。在F5低维多峰函数中,虽然ESSA算法的最优值相对于其他算法没有明显的优越性,但仍取得了显著的改进效果。
综上,可以得出ESSA算法在收敛结果和优化效果方面比其他算法表现出色,为进一步研究和应用ESSA算法提供了有力的支持。
2.2 预测准确率优化实验
2.2.1 实验设置
为了验证ESSA算法对LSTM模型参数改进的效果最佳,选取基本的SSA算法、GWO算法、MPA算法、WOA算法分别对LSTM模型改进进行对比。为了保证实验的准确性和算法对比结果的可信度,实验均采用Matlab 2019b软件进行,几种算法的实验参数一致,且每种算法的仿真实验独立运行10次然后取平均值。在实验过程中,每种算法中寻找最优解隐藏层的神经个数范围均设置为 0~100,初始速率范围设置为 0~0.2。
2.2.2 优化对比
表4是实验设置的主要参数,各算法独立运行10次优化后的平均预测率结果对比如表5所示。各算法实验得出的平均绝对误差 、均方根误差、平均绝对百分比误差和决定系数的评价指标如图7所示。根据实验所得的平均绝对百分比误差,得到ESSA、SSA、GWO、MPA和WOA算法优化的预测率收敛曲线如图8所示。
表 4 参数设置Table 4 Settings of parameters参数
Parameter取值
Value训练样本个数 Number of training samples 3 457 测试样本个数 Number of test samples 288 训练迭代次数 Number of training iterations 100 批处理大小个数 Number of batch sizes 32 种群大小 Population size 30 表 5 预测率结果对比Table 5 Comparison of predicted rate算法
Algorithm平均预测率 (10次)
Average rate of prediction (10 times)增强型麻雀搜索算法 ESSA 99.071% 麻雀搜索算法 SSA 96.929% 灰狼优化算法 GWO 92.418% 海洋捕食算法 MPA 92.389% 鲸鱼优化算法 WOA 91.357% 从表5可以直观的看出,在同样的实验环境和参数下,ESSA算法对于优化改进LSTM模型参数得出来的预测率具有明显的效果。ESSA算法优化改进效果略好于基本的SSA算法,但是由ESSA算法改进得出的预测率准确性远远大于GWO、MPA和WOA算法所改进的。
从图7可以清晰地看出,ESSA算法在MAE、RMSE和MAPE指标方面均表现出优于其他4种算法的结果。此外,虽然SSA拟合的效果所对应的
$ {R}^{2} $ 值也接近1,但是稍逊于ESSA,这进一步说明了ESSA算法在预测准确率方面表现出卓越的效果,即ESSA算法具有最佳的预测性能。由图8可以看出,ESSA算法在初始迭代阶段就避免了陷入局部最优解的问题,因为在初始迭代时预测准确率显著提高,且在大约10次迭代后便达到最优解。相比之下,SSA算法在初始迭代时的预测准确率虽然远高于MPA算法,但在第80次迭代时已经陷入了局部最优解。从GWO算法的结果可以看出,随着迭代次数的增加,预测准确率趋于稳定,约70次迭代后达到稳定效果。然而,与ESSA和SSA算法相比,GWO算法的最终结果并不理想。WOA算法在初始迭代阶段的预测效果较差,但随着迭代次数的增加也趋于稳定,明显可以观察到该算法的最终效果低于其他算法。
通过ESSA算法对LSTM模型参数寻优,可以得到隐藏层的神经个数为17,初始学习速率为0.008 6。多种算法优化曲线对比如图9-a和图9-b所示。
根据各种算法的仿真实验结果,可以看出ESSA算法具有收敛速度快、强大的跳出局部最优解的能力等特点。通过对比各种算法对LSTM模型优化后的曲线,可以发现ESSA算法优于其他算法,因此,将ESSA算法用于LSTM模型参数改进,不仅能很好地获得具体的参数值,而且能实现最佳的优化效果。这种方法避免了在模型中盲目设定参数而耗费大量时间和精力,并且当模型中特征向量的维度和数据量发生变化时,运用该算法对模型进行优化,可以迅速、准确地获取具体的模型参数,无需根据人工经验进行多次实验得出参数值再进行比较。将ESSA算法与LSTM模型相结合,可充分发挥两者的优势,实现更精确的预测。这种方法不仅能提高预测的准确率,还减少了参数调整的复杂性,使得模型的优化过程更加高效。对溶氧预测优化研究,不仅适用于养殖工船场景,也适用于如预测某水域某个水质参数的其他场景,可为其他领域的相关研究提供参考。
3. 结论
针对养殖工船水质参数溶解氧预测及其精确性,本文提出了一种利用ESSA算法对LSTM模型参数作优化改进的方法,以实现高效的预测效果。首先,该算法在基本的麻雀算法中加入了Circle混沌映射进行种群初始化,加入正弦余弦算法和Levy飞行策略分别对侦察者位置、跟踪者位置更新,提升了算法的全局探索和局部优化问题;然后,与其他改进算法对不同形态的基准函数进行优化求解,验证了ESSA算法在搜索能力、收敛性、稳定性等方面均具有明显的优势;最后,在相同的实验条件和参数下,将ESSA算法与其他算法进行对比实验。结果显示ESSA算法优化模型后的预测精度优于其他算法,能很快跳出局部最优,具有很好的稳定性。
采用ESSA算法对LSTM模型进行参数优化,明显缩减了参数设置的盲目性和时间成本,提高了优化效率和预测准确率。这种方法结合了现代优化算法与深度学习模型的优势,实现了对LSTM模型参数的寻优。在实验过程中,该算法在全局探索和局部优化方面均表现出色,达到了实验目的。不过,该算法在应用于特定类型数据方面的实验时存在一定的缺陷,后期需进一步考虑算法的通用性,以确保其在不同数据分布和特性下的鲁棒性;此外,还发现在获取实验结果过程中需要较长时间取得最优解,后续将持续优化算法,以提高计算效率,减少时间成本。
-
表 1 马氏珠母贝肉酶解产物及其超滤组分中蛋白质和小分子肽含量
Table 1 Content of protein and small molecule peptide in EP and its ultrafiltration components
$\overline {\mathit{\boldsymbol{X}}}{\bf \pm {{SD}}} $ 组别
Group蛋白质质量分数
Protein mass fraction/(g·kg−1)小分子肽质量分数
Small molecular peptide mass fraction/(g·kg−1)小分子肽占比
Proportion of small molecular peptides/%马氏珠母贝肉酶解产物组 EP group 616.65±50.15 451.14±15.39 73.38±4.11 EP-Ⅰ组 EP-Ⅰ group 575.67±36.22 91.91±17.89** 16.08±3.75** EP-Ⅱ组 EP-Ⅱ group 554.42±39.59 155.33±26.00** 28.13±5.09** EP-Ⅲ组 EP-Ⅲ group 490.29±24.35* 431.74±16.74 88.22±5.90** 注:与马氏珠母贝肉酶解产物组相比,*. P<0.05,**. P<0.01;表2同此 Note: Compared with the enzymatic hydrolysate of P. martensii; *. P < 0.05; **. P < 0.01. The same case in Table 2. 表 2 马氏珠母贝肉酶解产物及其超滤组分的氨基酸组成
Table 2 Amino acid composition of EP and its ultrafiltration components
g·kg−1 氨基酸种类
Amino acid species马氏珠母贝肉酶解产物
EP groupEP-Ⅰ组
EP-Ⅰ groupEP-Ⅱ组
EP-Ⅱ groupEP-Ⅲ组
EP-Ⅲ group天冬氨酸 Asp 49.10±0.35 31.80±0.41** 35.10±0.61** 46.10±0.61** 苏氨酸 Thr* 22.30±0.65 13.50±0.74** 21.10±0.45 19.10±0.66** 丝氨酸 Ser 18.10±0.28 11.20±0.31** 11.80±0.74** 10.40±0.54** 谷氨酸 Glu 55.20±0.71 58.90±0.11** 58.70±0.40** 72.20±0.34** 脯氨酸 Pro# 20.10±0.55 17.30±0.58 ** 17.90±0.27** 14.10±0.17** 甘氨酸 Gly# 35.30±0.82 29.00±0.34** 32.10±0.45** 36.00±0.45 丙氨酸 Ala# 30.00±0.41 26.20±0.62** 31.00±0.69 31.40±0.92 胱氨酸 Cys 0.70±0.04 0.50±0.07 0.60±0.16 1.50±0.10** 缬氨酸 Val*# 24.60±0.92 21.20±0.69* 25.20±0.24 24.70±0.75 蛋氨酸 Met* 11.60±0.54 12.80±0.20** 10.10±0.44** 10.00±0.45** 异亮氨酸 Ile*# 20.60±0.61 21.80±0.65 21.80±0.25 21.50±0.83 亮氨酸 Leu*# 31.10±0.44 28.50±0.74** 33.40±0.54** 36.50±0.38** 酪氨酸 Tyr 5.10±0.17 6.90±0.13** 9.10±0.18** 7.80±0.21** 苯丙氨酸 Phe*# 17.00±0.96 14.60±0.51** 18.90±0.17** 19.70±0.66** 赖氨酸 Lys* 37.40±0.83 23.20±0.45** 26.00±0.98** 32.60±0.88** 组氨酸 His 6.80±0.59 8.30±0.27* 8.90±0.25* 7.70±0.20 精氨酸 Arg 28.00±0.59 13.60±0.54** 20.10±0.88** 15.20±0.34** 氨基酸总和 Total amino acid, TAA 413.00±4.53 339.30±3.51 381.80±3.82 406.50±1.80 必需氨基酸 Essential amino acid, EAA 164.60±2.86 135.60±1.17 156.50±1.68 164.10±0.78 疏水性氨基酸 Hydrophobic amino acid, HAA 178.70±1.51 158.60±1.73 180.30±1.23 183.90±1.74 注:*. 必需氨基酸;#. 疏水性氨基酸 Note: *. Essential amino acid; #. Hydrophobic amino acid 表 3 马氏珠母贝肉酶解产物超滤组分的体外抗氧化能力
Table 3 Antioxidant capacity of EP ultrafiltration components in vitro
$\overline {\mathit{\boldsymbol{X}}}{\bf \pm {{SD}}} $ 组别
Group2,2-联氮-二(3-乙基-苯并噻唑-6-磺酸)
二铵盐自由基清除率
ABTS free radical scavenging activity/%2,2-联苯基-1-苦基肼基自由基清除率
DPPH free radical scavenging activity/%还原力A700 nm
Reducing power维生素C组 VC group 89.60±0.23 89.39±0.82 1.878±0.083 EP-Ⅰ组 EP-Ⅰ group 86.68±0.44** 84.54±1.47* 1.568±0.060** EP-Ⅱ组 EP-Ⅱ group 87.46±0.18** 85.13±2.57* 1.811±0.050 EP-Ⅲ组 EP-Ⅲ group 88.14±0.25** 86.34±1.21 1.811±0.046 注:与维生素C组相比,*. P < 0.05;**. P < 0.01 Note: Compared with VC group, *. P < 0.05; **. P < 0.01 表 4 马氏珠母贝肉酶解产物超滤组分对酒精性肝损伤小鼠体质量及其肝脏指数的影响
Table 4 Effects of EP ultrafiltration components on body mass and liver index of ALD mice
$\overline {\mathit{\boldsymbol{X}}}{\bf \pm {{SD}}} $ 组别
Group初次灌胃体质量
Body mass of primary gastric perfusion/g末次灌胃体质量
Body mass of last gastric perfusion/g肝脏指数
Liver index/(mg·g−1)空白对照组 Blank control group 23.63±1.05 25.31±0.71 3.55±0.21## 酒精模型组 Alcohol model group 24.37±0.83 26.35±0.87 4.40±0.26** 阳性对照组 Positive control group 23.35±1.40 25.40±1.29 3.97±0.37*# EP-Ⅰ组 EP-Ⅰ group 24.75±0.97 26.23±0.80 4.20±0.38** EP-Ⅱ组 EP-Ⅱ group 24.31±0.73 26.47±1.44 3.90±0.32## EP-Ⅲ组 EP-Ⅲ group 23.75±0.76 26.57±1.59 3.78±0.15## 注:与空白对照组相比,*. P<0.05,**. P<0.01;与酒精模型组相比,#. P<0.05,##. P<0.01;图2、图3、表5同此 Note: Compared with blank control group, *. P<0.05, **. P< 0.01; compared with alcohol model group, #. P < 0.05, ##. P<0.01. The same case in Figure 2–3 and Table 5. 表 5 马氏珠母贝肉酶解产物超滤组分对酒精性肝损伤小鼠肝脏中超氧化物歧化酶、谷胱甘肽、丙二醛和甘油三酯的影响
Table 5 Effects of EP ultrafiltration components on SOD, GSH, MDA and TG in liver of ALD mice
$\overline {\mathit{\boldsymbol{X}}}{\bf \pm {{SD}}} $ 组别
Group超氧化物歧化酶活力
SOD activity/(U·mg−1)谷胱甘肽
GSH/(nmol·mL−1)丙二醛
MDA/(nmol·mg−1)甘油三酯
TG/(nmol·mL−1)空白对照组 Blank control group 15.88±0.45## 2.64±0.27## 1.10±0.28## 0.02±0.01## 酒精模型组 Alcohol model group 14.75±0.40** 1.58±0.13** 1.90±0.60** 0.07±0.02** 阳性对照组 Positive control group 15.24±0.38 1.92±0.28** 1.19±0.24## 0.02±0.02## EP-Ⅰ组 EP-Ⅰ group 15.13±0.28** 2.07±0.29**## 1.44±0.28## 0.05±0.03** EP-Ⅱ组 EP-Ⅱ group 15.20±0.61** 2.14±0.38**## 1.22±0.35## 0.03±0.03## EP-Ⅲ组 EP-Ⅲ group 15.58±0.61## 2.25±0.45## 1.04±0.17## 0.03±0.01## -
[1] 章超桦, 吴红棉. 马氏珠母贝肉的营养成分及其游离氨基酸组成[J]. 水产学报, 2000, 24(2): 180-184. [2] 于志鹏, 武思佳, 赵文竹, 等. 海洋贝类蛋白源生物活性肽及肽组学的研究进展[J]. 食品工业科技, 2015, 36(22): 384-388. [3] 吴燕燕, 尚军, 李来好, 等. 合浦珠母贝肉短肽的分离及其抗氧化活性研究[J]. 食品工业科技, 2013, 33(7): 123-126. [4] 吴燕燕, 宫晓静, 李来好, 等. 风味蛋白酶水解合浦珠母贝肉制备抗菌肽人工神经网络法优化工艺[J]. 食品科学, 2011, 32(20): 63-68. [5] 章超桦, 刘亚, 杨萍, 等. 马氏珠母贝肉酶解蛋白抗疲劳功能的初步研究[J]. 中国海洋药物, 2006(4): 46-47. doi: 10.3969/j.issn.1002-3461.2006.04.011 [6] 曹文红, 吴红棉, 章超桦, 等. 马氏珠母贝肉酶解产物ACE抑制活性的研究[J]. 食品与发酵工业, 2008, 34(8): 60-64. [7] 韩丽娜, 秦小明, 林华娟, 等. 马氏珠母贝肉的醒酒作用机理初探[J]. 食品科技, 2010, 35(10): 180-183. [8] WU S, YANG C, XU N, et al. The protective effects of helix b surface peptide on experimental acute liver injury induced by carbon tetrachloride[J]. Digest Dis Sci, 2017, 62(6): 1537-1549. doi: 10.1007/s10620-017-4553-7
[9] MING C W, PEI L Z, CHANG X J, et al. Preliminary characterization, antioxidant activity in vitro and hepatoprotective effect on acute alcohol-induced liver injury in mice of polysaccharides from the peduncles of Hovenia dulcis[J]. Food Chem Toxicol, 2012, 50(9): 2964-2970. doi: 10.1016/j.fct.2012.06.034
[10] YAN S L, YANG H T, LEE H L, et al. Protective effects of maslinic acid against alcohol-induced acute liver injury in mice[J]. Food Chem Toxicol, 2014, 74(6): 1537-1549.
[11] BRANDON W E, SCHRUM L W, SCHMIDT C M, et al. Rodent models of alcoholic liver disease: of mice and men[J]. Alcohol, 2012, 46(8): 715-725. doi: 10.1016/j.alcohol.2012.08.004
[12] CRABB D W, GALLI A, FISCHER M, et al. Molecular mechanisms of alcoholic fatty liver: role of peroxisome proliferator-activated receptor alpha[J]. Alcohol, 2004, 34(1): 35-38. doi: 10.1016/j.alcohol.2004.07.005
[13] HETAL A K, SAMIR P. Management of alcoholic hepatitis: current concepts[J]. World J Hepatol, 2012, 4(12): 335-341. doi: 10.4254/wjh.v4.i12.335
[14] 胡滨, 李康林, 吴桥, 等. 猪血蛋白酶解物对小鼠急性酒精性肝损伤的保护作用[J]. 食品科学, 2018, 39(11): 185-190. doi: 10.7506/spkx1002-6630-201811029 [15] 刘鹏. 低分子量促乙醇代谢玉米肽制备与蛋白成分对活性贡献关系研究[D]. 广州: 华南理工大学, 2015: 17. [16] CAI X X, YAN A N, FU N Y, et al. In vitro antioxidant activities of enzymatic hydrolysate from Schizochytrium sp. and its hepatoprotective effects on acute alcohol-induced liver injury in vivo[J]. Mar Drugs, 2017, 15(4): 115-124. doi: 10.3390/md15040115
[17] UDENIGWE C C, ALUKO R E. Chemometric analysis of the amino acid requirements of antioxidant food protein hydrolysates[J]. Int J Mol Sci, 2011, 12(12): 3148-3161.
[18] 郭辉, 何慧, 韩樱, 等. 玉米肽对小鼠酒后肝脏乙醇脱氢酶活力的影响及醒酒机理[J]. 食品科学, 2011, 32(11): 265-269. [19] SARMADI B H, ISMAIL A. Antioxidative peptides from food proteins: a review[J]. Peptides, 2010, 31(10): 1949-1956.
[20] GIMÉNEZ B, ALEMÁN A, MONTERO P, et al. Antioxidant and functional properties of gelatin hydrolysates obtained from skin of sole and squid[J]. Food Chem, 2009, 114(3): 976-983. doi: 10.1016/j.foodchem.2008.10.050
[21] GUO C Y, JIANG T L, HUI H E, et al. Ultrafiltration preparation of potent bioactive corn peptide as alcohol metabolism stimulator in vivo and study on its mechanism of action[J]. J Food Biochem, 2013, 37(2): 161-167. doi: 10.1111/j.1745-4514.2011.00613.x
[22] 翟硕, 张海悦, 田田. 黑豆多肽的制备及其对乙醇脱氢酶活性的影响研究[J]. 食品安全质量检测学报, 2016, 7(12): 4864-4869. [23] 蒲月华, 邓旗, 杨萍, 等. 珍珠贝多肽体外抗氧化活性的研究[J]. 食品科技, 2016(11): 124-128. [24] 吴静, 胡晓, 杨贤庆, 等. 鸢乌贼酶解产物的抗氧化稳定性与功能特性[J]. 南方水产科学, 2016, 12(5): 105-111. doi: 10.3969/j.issn.2095-0780.2016.05.013 [25] 彭文锋, 钟政永. ADA与ALT、AST、GGT联合检测在肝脏疾病诊断中的意义[J]. 当代医学, 2011, 17(9): 4-6. doi: 10.3969/j.issn.1009-4393.2011.9.003 [26] 王佳佳, 赵莎莎, 杨最素, 等. 文蛤寡肽对小鼠急性肝损伤的保护作用[J]. 食品科学, 2017, 38(13): 190-195. doi: 10.7506/spkx1002-6630-201713031 [27] JE J Y, CHA J Y, CHO Y S, et al. Hepatoprotective effect of peptic hydrolysate from salmon pectoral fin protein byproducts on ethanol-induced oxidative stress in Sprague-Dawley rats[J]. Food Res Int, 2013, 51(2): 648-653. doi: 10.1016/j.foodres.2013.01.045
[28] 刘奇奇, 温久福, 区又君, 等. 急性操作胁迫对四指马鲅幼鱼肝脏组织结构和氧化应激的影响[J]. 南方水产科学, 2017, 13(5): 103-109. doi: 10.3969/j.issn.2095-0780.2017.05.014 [29] SCHEMITT E G, HARTMANN R M, COLARES J R, et al. Protective action of glutamine in rats with severe acute liver failure[J]. World J Hepatol, 2019, 11(3): 27-40.
[30] KOCH O R, PANI G, BORRELLO S, et al. Oxidative stress and antioxidant defenses in ethanol-induced cell injury[J]. Mol Aspects Med, 2004, 25(1/2): 191-198.
[31] LU J, LYU Y B, LI M T, et al. Alleviating acute alcoholic liver injury in mice with Bacillus subtilis co-expressing alcohol dehydrogenase and acetaldehyde dehydrogenase[J]. J Funct Foods, 2018, 49: 342-350. doi: 10.1016/j.jff.2018.09.006
[32] XIAO C Q, ZHOU F B, ZHAO M M, et al. Chicken breast muscle hydrolysates ameliorate acute alcohol-induced liver injury in mice through alcohol dehydrogenase (ADH) activation and oxidative stress reduction[J]. Food Funct, 2018, 9(2): 1-42.
-
期刊类型引用(4)
1. 路俊怡,姜立妍,龙凯,王韬,吴正理,李艳红. HcTLR1通过MyD88-NF-κB信号通路参与三角帆蚌抗菌免疫应答. 南方水产科学. 2024(06): 19-30 . 本站查看
2. 张林宝,孙伟,宁璇璇,蔡文贵,张喆,陈海刚,马胜伟,贾晓平. 菲律宾蛤仔脂多糖诱导的肿瘤坏死因子α VpLITAF基因克隆及其对微生物侵染的应答研究. 海洋与湖沼. 2014(04): 895-901 . 百度学术
3. 王志新,梁海鹰,杜晓东,朱家萍. 蛋白质组学在贝类研究中的应用. 生命科学研究. 2014(02): 184-188 . 百度学术
4. 张晗,姜敬哲,王江勇. 蛋白质组学研究的发展及其在经济贝类研究中的应用. 广东农业科学. 2013(22): 170-173+182 . 百度学术
其他类型引用(4)