同步辐射晶体学衍射.ppt

1、同步辐射晶体学衍射,董宇辉 BSRF，IHEP，CAS,什么是同步辐射？,接近光速运动的电子或正电子在改变运动方向时放出的电磁辐射，因为是在同步加速器上发现的，所以称为同步辐射。 这种辐射强度高、覆盖的频谱范围广，可以任意选择所需要的波长且连续可调，因此成为一种科学研究的新光源。,同步辐射装置示意图,Booster,储存环,三种发光元件,弯转磁铁,扭摆器(Wiggler),波荡器(undulator),同步辐射的作用,同步辐射的加入，大大提高了结构测定的速度。,PDB(Protein Data Bank),同步辐射促进了蛋白质结构的解析,1980-1990年，获得解析的蛋白质数量大幅度增加；

2、1980左右，同步辐射开始应用到蛋白质结构解析中，1990年各种技术趋于成熟，同步辐射大规模应用于生物大分子结构解析。,SR vs NMR,80%以上的蛋白质结构由同步辐射晶体学衍射方法解析； 大约15%的结构由NMR（核磁共振）方法得到； 同步辐射解析结构所占的份量在不断的增加。,生物学家的评论,生物大分子结构测定影响了生物学几乎所有的领域。 几乎每个星期都有激动人心的结构发表在如Science和Nature这样的杂志上。 生物大分子晶体学已经比其他任何学科更多地得益于同步辐射的应用。 -Steven E. Ealick, Cornell Univ.,生物大分子晶体,衍射能力很弱：大部分原子

3、为C、N、O、S； 晶体质量差； 相位问题解决困难。,同步辐射的高亮度、准直性好、能量 可调正好提供了适合的光源,为什么要使用同步辐射,高亮度：衍射弱的晶体可以得到好的信号； 准直性好：质量差的晶体可以进行实验； 能量可调：可以进行多波长反常衍射实验（解决全新结构的有效方法）。,衍射数据,MR: CNS,CCP4,SIR/MIR/SAD/MAD： Solve/Resolve, shelxd, SnB, CNS, Hyss, OASIS,CCP4,sharp,数据处理： HKL2000， Automar, mosflm, XDS,建模： ARP/wARP，Resolve O, Xtalview,

4、模型修正： CNS,模型检验： CNS CCP4,解析结构的流程,收集衍射数据的目的,确定出准确的结构因子F(hkl)。,收集单晶衍射数据的方法,回摆法：最常用 旋进法 魏森堡法,回摆法：rotation (oscillation) method,入射X光,单晶,绕某个轴旋转,探测器（CCD、IP）,晶体的安装方式,毛细管：无法冷却,loop：冷却,蛋白质晶体的衍射,单晶样品冷冻,使用液氮冷却的干燥空气（防止结冰）来冷冻样品。 实践证明可以大大延长蛋白质晶体耐受同步辐射辐照的时间。 对于花费时间长的实验（如MAD）是必不可少的。 一般用100K的温度就可满足要求。,收集衍射数据的策略,是否需要

5、冷冻？ 为了减少辐射衰减，需要冷冻。但是冷冻可能导致晶体损坏，往往需要防冻剂（有些晶体生长时使用的沉淀剂本身就有防冻效果）。 防冻剂需要尝试。（花费大量时间和精力的工作！） 不要破坏晶体，不出现冰的衍射（不结冰）。,选取冷冻条件,可能晶体本身就有抗冻能力。 安装晶体时尽量减少（生长晶体所用的）母液。 最简单的方法：把母液中部分水替换成甘油或蔗糖。 先找出没有结冰的替换条件，再尝试晶体是否不被破坏。 使用Hampton提供的防冻剂。 使用矿物油。,是否为需要的蛋白晶体,生长出来的晶体可能是其他成分，特别是结晶所使用的小分子试剂的晶体（盐晶）。 在衍射上可以判断：蛋白质晶体晶胞很大，起码有几十埃，

6、衍射点很多，盐晶最大只有十几埃，衍射点很少。 蛋白晶体偏光弱，感觉很软（像果冻）；盐晶偏光强，感觉很硬（像盐粒）。,晶体的单晶性,如果要顺利的解析出结构，晶体的单晶性要好； 仔细观察衍射点的形状和衍射的上限：衍射点要圆而锐利，能达到的衍射上限不要低于3埃，否则结构解析就难了。 出现两个点很靠近的情况要小心：这可能不是单晶（孪晶），也会是晶胞某个轴特别长（移动探测器距离，处理一下衍射图）。一定要排除孪晶的可能性。,适合衍射的晶体,得到适合衍射的晶体是一件相当花费时间和精力的任务。 首先要得到大量的纯净蛋白质（分子生物学表达重组蛋白，生化手段纯化）； 寻找可能的结晶条件：如用Hampton公司的筛

7、选试剂盒Index，Screen I和Screen II。,优化结晶条件：可能优化的条件很多，如沉淀剂浓度，盐浓度，离子强度，pH值，蛋白浓度等等等等。 得到大量的、衍射能力好的晶体。 筛选合适的防冻剂。 每一步都可能遇上困难。,探测器距离的设置,距离近可以收集到的衍射分辨率高，但是衍射点可能重叠（overlap）； 距离远可以使信噪比提高，衍射点分辨得好，但是能收集的分辨率降低。 经验：晶体中最长的晶胞长度数值（以埃为单位）=晶体到探测器距离数值（以毫米为单位）。 或者探测器最外面正好达到衍射分辨率上限。 高、低分辨率数据可能要分别收集（如果晶体分辨率极高）。,回摆角度、画面张数、曝光时间,

8、回摆角度一般取1度。但是晶胞很大的情况下需要减少角度。不要出现回摆角度过大导致衍射点重叠。 需要收集的画面张数取决于晶体的对称性，至少要覆盖一个独立衍射区。当然在时间许可的情况下，越多越好。 曝光时间取决于晶体衍射能力，尽量使探测器的动态范围得到利用，也要考虑晶体的辐射耐受能力。,波长选择,波长越长，衍射强度会大（l3）；但是吸收效应不好处理。 同步辐射上还是需要较短的波长，除非反常衍射需要调节特定的波长。 一般的束线能量优化在1埃左右，选择最优强度的波长最好。 短的波长使得衍射点间距减少，可以考虑探测器距离加大一些。,晶体尺寸和外形,同步辐射由于强度高，不需要大的晶体。其实太大的晶体不容易冷

9、冻。一般100-200微米即可。 好看的晶体衍射不一定好：衍射说了算。,安装晶体的方位,一般不需要考虑：由于处理技术的发展，随机取向的晶体就可以了。 选择特殊取向可以减少衍射盲区（即收集不到的衍射点区域），但是也会导致盲区增加。有时候反常衍射需要选择特殊的取向。 除非有精确调节晶体取向的衍射仪（三圆衍射仪），一圆衍射仪是随机取向装样的。经验表明，一旦晶体取向的某个轴和入射光太接近，会导致数据处理困难。,相位问题的解决,解决相位问题的主要方法,同晶差值付里叶法 分子置换法 同晶置换法 反常衍射法,反常衍射,如果得到的晶体中含有重金属（蛋白本身含有、硒代引入、浸泡重元素引入的都可以），采集数据时就

10、可以考虑作反常衍射。 现在利用单波长反常衍射（SAD）解析结构的数目已经接近或超过MAD。所以多数情况下一个波长的反常衍射已经足够了。,MAD,选择吸收边附近的几个波长，相当于几个不同的完全同晶的置换。 分子生物学提供了硒代的手段，对于解析全新结构非常有利。,波长选择,f1,f2,f3,f4,四个波长,f1：peak，f最大； f2：inflection，f最大； f3：high energy remote； f4：low energy remote。 流行做法：一边收集数据，一边解析结构。一般完成f1数据的收集，进行f2数据收集时就完成相位解析了。,数据,同晶置换：I_native（nati

11、ve晶体的衍射强度），I_h（同晶置换晶体的衍射强度）； 反常衍射：I+，I-（Friedel对的强度）。,寻找重原子位置,同晶置换法得到的晶体可以看成在原有的晶胞内的原子（位置在xi，yi，zi，一共N个）外，加入了某些重原子（位置在xhj，yhj，zhj，一共M个）。,AB*，A*B由于重原子和其他原子之间位置没有相干关系，这一项求和的结果接近于零。所以IH-Inative主要是重原子的贡献。 通过差值patterson图可以得到重原子位置，然后根据F的关系求出native的相位来。,反常衍射也是类似处理：I+和I-主要是反常散射原子的贡献。 得到反常散射中心的位置后，在求解相位。 反常衍

12、射可以看作完全同晶的置换。,总结：,同晶置换和反常衍射解相位的一般步骤： 差值Patterson图寻找重原子； 精修重原子位置； 解出相位。 Solve：自动化较好的程序，以上步骤自动完成； CNS：需要一步一步来进行。,S vs M,单对同晶置换或者单波长反常衍射都存在相位双解问题。以前是利用多对同晶置换或者多波长反常衍射解决的。 现在通过与直接法的结合，不需要这么做了。 直接法可以计算出两个解的概率，选取概率大的那个作为真实解。,Solve：双解中利用sim概率选一个可能性较大的解，得到电子密度图后，利用溶剂平滑改善。 OASIS：再加上cochran概率，选取的解更接近真实解，其结果比s

13、olve好。,例子：Se-MAD,大肠杆菌（escherichia coli，E. coli）的巯基化酶II（thioesterase ii）。 Se-MAD收集了peak，inflection，low energy remote，high energy remote的数据，分辨率为2.5A。 利用solve解决相位问题。,solve/resolve软件包,由LANL的Tom Terwilliger开发的解析SAD/MAD/SIR/MIR的软件包。 运行方式：编辑一个script（脚本），里面包含了需要读入的数据文件，格式，需要进行的运算等。,script,#!/bin/csh -f # se

14、tenv CCP4_OPEN UNKNOWN setenv SOLVEDIR /usr/local/lib/solve #solve solve.log title 4-wavelength MAD dataset ! a title for this dataset solve.setup ! get our standard information read in logfile mad.logfile ! write out most information to this file. ! summary info will be written to solve.prt mad_ato

15、m se ! the anomalously scattering atom is selenium refscattfactors ! do not refine scattering factors (you can if ! you want as long as the data is good) readdenzo ! data is from scalepack. Note: for best results ! use no merge original index in scalepack ! Alternative to readdenzo is readformatted

16、premerged ! data have not been merged to the asymm unit ! Alternative to unmerged is premerged,script,lambda 1 ! info on wavelength #1 follows label Peak Wavelength ! a label for this wavelength rawmadfile ./jia_peak.sca ! datafile for lambda 1 ; raw Intensities wavelength 0.9787 ! wavelength value

17、fprimv_mad -5.9 ! f value at this wavelength fprprv_mad 4.1 ! f value at this wavelength 四个波长的数据都输入。 f和f的值不用很准确，solve会拟合这两个值。,script,nres 570 !approx # of residues of protein in asymmetric unit nanomalous 8 !approx # of anomalously scattering atoms in a.u. SCALE_MAD ! read in and localscale the data

18、 ANALYZE_MAD ! run MADMRG and MADBST and analyze all the Pattersons SOLVE ! Solve the structure Exit 这里告诉solve程序的内容包括晶体的空间群（在solve.setup中），数据格式（denzo处理的数据，已经scale过了），数据文件名，蛋白中的氨基酸数目，反常散射原子数目等。 然后运行solve。,结果：重原子位置,-TOP SOLUTION FOUND BY SOLVE ( = 0.56; score = 48.66) - X Y Z OCCUP B HEIGHT/SIGMA 1 0.239 0.094 0.150 0.536 19.1 27.6 1 0.807 0.423 0.156 0.612 25.4 28.3 1 0.244 0.080 0.098 0.689 32.4 28.1 1 0.800 0.433 0.102 0.646 30.2 27.9 1 0.156 0.410 0.243 0.669 36.3 24.4 1 0.198 0.118