The DLPNO-CCSD(T) method is designed to study large molecular systems at significantly reduced cost relative to its canonical counterpart. However, the error in this approach is also size-extensive and relies on cancellation of errors for the calculation of relative energies. This work provides a direct comparison of canonical CCSD(T) and TightPNO DLPNO-CCSD(T) calculations of reaction energies and barriers of a broad range of chemical reactions. The dataset includes acidities, anion binding affinities, enolization, Diels−Alder, nucleophilic substitution, and atom transfer reactions and complements existing theoretical datasets in terms of system size as well as new reaction types (e.g., anion binding affinities and chlorine atom transfer reactions). The performance of DLPNO-CCSD(T) was further examined with respect to systematic variation of basis set and system size and amounts of nonbonded interaction present in the system. The errors in the DLPNO-CCSD(T) were found to be relatively insensitive to the choice of basis set for small systems but increase monotonically with system size. Additionally, calculations of barriers appear to be more challenging than reaction energies with errors exceeding 5 kJ mol −1 for many Diels−Alder reactions. Further tests on three realistic organic reactions reveal the impact of the DLPNO approximation in calculating absolute and relative barriers that are important for predictions such as stereoselectivity.

Solvents are one of the key variables in the optimization of a synthesis yield or properties of a synthesis product. In this paper, contemporary solvent models are applied to predict the rates of SN2 reactions in a range of aqueous and non-aqueous solvents. High-level CCSD(T)/CBS//M06-2X/6-31+G(d) gas phase energies were combined with solvation free energies from SMD, SM12, and ADF-COSMO-RS continuum solvent models, as well as molecular mechanics (MM) explicit solvent models with different atomic charge schemes to predict the rate constants of three SN2 reactions in eight protic and aprotic solvents. It is revealed that the prediction of rate constants in organic solvents is not necessarily less challenging than in water and popular solvent models struggle to predict their rate constants to within 3 log units of experimental values. Among the continuum solvent models, the ADF-COSMO-RS model performed the best in predicting absolute rate contants while the SM12 model was best at predicting relative rate constants with an average accuracy of about 1.5 and 0.8 log units, respectively. The use of computationally more demanding MM explicit solvent models did not translate to improvements in absolute rate constants but was quite effective at predicting relative rate constants due to systematic error cancellation. Free energy barriers obtained from umbrella sampling with explicit solvent QM/MM simulations led to excellent agreement with experimental values, provided that a validated level of theory is used to treat the QM region.

Carbonic anhydrase is an attractive drug target for the treatment of many diseases. This paper examines the ability of end-state MM/GBSA methods to rank inhibitors of carbonic anhydrase in terms of their binding affinities. The MM/GBSA binding energies were evaluated using different atomic charge schemes (Mulliken, ESP and NPA) at different levels of theories, including Hartree–Fock, B3LYP-D3(BJ), and M06-2X with the 6–31G(d,p) basis set. For a large test set of 32 diverse inhibitors, the use of B3LYP-D3(BJ) ESP atomic charges yielded the strongest correlation with experiment (R2 = 0.77). The use of the recently enhanced Autodock Vina and zinc optimised AD4Zn force field also predicted ligand binding affinities with moderately strong correlation (R2 = 0.64) at significantly lower computational cost. However, the docked poses deviate significantly from crystal structures. Overall, this study demonstrates the applicability of docking to estimate ligand binding affinities for a diverse range of CA inhibitors, and indicates that more theoretically robust MM/GBSA simulations show promise for improving the accuracy of predicted binding affinities, as long as a validated set of parameters is used. Graphical abstract

Calculation of molecular geometries and harmonic vibrational frequencies are pre-requisites for thermochemistry calculations. Contrary to conventional wisdom, this paper demonstrates that quantum chemical predictions of the thermochemistry of many gas and solution phase chemical reactions appear to be very insensitive to the choice of basis sets. For a large test set of 80 diverse organic and transition-metal-containing reactions, variations in reaction free energy based on geometries and frequencies calculated using a variety of double and triple-zeta basis sets from the Pople, Jensen, Ahlrichs, and Dunning families are typically less than 4 kJ mol–1, especially when the quasiharmonic oscillator correction is applied to mitigate the effects of low-frequency modes. Our analysis indicates that for many organic molecules and their transition states, high-level revDSD-PBEP86-D4 and DLPNO-CCSD(T)/(aug-)cc-pVTZ single-point energies usually vary by less than 2 kJ mol–1 on density functional theory geometries optimized using basis sets ranging from 6-31+G(d) to aug-pcseg-2 and aug-cc-pVTZ. In cases where these single-point energies vary significantly, indicating sensitivity of molecular geometries to the choice of basis set, there is often substantial cancellation of errors when the reaction energy or barrier is calculated. The study concludes that the choice of basis set for molecular geometry and frequencies, particularly those considered in this study, is not critical for the accuracy of thermochemistry calculations in the gas or solution phase.

Solvents are one of the key variables in the optimisation of a synthesis yield or properties of a synthesis product. In this paper, contemporary solvent models are applied to predict the rates of SN2 reactions in a range of aqueous and non-aqueous solvents. High-level CCSD(T)/CBS//M06-2X/6-31+G(d) gas phase energies were combined with solvation free energies from SMD, SM12 and ADF-COSMO-RS continuum solvent models as well as molecular mechanics (MM) explicit solvent models with different atomic charge schemes to predict the rate constants of three SN2 reactions in eight protic and aprotic solvents. It is revealed that popular solvent models struggle to predict their rate constants to within 3 log units of experimental values and deviations as large as 7.6 log units were observed. Amongst the implicit solvent models, the ADF-COSMO-RS model performed the best in predicting absolute rate constants with an average accuracy of 1.5 log units while the SM12 and CGenFF/TIP3P MM explicit solvent models were most accurate in the prediction of relative rate constants in different solvents due to systematic error cancellation. Free energy barriers obtained from umbrella sampling with explicit solvent QM/MM simulations led to excellent agreement with experiment provided that a validated level of theory is used to treat the QM region.