细胞培养、克隆计数与平均存活分数(SF)的计算将接受照射后的细胞
放回细胞培养箱继续培养,每三天更换新鲜培养液一次。
待培养至10-14天后
弃去培养基,用0.5%结晶紫染色贴壁细胞。
计数并记录不同照射剂量下各组培养瓶中细胞存活克隆个数。
根据假照射处理组细胞克隆形成数目,计算细胞接种效率(PE;各照射处理组细胞存活分数(SF)为各组细胞克隆形成的平均数除以相应细胞接种数与接种效率的乘积。
细胞存活分数(SF)与辐射剂量之间数学模型的建立
哺乳细胞接受不同剂量的电离辐射照射后,细胞存活的规律大多是符合单击多靶细胞存活模型,即SF=1-[1-EXP(-k*D)**n]。
根据实验实测细胞存活数
据进行该数学模型的非线性回归拟合分析,可获得单击多靶细胞存活模型的一系列有意义的参数,如k、n、Do、D37和De。
其中,k为细胞存活曲线的钝化常数,其值可由拟合曲线方程直接给出;n称外推数,代表理论上的靶数目,数值大小也是由拟合方程直接给出。
Do、D37和Dq是根据k、n值演算出来的,Db称为平均致死剂量,是理论上使每个细胞平均被击中一次所需要的辐射剂量,数值上Do
=1/k;Dq称为准阈剂量,表征修复亚致死损伤的能力,即克服单击多靶曲线中
“肩宽”的剂量,数值上Dq=InnxDo;6是细胞存活分数为37%寸相应的辐射剂量,在单击多靶模型中D37=D0+Dq。
细胞存活单击多靶模型(下图)的表
达式为:
SPSS软件对单击多靶模型的数学拟合,求出模型中k和n参数值大致步骤如下:
计算、整理不同辐射吸收剂量D下细胞的平均存活分数SF(校正接种效率后的平均存活分数,这是第一步),制成剂量D与存活分数SF的二维表
格;启动SPSS软件,定义变量名D和SF、字段宽度以及小数点后位数,将SF显示精度定为小数点后4位数,并输入D和SF的原始数据;选择SPSS软件中非线性回归(Nonlinearregression)窗口,将SF设置为因变量(Dependent),数学公式表达模型中输入1-[1-EXP[-k*D]**n],同时将模型参数k和n迭代的起始值分别设置为0.0001,同时将模型参数k和n数值各自限制在大于或等于0.0001以内,保存上述分析条件;迭代拟合的算法方式选择levenburg-
Marquadt,为绘制存活曲线图方便,显示项中选择根据拟合方程计算出的细胞存活分数的理论预测值;执行运算并显示出SPSSt有的结果文件。
在运算结果显示文件中有不同迭代过程剩余的残差平方和、根据1-[1-EXP[-k*D]**n]拟合后
的最终优选k和n参数值、以及k和n的统计学显著性检验概率等。
最终优选的k和n参数值最为重要,是我们模拟计算的目的。
此节是本次实验的重点与难点,有关详细过程可参考六、SPSS软件使用指南参考。
单击多靶模型中DO、Dq和D37值,以及增敏比SER勺计算
根据最终优选的k和n参数值,我们可以进一步计算出DO,D37和Dq,其相互之间的关系式如下:
sf=丄一a—&屁口)直.
5-
A:
'-1-
——1±X11XZT口
=S+Dm
同一细胞株接受不同因素(如乏氧、增敏剂等)处理后再行照射,细胞存
活曲线及其参数会发生变化,根据这些参数的变化我们还可进一步计算出增敏比
(sensitivityenhaneeratio,SER,SER增敏比反映细胞放射敏感性的不同
及其变化。
SEF一般大于或等于1,SEF越大说明经该试剂处理或缺氧后对辐射表现出更加敏感。
SER=Dq(未接受化学处理或有氧情况下对照组)/Dq(接受化学处理或无氧情况下处理组)
有研究表明,实体瘤细胞因缺血缺氧而表现出肿瘤细胞的内质网应激。
Tun
(衣霉素)具有抑制细胞内糖蛋白的糖基化修饰,从而导致内质网应激的发生。
本次实习是利用Tun和另一种内质网应激试剂DTT—起对肿瘤细胞株HeLa细胞进行预处理,再行电离辐射照射。
因此本次实验我们要求同学利用附后的实验结果,对内质网应激是否具有增强辐射敏感性的作用提出评价,另外如若要做出此
实验结论时您还希望要做的进一步补充实验是什么,请说明为什么。
利用EXCEL软件对细胞存活曲线的绘制最终优选的k和n参数值后除
了可以进一步计算单击多靶模型中DO、Dq和D37值以及增敏比SER外,我们还要求同学利用k和n参数值回代到该模拟公式中,对不同辐射剂量D下的存活分数
进行理论预测,然后利用理论预测的存活分数SF值与实际给予的辐射剂量D之
间绘制一条平滑的存活曲线。
绘制过程大致步骤如下:
启动EXCELS件;输入不同剂量值D;fx栏输入=1-P0WER[(1-EXP(-k*D,n]公式,其中k和n值换成来自SPSS模拟方程的计算结果,运算各剂量组细胞存活分数预测值;根据各组剂量值D与细胞存活分数
SF预测值建立数据二维表;启动绘图XY散点标准模式,输入标题、纵横轴名称
【如HeLa细胞伽玛照射后存活曲线单击多靶模型图、SF、Dose(Gy)】,去掉网
格线、图背景阴影、图例和线框;点击丫轴坐标轴,选刻度项,修改最小刻度0.01,最大刻度10,主要刻度10,次要刻度10,数轴交叉于0.001,选对数刻度;点击各组各点,选数据系列格式,线型选自定义,线形平滑;图幅大小设置5X20。
获得曲线图之后,我们也可以在不同预测曲线上添加出实测存活分数值,同学们由此可观察到理论预测值与实际值之间的真实差异,这种差异也可以反映在运算结果显示文件中显著性检验分析上。
四、实习作业
演算练习,根据以下Tun和DTT组数据计算各自D增敏比(SER并作图
cell
seeding
CTR
survival
cell
CTRmeancell
Tun
survival
cell
Tunmeancell
DTT
survival
cell
DTT
survival
cell
Dose
number
number
number
number
number
number
mean
0
600
411
390
300
0
600
427
427
365
0
600
442
427
442
420
432
366
0.5
600
283
362
293
0.5
600
324
349
316
0.5
600
217
275
280
330
290
300
1
1000
322
344
353
1
1000
202
380
301
1
1000
288
271
327
350
324
326
2
4000
344
492
511
2
4000
357
408
458
2
4000
342
348
601
500
546
505
6
8000
260
305
300
6
8000
233
324
298
6
8000
203
232
281
303
363
320
8
10000
88
122
169
五、SPSS软件使用指南参考:
下面以原始数据进行1-[1-EXP[-k*D]**n]非线性回归方程拟合为例详细讲
解SPSS软件的使用。
照射剂量与平均细胞存活分数原始数据整理如下:
组别
Dose(Gy)
平均(SF)
0
0.00
:
计算PE用
1
0.00
1.00
2
2.00
[0.70
3
4.00
0.30
4
6.00
[0.10
5
8.00
0.02
第一步:
启动SPSS定义变量名,并把SF显示为精度4位数,便于和拟合数据比较!
!
其中定义变量名非常关键
两Untitled-SPSSEditor
3凶
File£jdit
Vtew巫恥
JjransformAnalyE*GraphsUtilitiefWindowyelp
胡翔國闸州旳垃|@|gb|g|
p:
'I
VARD00011VAR00002]
var
I
A.
—
1
00
A1-00
2
2.00
7g
3
400
丄.30
4
6.00
JO
5
8.00
r~
-I
*I>l\DataViewk
Vaiiable
HLd
;PS5Processor
isready
I
e旧1圖團碉iWIm|SBHHi
引QI
9:
SP
n
CSF>
var
rar
V3r
var
var斗
1
.00
I1.0000
―1
2
200
.7000
3
4.00
.3000
4
6阿
~~J000
5
8.00(
0200
G
.—r
itaView/VariableViewf
[
』
SPSSProcessorisready
I
/
FfeE(ftView国沽Traf^oniflAnalyze&aphsUtilities塑rxfow曲p
园UrMed-SPSSDataEditor
fQlxl
第二步:
激活Analysis菜单选Regression中Nonlinear...项,弹出NonlinearRegression对话框,依次进行如下图中操作,其中ModelExpression的选项非常关键,本例为1-(1-EXP(-k*D))**N
第三步:
点击上图中Parameters,出现下面的对话框,进行如下图的设置
NonlineiarRegrerssioin:
ParameterConstraiints
第四步:
打开限定条件对话框后,依次进行如下图的设置,至于为什么要进行这种限定,以及这样限定的合理性,简单解释:
k肯定大于0,N也肯定大于0,所以就
设置他们都大于一个很小的数字吧,这样设置的巧妙性在于避免SPSS拟合过程中出现诸如(-3)9.5这样的死步骤。
给好限定条件后,会跳出一个框子提醒你什么什么什么,不管他,继续望下走
Functionprecision
第五步:
回到拟合界面对话框后,依次进行下图操作,目的是根据拟合方程显示理论预测值,其中拟合算法很关键,请一定要按照图中所示进行,最后点击0K—切结果就显示出来啦!
!
ConstrainedNonlinearRegression
Allthederivativeswillbecalculatednumerically.Thefollowingnewvariablesarebeingcreated:
NameLabel
PRED_PredictedValues
IterationResidualSSkN
0.1
.5885643726.000100000.000100000
1.1
.4326928450.000100000.010267659
2.1
.2644240765.000100000.032263147
3.2
.2460533309.000100000.040666221
4.1
.2353237999.000550799.058576000
5.1
.2343114615.000749210.065273484
6.1
.2310963020.001052919.068291791
7.1
.2267676825.001668720.074229208
8.1
.2253204323.001932692.071071008
9.1
.2252782354.010525019.155223311
10.1
.2114927788.007972194.123511525
11.1
.2072941880.007258942.105891983
12.1
.2037004224.008892185.114062200
13.1
.1849657836.020745415.163259238
14.1
.1678978737.036031833.217246455
15.1
.1359764487.073514715.341227935
16.1
.0696285850.177899360.683299005
17.1
.0193661970.3970542741.40880796
18.1
.0192417786.3931284801.39728887
19.1
.0187180868.3790873691.36145913
20.1
.0180956915.3717613231.35257103
21.1
.0160190997.3640753701.38506874
22.1
.0126064378.3751342921.51881114
23.1
.0059547114.4289360571.90344372
24.1
.0019056787.5007576152.39973577
25.1
.0011632693.5290739082.63477892
26.1
.0006871616.5469015662.82717444
27.1
.0003546795.5574224162.99265985
28.1
.0002842715.5604100253.06737893
29.1
.0002815003.5614781273.08761846
30.1
.0002814377.5618414173.09187019
31.1
.0002814332.5619615613.09297517
32.1
.0002814332.5619699613.09301099
Runstoppedafter32majoriterations.
Optimalsolutionfound.
NonlinearRegressionSummaryStatisticsDependentVariableSF
Source
DFSumofSquaresMeanSquare
Regression21.59012.79506
Residual32.814332E-049.381106E-05
UncorrectedTotal51.59040
(CorrectedTotal)4.69152
Rsquared=1-ResidualSS/CorrectedSS=.99959
Asymptotic95%
AsymptoticConfidenceInterval
ParameterEstimateStd.ErrorLowerUpper
k.561969961.019670604.499369321.624570602
N3.093010989.2050164822.4405570433.745464936
AsymptoticCorrelationMatrixoftheParameterEstimates
kNk1.0000.9452
N.94521.0000
第六步:
阅读结果文件。
结果文件详细显示了叠代步骤和最终优选的参数值,以及残差平方和等等,本例结果如下:
结果分析如下:
SF=1-(1-EXP(-0.562*D))**3.093
根据上式k=0.562和N=3.093可以演算出以下参数:
D0=1/k=1/0.562=1.779
Dq=lnn*D0=ln3.093*1.779=2.009
D37=D0+Dq=1.779+2.009=3.788