Halogen Bonding—A Novel Interaction for Rational Drug Design?

Can an entirely negative fluorine in a molecule, viz. perfluorobenzene, interact attractively with the entirely negative site(s) on another molecule(s)? Like liking like!

The study presents the possibility for the formation of attractive intermolecular interactions between various entirely negative Lewis bases, leading to the formation of the thirteen isolated binary complexes examined.

A halogen bond does not dictate the conformational preferences of cis-1,3-disubstituted cyclohexanes

A halogen bond is unlikely to control the conformation of cis-1,3-disubstituted cyclohexanes. In addition, polarization rather than the σ-hole effect describes this weak interaction more appropriately.

Predictive Models for Halogen-bond Basicity of Binding Sites of Polyfunctional Molecules

Halogen bonding (XB) strength assesses the ability of an electron-enriched group to be involved in complexes with polarizable electrophilic halogenated or diatomic halogen molecules. Here, we report QSPR models of XB of particular relevance for an efficient screening of large sets of compounds. The basicity is described by pKBI2 , the decimal logarithm of the experimental 1 : 1 (B : I2 ) complexation constant K of organic compounds (B) with diiodine (I2 ) as a reference halogen-bond donor in alkanes at 298 K. Modeling involved ISIDA fragment descriptors, using SVM and MLR methods on a set of 598 organic compounds. Developed models were then challenged to make predictions for an external test set of 11 polyfunctional compounds for which unambiguous assignment of the measured effective complexation constant to specific groups out of the putative acceptor sites is not granted. At this stage, developed models were used to predict pKBI2 of all putative acceptor sites, followed by an estimation of the predicted effective complexation constant using the ChemEqui program. The best consensus models perform well both in cross-validation (root mean squared error RMSE=0.39-0.47 logKBI2 units) and external predictions (RMSE=0.49). The SVM models are implemented on our website (http://infochim.u-strasbg.fr/webserv/VSEngine.html) together with the estimation of their applicability domain and an automatic detection of potential halogen-bond acceptor atoms.

Synergistic and diminutive effects between halogen bond and lithium bond in complexes involving aromatic compounds

Quantum chemical calculations have been performed to study the interplay between halogen bond and lithium bond in the ternary systems FX-C6H5CN-LiF, FLi-C6H5CN-XF, and FLi-C6H5X-NH3 (X = Cl, Br, and I) involving aromatic compounds. This effect was studied in terms of interaction energy, electron density, charge transfer, and orbital interaction. The results showed that both FX-C6H5CN-LiF and FLi-C6H5CN-XF exhibit diminutive effects with the weakening of halogen bond and lithium bond, while FLi-C6H5X-NH3 displays synergistic effects with the strengthening of halogen bond and lithium bond. The nature of halogen bond and lithium bond in the corresponding binary complexes has been unveiled by the quantum theory of atoms in molecules methodology and energy decomposition analysis.

Synthesis and In Vitro Activity of Polyhalogenated 2-phenylbenzimidazoles as a New Class of anti-MRSA and Anti-VRE Agents

A series of novel polyhalogenated 2-phenylbenzimidazoles have been synthesized and evaluated for in vitro antistaphylococcal activity against drug-resistant bacterial strains (methicillin-resistant Staphylococcus aureus, and vancomycin-resistant Enterococcus faecium. Certain compounds inhibit bacterial growth perfectly. 11 was active than vancomycin (0.78 μg/mL) with the lowest MIC values with 0.19 μg/mL against methicillin-resistant Staphylococcus aureus, 8 and 35 exhibited best inhibitory activity against vancomycin-resistant Enterococcus faecium (1.56 μg/mL). The mechanism of action for this class of compounds appears to be different than clinically used antibiotics. These polyhalogenated benzimidazoles have potential for further investigation as a new class of potent anti-methicillin-resistant Staphylococcus aureus and anti-vancomycin-resistant Enterococcus faecium agents.

Effect of cooperativity in lithium bonding on the strength of halogen bonding and tetrel bonding: (LiCN)n···ClYF3 and (LiCN)n···YF3Cl (Y = C, Si and n = 1-5) complexes as a working model

This paper reports results of cooperativity in lithium bonding on the strength of halogen bonding and tetrel bonding in complexes pairing CF3Cl and SiF3Cl with (LiCN)n complexes, where n varies from 1 to 5. Molecular geometries and stabilization energies of title complexes are calculated at the MP2 level with 6-311++G(d,p) basis set. Cooperative effects are found in terms of structural and energetic properties when lithium, halogen, and tetrel bonds are present in these complexes simultaneously. Our results reveal that strength of halogen and tetrel bondings are enhanced due to cooperativity of Li···N interactions in lithium bonded complexes. Good linear correlations between cooperativity parameters and electronic properties of complexes were established in the present study.

Novel benzidine and diaminofluorene prolinamide derivatives as potent hepatitis C virus NS5A inhibitors

Our study describes the discovery of a series of highly potent hepatitis C virus (HCV) NS5A inhibitors based on symmetrical prolinamide derivatives of benzidine and diaminofluorene. Through modification of benzidine, l-proline, and diaminofluorene derivatives, we developed novel inhibitor structures, which allowed us to establish a library of potent HCV NS5A inhibitors. After optimizing the benzidine prolinamide backbone, we identified inhibitors embedding meta-substituted benzidine core structures that exhibited the most potent anti-HCV activities. Furthermore, through a battery of studies including hERG ligand binding assay, CYP450 binding assay, rat plasma stability test, human liver microsomal stability test, and pharmacokinetic studies, the identified compounds 24, 26, 27, 42, and 43 are found to be nontoxic, and are expected to be effective therapeutic anti-HCV agents.

Uses of 3-(2-Bromoacetyl)-2H-chromen-2-one in the Synthesis of Heterocyclic Compounds Incorporating Coumarin: Synthesis, Characterization and Cytotoxicity

In this work, 3-bromoacetylcoumarin was used as the key starting material for the synthesis of pyran, pyridine, thiophene, thiazole and pyrazole derivatives through its reaction with different reagents. The structures of the newly synthesized compounds were confirmed on the basis of their spectral data and elemental analyses. All of the synthesized compounds were screened for their in vitro anticancer activity against six human cancer cell lines, namely: human gastric cancer (NUGC), human colon cancer (DLD1), human liver cancer (HA22T and HEPG2), nasopharyngeal carcinoma (HONE1), human breast cancer (MCF) and normal fibroblast cells (WI38). The IC₅₀ values (the sample concentration that produces 50% reduction in cell growth) in nanomolars (nM)) showed most of the compounds exhibited significant cytotoxic effect. Among these derivatives, compound 6d showed almost equipotent cytotoxic activity against NUGC (IC₅₀ = 29 nM) compared to the standard CHS 828 (IC₅₀ = 25 nM).

Contrasting Anticancer Activity of Half-Sandwich Iridium(III) Complexes Bearing Functionally Diverse 2-Phenylpyridine Ligands

We report the synthesis, characterization, and antiproliferative activity of 15 iridium(III) half-sandwich complexes of the type [(η5-Cp*)Ir(2-(R'-phenyl)-R-pyridine)Cl] bearing either an electron-donating (-OH, -CH2OH, -CH3) or electron-withdrawing (-F, -CHO, -NO2) group at various positions on the 2-phenylpyridine (2-PhPy) chelating ligand giving rise to six sets of structural isomers. The X-ray crystal structures of [(η5-Cp*)Ir(2-(2'-fluorophenyl)pyridine)Cl] (1) and [(η5-Cp*)Ir(2-(4'-fluorophenyl)pyridine)Cl] (2) exhibit the expected "piano-stool" configuration. DFT calculations showed that substituents caused only localized effects on the electrostatic potential surface of the chelating 2-PhPy ligand of the complexes. Hydrolysis of all complexes is rapid, but readily reversed by addition of NaCl. The complexes show preferential binding to 9-ethylguanine over 9-methyladenine and are active catalysts for the oxidation of NADH to NAD+. Antiproliferative activity experiments in A2780 ovarian, MCF-7 breast, A549 lung, and HCT116 colon cancer cell lines showed IC50 values ranging from 1 to 89 μM, with the most potent complex, [(η5-Cp*)Ir(2-(2'-methylphenyl)pyridine)Cl] (13) (A2780 IC50 = 1.18 μM), being 10× more active than the parent, [(η5-Cp*)Ir(2-phenylpyridine)Cl], and 2× more active than [(η5-CpxPh)Ir(2-phenylpyridine)Cl]. Intriguingly, contrasting biological activities are observed between structural isomers despite exhibiting similar chemical reactivity. For pairs of structural isomers both the nature and position of the functional group can affect the hydrophobicity of the complex. An increase in hydrophobicity resulted in enhanced cellular-iridium accumulation in A2780 ovarian cells, which generally gave rise to an increase in potency. The structural isomers [(η5-Cp*)Ir(2-(4'-fluorophenyl)pyridine)Cl] (2) and [(η5-Cp*)Ir(2-phenyl-5-fluoropyridine)Cl] (4) preferentially localized in the cytosol > membrane and particulate > nucleus > cytoskeleton. This work highlights the strong dependence of biological behavior on the nature and position of the substituent on the chelating ligand and shows how this class of organometallic anticancer complexes can be fine-tuned to increase their potency without using extended cyclopentadienyl systems.

Substituted Bis-THF Protease Inhibitors with Improved Potency against Highly Resistant Mature HIV-1 Protease PR20

An extremely drug resistant mutant of HIV-1 protease (PR) bearing 20 mutations (PR20) has been studied with two potent antiviral investigational inhibitors. GRL-5010A and GRL-4410A were designed to introduce hydrogen bond interactions with the flexible flaps of the PR by incorporating gem-difluorines and alkoxy, respectively, at the C4 position of the bis-THF of darunavir. PR20 provides an excellent model for high level resistance, since clinical inhibitors are >1000-fold less active on PR20 than on wild-type enzyme. GRL-5010A and GRL-4410A show inhibition constants of 4.3 ± 7.0 and 1.7 ± 1.8 nM, respectively, for PR20, compared to the binding affinity of 41 ± 1 nM measured for darunavir. Crystal structures of PR20 in complexes with the two inhibitors confirmed the new hydrogen bond interactions with Gly 48 in the flap of the enzyme. The two new compounds are more effective than darunavir in inhibiting mature PR20 and show promise for further development of antiviral agents targeting highly resistant PR mutants.

The enhancing effects of group V σ-hole interactions on the F···O halogen bond

The σ-hole interaction, which occurs between the covalent IV-VII atoms and nucleophilic substances, has become a hot issue of weak interaction. In this work, NCF···O=PX3···(NCF)n (X = F, Cl, Br, H, CH3·; n = 0, 1, 2) complexes were constructed and studied based on the second-order Møller-Plesset perturbation theory (MP2) calculations to investigate the enhancing effects of group V σ-hole interactions on the F···O halogen bond. With increasing n, the FO halogen bond becomes stronger, indicating that the group V σ-hole interactions could enhance the F···O halogen bond. As the capacity of donating electrons of X increases, the most negative electrostatic potentials outside the oxygen atom of O=PX3···(NCF)n (n = 0, 1, 2) become more negative, resulting in a stronger F···O halogen bond. In the formation of a F···O halogen bond, along the sequence of X = F, Cl, Br, H, CH3 of the negative sites O=PX3, the electric field of the lone pair of oxygen becomes greater and causes a larger decrease in electron density outside the fluorine atom. On the other hand, with increasing n from 0 to 2, the group V σ-hole interactions also increase the electric field of lone pair of oxygen and results in a larger decrease in electron density outside the fluorine atom.

Extended halogen bonding between fully fluorinated aromatic molecules

Halogen bonding is a noncovalent interaction where an electrophilic cap on a halogen atom, the so-called σ-hole, attracts a nucleophilic site on an adjacent molecule. The polarizability of halogens relates to the strength of the σ-hole, and accordingly the halogen-halogen distance becomes shorter in the order of Cl, Br, and I. Fully fluoro-substituted aromatic molecules, on the contrary, are generally believed not to form halogen bonds due to the absence of a σ-hole. Here, we study atomic-scale in-plane F-F contacts with high-resolution force microscopy. Our ab initio calculations show that the attractive dispersion forces can overcome the electrostatic repulsion between the fluorine atoms, while the anisotropic distribution of the negative electrostatic potential leads the directional bond and even changes the gap. The coexistence of these two competing forces results in the formation of a "windmill" structure, containing three C-F···F bonds among neighboring molecules. While the σ-hole is absent, the scheme of the C-F···F bonding has a high similarity to halogen bonding.

An interaction map of small-molecule kinase inhibitors with anaplastic lymphoma kinase (ALK) mutants in ALK-positive non-small cell lung cancer

Human anaplastic lymphoma kinase (ALK) has become a well-established target for the treatment of ALK-positive non-small cell lung cancer (NSCLC). Here, we have profiled seven small-molecule inhibitors, including 2 that are approved drugs, against a panel of clinically relevant mutations in ALK tyrosine kinase (TK) domain, aiming at a comprehensive understanding of molecular mechanism and biological implication underlying inhibitor response to ALK TK mutation. We find that (i) the gatekeeper mutation L1196M causes crizotinib resistance by simultaneously increasing and decreasing the binding affinities of, respectively, ATP and inhibitor to ALK, whereas the secondary mutation C1156Y, which is located far away from the ATP-binding site of ALK TK domain, causes the resistance by inducing marked allosteric effect on the site, (ii) the 2nd and 3rd generation kinase inhibitors exhibit relatively high sensitivity towards ALK mutants as compared to 1st generation inhibitors, (iii) the pan-kinase inhibitor staurosporine is insensitive for most mutations due to its high structural compatibility, and (iv) ATP affinity to ALK is generally reduced upon most clinically relevant mutations. Furthermore, we also identify six novel mutation-inhibitor pairs that are potentially associated with drug resistance. In addition, the G1202R and C1156Y mutations are expected to generally cause resistance for many existing inhibitors, since they can address significant effect on the geometric shape and physicochemical property of ALK active pocket.

Comparing sixteen scoring functions for predicting biological activities of ligands for protein targets

Accurately predicting relative binding affinities and biological potencies for ligands that interact with proteins remains a significant challenge for computational chemists. Most evaluations of docking and scoring algorithms have focused on enhancing ligand affinity for a protein by optimizing docking poses and enrichment factors during virtual screening. However, there is still relatively limited information on the accuracy of commercially available docking and scoring software programs for correctly predicting binding affinities and biological activities of structurally related inhibitors of different enzyme classes. Presented here is a comparative evaluation of eight molecular docking programs (Autodock Vina, Fitted, FlexX, Fred, Glide, GOLD, LibDock, MolDock) using sixteen docking and scoring functions to predict the rank-order activity of different ligand series for six pharmacologically important protein and enzyme targets (Factor Xa, Cdk2 kinase, Aurora A kinase, COX-2, pla2g2a, β Estrogen receptor). Use of Fitted gave an excellent correlation (Pearson 0.86, Spearman 0.91) between predicted and experimental binding only for Cdk2 kinase inhibitors. FlexX and GOLDScore produced good correlations (Pearson>0.6) for hydrophilic targets such as Factor Xa, Cdk2 kinase and Aurora A kinase. By contrast, pla2g2a and COX-2 emerged as difficult targets for scoring functions to predict ligand activities. Although possessing a high hydrophobicity in its binding site, β Estrogen receptor produced reasonable correlations using LibDock (Pearson 0.75, Spearman 0.68). These findings can assist medicinal chemists to better match scoring functions with ligand-target systems for hit-to-lead optimization using computer-aided drug design approaches.

Spectral and thermal characterization of halogen-bonded novel crystalline oligo(p-bromoacetophenone formaldehyde)

A novel oligomer p-bromoacetophenone-formaldehyde (OPBAF) was prepared by condensation polymerization in the presence of an acid as catalyst. It was characterized by FT-IR, NMR, pyrolysis GC/MS, XRD, GPC, and TG-DTG. The crystallographic parameters and space group for hexagonal OPBAF were a = b = 2.0810 Å and c = 9.2340 Å and P3̅m1, respectively. The degradation activation energy of the oligomer was studied by the Kissinger method. The kinetic parameters were also obtained. Halogen bonding interactions in the crystalline oligomers are identified between halogen···carbonyl and halogen···halogen. Little correlation was found in the halogen bonding motifs exhibited as a function of bromine present in this oligomer, and a unique bifurcated Br···Br/Br···O═C halogen bonding synthon was identified. This newly developed oligomer may be used as an interesting material for the development of 3D-designed structural products.

Molecular recognition in chemical and biological systems

Structure-based ligand design in medicinal chemistry and crop protection relies on the identification and quantification of weak noncovalent interactions and understanding the role of water. Small-molecule and protein structural database searches are important tools to retrieve existing knowledge. Thermodynamic profiling, combined with X-ray structural and computational studies, is the key to elucidate the energetics of the replacement of water by ligands. Biological receptor sites vary greatly in shape, conformational dynamics, and polarity, and require different ligand-design strategies, as shown for various case studies. Interactions between dipoles have become a central theme of molecular recognition. Orthogonal interactions, halogen bonding, and amide⋅⋅⋅π stacking provide new tools for innovative lead optimization. The combination of synthetic models and biological complexation studies is required to gather reliable information on weak noncovalent interactions and the role of water.

2-Amino-4-arylthiazole Derivatives as Anti-giardial Agents: Synthesis, Biological Evaluation and QSAR Studies

A series of seven 2-amino-4-arylthiazoles were prepared following Hantzsch’s modified method under microwave irradiation. A set of 50 derivatives was obtained and the in vitro activity against Giardia intestinalis was evaluated. The results on the biological activity revealed that, in general, the N-(5-bromo-4-aryl-thiazol-2-yl)-acetamide scaffold showed high bioactivity. In particular, compounds 6e (IC50 = 0.39 μM) and 6b (IC50 = 0.87 μM) were found to be more potent than the positive control metronidazole. Citoxicity and acute toxicity tests performed showed low toxicity and high selectivity of the most active compounds (6e SI = 139, 6b SI = 52.3). A QSAR analysis was applied to a data set of 37 obtained 2-amino-4-arylthiazoles derivatives and the best model described a strongly correlation between the anti-giardiasic activity and molecular descriptors as E2M, RDF115m, F10, MATS6v, and Hypnotic-80, with high statistical quality. This finding indicates that N-substituted aminothiazole scaffold should be investigated for the development of highly selective anti-giardial agent.

Theoretical calculations of π-type pnicogen bonds in the triad intermolecular complexes

The pnicogen bonding interactions of PCl3and π-electron systems (acetylene, ethylene, benzene) were calculated by using MP2/aug-cc-pVDZ method and the effect of hydrogen bond on pnicogen bond systems were investigated. It has been indicated that the hydrogen bonding and the pnicogen bonding interactions have influence on each other and the positively cooperative effect has been detected. The interaction energies of pnicogen bonded supramolecular system were also calculated by using DFT method (M06-2X) and some simple comparisons with those by using MP2 method were made. It has been disclosed from natural bond orbitals (NBO) analysis that more the amount of charge transfer of pnicogen bonding interaction, the greater the stability of the corresponding complex. Through AIM topological analysis, it has been revealed that the electron density of pnicogen bond BCP point is positively correlated with the stability of trimeric complex. Electron localization function (ELF) was also adopted to analyze the nature of pnicogen bonding interactions. Furthermore, density difference function (DDF) method was adopted to analyze the variation of electron density of pnicogen bond system because of hydrogen bond.

Combining coordination and supramolecular chemistry to explore uranyl assembly in the solid state

Supramolecular assembly of uranyl species via halogen–oxo and halogen–halogen interactions is explored in the solid state.

Ether modifications to 1-[2-(3,4-dimethoxyphenyl)ethyl]-4-(3-phenylpropyl)piperazine (SA4503): effects on binding affinity and selectivity for sigma receptors and monoamine transporters

Two series of novel ether analogs of the sigma (σ) receptor ligand 1-[2-(3,4-dimethoxyphenyl)ethyl]-4-(3-phenylpropyl)piperazine (SA4503) have been prepared. In one series, the alkyl portion of the 4-methoxy group was replaced with allyl, propyl, bromoethyl, benzyl, phenethyl, and phenylpropyl moieties. In the second series, the 3,4-dimethoxy was replaced with cyclic methylenedioxy, ethylenedioxy and propylenedioxy groups. These ligands, along with 4-O-des-methyl SA4503, were evaluated for σ1 and σ2 receptor affinity, and compared to SA4503 and several known ether analogs. SA4503 and a subset of ether analogs were also evaluated for dopamine transporter (DAT) and serotonin transporter (SERT) affinity. The highest σ1 receptor affinities, Ki values of 1.75-4.63 nM, were observed for 4-O-des-methyl SA4503, SA4503 and the methylenedioxy analog. As steric bulk increased, σ1 receptor affinity decreased, but only to a point. Allyl, propyl and bromoethyl substitutions gave σ1 receptor Ki values in the 20-30 nM range, while bulkier analogs having phenylalkyl, and Z- and E-iodoallyl, ether substitutions showed higher σ1 affinities, with Ki values in the 13-21 nM range. Most ligands studied exhibited comparable σ1 and σ2 affinities, resulting in little to no subtype selectivity. SA4503, the fluoroethyl analog and the methylenedioxy congener showed modest six- to fourteen-fold selectivity for σ1 sites. DAT and SERT interactions proved much more sensitive than σ receptor interactions to these structural modifications. For example, the benzyl congener (σ1Ki=20.8 nM; σ2Ki=16.4 nM) showed over 100-fold higher DAT affinity (Ki=121 nM) and 6-fold higher SERT affinity (Ki=128nM) than the parent SA4503 (DAT Ki=12650 nM; SERT Ki=760 nM). Thus, ether modifications to the SA4503 scaffold can provide polyfunctional ligands having a broader spectrum of possible pharmacological actions.

Unusual bonding modes of perfluorobenzene in its polymeric (dimeric, trimeric and tetrameric) forms: entirely negative fluorine interacting cooperatively with entirely negative fluorine

The Laplacian of the charge density plot for the tetramer of perfluorobenzene, displaying the attraction between entirely negative fluorine atoms, as well as the windmill type intermolecular topology between those atoms.

Unstable, metastable, or stable halogen bonding interaction involving negatively charged donors? A statistical and computational chemistry study

The noncovalent halogen bonding could be attributed to the attraction between the positively charged σ-hole and a nucleophile. Quantum mechanics (QM) calculation indicated that the negatively charged organohalogens have no positively charged σ-hole on their molecular surface, leading to a postulation of repulsion between negatively charged organohalogens and nucleophiles in vacuum. However, PDB survey revealed that 24% of the ligands with halogen bonding geometry could be negatively charged. Moreover, 36% of ionizable drugs in CMC (Comprehensive Medicinal Chemistry) are possibly negatively charged at pH 7.0. QM energy scan showed that the negatively charged halogen bonding is probably metastable in vacuum. However, the QM calculated bonding energy turned negative in various solvents, suggesting that halogen bonding with negatively charged donors should be stable in reality. Indeed, QM/MM calculation on three crystal structures with negatively charged ligands revealed that the negatively charged halogen bonding was stable. Hence, we concluded that halogen bonding with negatively charged donors is unstable or metastable in vacuum but stable in protein environment, and possesses similar geometric and energetic characteristics as conventional halogen bonding. Therefore, negatively charged organohalogens are still effective halogen bonding donors for medicinal chemistry and other applications.