Role and nature of halogen bonding in inhibitor⋅⋅⋅receptor complexes for drug discovery: casein kinase-2 (CK2) inhibition as a case study

Halogenated Antimicrobial Agents to Combat Drug-Resistant Pathogens

Antimicrobial resistance presents us with a potential global crisis as it undermines the abilities of conventional antibiotics to combat pathogenic microbes. The history of antimicrobial agents is replete with examples of scaffolds containing halogens. In this review, we discuss the impacts of halogen atoms in various antibiotic types and antimicrobial scaffolds and their modes of action, structure-activity relationships, and the contributions of halogen atoms in antimicrobial activity and drug resistance. Other halogenated molecules, including carbohydrates, peptides, lipids, and polymeric complexes, are also reviewed, and the effects of halogenated scaffolds on pharmacokinetics, pharmacodynamics, and factors affecting antimicrobial and antivirulence activities are presented. Furthermore, the potential of halogenation to circumvent antimicrobial resistance and rejuvenate impotent antibiotics is addressed. This review provides an overview of the significance of halogenation, the abilities of halogens to interact in biomolecular settings and enhance pharmacological properties, and their potential therapeutic usages in preventing a postantibiotic era. SIGNIFICANCE STATEMENT: Antimicrobial resistance and the increasing impotence of antibiotics are critical threats to global health. The roles and importance of halogen atoms in antimicrobial drug scaffolds have been established, but comparatively little is known of their pharmacological impacts on drug resistance and antivirulence activities. This review is the first to extensively evaluate the roles of halogen atoms in various antibiotic classes and pharmacological scaffolds and to provide an overview of their ability to overcome antimicrobial resistance.

External Electric Field Effect on the Strength of σ-Hole Interactions: A Theoretical Perspective in Like⋯Like Carbon-Containing Complexes

For the first time, σ-hole interactions within like⋯like carbon-containing complexes were investigated, in both the absence and presence of the external electric field (EEF). The effects of the directionality and strength of the utilized EEF were thoroughly unveiled in the (F-C-F₃)₂, (F-C-H₃)₂, and (H-C-F₃)₂ complexes. In the absence of the EEF, favorable interaction energies, with negative values, are denoted for the (F-C-F₃)₂ and (H-C-F₃)₂ complexes, whereas the (F-C-H₃)₂ complex exhibits unfavorable interactions. Remarkably, the strength of the applied EEF exhibits a prominent role in turning the repulsive forces within the latter complex into attractive ones. The symmetrical nature of the considered like⋯like carbon-containing complexes eradicated the effect of directionality of the EEF. The quantum theory of atoms in molecules (QTAIM), and the noncovalent interaction (NCI) index, ensured the occurrence of the attractive forces, and also outlined the substantial contributions of the three coplanar atoms to the total strength of the studied complexes. Symmetry-adapted perturbation theory (SAPT) results show the dispersion-driven nature of the interactions.

σ-Hole and LP-Hole Interactions of Pnicogen···Pnicogen Homodimers under the External Electric Field Effect: A Quantum Mechanical Study

σ-Hole and lone-pair (lp)-hole interactions within σ-hole···σ-hole, σ-hole···lp-hole, and lp-hole···lp-hole configurations were comparatively investigated on the pnicogen···pnicogen homodimers (PCl3)2, for the first time, under field-free conditions and the influence of the external electric field (EEF). The electrostatic potential calculations emphasized the impressive versatility of the examined PCl3 monomers to participate in σ-hole and lp-hole pnicogen interactions. Crucially, the sizes of σ-hole and lp-hole were enlarged under the influence of the positively directed EEF and decreased in the case of reverse direction. Interestingly, the energetic quantities unveiled more favorability of the σ-hole···lp-hole configuration of the pnicogen···pnicogen homodimers, with significant negative interaction energies, than σ-hole···σ-hole and lp-hole···lp-hole configurations. Quantum theory of atoms in molecules and noncovalent interaction index analyses were adopted to elucidate the nature and origin of the considered interactions, ensuring their closed shell nature and the occurrence of attractive forces within the studied homodimers. Symmetry-adapted perturbation theory-based energy decomposition analysis alluded to the dispersion force as the main physical component beyond the occurrence of the examined interactions. The obtained findings would be considered as a fundamental underpinning for forthcoming studies pertinent to chemistry, materials science, and crystal engineering.

Type I-IV Halogen⋯Halogen Interactions: A Comparative Theoretical Study in Halobenzene⋯Halobenzene Homodimers

In the current study, unexplored type IV halogen⋯halogen interaction was thoroughly elucidated, for the first time, and compared to the well-established types I−III interactions by means of the second-order Møller−Plesset (MP2) method. For this aim, the halobenzene⋯halobenzene homodimers (where halogen = Cl, Br, and I) were designed into four different types, parodying the considered interactions. From the energetic perspective, the preference of scouted homodimers was ascribed to type II interactions (i.e., highest binding energy), whereas the lowest binding energies were discerned in type III interactions. Generally, binding energies of the studied interactions were observed to decline with the decrease in the σ-hole size in the order, C6H5I⋯IC6H5 > C6H5Br⋯BrC6H5 > C6H5Cl⋯ClC6H5 homodimers and the reverse was noticed in the case of type IV interactions. Such peculiar observations were relevant to the ample contributions of negative-belt⋯negative-belt interactions within the C6H5Cl⋯ClC6H5 homodimer. Further, type IV torsional trans → cis interconversion of C6H5X⋯XC6H5 homodimers was investigated to quantify the π⋯π contributions into the total binding energies. Evidently, the energetic features illustrated the amelioration of the considered homodimers (i.e., more negative binding energy) along the prolonged scope of torsional trans → cis interconversion. In turn, these findings outlined the efficiency of the cis configuration over the trans analog. Generally, symmetry-adapted perturbation theory-based energy decomposition analysis (SAPT-EDA) demonstrated the predominance of all the scouted homodimers by the dispersion forces. The obtained results would be beneficial for the omnipresent studies relevant to the applications of halogen bonds in the fields of materials science and crystal engineering.

R•-hole interactions of group IV-VII radical-containing molecules: A comparative study

For the first time, the potentiality of the sp²-hybridized group IV-VII radical (R^•)-containing molecules to participate in R^•-hole interactions was comparatively assessed using ^•SiF_3,^•POF₂, ^•SO₂F, and ^•ClO₃ models in the trigonal pyramidal geometry. In that spirit, a plethora of quantum mechanical calculations was performed at the MP2/aug-cc-pVTZ level of theory. According to the results, all the investigated R^•-containing molecules exhibited potent versatility to engage in R^•-hole … Lewis base interactions with significant negative binding energies for the NCH-based complexes. The strength of R^•-hole interactions was perceived to obey the ^•ClO₃ … > ^•SO₂F … > ^•POF₂ … > ^•SiF₃ … Lewis base order, outlining an inverse correlation between the binding energy and the atomic size of the R^•-hole donor. Benchmarking of the binding energy at the CCSD/CBS(T) computational level was executed for all the explored interactions and addressed an obvious similarity between the MP2 and CCSD energetic findings. QTAIM analysis critically unveiled the closed-shell nature of the explored R^•-hole interactions. SAPT-EDA proclaimed the reciprocal contributions of electrostatic and dispersion forces to the total binding energy. These observations demonstrate in better detail the nature of R^•-hole interactions, leading to a convincing amelioration for versatile fields relevant to materials science and drug design.

Effect of External Electric Field on Tetrel Bonding Interactions in (FTF3···FH) Complexes (T = C, Si, Ge, and Sn)

A quantum chemical study was accomplished on the σ-hole interactions of the barely explored group IV elements, for the first time, in the absence and presence of the positively and negatively directed external electric field (EEF). The analyses of molecular electrostatic potential addressed the occurrence of the σ-hole on all the inspected tetrel atoms, confirming their salient versatility to engage in σ-hole interactions. MP2 energetic findings disclosed the occurrence of favorable σ-hole interactions within the tetrel bonding complexes. The tetrel bonding interactions became stronger in the order of C < Si < Ge < Sn for F-T-F3···FH complexes with the largest interaction energy amounting to -19.43 kcal/mol for the optimized F-Sn-F3···FH complex under the influence of +0.020 au EEF. The interaction energy conspicuously evolved by boosting the magnitude of the positively directed EEF value and declining the negatively directed EEF one. The decomposition analysis for the interaction energies was also executed in terms of symmetry-adapted perturbation theory, illuminating the dominant electrostatic contribution to all the studied complexes' interactions except carbon-based interactions controlled by dispersion forces. The outcomes that emerged from the current work reported significantly how the direction and strength of the EEF affect the tetrel-bonding interactions, leading to further improvements in the forthcoming studies of supramolecular chemistry and materials science.

On the Potentiality of X-T-X3 Compounds (T = C, Si, and Ge, and X = F, Cl, and Br) as Tetrel- and Halogen-Bond Donors

The versatility of the X-T-X3 compounds (where T = C, Si, and Ge, and X = F, Cl, and Br) to participate in tetrel- and halogen-bonding interactions was settled out, at the MP2/aug-cc-pVTZ level of theory, within a series of configurations for (X-T-X3)2 homodimers. The electrostatic potential computations ensured the remarkable ability of the investigated X-T-X3 monomers to participate in σ-hole halogen and tetrel interactions. The energetic findings significantly unveil the favorability of the tetrel···tetrel directional configuration with considerable negative binding energies over tetrel···halogen, type III halogen···halogen, and type II halogen···halogen analogs. Quantum theory of atoms in molecules and noncovalent interaction analyses were accomplished to disclose the nature of the tetrel- and halogen-bonding interactions within designed configurations, giving good correlations between the total electron densities and binding energies. Further insight into the binding energy physical meanings was invoked through using symmetry-adapted perturbation theory-based energy decomposition analysis, featuring the dispersion term as the most prominent force beyond the examined interactions. The theoretical results were supported by versatile crystal structures which were characterized by the same type of interactions. Presumably, the obtained findings would be considered as a solid underpinning for future supramolecular chemistry, materials science, and crystal engineering studies, as well as a fundamental linchpin for a better understanding of the biological activities of chemicals.

Optimized Halogen Atomic Radii for PBSA Calculations Using Off-Center Point Charges

In force-field methods, the usage of off-center point charges, also called extra points (EPs), is a common strategy to tackle the anisotropy of the electrostatic potential of covalently bonded halogens (X), thus allowing the description of halogen bonds (XBs) at the molecular mechanics/molecular dynamics (MM/MD) level. Diverse EP implementations exist in the literature differing on the charge sets and/or the X-EP distances. Poisson-Boltzmann and surface area (PBSA) calculations can be used to obtain solvation free energies (ΔG_solv) of small molecules, often to compute binding free energies (ΔG_bind) at the MM-PBSA level. This method depends, among other parameters, on the empirical assignment of atomic radii (PB radii). Given the multiplicity of off-center point-charge models and the lack of specific PB radii for halogens compatible with such implementations, in this work, we assessed the performance of PBSA calculations for the estimation of ΔG_solv values in water (ΔG_hyd), also conducting an optimization of the halogen PB radii (Cl, Br, and I) for each EP model. We not only expand the usage of EP models in the scope of the general AMBER force field (GAFF) but also provide the first optimized halogen PB radii in the context of the CHARMM general force field (CGenFF), thus contributing to improving the description of halogenated compounds in PBSA calculations.

Cospatial σ-Hole and Lone Pair Interactions of Square-Pyramidal Pentavalent Halogen Compounds with π-Systems: A Quantum Mechanical Study

In the spirit of the mounting interest in noncovalent interactions, the present study was conducted to scrutinize a special type that simultaneously involved both σ-hole and lone pair (lp) interactions with aromatic π-systems. Square-pyramidal pentavalent halogen-containing molecules, including X-Cl-F4, F-Y-F4, and F-I-X4 compounds (where X = F, Cl, Br, and I and Y = Cl, Br, and I) were employed as σ-hole/lp donors. On the other hand, benzene (BZN) and hexafluorobenzene (HFB) were chosen as electron-rich and electron-deficient aromatic π-systems, respectively. The investigation relied upon a variety of quantum chemical calculations that complement each other. The results showed that (i) the binding energy of the X-Y-F4···BZN complexes increased (i.e., more negative) as the Y atom had a larger magnitude of σ-hole, contrary to the pattern of X-Y-F4···HFB complexes; (ii) the interaction energies of X-Y-F4···BZN complexes were dominated by both dispersion and electrostatic contributions, while dispersive interactions dominated X-Y-F4···HFB complexes; and (iii) the X4 atoms in F-I-X4···π-system complexes governed the interaction energy pattern: the larger the X4 atoms were, the greater the interaction energies were, for the same π-system. The results had illuminating facets in regard to the rarely addressed cases of the σ-hole/lp contradictory scene.

Comparison of ±σ-hole and ±R˙-hole interactions formed by tetrel-containing complexes: a computational study

For the first time, unconventional ±R˙-hole interactions were unveiled in tetrel-containing complexes. The nature and characteristics of ±R˙-hole interactions were explored relative to their ±σ-hole counterparts for ˙TF3⋯ and W-T-F3⋯B/R˙/A complexes (where T = C, Si, and Ge, W = H and F, B = Lewis bases, R˙ = free radicals, and A = Lewis acids). In an effort to thoroughly investigate such interactions, a plethora of quantum mechanical calculations, including molecular electrostatic potential (MEP), maximum positive electrostatic potential (V s,max), point-of-charge (PoC), interaction energy, symmetry adapted perturbation theory (SAPT), and reduced density gradient-noncovalent interaction (RDG-NCI) calculations, were applied. The most notable findings to emerge from this study are that (i) from the electrostatic perspective, the molecular stabilization energies of ˙TF3 and W-T-F3 monomers became more negative as the Lewis basicity increased, (ii) the most stable complexes were observed for the ones containing Lewis bases, forming -σ-hole and -R˙-hole interactions, and the interaction energies systematically increased in the order H-T-F3⋯B < ˙TF3⋯B < F-T-F3⋯B, (iii) contrariwise, the +σ-hole and +R˙-hole interactions with Lewis acids are more energetically favorable in the order F-T-F3⋯A < ˙TF3⋯A < H-T-F3⋯A, and (iv) generally, the dispersion force plays a key role in stabilizing the tetrel-containing complexes, jointly with the electrostatic and induction forces for the interactions with Lewis bases and acids, respectively. Concretely, the findings presented in this paper add to our understanding of the characteristics and nature of such intriguing interactions.

External electric field effects on the σ-hole and lone-pair hole interactions of group V elements: a comparative investigation

σ-hole and lone-pair (lp) hole interactions of trivalent pnicogen-bearing (ZF3) compounds were comparatively scrutinized, for the first time, under field-free and external electric field (EEF) conditions. Conspicuously, the sizes of the σ-hole and lp-hole were increased by applying an EEF along the positive direction, while the sizes of both holes decreased through the reverse EEF direction. The MP2 energetic calculations of ZF3⋯FH/NCH complexes revealed that σ-holes exhibited more impressive interaction energies compared to the lp-holes. Remarkably, the strengths of σ-hole and lp-hole interactions evolved with the increment of the positive value of the considered EEF; i.e., the interaction energy increased as the utilized EEF value increased. Unexpectedly, under field-free conditions, nitrogen-bearing complexes showed superior strength for their lp-hole interactions than phosphorus-bearing complexes. However, the reverse picture was exhibited for the interaction energies of nitrogen- and phosphorus-bearing complexes interacting within lp-holes by applying the high values of a positively directed EEF. These results significantly demonstrate the crucial influence of EEF on the strength of σ-hole and lp-hole interactions, which in turn leads to an omnipresent enhancement for variable fields, including biological simulations and material science.

Unconventional Type III Halogen···Halogen Interactions: A Quantum Mechanical Elucidation of σ-Hole···σ-Hole and Di-σ-Hole Interactions

Herein, two unconventional type III halogen···halogen interactions, namely, σ-hole···σ-hole and di-σ-hole interactions, were reported in a series of halogenated complexes. In type III, the A-halogen···halogen angles are typically equal to 180°, and the occurrence of σ-hole on halogen atoms is mandatory. Using diverse quantum mechanical calculations, it was demonstrated that the occurrence of such interactions with binding energies varied from -0.35 to -1.30 kcal/mol. Symmetry-adapted perturbation theory-based energy decomposition analysis (SAPT-EDA) revealed that type III interactions are dominated by dispersion forces, while electrostatic forces are unfavorable. Cambridge Structure Database (CSD) survey unveiled the experimental evidence for the manifestation of σ-hole···σ-hole interactions in crystal structures. This work might be deemed as a foundation for a vast number of forthcoming crystal engineering and materials science studies.

σ-Hole and Lone-Pair Hole Interactions in Chalcogen-Containing Complexes: A Comparative Study

The potentiality of sp3-hybridized chalcogen-containing molecules to participate in lone-pair (lp) hole interactions was reported for the first time. lp hole interactions were characterized and compared to σ-hole ones for OF2 and SF2 molecules as a case study. Various quantum mechanical calculations, including molecular electrostatic potential (MEP), maximum positive electrostatic potential (V s,max), point of charge (PoC), symmetry-adapted perturbation theory (SAPT), quantum theory of atoms in molecule (QTAIM), and reduced density gradient-noncovalent interaction (RDG-NCI) calculations, were carried out. The more significant findings to emerge from this study are the following: (i) the V s,max calculation was proved to be an unreliable method to determine the precise σ-hole and lp hole locations. (ii) The maximum positive electrostatic potential of the σ hole and lp hole was found to be at the F-Chal···PoC angle (θ) of 180° and at the centroid of XYlp plane, respectively. (iii) Lewis basicity has a significant effect on the strength of σ-hole and lp hole interactions. (iv) The studied molecules more favorably interact with Lewis bases via the σ hole compared to the lp hole, and (v) stabilization of the σ-hole and lp hole interactions stems from the electrostatic and dispersion forces, respectively.

Thermodynamic contribution of iodine atom to the binding of heterogeneously polyhalogenated benzotriazoles by the catalytic subunit of human protein kinase CK2

A series of novel benzotriazole derivatives containing iodine atom(s) were synthesized. The binding of these compounds to the catalytic subunit of human protein kinase CK2 was evaluated using differential scanning fluorimetry. The obtained thermodynamic data were compared with those determined previously for the brominated and chlorinated benzotriazole analogues to get a deeper insight into the thermodynamic contribution of iodine substitution to the free energy of ligand binding. We have shown that iodine atom(s) attached to the benzene ring of benzotriazole enhance(s) its binding by the target protein. This effect is the strongest when two iodine atoms are attached at positions peripheral to the triazole ring, which according to the structures deposited in protein data bank may be indicative for the formation of the halogen bond between iodine and carbonyl groups of residues located in the hinge region of the protein. Finally, quantitative structure-activity relationship analysis pointed the solute hydrophobicity as the main factor contributing to the binding affinity.

Tackling Halogenated Species with PBSA: Effect of Emulating the σ-Hole

To model halogen-bond phenomena using classical force fields, an extra point (EP) of charge is frequently introduced at a given distance from the halogen (X) to emulate the σ-hole. The resulting molecular dynamics (MD) trajectories can be used in subsequent molecular mechanics (MM) combined with Poisson-Boltzmann and surface area calculations (PBSA) to estimate protein-ligand binding free energies (Δ G_bind). While EP addition improves the MM/MD description of halogen-containing systems, its effect on the calculation of solvation free energies (Δ G_solv) using the PBSA approach is yet to be assessed. As the PBSA calculations depend, among other parameters, on the empirical assignment of radii (PB radii), a problematic issue arises, since standard halogen radii are smaller than the typical X···EP distances, thus placing the EP within the solvent dielectric. Herein, we took a common literature EP parametrization scheme, which uses X···EP = R_min and RESP charges in the context of GAFF, and performed a comprehensive study on the performance of PBSA (using three different setups) in the calculation of Δ G_solv values for 142 halogenated compounds (bearing Cl, Br, or I) for which the experimental values are known. By conducting an optimization (minimizing the error against experimental values), we provide a new optimized set of halogen PB radii, for each PBSA setup, that should be used in the context of the aforementioned scenario. A simultaneous optimization of PB radii and X···EP distances shows that a wide range of distance/radius pairs can be used without significant loss of accuracy, therefore laying the basis for expanding this halogen radii optimization strategy to other force fields and EP implementations. As ligand Δ G_solv estimation is an important term in the determination of protein-ligand Δ G_bind, this work is particularly relevant in the framework of structure-based virtual screening and related computer-aided drug design routines.

Halogen bonding in halocarbon-protein complexes and computational tools for rational drug design

Introduction: Halogens have a prominent role in drug design. Often used as a mean to improve ADME properties, they are also becoming a tool in protein-ligand recognition given their ability to form a non-covalent interaction, termed halogen bond, where halogens act as electrophilic species interacting with electron-rich partners. Rational drug design of halogen-bonding lead molecules requires an accurate description of halocarbon-protein complexes by computational tools though not all methods are able to tackle this non-covalent interaction. Areas covered: The authors present a review of computational methodologies that can be used to properly describe halogen bonds in the context of protein-ligand complexes, providing also insights on how these methods can be used in the context of computer-aided drug design. Expert opinion: Although in the last few years many computational tools, ranging from fast screening methods to the more expensive QM calculations, have been developed to tackle the halogen bonding phenomenon, they are not yet standard in the literature. This will eventually change as official software distributions are including support for halogen bonding in their methods. Tackling desolvation of halogenated species seems to be a good strategy to improve the accuracy of computational methods, that will be more commonly used prior to laboratory work in the future.

Comparative investigation of interactions of hydrogen, halogen and tetrel bond donors with electron-rich and electron-deficient π-systems

Recently, noncovalent interactions in complexes and crystals have attracted considerable interest. The current study was thus designed to gain a better understanding of three seminal types of noncovalent interactions, namely: hydrogen, halogen and tetrel interactions with π-systems. This study was performed on three models of Lewis acids: X₃–C–H, F₃–C–X and F–T–F₃ (where X = F, Cl, Br and I; and T = C, Si, Ge and Sn) and three π-systems as Lewis bases: benzene (BZN), 1,3,5-trifluorobenzene (TFB) and hexafluorobenzene (HFB). Quantum mechanical calculations, including geometrical optimization, molecular electrostatic potential (MEP), maximum positive electrostatic potential (V_s,max), Point-of-Charge (PoC), potential energy surface (PES), quantum theory of atoms in molecules (QTAIM) and noncovalent interaction (NCI) calculations, were carried out at the MP2/aug cc-pVDZ level of theory. The binding energies were additionally benchmarked at the CCSD(T)/CBS level. The results showed that: (i) the binding energies of the X₃–C–H⋯π-system complexes were unexpectedly inversely correlated with the V_s,max values on the hydrogen atom but directly correlated with the X atomic sizes; (ii) the binding energies for the F₃–C–X⋯π-system and F–T–F₃⋯π-system complexes were correlated with the σ-hole magnitudes of the X and T atoms, respectively; and (iii) for the F₃–C–F⋯π-system complexes, the binding energy was as strong as the π-system was electron-deficient, indicating the dominating nucleophilic character of the fluorine atom. NCI analysis showed that the unexpected trend of X₃–C–H⋯π-system binding energies could be attributed to additional attractive interactions between the X atoms in the X₃–C–H molecule and the carbon atoms of the π-system. Furthermore, the I₃–Sn–H molecule was employed as a case study of hydrogen, halogen and tetrel interactions with π-systems. It was found that hydrogen and halogen interactions of the I₃–Sn–H molecule correlated with the electron-richness of the π-system. In contrast, tetrel interactions correlated with the electron deficiency of the π-system.

Three seminal types of noncovalent interaction, namely: hydrogen, halogen and tetrel interactions with π-systems, were investigated using quantum mechanical calculations.

Improved Modeling of Halogenated Ligand-Protein Interactions Using the Drude Polarizable and CHARMM Additive Empirical Force Fields

Halogenated ligands can participate in nonbonding interactions with proteins via halogen bond (XB) or halogen-hydrogen bond donor (X-HBD) interactions. In the context of molecular dynamics (MD) simulations, the accuracy of the simulations depends strongly on the force field (FF) used. To ensure good reproduction of XB and X-HBD interactions with proteins, we optimized the previously developed additive CHARMM36/CHARMM General force field (CGenFF) and Drude polarizable force field by including atom pair-specific Lennard-Jones parameters for aromatic halogen-protein interactions. The optimization targeted quantum mechanical interaction energy surfaces with the developed parameters then examined for their ability to reproduce experimental halogen-containing ligand-protein interactions in MD simulations. The calculated halogenated ligand interaction geometries were in good overall agreement with the experimental crystal data for both the polarizable and additive FFs, showing that these models can accurately treat both XB and X-HBD interactions. Analysis of the ligand-protein interactions shows significant contributions of polarizability to binding occurring in the Drude FF, with self-polarization energy making both favorable and unfavorable contributions to binding. Further analysis of the dipole moments from aqueous solution to protein indicates the polarizable FF accounts for subtle changes of the environment of the ligands that can impact binding. The present work demonstrates the utility of the updated additive CHARMM36/CGenFF and polarizable Drude FFs for the study of halogenated ligand-protein interactions in computer-aided drug design.

Quantum-mechanical investigation of tetrel bond characteristics based on the point-of-charge (PoC) approach

The point-of-charge (PoC) approach was employed to investigate the characteristics of the tetrel bond from an electrostatic perspective. W-T-XYZ···B nomenclature was suggested where T is a tetrel atom, W is the atom along the σ-hole extension, B is a Lewis base, and X, Y, and Z are three atoms on the same side of the σ-hole. Quantum-mechanical calculations were carried out on F-T-F₃ systems (where T = C, Si, Ge, or Sn) at the MP2/aug-cc-pVTZ level of theory, with PP functions for Ge and Sn atoms. The tetrel bond strength was estimated via the molecular stabilization energy. Tetrel bond strength was found to increase with increasing PoC negativity (i.e., Lewis basicity) and the electronegativity of the W atom. Moreover, the effects of the T···PoC distance, the W-T···PoC angle, and the aqueous medium on the tetrel bond strength were also investigated. Correlations between tetrel bond strength and several atomic and molecular descriptors such as the natural charge on the tetrel atom, E_HOMO, and the p-orbital contribution to W-T bond hybridization were observed. Contrary to expectations, the tetrel bond strength in F-C-X₃ increased as the electronegativity of X decreased. The σ-node criteria for the studied molecules were also introduced and discussed. The ability of these molecules to simultaneously form more than one tetrel bond was examined via the σⁿ-hole test. In conclusion, the tetrel bond strength was found to be governed by the strengths of (i) the attractive electrostatic interaction of the Lewis base with the σ-hole, (ii) the attractive/repulsive interaction between the Lewis base and the X, Y, and Z atoms, and (iii) the van der Waals interaction between the Lewis base and the X, Y, and Z atoms. Graphical Abstract Characterization of tetrel bond using the Point-of-Charge (PoC) approach.