On the origin of the substantial stabilisation of the electron-donor 1,3-dithiole-2-thione-4-carboxyclic acid···I2 and DABCO···I2 complexes

The stabilisation energies of the crystal structures of 1,3-dithiole-2-thione-4-carboxyclic acid···I2 and DABCO···I2 complexes determined by the CCSD(T)/CBS method are very large and exceed 8 and 15 kcal mol(-1), respectively. The DFT-D method (B97-D3/def2-QZVP) strongly overestimates these stabilisation energies, which support the well-known fact that the DFT-D method is not very applicable to the study of charge-transfer complexes. On the other hand, the M06-2X/def2-QZVP method provides surprisingly reliable energies. A DFT-SAPT analysis has shown that a substantial stabilisation of these complexes arises from the charge-transfer energy included in the induction energy and that the respective induction energy is much larger than that of other non-covalently bound complexes. The total stabilisation energies of the complexes mentioned as well as of those where iodine has been replaced by lighter halogens (Br2 and Cl2) or by hetero systems (IF, ICH3, N2) correlate well with the magnitude of the σ-hole (Vs,max value) as well as with the LUMO energy. The nature of the stabilisation of all complexes between both electron donors and X2 (X = I, Br, Cl, N) systems is explained by the magnitude of the σ-hole but surprisingly also by the values of the electric quadrupole moment of these systems. Evidently, the nature of the stabilisation of halogen-bonded complexes between electron donors and systems where the first non-zero electric multipole moment is the quadrupole moment can be explained not only by the recently introduced concept of the σ-hole but also by the classical concept of electric quadrupole moments.

Selective induced polarization through electron transfer in acetone and pyrazole ester derivatives via C–H⋯OC interaction

A set of organic compounds (pyrazole ester derivatives,viz.5-[3-(substituted)-propoxy]-3-methyl-1-phenyl-1H-pyrazole-4-carboxylic acid methyl ester and 5-[2-(substituted)-ethoxy]-3-methyl-1-phenyl-1H-pyrazole-4-carboxylic acid methyl ester) was synthesized and their affinity and stability towards the acetone molecule were tested by NMR.

The nature of bonding and electronic properties of graphene and benzene with iridium adatoms

The strong correlation effects contribute to the opening of the band gap in graphene covered with the Ir adatoms.

Evaluation of composite schemes for CCSDT(Q) calculations of interaction energies of noncovalent complexes

Evaluation of composite schemes for CCSDT(Q) calculations of interaction energies of noncovalent complexes.

Statistical analysis of σ-holes: a novel complementary view on halogen bonding

Halogen σ-hole characteristics are defined and their typical values are provided based on an extensive analysis of more than 2500 molecules from the ZINC database.

Quantum mechanics-based scoring rationalizes the irreversible inactivation of parasitic Schistosoma mansoni cysteine peptidase by vinyl sulfone inhibitors

The quantum mechanics (QM)-based scoring function that we previously developed for the description of noncovalent binding in protein-ligand complexes has been modified and extended to treat covalent binding of inhibitory ligands. The enhancements are (i) the description of the covalent bond breakage and formation using hybrid QM/semiempirical QM (QM/SQM) restrained optimizations and (ii) the addition of the new ΔG(cov)' term to the noncovalent score, describing the "free" energy difference between the covalent and noncovalent complexes. This enhanced QM-based scoring function is applied to a series of 20 vinyl sulfone-based inhibitory compounds inactivating the cysteine peptidase cathepsin B1 of the Schistosoma mansoni parasite (SmCB1). The available X-ray structure of the SmCB1 in complex with a potent vinyl sulfone inhibitor K11017 is used as a template to build the other covalently bound complexes and to model the derived noncovalent complexes. We present the correlation of the covalent score and its constituents with the experimental binding data. Four outliers are identified. They contain bulky R1' substituents structurally divergent from the template, which might induce larger protein rearrangements than could be accurately modeled. In summary, we propose a new computational approach and an optimal protocol for the rapid evaluation and prospective design of covalent inhibitors with a conserved binding mode.

Modulation of aldose reductase inhibition by halogen bond tuning

In this paper, we studied a designed series of aldose reductase (AR) inhibitors. The series was derived from a known AR binder, which had previously been shown to form a halogen bond between its bromine atom and the oxygen atom of the Thr-113 side chain of AR. In the series, the strength of the halogen bond was modulated by two factors, namely bromine-iodine substitution and the fluorination of the aromatic ring in several positions. The role of the single halogen bond in AR-ligand binding was elucidated by advanced binding free energy calculations involving the semiempirical quantum chemical Hamiltonian. The results were complemented with ultrahigh-resolution X-ray crystallography and IC50 measurements. All of the AR inhibitors studied were shown by X-ray crystallography to bind in an identical manner. Further, it was demonstrated that it was possible to decrease the IC50 value by about 1 order of magnitude by tuning the strength of the halogen bond by a monoatomic substitution. The calculations revealed that the protein-ligand interaction energy increased upon the substitution of iodine for bromine or upon the addition of electron-withdrawing fluorine atoms to the ring. However, the effect on the binding affinity was found to be more complex due to the change of the solvation/desolvation properties within the ligand series. The study shows that it is possible to modulate the strength of a halogen bond in a protein-ligand complex as was designed based on the previous studies of low-molecular-weight complexes.

The accuracy of quantum chemical methods for large noncovalent complexes

We evaluate the performance of the most widely used wavefunction, density functional theory, and semiempirical methods for the description of noncovalent interactions in a set of larger, mostly dispersion-stabilized noncovalent complexes (the L7 data set). The methods tested include MP2, MP3, SCS-MP2, SCS(MI)-MP2, MP2.5, MP2.X, MP2C, DFT-D, DFT-D3 (B3-LYP-D3, B-LYP-D3, TPSS-D3, PW6B95-D3, M06-2X-D3) and M06-2X, and semiempirical methods augmented with dispersion and hydrogen bonding corrections: SCC-DFTB-D, PM6-D, PM6-DH2 and PM6-D3H4. The test complexes are the octadecane dimer, the guanine trimer, the circumcoronene…adenine dimer, the coronene dimer, the guanine-cytosine dimer, the circumcoronene…guanine-cytosine dimer, and an amyloid fragment trimer containing phenylalanine residues. The best performing method is MP2.5 with relative root mean square deviation (rRMSD) of 4 %. It can thus be recommended as an alternative to the CCSD(T)/CBS (alternatively QCISD(T)/CBS) benchmark for molecular systems which exceed current computational capacity. The second best non-DFT method is MP2C with rRMSD of 8 %. A method with the most favorable "accuracy/cost" ratio belongs to the DFT family: BLYP-D3, with an rRMSD of 8 %. Semiempirical methods deliver less accurate results (the rRMSD exceeds 25 %). Nevertheless, their absolute errors are close to some much more expensive methods such as M06-2X, MP2 or SCS(MI)-MP2, and thus their price/performance ratio is excellent.

Assessing the accuracy and performance of implicit solvent models for drug molecules: conformational ensemble approaches

The accuracy and performance of implicit solvent methods for solvation free energy calculations were assessed on a set of 20 neutral drug molecules. Molecular dynamics (MD) provided ensembles of conformations in water and water-saturated octanol. The solvation free energies were calculated by popular implicit solvent models based on quantum mechanical (QM) electronic densities (COSMO-RS, MST, SMD) as well as on molecular mechanical (MM) point-charge models (GB, PB). The performance of the implicit models was tested by a comparison with experimental water-octanol transfer free energies (ΔG(ow)) by using single- and multiconformation approaches. MD simulations revealed difficulties in a priori estimation of the flexibility features of the solutes from simple structural descriptors, such as the number of rotatable bonds. An increasing accuracy of the calculated ΔG(ow) was observed in the following order: GB1 ~ PB < GB7 ≪ MST < SMD ~ COSMO-RS with a clear distinction identified between MM- and QM-based models, although for the set excluding three largest molecules, the differences among COSMO-RS, MST, and SMD were negligible. It was shown that the single-conformation approach applied to crystal geometries provides a rather accurate estimate of ΔG(ow) for rigid molecules yet fails completely for the flexible ones. The multiconformation approaches improved the performance, but only when the deformation contribution was ignored. It was revealed that for large-scale calculations on small molecules a recent GB model, GB7, provided a reasonable accuracy/speed ratio. In conclusion, the study contributes to the understanding of solvation free energy calculations for physical and medicinal chemistry applications.

On the importance and origin of aromatic interactions in chemistry and biodisciplines

Aromatic systems contain both σ- and π-electrons, which in turn constitute σ- and π-molecular orbitals (MOs). In discussing the properties of these systems, researchers typically refer to the highest occupied and lowest unoccupied MOs, which are π MOs. The characteristic properties of aromatic systems, such as their low ionization potentials and electron affinities, high polarizabilities and stabilities, and small band gaps (in spectroscopy called the N → V1 space), can easily be explained based on their electronic structure. These one-electron properties point to characteristic features of how aromatic systems interact with each other. Unlike hydrogen bonding systems, which primarily interact through electrostatic forces, complexes containing aromatic systems, especially aromatic stacked pairs, are predominantly stabilized by dispersion attraction. The stabilization energy in the benzene dimer is rather small (~2.5 kcal/mol) but strengthens with heteroatom substitution. The stacked interaction of aromatic nucleic acid bases is greater than 10 kcal/mol, and for the most stable stacked pair, guanine and cytosine, it reaches approximately 17 kcal/mol. Although these values do not equal the planar H-bonded interactions of these bases (~29 kcal/mol), stacking in DNA is more frequent than H-bonding and, unlike H-bonding, is not significantly weakened when passing from the gas phase to a water environment. Consequently, the stacking of aromatic systems represents the leading stabilization energy contribution in biomacromolecules and in related nanosystems. Therefore stacking (dispersion) interactions predominantly determine the double helical structure of DNA, which underlies its storage and transfer of genetic information. Similarly, dispersion is the dominant contributor to attractive interactions involving aromatic amino acids within the hydrophobic core of a protein, which is critical for folding. Therefore, understanding the nature of aromatic interactions, which depend greatly on quantum mechanical (QM) calculations, is of key importance in biomolecular science. This Account shows that accurate binding energies for aromatic complexes should be based on computations made at the (estimated) CCSD(T)/complete basis set limit (CBS) level of theory. This method is the least computationally intensive one that can give accurate stabilization energies for all common classes of noncovalent interactions (aromatic-aromatic, H-bonding, ionic, halogen bonding, charge-transfer, etc.). These results allow for direct comparison of binding energies between different interaction types. Conclusions based on lower-level QM calculations should be considered with care.

Quantification of the interaction forces between metals and graphene by quantum chemical calculations and dynamic force measurements under ambient conditions

The two-dimensional material graphene has numerous potential applications in nano(opto)electronics, which inevitably involve metal graphene interfaces.Theoretical approaches have been employed to examine metal graphene interfaces, but experimental evidence is currently lacking. Here, we combine atomic force microscopy (AFM) based dynamic force measurements and density functional theory calculations to quantify the interaction between metal-coated AFM tips and graphene under ambient conditions. The results show that copper has the strongest affinity to graphene among the studied metals (Cu, Ag, Au, Pt, Si), which has important implications for the construction of a new generation of electronic devices. Observed differences in the nature of the metal-graphene bonding are well reproduced by the calculations, which included nonlocal Hartree-Fock exchange and van der Waals effects.