全 文 :生命科学
Chinese Bulletin of Life Sciences
第22卷 第8期
2010年8月
Vol. 22, No. 8
Aug., 2010
文章编号 :1004-0374(2010)08-0817-06
原子力显微镜研究细胞弹性的数据分析方法
叶志义*, 张 丽
(重庆大学生物工程学院生物流变科学与技术教育部重点实验室, 重庆 400044)
摘 要:细胞的力学性质与生命体的功能和健康都是息息相关的,而原子力显微镜(atomic fo r c e
microscope,AFM)是研究细胞力学性质最好的仪器之一。该文较为详细叙述了 AFM 力曲线测量原理,
以及分析 AF M 实验数据的常用模型。作者用该文所提到的模型分析了肺癌细胞的弹性,结果显示在加
载速率比较低的情况下,三种计算方法计算得到的结果没有显著差别,而在加载速率高于8 μm/s 时,结
果则有显著差异。
关键词:原子力显微镜;力曲线;杨氏模量;模型
中图分类号:Q-336;Q2-33 文献标识码:A
Elasticity measurement of cells with AFM: the data processing
YE Zhi-yi*, ZHANG Li
(Key Laboratory of Biorheological Science and Technology, Ministry of Education, Bioengineering College,
Chongqing University, Chongqing 400044, China)
Abstract: Mechanical properties of individual living cells which can be probed by AFM, one of the best tools,
are known to be closely related to both health and function of the human body. Measurement principle of force
curve of AFM, as well as several models applied to analyze the experimental data, is introduced in this review.
Furthermore, elasticity of lung cancer cells is calculated according to the models proposed in this paper, which
indicates that corresponding Youngs moduluses are approximate under slow loading rate, but different signifi-
cantly when loading rate is greater than 8 μm/s.
Key words: AFM; force curve; Youngs modulus; model
细胞的物理性质,如细胞的硬度、弹性以及
黏性等,在生物体生命活动中起着非常重要的作
用。这些物理性质的改变既可能改变细胞的生理功
能,使得生命体产生疾病,同时,细胞功能的异
常也可以通过这些物理性质表征出来。因此,研究
细胞的物理性质对探索生命体的生理功能有着非常
直接的意义。原子力显微镜具有高分辨率和能在生
理状态下研究生物样品的优点,已经广泛应用于细
胞黏弹性、细胞硬度以及细胞表面分子间相互作用
的研究。Yang等[1]研究了细胞表面黏附分子整合素
Mac-1 和 ICAM-1 相互作用。AFM 单分子力谱结果
表明,Mac-1 在静息状态下与ICAM-1 有一定结合,
被活化后对ICAM-1 结合力增大。Bálint 等[2]用 AFM
研究甘露醇对脑血管内皮细胞的形状和力学性质的
影响。通过 AFM 力曲线测量细胞的弹性,发现未
经甘露醇处理细胞的杨氏弹性模量是(8.04±0.12) kPa,
而经过甘露醇处理后细胞变软了,其杨氏弹性模量
是(0.93±0.04) kPa。Kasas等[3]用AFM研究细胞的膜
下蛋白和深层微管的力学性质,这一研究表明,
AFM 成为研究细胞表层和深层的生物力学性质的一
个重要工具。Cross 等[4]用 AFM测量来自肺癌、乳
收稿日期:2009-12-30;修回日期:2010-01-28
基金项目:国家自然科学基金项目(30570453); 高等学
校学科创新引智计划(“111”计划)(B06023)
*通讯作者:E-mail: yzybioafm@163.com; Tel:023-
65112452
· 技术与应用 ·
818 生命科学 第22卷
腺癌和胰腺癌患者胸膜转移癌细胞的硬度,发现转
移癌细胞的细胞硬度比良性细胞软70%。被氧化的
低密度脂蛋白(ox-LDL) 对内皮细胞的伤害起着很重
要的作用,Chouinard 等[5]用 AFM 对 HUVEC 细胞
(human umbilical vein endothelial cells,人脐静脉血
管内皮细胞)进行力曲线的测量,结果表明ox-LDL
对 HUVEC 的硬度有影响,而 LDL 对 HUVEC 的硬
度基本没影响,该研究在一定程度上说明细胞物理
性质能够反映出细胞的功能。利用AFM 分析细胞力
学性质的研究已经是非常热门的一个领域,要想更
好地利用AFM 研究细胞的物理性质,必须要深刻认
识AFM力曲线的测量原理和数据分析方法。目前分
析力曲线的最经典模型是赫兹模型。对于细胞而
言,虽然并不完全符合以上的假设,但是,赫兹
模型对于研究细胞的弹性仍然是非常有用的。赫兹
模型把细胞考虑成一个弹性体,忽略了细胞的黏
性。此外Sirghi等[6]利用力曲线回针曲线分析细胞的
物理性质(该理论以下简称为回针曲线模型),这一
理论考虑了细胞与探针之间的黏附力,而且还能够
获取细胞与探针之间的黏附力。本文将介绍AFM力
曲线测量原理和三种分析力曲线的方法,并且用这
三种方法进行分析,并以此进一步分析这三种方法
之间的异同。
1 AFM 力曲线测量原理
A F M 探针在靠近、接触和压入细胞这一过程
中,探针微悬臂由于受到细胞的阻碍而弯曲变形,
使得通过微悬臂反射得到光敏二极管阵列的激光信
号发生变化而记录下探针所受到的力变化过程。
AFM 力曲线是将 AFM 探针与细胞表面之间的距离
(distance)作为横坐标,将探针微悬臂的变形(cantilever
deflection)作为纵坐标而得到的曲线,根据胡克定律,
探针所受到的作用力(force)与探针微悬臂变形量之
间的关系:力= 微悬臂形变量 × 探针弹性系数,其
中探针弹性系数可以由热振动得到。力曲线测量过
程示意图如图 1 所示。
从图 1 可以看到,在探针没有接触到细胞之
前,探针微悬臂没有变化,此时,力曲线是一条
几乎平行于横坐标的直线;当探针接触细胞,并继
续压细胞时,探针由于受到细胞对它的反作用力而
使得微悬臂变形,力曲线开始慢慢上走。细胞膜性
质不一样,探针受到的反作用力就会不一样,所得
到的力曲线也会有差别。
2 AFM 数据分析方法
2.1 赫兹模型
赫兹模型[7]是赫兹在1882年建立的模型,主要
用于计算弹性材料的力学性质。赫兹模型作了一些
假设,比如:表面是连续的、光滑的,变形是相
对很小的。对于细胞而言,虽然并不完全符合以上
的假设,但是,赫兹模型对于研究细胞的弹性仍然
是非常有用的,在细胞的变形小于细胞高度的一半
时[8],认为细胞的变形规律符合弹性体。因此,赫
兹模型广泛应用于细胞力学性质的计算,是现在应
用最为广泛的一个模型。
根据 AF M 探针形状的不同,赫兹模型有两个
图1 AFM 力曲线测量过程示意图
819第8期 原子力显微镜研究细胞弹性的数据分析方法
表达式,当探针形状为锥形时,赫兹模型方程为:
2
2
2 tan
1cone
Ef δπ ν= ∂− ............................... (1)
当探针形状是球形时,赫兹模型表达式为:
3
2
2
4
3 1sphere
Ef Rδν= − .................. (2)
其中,f 是探针对细胞所施加的压力,E 是杨
氏模量,v 是泊松比,生物样本取值为0.5[9,10],∂
是探针半开角,R 是探针半径,δ 是压痕深度。下
面我们列出两种常见的计算方法。
2.1.1 斜率法[11]
我们将式子(1)(2)等号两边同时乘方,线性化
后分别如下:
1 1
2 2
2
2( ) ( tan )
1cone
Ef δπ ν= ∂− ............. (3)
2 2
3 3
2
4( ) ( )
3 1sphere
Ef R δν= − .............. (4)
从(3)(4)两式可以看到,通过变形之后,fcone
1
2
——δ (或者(fsphere 23 ——δ )呈一次关系式,通过对
力曲线进行拟合之后可以得到杨氏模量E,分别如
以下两式所述:
1
1 2
2
2
( )2( tan )
1
conefE slopπ υ δ
Δ∂ = =− Δ ....... (5)
2
2 3
3
2
( )4( )
3 1
spherefE R slopν δ
Δ= =− Δ ........ (6)
其中slop 是斜率,所以,杨氏模量的计算方
法如下:
1
2 22
2 2( ) 1 1( ) ( )
2 tan 2 tan
conefE slopπ ν π νδ
Δ − −= =Δ ∂ ∂ ... (7)
3 32 2
2 2
3
2
( )3 1 3 1( ) ( )
4 4
spherefE slop
R R
ν ν
δ
Δ − −= =
Δ (8)
Carl和Schillers [11]应用这种方法分析了中国仓
鼠卵巢细胞在生理条件下的力学性质,得到了细胞
不同层面所表现的弹性不同,如图2所示。他们的
结果显示:在膜下500 nm 时,细胞表现出来的弹
性约340 Pa,在膜下约1 000 nm的压痕下表现出来
的弹性是525 Pa,造成这种差异的原因可能是由于
细胞内部骨架分布不同,前者只要是表现出膜的性
质,如膜蛋白、膜下骨架的性质;后者是由细胞
内的骨架结构引起的,如微管、微丝等网状结构。
线性化的第二种处理方法[12]是直接将式子(1)(2)
图2 赫兹模型斜率法示意图[11]
将力曲线的进针曲线变形得到
2
3f δ− 的曲线关系,再将该曲线线性拟合得到两段直线,其斜率即是该压痕下细胞膜的杨氏
模量 E 。插入图是一条完整的力曲线。
820 生命科学 第22卷
分别看成(fcone)—— δ 2、(fsphere —— δ
2
3 的一次关系
式,得到两者的关系图后,拟合得到曲线斜率
slop,这样杨氏模量就可以计算得到:
2 2
2 2
2
( ) 1 1( ) ( )
2 tan 2 tan
conefE slopπ ν π νδ
Δ − −= =Δ ∂ ∂ . (9)
3 32 2
2 2
3
2
( )3 1 3 1( ) ( )
4 4
spherefE slop
R R
ν ν
δ
Δ − −= =
Δ (10)
Tamara等[12]用这种处理方法得到了人内皮细胞
的弹性,如图3 所示,同样内皮细胞在不同层面所
表现的弹性也有所差异。
线性化方法分析数据的好处在于不用找接触
点,直接将原始数据乘方后作图拟合得到斜率即
可,而且还能同时得到不同深度压痕下细胞的弹性
差异。但是这种方法要求力曲线比较平滑,而且测
量环境中的液体黏性不能过大,所以在AFM实验过
程中要求仪器稳定性好,外界噪音很小,否则得到
的力曲线所拟合成的直线线性关系不好,导致难以
确定直线的斜率而无法很好地得到杨氏模量。
2.1.2 两点法[13]
由(1)(2)两个式子可以看出,只要知道f和 δ,
就可以求出样本杨氏模量E。如果我们找到探针与
细胞的接触点,在力曲线中得到接触点所对应的压
力fcontact和距离dcontact,以及我们所测量压痕深度对
应的结束点所对应的压力fend和距离dend,就可以得
到在这一深度压痕处探针施加给细胞的压力f和细胞
压痕 δ,再利用赫兹模型求出细胞的杨氏模量 E。
下面给出这种方法的计算过程:
| |end contactf f f= − ..................... (11)
| |end contactd dδ = − ..................... (12)
将上面两个式子代入式子(1)(2)分别得到:
2
2
2| | tan ( )
1end contact end contact
Ef f f d dπ ν− = ∂ −− (13)
3
2
2
4| | ( )
3 1end contact end contact
Ef f R d dν− = −− . (14)
其中,fend、fcontact分别是力曲线中结束点和接
触点所对应的力,dend、dcontact 分别是力曲线中结束
点和接触点所对应的距离值,如图4所示。所以我
们确定接触点和压痕深度之后,代入相应数据就可
以求解出样本杨式模量E。Charles等[13]用这种方法
分析了活细胞膜的物理性质。这种方法计算简单,
但是因为AFM实验过程中受到很多因素影响[14],比
如探针移动过程中液体对探针的作用力、细胞本身
的黏性,以及实验过程中环境的噪音和系统本身的
不稳定性等,所以很难准确地确定接触点。
2.2 回针曲线模型[15]
赫兹模型是通过分析力曲线中进针曲线来得到
图3 赫兹模型斜率法示意图[12]
将力曲线的进针曲线变形得到
3
2f δ− 的曲线关系,再将该曲
线线性拟合得到两段直线,其斜率即是该压痕下细胞膜的杨
氏模量 E
图4 两点法示意图
如图是一条力曲线的进针曲线,在图上已经标明接触点
(contact point)和某一深度压痕处的结束点(end point),(我们
以压痕500 nm 为例),将这两个点分别对应的距离dcontact、
dend和压力值fcontact、fend代入式子(13)或者式子(14)就可以求
得杨氏模量 E 。
821第8期 原子力显微镜研究细胞弹性的数据分析方法
杨氏模量,Oliver和Pharr[15]发现从回针曲线中也能
得到杨氏模量,而且还能得到所研究样本的黏性。
Sirghi等[6]利用AFM压痕曲线的回针曲线分析得到了
细胞膜的弹性和黏性,其表达式如下:
*
2
2
324 tan tan( )
cos
aEP h h hγα ααππ π= − ………………… (15)
其中P(h)表示微悬臂变形为h时探针施加的压
力,∂ 是探针半开角,h 是微悬臂变形量,E* 是杨
氏模量 E 的近似值,两者的关系是
E*=E/(1-v 2) ......................... (16)
ra 是细胞黏性系数,它与探针抵抗细胞黏性做
功 W a 的关系是:
Wa=-γaA c ................................................ (17)
Ac 是探针与细胞的接触面积
该模型示意图如图 5 所示:
由于AFM所获得的力曲线通常会受到黏附力的
影响,探针与细胞之间的黏附力会同时影响压痕实
验中加载和非加载部分,在加载区域它会加速探针
的位移,在非加载区域它会阻止探针的回拉。因
此,若是忽略探针与细胞之间的黏附力,将会给最
后结果带来系统性的误差[6]。回针曲线模型是Sirghi
等[6]建立的理论分析方法,考虑AFM 所得到的力曲
线受到探针与细胞之间黏附力的影响,而且利用该
模型可以获取细胞表面的弹性信息以及探针与细胞
膜表面之间黏附力信息。黏附力可以用探针和细胞
接触面积上的黏附能量来表示,而且黏附能量和接
触面积在整个压痕过程中是处于一个连续的变化状
态。探针和细胞之间的黏附力同时影响压痕过程中
的进针和回针,但是在非加载(回针)时间段可以直
接获得黏附力的大小。当黏附力是探针和活体细胞
接触区域上大量分子相互作用的结果时,需要引用
一个更适合的参数 γa,它指的是单位面积上分开两
个物体所需要的能量值。
这个理论分析方法能够较为准确地计算出探针
与细胞发生力学接触的面积,在压痕的非加载过程
中细胞骨架弹性会发生松弛,并且细胞和探针接触
面积会连续地变小,探针和细胞之间接触面积的连
续变化决定了两者之间黏附能量和黏附力的变化。
从锥形探针对活体细胞进行的压痕实验来看,探针
和细胞膜之间存在有限的黏附力。通过该理论模型
拟合力曲线可以获取细胞表面弹性以及细胞和探针
接触面上黏附能量的信息。
3 讨论
本文较为详细地介绍了分析AFM力曲线求解细
胞弹性的常用模型,下面分别用以上三种方法分析
实验数据三种方法的异同。用AFM在不同加载速率
下对肺癌细胞 A549 做力曲线,分别用以上三个方
法得到杨氏模量,其结果如图 6 所示。
从结果可以看到,随着加载速率的增加,杨
氏模量也随之增大,这与很多研究是一致[14],这可能
是因为受到细胞黏性的影响。在加载速率较小时(小
于8 μm/s),三种方法计算得到的杨氏模量没有显著性
差异;在加载速率较大的时候(大于8 μm/s),从图中
可以看到,回针曲线模型得到的杨氏模量更小,这
可能是因为回针曲线模型考虑了细胞的黏性影响。
有研究表明,在加载速率小于10 μm/s 的情况下,
液体对探针的影响很小[11],加载速率越大,细胞黏
性对AFM 杨氏模量的测量结果影响越大[16]。所以,
在加载速率较大的情况下,用回针曲线模型可能会
更为符合实际,在加载速率较小的情况下,赫兹模
图5 探针与细胞接触示意图[6]
a 图是坚硬的锥形探针对一个半空间软的弹性材料的压痕示
意图(其中 α 为锥形探针的半开角,he 表征的是材料表面的
变形,hc 为探针接触深度,h 为探针向细胞内移动的整个距
离深度); b图为与细胞完全接触的那部分探针的立体示意图
(其中 A 为接触面的投影面积,lc 接触面的边长)
822 生命科学 第22卷
型计算更简单一些。
文献调查发现,不同研究组所测到的细胞杨氏
模量差异非常显著[17-20],一方面这是由不同细胞本
身性质导致的,另一方面,使用不同模型进行数据
分析,使用不同仪器进行实验,这些对测量结果的
影响可能也是非常大的。所以如果要比较细胞性质
差异,我们建议用同一仪器实验和同一种模型进行
数据分析,只有在相同条件下所得到的结果才能更
好反映细胞之间弹性的差异。
[参 考 文 献]
[1] Yang H, Yu J, Fu G, et al. Interaction between single mol-
ecules of Mac-1 and ICAM-1 in living cells: an atomic force
microscopy study. Exp Cell Res, 2007, 313(16): 3497-504
[2] Bálint Z, Krizbai IA, Wilhelm I, et al. Changes induced by
hyperosmotic mannitol in cerebral endothelial cells: an atomic
force microscopic study. Eur Biophys J, 2007, 36(2): 113-
20
[3] Kasas S, Wang X, Hirling H, et al. Superficial and deep
changes of cellular mechanical properties following cytosk-
eleton disassembly. Cell Motil Cytoskeleton, 2005, 62(2):
124-32
[4] Cross SE, Jin YS, Rao JY, et al. Nanomechanical analysis of
cells from cancer patients. Nat Nanotechnol, 2007, 2(12):
780-3
[5] Chouinard JA, Grenier G, Khalil A, et al. Oxidized-LDL
induce morphological changes and increase stiffness of endot-
helial cells. Exp Cell Res, 2008, 314(16): 3007-16
[6] Sirghi L, Ponti J, Broggi F, et al. Probing elasticity and adhe-
sion of live cells by atomic force microscopy indentation.
Eur Biophys J, 2008, 37(6): 935-45
[7] Hertz H. über den kontakt elastischer korper. J Reome
Amgew Mathematik, 1882, 92: 156-71
[8] Zhao MH, Charudharshini S, Diane JB, et al. Rate- and
depth-dependent nanomechanical behavior of individual liv-
ing Chinese hamster ovary cells probed by atomic force
microscopy. J Mater Res, 2006, 21(8): 1906-12
[9] Costa KD. Single-cell elastography: probing for disease with
the atomic force microscope. Dis Markers, 2003, 2(3): 139-
54
[10] Radmacher M, Fritz M, Kacher CM, et al. Measuring the
viscoelastic properties of human platelets with the atomic
force microscope. Biophys J, 1996, 70(1): 556-67
[11] Carl P, Schillers H. Elasticity measurement of living cells
with an atomic force microscope: data acquisition and
processing. Eur J Physiol, 2008, 457(2): 551-9
[12] Tamara KB, Craig DW, Igor S. Human epithelial cells in-
crease their rigidity with ageing in vitro: direct measurements.
Phys Med Biol, 2005, 50(1): 81-92
[13] Charles RF, Gisouvander G, Paolo DLR, et al. Elastic mem-
brane heterogeneity of living cells revealed by stiff nanoscale
membrane domains. Biophys J, 2008, 94(4): 1521-32
[14] Li QS, Lee GY, Ong CN, et al. AFM indentation study of
breast cancer cell. Biochem Biophys Res Commun 2008,
374 (4): 609-13
[15] Oliver WC, Pharr GM. An improved technique for deter-
mining hardness and elastic modulus using load and dis-
placement sensing indentation experiments. J Mater Res,
1992, 7(6): 1564-83
[16] Bremmell KE, Evans A, Prestidge CA. Deformation and
nano-rheology of red blood cells: an AFM investigation.
Colloids Surf B: Biointerfaces, 2006, 50(1):43-8
[17] Berry MF. Mesenchymal stem cell injection after myocar-
dial infarction improves myocardial compliance.Am J Physiol
Heart Circ Physiol, 2006, 290 (6): 2196-203
[18] Yoshikawa Y, Yasuike T, Yagi A, et al. Transverse elasticity
of myofibrils of rabbit skeleta muscle studied by atomic
force microscopy. Biochem Biophy Res Commun, 1999,
256 (1): 13-9
[19] Engler AJ. Myotubes differentiate optimally on substrates
with tissue-like stiffness: pathological implications for soft
or stiff microenvironments. J Cell Biol, 2004, 166 (6): 877-
87
[20] Lieber SC, Aubry N, Pain J, et al. Aging increases stiffness
of cardiac myocytes measured by atomic force microscopy
nanoindentation. Am J Physiol Heart Circ Physiol, 2004,
287(2): H645-51
图6 在三种加载速率下,分别用上述三种方法所得到
的杨氏模量E