Perspectives of Halogen Bonding Description in Scoring Functions and QSAR/QSPR: Substituent Effects in Aromatic Core

Assessing How Residual Errors of Scoring Functions Correlate to Ligand Structural Features

Scoring functions (SFs) are ubiquitous tools for early stage drug discovery. However, their accuracy currently remains quite moderate. Despite a number of successful target-specific SFs appearing recently, up until now, no ideas on how to systematically improve the general scope of SFs have been formulated. In this work, we hypothesized that the specific features of ligands, corresponding to interactions well appreciated by medicinal chemists (e.g., hydrogen bonds, hydrophobic and aromatic interactions), might be responsible, in part, for the remaining SF errors. The latter provides direction to efforts aimed at the rational and systematic improvement of SF accuracy. In this proof-of-concept work, we took a CASF-2016 coreset of 285 ligands as a basis for comparison and calculated the values of scores for a representative panel of SFs (including AutoDock 4.2, AutoDock Vina, X-Score, NNScore2.0, ΔVina RF20, and DSX). The residual error of linear correlation of each SF value, with the experimental values of affinity and activity, was then analyzed in terms of its correlation with the presence of the fragments responsible for certain medicinal chemistry defined interactions. We showed that, despite the fact that SFs generally perform reasonably, there is room for improvement in terms of better parameterization of interactions involving certain fragments in ligands. Thus, this approach opens a potential way for the systematic improvement of SFs without their significant complication. However, the straightforward application of the proposed approach is limited by the scarcity of reliable available data for ligand-receptor complexes, which is a common problem in the field.

Neuropsychopharmacology of Emerging Drugs of Abuse: meta- and para-Halogen-Ring-Substituted α-PVP (“flakka”) Derivatives

Changes in the molecular structure of synthetic cathinones has led to an increase in the number of novel emerging drugs in the illicit drug market at an unprecedented rate. Unfortunately, little is known about the neuropsychopharmacology of recently emerged halogen-substituted α-PVP derivatives. Thus, the aim of this study was to investigate the role of para- and meta-halogen (F-, Cl-, and Br-) substitutions on the in vitro, in silico, and in vivo effects of α-pyrrolidinopentiophenone (α-PVP) derivatives. HEK293 cells expressing the human dopamine or serotonin transporter (hDAT and hSERT) were used for the uptake inhibition and transporter affinity assays. Molecular docking was used to model the interaction mechanism against DAT. Swiss CD-1 mice were used for the horizontal locomotor activity, open field test, and conditioned place preference paradigm. All compounds demonstrated potent DA uptake inhibition and higher DAT selectivity than cocaine. Meta-substituted cathinones showed higher DAT/SERT ratios than their para- analogs, which correlates with an increased psychostimulant effect in vivo and with different meta- and para-in silico interactions at DAT. Moreover, all compounds induced rewarding and acute anxiogenic effects in mice. In conclusion, the present study demonstrates the role of meta- and para-halogen substitutions in the mechanism of action and provides the first evidence of the rewarding and anxiety-like properties of halogenated α-PVP derivatives.

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.

Multipole models of sulphur for accurate anisotropic electrostatic interactions within force fields

Nowadays, as computing has become much more available, a fresh momentum has been observed in the field of re-visioning and re-parameterizing the usual tools, as well as estimating for the incorporation of new qualitative capabilities, aimed at making more accurate and reliable predictions in drug discovery processes. Inspired by the success of modelling the electrostatic part of the halogen bonding (XB) by means of the distributed multipole expansion, a study is presented which attempts to extend this approach to a tougher case of σ-hole interaction: sulphur-based chalcogen bonding. To that end, 11 anisotropic models have been derived and tested for their performance in the reproduction of reference ab initio molecular electrostatic potential. A careful examination resulted in three models which have been selected for further examination as a part of the molecular mechanics force field (GAFF). The combined force field was used to estimate inter- and intra-molecular interactions for the molecular systems, capable of differentiating the binding from the σ-hole and other directions. The anisotropic models proposed were generally able to correct the wrong predictions of the sulphur models based only on isotropic charges and, thus, are a promising direction for further development of the refined electrostatics force fields.

Repositioning organohalogen drugs: a case study for identification of potent B-Raf V600E inhibitors via docking and bioassay

Drug repositioning has been attracting increasingly attention for its advantages of reducing costs and risks. Statistics showed that around one quarter of the marketed drugs are organohalogens. However, no study has been reported, to the best of our knowledge, to aim at efficiently repositioning organohalogen drugs, which may be attributed to the lack of accurate halogen bonding scoring function. Here, we present a study to show that two organohalogen drugs were successfully repositioned as potent B-Raf V600E inhibitors via molecular docking with halogen bonding scoring function, namely D(3)DOCKxb developed in our lab, and bioassay. After virtual screening by D(3)DOCKxb against the database CMC (Comprehensive Medicinal Chemistry), 3 organohalogen drugs that were predicted to form strong halogen bonding with B-Raf V600E were purchased and tested with ELISA-based assay. In the end, 2 of them, rafoxanide and closantel, were identified as potent inhibitors with IC50 values of 0.07 μM and 1.90 μM, respectively, which are comparable to that of vemurafenib (IC50: 0.17 μM), a marketed drug targeting B-Raf V600E. Single point mutagenesis experiments confirmed the conformations predicted by D(3)DOCKxb. And comparison experiment revealed that halogen bonding scoring function is essential for repositioning those drugs with heavy halogen atoms in their molecular structures.

Halogen bonding in drug-like molecules: a computational and systematic study of the substituent effect

Systematic study of the effect of fourteen chemical groups at the ortho, para and meta positions of NMA⋯halobenzene complexes showed a significant influence on halogen bonding, and also non-additive effects. A comprehensive description is reported.