Assessment of Halogen Off-Center Point-Charge Models Using Explicit Solvent Simulations

Compounds containing halogens can form halogen bonds (XBs) with biological targets such as proteins and membranes due to their anisotropic electrostatic potential. To accurately describe this anisotropy, off-center point-charge (EP) models are commonly used in force field methods, allowing the description of XBs at the molecular mechanics and molecular dynamics level. Various EP implementations have been documented in the literature, and despite being efficient in reproducing protein-ligand geometries and sampling of XBs, it is unclear how well these EP models predict experimental properties such as hydration free energies (ΔGhyd), which are often used to validate force field performance. In this work, we report the first assessment of three EP models using alchemical free energy calculations to predict ΔGhyd values. We show that describing the halogen anisotropy using some EP models can lead to a slight improvement in the prediction of the ΔGhyd when compared with the models without EP, especially for the chlorinated compounds; however, this improvement is not related to the establishment of XBs but is most likely due to the improvement of the sampling of hydrogen bonds. We also highlight the importance of the choice of the EP model, especially for the iodinated molecules, since a slight tendency to improve the prediction is observed for compounds with a larger σ-hole but significantly worse results were obtained for compounds that are weaker XB donors.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Enantioseparation of planar chiral ferrocenes on cellulose-based chiral stationary phases: Benzoate versus carbamate pendant groups

In this study, the enantioseparation of 14 planar chiral ferrocenes containing halogen atoms, and methyl, iodoethynyl, phenyl, and 2-naphthyl groups, as substituents, was explored with a cellulose tris(4-methylbenzoate) (CMB)-based chiral column under multimodal elution conditions. n-Hexane/2-propanol (2-PrOH) 95:5 v/v, pure methanol (MeOH), and MeOH/water 90:10 v/v were used as mobile phases (MPs). With CMB, baseline enantioseparations were achieved for nine analytes with separation factors (α) ranging from 1.24 to 1.77, whereas only three analytes could be enantioseparated with 1.14 ≤ α ≤ 1.51 on a cellulose tris(3,5-dimethylphenylcarbamate) (CDMPC)-based column, used as a reference for comparison, under the same elution conditions. Pendant group-dependent reversal of the enantiomer elution order was observed in several cases by changing CMB to CDMPC. The impact of analyte and chiral stationary phase (CSP) structure, and MP polarity on the enantioseparation, was evaluated. The two cellulose-based CSPs featured by different pendant groups were also compared in terms of thermodynamics. For this purpose, enthalpy (ΔΔH°), entropy (ΔΔS°) and free energy (ΔΔG°) differences, isoenantioselective temperatures (T_iso ), and enthalpy/entropy ratios (Q), associated with the enantioseparations, were derived from van 't Hoff plots by using n-hexane/2-PrOH 95:5 v/v and methanol/water 90:10 v/v as MPs. With the aim to disclose the functions of the different substituents in mechanisms and noncovalent interactions underlying analyte-selector complex formation at molecular level, electrostatic potential (V) analysis and molecular dynamics simulations were used as computational techniques. On this basis, enantioseparations and related mechanisms were investigated by integrating theoretical and experimental data.

Targeting Protein Pockets with Halogen Bonds: The Role of the Halogen Environment

Protein pockets that form a halogen bond (X-bond) with a halogenated ligand molecule simultaneously form other (mainly hydrophobic) interactions with the halogen atom that can be considered as its "X-bond environment" (XBenv). Most studies in the field have focused on the X-bond, with the properties of the XBenv usually overlooked. In this work, we derived a protocol that evaluates the XBenv strength as a measure of the propensity of a protein pocket to host an X-bond. The charge density-based topological descriptors in combination with machine learning tools were employed to predict formation and strength of the interactions that conform the XBenv as a function of their geometrical parameters. On the basis of these results, we propose that the XBenv can be used as a footprint to judge the chance of a protein pocket to form an X-bond.

In silico investigations of some Cyperus rotundus compounds as potential anti-inflammatory inhibitors of 5-LO and LTA4H enzymes

The present study aimed to experimentally identify the essential oil of Algerian Cyperus rotundus L. and to model the interaction of some known anti-inflammatory molecules with two key enzymes involved in inflammation, 5-Lypoxygenase (5-LO) and leukotriene A4 hydrolase (LTA4H). Gas chromatography/gas chromatography-mass spectrometry (GC/GC-MS) revealed that 92.7% of the essential oil contains 35 compounds, including oxygenated sesquiterpenes (44.2%), oxygenated monoterpenes (30.2%), monoterpene hydrocarbons (11.8%) and sesquiterpene hydrocarbons (6.5%). The major identified oxygenated terpenes are humulene oxide II, caryophyllene oxide, khusinol, agarospirol, spathulinol and trans-pinocarveol Myrtenol and α-terpineol are known to exhibit anti-inflammatory activities. Several complexes obtained after docking the natural terpenes with 5-LO and LTA4H have shown strong hydrogen bonding interactions. The best docking energies were found with α-terpineol, Myrtenol and khusinol. The interaction between the natural products and amino-acid residues HIS367, ILE673 and GLN363 appears to be critical for 5-LO inhibition, while the interaction with residues GLU271, HIS295, TYR383, TYR378, GLU318, GLU296 and ASP375 is critical for LTA4H inhibition. Molecular dynamics (MD) trajectories of the selected docked complexes showed stable backbone root mean square deviation (RMSD), supporting the stability of the natural product-enzyme interaction.Communicated by Ramaswamy H. Sarma.

1,5-Dichloroethanoanthracene Derivatives As Antidepressant Maprotiline Analogs: Synthesis, DFT Computational Calculations, and Molecular Docking

The chlorinated tetracyclic 1,5-dichloro-9,10-dihydro-9,10-ethanoanthracen-12-yl)-N-methylmethanamine 1, a maprotiline analog, has been synthesized via reduction and the Diels–Alder reaction followed by reductive amination of aldehyde 2.1D-NMR (DEPT) and 2D-NMR (HSQC, DQF-COSY) techniques were recruited for structural elucidation in addition to HRMS. Density functional theory calculations were performed to identify the possible isomers of the intermediate compound aldehyde 2; these calculations were in good agreement with experimental results where aldehyde 2 could exist in three isomers with comparable energies. In addition, the side chain of this aldehyde 2 was extended via the Wittig reaction to obtain the unsaturated ester 5 that was subjected to selective olefinic catalytic hydrogenation to obtain the corresponding saturated ester 6. Molecular docking simulation showed that all the compounds (1, 2, 5, and 6) have high antidepressant activities and form stable complexes with LeuT by inhibiting the neurotransmitter reuptake at the synapse and hence are good candidates as antidepressant drugs.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

Stereoselective Processes Based on σ-Hole Interactions

The σ-hole interaction represents a noncovalent interaction between atoms with σ-hole(s) on their surface (such as halogens and chalcogens) and negative sites. Over the last decade, significant developments have emerged in applications where the σ-hole interaction was demonstrated to play a key role in the control over chirality. The aim of this review is to give a comprehensive overview of the current advancements in the use of σ-hole interactions in stereoselective processes, such as formation of chiral supramolecular assemblies, separation of enantiomers, enantioselective complexation and asymmetric catalysis.

The effect of off-center σ -hole on the atom-centered partial charges in halogenated molecules

Partial atomic charges belong to key concepts of computational chemistry. In some cases, however, they fail in describing the electrostatics of molecules. One such example is the σ -hole, a region of positive electrostatic potential located on halogens and other atoms. In molecular mechanics, the σ -hole is often modeled as a pseudo-atom with a positive partial charge located off the halogen nucleus. Here we address a question, to what extent the pseudo-atom affects partial charges of other atoms in the molecule. To this aim, we have thoroughly analyzed partial charges of over 2300 halogenated molecules from the ZINC database calculated by the restricted electrostatic potential (RESP) method and compared them with the charges fitted by RESP including the pseudo-atom. We show that the pseudo-atom improves charge fitting for a vast majority of molecules. The σ -hole, modeled as the off-center charge, affects the atoms within three covalent bonds from the halogen.

Exploring Orthogonality between Halogen and Hydrogen Bonding Involving Benzene

The concept of orthogonality between halogen and hydrogen bonding, brought out by Ho and coworkers some years ago, has become a widely accepted idea within the chemists' community. While the original work was based on a common carbonyl oxygen as acceptor for both interactions, we explore here, by means of M06-2X, M11, ωB97X, and ωB97XD/aug-cc-PVTZ DFT calculations, the interdependence of halogen and hydrogen bonding with a shared π-electron system of benzene. The donor groups (specifically NCBr and H2O) were placed on either or the same side of the ring, according to a double T-shaped or a perpendicular geometry, respectively. The results demonstrate that the two interactions with benzene are not strictly independent on each other, therefore outlining that the orthogonality between halogen and hydrogen bonding, intended as energetical independence between the two interactions, should be carefully evaluated according to the specific acceptor group.

Biological halogen bonds in protein-ligand complexes: a combined QTAIM and NCIPlot study in four representative cases

In this study, the PDB has been manually scrutinized by using a subset of all PDB entries containing organic iodinated ligands. Four structures exhibiting short IA halogen bonding (HaB) contacts (A stands for the σ-hole acceptor) have been selected and further analysed. In most hits, the sigma-hole acceptor corresponds to an O-atom of the amido group belonging to the protein backbone. In a minority of hits, the electron donors are O, S, Se or π-systems of the amino-acid side chains. A judicious selection of four PDB structures presenting all four types of HaB interactions (C-IA, A = O, S, Se, π) has been performed. For these selected structures, a comprehensive RI-MP2/def2-TZVP study has been carried out to evaluate the HaB energetically. Moreover, the interactions have been characterized by combining the quantum theory of "atoms-in-molecules" (QTAIM) and the noncovalent interaction plot (NCIPlot) and rationalized using the molecular electrostatic potential (MEP) surface.

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.

Enantioseparations of polyhalogenated 4,4'-bipyridines on polysaccharide-based chiral stationary phases and molecular dynamics simulations of selector-selectand interactions

2'-(4-Pyridyl)- and 2'-(4-hydroxyphenyl)-TCIBPs (TCIBP = 3,3',5,5'-tetrachloro-2-iodo-4,4'-bipyridyl) are chiral compounds that showed interesting inhibition activity against transthyretin fibrillation in vitro. We became interested in their enantioseparation since we noticed that the M-stereoisomer is more effective than the P-enantiomer. Based thereon, we recently reported the enantioseparation of 2'-substituted TCIBP derivatives with amylose-based chiral columns. Following this study, herein we describe the comparative enantioseparation of both 2'-(4-pyridyl)- and 2'-(4-hydroxyphenyl)-TCIBPs on four cellulose phenylcarbamate-based chiral columns aiming to explore the effect of the polymer backbone, as well as the nature and position of substituents on the side groups on the enantioseparability of these compounds. In the frame of this project, the impact of subtle variations of analyte and polysaccharide structures, and mobile phase (MP) polarity on retention and selectivity was evaluated. The effect of temperature on retention and selectivity was also considered, and overall thermodynamic parameters associated with the analyte adsorption onto the CSP surface were derived from van 't Hoff plots. Interesting cases of enantiomer elution order (EEO) reversal were observed. In particular, the EEO was shown to be dependent on polysaccharide backbone, the elution sequence of the two analytes being P-M and M-P on cellulose and amylose tris(3,5-dimethylphenylcarbamate), respectively. In this regard, a theoretical investigation based on molecular dynamics (MD) simulations was performed by using amylose and cellulose tris(3,5-dimethylphenylcarbamate) nonamers as virtual models of the polysaccharide-based selectors. This exploration at the molecular level shed light on the origin of the enantiodiscrimination processes.

New Aspects of Monoamine Oxidase B Inhibitors: The Key Role of Halogens to Open the Golden Door

A large plethora of drugs and promising lead compounds contain halogens in their structures. The introduction of such moieties strongly modulates their physical-chemical features as well as pharmacokinetic and pharmacodynamic profile. The most important outcome was shown to be the ability of these halogens to favourably influence the drug-target interaction and energetic stability within the active site by the establishment of halogen bonds. This review attempted to demonstrate the key role exerted by these versatile moieties when correctly located in an organic scaffold to display Monoamine Oxidase (MAO) inhibition and selectivity towards the B isoform of this important enzyme. Human MAOs are well-recognized as therapeutic targets for mood disorders and neurodegenerative diseases and medicinal chemists were prompted to discover the structural requirements crucial to discriminate the slight differences between the active sits of the two isoforms (MAO-A and MAOB). The analysis of the structure-activity relationships of the most important scaffolds (hydrazothiazoles, coumarins, chromones, chalcones, pyrazolines) and the impact of halogen (F, Cl, Br and I) insertion on this biological activity and isozyme selectivity have been reported being a source of inspiration for the medicinal chemists.

Surface-controlled reversal of the selectivity of halogen bonds

Intermolecular halogen bonds are ideally suited for designing new molecular assemblies because of their strong directionality and the possibility of tuning the interactions by using different types of halogens or molecular moieties. Due to these unique properties of the halogen bonds, numerous areas of application have recently been identified and are still emerging. Here, we present an approach for controlling the 2D self-assembly process of organic molecules by adsorption to reactive vs. inert metal surfaces. Therewith, the order of halogen bond strengths that is known from gas phase or liquids can be reversed. Our approach relies on adjusting the molecular charge distribution, i.e., the σ-hole, by molecule-substrate interactions. The polarizability of the halogen and the reactiveness of the metal substrate are serving as control parameters. Our results establish the surface as a control knob for tuning molecular assemblies by reversing the selectivity of bonding sites, which is interesting for future applications.

Understanding noncovalent bonds and their controlling forces

The fundamental underpinnings of noncovalent bonds are presented, focusing on the σ-hole interactions that are closely related to the H-bond. Different means of assessing their strength and the factors that control it are discussed. The establishment of a noncovalent bond is monitored as the two subunits are brought together, allowing the electrostatic, charge redistribution, and other effects to slowly take hold. Methods are discussed that permit prediction as to which site an approaching nucleophile will be drawn, and the maximum number of bonds around a central atom in its normal or hypervalent states is assessed. The manner in which a pair of anions can be held together despite an overall Coulombic repulsion is explained. The possibility that first-row atoms can participate in such bonds is discussed, along with the introduction of a tetrel analog of the dihydrogen bond.

Rational Design, Synthesis, Characterization and Evaluation of Iodinated 4,4'-Bipyridines as New Transthyretin Fibrillogenesis Inhibitors

The 3,3',5,5'-tetrachloro-2-iodo-4,4'-bipyridine structure is proposed as a novel chemical scaffold for the design of new transthyretin (TTR) fibrillogenesis inhibitors. In the frame of a proof-of-principle exploration, four chiral 3,3',5,5'-tetrachloro-2-iodo-2'-substituted-4,4'- bipyridines were rationally designed and prepared from a simple trihalopyridine in three steps, including a Cu-catalysed Finkelstein reaction to introduce iodine atoms on the heteroaromatic scaffold, and a Pd-catalysed coupling reaction to install the 2'-substituent. The corresponding racemates, along with other five chiral 4,4'-bipyridines containing halogens as substituents, were enantioseparated by high-performance liquid chromatography in order to obtain pure enantiomer pairs. All stereoisomers were tested against the amyloid fibril formation (FF) of wild type (WT)-TTR and two mutant variants, V30M and Y78F, in acid mediated aggregation experiments. Among the 4,4'-bipyridine derivatives, interesting inhibition activity was obtained for both enantiomers of the 3,3',5,5'-tetrachloro-2'-(4-hydroxyphenyl)-2-iodo-4,4'-bipyridine. In silico docking studies were carried out in order to explore possible binding modes of the 4,4'-bipyridine derivatives into the TTR. The gained results point out the importance of the right combination of H-bond sites and the presence of iodine as halogen-bond donor. Both experimental and theoretical evidences pave the way for the utilization of the iodinated 4,4'-bipyridine core as template to design new promising inhibitors of TTR amyloidogenesis.

A Halogen Bonding Perspective on Iodothyronine Deiodinase Activity

Iodothyronine deiodinases (Dios) are involved in the regioselective removal of iodine from thyroid hormones (THs). Deiodination is essential to maintain TH homeostasis, and disruption can have detrimental effects. Halogen bonding (XB) to the selenium of the selenocysteine (Sec) residue in the Dio active site has been proposed to contribute to the mechanism for iodine removal. Polybrominated diphenyl ethers (PBDEs) and polychlorinated biphenyls (PCBs) are known disruptors of various pathways of the endocrine system. Experimental evidence shows PBDEs and their hydroxylated metabolites (OH-BDEs) can inhibit Dio, while data regarding PCB inhibition are limited. These xenobiotics could inhibit Dio activity by competitively binding to the active site Sec through XB to prevent deiodination. XB interactions calculated using density functional theory (DFT) of THs, PBDEs, and PCBs to a methyl selenolate (MeSe⁻) arrange XB strengths in the order THs > PBDEs > PCBs in agreement with known XB trends. THs have the lowest energy C–X*-type unoccupied orbitals and overlap with the Se lp donor leads to high donor-acceptor energies and the greatest activation of the C–X bond. The higher energy C–Br* and C–Cl* orbitals similarly result in weaker donor-acceptor complexes and less activation of the C–X bond. Comparison of the I···Se interactions for the TH group suggest that a threshold XB strength may be required for dehalogenation. Only highly brominated PBDEs have binding energies in the same range as THs, suggesting that these compounds may inhibit Dio and undergo debromination. While these small models provide insight on the I···Se XB interaction itself, interactions with other active site residues are governed by regioselective preferences observed in Dios.

The Significance of Halogen Bonding in Ligand-Receptor Interactions: The Lesson Learned from Molecular Dynamic Simulations of the D4 Receptor

Recently, a computational approach combining a structure-activity relationship library containing pairs of halogenated ligands and their corresponding unsubstituted ligands (called XSAR) with QM-based molecular docking and binding free energy calculations was developed and used to search for amino acids frequently targeted by halogen bonding, also known as XB hot spots. However, the analysis of ligand-receptor complexes with halogen bonds obtained by molecular docking provides a limited ability to study the role and significance of halogen bonding in biological systems. Thus, a set of molecular dynamics simulations for the dopamine D4 receptor, recently crystallized with the antipsychotic drug nemonapride (5WIU), and the five XSAR sets were performed to verify the identified hot spots for halogen bonding, in other words, primary (V5x40), and secondary (S5x43, S5x461 and H6x55). The simulations confirmed the key role of halogen bonding with V5x40 and H6x55 and supported S5x43 and S5x461. The results showed that steric restrictions and the topology of the molecular core have a crucial impact on the stabilization of the ligand-receptor complex by halogen bonding.

Could Quantum Mechanical Properties Be Reflected on Classical Molecular Dynamics? The Case of Halogenated Organic Compounds of Biological Interest

Essential to understanding life, the biomolecular phenomena have been an important subject in science, therefore a necessary path to be covered to make progress in human knowledge. To fully comprehend these processes, the non-covalent interactions are the key. In this review, we discuss how specific protein-ligand interactions can be efficiently described by low computational cost methods, such as Molecular Mechanics (MM). We have taken as example the case of the halogen bonds (XB). Albeit generally weaker than the hydrogen bonds (HB), the XBs play a key role to drug design, enhancing the affinity and selectivity toward the biological target. Along with the attraction between two electronegative atoms in XBs explained by the σ-hole model, important orbital interactions, as well as relief of Pauli repulsion take place. Nonetheless, such electronic effects can be only well-described by accurate quantum chemical methods that have strong limitations dealing with supramolecular systems due to their high computational cost. To go beyond the poor description of XBs by MM methods, reparametrizing the force-fields equations can be a way to keep the balance between accuracy and computational cost. Thus, we have shown the steps to be considered when parametrizing force-fields to achieve reliable results of complex non-covalent interactions at MM level for In Silico drug design methods.

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.

Halogenated derivatives of methotrexate as human dihydrofolate reductase inhibitors in cancer chemotherapy

Methotrexate is a widely used anti-metabolite in cancer chemotherapy. A series of halogenated drugs is designed from Methotrexate to assess their interactions with human dihydrofolate reductase. The aim of this study is to evaluate the performance of the modified drugs compared to the parent Methotrexate. Density Functional Theory is employed to optimize these drugs. Molecular docking calculation of these optimized drugs against dihydrofolate reductase is performed to find out binding affinity. In addition, molecular dynamics simulation is considered for the complexes of best two modified drugs with their receptors. Modifications by the halogens show significant changes in the charge distribution, dipole moment, thermodynamic stability, enthalpy and free energy. The highest binding affinity value (-36.401 KJ/mol) was obtained for M14. Hybrid quantum mechanics/molecular mechanics calculation shows a binding energy of -255.140 KJ/mol. Modified drugs have significant hydrogen and non-covalent bonding interactions with amino acids of the receptor. Molecular dynamics simulation disclosed that the root-mean-square-deviation of the alpha carbon associated with M6-1KMV and M14-1KMV complexes is 2.367 Å and 2.622 Å, respectively. Moreover, the interactions between modified drugs and receptor are mostly persevered in 25 nanosecond molecular dynamics simulation. Ensemble-based docking also confirmed that modified drugs show strong non-bonding interactions with different crystallographic and molecular dynamics based conformers. The best scored drugs show considerable pharmacokinetic properties. Modified derivatives M5, M6, M8, M10, M13 and M14 show the better binding affinity and a good number of hydrogen and other non-bonding interactions with the target protein which are similar to other anticancer drugs.Communicated by Ramaswamy H. Sarma.

Recent studies of docking and molecular dynamics simulation for liquid-phase enantioseparations

Liquid-phase enantioseparations have been fruitfully applied in several fields of science. Various applications along with technical and theoretical advancements contributed to increase significantly the knowledge in this area. Nowadays, chromatographic techniques, in particular HPLC on chiral stationary phase, are considered as mature technologies. In the last thirty years, CE has been also recognized as one of the most versatile technique for analytical scale separation of enantiomers. Despite the huge number of papers published in these fields, understanding mechanistic details of the stereoselective interaction between selector and selectand is still an open issue, in particular for high-molecular weight chiral selectors like polysaccharide derivatives. With the ever growing improvement of computer facilities, hardware and software, computational techniques have become a basic tool in enantioseparation science. In this field, molecular docking and dynamics simulations proved to be extremely adaptable to model and visualize at molecular level the spatial proximity of interacting molecules in order to predict retention, selectivity, enantiomer elution order, and profile noncovalent interaction patterns underlying the recognition process. On this basis, topics and trends in using docking and molecular dynamics as theoretical complement of experimental LC and CE chiral separations are described herein. The basic concepts of these computational strategies and seminal studies performed over time are presented, with a specific focus on literature published between 2015 and November 2018. A systematic compilation of all published literature has not been attempted.

Synthesis and Evaluation of the 4-Substituted 2-Hydroxy-5-Iodochalcones and Their 7-Substituted 6-Iodoflavonol Derivatives for Inhibitory Effect on Cholinesterases and β-Secretase

A series of 2-aryl-3-hydroxy-6-iodo-4H-chromen-4-ones substituted at the 7-position with a halogen atom (X = F, Cl and Br) or methoxy group and their corresponding 4-substituted 2-hydroxy-5-iodochalcone precursors were evaluated in vitro for inhibitory effect against acetylcholinesterase (AChE), butyrylcholinesterase (BChE) and β-secretase (BACE1) activities. Although moderate inhibitory effect was observed for the chalcones against AChE, derivatives 2h, 2j and 2n exhibited significant inhibitory effect against BChE and BACE-1. The 2-aryl-7-fluoro-8-iodoflavonols 3b and 3c, on the other hand, exhibited increased activity and selectivity against AChE and reduced effect on BACE-1. The flavonols 3h, 3i, 3k, 3l and 3p exhibited moderate inhibitory effect against AChE, but significant inhibition against BChE. Compounds 2j and 3l exhibited non-competitive mode of inhibition against BACE-1. Molecular docking predicted strong interactions with the protein residues in the active site of BACE-1 implying these compounds bind with the substrate. Similarly docking studies predicted interaction of the most active compounds with both CAS and PAS of either AChE or BChE with mixed type of enzyme inhibition confirmed by kinetic studies.

Biomolecular Simulations of Halogen Bonds with a GROMOS Force Field

Halogen bonds (XBs) are non-covalent interactions in which halogens (X), acting as electrophiles, interact with Lewis bases. XBs are able to mediate protein-ligand recognition and therefore play an important role in rational drug design. In this context, the development of molecular modeling tools that can tackle XBs is paramount. XBs are predominantly explained by the existence of a positive region on the electrostatic potential of X named the σ-hole. Typically, with molecular mechanics force fields, this region is modeled using a charged extra point (EP) linked to X along the R-X covalent bond axis. In this work, we developed the first EP-based strategy for GROMOS force fields (specifically GROMOS 54A7) using bacteriophage T4 lysozyme in complex with both iodobenzene and iodopentafluorobenzene as a prototype system. Several EP parametrization schemes were tested by adding a virtual interaction site to ligand topologies retrieved from the Automated Topology Builder (ATB) and Repository. Contrary to previous approaches using other force fields, our analysis is based on the capability of each parametrization scheme to sample XBs during MD simulations. Our results indicate that the implementation of an EP at a distance from iodine corresponding to R_min provides a good qualitative description of XBs in MD simulations, supporting the compatibility of our approach with the GROMOS 54A7 force field.

Halogen bond in high-performance liquid chromatography enantioseparations: Description, features and modelling

Halogen bond (XB)-driven enantioseparations involve halogen-centred regions of electronic charge depletion (σ-hole) as electrophilic recognition sites. The knowledge in this field is still in its infancy. Indeed, although the influence of halogens on enantioseparation have been often considered, only recently the function of electrophilic halogens (Cl, Br, I) as enantioseparations 'drivers' has been demonstrated by our groups. Further to these studies, in this paper we focus on some unexplored issues. First, as XB-driven chiral recognition mechanisms are at an early stage of comprehension, a theoretical investigation based on a series of 32 molecular dynamic (MD) simulations was performed by using polyhalogenated 4,4'-bipyridines and polysaccharide-based polymers as ligands and receptors, respectively. Enantiomer elution orders (EEOs) were derived from calculations and the theoretical model accounted for some analyte- and chiral stationary phase (CSP)-dependent experimental EEO inversions. Then, the function of halogen-centred σ-holes in competitive systems, presenting also hydrogen bond (HB) centres as recognition sites, was considered. In this regard, Pirkle's enantioseparations of halogenated compounds performed on Whelk-O1 were theoretically re-examined and electrostatic potentials (EPs) associated with both σ-holes on halogens and HB centres were computed and compared. Then, the enantioseparation of halogenated 2-nitro-1-arylethanols was performed on cellulose tris(3,5-dimethylphenylcarbamate) (CDMPC) and the influence of halogen substituents on the chromatographic results was evaluated by correlating theoretical and experimental data.

Polysaccharide-based chiral stationary phases as halogen bond acceptors: A novel strategy for detection of stereoselective σ-hole bonds in solution

In the last few years, halogen bonds have been exploited in a variety of research areas both in the solid state and in solution. Nevertheless, several factors make formation and detection of halogen bonds in solution challenging. Moreover, to date, few chiral molecules containing electrophilic halogens as recognition sites have been reported. Recently, we described the first series of halogen-bond-driven enantioseparations performed on cellulose tris(3,5-dimethylphenylcarbamate) by high-performance liquid chromatography. Herein the performances of amylose tris(3,5-dimethylphenylcarbamate) as halogen bond acceptor were also investigated and compared with respect to cellulose tris(3,5-dimethylphenylcarbamate). With the aim to explore the effect of polysaccharide backbone on the enantioseparations, the thermodynamic parameters governing the halogen-dependent enantioseparations on both cellulose and amylose polymers were determined by a study at variable temperature and compared. Molecular dynamics were performed to model the halogen bond in polysaccharide-analyte complexes. Chiral halogenated 4,4'-bipyridines were used as test compounds (halogen bond donors). On this basis, a practical method for detection of stereoselective halogen bonds in solution was developed, which is based on the unprecedented use of high-performance liquid chromatography as technical tool with polysaccharide polymers as molecular probes (halogen bond acceptors). The analytical strategy showed higher sensitivity for the detection of weak halogen bonds.

Halogen-directed drug design for Alzheimer's disease: a combined density functional and molecular docking study

A series of halogen-directed donepezil drugs has been designed to inhibit acetyl cholinesterase (AChE). Density Functional theory (DFT) has been employed to optimize the chair as well as boat conformers of the parent drug and modified ligands at B3LYP/MidiX and B3LYP/6-311G + (d,p) level of theories. Charge distribution, dipole moment, enthalpy, free energy and molecular orbitals of these ligands are also investigated to understand how the halogen-directed modifications impact the ligand structure and govern the non-bonding interactions with the receptors. Molecular docking calculation has been performed to understand the similarities and differences between the binding modes of unmodified and halogenated chair-formed ligands. Molecular docking indicated donepezil and modified ligands had non-covalent interactions with hydrophobic gorges and anionic subsites of AChE. The -CF3-directed ligand possessed the most negative binding affinity. Non-covalent interactions within the ligand-receptor systems were found to be mostly hydrophobic and π- stacking type. F, Cl and -CF3 containing ligands emerge as effective and selective AChE inhibitors, which can strongly interact with the two active sites of AChE. In addition, we have also investigated selected pharmacokinetic parameters of the parent and modified ligands.

Insights into halogen bond-driven enantioseparations

Although the halogen bond (XB) has been so far mainly studied in silico and in the solid state, its potential impact in solution is yet to be fully understood. In this study, we describe the first systematic investigation on the halogen bond in solvated environment by high-performance liquid chromatography (HPLC). Thirty three atropisomeric polyhalogenated-4,4'-bipyridines (HBipys), containing Cl, Br and I as substituents, were selected and used as potential XB donors (XBDs) on two cellulose-based chiral stationary phases (CSPs) containing potential XB acceptors (XBAs). The impact of the halogens on the enantiodiscrimination mechanism was investigated and iodine showed a pivotal role on the enantioseparation in non-polar medium. Electrostatic potentials (EPs) were computed to understand the electrostatic component of CSP-analyte interaction. Moreover, van't Hoff studies for ten HBipys were performed and the thermodynamic parameters governing the halogen-dependent enantioseparations are discussed. Finally, a molecular dynamic (MD) simulation is proposed to model halogen bond in polysaccharide-analyte complexes by inclusion of a charged extra point to represent the positive 'σ-hole' on the halogen atom. On the basis of both experimental results and theoretical data, we have profiled the halogen bond as a chemo-, regio-, site- and stereoselective interaction which can work in HPLC environment besides other known interactions based on the complementarity between selector and selectand.

AutoDock VinaXB: implementation of XBSF, new empirical halogen bond scoring function, into AutoDock Vina

Background

Halogen bonding has recently come to play as a target for lead optimization in rational drug design. However, most docking program don’t account for halogen bonding in their scoring functions and are not able to utilize this new approach. In this study a new and improved halogen bonding scoring function (XBSF) is presented along with its implementation in the AutoDock Vina molecular docking software. This new improved program is termed as AutoDock VinaXB, where XB stands for the halogen bonding parameters that were added.

Results

XBSF scoring function is derived based on the X···A distance and C–X···A angle of interacting atoms. The distance term was further corrected to account for the polar flattening effect of halogens. A total of 106 protein-halogenated ligand complexes were tested and compared in terms of binding affinity and docking poses using Vina and VinaXB. VinaXB performed superior to Vina in the majority of instances. VinaXB was closer to native pose both above and below 2 Å deviation categories almost twice as frequently as Vina.

Conclusions

Implementation of XBSF into AutoDock Vina has been shown to improve the accuracy of the docking result with regards to halogenated ligands. AutoDock VinaXB addresses the issues of halogen bonds that were previously being scored unfavorably due to repulsion factors, thus effectively lowering the output RMSD values.

Electronic supplementary material

The online version of this article (doi:10.1186/s13321-016-0139-1) contains supplementary material, which is available to authorized users.

Asymmetric bifurcated halogen bonds

Halogen bonding (XB) is being extensively explored for its potential use in advanced materials and drug design. Despite significant progress in describing this interaction by theoretical and experimental methods, the chemical nature remains somewhat elusive, and it seems to vary with the selected system. In this work we present a detailed DFT analysis of three-center asymmetric halogen bond (XB) formed between dihalogen molecules and variously 4-substituted 1,2-dimethoxybenzene. The energy decomposition, orbital, and electron density analyses suggest that the contribution of electrostatic stabilization is comparable with that of non-electrostatic factors. Both terms increase parallel with increasing negative charge of the electron donor molecule in our model systems. Depending on the orientation of the dihalogen molecules, this bifurcated interaction may be classified as 'σ-hole - lone pair' or 'σ-hole - π' halogen bonds. Arrangement of the XB investigated here deviates significantly from a recent IUPAC definition of XB and, in analogy to the hydrogen bonding, the term bifurcated halogen bond (BXB) seems to be appropriate for this type of interaction.

Synthesis of Chlorinated Tetracyclic Compounds and Testing for Their Potential Antidepressant Effect in Mice

The synthesis of the tetracyclic compounds 1-(4,5-dichloro-9,10-dihydro-9,10-ethanoanthracen-11-yl)-N-methylmethanamine (5) and 1-(1,8-dichloro-9,10-dihydro-9,10-ethanoanthracen-11-yl)-N-methylmethanamine (6) as a homologue of the anxiolytic and antidepressant drugs benzoctamine and maprotiline were described. The key intermediate aldehydes (3) and (4) were successfully synthesized via a [4 + 2] cycloaddition between acrolein and 1,8-dichloroanthracene. The synthesized compounds were investigated for antidepressant activity using the forced swimming test. Compounds (5), (6) and (3) showed significant reduction in the mice immobility indicating significant antidepressant effects. These compounds significantly reduced the immobility times at a dose 80 mg/kg by 84.0%, 86.7% and 71.1% respectively.

Antileishmanial, DNA Interaction, and Docking Studies of Some Ferrocene-Based Heteroleptic Pentavalent Antimonials

A series of ferrocenyl pentavalent antimonials (1-8) were synthesized and characterized by elemental analysis, FT-IR, and multinuclear ((1) H and (13) C) NMR spectroscopy. These antimonials were evaluated for their antileishmanial potential against Leishmania tropica KWH23, and by biocompatibility and membrane permeability assays. Moreover, mechanistic studies were carried out, mediated by DNA targeting followed by computational docking of ferrocenyl antimonials against the leishmanial trypanothione reductase enzyme. It was observed that the antimonials 1-8 were 390-fold more efficacious (IC50 ) as compared with the standard antimonial drug used. Cytotoxicity results showed that these antimonials are highly active even at low concentrations and are biocompatible with human macrophages. Antimonials 1-8 exhibited extensive intercalation with DNA and, furthermore, docking interactions highlighted the potential interactive binding of the anitimonials within the trypanothione reductase active site, with van der Waals interactions contributing significantly to the process. Hence, it is suggested that the reported antimonials demonstrate high efficacy, less toxicity, and target multiple sites of the Leishmania parasite.

Molecular dynamics simulation of halogen bonding in Cl2, BrCl, and mixed Cl2/Br2 clathrate hydrates

Clathrate hydrate phases of dihalogen molecules have properties that differ from those of other guest molecules of similar size. The water oxygen–chlorine distances in the structure I (sI) Cl2 hydrate are smaller than the sum of the van der Waals radii of oxygen and chlorine. Bromine hydrate forms a unique clathrate hydrate structure that is not seen in other guest substances. In mixed Cl2/Br2 structure I hydrate, the water oxygen–bromine distances are also smaller than the sum of the oxygen and bromine van der Waals radii. We previously studied the structure of three dihalogen clathrate hydrates using single crystal X-ray diffraction and described these structural features in terms of halogen bonding between the dihalogen and water molecules. In this work, we perform molecular dynamics simulations of cubic sI Cl2, mixed Cl2/Br2, and BrCl clathrate hydrate phases. We perform quantum chemical computations on the dihalogen molecules to determine the nature of σ-hole near the halogen atoms. We fit the electrostatic potential of the molecules to point charge models including dummy atoms that represent σ-holes adjacent to the halogen molecules. Molecular dynamics simulations are used to determine the lattice constants, radial distribution functions, and guest dynamics in these phases. We determine the effect of guest size and difference in halogen bonding on the properties of the clathrate hydrate phase. Simulations for the Cl2, BrCl, and mixed Cl2/Br2 hydrates are performed with small cages of the sI clathrate hydrate phases completely full or filled with experimental occupancies with Cl2 guests.

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.