Institut für Anorganische Chemie, Universität Göttingen, Göttingen, Germany.
Departamento de Química, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil.
J Comput Chem. 2023 Mar 5;44(6):745-754. doi: 10.1002/jcc.27037. Epub 2022 Nov 25.
Since atomic or functional-group properties in the bulk are generally not available from experimental methods, computational approaches based on partitioning schemes have emerged as a rapid yet accurate pathway to estimate the materials behavior from chemically meaningful building blocks. Among several applications, a comprehensive and systematically built database of atomic or group polarizabilities and related opto-electronic quantities would be very useful not only to envisage linear or non-linear optical properties of biomacromolecules but also to improve the accuracy of classical force fields devoted to simulate biochemical processes. In this work, we propose the first entries of such database that contains distributed polarizabilities and dipole moments extracted from fragments of peptides. Twenty three prototypical conformers of the dipeptides alanine-alanine and glycine-glycine were used to extract functional groups such as CH , CHCH , NH , COOH, CONH, thus allowing construction of a diversity of chemically relevant environments. To evaluate the accuracy of our database, reconstructed properties of larger peptides containing up to six residues of alanine and glycine were tested against density functional theory calculations at the M06-HF/aug-cc-pVDZ level of theory. The procedure is particularly accurate for the diagonal components of the polarizability tensor with errors up to 15%. In order to include solvent effects explicitly, the peptides were also surrounded by a box of water molecules whose distribution was optimized using the CHARMM force field. Solvent effects introduced by a classical dipole-dipole interaction model were compared to those obtained from polarizable-continuum model calculations.
由于从实验方法中通常无法获得体相的原子或官能团性质,因此基于划分方案的计算方法已成为从具有化学意义的构建块估计材料行为的快速而准确的途径。在几种应用中,原子或基团极化率和相关光电量的综合且系统构建的数据库不仅对于设想生物大分子的线性或非线性光学性质非常有用,而且对于提高用于模拟生化过程的经典力场的准确性也非常有用。在这项工作中,我们提出了此类数据库的第一个条目,其中包含从肽片段中提取的分布式极化率和偶极矩。使用二十三个丙氨酸-丙氨酸和甘氨酸-甘氨酸二肽的原型构象来提取官能团,例如CH ,CHCH ,NH, COOH,CONH,从而允许构建各种具有化学相关性的环境。为了评估我们的数据库的准确性,使用 M06-HF/aug-cc-pVDZ 理论水平的密度泛函理论计算测试了包含多达六个丙氨酸和甘氨酸残基的更大肽的重建性质。对于极化率张量的对角分量,该方法的准确性特别高,误差可达 15%。为了明确包括溶剂效应,肽也被水分子盒包围,其分布使用 CHARMM 力场进行了优化。经典偶极-偶极相互作用模型引入的溶剂效应与从极化连续体模型计算获得的溶剂效应进行了比较。