NMR方法解析蛋白质结构

NMR方法解析蛋白质结构冯银刚北京核磁共振中心()北京大学化学与分子工程学院[e-mail]@gmail.com2006.4.12蛋白质结构层次氨基酸通过肽键形成的生物高分子一级结构、二级结构、三级结构、四级结构•肽键具有双键性质而不能任意旋转•主链可旋转的二面角，20种常见的氨基酸残基通常NMR中一个自旋系统是指一个氨基酸残基上的所有原子核磁共振方法解析蛋白质溶液结构测定原子（氢原子）之间的距离信息和其他约束信息，得到空间结构模型化学结构（氨基酸序列，即一级结构）已知，测定空间结构（三级结构，四级结构）适用范围——分子量分子量受限：谱峰重叠，分子增大造成横向弛豫时间减少（线宽增加）同核二维1H谱：80氨基酸异核多维（二维、三维）：常规三共振（1H,15N,13C）200氨基酸氘标记，特殊标记已报道有700氨基酸同核样品可以是天然提取、化学合成或分子生物学方法培养制备标记样品通常只能通过分子生物学方法培养制备常规样品浓度：毫摩尔每升500uL10kDa5毫克每样品均一没有聚合水溶液（10%重水锁场）pH7.5常用缓冲体系：Tris-HCl,磷酸盐，醋酸盐稳定性：室温下大于一周NaN3蛋白酶抑制剂主链HN侧链HN芳环HA脂肪族氨基酸侧链甲基H2ODSSDTT蛋白质的一维谱同核方法指认方法不同：2D1H-1H两种样品：H2O,D2O实验TOCSY,COSY,NOESY自旋系统指认，序列指认指纹区：HN-HA每个残基一个峰同核方法指认基于相邻序列间的NOE三共振方法—以杨树谷氧还蛋白C1为例实验之前——蛋白质基本性质杨树谷氧还蛋白C1(Grx-C1)分子量：117个氨基酸残基，12.5kDa等电点：8.49同源性：~35%功能：谷胱甘肽依赖的氧化还原酶稳定、无聚合MASKQELDAALKKAKELASSAPVVVFSKTYCGYCNRVKQLLTQVGASYKVVELDELSDGSQLQSALAHWTGRGTVPNVFIGGKQIGGCDTVVEKHQRNELLPLLQDAAATAKNPAQL样品制备——分子生物学方法基因克隆：cDNA质粒载体蛋白质表达纯化大肠杆菌E.coliBL21(DE3)离子交换色谱，凝胶过滤纯度：SDS-PAGE电泳上为单一条带同位素标记使用M9培养基：无机盐，葡萄糖15NNH4Cl,13C葡萄糖样品~1mM蛋白质磷酸钾缓冲液，pH6.4，90%H2O/10%D2O0.01%DSS（化学位移校准）40mMDTT，0.01%NaN3，蛋白酶抑制剂数据采集谱仪Bruker500MHz,600MHz三共振探头采集用于解析一个蛋白质所需要的全部谱图时间约1~2月数据处理处理的内容：线性预测－加窗函数－充零－傅立叶变换－相位校正－基线校正在线与离线处理二维谱：每个约几兆三维谱：每个几十兆到几百兆nmrPipe多维谱的处理：每一维都要分别进行傅立叶变换化学位移校准重要性1H使用DSS(2,2-dimethylsilapentane-5-sulfonicacidsodiumsalt)异核：间接校准以DSS的0ppm处的频率乘以换算因子得到异核的0ppm的频率NMR初步鉴定蛋白质折叠程度H2ODSSDTTNMR初步鉴定1H-15NHSQC分散程度：折叠程度确定谱宽指认基础：每个氨基酸（除脯氨酸外）都在HSQC上有一个峰三维实验（双共振、三共振）减少谱峰堆积利用异核之间较大的耦合常数基于J耦合进行共振指认，减少对空间构象的依赖三维实验绝大多数异核多维实验可以类比于二维COSY,TOCSY,NOESY及其组合主链指认实验侧链指认实验NOE实验123HN---NH---Ca(i,i-1)HNCaNH化学位移指认化学位移：每个原子的身份证主链指认序列连接氨基酸类型判断主链指认HNCACB-CBCA(CO)NH化学位移指认侧链指认3DTOCSY或COSYH(C)CH-TOCSY芳环(HB)CB(CGCD)HD(HB)CB(CGCDCE)HE通过NOE确认指认主链侧链用于Grx-C1指认的实验主链指认：2D1H-15NHSQC,3DHNCA,HNCACB,CBCA(CO)NH,HNCO,HN(CA)CO~10天侧链指认：2D1H-13CHSQC,3DHBHA(CBCA)(CO)NH,(H)C(CO)NH,H(C)(CO)NH,1H-15NTOCSY-HSQCHCCH-TOCSY,CCH-TOCSY,HCCH-COSY,CCH-COSY,2D(HB)CB(CGCD)HD,(HB)CB(CGCDCE)HE~20天NOE：3D15NNOESY-HSQC,13CNOESY-HSQC(foraliphaticregion),13CNOESY-HSQC(foraromaticregion)~10天Grx-C1指认结果1H,15N,13C指认率：98%1H-15NHSQC上所有谱峰得到了指认J.Biomol.NMR2005,31:263-264BMRB6410由化学位移得到二级结构的信息二级结构CSI：化学位移与无规卷曲的多肽化学位移之差TALOS：基于数据库比对预测二面角可用于结构计算和分析结构计算基本方法由实验得到各种构象约束信息：距离，二面角约束信息是不完备的约束信息是不精确的计算满足这些约束条件的构象距离几何(DistanceGeometry)约束条件下的分子动力学模拟(RestrainedMolecularDynamicsSimulation)约束信息距离约束1H-1HNOE氢键二面角约束主链侧链肽键：反式－顺式其他手性L-氨基酸约束生成NOE指认Unique：只有一种可能Ambiguous：有多种可能转化为距离实验通常可检测5Å距离是不精确的：距离范围距离是大量的精确结构•其中一种可能是正确的•多种可能都是正确的（谱峰重叠）二面角约束a-helix:f=-100to-20,=-70to-10b-sheet:f=-180to-60,=20to180二面角预测：CSI或TALOS实验：HNHA,HNHB等,对特定的二级结构有一定的分布拉氏图(Ramachandranplot)最佳区次佳区一般区不允许区二面角约束二级结构判断：特征的NOE信号氢键约束α-HelixParallelβ-StrandsAnti-parallelβ-StrandsH----O2.0ÅN----O3.0ÅN-H----O•根据二级结构或实验判断氢键分子动力学模拟（结构计算）模拟分子随机运动，使其达到能量的最小值约束条件转化为MD中的能量项Vtotal=Vbond+Vangle+Vdihedr+Vvdw+Vcoulomb+VNMR模拟退火：克服局部的能量最小点计算多个结构，取能量较低的若干结构作为结果Ensemble(NMR信息的不完备和不精确)KNOE(rij-rij1)2ifrijrijlKNOE(rij-riju)2ifrijriju0ifrijlrijrijuVNOE=结构计算结构检查违约分析：查错化学位移指认NOE指认其他效应对NOE的影响结构评价反复计算若干循环违约：结构计算是计算能量最低的结果，因此最后的结构可能不完全符合每一个约束。违约产生原因是约束中存在错误，造成约束之间或约束和化学结构之间存在矛盾。结构评价能量二级结构拉氏图尽可能少的氨基酸残基处于不允许区RMSD（均方根偏差）表明结构的收敛程度自动化结构计算NOE自动指认CANDID/CYANA只需提供指认的化学位移列表和NOE谱峰列表自动进行NOE指认自动通过7个结构计算循环（100个结构取二十个），逐步优化指认结果自动计算结果正确性依赖于原子化学位移指认的比例和正确性（至少90%）SANE依赖初始结构的自动NOE指认Grx-C1的结构计算CYANA结构初步优化力场相对简单运行速度快无需初始结构结构相对较为粗略Amber结构精修具有更精细的力场参数，使用溶剂化模型（或显式加溶剂）从而获得更加合理的局部构象需要整体折叠正确的初始结构运算量大，速度慢Grx-C1的结构计算CANDID/CYANA得到初始结构（全自动）2CPU,~8小时SANE-CYANA循环，进行初步优化（半自动）手工分析违约和未指认的NOE每个循环2CPU～1小时～20-40个循环SANE-AMBER循环，进行结构精修（半自动）每个循环20CPU～15小时～10~30个循环结构计算结果PDB1Z7P(ensemble),1Z7R(mean)结构计算结果约束统计NOE：4845二面角：160氢键：47手性：287违约状况距离无0.2Å二面角无结构计算结果结构评价Mostfavoredregions(%)88.8Additionallyallowedregions(%)10.7Generouslyallowedregions(%)0.5Disallowedregions(%)0.0RMSDAllresiduesRegularsecondarystructureBackboneheavyatoms0.880.32Allheavyatoms1.130.68三共振方法所需时间提取样品——2～4周数据采集——4～8周数据处理——1～4小时每实验共振指认——2～4周结构计算——2～4周数据处理通常在数据采集后立即进行，在总时间中可忽略不计数据采集与共振指认的时间可部分重叠一些其他因素水峰压制小分子干扰对异核实验影响不大温度：控温装置分辨率、信噪比与时间分辨率：间接维点数实验时间与间接维点数成正比信噪比与累加次数的平方根成正比实验时间与累加次数成正比最新的进展谱仪：超低温探头提高灵敏度2－4倍TROSY（横向弛豫优化谱）在700MHz以上高场谱仪中可以显著减小线宽，从而可用于测定更大分子量的蛋白质RDC（残余偶极耦合）得到化学键的空间相对取向信息自动化方法显著加速结构计算进程使用的软件NMRPipe数据处理NMRView指认分析CYANA结构计算500EuroTALOS基于化学位移预测主链二面角(NMRPipe的一部分)SANE基于结构的NOE自动指认JBiomolNMR,200119(4)321-9Amber分子动力学模拟。用于结构优化$400PROCHECK-NMR结构分析与评价~roman/procheck_nmr/procheck_nmr.htmlMOLMOL结构分析与绘图其他软件数据处理Felix$????AZARAFreePROSA(Free?)