空气辅助雾化的细胞递送参数的优化

作者：严心涛王策王耀宋飞飞武晓东来源：《光学精密工程》日期：2022-09-22人气：588

大面积皮肤伤口与多种突发事件、重大疾病紧密相关，是临床上常见病症，其迁延不愈的特性严重降低了患者的生活质量，甚至危及生命，并且在全球范围内造成了巨大的经济负担^［1-2］。由大面积皮肤创伤而引发的一系列严重并发症和感染是威胁患者生命的两个重要因素^［3］。因此，除了给患者静脉注射液体和营养素以抵消脱水和补充失去的蛋白质外，采取有效的治疗方法以减轻患者疼痛和挽回生命显得至关重要。目前临床实践的有效治疗方法主要为自体皮肤移植术^［4-5］，但对于自体皮源本已稀缺的患者，该方法仍面临着巨大挑战。

细胞移植术是解决这个问题的有效方法^［6-8］，即从患者身上采集自体皮肤细胞（通常是角质形成细胞）后直接递送到伤口，或经过细胞培养传代后悬浮，再递送到伤口。事实证明，细胞移植术可以显著影响细胞在伤口中的分布和增殖^［9-11］。近些年，细胞移植术已经从静脉注射、肌肉注射、注射器或移液枪的局部滴/注^［12-13］发展到雾化式喷洒^［14-16］。尽管静脉注射、肌肉注射或局部细胞给药均显示出治疗效果，但前两种给药方法会导致作用于伤口部位的细胞浓度低，治疗效果不佳^［12］；通过移液枪或注射器进行局部细胞给药虽可减少患者疼痛，但细胞递送效率低，且分布不均匀^［13］。细胞悬浮液喷雾自体移植用自体表皮衍生细胞，使用更小的自体部位（供体部位与烧伤面积比例ε可达1∶100）即可在短时间内完成再上皮化^［16-17］，对于深部、大面积和不规则创面，该方法能实现细胞的快速和高效的递送。目前涉及组织修复与再生应用的空气喷枪或雾化装置仍存在一些缺陷，如生成的雾滴大小不一致、分布不均匀^［18］，雾化参数对雾化特性的精确控制尚不清楚^［9］，在辅助空气压力为69 kPa下，细胞活化率仅为65%^［19］，有相关文献报道^［20］对细胞雾化递送4天后，细胞增殖能力大受影响，细胞活化率降为45%。如何实现药物的均匀化雾化、不同面积的雾化或较高的细胞活化率雾化是细胞给药雾化递送装置开发中需面临的几个关键技术问题。

空气辅助雾化方法具有仅需较低空气压力和射流压力即可实现较好的雾化效果，以及较大的射流通道可避免堵塞等特点^［21］，本文构建了基于该雾化方法的可用于细胞递送的雾化装置，通过对其雾化特性测试方法研究，侧重分析雾化高度、辅助气体流量和射流流量等雾化参数对雾化特性和生物墨水（Bio-ink）递送的影响，以指导雾化参数的优化设计和提高生物墨水的雾化喷涂效率。

2 装置设计及试剂材料

2.1　装置原理及设计

空气辅助雾化的细胞递送装置（Air-Assisted Atomization Device，AAAD）的工作原理如图1（a）所示，空气泵输出的气体通过空气滤芯过滤后分两个流路，一路经气压调节阀后在雾化喷嘴出口处形成稳定的高速涡旋气流，另一路也经气压调节阀以稳定压力驱动无菌储液管中细胞悬浮液，并在喷嘴的针管出口形成液柱或液膜。在喷嘴出口处，气体与生物墨水之间会形成很高的相对速度，从而破坏液柱或液膜，改变生物墨水的流动性、表观粘性和表面张力，最终使连续相的生物墨水变成分散相以达到雾化效果。其中，两个流路中气压调节阀分别用于调节辅助空气气压以控制气体流量和调节无菌储液管的驱动压力以控制射流流量。为分析雾化高度对雾化性能的影响，雾化喷嘴可沿导杆做上下滑动。

图1 AAAD原理及装置

Fig.1 Principle and setup for air-assisted atomization device

AAAD主要包含气动控制模块和雾化喷嘴模块两部分，实物分别如图1（b）和图1（c）所示。装置主要组成部件包含流量为10 L/min的微型空气压缩泵（中国，SKOOCOM），100 kPa与200 kPa量程且灵敏度均优于0.2%的精密调压阀（日本，SMC），储液筒（日本，MUSASHI），27G雾化针管（日本，MUSASHI），雾化针管的内径为200~250 μm，喷嘴的内径为0.7~1.0 mm。采用量程为5 L/min且精度优于1.5%的微型空气流量计（中国，SENLOD）和量程为40 mL/min且精度优于5%的微型液体流量计（瑞士，SENSIRION）来测量辅助气体流量和射流流量。

2.2　试剂材料

人永生化表皮细胞（HaCaT，中国，COBIOER）用含10%胎牛血清（FBS，中国，TIANHANG）和1%青霉素/链霉素（A5955，Sigma）的高糖培养基（DMEM，Gibco）在T75 cm培养瓶（美国，CORNING）中培养，孵育条件为37 ℃、CO₂的体积分数为5%。每4天传代一次，每隔一天更换一次培养基。当细胞达到90%以上汇合时，用胰蛋白酶（15090046，Gibco）收集细胞，将细胞悬液在室温下以1 000 r/min的转速离心5 min，弃去上清液，再将细胞重悬于无菌生理盐水（ST341，中国，Beyotime）中，配制成细胞密度约为10⁶ cells/mL的生物墨水。在雾化基本特性测试中，均以无菌生理盐水为研究对象。

3 雾化特性测试方法

3.1　雾滴粒径分布及均匀度分析

本文采用激光粒度仪（JL-3000，中国，JNGX）对雾滴粒径进行测量。该仪器主要基于夫琅和费（Fraunhofer）衍射理论测量雾滴尺寸及对应的雾滴个数。雾化喷嘴产生的雾滴是由大小不等的雾滴群组成，本文雾滴粒径分布采用体积累积分布表示，平均粒径采用体积平均粒径D［4，3］计算，即：

（1）

经换算可得：

（2）

其中： $n$ 是各个单一粒径对应的雾滴个数； $D$ 是雾滴的粒径，单位为μm； $f$ 是各单一粒径的雾滴群对应的体积频度。

通常，雾滴粒径是不连续的，但当雾滴数量足够多时，可假设雾滴的粒径累计分布是连续的，通常其符合Rosin-Rammler分布（简称“R-R”分布）规律^［22］，即：

（3）

其中： $F (D_{i})$ 表示雾滴粒径小于 $D_{i}$ 的体积累计百分比； $\bar{D}$ 为雾滴的特征尺寸，即 $F (D_{i}) = 0.632$ 对应的雾滴粒径； $k$ 为均匀度指数， $k$ 越大，表示雾滴粒径分布的均匀性越好。对于符合“R-R”分布规律的雾化，若粒径 $D_{i}$ 和 $D_{j}$ 对应的累计百分比分别为 $F (D_{i})$ 和 $F (D_{j})$ ，则由式（3）不难推导出雾化均匀度指数 $k [D_{i}, D_{j}]$ ，即：

（4）

3.2　雾化角测试及雾化面积评估

雾化过程中雾滴群在空间呈圆锥状分布，该圆锥的张角即为雾化角，其大小直接影响雾化面积A。本文采用工业相机（DMK 23G445，德国，THE IMAGING SOURCE）对雾化状态摄像后，在Matlab（美国，MathWorks）中计算雾化角。使用宽工作距离范围的物镜（BLADE 2.4/12， Korea， FIM optcis），摄像分辨率为640×480像素，其中每个像素对应物方尺寸为139 μm。

由于单帧图像信噪比（Signal-to-Noise Ratio，SNR）较低，边界提取不准确，而随着图像叠加数量的增加，SNR也将增大。经测试，单帧图像的SNR仅为11.8 dB，用于雾化边界提取误差较大，通过对多帧图像进行叠加，可显著提高图像的SNR。当叠加图像达到227帧时，图像SNR为50 dB，对叠加后的高SNR图像进行二值化处理，能够有效保留雾化边界，且边界更加清晰。随着图像的叠加数量增多，SNR将进一步增大，如当叠加视频中所有图像时（总帧数为238），图像SNR可达65.5 dB。本文计算雾化角的方法是：先对雾化过程进行连续视频拍摄，分析时，叠加视频的所有帧的图像，以提高图像SNR，如图2（a）为单帧图像，图2（b）为视频叠加图像；对雾化角图像进行二值化，采用局部自适应阈值分割算法，该算法可根据局部光强分布自动调整二值化阈值，对于边界的识别更加准确和稳定，如图2（c）；图2（d）为二值化图像的边界数据点，考虑到靠近喷嘴和远离喷嘴处的雾化角度有所变化，因此数据点的分析也分为两部分，采用最小二乘法拟合数据点，得到边界直线方程，进而计算得到靠近喷嘴的雾化角 $θ_{1}$ 和远离喷嘴处的雾化角 $θ_{2}$ 。如图2（d）所示，计算出雾化角 $θ_{1}$ 和 $θ_{2}$ 后，可进一步估算出生物墨水喷洒在基板上的面积A，即：

图2 雾化角测试方法

Fig.2 Measurement method of atomization angle

（5）

其中： $h_{N}$ 为最小二乘法拟合直线断点位置距离喷嘴出口位置的高度， $h_{F}$ 为雾化基板距离喷嘴出口位置的高度，单位均为mm。

3.3　雾滴沉降速度

由于雾滴沉降速度会影响细胞活性^［23］，本文利用高速相机（HS1，德国，PCO）的双快门捕获功能来测算雾滴沉降速度，相机单次曝光时间设置为10 μs，拍摄时间间隔T为3.58 μs。采用10×远距离物镜（375-039，日本，三丰），视场为1.65×2.75 mm²。由于曝光时间极短，拍摄图像SNR较低，如图3（a），很难对视场中的雾滴进行准确定位。本文采用BM3D去噪算法，可有效提高图像质量，如图3（b）；利用二值化算法对图像中的雾滴进行识别后，如图3（c），根据雾滴光强分布计算得到雾滴质心坐标 $(x_{T 1}, z_{T 1})$ 和 $(x_{T 2}, z_{T 2})$ ，高速相机双快门捕获图像处理前与处理后分别如图3（d）和图3（e）所示；再根据相邻两帧图像中雾滴的坐标，分别计算出雾滴的横向运动速度 $v_{x} = |x_{T 2} - x_{T 1}| / T$ ，纵向运动速度 $v_{z} = |z_{T 2} - z_{T 1}| / T$ ，以及合速度 $v_{r} = \sqrt[]{{v_{x}}^{2} + {v_{z}}^{2}}$ 。

图3 雾滴沉降速度测试方法

Fig.3 Measurement method of droplet settling velocity

3.4　细胞活化率测试

将配制好的生物墨水一部分用于无喷涂的阴性对照，采用注射器（2 mL，中国，KINDLY）挤压递送；另一部分用AAAD喷洒。为了评估雾化后活细胞数与死细胞数比例，采用Calcein/PI细胞活性与细胞毒性检测试剂盒（C2015M，中国，Beyotime）在37 ℃避光孵育30 min，再用倒置荧光显微镜（IX73，日本，OLYMPUS）拍摄荧光图片，并利用ImageJ（美国，National Institutes of Health）对荧光图像进行自动计数和Merge处理。细胞活化率 $η$ 由活细胞总数占细胞总数（含活细胞与死细胞）的百分比来表示。采用CCK-8试剂盒（C0043，中国，Beyotime）检测喷涂或未喷涂的细胞在24 h、48 h及72 h后的活性，来分析AAAD喷涂对细胞增殖性能的影响。

公开的文献［15］中已分析了辅助空气压力 $P_{A}$ 对 $η$ 的影响，即在 $h_{F}$ 为 $100 m m$ ， $P_{A} \leq 70 k P a$ 时， $P_{A}$ 对 $η$ 无显著差异。在本文中，重点分析 $h_{F}$ 和辅助空气流量 $Q_{A}$ 对 $η$ 的影响。

3.5　统计方法

数据表示为平均值±标准偏差（Standard Deviation，SD），并使用Origin（美国，OriginLab）做非配对t检验（Unpaired t test）、单因素方差分析（One-Way ANOVA）或双因素方差分析（Two-way ANOVA），以评估各实验组之间的差异性。

4 结果与讨论

4.1　射流流量和雾化高度对雾滴粒径分布的影响

图4为雾滴粒径分布及雾化参数对其影响。 $Q_{A}$ 设定为3 L/min，在不同 $Q_{L}$ 和不同 $h_{F}$ 的雾化条件下，经激光粒度仪得到不同雾滴粒径及对应的体积频度 $f$ ，如图4（a）所示，由式（2）可计算出D［4，3］，根据式（3）可以得到“R-R”分布曲线，在该曲线上分别选取25%和75%累计百分比对应下的粒径大小D25和D75，一般选取D25和D75之间的分布情况来说明颗粒的集中情况^［15］，利用式（4）计算出 $k [D 25, D 75]$ 来定性分析雾滴的均匀度程度。

图4 雾滴粒径分布及雾化参数对其影响，其中，H1~H4分别表示雾化高度为25、50、75及100 mm

Fig.4 Droplet size distribution and the influence of spray parameters on it，H1-H4 represented the spray height of 25， 50， 75 and 100 mm， respectively

通过双因素方差分析，可得出 $h_{F}$ 、 $Q_{L}$ 以及两者间的交互作用均对D［4，3］有显著影响。利用单因素方差分析来具体探究它们之间的关系，如图4（b），不难得到：（1）在同等 $h_{F}$ 下， $Q_{L}$ 对D［4，3］有显著影响，表现为 $Q_{L}$ 越大，D［4，3］越大；（2） $Q_{L} = 0.2 m L / m i n$ 时， $h_{F}$ 对D［4，3］无显著影响；（3）以高度为25 mm实验组为参考组，在相同 $Q_{L}$ 条件下，50 mm实验组与参考组差异不大，但75 mm实验组和100 mm实验组与参考组相比均显著增大，这说明在雾化过程中，当 $h_{F}$ 大于一定值后，随着其进一步增大，雾滴与雾滴间存在融合现象，使得D［4，3］有显著增大趋势，但当 $Q_{L}$ 过低时， $h_{F}$ 对D［4，3］无显著影响。

封装细胞的雾滴与基板发生撞击中可能造成细胞损失，一般而言，D［4，3］越小，细胞活化率 $η$ 越低^［23］；但D［4，3］过大时，细胞位于雾滴边缘概率会增加，雾滴与基板撞击过程中，细胞受到的剪切应力会增大， $η$ 降低^［24］。因此， $h_{F}$ 与 $Q_{L}$ 在细胞递送中是重要雾化参数。细胞传送装置在实际使用中 $Q_{L}$ 一般低于5 mL/min^{［14，16，25］}，否则会影响生物墨水的利用率。

经测量，当 $Q_{A} = 3 L / m i n$ 时， $P_{A} \approx 20 k P a$ 。利用单因素方差分析可以得到：在H2实验组和H4实验组中， $Q_{L}$ 对 $k [D 25, D 75]$ 有显著不同的影响，但影响程度不一致，如图4（c）所示，说明这两组实验的雾滴生成状态不稳定；但在H1~H4四组实验中，当 $Q_{L}$ 相同时， $h_{F}$ 对 $k [D 25, D 75]$ 影响不大。为进一步分析喷嘴出口的辅助空气流场对雾化的影响，研究了20~100 kPa辅助气体压力下的雾滴分布状态，如图4（d）所示，可得到：在 $P_{A} \leq 50 k P a$ 时，各实验组的 $k [D 25, D 75]$ 相差不显著；而在 $P_{A} \geq 60 k P a$ 时，各实验组的 $k [D 25, D 75]$ 显著提高，尤其在80 kPa和90 kPa实验组中，相对于前一个梯度实验组而言， $k [D 25, D 75]$ 提高地非常显著。虽然 $P_{A}$ 的增大能提高 $k [D 25, D 75]$ ，但同样会影响细胞递送中 $η$ 值，在文献［15］中已做了相关分析。因此，在一些生物材料需均匀化喷涂而不考虑 $η$ 指标时，可以重点关注 $P_{A}$ 这项雾化参数。

4.2　雾化角及雾化面积的影响因素

图5为雾化参数对雾化角及雾化面积的影响。经最小二乘法拟合边界直线方程，在不同 $Q_{L}$ 下可计算得到 $θ_{1}$ 和 $θ_{2}$ ，且 $h_{N} \approx 25 m m$ 。从图5（a）中可看到：在 $Q_{A}$ 为 $3 L / m i n$ 的实验组中，随着 $Q_{L}$ 增大， $θ_{1}$ 先从26.44°增加到30.81°之后一直递减到8.04°，而 $θ_{2}$ 在15.55°~25.26°间上下波动；当 $Q_{L} \leq 7 m L / m i n$ 时， $θ_{1} > θ_{2}$ ；而在 $Q_{L}$ 为8 mL/min或9 mL/min时， $θ_{1} < θ_{2}$ 。结果表明，在 $Q_{A}$ 恒定条件下， $Q_{L}$ 对 $θ_{1}$ 的影响比对 $θ_{2}$ 的影响显著。利用式（5）可计算出在不同 $Q_{L}$ 下的雾化面积 $A$ ，其中取 $h_{N} = 25 m m$ 。不难发现在固定 $h_{F}$ 下调节 $Q_{L}$ 时，在 $Q_{L} = 3 m L / m i n$ 条件下的 $A$ 最大，且 $h_{F}$ 越大，差异越显著；在 $h_{F} = 25 m m$ 时，雾化面积可控制在100 mm²内，可满足小面积喷洒的工况，如图5（b）所示。

图5 雾化参数对雾化角及雾化面积的影响

Fig.5 Effect of atomization parameters on spray cone angle and spray area

如图5（c），在 $Q_{L}$ 为为 $3 m L / m i n$ 的实验组中，当 $Q_{A} \leq 1.5 L / m i n$ 时，雾化不充分，且 $θ_{1}$ 略小于 $θ_{2}$ ；而当 $Q_{A} \geq 2 L / m i n$ 时，雾化效果良好， $θ_{1}$ 和 $θ_{2}$ 会显著增大，且趋于平稳， $θ_{1} > θ_{2}$ 。结果表明，在 $Q_{L}$ 恒定条件下，为形成良好的雾化效果， $Q_{A}$ 需足够大，建议不低于2 L/min。其中，部分雾化图像如图5（d）所示。

4.3　气体流量和雾化高度对雾滴沉降速度的影响

图6（a）是雾滴沉降速度测试平台；在 $Q_{L}$ 为 $3 m L / m i n$ ， $Q_{A} \leq 5 L / m i n$ 时，如图6（b）所示，H1实验组的 $v_{x}$ 最高，随着雾化高度 $h_{F}$ 变大， $v_{x}$ 会降低，但H3与H4实验组的 $v_{x}$ 变化不大；在各个实验组中， $Q_{A}$ 对 $v_{x}$ 影响不大；由图6（c）和图6（d）可知， $v_{z}$ 和 $v_{r}$ 均比 $v_{x}$ 大的多；在4个实验组中，H2实验组的 $v_{z}$ 和 $v_{r}$ 最大；在各个实验组中， $Q_{A}$ 对 $v_{z}$ 和 $v_{r}$ 影响也不显著。

图6 雾滴沉降速度测试及影响分析

Fig.6 Measurement and influence analysis of drop velocities

结果表明，在辅助空气低流量下，雾滴的运动速度受 $h_{F}$ 的影响较大；辅助气体流场主要在喷嘴出口的小范围内影响 $v_{x}$ ，随着 $h_{F}$ 变大，雾滴在空气阻力作用下， $v_{x}$ 会逐渐减小；雾滴的 $v_{z}$ 和 $v_{r}$ 在辅助气体的作用下，有一段加速积累的过程，需运行一定距离后才会逐渐降低。雾滴的运动速度直接影响其与雾化基板之间的剪切应力，为保证细胞雾化递送中的细胞活化率 $η$ ，雾滴的运动速度应尽可能低。因此，对于本文设计的雾化喷嘴，应避免在50 mm的雾化高度喷涂。

4.4　生物墨水雾化的细胞活化率影响因素

经倒置荧光显微镜拍摄和Image J软件分析后，可得出阴性对照组和不同 $h_{F}$ 下实验组中 $η$ ，如图7（a）。经单因素方差分析，得到 $h_{F}$ 对 $η$ 有显著的影响（p<0.001）；将各个实验组与阴性对照组作非配对t检验，可得到在50 mm的雾化高度进行细胞传递时，确实会显著降低 $η$ （p<0.001），如图7（b）所示，该结果同时证明了 $η$ 与雾滴粒径分布和雾滴沉降速度均存在一定关联。图7（c）结果表明，当 $Q_{A} \leq 3 L / m i n$ 下， $Q_{A}$ 对 $η$ 无显著影响；而在 $Q_{A} \geq 4 L / m i n$ 下， $Q_{A}$ 对 $η$ 会有显著负作用。采用CCK-8试剂盒对HaCaT细胞在24 h、48 h以及72 h后增殖检测的结果显示，AAAD喷涂对细胞增殖性能无显著影响，如图8所示。

图7 AAAD喷涂的HaCaT细胞活性

Fig.7 Viability of cells sprayed by AAAD

图8 细胞增殖不受AAAD喷涂的影响

Fig.8 Proliferation of cells was not impaired by AAAD sprayed

5 结论

本文构建了一种基于空气辅助雾化技术的可用于细胞递送的雾化装置；系统地研究了雾滴粒径分布、雾化角、雾滴沉降速度等雾化特性的光学测试方法、图像优化处理算法以及雾滴均匀度和雾化面积的定量评估方法；并利用统计学方法分析了多种雾化参数对雾化特性或细胞活化率的影响。结果表明：高于60 kPa的辅助气体压力与雾滴均匀度呈正相关关系，但会降低细胞活化率；通过控制射流流量和雾化高度，可实现50~1 800 mm²大宽幅的雾化面积调节；雾滴的运动速度受雾化高度影响较大，且在50 mm附近的运动速度最大，均值可达14 m/s；雾化高度显著影响着雾滴的粒径分布和运行速度，进而影响到细胞活化率，结果显示雾化高度为50 mm时对细胞活化率有非常显著影响，阴性对照细胞活化率均值为92.98%，在50 mm雾化高度时其仅为64.01%，而在100 mm雾化高度时其可高达90.24%。本文研制的雾化装置为细胞或生物材料提供了一种递送手段，雾化特性测试方法的研究也对雾化参数的优化设计和细胞或生物材料的喷涂效率提升具有指导意义。此外，雾化装置所需气源的输出流量不超过10 L/min即可实现较好的雾化效果，这为未来开发一体化手持式的细胞雾化喷枪创造了有利条件，在细胞移植术领域将具有广泛的应用前景。

关键字：优秀论文

上一篇：压电定位平台Hammerstein建模与反馈线性化控制
下一篇：面向可见光和SAR影像配准的特征点检测

栏目分类

热门排行

推荐信息

期刊知识

空气辅助雾化的细胞递送参数的优化

2 装置设计及试剂材料

2.1 装置原理及设计

2.2 试剂材料

3 雾化特性测试方法

3.1 雾滴粒径分布及均匀度分析

3.2 雾化角测试及雾化面积评估

3.3 雾滴沉降速度

3.4 细胞活化率测试

3.5 统计方法