贺 越 赵圣国 张晓音 郑 楠 王加启*
(1.中国农业科学院北京畜牧兽医研究所,动物营养学国家重点实验室,北京 100193;2.中国农业科学院北京畜牧兽医研究所,农业农村部奶产品质量安全风险评估实验室(北京),北京 100193)
反刍动物瘤胃微生物能利用尿素分解产生的氨合成机体可利用的微生物蛋白,所以尿素可作为一种经济的蛋白质饲料替代品[1-2]。瘤胃尿素的分解依赖于瘤胃细菌脲酶,脲酶催化尿素水解的速度是尿素自然分解速度的1014倍[3],尿素饲喂不当容易造成动物氨中毒[4],这严重制约了尿素在动物生产上的广泛应用。脲酶抑制剂是脲酶活性调控研究的前沿。
细菌脲酶是一种镍金属酶[5],通常由结构蛋白(UreA、UreB、UreC)和辅助蛋白(UreD/H、UreF、UreG、UreE)组成[6]。不同微生物来源的脲酶的结构蛋白具有差异,如产气克雷伯氏菌脲酶由UreA、UreB、UreC 3个亚基聚合形成(UreABC)3[7],耶尔森氏菌脲酶为[(UreABC)3]4[8],幽门螺杆菌脲酶为[(UreAB)3]4[9]。细菌脲酶的活性中心位于序列高度保守的UreC亚基[10],每1个活性中心包含2个镍离子,2个镍离子由1个羧化的赖氨酸连接[11],活性中心底物的进入和产物的排出由1个柔性的flap区域调控[12]。脲酶活性中心及flap区域是活性研究的重要靶点。
目前已经报道了多种脲酶抑制剂,其中多为化学合成的化合物。乙酰氧肟酸是目前唯一被批准商用的脲酶抑制剂[13],但是在瘤胃中稳定性差、对动物的健康有影响[14],所以开发安全、高效的脲酶抑制剂是目前研究的重点。自然界蕴含丰富的天然产物资源,其中包括大量具有生物活性、化学结构多样的植物来源的天然化合物,这些化合物在疾病治疗、粮食资源等方面应用广泛[15]。在这种背景下,研究植物来源的天然产物作为脲酶抑制剂调控脲酶活性在动物生产中具有重要价值。本研究拟运用分子对接的方法对植物源化合物库进行筛选,并进一步结合瘤胃细菌脲酶活性试验进行评价,以期获得对瘤胃脲酶具有抑制效果的天然化合物,为天然化合物抑制瘤胃细菌脲酶的研究提供科学依据。
试验用瘤胃液来自中国农业科学院昌平基地荷斯坦奶牛[泌乳中期,体重(550±23) kg];试验用植物源化合物库信息及植物源天然化合物来源于上海陶素生化科技有限公司,纯度>98%;刀豆脲酶、尿素、4-羟乙基哌嗪乙磺酸(HEPES)粉末、四水合酒石酸钾钠、还原性谷胱甘肽(GSH)、L-半胱氨酸(L-Cys)、硼酸(BA)购于西格玛奥德里奇(上海)贸易有限公司,分析纯;偏磷酸购于国药集团化学试剂有限公司,分析纯;纳氏试剂购于上海阿拉丁生化科技股份有限公司,分析纯;二甲基亚砜(DMSO)购于北京索莱宝科技有限公司;尿素氮含量测试试剂盒购于南京建成生物工程研究所;试验用气体购于北京市北温气体制造厂;试验均使用超纯水。
试验用到的仪器和软件包括Auto Dock vina软件、Discovery Studio 4.5 Client软件、超纯水仪(F6BA99920)、pH计(SG2)、台式高速冷冻离心机(sigma 3-18k)、超声波细胞粉碎机(UP-250)、生化振荡培养箱(HZQ-F100)、磁力搅拌器(IKA big squid)、酶标仪(Thermo Electron Varioskan Flash)、厌氧培养箱(PLAS-LAB 855-AC & 855-ACB)。
50 mmol/L HEPES缓冲液:称取11.9 g HEPES粉末,适量去离子水溶解后,调节pH至7.5,定容至1 L。
50 mmol/L尿素缓冲液:称取0.3 g尿素,用HEPES缓冲液溶解,定容至100 mL。
瑞香素溶液:10 mg瑞香素中加入561.4 μL DMSO,即配制成浓度为100 mmol/L的瑞香素母液,后续使用根据所需浓度进行稀释。
DMSO溶液:DMSO溶液中DMSO的含量与对应试验组化合物中DMSO含量一致。
硫醇试剂溶液:称取L-Cys 1.21 g,HEPES缓冲液定容至100 mL,称取GSH 3.07 g,HEPES缓冲液定容至100 mL,2种硫醇试剂母液浓度均为100 mmol/L,后续使用根据所需进行稀释。
BA溶液:称取BA 0.62 g,HEPES缓冲液定容至100 mL,BA溶液母液浓度为100 mmol/L,后续使用根据所需进行稀释。
本试验使用的瘤胃细菌脲酶结构参考实验室之前的同源建模[16]。利用Auto Dock vina软件对瘤胃细菌脲酶进行去除水分子、加氢及计算电荷等处理,对植物源化合物库中的化合物进行同样处理,而后选择细菌脲酶UreC亚基的第308位氨基酸残基至第336位氨基酸残基的区域和活性中心区域进行对接,选择最佳结合模型。进一步测定对接结果较好的4种天然化合物对瘤胃细菌脲酶的抑制效果。
2.2.1 瘤胃细菌脲酶的提取
挑选3头荷斯坦瘘管奶牛,采集瘤胃内容物100 mL,通过4层无菌纱布过滤,将滤液(瘤胃液)分装到50 mL无菌离心管中,并置于液氮中保存。
瘤胃液用流水解冻后,300×g离心5 min,收集上清液到新的离心管中,12 000×g离心10 min,弃上清,用25 mL HEPES缓冲液重悬沉淀,12 000×g离心10 min,弃上清,用25 mL HEPES缓冲液重悬沉淀,将其置于冰上,用超声波细胞粉碎机进行破碎(150 W,超声5 s,停5 s,超声10 min),12 000×g离心10 min,收集上清,即瘤胃细菌脲酶。离心过程均为4 ℃。
2.2.2 4种天然化合物对瘤胃细菌脲酶的抑制效果的测定
本试验测定脲酶活性采用纳氏试剂法。试验组中,取40 μL瘤胃细菌脲酶于96孔板中,加入40 μL化合物溶液(终浓度10 μmol/L)和100 μL尿素缓冲液(终浓度28 mmol/L);试验阴性对照组中,取40 μL瘤胃细菌脲酶于96孔板中,加入40 μL化合物溶液(终浓度10 μmol/L)和100 μL HEPES缓冲液;脲酶组中,取40 μL瘤胃细菌脲酶于96孔板中,加入40 μL DMSO溶液和100 μL尿素缓冲液;脲酶阴性对照组中,取40 μL瘤胃细菌脲酶于96孔板中,加入40 μL DMSO溶液和100 μL HEPES缓冲液。以上每组设置3个重复,样品加好后,置于生化振荡培养箱中,37 ℃孵育30 min,加入10 μL酒石酸钾钠溶液,混合均匀,加入10 μL纳氏试剂溶液,混合均匀后,室温放置10 min。而后置于酶标仪中,设置420 nm波长,测定吸光度(OD)值。剩余活性计算公式如下:
剩余活性(%)=[(试验组OD值-试验阴性
对照组OD值)/(脲酶组OD值-脲酶
阴性对照组OD值)]×100。
根据4种天然化合物对瘤胃细菌脲酶的抑制结果,对其中抑制效果最好的化合物瑞香素作进一步的评价。
测定瑞香素对瘤胃细菌脲酶的IC50值采用纳氏试剂法,并且以刀豆脲酶作为对照,瑞香素对瘤胃细菌脲酶IC50值测定浓度设置为0、10、20、80、100、160、200 μmol/L,瑞香素对刀豆脲酶IC50值测定浓度设置为0、2、10、20、30、40、45 μmol/L,每个浓度3个重复,测定体系及相应对照设置、剩余活性计算同2.2.2。通过GraphPad Prism 6 进行曲线拟合,得到IC50值。
瑞香素对瘤胃脲酶动力学测定采用纳氏试剂法,瑞香素设0、80、160、200 μmol/L 4个终浓度,尿素设0、0.875、1.750、3.500、7.000、14.000 mmol/L 6个终浓度,每个处理3个重复,其中尿素0 mmol/L为试验阴性对照组,测定体系同2.2.2。根据测得的OD值计算出相应处理下的反应速率(V),在GraphPad Prism 6中进行反应速率的倒数关于尿素浓度(S)的倒数的线性拟合,进一步根据米氏方程计算得到酶动力学参数——米氏常数(Km)和最大反应速率(Vmax),方程如下:
本试验选择L-Cys、GSH 2种硫醇试剂和BA这种脲酶Ni2+竞争性抑制剂。脲酶活性测定选择纳氏试剂法,每个处理中添加40 μL瘤胃细菌脲酶和100 μL尿素缓冲液,其中对照组中添加40 μL DMSO溶液,瑞香素组中添加40 μL瑞香素溶液,L-Cys、GSH、BA组中添加相应的溶液40 μL,L-Cys/GSH/BA+瑞香素组中各添加20 μL,瑞香素终浓度为200 μmol/L,L-Cys、GSH、BA终浓度都为0.25 mmol/L,用HEPES缓冲液代替尿素缓冲液作为相对应的阴性对照组,计算方法同2.2.2。利用GraphPad Prism 6进行作图分析。
利用Discovery Studio 4.5 Client软件对瑞香素与瘤胃脲酶活性中心对接的结果进行作图,分析瑞香素与瘤胃脲酶相互结合的氨基酸位点及作用力。
瘤胃液的采集同2.2.1,经过纱布过滤后,用20 mL无菌注射器注射入进过厌氧及灭菌后的血清瓶中,冰上保存。
厌氧稀释液及厌氧培养基配制参考实验室之前的配方[17]。瑞香素溶液和尿素溶液用厌氧稀释液配制,浓度分别为3.8 mmol/L和0.4 mol/L。在厌氧培养箱中将20 mL新鲜瘤胃液与40 mL厌氧培养基混合均匀后分装到螺口培养管中,每管9 mL,试验组中加入500 μL瑞香素溶液,对照组中则加入500 μL DMSO溶液(DMSO含量与瑞香素溶液一致,用厌氧稀释液配制),每组设置3个重复,室温下反应30 min后,加入500 μL尿素溶液,置于厌氧培养箱中39 ℃培养,分别于0、2、4、6、12 h采集200 μL培养液,采集的样品12 000×g,4 ℃离心5 min,收集上清,加入400 μL 25%偏磷酸,于-20 ℃保存。尿素氮含量利用尿素氮测试盒(二乙酰肟比色法)进行测定。
试验数据利用SPSS 22.0统计软件进行统计学分析,结果以平均值±标准差表示。在统计分析前,先利用软件中分析-探索-Shapiro-Wilk检验对数据进行正态性检验,服从正态分布的2组数据间的比较采用分析-比较均值-独立样本t检验进行差异性比较;服从正态分布的多组数据间的比较采用分析-比较均值-单因素ANOVA,方差齐时用LSD检验进行多重比较。P<0.05表示结果差异显著。
通过分子对接对植物来源化合物库进行多轮筛选,最后得到4种天然化合物作为瘤胃脲酶候选抑制剂,表1列出了这4种化合物分子对接得分以及化学结构式。
表1 4种天然化合物分子对接得分及化学结构式信息Table 1 Molecular docking score and chemical formula information of four natural compounds
测定分子对接筛选出的4种天然化合物对瘤胃脲酶的抑制效果,结果如图1所示,在浓度一致的情况下,瑞香素对瘤胃脲酶的抑制效果最好。
图1 4种天然化合物作用后瘤胃细菌脲酶的剩余活性Fig.1 Residual activity of ruminal microbial urease interacted by four natural compounds
为了表征瑞香素对脲酶的抑制效果,通过测定不同浓度的瑞香素对瘤胃脲酶和刀豆脲酶的抑制率,并计算瑞香素对2种脲酶的IC50值。结果如图2所示,瑞香素对瘤胃脲酶的IC50值为(93.43±3.87) μmol/L(图2-A),瑞香素对刀豆脲酶的IC50值为(16.47±0.79) μmol/L(图2-B)。
图2 不同瑞香素浓度作用后瘤胃细菌脲酶(A)、刀豆脲酶(B)剩余活性Fig.2 Residual activity of ruminal microbial urease (A) and jack been urease (B) with different concentrations of daphnetin
表2 不同瑞香素浓度对瘤胃细菌脲酶动力学参数的影响Table 2 Effects of kinetic parameters of rumen bacterial urease by different concentrations of daphnetin
图3 Lineweaver-Burk双倒数作图法测定瑞香素对瘤胃细菌脲酶的酶动力学参数Fig.3 Kinetic parameters of rumen bacterial urease tested through Lineweaver-Burk plots of different concentrations of daphnetin
2种硫醇试剂L-Cys、GSH以及Ni2+竞争性抑制剂BA对瘤胃脲酶均具有一定的抑制效果,但是抑制效果均显著低于瑞香素(图4)。当瑞香素与L-Cys或GSH联合作用时,瑞香素对瘤胃脲酶的抑制效果显著降低(P<0.05)(图4-A和图4-B),当瑞香素与BA联合作用时,瑞香素对瘤胃脲酶的抑制效果变化不显著(P>0.05)(图4-C)。这一结果说明,瑞香素对瘤胃脲酶的抑制可能是通过与活性中心周围的巯基相互作用。
通过Discovery Studio 4.5 Client软件分析得到瑞香素与瘤胃脲酶的最佳结合模式(图5)。结果显示,瑞香素上的2个羟基分别与氨基酸残基Ala-166和His-245形成氢键,其中1个羟基还与氨基酸残基Lys-216产生离子相互作用(盐桥),瑞香素的苯环结构与氨基酸残基Ala-362形成疏水结合Pi-Alkyl作用,瑞香素的杂环结构能与氨基酸残基Ala-362和Cys-318形成疏水结合Pi-Alkyl作用,并且还能与Ala-362形成碳氢键。
由图6结果可知,在瘤胃微生物体外厌氧培养试验中,添加瑞香素与只添加等量DMSO的对照组相比,0 h时,2组培养体系中尿素含量差异不显著(P>0.05),而在2、4、6 h时,瑞香素体系中尿素含量显著高于对照组(P<0.05)。这一结果说明,添加瑞香素在一定时间内能有效地降低瘤胃微生物对尿素的分解。
标注相同小写字母表示差异不显著(P>0.05),标注不同小写字母表示差异显著(P<0.05)
A图中深蓝色球体为瘤胃细菌优势脲酶活性中心镍离子,玫红色结构为瑞香素。B图中圆形球体表示脲酶氨基酸残基;实线表示瑞香素;虚线表示瑞香素与氨基酸之间的相互作用,其中深绿色表示氢键,浅绿色表示碳氢键,橙色表示离子相互作用,玫红色表示Pi-Alkyl作用。
标注*表示差异显著(P<0.05),未标注*表示差异不显著(P>0.05)。
分子对接是目前药物筛选常用的方法[18],这一方法在对脲酶有抑制效果的天然化合物研究中也有诸多应用。如Kataria等[19]利用刀豆脲酶作为目标蛋白,通过分子对接技术从多种天然酚类化合物中筛选出了桑色素、白藜芦醇等化合物。有研究发现,桑色素及其类似物对幽门螺杆菌脲酶具有较好的抑制效果[20],白藜芦醇能通过抑制幽门螺杆菌脲酶而影响幽门螺杆菌在胃肠环境的生存和定植[21],也有研究表明,白藜芦醇的衍生物对刀豆脲酶具有抑制作用[22]。这些研究提示在众多天然化合物种中筛选有效的脲酶抑制剂,分子对接技术是一种有效的方法。虽然脲酶的活性中心具有高度保守的特点,但是同一种化合物对不同微生物脲酶的抑制效果具有一定的差异。Li等[23]研究发现,小檗碱对幽门螺杆菌脲酶的IC50值为2.40 mg/mL,而对刀豆脲酶的IC50值大于3.06 mg/mL,这提示我们靶向性的筛选脲酶抑制剂很有必要。本试验利用瘤胃优势脲酶基因簇的结构基因[24]进行同源建模,将得到的蛋白结构通过分子对接对植物来源化合物库进行筛选,并测定了天然化合物对瘤胃细菌脲酶的抑制效果,最终确定瑞香素对瘤胃细菌脲酶抑制效果最好。
香豆素及其衍生物是分布广泛的天然化合物,具有多种生物活性,一直以来备受化学家和药物化学家的青睐。研究发现,这些天然化合物中有多种具有脲酶抑制效果[25]。瑞香素是一种天然香豆素衍生物,具有抗菌、抗炎、抗肿瘤等特性[26-27]。幽门螺杆菌产生的脲酶有利于幽门螺杆菌在人体胃肠环境中的生存,从而引发多种疾病,瑞香素能有效降低幽门螺杆菌活性[28],并且也有研究发现瑞香素能抑制幽门螺杆菌脲酶活性[29]。本试验研究发现,瑞香素对瘤胃细菌脲酶以及刀豆脲酶均具有抑制效果,并且对刀豆脲酶的抑制效果更好,猜测可能和本试验所用瘤胃细菌脲酶不纯有关。这一现象也出现在天然化合物黄连碱的研究中,Li等[30]研究显示,黄连碱对刀豆脲酶的IC50值低于对幽门螺杆菌脲酶。在后续天然化合物研究中可以考虑对瘤胃细菌脲酶进行提纯后分析。
随着脲酶结构以及抑制剂-脲酶结构的解析,脲酶抑制剂对脲酶产生抑制作用的结构机制有:一方面是脲酶抑制剂进入脲酶活性口袋,与活性中心的镍离子发生相互作用,阻碍底物与镍离子的结合[31];另一方面则是脲酶抑制剂与flap区域发生相互作用,限制这个皮瓣结构的运动而产生抑制作用[32]。有研究表明,含硫醇化合物是一种脲酶活化剂,它能与脲酶氨基酸中的巯基发生相互作用[33],脲酶flap区域的关键氨基酸是1个半胱氨酸,当硫醇化合物和脲酶抑制剂同时存在时,硫醇可通过与半胱氨酸的巯基结合而阻碍脲酶抑制剂与半胱氨酸的相互作用,从而降低脲酶抑制剂对脲酶的抑制效果[34]。本试验结果显示,存在硫醇试剂时,瑞香素的抑制效果降低,说明瑞香素可能是通过与瘤胃细菌脲酶flap区域的半胱氨酸发生相互作用而产生抑制效果,同时分子对接结果也显示瑞香素能与半胱氨酸形成比较强的相互作用。氨基酸残基Lys-216是连接活性中心2个镍离子的关键氨基酸,氨基酸残基His245与其中的1个镍离子发生配位[16]。本试验分子对接结果显示了瑞香素结构中的1个羟基能与这2个关键氨基酸残基发生相互作用,而BA作用下瑞香素对瘤胃脲酶的抑制效果没有明显降低,猜测可能受到瘤胃脲酶纯度较低的影响。很多属于天然化合物的脲酶抑制剂也能通过与脲酶flap区域的半胱氨酸以及与活性中心的关键氨基酸残基发生相互作用而对脲酶产生抑制效果,如槲皮素[35]、芦荟树脂和枸杞中的天然脲酶抑制剂成分[36]、两面针中对脲酶具有抑制效果的天然化合物[37]等。
有研究发现,瘤胃微生物产生的脲酶水解尿素是一个很快的过程,细菌脲酶能在30 min到2 h内快速水解来自饲粮或者血液循环而来的尿素[38],本研究对瘤胃微生物体外厌氧培养的结果显示,培养体系中尿素含量在2 h显著降低。一些脲酶抑制剂如乙酰氧肟酸、正丁基硫代磷酸三胺、铋盐类化合物[39]等能降低尿素的分解,但是这些化合物也存在诸多问题,如乙酰氧肟酸在瘤胃内易降解,瘤胃微生物对正丁基硫代磷酸三胺能产生适应性,还有些化合物对人和动物机体有损害。本试验得到的天然化合物瑞香素在降低尿素分解方面具有很大的潜力,但是瑞香素在瘤胃内的稳定性、微生物对瑞香素是否具有适应性等还需要进一步的研究。
① 植物源化合物库中对瘤胃细菌脲酶抑制效果最好的化合物是瑞香素,该化合物对瘤胃细菌脲酶的IC50值为93.43 μmol/L,是瘤胃细菌脲酶的混合型抑制剂。
② 瑞香素能通过与瘤胃细菌脲酶活性中心的多个氨基酸残基以及flap区域的半胱氨酸相互作用而阻碍脲酶对尿素的分解。
③ 瑞香素能降低瘤胃微生物体外厌氧培养过程中尿素的分解。