Jindřich Fanfrlík

Halogen bond tunability I: the effects of aromatic fluorine substitution on the strengths of halogen-bonding interactions involving chlorine, bromine, and iodine

AbstractIn the past several years, halogen bonds have been shown to be relevant in crystal engineering and biomedical applications. One of the reasons for the utility of these types of noncovalent interactions in the development of, for example, pharmaceutical ligands is that their strengths and geometric properties are very tunable. That is, substitution of atoms or chemical groups in the vicinity of a halogen can have a very strong effect on the strength of the halogen bond. In this study we investigate halogen-bonding interactions involving aromatically-bound halogens (Cl, Br, and I) and a carbonyl oxygen. The properties of these halogen bonds are modulated by substitution of aromatic hydrogens with fluorines, which are very electronegative. It is found that these types of substitutions have dramatic effects on the strengths of the halogen bonds, leading to interactions that can be up to 100% stronger. Very good correlations are obtained between the interaction energies and the magnitudes of the positive electrostatic potentials (σ-holes) on the halogens. Interestingly, it is seen that the substitution of fluorines in systems containing smaller halogens results in electrostatic potentials resembling those of systems with larger halogens, with correspondingly stronger interaction energies. It is also shown that aromatic fluorine substitutions affect the optimal geometries of the halogen-bonded complexes, often as the result of secondary interactions. FigureSchematic models of halogen bonding complexes studied in this work

Semiempirical Quantum Chemical PM6 Method Augmented by Dispersion and H-Bonding Correction Terms Reliably Describes Various Types of Noncovalent Complexes.

Because of its construction and parametrization for more than 80 elements, the semiempirical quantum chemical PM6 method is superior to other similar methods. Despite its advantages, however, the PM6 method fails for the description of noncovalent interactions, specifically the dispersion energy and H-bonding. Upon inclusion of correction terms for dispersion and H-bonding, the performance of the method was found to be dramatically improved. The former correction included two parameters in the damping function that were parametrized to reproduce the benchmark interaction energies [CCSD(T)/complete basis set (CBS) limit] of the dispersion-bonded complexes from the S22 data set. The latter correction was parametrized on an extended set of H-bonded stabilization energies determined at the MP2/cc-pVTZ level. The resulting PM6-DH method was tested on the S22 data set, for which chemical accuracy (error < 1 kcal/mol) was achieved, and also on the JSCH2005 set, for which significant improvement over the original PM6 method was also obtained. Implementation of analytical gradients allows very efficient geometry optimization, which, for all complexes, provides better agreement with the benchmark data. Excellent results were also achieved for small peptides, and here again, chemical accuracy was obtained (i.e., the error with respect to CCSD(T)/CBS results was smaller than 1 kcal/mol). The performance of the technique was finally demonstrated on extended complexes, namely, the porphine dimer and various graphene models with DNA bases and base pairs, where the PM6-DH stabilization energies agree very well with available benchmark data obtained with DFT-D, SCS-MP2, and MP2.5 methods. The PM6-DH calculations are very efficient and can be routinely applied for systems of up to 1000 atoms. For nonaromatic systems, the use of a linear scaling version of the SCF procedure based on localized orbitals speeds up the method significantly and allows one to investigate systems with several thousand atoms. The method can thus replace force fields, which face basic problems for the description of quantum effects, in many applications.

Halogen bond tunability II: the varying roles of electrostatic and dispersion contributions to attraction in halogen bonds

AbstractIn a previous study we investigated the effects of aromatic fluorine substitution on the strengths of the halogen bonds in halobenzene…acetone complexes (halo = chloro, bromo, and iodo). In this work, we have examined the origins of these halogen bonds (excluding the iodo systems), more specifically, the relative contributions of electrostatic and dispersion forces in these interactions and how these contributions change when halogen σ-holes are modified. These studies have been carried out using density functional symmetry adapted perturbation theory (DFT-SAPT) and through analyses of intermolecular correlation energies and molecular electrostatic potentials. It is found that electrostatic and dispersion contributions to attraction in halogen bonds vary from complex to complex, but are generally quite similar in magnitude. Not surprisingly, increasing the size and positive nature of a halogen’s σ-hole dramatically enhances the strength of the electrostatic component of the halogen bonding interaction. Not so obviously, halogens with larger, more positive σ-holes tend to exhibit weaker dispersion interactions, which is attributable to the lower local polarizabilities of the larger σ-holes. FigureIn this work we investigate the roles played by electrostatic and dispersion forces in stabilizing halogen bonding interactions.

The Dominant Role of Chalcogen Bonding in the Crystal Packing of 2D/3D Aromatics

The chalcogen bond is a nonclassical σ-hole-based noncovalent interaction with emerging applications in medicinal chemistry and material science. It is found in organic compounds, including 2D aromatics, but has so far never been observed in 3D aromatic inorganic boron hydrides. Thiaboranes, harboring a sulfur heteroatom in the icosahedral cage, are candidates for the formation of chalcogen bonds. The phenyl-substituted thiaborane, synthesized and crystalized in this study, forms sulfur⋅⋅⋅π type chalcogen bonds. Quantum chemical analysis revealed that these interactions are considerably stronger than both in their organic counterparts and in the known halogen bond. The reason is the existence of a highly positive σ-hole on the positively charged sulfur atom. This discovery expands the possibilities of applying substituted boron clusters in crystal engineering and drug design.

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.

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.

QM/MM calculations reveal the different nature of the interaction of two carborane-based sulfamide inhibitors of human carbonic anhydrase II.

The crystal structures of two novel carborane-sulfamide inhibitors in the complex with human carbonic anhydrase II (hCAII) have been studied using QM/MM calculations. Even though both complexes possess the strongly interacting sulfamide···zinc ion motif, the calculations have revealed the different nature of binding of the carborane parts of the inhibitors. The neutral closo-carborane cage was bound to hCAII mainly via dispersion interactions and formed only very weak dihydrogen bonds. On the contrary, the monoanionic nido cage interacted with the protein mainly via electrostatic interactions. It formed short and strong dihydrogen bonds (stabilization of up to 4.2 kcal/mol; H···H distances of 1.7 Å) with the polar hydrogen of protein NH2 groups. This type of binding is unique among all of the classical organic and inorganic inhibitors of hCAII. Virtual glycine scanning allowed us to identify the amino-acid side chains, which made important contributions to ligand-binding energies. In summary, using QM/MM calculations, we have provided a detailed understanding of the differences between the interactions of two carborane sulfamides, identified the amino acids of hCAII with which they interact, and thus paved the way for the computer-aided rational design of selective boron-cluster-containing hCAII inhibitors.

The Effect of Halogen-to-Hydrogen Bond Substitution on Human Aldose Reductase Inhibition

The effect of halogen-to-hydrogen bond substitution on the binding energetics and biological activity of a human aldose reductase inhibitor has been studied using X-ray crystallography, IC50 measurements, advanced binding free energy calculations, and simulations. The replacement of Br or I atoms by an amine (NH2) group has not induced changes in the original geometry of the complex, which made it possible to study the isolated features of selected noncovalent interactions in a biomolecular complex.

Interactions of Boranes and Carboranes with Aromatic Systems: CCSD(T) Complete Basis Set Calculations and DFT-SAPT Analysis of Energy Components

The noncovalent interactions of heteroboranes with aromatic systems have only recently been acknowledged as a source of stabilization in supramolecular complexes. The physical basis of these interactions has been studied in several model complexes using advanced computational methods. The highly accurate CCSD(T)/complete basis set (CBS) value of the interaction energy for the model diborane···benzene complex in a stacking geometry exhibiting a B(2)H···π hydrogen bond was calculated to be -4.0 kcal·mol(-1). The DFT-SAPT/CBS approach, which is shown to reproduce the CCSD(T)/CBS data reliably asserted that the major stabilizing component was dispersion, followed by electrostatics. Furthermore, the effect of the benzene heteroatom- and exosubstitutions was studied and found to be small. Next, when aromatic molecules were changed to cyclic aliphatic ones, van der Waals complexes stabilized by the dispersion term only were formed. As the last step, interactions of two larger icosahedral borane cages with benzene were explored. The complex of the monoanionic CB(11)H(12)(-) exhibited two minima: the first stacked above the plane of the benzene ring with a C-H···π hydrogen bond and the second planar, in which the carborane cage bound to benzene via five B-H···H-C dihydrogen bonds. The DFT-SAPT/CBS calculations revealed that both of these binding motifs were stabilized by dispersion followed by electrostatic terms, with the planar complex being 1.4 kcal·mol(-1) more stable than the stacked one. The dianionic B(12)H(12)(2-) interacted with benzene only in the planar geometry, similarly as smaller anions do. The large stabilization energy of 11.0 kcal·mol(-1) was composed of dominant attractive dispersion and slightly smaller electrostatic and induction terms. In summary, the borane/carborane···aromatic interaction is varied both in the complex geometries and in the stabilizing energy components. The detailed insight derived from high-level quantum chemical computations can help us understand such important processes as host-guest complexation or carborane···biomolecule interactions.

Chalcogen and Pnicogen Bonds in Complexes of Neutral Icosahedral and Bicapped Square-Antiprismatic Heteroboranes

A systematic quantum mechanical study of σ-hole (chalcogen, pnicogen, and halogen) bonding in neutral experimentally known closo-heteroboranes is performed. Chalcogens and pnicogens are incorporated in the borane cage, whereas halogens are considered as exo-substituents of dicarbaboranes. The chalcogen and pnicogen atoms in the heteroborane cages have partial positive charge and thus more positive σ-holes. Consequently, these heteroboranes form very strong chalcogen and pnicogen bonds. Halogen atoms in dicarbaboranes also have a highly positive σ-hole, but only in the case of C-bonded halogen atoms. In such cases, the halogen bond of heteroboranes is also strong and comparable to halogen bonds in organic compounds with several electron-withdrawing groups being close to the halogen atom involved in the halogen bond.