MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Szeged, Hungary.
J Comput Chem. 2022 Jan 5;43(1):19-28. doi: 10.1002/jcc.26767. Epub 2021 Oct 22.
We determine the proton affinity (PA) and gas-phase basicity (GB) of amino acid α-alanine at a chemically accurate level by performing explicitly-correlated CCSD(T)-F12b/aug-cc-pVDZ geometry optimizations and normal mode vibrational frequency calculations as well as CCSD(T)-F12b/aug-cc-pVTZ energy computations at the possible neutral and protonated geometries. Temperature effects at 298.15 K considering translational, rotational, and vibrational enthalpy and entropy corrections are obtained via standard statistical mechanics utilizing the molecular geometries and the harmonic vibrational energy levels. Both the amino nitrogen (N) and the carbonyl oxygen (O) atoms are proven to be potential protonation sites and a systematic conformational search reveals 3 N- and 9 O-protonated conformers in the 0.00-7.88 and 25.43-30.43 kcal/mol energy ranges at 0 K, respectively. The final computed PA and GB values at (0)298.15 K in case of N-protonation are (214.47)216.80 and 207.07 kcal/mol, respectively, whereas the corresponding values for O-protonation are (189.04)190.63 and 182.31 kcal/mol. The results of the benchmark high-level coupled-cluster computations are utilized to assess the accuracy of several lower-level cost-effective methods such as MP2 and density functional theory with various functionals (SOGGA11-X, M06-2X, PBE0, B3LYP, M06, TPSS).
我们通过执行明确相关的 CCSD(T)-F12b/aug-cc-pVDZ 几何优化和正则振动频率计算,以及 CCSD(T)-F12b/aug-cc-pVTZ 能量计算,在可能的中性和质子化几何结构上,以化学精确的水平确定氨基酸α-丙氨酸的质子亲和能 (PA) 和气相碱度 (GB)。在 298.15 K 下,考虑到平移、旋转和振动焓和熵校正的温度效应,通过标准统计力学利用分子几何形状和简谐振动能级来获得。氨基氮 (N) 和羰基氧 (O) 原子都被证明是潜在的质子化位点,系统的构象搜索在 0 K 时分别在 0.00-7.88 和 25.43-30.43 kcal/mol 能量范围内发现 3 个 N-和 9 个 O-质子化构象。在 N-质子化的情况下,最终计算得到的 298.15 K 时的 PA 和 GB 值分别为 (214.47)216.80 和 207.07 kcal/mol,而 O-质子化的相应值分别为 (189.04)190.63 和 182.31 kcal/mol。基准高精度耦合簇计算的结果被用来评估几种低成本效益方法的准确性,如 MP2 和各种函数 (SOGGA11-X、M06-2X、PBE0、B3LYP、M06、TPSS) 的密度泛函理论。