Bio-X Lab, Department of Physics, and Soft Matter Research Center, Zhejiang University, Hangzhou 310027, China.
J Chem Phys. 2012 Jan 14;136(2):025103. doi: 10.1063/1.3675486.
Understanding the interaction between carbon nanotubes (CNTs) and biomolecules is essential to the CNT-based nanotechnology and biotechnology. Some recent experiments have suggested that the π-π stacking interactions between protein's aromatic residues and CNTs might play a key role in their binding, which raises interest in large scale modeling of protein-CNT complexes and associated π-π interactions at atomic detail. However, there is concern on the accuracy of classical fixed-charge molecular force fields due to their classical treatments and lack of polarizability. Here, we study the binding of three aromatic residue analogues (mimicking phenylalanine, tyrosine, and tryptophan) and benzene to a single-walled CNT, and compare the molecular mechanical (MM) calculations using three popular fixed-charge force fields (OPLSAA, AMBER, and CHARMM), with quantum mechanical (QM) calculations using the density-functional tight-binding method with the inclusion of dispersion correction (DFTB-D). Two typical configurations commonly found in π-π interactions are used, one with the aromatic rings parallel to the CNT surface (flat), and the other perpendicular (edge). Our calculations reveal that compared to the QM results the MM approaches can appropriately reproduce the strength of π-π interactions for both configurations, and more importantly, the energy difference between them, indicating that the various contributions to π-π interactions have been implicitly included in the van der Waals parameters of the standard MM force fields. Meanwhile, these MM models are less accurate in predicting the exact structural binding patterns (matching surface), meaning there are still rooms to be improved. In addition, we have provided a comprehensive and reliable QM picture for the π-π interactions of aromatic molecules with CNTs in gas phase, which might be used as a benchmark for future force field developments.
理解碳纳米管(CNT)与生物分子之间的相互作用对于基于 CNT 的纳米技术和生物技术至关重要。一些最近的实验表明,蛋白质芳香族残基与 CNT 之间的π-π堆积相互作用可能在它们的结合中起着关键作用,这引起了人们对大规模建模蛋白质-CNT 复合物和相关 π-π相互作用的原子细节的兴趣。然而,由于经典处理和缺乏极化率,经典固定电荷分子力场的准确性受到关注。在这里,我们研究了三个芳香族残基类似物(模拟苯丙氨酸、酪氨酸和色氨酸)和苯与单壁 CNT 的结合,并比较了三种流行的固定电荷力场(OPLSAA、AMBER 和 CHARMM)的分子力学(MM)计算,以及使用包含色散校正的密度泛函紧束缚方法(DFTB-D)的量子力学(QM)计算。使用两种常见的π-π相互作用构型,一种是芳环与 CNT 表面平行(平面),另一种是垂直(边缘)。我们的计算表明,与 QM 结果相比,MM 方法可以适当再现两种构型的π-π相互作用的强度,更重要的是,它们之间的能量差异,表明各种对π-π相互作用的贡献已经隐含在标准 MM 力场的范德华参数中。同时,这些 MM 模型在预测精确的结构结合模式(匹配表面)方面不太准确,这意味着仍有改进的空间。此外,我们还为气相中芳香族分子与 CNT 之间的π-π相互作用提供了全面可靠的 QM 图像,这可能作为未来力场发展的基准。