Non-covalent bonds are an important class of intermolecular interactions, which result in the ordering of atoms and molecules on the supramolecular scale. One such type of interaction is brought about by the bond formation between a region of positive electrostatic potential (σ-hole) interacts and a Lewis base. Previously, the halogen bond has been extensively studied as an example of a σ-hole interaction, where the halogen atom acts as the bond donor. Similarly, carbon, and the other tetrel elements can participate in σ-hole bonds. This thesis explores the nature of the carbon tetrel bond through the use of computational chemistry and solid state nuclear magnetic resonance (NMR) spectroscopy. The results of calculations of interaction energies and NMR parameters are reported for a series of model compounds exhibiting tetrel bonding from a methyl carbon to the oxygen and nitrogen atoms in a range of functional groups. The ¹³C chemical shift (𝛿iso) and the ¹ᶜ𝐽(¹³C,¹⁷O/¹⁵N) coupling across the tetrel bond are recorded as a function of geometry. The sensitivity of the NMR parameters to the non-covalent interaction is demonstrated via an increase in 𝛿iso and in |¹ᶜ𝐽(¹³C,¹⁷O/¹⁵N)| as the tetrel bond strengthens. There is no direct correlation between the NMR trends and the interaction energy curves; the energy minimum does not appear to correspond to a maximum or minimum chemical shift or J-coupling value. Gauge-including projector-augmented wave density functional theory (DFT) calculations of 𝛿iso are reported for crystals which exhibit tetrel bonding in the solid state. Experimental 𝛿iso values for sarcosine, betaine and caffeine and their tetrel-bonded salts generally corroborate the computational findings. This work offers new insights into tetrel bonding and facilitates the incorporation of tetrel bonds as restraints in NMR crystallographic structure refinement.

The aim of this thesis is to investigate strength and origin of the stabilization for various types of noncovalent interactions. As this knowledge could lead to a deeper understand- ing and rationalization of the binding phenomena. Further, to participate on the de- velopment of new noncovalent data sets, which are nowadays inevitable in the process of parametrization and validation of new computational methods. In all the studies, different binding motifs of model complexes, which represent usually crystal structures, structures from unrelaxed scans or the local minima, were investi- gated. The calculations of the reference stabilization energies were carried out at ab initio level (e.g. CCSD(T)/CBS, QCISD(T)/CBS). Further, the accuracy of more ap- proximate methods (e.g. MP2.5, DFT-D or SQM methods) toward reference method, was tested. In order to obtain the nature of the stabilization the DFT-SAPT decompo- sition was frequently utilized. In the first part of the thesis, the importance and basic characteristics of different types of noncovalent interactions (e.g. halogen bond, hydrogen bond, π· · · π interaction etc.), are discussed. The second part provides the description of computational methods which were essential for our investigation. The third part of the thesis provides an overview for part...

Noncovalent interactions are of great importance in studies on crystal design and drug discovery. One such noncovalent interaction, halogen bonding, is present between a covalently bound halogen atom and a Lewis base. A halogen bond is a directional interaction caused by the anisotropic distribution of charge on a halogen atom X covalently bound to A, which in turn forms a positive region called σ-hole on the A–X axis. Utilization of halogen bonds in lead optimization have been rarely considered in drug discovery until recently and yet more than 50% of the drug candidates are halogenated. To date, the halogen bond has not been subjected to practical molecular mechanical-molecular dynamics (MM-MD) study, where this noncovalent interaction cannot be described by conventional force fields because they do not account for the anisotropic distribution of the charge density on the halogen atoms. This problem was solved by the author and, for the first time, an extra-point of positive charge was used to represent the σ-hole on the halogen atom. This approach is called positive extra-point (PEP) approach. Interestingly, it was found that the performance of the PEP approach in describing halogen bond was better than the semiempirical methods including the recent halogen-bond corrected PM6 (PM6-DH2X) method. The PEP approach also gave promising results in describing other noncovalent halogen interactions, such as C–X···H and C–X···π-systems. The PEP resulted in an improvement in the accuracy of the electrostatic-potential derived charges of halogen-containing molecules, giving in turn better dipole moments and solvation free energies compared to high-level quantum mechanical and experimental data.With the aid of our PEP approach, the first MM-molecular dynamics (MM-MD) study of inhibitors that form a halogen bond with a receptor was performed for tetrahalobenzotriazole inhibitors complexed to cyclin-dependent protein kinase (CDK2). When the PEP approach was used, the calculated MM-generalized Born surface area (MM-GBSA)//MM-MD binding energies for halobenzimidazole and halobenzotriazole inhibitors complexed with protein kinase CK2 were found to correlate well with the corresponding experimental data, with correlation coefficients R2 of greater than 0.90. The nature and strength of halogen bonding in halo molecule···Lewis base complexes were studied in terms of molecular mechanics using our PEP approach. The contributions of the σ-hole (i.e., positively charged extra-point) and the halogen atom to the strength of this noncovalent interaction were clarified using the atomic parameter contribution to the molecular interaction approach. The molecular mechanical results revealed that the halogen bond is electrostatic and van der Waals in nature. The strength of the halogen bond increases with increasing the magnitude of the extra-point charge. The van der Waals interaction’s contribution to the halogen bond strength is most favorable in chloro complexes, whereas the electrostatic interaction is dominant in iodo complexes.The failure of the PM6 semiempirical method in describing noncovalent halogen interactions —not only halogen bonds, but also hydrogen bonds involving halogen atoms— was reported and corrected by the introduction of a second and third generation of noncovalent halogen interactions correction. The developed correction yielded promising results for the four examined noncovalent halogen interactions, namely: C–X···O, C–X···N, C–X···π-system, and C–X···H interactions.