一种热稳定的脂肪酶及制备方法与应用
【专利摘要】本发明公开了一种热稳定的脂肪酶及制备方法与应用。该脂肪酶为氨基酸序列如SEQ ID NO.1所示的S2?210脂肪酶、氨基酸序列如SEQ ID NO.2所示的S8?214脂肪酶、氨基酸序列如SEQ ID NO.3所示的S14?216脂肪酶或氨基酸序列如SEQ ID NO.4所示的S191?241脂肪酶。本发明通过长时间的高温分子动力学模拟模拟耶氏解脂酵母脂肪酶2的解折叠过程，分析了蛋白解折叠的主要区域和湿熔球态形成的关键步骤，有效的筛选出动力学稳定性改造的关键位点。该热稳定的脂肪酶耐热，半衰期长，特别适合在工业上进行应用。
【专利说明】
一种热稳定的脂肪酶及制备方法与应用
技术领域
[0001] 本发明属于酶工程领域，特别涉及一种热稳定的脂肪酶及制备方法与应用。
【背景技术】
[0002] 脂肪酶(EC 3.1.1.3)全称三酰甘油水解酶，被广泛应用于饲料工业、食品加工、化 妆品、洗涤剂、生物医学和生物能源等方面。
[0003] 耶氏解脂假丝酵母(Yarrowia lipolytica)是一种非常规酵母，属于食品安全型 酵母，目前通过连续基因干扰分析，，已发现其能编码16种脂肪酶Lip2、Lip4、Lip5、Lip7~ 19。这些脂肪酶同源性大小不一，且性质差异较大。其中，耶氏解脂酵母脂肪酶2(Lip2)是主 要的胞外分泌脂肪酶，其对中长链脂肪酸甘油三酯(C12~C16)具有较高的催化活性，被广 泛应用于油脂水解、污水处理、食品加工、生物能源、化学合成以及胰腺缺乏症的治疗等多 个领域，具有广泛的应用前景，但当环境温度超过40°C时，该脂肪酶迅速失活，在生产加工、 储存、运输等过程中造成较大的损失，成为其在工业上应用的主要瓶颈。因此，有必要对耶 氏解脂酵母脂肪酶2改造，提高其热稳定性，以拓宽其应用范围。
[0004] 在天然蛋白结构中引入新的二硫键是改造蛋白热稳定性的一种重要方法，被广泛 应用于各类蛋白的改造，特别是在药物储存和工业应用改造领域。随着预测二硫键算法的 不断发展，其中应用较为广泛的是DbD和M0DIP算法，通过对C e和SY原子进行定位，快速搜索 并精确地估算X3二面角，筛选符合条件的二硫键。这两种算法都是基于静态晶体结构的计 算，对二硫键的预测具有高准确性，并且已经成功的应用于米黑根毛霉脂肪酶、南极假丝酵 母脂肪酶B、华根霉脂肪酶和铜绿假单胞菌脂肪酶等多种脂肪酶的热稳定性设计，但是对筛 选热稳定型二硫键的准确度仍然十分有限。
[0005] 目前对蛋白质热稳定性改造主要分为定向进化、半理性设计和理性设计，其中定 向进化和半理性设计需要大量的人力物力，而理性设计需要对蛋白质结构与功能有全面的 理解。目前理性设计主要根据蛋白质的晶体静态结构的计算，目前主要利用一下2类筛选方 法：
[0006] (1)基于蛋白质柔性区域的设计，如B-FITTER、FIRST、分子动力学模拟筛选波动较 大的区域。
[0007] (2)基于经验力场对蛋白热力学自由能进行计算，筛选蛋白质热力学稳定的突变 位点，如F0LDX力场等。
[0008] 以上均没有深刻理解蛋白质解折叠动力学过程，亦无涉及蛋白质动力学稳定性的 设计。
【发明内容】
[0009] 本发明的首要目的在于克服现有技术的缺点与不足，利用分子动力学模拟筛选出 脂肪酶解折叠和湿熔球态形成的关键区域，并分析湿熔球态的形成，有针对性地引入二硫 键突变，抑制蛋白湿恪球态和共价聚沉的形成，提尚热稳定型^?硫键筛选的精确性，并提供 一种热稳定的脂肪酶。
[0010] 本发明的另一目的在于提供所述热稳定的脂肪酶的制备方法。
[0011] 本发明的再一目的在于提供所述热稳定的脂肪酶的应用。
[0012] 本发明的目的通过下述技术方案实现:一种热稳定的脂肪酶，是在耶氏解脂假丝 酵母脂肪酶2的解折叠和湿熔球态形成的关键区域引入二硫键突变；
[0013] 所述的关键区域为耶氏解脂假丝酵母脂肪酶2的第180位~280位氨基酸序列。
[0014] 所述热稳定的脂肪酶优选为S2-210脂肪酶、S8-214脂肪酶、S14-216脂肪酶或 S191-241脂肪酶；
[0015] S2-210脂肪酶的氨基酸序列如下：
[0017] S8-214脂肪酶的氨基酸序列如下：
[0019] S14-216脂肪酶的氨基酸序列如下：
[0021] S191-241脂肪酶的氨基酸序列如下：
[0023]编码所述热稳定的脂肪酶的核苷酸序列，如下所示：
[0024] S2-210脂肪酶的核苷酸序列如下：
[0026] S8-214脂肪酶的核苷酸序列如下：
[0028] S14-216脂肪酶的核苷酸序列如下：
[0030] S191-241脂肪酶的核苷酸序列如下：
[0032]所述热稳定的脂肪酶的制备方法，包括如下步骤：
[0033] (1)通过Gromacs分子动力学软件模拟耶氏解脂酵母脂肪酶2的解折叠过程，通过 对所得轨迹的波动性分析与蛋白表面水分子的统计，分析出脂肪酶解折叠和湿熔球态形成 的关键区域，并进一步对β折叠尾巴区域分区，分析湿熔球态形成的步骤，最后使用 Disulfide byDesign软件算法筛选潜在二硫键突变位点；
[0034] (2)通过反向PCR突变，将所选出的氨基酸位点突变为半胱氨酸，并转入大肠杆菌 工程菌中，进行扩增培养与质粒提取、测序；
[0035] (3)将测序正确的突变质粒以Pmel限制性内切酶线性化处理，并电击转化入感受 态毕赤酵母X33中，进一步通过博来霉素抗性平板筛选和BMMY-罗丹明B产酶平板筛选，得到 对应的突变工程菌；
[0036] (4)将突变工程菌在YH)液体培养基中进行扩繁培养后转接至BMGY液体培养基进 行去抑制培养，最后接种BMMY液体培养基进行发酵，将菌液离心获取上清粗酶液
[0037] (5)使用超滤管将粗酶液进行超滤浓缩后，使用镍柱一步法纯化，分离出带组氨酸 标签的脂肪酶蛋白，并以还原性SDS-PAGE检测蛋白纯度，获取纯化后的脂肪酶；
[0038] (6)通过DTNB法测定蛋白质中的游离巯基浓度和总巯基浓度，计算蛋白分子内的 二硫键数目，检测引入的二硫键突变是否成键；
[0039] (7)通过p-NPP比色法和DSF荧光检测测定突变脂肪酶的热稳定性指标:保温15min 后残余50%活性的温度(T5Q)、5(TC下的半衰期（t1/2)、最适反应温度(T Qpt)以及蛋白熔解温 度(Tm);
[0040] (8)通过不同尿素诱导解折叠测定湿熔球态的形成所需要的浓度，确定突变脂肪 酶动力学稳定性大小；
[0041 ] (9)将突变脂肪酶在50°C下保温不同的时间，以非还原性SDS-PAGE检测突变脂肪 酶的抗聚沉能力大小，得到热稳定的脂肪酶。
[0042] 步骤（1)中所述的脂肪酶解折叠和湿熔球态形成的关键区域为耶氏解脂假丝酵母 脂肪酶2的第180位~280位氨基酸序列，即β折叠尾巴区域。
[0043] 步骤（1)中所述的湿熔球态形成的步骤包括初级步骤和关键步骤；耶氏解脂假丝 酵母脂肪酶2底部无规卷曲1〇〇ρ4(氨基酸207位至221位）与Ν端区域(氨基酸1位至13位）的 解离是湿熔球态形成的初级步骤，1〇〇ρ3结构(氨基酸228位至246位）的解离是湿熔球态形 成的关键步骤。
[0044] 步骤(1)中所述的解折叠的区域优选为β折叠尾巴区域(氨基酸180位至280位）。
[0045] 步骤(2)中所述的大肠杆菌工程菌优选为大肠杆菌Τ0Ρ10。
[0046] 步骤(4)中所述的培养的条件优选为摇瓶培养96小时。
[0047] 所述热稳定型脂肪酶耐热，半衰期长，特别适合在工业上进行应用。
[0048] 本发明相对于现有技术具有如下的优点及效果：
[0049] 本发明利用生物信息学对蛋白质改造属于理性设计的方法。该方法相对于传统的 非理性改造方法而言，能节约大量的人力物力，但需要对蛋白质分子进行全面的分析，方能 找筛选出有效的改造点。
[0050] 本发明通过长时间的高温分子动力学模拟筛选出耶氏解脂酵母脂肪酶2解折叠和 湿熔球态形成的关键区域，并分析湿熔球态的形成，关注蛋白质在动力学过程中的变化，有 针对性地引入二硫键突变，大大提高热稳定型二硫键筛选的准确性。
[0051] 1)本发明通过在脂肪酶底部无规卷曲1〇〇ρ4(氨基酸207位至221位)与Ν端区域(氨 基酸1位至13位)之间引入二硫键，加强了底部无规卷曲1 〇〇ρ4与Ν端的联系，稳定了 β折叠尾 部区域(氨基酸180位至280位)的稳定，相对于亲本脂肪酶，S2-210、S8-214、S14-216突变体 脂肪酶的T m值分别提升了 10.48 °C、7.50 °C和5.67 °C，50 °C半衰期分别达到198、43.87、 17.24 分钟，分别提升了 120、26.58、10.44 倍；
[0052] 2)本发明通过在脂肪酶β7、β8折叠之间的无规卷曲1〇〇ρ3(氨基酸228-246位)与β6 折叠和α5螺旋间的无规卷曲（氨基酸188位至氨基酸198位)之间引入的二硫键突变，相对于 亲本脂肪酶，S191-241突变体脂肪酶的1?值提升了5.63°C，50°C半衰期达到14.87分钟，提 升了 9倍。
[0053]以上所述的二硫键突变均有效地抑制了蛋白解折叠、减缓湿熔球态的形成，并减 少了脂肪酶之间的共价聚沉反应。
【附图说明】
[0054] 图1是不同温度下氨基酸的均方根涨幅曲线图。
[0055] 图2是熔球态形成区域分析图；其中：图a为盖子1区域、β折叠尾巴区域、财斤叠头部 区域4 Α范围内水分子数量随时间的变化曲线；图b为β折叠尾巴区域的4个亚区4 Α范围内 的水分子数量随时间变化曲线。
[0056]图3是β折置尾巴区域解折置波动曲线图；其中，图a为06、07、08、09折置的均方根 波动随模拟时间的变化曲线，图b为1〇〇ρ2、1〇〇ρ3、1〇〇ρ4的均方根波动随模拟时间的变化曲 线。
[0057] 图4是BMMY-罗丹明B产酶平板筛选图；其中，A-F分别为1^?2、38-214、32-210、314- 216、S191-241、S4-266在紫外光下脂肪酶的水解光圈。
[0058]图5是还原性SDS-PAGE检测蛋白纯度的结果图；其中：泳道1-7分别为1^?2、52_ 210、8-214、514-216、5191-241、5198-242、54-266突变脂肪酶的蛋白条带。
[0059] 图6是尿素诱导解折叠 bis-ANS荧光曲线图。
[0060] 图7是非还原性SDS-PAGE蛋白检测的结果图；其中：图&、13、(：、(1、6分别代表1^2、 52-210、58-214、514-216、5191-241在50°(：下热处理不同时间。
【具体实施方式】
[0061] 下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限 于此。
[0062] 材料与试剂:pPICZaA_Lip2毕赤酵母表达载体由上海金斯瑞生物公司全基因合成 与构建;质粒提取试剂盒购自Omega贸易有限公司，K0D-PLUS突变试剂盒购自东洋纺公司， ProteinThermal Shift筛选试剂盒购自Thermo公司；T0P10大肠杆菌感受态细胞购自天根 生物公司，突变引物由上海生工生物工程公司合成；Pmel限制性内切酶购自New England Biolabs公司；PCR产物纯化回收试剂盒购自大连宝生物公司；电转仪购自Bio-Rad公司； 0^、1^^+26〇(^11、￥?0、81?^、81?^培养基均按照1鮮1杜(^611毕赤酵母表达试剂盒操作手册 配制，镍柱纯化试剂盒购自上海生工生物工程公司，其余试剂均为国内外购买的分析纯级 别。
[0063] 实施例1耶氏解脂酵母脂肪酶2的解折叠模型建立与二硫键的筛选。
[0064] (1)从RCSB PDB晶体数据库中搜索下载Lip2的晶体结构(PDB ID:300D)，精度为 1.7 A，晶体为七聚体，每个聚体各有不同程度的原子缺失，因此选取骨架原子完整的A链作 为初始模型，使用Swiss-pdb Viewer软件对A链缺失的侧链氨基酸进行修复，并使用 Disulfide byDesign软件，将A链中Cysl20_Cysl23二硫键重构，在模拟中保留A链的所有结 晶水，除了三联体中心的HIS289使用HID质子化状态，其余组氨酸均设置为Gromacs默认下 pH=7.00的质子化状态，得到预处理模型；
[0065] (2)使用6如1^03 5.04软件对预处理模型进行三维建模，使用4!1^6"€9938-匪1?- ILDN力场的参数，将蛋白放置于十二面体的菱形盒子中央并填充TIP3P水分子与10个Na+中 和的体系电荷，并以Steepest算法对整个体系进行能量优化。随后在等温等压体系(NPT)体 系下，限制蛋白骨架重原子的运动，在500ps内缓慢地从0K升温至303K，充分平衡酶周围的 水分子，体系的温度和压力偶联使用Berendsen算法，得到预平衡模型；
[0066] (3)将预平衡模型在2ns内从303K缓慢地升温至目标温度，然后在等温等压体系 (NPT)下进行模拟。范德华力截断半径、静电相互作用力截断半径、邻近搜索截断半径均根 据AMBER力场推荐，设置为0.9nm。所有的键长采用LINCS算法约束，长程静电相互作用力采 用PME算法，计算步长设置为2f s;
[0067] (4)对分子模拟所得轨迹进行分析：
[0068] 1)通过波动性分析解折叠的关键区域:通过Gromacs自带的分析工具gmx_rmsf计 算C。原子的均方根涨幅随解折叠进程的变化，分析骨架氨基酸均方根涨幅曲线，结果如图1 所示，在450K温度下底部环形区1〇〇ρ4(氨基酸207位至221位)浮动最大，该区域为容易运动 的结构;在473K温度下N端(氨基酸1位至13位)和底部环形区1 〇〇ρ4(氨基酸207位至221位） 的均方根涨幅大幅增加，该区域可能是脂肪酶发生初步解折叠的区域;在500K温度下脂肪 酶整体的均方根涨幅显著提升，主要集中在β折叠头部区域(氨基酸30位至90位）、β折叠尾 巴区域、盖子1区域(氨基酸91位至130位）。
[0069] 2)分析湿熔球态形成的关键区域:根据波动性分析的结果，针对Lip2脂肪酶主要 的三个波动区域进行熔球态形成分析，通过在VMD中统计氨基酸4A范围内的水分子数目随 解折叠进程的变化，分析熔球态形成的关键区域和位点。结果如图2所示，水分子进入Lip2 脂肪酶的主要区域为β折叠尾巴区域。并且通过进一步将β折叠尾巴区域分为4个亚区进行 水分子数目统计，其中亚区1(氨基酸251位至279号)包含盖子2铰链区、β9折叠;亚区2(氨基 酸221位至250号）包含β7、8折叠和之间的无规卷曲1οορ3;亚区3(氨基酸194位至220号）包 括底部环形区1〇〇ρ4与α5螺旋;亚区4(氨基酸180位至193号)包含β6折叠与α5螺旋间的无规 卷曲，结果显示亚区2周围的水分子数目在解折叠过程中增长最为明显，是湿熔球态形成的 主要区域。
[0070] 3)分析湿恪球态形成的过程:通过Gromacs自带的分析工具gmx_rmsf计算β折叠尾 巴区域中每段结构的Ca原子RMSD随时间的变化曲线，以分析熔球态形成的关键结构的波动 与步骤。结果如图3所示，在湿恪球态形成的时间区间内1οορ3、1οορ4、β7和β8折叠的RMSD急 剧上升，暗示着这些区域的运动与恪球态形成直接相关，并且β9折叠与1οορ4的运动发生在 早期，可能是熔球态形成的诱因，但β6折叠在解折叠运动中所受影响不大。
[0071]综上所述方法分析出脂肪酶在热变性的过程中解折叠和湿熔球态形成的关键区 域为β折叠尾巴区域，其中湿熔球态的形成可分为2个步骤，初级步骤为脂肪酶的底部无规 卷曲1〇〇ρ4(氨基酸207位至221位)与Ν端区域(氨基酸1位至13位)的解离，关键步骤为β7、β8 折叠之间的无规卷曲区域(氨基酸228位至246位)结构的解离。
[0072] (7)通过Disulfide by Design软件筛选存在于底部无规卷曲1οορ4(氨基酸207位 至221位)与Ν端区域(氨基酸1位至13位）、β7折叠至β8折叠之间的无规卷曲1〇〇ρ3区域(氨基 酸228位至246位)与周围之间的所有二硫键潜在位点。筛选结果显示，筛选出六对潜在的二 硫键位点，分别为 32-210、58-214、514-216、5191-241、5198-242、54-266。具体信息如表1所 不。
[0073] 表1二硫键筛选位点
[0075]其中，6段氨基酸序列如下所示：
[0088] 编码上述序列的脂肪酶的核苷酸序列，如下所示：
[0089] S2-210：
[0090]
[0101 ]实施例2脂肪酶突变体表达质粒的构建
[0102] 以耶氏解脂酵母脂肪酶2序列（Genbank ID:AJ012632.1)为目的片段，以EcoRI与 Notl为限制性酶切位点，全基因合成并构建pPICZaA-Lip2,通过两次连续的反向突变PCR方 法引入二硫键突变，点突变所用引物如表2所示。
[0103] 表2突变引物汇总
[0106] 注:斜线加粗处为突变位点
[0107] PCR 扩增条件为：941€2111丨11;94°(：1〇8、661€3〇8、68 1€51^11，10个循环。反应体系如下 表3所示。
[0108] 表3 PCR反应体系
[0110]扩增产物以Dnpl酶消化模板，在琼脂糖凝胶电泳检测突变条带大小后，用T4连接 酶连接环化过夜，随后使用热激法将突变质粒转入T0P10大肠杆菌感受态细胞，并涂布于 LLB+Ze〇Cin(Ze〇Cin浓度为25μg/ml)平板37°C过夜培养，挑选阳性转化子进行质粒的测序。
[0111] 实施例3:线性化质粒电转化毕赤酵母、转化子筛选与产酶筛选
[0112] 将测序正确的阳性转化子于LLB液体培养基过夜扩培后，提取质粒，以Pmel线性化 处理并纯化回收，以总量为5yg的质粒线性化产物与X33毕赤酵母感受态混合电击转化。毕 赤酵母感受态制备参照Invitrogen公司操作手册。电转程序按照Bio-Rad公司推荐参数设 置。
[0113] 电转完毕立刻加入lmL lmol/L山梨醇溶液，将菌液在30°C孵育复苏1小时后，均匀 涂布于￥?03+26〇〇;[11(2600;[11浓度为20(^8/1111)抗性平板上筛选 ;培养3天后，将生长出的单 菌落挑取至BMMY-罗丹明B平板上进行产酶筛选，每12小时添加100yL甲醇诱导;诱导3天后 观察脂肪酶水解圈，筛选出正常产酶的基因工程菌。如图4所示，所有的突变体脂肪酶均产 生了油脂水解圈，二硫键的引入没有影响工程菌的产酶。
[0114] 实施例4:工程菌摇床发酵
[0115]参考Invitrogen公司毕赤酵母表达试剂盒操作手册并稍作修改，修改内容如下： 把工程菌株单菌落接种到2mL YPDS-Zeocin(Zeocin浓度为200yg/mL)液体培养基中进行纯 化培养过夜，离心将细胞以BMGY液体培养基重悬并过夜培养，再接种至30mL BMMY液体培养 基，以25°C、300r/min培养96小时，每天补充甲醇至终浓度为1 %。
[0116] 实施例5:脂肪酶的分离和纯化
[0117] 1)超滤浓缩
[0118] 将100mL发酵液于4°C5000r/min离心5分钟后吸取上清液并用10kDa超滤管于 5000r/min 4°C离心50min，收集浓缩酶液。
[0119] 2) -步法镍柱纯化
[0120] ①用5mL的含咪唑10mM的Binding Buffer平衡镍柱，充分除去残留的乙醇；
[0121 ] ②浓缩酶液与含咪唑120mM的Binding Buffer按照1:1比例混合后添加至镍柱中 结合；
[0122] ③用15mL含咪唑60mM的Washing buffer充分洗脱杂蛋白；
[0123] ④用15mL含咪唑300mM的Elution Buffer洗脱目标蛋白；
[0124] ⑤纯化酶液按照上述条件超滤浓缩；
[0125] 得到纯化的突变脂肪酶，最后以还原性SDS-PAGE垂直电泳检测酶纯度，纯度均在 99%以上，结果如图5所示。
[0126] 实施例6:二硫键测定
[0127] 1)将纯化蛋白样品以含8M尿素的Tris-Gly缓冲液(50mM，pH=8 · 00)稀释至0 · 5mg/ mL。混合均匀后，分装2mL至两支4mL离心管中，A管用于测定游离巯基的浓度，B管用于测定 总半胱氨酸含量，并在B管中加入50μ?3-巯基乙醇还原分子内的二硫键，于37 °C保温1小时。
[0128] 2)在B管中加入2mL 30%(w/v)三氯乙酸溶液，37°C水浴1小时，充分沉淀蛋白，随 后以12000r/min离心30min，弃上清，沉淀使用2mL 30%三氯乙酸溶液重悬润洗并再次离 心，晾干5min，挥发残余的β-巯基乙醇。最后用2mL含8M尿素的Tris-Gly缓冲液（50mM，pH = 8.00)重新溶解沉淀。每管吸取lmL做为空白组，并在剩余反应液中加入10yL DTNB(4mg/mL) 显色。孵育30min后以12000r/min离心10min，除去沉淀，吸取上清，并稀释一定的倍数，于 412nm测定吸光值，保证其吸光值在0.2-0.8之间。二硫键数量由以下公式计算：
[0129] SH=73.53XA4i2XD/V----------------(公式 1)
[0130] SH3 = SH1-SH2----------------(公式 2)
[0131] SS = N/2X(SH3/SH1)----------------(公式 3)
[0132] 其中SH为巯基浓度(ymol/mL)，SH3为成键半胱氨酸巯基浓度(ymol/mL)，SH2为游 离半胱氨酸疏基浓度(μπιο 1 /mL)，SH1为总半胱氨酸疏基浓度(μπιο 1 /mL)，A4i2为除去空白样 品后的吸光光度值，D为稀释倍数，V为溶液体积(mL)，SS为蛋白质内部二硫键的数目。
[0133] 经过二硫键测定，结果如表4所示，除S198-242突变脂肪酶未形成二硫键，其余新 引入的二硫键均正确成键。
[0134] 表4二硫键数目测定
[0136] 实施例7:脂肪酶突变体的热稳定性与催化性质的测定。
[0137] 利用荧光定量PCR仪，按照Protein Thermal Shift试剂盒推荐反应程序测定突变 脂肪酶的1值，结果如表5所示，除了 S4-266突变脂肪酶，其余引入的二硫键均大幅度提升 了脂肪酶的Tm，引入的^硫键属于热稳定型^硫键，基于解折置t吴型设计^?硫键大大提升 了设计的成功率。
[0138] 表5 DSF测定结果
[0140] T5Q测定方法：以0.1mg/mL的纯化蛋白溶液在PCR仪中精确保温15min，以50mM含 40mM的p-NPP的Tris-Hcl缓冲液为反应体系(pH=7.50)，精确反应lOmin后，添加20%(w/v) 的三氯乙酸终止反应5min，在20 % (w/v)的碳酸钠溶液显色，测定41 Onm处吸光值，计算不同 温度保温后，脂肪酶残余的相对酶活。
[0141 ] t1/2测定方法：以0. lmg/mL的纯化蛋白溶液在PCR仪中50 °C精确保温0、5、10、15、 30、45、60min，以50mM含40mM的p-NPP的Tris-Hcl缓冲液为反应体系（pH= 7.50)，精确反应 10min后，添加20% (w/v)的三氯乙酸终止反应，20% (w/v)碳酸钠溶液显色，测定410nm处吸 光值，计算不同温度保温后，脂肪酶残余的相对酶活。
[0142] 最适反应温度测定方法:将0.1mg/mL的纯化蛋白溶液加入在30、35、40、45、50、55 °C条件下预热的50mM含40mM的p-NPP的Tris-Hcl缓冲液(pH = 7.50)，精确反应10min后，添 加20% (w/v)的三氯乙酸终止反应，20% (w/v)碳酸钠溶液显色，测定410nm处吸光值，计算 不同温度下脂肪酶的相对酶活。以上结果如表6所示：
[0143] 表6热稳定性测定结果
[0146] 可见，S2-210、S8-214、S14-216和S191-241突变脂肪酶的热稳定性均显著地提高， 其中S2-210在50°C下具有较长的半衰期，并且除S191-241突变脂肪酶外，其余突变脂肪酶 的最适反应温度均提升至40°C。
[0147] 实施例8:尿素解折叠检测湿熔球态的形成。
[0148] 将O.lmg/mL纯化蛋白溶液与含0、0·25、0·5、0·75、1、1·5、2、2·5、3、4、5、6Μ 尿素的 Tris-Hcl缓冲液(50mM，ρΗ = 7.50)均匀混合，常温下静止12个小时，使体系充分平衡后，加 入bis-ANS荧光剂（终浓度为50μΜ)，常温结合1小时后，加入酶标板，在350nm波长激发，以 492nm波长采集荧光信号，以荧光信号强度对尿素浓度作图。结果如图6所示，所有突变脂肪 酶需要更高的尿素浓度方能完全进入湿熔球态，其中S2-210突变脂肪酶进入湿熔球态所需 要的尿素浓度最高，为2.5M，在脂肪酶底部无规卷曲1 〇〇ρ4(氨基酸207位至221位)与N端区 域(氨基酸1位至13位)之间引入二硫键突变，均抑制了湿熔球态的形成。并且S191-241突变 脂肪酶的荧光峰值显著低于其他突变脂肪酶，二硫键显著抑制了湿熔球态进一步地形成， 证实β7、β8折叠之间的无规卷曲区域1 〇〇ρ3(氨基酸228位至246位）的解离的为湿熔球态形 成的关键步骤，与分子模拟结果一致。
[0149] 实施例9:非还原性SDS-PAGE检测脂肪酶共价聚沉
[0150] 将50yL 0.5mg/mL纯化蛋白溶液在PCR仪上50°C精确保温5、10、15、30、45、60min。 以非还原性SDS-PAGE电泳检测脂肪酶的共价聚沉。结果如图7所示，Lip2脂肪酶在50°C处理 5min后，发生了严重聚沉，共价聚体分子量集中在75kDa、1 lOkDa和140kDa左右，随着热处理 的时间增加，单体与二聚体均向更高聚体的转变。在引入二硫键突变后，所有的突变脂肪酶 的共价聚沉均受到了抑制，其中32-210、38-214、314-216、3191-241突变脂肪酶的聚沉抗性 依次递减，其中S2-210、S8-214突变脂肪酶的低聚聚沉主要发生在10~15min的热处理期 间，随后以高聚体聚沉为主。S14-216突变脂肪酶、S191-241突变脂肪酶和Lip2脂肪酶的低 聚聚沉主要发生在热处理的前5min内，并且S191-241二硫键对1οορ3结构的限制并没有完 全抑制共价聚沉，因此的折叠与1οορ4的解离是Cys244游离半胱氨酸暴露的关键原因。
[0151] 上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的 限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化， 均应为等效的置换方式，都包含在本发明的保护范围之内。
【主权项】
1. 一种热稳定的脂肪酶，其特征在于:是在耶氏解脂假丝酵母脂肪酶2的解折叠和湿熔 球态形成的关键区域引入二硫键突变； 所述的关键区域为耶氏解脂假丝酵母脂肪酶2的第180位~280位氨基酸序列。2. 根据权利要求1所述热稳定的脂肪酶，其特征在于:所述热稳定的脂肪酶为S2-210脂 肪酶、S8-214脂肪酶、S14-216脂肪酶或S191-241脂肪酶； S2-210脂肪酶的氨基酸序列如SEQ ID N0.1所示； S8-214脂肪酶的氨基酸序列如SEQ ID NO.2所示； S14-216脂肪酶的氨基酸序列如SEQ ID NO.3所示； S191-241脂肪酶的氨基酸序列如SEQIDN0.4所示。3. 编码权利要求2所述热稳定的脂肪酶的核苷酸序列，其特征在于： 所述的S2-210脂肪酶的编码核苷酸序列如SEQ ID NO.5所示； 所述的S8-214脂肪酶的编码核苷酸序列如SEQ ID NO.6所示； 所述的S14-216脂肪酶的编码核苷酸序列如SEQ ID NO.7所示； 所述的S191-241脂肪酶的编码核苷酸序列如SEQ ID NO.8所示。4. 权利要求1所述热稳定的脂肪酶的制备方法，其特征在于包括如下步骤： (1) 通过Gromacs分子动力学软件模拟耶氏解脂酵母脂肪酶2的解折叠过程，通过对所 得轨迹的波动性分析与蛋白表面水分子的统计，分析出脂肪酶解折叠和湿熔球态形成的关 键区域，并进一步对β折叠尾巴区域分区，分析湿熔球态形成的步骤，最后使用Disulfide by Design软件算法筛选潜在二硫键突变位点； (2) 通过反向PCR突变，将所选出的氨基酸位点突变为半胱氨酸，并转入大肠杆菌工程 菌中，进行扩增培养与质粒提取、测序； (3) 将测序正确的突变质粒以Pme I限制性内切酶线性化处理，并电击转化入感受态毕 赤酵母X33中，进一步通过博来霉素抗性平板筛选和BMMY-罗丹明B产酶平板筛选，得到对应 的突变工程菌； (4) 将突变工程菌在YPD液体培养基中进行扩繁培养后转接至BMGY液体培养基进行去 抑制培养，最后接种BMMY液体培养基进行发酵，将菌液离心获取上清粗酶液； (5) 使用超滤管将粗酶液进行超滤浓缩后，使用镍柱一步法纯化，分离出带组氨酸标签 的脂肪酶蛋白，并以还原性SDS-PAGE检测蛋白纯度，获取纯化后的脂肪酶； (6) 通过DTNB法测定蛋白质中的游离巯基浓度和总巯基浓度，计算蛋白分子内的二硫 键数目，检测引入的二硫键突变是否成键； (7) 通过p-NPP比色法和DSF荧光检测测定突变脂肪酶的热稳定性指标:保温15min后残 余50%活性的温度、50°C下的半衰期、最适反应温度以及蛋白熔解温度； (8) 通过不同尿素诱导解折叠测定湿熔球态的形成所需要的浓度，确定突变脂肪酶动 力学稳定性大小； (9) 将突变脂肪酶在50°C下保温不同的时间，以非还原性SDS-PAGE检测突变脂肪酶的 抗聚沉能力大小，得到热稳定的脂肪酶。5. 根据权利要求4所述热稳定的脂肪酶的制备方法，其特征在于：步骤（1)中所述的脂 肪酶解折叠和湿熔球态形成的关键区域为耶氏解脂假丝酵母脂肪酶2的第180位~280位氨 基酸序列。6. 根据权利要求4所述热稳定的脂肪酶的制备方法，其特征在于：步骤（I)中所述的湿 熔球态形成的步骤包括初级步骤和关键步骤； 所述的初级步骤为耶氏解脂假丝酵母脂肪酶2第207~221位氨基酸序列与第1~13位 氨基酸序列的解离； 所述的关键步骤是耶氏解脂假丝酵母脂肪酶2第228~246位氨基酸序列的解离。7. 权利要求1所述的热稳定型脂肪酶在工业中的应用。
【文档编号】C12N9/20GK105950585SQ201610279266
【公开日】2016年9月21日
【申请日】2016年4月29日
【发明人】管武太, 吴炜坤, 李力浪
【申请人】华南农业大学