-
甘油脂的从头生物合成途径(de novo glycerolipid biosynthesis)是细胞中最基本的代谢过程。3-磷酸甘油酰基转移酶(glycerol-3-phosphate acyltransferase, GPAT)催化甘油脂从头合成的初始步骤,生成的溶血磷脂酸(LPA)在LPA酰基转移酶(lysophosphatidic acid acyltransferase, LPAAT)的作用下转化为磷脂酸(PA)[1−6]。PA是调节真核生物中多种细胞过程的重要信号分子,也是极性甘油磷脂与中性三酰甘油(TAG)的生物合成前体。在磷酸酶的催化下PA转化为二酰甘油,后者可经脂酰辅酶A-依赖型二酰甘油酰基转移酶(acyl-CoA: diacylglycerol acyltransferease, DGAT)和(或)磷脂-依赖型二酰甘油酰基转移酶(phospholipid-dependent diacylglycerol acyltransferase, PDAT)作用生成TAG[7−10]。TAG合成能力是油料作物的关键性状,同时与人类肥胖症等疾病密切相关。因而,运用遗传或化学遗传方法操控TAG生物合成,提高油料作物含油量,降低与肥胖症相关联的人类疾病,具有重要实践意义[11−16]。
脂酰基转移酶在TAG生物合成过程中发挥重要作用,但目前对于这类酶的结构与功能的内在关系知之甚少,参与第1步酰化反应的GPAT亦不例外,仅有少量关于其结构与功能关系的报道[17−18]。已知GPAT、LPAAT、磷酸二羟丙酮酰基转移酶(dihydroxyacetone-phosphate acyltransferase, DHAPAT)等脂酰基转移酶均含4个高度保守的结构域,结构域Ⅰ的组氨酸(H)、天冬氨酸(D),结构域Ⅲ的甘氨酸(G)和结构域Ⅳ的脯氨酸(P)是GPAT催化所必需的;而结构域Ⅱ的精氨酸(R)与结构域Ⅲ的谷氨酸(E)在结合底物3-磷酸甘油中起作用[18−20]。迄今为止,对于保守结构域外的其他氨基酸残基在酰基转移酶活性调控中的作用,及与保守结构域中的氨基酸残基存在的潜在互作关系的了解非常有限。
GPAT9位于植物细胞的内质网,参与膜脂和TAG的生物合成,其功能缺失会导致种子发育异常、油脂合成减弱[20−22]。本实验室前期的酵母遗传互补研究也显示:油菜Brassica napus BnGPAT9的异源表达能够恢复酵母条件致死型双敲除突变体(ZAFU1)因GPAT酶活性缺失引起的生长缺陷。然而,拟南芥Arabidopsis thaliana AtGPAT9却不具备这种互补能力[23−25],尽管AtGPAT9与BnGPAT9的进化关系密切[26−27],两者氨基酸序列的一致性高达 94.1%,且两者在4个酰基转移酶保守结构域的氨基酸残基完全一致。因此,可以假设保守结构域之外的某些氨基酸残基对GPAT9的活性起着重要的调节作用。本研究充分利用AtGPAT9和BnGPAT9在酵母异源系统中表现出的不同性质,并结合定点突变与酵母遗传互补技术,剖析单个和多个氨基酸残基改变对GPAT9酶活性的影响,鉴定新的关键活性位点,以深化对脂酰基转移酶结构与功能内在关系的认知,为酰基转移酶的分子改造与结构优化、真核生物中TAG合成途径的改良以及全新TAG从头合成途径的构建提供理论基础。
-
通过Vector NTI 11.5.4软件对AtGPAT9 (Genebank 登录号: ACT32031.1)和BnGPAT9 (Genebank 登录号: ANV28166.1)的氨基酸序列进行比对。使用TMHMM 2.0和Protter进行跨膜结构域和蛋白质拓扑异构模型预测[28-29]。由I-TASSER预测三维结构[30]。
-
将AtGPAT9和BnGPAT9编码序列分别通过BamH I/XhoⅠ和BamH I/EcoR I双酶切位点克隆至pMD19-T载体,得到新的质粒pMD19-T-AtGPAT9和pMD19-T-BnGPAT9;以之为模板,对AtGPAT9和BnGPAT9进行定点突变。具体方法如下:利用包含突变位点的引物(表1),PCR扩增整个质粒;质粒DNA经Dpn I消化与纯化后,转入大肠埃希菌Escherichia coli并进行测序分析;序列正确的质粒经BamH I和Xho Ⅰ双酶切后,连接到经相同酶切处理的pYES2-yADH1-Kan V2酵母表达载体,并对产生的重组质粒再次进行DNA测序分析,以确保突变位点的正确性。需要说明的是,选择pMD19-T质粒作为定点突变过程中的中间质粒,而非直接在pYES2-yADH1-Kan V2质粒上进行GPAT9基因的定点突变,是因为后者DNA长度(6 998 bp)是前者(2 660 bp)的2.6倍,采取这样的策略可以降低因PCR扩增时间延长导致潜在的错误碱基出现频率。
突变位点 引物序列 (5′→3′) 突变位点 引物序列 (5′→3′) BnGPAT9(R40S) AGCCTCGTGGCAAGCTCAGCCTGCGTGATTTGCTAGACAT AtGPAT9(N119H) TTTCATTGTTTATCCCTGTACACGCGTTGCTGAAAGGTCAAG BnGPAT9(W85Y) TCTACTTGTTTCCTTTATACTGCTGTGGTGTTGTTGTTAG AtGPAT9(D230N) TTGTAGCAAAAAAGTTAAGGAACCATGTCCAAGGAGCTGAC BnGPAT9(C87F) TTGTTTCCTTTATGGTGCTTTGGTGTTGTTGTTAGATACT AtGPAT9(A235T) TAAGGGACCATGTCCAAGGAACTGACAGTAATCCTCTTCTC BnGPAT9(I102F) TTCTCTTTCCCTTGAGGTGCTTCACTTTAGCTTTTGGATG AtGPAT9(S237N) ACCATGTCCAAGGAGCTGACAATAATCCTCTTCTCATATTTCC BnGPAT9(F109I) CATCACTTTAGCTTTTGGATGGATTATTTTCCTTTCAACG AtGPAT9(D230N/A235T) TTGTAGCAAAAAAGTTAAGGAACCATGTCCAAGGAACTGAC BnGPAT9(T114L) TGGTTTATTTTCCTTTCATTGTTTATCCCTGTACACTCTC AtGPAT9(A235T/ S237N) TAAGGGACCATGTCCAAGGAACTGACAATAATCCTCTTCTCATATTTC BnGPAT9(H119N) TTCAACGTTTATCCCTGTAAATTCTCTCCTGAAAGGTCAG AtGPAT9(D230N/A235T/
S237N)TTGTAGCAAAAAAGTTAAGGAACCATGTCCAAGGAACTGAC BnGPAT9(N230D) GTAGCAAGAAAGTTAAGGGACCATGTTCAAGGAACTGACA AtGPAT9(G332A) CATAAGGCCCGGTGAAACAGCAATTGAATTTGCAGAGAGGG BnGPAT9(T235A) TAAGGAACCATGTTCAAGGAGCTGACAATAACCCTCTTCT AtGPAT9(L335H) GGTCAGAGACATGATATCTCATCGGGCGGGTCTCAAAAAGG BnGPAT9(N237S) CATGTTCAAGGAACTGACAGTAACCCTCTTCTTATATTTC AtGPAT9(P355S) TGAAGTATTCGAGACCAAGCTCCAAGCATAGTGAACGCAAG BnGPAT9(A322G) AAGGCCTGGTGAAACAGGAATTGAGTTTGCAGAGAGGGTC AtGPAT9(T10A) GTACGGCAGGGAGGCTCGTGGCTTCAAAATCCGAGCTTGAC AtGPAT9(S40R) ATGAACCTCGCGGCAAGCTCCGCCTGCGTGATTTGCTAGA AtGPAT9(S11A) CGGCAGGGAGGCTCGTGACTGCAAAATCCGAGCTTGACCTC AtGPAT9(Y85W) ATTTACTTATTCCCACTATGGTGCTTTGGGGTTGTTGTTAG AtGPAT9(S13A) GGAGGCTCGTGACTTCAAAAGCCGAGCTTGACCTCGATCAC AtGPAT9(F87C) CTTATTCCCACTATACTGCTGTGGGGTTGTTGTTAGATACT AtGPAT9(S28A) AACATCGAAGATTACCTTCCTGCTGGTTCTTCCATCAATGAAC AtGPAT9(F102I) TCCTCTTTCCCTTGAGGTGCATCACTTTAGCTTTTGGGTGG AtGPAT9(S30A) GAAGATTACCTTCCTTCTGGTGCTTCCATCAATGAACCTCGCG AtGPAT9(I109F) TCACTTTAGCTTTTGGGTGGTTTATTTTCCTTTCATTGTTT AtGPAT9(S31A) GATTACCTTCCTTCTGGTTCTGCCATCAATGAACCTCGCGGCA AtGPAT9(L114T) GGGTGGATTATTTTCCTTTCAACGTTTATCCCTGTAAATGCG Table 1. Sequences of the primers used for site-directed mutagenesis of A. thaliana AtGPAT9 and B. napus BnGPAT9
-
条件致死型酵母双突变体ZAFU1[BY4742, gat1Δgat2Δ+(pGAL1::AtGPAT1 Leu2)] [25, 31],可在半乳糖的培养基上生长,但在葡萄糖培养基上丧失了生长能力。基于菌株ZAFU1建立的酵母遗传互补法对GPAT的鉴定具有很强的专一性[25, 31],本研究运用它鉴定不同氨基酸残基突变对GPAT9活性的影响。
使用基于醋酸锂的标准方法将重组酵母表达质粒导入菌株ZAFU1感受态细胞[25, 31],复苏4 h,分别涂布转化液于以葡萄糖(Glu)或半乳糖(Gal)为碳源,不含尿嘧啶、组氨酸和亮氨酸的培养基(SC-Ura-His-Leu)上, 30 ℃下培养3~5 d。为了准确地比较不同氨基酸残基突变对酶活性的影响,挑取在半乳糖培养基上生长的不同单菌落酵母进行浓度梯度稀释培养实验:从半乳糖培养基上随机挑选生长良好的单菌落至SC-Ura-His-Leu+Gal液体培养基,30 ℃振荡培养1~2 d至光密度[D(600)]为2.000 0~3.000 0,稀释菌液浓度至D(600)为1.000 0、0.200 0、0.040 0、0.008 0和0.001 6,取5 µL接种于SC-Ura-His-Leu+Glu和SC-Ura-His-Leu+Gal固体培养基上,30 ℃培养3~5 d。
-
30 ℃下,将表达不同GPAT9突变基因的ZAFU1细胞在SC-Ura-His-Leu+Gal液体培养基中培养至D(600)为3.000 0~4.000 0,稀释接种于SC-Ura-His-Leu+Glu液体培养基至D(600)为0.100 0,振荡培养并定时记录D(600)。
将平台生长期收获的细胞在真空冷冻干燥机中干燥,提取酵母总脂质[32],点样于硅胶板,通过薄层色谱法分离总脂质。喷洒质量浓度为0.05%樱草黄显色剂,在紫外灯下观察板上的脂质,从硅胶中提取TAG,并通过气相色谱法定量分析TAG含量[32]。
-
早期研究发现:在酵母中异源表达时,拟南芥AtGPAT9与油菜BnGPAT9表现出不同的活性[23-24]。为了找出潜在的关键活性调控位点,对AtGPAT9与BnGPAT9进行了序列比对。AtGPAT9和BnGPAT9均由376个氨基酸组成,两者序列一致性高达94.1%,4个保守的酰基转移酶结构域和C端内质网定位必需的疏水五肽结构域(−ILARL−)的氨基酸残基完全一致[18, 20],且在N端都含多个潜在的磷酸化位点,主要由丝氨酸和苏氨酸组成(图1)。它们之间共有22个不同的氨基酸残基,其中11个具有相似的性质。基于TMHMM和Protter的跨膜结构域预测显示:跨膜区结构具有较高相似性,但3个跨膜区域中有7个不同的氨基酸残基(图1)。另外,基于I-TASSER的三维结构预测发现:AtGPAT9和BnGPAT9在40、109、114、119、230、235、237和322位氨基酸残基的不同,可能会引起两者三维空间结构的差异,因此它们成为本研究重点剖析的位点(图2)。
-
为了明确哪些氨基酸位点对GPAT酶活性起着重要调控作用,据上述AtGPAT9和BnGPAT9中存在的氨基酸残基差异,运用定点突变技术对两者相应位置的氨基酸位点进行相互替代,即将AtGPAT9中的单个或多个氨基酸残基同时替换成与BnGPAT9中相应位置完全相同的氨基酸残基,反之亦然。本研究共构建了58种不同的GPAT9突变基因(表2)。
单个氨基酸残基突变 AtGPAT9氨基酸残基突变组合 BnGPAT9 AtGPAT9 BnR40S** AtT10A AtF102I/S237N AtY85W/D230N BnW85Y AtS11A AtI109F/S237N AtY85W/A235T BnC87F** AtS13A AtD230N/A235T AtY85W/S237N** BnI102F** AtS28A AtD230N/S237N AtY85W/D230N/A235T BnF109I** AtS30A AtA235T/S237N AtY85W/D230N/S237N BnT114L*** AtS31A AtD230N/A235T/S237N AtY85W/A235T/S237N BnH119N AtS40R AtS237N/G322A AtY85W/D230N/A235T/S237N BnN230D*** AtY85W AtY85W/N119H*** AtS40R/Y85W/S237N** BnT235A** AtF87C AtY85W/L114T/N119H* AtN119H/D230N BnN237S* AtF102I AtY85W/N119H/S237N*** AtN119H/A235T BnA322G* AtI109F AtY85W/L114T/N119H/S237N*** AtN119H/S237N*** AtL114T AtY85W/N119H/D230N** AtN119H/D230N/A235T AtN119H* AtY85W/N119H/A235T** AtN119H/D230N/S237N** AtD230N AtN119H/A235T/S237N*** AtA235T AtN119H/D230N/A235T/S237N AtS237N AtG322A AtL335H AtP355S 说明:每种突变以物种的首字母缩写和突变前后的氨基酸残基缩写表示,如AtS40R/S237N代表AtGPAT9的40位由丝氨酸(S)变为精氨酸(R),237位由丝氨酸(S)变为天冬酰胺(N)。*代表基因的异源表达能够恢复酵母双突变体ZAFU1的生长缺陷;*数目代表恢复能力的大小,数目越多,能力越强。 Table 2. Site-directed mutagenesis of amino acid residues at single and multiple sites in AtGPAT9 and BnGPAT9
-
将不同GPAT9突变基因克隆至带有葡萄糖诱导启动子(ADH1)的质粒中,并以空载体与含野生型AtGPAT9的质粒为阴性对照,以含野生型BnGPAT9的质粒为阳性对照,将这些重组质粒导入到ZAFU1菌株中,测定不同GPAT9突变基因恢复ZAFU1在葡萄糖培养基上的生长缺陷能力。根据转化酵母细胞的生长速率,可以比较直观地评估不同突变位点对GPAT9酶活性的影响。结果显示:当ZAFU1突变体在葡萄糖培养基上培养时,W85Y或H119N的单位点替换导致BnGPAT9丧失对突变体生长缺陷的恢复能力(图3A),说明W85和H119是BnGPAT9正常功能所必需的。另外,与野生型BnGPAT9相比,含N237S或A322G突变位点的BnGPAT9对ZAFU1菌株生长的促进作用下降(图3A),表明这2个位点亦参与BnGPAT9活性的调节。相反,其他5个单位点替换(R40S、C87F、I102F、F109I、T235A)对BnGPAT9的活性不产生明显影响(图3A)。特别是N230D和T114L单位点替换增强了BnGPAT9活性,这种上调作用对于BnT114L而言尤为突出。如图4所示:与野生型BnGPAT9相比,含T114L突变位点的BnGPAT9在酵母突变体中的表达能促使细胞生长速率大幅度提高,表现为经2 d培养,表达BnGPAT9 (T114L)的菌落在葡萄糖培养基上生长的数量为野生型BnGPAT9的4倍,且每个单菌落表面积更大。
Figure 4. Effects of the substitution of threonine for leucine at residue 114 (T114L) on BnGPAT9 activity
类似地,对19个不同的AtGPAT9突变基因进行了酵母遗传互补鉴定,其中的6个编码蛋白分别在N端的潜在磷酸化位点发生T10A、S11A、S13A、S28A、S30A和S31A 替换,另外13个分别发生了S40R、Y85W、F87C、F102I、I109F、L114T、N119H、D230N、A235T、S237N、G322A、L335H和P355S替换。结果显示:N端6个潜在磷酸化位点分别替换成中性氨基酸残基并不能改善AtGPAT9活性(图3B)。除了N119H,其他的单位点突变亦对AtGPAT9在酵母异源系统中的活性不产生可见影响(表2)。在酵母菌浓度梯度稀释培养实验中,N119H替换能使AtGPAT9恢复ZAFU1突变体在葡萄糖上的生长缺陷,但这种作用相对较弱(表2,图5D)。
-
基于相邻与非相邻氨基酸残基之间均可能对酶活性产生某种特定的互作效应的假设,根据AtGPAT9和BnGPAT9之间存在的氨基酸差异,进一步构建了28个含有2~4个氨基酸残基替换的AtGPAT9突变酶。与W85、H119和N237位点对BnGPAT9活性产生重要调节作用一致,同步替换AtGPAT9上2或3个相应位点的氨基酸残基(Y85W、N119H、S237N)均能大幅提高AtGPAT9酶活性,表现为携带Y85W/N119H、Y85W/S237N、N119H/S237N、Y85W/N119H/S237N突变的AtGPAT9均能互补菌株ZAFU1的生长缺陷(图5A);不过,单位点突变对AtGPAT9活性的影响十分有限(表2)。进一步调查发现:携带Y85W/N119H或Y85W/N119H/S237N突变组合的AtGPAT9比野生型BnGPAT9更能促进菌株ZAFU1的生长,且与携带T114L突变的BnGPAT9具有相近的作用效果(图3A,图5A)。
为了更好地剖析114位氨基酸性质对GAPT9活性的调节作用以及与其他氨基酸相互作用产生的效应,将双位点突变酶AtGPAT9 (Y85W/N119H)和三位点突变酶AtGPAT9 (Y85W/N119H/S237N)的114位亮氨酸(L114)替换为BnGPAT9相应位置存在的苏氨酸(T)。结果显示:生成的三位点突变酶AtGPAT9 (Y85W/L114T/N119H)的活性明显弱于AtGPAT9 (Y85W/N119H),但四位点突变酶AtGPAT9(Y85W/L114T/N119H/S237N)仍与AtGPAT9 (Y85W/N119H/S237N)的活性相当(图5B)。对这种现象的可能解释是,当2个潜在磷酸化位点即T114和S237同时出现于AtGPAT9 (Y85W/L114T/N119H) (其中含S237)或BnGPAT9 (N237S)(其中含T114)时,GPAT9的磷酸化程度可能加剧,从而抑制酶的活性。因此,推测植物GPAT9活性可能受磷酸化和非磷酸化机制调节,野生型GPAT9中存在的2个潜在磷酸化位点T114和S237可能对酰基转移酶活性产生负面效应。
进一步研究发现:230位氨基酸残基能与85、119位氨基酸残基产生互作效应而影响GPAT活性。虽然D230N本身或与A235T、S237N组合突变未能对AtGPAT9酶活性产生明显影响(表2),但当与N119H、Y85W/N119H或Y85W/S237N结合时, D230N突变能对AtGPAT9活性产生抑制作用。这是因为,与表达含N119H、Y85W/N119H或Y85W/S237N突变位点的酶的菌株相比,表达AtGPAT9 (N119H/D230N)、AtGPAT9 (Y85W/N119H/D230N)或AtGPAT9 (Y85W/D230N/S237N)的菌株ZAFU1在葡萄糖培养基上呈现生长速率明显减弱或不能生长的现象,暗示着D230N的替换不利于AtGPAT9活性(表2,图5B~D)。这一推测得到下述结果的支持,即N230D替换能提高BnGPAT9活性(图3A)。但是AtGPAT9(N119H/D230N/S237N)与AtGPAT9 (N119H/S237N)的活性相当,这说明D230N的负效应依赖于其他氨基酸的互作关系(表2,图5D)。
另外,N119H/D230N/S237N和N119H/A235T/S237N三突变组合均能提高AtGPAT9活性,使之具有拯救ZAFU1生长缺陷的能力,但在前者和后者中分别添加A235T和D230N得到的N119H/D230N/A235T/S237N四突变组合,能使相应蛋白丧失GPAT活性,如野生型AtGPAT9一样,无法恢复ZAFU1的生长缺陷(表2,图5D)。此外,无论是S237N和G322A单或双替换,均不能增强AtGPAT9在酵母中表达时的活性(表2),这与N237S或A322G的替换导致BnGPAT9活性下降现象不一致。由此可见,尽管在237位保留非磷酸化氨基酸(天冬酰胺,N)对BnGPAT9活性有利,但其效应取决于其他位置的氨基酸性质。
综上所述,GPAT9的酶活性受85、114、119、230、237和322位氨基酸残基性质的影响,它们之间存在互作效应;当这些位置的氨基酸残基分别为色氨酸(W)、亮氨酸(L)、组氨酸(H)、天冬氨酸(D)、天冬酰胺(N)和丙氨酸(A)时,酰基转移酶的活性较高。上述多位点突变增强AtGPT9酶活性的事实说明:运用分子设计能有效地改造和优化酰基转移酶的结构,并使之产生新的特性,这对将来人为操控甘油脂的从头合成十分有益。
-
为进一步探究氨基酸位点突变对GPAT9功能的影响,选取了3个活性程度不同的AtGPAT9突变酶,分别含N119H、Y85W/N119H 和Y85W/N119H/S237N突变,使其进行异源表达,测定其表达对酵母细胞中TAG合成的影响;阳性对照为野生型BnGPAT9。当相应酵母细胞的生长进入平台生长期时,提取总脂质,经薄层层析分离后,采用气相色谱法测定TAG含量。
需要指出的是,因不同突变酶活性的差异,相应酵母菌到达平台生长期所需的培养时间不一,且平台生长期时的细胞密度也不完全相同。譬如,表达BnGPAT9、AtGPAT9 (Y85W/N119H)和AtGPAT9 (Y85W/N119H/S237N)的酵母细胞在平台生长期时的细胞密度,以D(600)表示,分别为3.72、4.40、5.32 (图6)。总体而言,酶活性愈大,细胞生长越快、密度越高。
Figure 6. Effects of expression of AtGPAT9 bearing different mutation sites on cell density of the yeast strain ZAFU1 in different periods of growth
脂质分析显示:AtGPAT9 (Y85W/N119H/S237N)的表达使突变体酵母ZAFU1细胞中的TAG含量达0.51%,而BnGPAT9和AtGPAT9 (Y85W/N119H)的表达则分别产生0.35%和0.39%的含油量,前者比后两者分别增加了45.7%和 30.8% (图7)。与AtGPAT9 (N119H)具有相对较低的活性一致,表达此酶的酵母细胞,其TAG含量仅为0.22%,显著低于表达BnGPAT9的酵母细胞(0.35%)。以上结果表明:酵母细胞中TAG的合成能力受GPAT活性的调节,而GPAT活性大小与特定位置的氨基酸性质密切相关,因此通过分子设计优化GPAT9的氨基酸组成将有助于修饰细胞中TAG的合成能力。
Identification of key amino acid residues controlling the activities of glycerol-3-phosphate acyltransferases in Arabidopsis thaliana and Brassica napus
doi: 10.11833/j.issn.2095-0756.20220764
- Received Date: 2022-12-15
- Rev Recd Date: 2023-03-02
- Available Online: 2023-07-13
- Publish Date: 2023-08-20
-
Key words:
- Arabidopsis thaliana /
- Brassica napus /
- glycerol-3-phosphate acyltransferase /
- structure-function relationship /
- site-directed mutagenesis /
- yeast genetic complementation
Abstract:
Citation: | LIN Yixin, CHEN Dandan, LIU Hongbo, KE Xingxing, ZHENG Yueping, ZHENG Zhifu. Identification of key amino acid residues controlling the activities of glycerol-3-phosphate acyltransferases in Arabidopsis thaliana and Brassica napus[J]. Journal of Zhejiang A&F University, 2023, 40(4): 695-706. doi: 10.11833/j.issn.2095-0756.20220764 |