Assessment of various density functionals for intermolecular N→Sn interactions: The test case of trimethyltin cyanide dimer

Theoretical insight into the interaction between SnX2 (X = H, F, Cl, Br, I) and benzene

For a series of five model complexes composed of a singlet SnX₂ molecule (X = H, F, Cl, Br, I) and a benzene molecule, the first-principles calculations of their energetics and the analysis of their electron density topology have been performed. The CCSD(T)/CBS interaction energy between SnX₂ and C₆H₆ fall into the range between −10.0 and −11.2 kcal/mol, which indicates that the complexes are rather weakly bound. The relevant role of electrostatic and dispersion contributions to the interaction energy between SnX₂ and C₆H₆ is highlighted in the results obtained from the symmetry-adapted perturbation theory (SAPT). The electron density topological analysis has been carried out using the quantum theory of atoms in molecules (QTAIM) and the noncovalent interactions (NCI) visualization index. Both QTAIM and NCI prove the closed-shell, noncovalent and attractive character of the interaction. A very small charge transfer from C₆H₆ to SnX₂ has been detected. The formation of the five complexes is accompanied by the electron density deformations that are spatially restricted mostly to the region around the Sn atom and its adjacent C atom. The results presented in this work shed some light on the nature of the interactions associated with crystalline structural motifs involving low-valent tin complexed with neutral aryl rings.

Binding energies of tyrosine kinase inhibitors: Error assessment of computational methods for imatinib and nilotinib binding

The binding energies of imatinib and nilotinib to tyrosine kinase have been determined by quantum mechanical (QM) computations, and compared with literature binding energy studies using molecular mechanics (MM). The potential errors in the computational methods include these critical factors. Errors in X-ray structures such as structural distortions and steric clashes give unrealistically high van der Waals energies, and erroneous binding energies.MM optimization gives a very different configuration to the QM optimization for nilotinib, whereas the imatinib ion gives similar configurations. Solvation energies are a major component of the overall binding energy. The QM based solvent model (PCM/SMD) gives different values from those used in the implicit PBSA solvent MM models. A major error in inhibitor—kinase binding lies in the non-polar solvation terms. Solvent transfer free energies and the required empirical solvent accessible surface area factors for nilotinib and imatinib ion to give the transfer free energies have been reverse calculated. These values differ from those used in the MM PBSA studies.An intertwined desolvation—conformational binding selectivity process is a balance of thermodynamic desolvation and intramolecular conformational kinetic control.The configurational entropies (TΔS) are minor error sources.

Performance of Møller-Plesset second-order perturbation theory and density functional theory in predicting the interaction between stannylenes and aromatic molecules

The performances of Møller-Plesset second-order perturbation theory (MP2) and density functional theory (DFT) have been assessed for the purposes of investigating the interaction between stannylenes and aromatic molecules. The complexes between SnX2 (where X = H, F, Cl, Br, and I) and benzene or pyridine are considered. Structural and energetic properties of such complexes are calculated using six MP2-type and 14 DFT methods. The assessment of the above-mentioned methods is based on the comparison of the structures and interaction energies predicted by these methods with reference computational data. A very detailed analysis of the performances of the MP2-type and DFT methods is carried out for two complexes, namely SnH2-benzene and SnH2-pyridine. Of the MP2-type methods, the reference structure of SnH2-benzene is reproduced best by SOS-MP2, whereas SCS-MP2 is capable of mimicking the reference structure of SnH2-pyridine with the greatest accuracy. The latter method performs best in predicting the interaction energy between SnH2 and benzene or pyridine. Among the DFT methods, ωB97X provides the structures and interaction energies of the SnH2-benzene and SnH2-pyridine complexes with good accuracy. However, this density functional is not as effective in reproducing the reference data for the two complexes as the best performing MP2-type methods. Next, the DFT methods are evaluated using the full test set of SnX2-benzene and SnX2-pyridine complexes. It is found that the range-separated hybrid or dispersion-corrected density functionals should be used for describing the interaction in such complexes with reasonable accuracy.