下载:
-
纳米级别构筑单元的尺寸效应赋予纳米粒子材料独特的电学、光学、磁学和力学性能[1]。自组装是一种“自下而上”制备纳米有序结构材料的有效手段。根据Watson-Crick原则，互补DNA序列可以通过氢键形成热力学可逆的碱基对，使DNA成为纳米粒子间理想的“黏合剂”[2]。Mirkin等[3]于1996年将末端带有硫醇基团的DNA链与球形金纳米粒子混合，首次制备了表面DNA链均匀分布的纳米粒子，即DNA均匀功能化纳米粒子。此后，研究人员利用DNA功能化纳米粒子的自组装行为得到了线性、薄膜状、超晶格等多种超结构[4-6]。
DNA功能化纳米粒子的自组装行为与原子间形成共价键的过程有很多相似之处，因而DNA功能化纳米粒子被称为可编程原子等价物[7]。但是，DNA均匀功能化纳米粒子难以控制其形成杂化碱基对的数量和方向，在制备具有特殊结构的纳米粒子材料方面具有一定局限性。近来，仅在表面特定位置分布DNA链的纳米粒子，即DNA非均匀功能化纳米粒子为精确控制自组装结构提供了新方法[8-11]，纳米粒子的结构可设计性和序列可编程性使其自组装行为更加复杂。定量描述DNA非均匀功能化纳米粒子的自组装动力学行为对获得结构可控的自组装体具有十分重要的意义。
计算机模拟在示踪粒子运动轨迹、阐明实验机理等方面具有极大优势，已成为研究纳米粒子自组装行为和结构的重要方法。本文以DNA非均匀功能化纳米粒子为研究对象，通过粗粒化分子动力学模拟纳米粒子可编程自组装过程，探索自组装的结构与动力学。
-
基于Travesset课题组[12-13]发展的DNA均匀功能化纳米粒子的粗粒化模型，本文建立了DNA非均匀功能化纳米粒子模型（如图1所示）。设球形纳米粒子的半径（Rs）为3σ，其中σ为模拟的长度约化单位（σ =2.0 nm）。3条DNA链等间距地锚定于球面赤道位置。每条DNA链含有中性部分（Spacer）和DNA黏性端（Linker），分别含有ns和nl个粗粒化珠子。黏性端中1个粗粒化珠子代表3～7个碱基；若nl=3时，粗粒化DNA链的黏性端相当于序列长度约为15个碱基的寡聚核酸。模型中，DNA链的尺寸与纳米粒子直径相近[1, 14, 15]。为了确保DNA链完全互补，黏性端序列分别设置为-ACK和-FGT。需要说明的是A-T和C-G为互补碱基对，而K-F互补碱基对的引入是为了防止DNA链自杂化。含有-ACK序列的功能化纳米粒子标记为V型纳米粒子；相应地，含有-FGT序列的纳米粒子标记为W型纳米粒子。
图 1 （a）DNA非均匀功能化纳米粒子的粗粒化模型；（b）互补DNA链的黏性端杂化示意图
Figure 1. （a）Coarse-grained model of non-uniformly DNA-functionalized nanoparticles;（b）Hybridization between complementary DNA strands containing ns spacer beads and nl linker beads
模型中，相邻珠子的键连作用为谐振子势，包括键长作用势（Ubond）与键角作用势（Uangle）
$ \begin{array}{l} {U_{\rm{bond}}}\left( R \right) = \dfrac{1}{2}{k_{\rm{s}}}{\left( {R - {R_0}} \right)^2}; \quad{U_{\rm{angle}}}\left( \theta \right) = \dfrac{1}{2}{k_\theta }{\left( {\theta - {\theta _0}} \right)^2} \end{array} $ 式中，R表示2个粗粒化珠子的间距，其平衡值R0 = 0.84σ，系数ks = 330ε/σ2（ε为能量约化单位）。θ表示3个连续粗粒化珠子之间的角度，其平衡值θ0 = π，系数kθ = 30ε。互补碱基对粗粒化珠子间的非键相互作用由Lennard-Jones势能表示：
$U\left( r \right) = 4{\varepsilon _{\rm{bp}}}\left[ {{{\left( {\dfrac{\sigma }{R}} \right)}^{12}} - {{\left( {\dfrac{\sigma }{R}} \right)}^6}} \right] $ 。式中，作用强度εbp = 10ε，截断半径Rc = 2.75σ。其他粗粒化珠子之间的相互作用为纯排斥力，由Weeks-Chandler-Andersen势能表示[16]：$U\left( r \right) = 4\varepsilon {\left( {\dfrac{\sigma }{R}} \right)^{12}} $ ，此处Rc=21/6σ。分子动力学模拟在正则系综（NVT）下进行。体系一共含有N=200个DNA非均匀功能化纳米粒子。V和W型纳米粒子的个数分别为NV和NW（N =NV+NW），化学计量比（f）定义为f =NV/NW，f≥1。模拟盒子满足周期性边界条件，其尺寸为L×L×L。DNA功能化纳米粒子的体积分数（
$\varphi $ ）定义为$\varphi \equiv 4\pi N{\left( {{R_{\rm{s}}} + {R_{\rm{e}}}} \right)^3}/3{L^3}$ ，其中Re是DNA链特征尺寸，${R_{\rm{e}}} = \sigma {\left( {{n_{\rm{s}}} + {n_{\rm{l}}}} \right)^{1/2}}$ 。分子动力学模拟的时间步长（t）设置为t = 0.001τ，其中τ为时间约化单位。为了确保体系达到平衡状态，每次分子动力学模拟约进行4×108步。 -
在分子聚合反应中，多官能度单体聚合可以得到非线性拓扑的聚合物，如支化和网络状聚合物。通过调控单体的化学计量比可以得到不同的结构。类比于多官能度单体的分子聚合反应，当纳米粒子功能化3条DNA链时，每一个通过DNA序列杂化加入到聚集体中的纳米粒子将带来2个新的“生长点”（DNA黏性端）；在此驱动下，游离的纳米粒子自组装形成支化和三维网络状超结构[17]。同高分子链拓扑结构的调控机制类似，自组装超结构的拓扑依赖于V和W型纳米粒子的化学计量比f。
为了证实上述推测，构建了不同数目的V和W型纳米粒子体系，并进行了分子动力学模拟。V和W型纳米粒子的数目之比（NV/NW）分别为100/100、125/75、150/50和170/30，对应的化学计量比分别为1.0、1.7、3.0和5.7。图2展示了不同化学计量比之下，纳米粒子通过互补DNA序列杂化所形成的自组装超结构。当f=1.0时，在互补DNA序列杂化的引导下，纳米粒子自组装形成三维网络状超结构。需要说明的是，此网络状超结构在各个方向上是贯穿的，为逾渗网络。随着f的增大，过量的DNA序列参与杂化的可能性减小，网络状超结构变稀疏（f =1.7），进而转变为支化超结构（f =3.0）。当f=5.7时，数量较少的W型纳米粒子完全被V型纳米粒子包围，只能够形成离散分布的小聚集体；大部分V型纳米粒子仍处于游离态。由此可见，通过改变体系的化学计量比，DNA链非均匀功能化纳米粒子从三维网络状超结构转变为支化超结构，甚至是离散分布的小聚集体。
图 2 化学计量比f对自组装超结构的影响（图中只显示通过互补DNA序列杂化所形成的最大纳米粒子超结构，红线表示由互补DNA序列杂化形成的纳米粒子间连接）
Figure 2. Effect of the stoichiometric ratio on the superstructure of nanoparticles（Note that only the largest superstructures are shown; Red lines represent the connectivity between various nanoparticles through the hybridization of complementary DNA strands）
为了定量表征纳米粒子自组装结构，我们构建了自组装结构的几何模型（图3（a））。在分子动力学模型中，DNA链有一定刚性，并且等间距地锚定于球面赤道位置。通过互补DNA序列杂化（实线），最邻近纳米粒子间距离为r0，即为典型“键合”距离[18]。次邻近纳米粒子间距离为
$\sqrt 3 {r_0}$ （虚线）。纳米粒子自组装结构特征可通过粒子间径向分布函数g（r）表征，其中径向分布函数特征峰的位置和强度可以表征自组装超结构中纳米粒子之间的相对位置及其分布[19-22]。
图 3 （a）纳米粒子自组装结构的几何模型（实线为通过互补DNA序列杂化所形成的纳米粒子间连接，虚线连接的是次邻近纳米粒子）；（b～d）不同化学计量比（f）和中性部分长度（ns）下纳米粒子径向分布函数g（r）（约化温度T *= 0.90，纳米粒子体积分数φ = 0.10）
Figure 3. （a）Geometrical representation of self-assembled superstructures of nanoparticles（Solid lines denote the connectivity between nanoparticles, and dashed lines denote the connectivity of nanoparticles with sub-nearest distance）；（b—d）Radial distribution function g（r）of nanoparticles for a range of stoichiometric ratios （f） and the lengths of spacer（ns）（The reduced temperature is set as T *= 0.90, and the packing fraction is set as φ = 0.10）
图3（b～d）展示了不同化学计量比f和中性部分长度（即中性珠子个数，ns）下的径向分布函数g（r）。图中虚线从左至右分别对应于特征峰的理论预测位置r/r0为1和
$\sqrt 3 $ 。当f=1.0时，带有互补序列的纳米粒子相互连接形成三维网络状结构，纳米粒子间相对位置关系符合图3（a）的预测（如图3（b）所示，在r/r0为1和$\sqrt 3 $ 处特征峰的峰值较大，峰形较为显著）。随着f的增大，体系中存在多个松散分布的支化聚集体，使得特征峰强度减弱（图3（c））。当ns增加时，径向分布函数g（r）的特征峰进一步减弱（图3（d））。此外，当ns=4时，径向分布函数g（r）在r/r0 = 1之前没有明显的峰（图3（c））；而当ns = 10时，径向分布函数g（r）在r/r0 = 1之前出现多个明显的峰（图3（d））。其原因是，较长DNA链占据的空间较大，使得相连的V和W型粒子间空隙变大，从而形成互穿支化结构[18]。 -
化学计量比不仅影响自组装超结构，还会改变自组装动力学。为了定量描述自组装过程中聚集体的生长过程和DNA黏性端的杂化程度，模拟中统计如下2个物理量：聚集体数量（Na），它表示包含2个及以上纳米粒子的自组装超结构数量；碱基对数量（Nh），它表示通过互补DNA序列杂化生成的碱基对数量。图4是聚集体数量和碱基对数量随时间的变化曲线。在给出的化学计量比之下，Nh随着模拟时间的延长而增加，最终达到平衡值。Na则随模拟时间的延长表现出非单调变化。在f =1.0、1.7、3.0时，Na先增加后减小，直至达到平衡值。这说明在自组装早期，纳米粒子先形成较小的聚集体；随后，在聚集体外围DNA黏性端的引导下，小聚集体之间相互融合形成更大的聚集体，甚至网络状结构，导致聚集体数量下降。当f =5.7时，过量的V型纳米粒子分布在聚集体外围，消耗并屏蔽了非过量的W型纳米粒子，使得聚集体数量变化不再显著。
图 4 不同化学计量比下聚集体数量（Na）与DNA杂化碱基对数量（Nh）随时间的变化曲线（约化温度T *= 0.90，纳米粒子体积分数φ = 0.10）
Figure 4. Temporal evolution of numbers of aggregates （ Na） and numbers of hybridizations （Nh） at various stoichiometric ratios （The reduced temperature is set as T *= 0.90, and the packing fraction φ = 0.10）
化学计量比还影响了自组装结构的动力学行为。本文考察了纳米粒子的均方位移（MSD）以及最大聚集体所包含的纳米粒子数目与纳米粒子总数之比（η）。图5（a）对比了不同化学计量比下的均方位移随时间的变化趋势[12, 23]。图5（b）为自组装体系中η随时间的变化曲线。当f =5.7时，体系中只能形成非常小的聚集体，同时存在未发生DNA序列杂化的纳米粒子；纳米粒子的MSD接近于线性增长。随着f的减小，MSD曲线在经历自组装初期的快速增长后，增幅出现不同程度的减弱。特别是，在f =1.0时，自组装时间t超过5×104τ之后，MSD达到平衡，不再增加。由图5（b）可知，此时纳米粒子自组装形成一个大的聚集体（网络状结构），导致纳米粒子扩散缓慢。
图 5 不同纳米粒子化学计量比下（a）MSD和（b）η随时间的变化曲线（约化温度T *= 0.90，纳米粒子体积分数φ = 0.10）
Figure 5. Effects of stoichiometric ratios on （a）mean square displacement （MSD） and（b）the nanoparticles’ fraction of largest aggregates （η） in terms of the assembly time（The reduced temperature is set as T *= 0.90, and the packing fraction φ = 0.10）
-
本文构建了DNA非均匀功能化纳米粒子的粗粒化模型，采用粗粒化分子动力学方法模拟了纳米粒子的可编程自组装过程，研究了化学计量比对自组装行为的影响。模拟发现：通过改变体系的化学计量比，DNA功能化纳米粒子从三维网络状超结构转变为支化超结构；构建了纳米粒子结构的几何模型，正确地预测了纳米粒子之间的相对位置及其分布；化学计量比显著地改变自组装动力学行为和自组装结构的动力学行为。
DNA非均匀功能化纳米粒子的可编程自组装行为
Programmable Self-Assembly of Non-Uniformly DNA-Functionalized Nanoparticles
-
摘要: DNA非均匀功能化纳米粒子作为一种可编程原子等价物，在多层次自组装结构领域具有重要的应用前景。构建了DNA非均匀功能化纳米粒子的粗粒化模型，并利用分子动力学模拟其自组装过程。通过计算机模拟发现，互补DNA序列间发生杂化反应，纳米粒子形成三维网络状和支化超结构；通过构建纳米粒子自组装结构的几何模型，能够正确预测纳米粒子之间的相对位置及其分布；通过调节互补的DNA功能化纳米粒子的化学计量比，显著地改变了自组装超结构和动力学行为。Abstract: DNA-functionalized nanoparticles, regarded as the programmable atom equivalent, enable the realization of hierarchically self-assembled superstructures. The self-assembled superstructures possessing unique mechanic, optical and electronic properties have prospective applications in the field of energy conservation, catalysis and medical diagnostics. With the development of nanotechnology, non-uniformly DNA-functionalized nanoparticles with DNA strands regioselectively distributed on the surfaces have been successfully synthesized. Recently, by utilizing the non-uniformly DNA-functionalized nanoparticles, researchers have created nanoparticle superstructures with complex architecture in the lab, such as discrete planet-satellite nanostructures, one-dimensional nanoparticle chains, and even three-dimensional networks. However, the mechanism and design rules of self-assembly of non-uniformly DNA-functionalized nanoparticles remain to be explored. Herein, the coarse-grained model of non-uniformly DNA-functionalized nanoparticles is constructed, and molecular dynamics is utilized to simulate the self-assembly process of DNA-programmable nanoparticles. It is demonstrated that the non-uniformly DNA-functionalized nanoparticles self-assemble into branched or even network-like superstructures through the hybridization of complementary DNA strands. The geometrical model of self-assembled superstructures is proposed to predict the relative position and distribution of nanoparticles inside the superstructures. Sparked by the molecular polymerization, we further explored the effect of stoichiometric ratio on the self-assembly of nanoparticles. The stoichiometric ratio of nanoparticles has remarkable effects on both the architecture of superstructures and the kinetics of DNA-programmable self-assembly of nanoparticles. As the stoichiometric ratio increased from 1.0 to 5.7, the self-assembled superstructures switch from the spanning networks to discrete branched clusters.
-
图 1 （a）DNA非均匀功能化纳米粒子的粗粒化模型；（b）互补DNA链的黏性端杂化示意图
Figure 1. （a）Coarse-grained model of non-uniformly DNA-functionalized nanoparticles;（b）Hybridization between complementary DNA strands containing ns spacer beads and nl linker beads
图 2 化学计量比f对自组装超结构的影响（图中只显示通过互补DNA序列杂化所形成的最大纳米粒子超结构，红线表示由互补DNA序列杂化形成的纳米粒子间连接）
Figure 2. Effect of the stoichiometric ratio on the superstructure of nanoparticles（Note that only the largest superstructures are shown; Red lines represent the connectivity between various nanoparticles through the hybridization of complementary DNA strands）
图 3 （a）纳米粒子自组装结构的几何模型（实线为通过互补DNA序列杂化所形成的纳米粒子间连接，虚线连接的是次邻近纳米粒子）；（b～d）不同化学计量比（f）和中性部分长度（ns）下纳米粒子径向分布函数g（r）（约化温度T *= 0.90，纳米粒子体积分数φ = 0.10）
Figure 3. （a）Geometrical representation of self-assembled superstructures of nanoparticles（Solid lines denote the connectivity between nanoparticles, and dashed lines denote the connectivity of nanoparticles with sub-nearest distance）；（b—d）Radial distribution function g（r）of nanoparticles for a range of stoichiometric ratios （f） and the lengths of spacer（ns）（The reduced temperature is set as T *= 0.90, and the packing fraction is set as φ = 0.10）
图 4 不同化学计量比下聚集体数量（Na）与DNA杂化碱基对数量（Nh）随时间的变化曲线（约化温度T *= 0.90，纳米粒子体积分数φ = 0.10）
Figure 4. Temporal evolution of numbers of aggregates （ Na） and numbers of hybridizations （Nh） at various stoichiometric ratios （The reduced temperature is set as T *= 0.90, and the packing fraction φ = 0.10）
图 5 不同纳米粒子化学计量比下（a）MSD和（b）η随时间的变化曲线（约化温度T *= 0.90，纳米粒子体积分数φ = 0.10）
Figure 5. Effects of stoichiometric ratios on （a）mean square displacement （MSD） and（b）the nanoparticles’ fraction of largest aggregates （η） in terms of the assembly time（The reduced temperature is set as T *= 0.90, and the packing fraction φ = 0.10）
-
[1] LI T I N G, SKNEPNEK R, OLVERA de La Cruz M. Thermally active hybridization drives the crystallization of DNA-functionalized nanoparticles [J]. Journal of the American Chemical Society,2013,135(23):8535-8541. doi: 10.1021/ja312644h [2] KIM A J, BIANCANIELLO P L, CROCKER J C. Engineering DNA-mediated colloidal crystallization [J]. Langmuir,2006,22(5):1991-2001. doi: 10.1021/la0528955 [3] MIRKIN C A, LETSINGER R L, MUCIC R C, et al. A DNA-based method for rationally assembling nanoparticles into macroscopic materials [J]. Nature,1996,382(6592):607-609. doi: 10.1038/382607a0 [4] ZHANG Y, LU F, YAGER K G, et al. A general strategy for the DNA-mediated self-assembly of functional nanoparticles into heterogeneous systems [J]. Nature Nanotechnology,2013,8(11):865-872. doi: 10.1038/nnano.2013.209 [5] GANG O, TKACHENKO A V. DNA-programmable particle superlattices: Assembly, phases, and dynamic control [J]. Mrs Bulletin,2016,41(5):381-387. doi: 10.1557/mrs.2016.92 [6] MACFARLANE R J, LEE B, JONES M R, et al. Nanoparticle superlattice engineering with DNA [J]. Science,2011,334(6053):204-208. doi: 10.1126/science.1210493 [7] MACFARLANE R J, THANER R V, BROWN K A, et al. Importance of the DNA “bond” in programmable nanoparticle crystallization [J]. Proceedings of the National Academy of Sciences,2014,111(42):14995-15000. doi: 10.1073/pnas.1416489111 [8] WANG Y, WANG Y, BREED D R, et al. Colloids with valence and specific directional bonding [J]. Nature,2012,491(7422):51-55. doi: 10.1038/nature11564 [9] MCMILLAN J R, MIRKIN C A. DNA-functionalized, bivalent proteins [J]. Journal of the American Chemical Society,2018,140(22):6776-6779. doi: 10.1021/jacs.8b03403 [10] ZHU G, XU Z, YANG Y. Hierarchical crystals formed from DNA-functionalized janus nanoparticles [J]. ACS Nano,2018,12(9):9467-9475. doi: 10.1021/acsnano.8b04753 [11] KASHIWAGI D, SIM S, NIWA T. Protein nanotube selectively cleavable with DNA: Supramolecular polymerization of “DNA-appended molecular chaperones” [J]. Journal of the American Chemical Society,2017,140(1):26-29. [12] KNOROWSKI C, BURLEIGH S, TRAVESSET A. Dynamics and statics of DNA-programmable nanoparticle self-assembly and crystallization [J]. Physical Review Letters,2011,106(21):215501. doi: 10.1103/PhysRevLett.106.215501 [13] KNOROWSKI C, TRAVESSET A. Materials design by DNA programmed self-assembly [J]. Current Opinion in Solid State and Materials Science,2011,15(6):262-270. doi: 10.1016/j.cossms.2011.07.002 [14] YU Q, HU J, HU Y, et al. Significance of DNA bond strength in programmable nanoparticle thermodynamics and dynamics [J]. Soft Matter,2018,14(14):2665-2670. doi: 10.1039/C7SM02456H [15] LYTTON-JEAN A K, MIRKIN C A. A thermodynamic investigation into the binding properties of DNA functionalized gold nanoparticle probes and molecular fluorophore probes [J]. Journal of the American Chemical Society,2005,127(37):12754-12755. doi: 10.1021/ja052255o [16] KNOROWSKI C, TRAVESSET A. Dynamics of DNA-programmable nanoparticle crystallization: Gelation, nucleation and topological defects [J]. Soft Matter,2012,8(48):12053-12059. doi: 10.1039/c2sm26832a [17] FLORY P J. Principles of polymer chemistry[M]. Ithaca, New York: Cornell University Press, 1953. [18] DAI W, HSU C W, SCIORTINO F, et al. Valency dependence of polymorphism and polyamorphism in DNA-functionalized nanoparticles [J]. Langmuir,2009,26(5):3601-3608. [19] KNOROWSKI C, TRAVESSET A. Self-assembly and crystallization of hairy (f-star) and DNA-grafted nanocubes [J]. Journal of the American Chemical Society,2014,136(2):653-659. doi: 10.1021/ja406241n [20] ZHANG R, SHI T, AN L, et al. Conformational study on sol-gel transition in telechelic polyelectrolytes solutions [J]. The Journal of Physical Chemistry B,2010,114(10):3449-3456. doi: 10.1021/jp9092404 [21] LI Y, SUN Z, SHI T, et al. Conformation studies on sol-gel transition in triblock copolymer solutions [J]. The Journal of Chemical Physics,2004,121(2):1133-1140. doi: 10.1063/1.1758938 [22] BHATTACHARYA S, KIEFFER J. Molecular dynamics simulation study of growth regimes during polycondensation of silicic acid: From silica nanoparticles to porous gels [J]. The Journal of Physical Chemistry C,2008,112(6):1764-1771. doi: 10.1021/jp073808f [23] WANG K W, BETANCOURT T, HALL C K. Computational study of DNA-cross-linked hydrogel formation for drug delivery applications [J]. Macromolecules,2018,51(23):9758-9768. doi: 10.1021/acs.macromol.8b01505 -
点击查看大图
计量
- 文章访问数: 5184
- HTML全文浏览量: 811
- PDF下载量: 27
- 被引次数: 0
