On Extension of the Current Biomolecular Empirical Force Field for the Description of Halogen Bonds

Interplay of halogen bonding and solvation in protein-ligand binding

Halogen bonding is increasingly utilized in efforts to achieve high affinity and selectivity of molecules designed to bind proteins, making it paramount to understand the relationship between structure, dynamics, and thermodynamic driving forces. We present a detailed analysis addressing this problem using a series of protein-ligand complexes involving single halogen substitutions - F, Cl, Br, and I - and nearly identical structures. Isothermal titration calorimetry reveals an increasingly favorable binding enthalpy from F to I that correlates with the halogen size and σ-hole electropositive character, but is partially counteracted by unfavorable entropy, which is constant from F to Cl and Br, but worse for I. Consequently, the binding free energy is roughly equal for Cl, Br, and I. QM and solvation-free-energy calculations reflect an intricate balance between halogen bonding, hydrogen bonds, and solvation. These advances have the potential to aid future drug design initiatives involving halogenated compounds.

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.

A General Picture of Cucurbit[8]uril Host–Guest Binding: Recalibrating Bonded Interactions

Atomic-level understanding of the dynamic feature of host–guest interactions remains a central challenge in supramolecular chemistry. The remarkable guest binding behavior of the Cucurbiturils family of supramolecular containers makes them promising drug carriers. Among Cucurbit[n]urils, Cucurbit[8]uril (CB8) has an intermediate portal size and cavity volume. It can exploit almost all host–guest recognition motifs formed by this host family. In our previous work, an extensive computational investigation of the binding of seven commonly abused and structurally diverse drugs to the CB8 host was performed, and a general dynamic binding picture of CB8-guest interactions was obtained. Further, two widely used fixed-charge models for drug-like molecules were investigated and compared in great detail, aiming at providing guidelines in choosing an appropriate charge scheme in host-guest modelling. Iterative refitting of atomic charges leads to improved binding thermodynamics and the best root-mean-squared deviation from the experimental reference is 2.6 kcal/mol. In this work, we focus on a thorough evaluation of the remaining parts of classical force fields, i.e., the bonded interactions. The widely used general Amber force fields are assessed and refitted with generalized force-matching to improve the intra-molecular conformational preference, and thus the description of inter-molecular host–guest interactions. The interaction pattern and binding thermodynamics show a significant dependence on the modelling parameters. The refitted system-specific parameter set improves the consistency of the modelling results and the experimental reference significantly. Finally, combining the previous charge-scheme comparison and the current force-field refitting, we provide general guidelines for the theoretical modelling of host–guest binding.

Substituent Effect and Its Halogen-Atom Dependence of Halogen Bonding Viewed through Electron Density Changes

Elucidating how the halogen-bonding ability and strength are controlled by the substituent effect and how this control depends on halogen atom will be essential for finely-tuned design of functionally important molecules. Here, this problem is tackled by analyzing the electron density differences/changes for variously substituted halobenzenes. It is shown that the anisotropy of the electron distribution around the halogen atom, which is an important factor for halogen-bonding ability, is not much affected by the substituent effect and rather simply depends on the halogen atom, while the partial charge on the halogen atom, which is related to the bond dipole of the C-X bond, is significantly modulated by the substituent effect and gives rise to enhancement of the electrostatic potential on the line extended from the C-X bond. The properties related to the polarization effect are also discussed.

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.

Efficient Computation of the Interaction Energies of Very Large Non-covalently Bound Complexes

We present a new benchmark set consisting of 16 large non-covalently bound systems (LNCI16) ranging from 380 up to 1988 atoms and featuring diverse interaction motives. Gas-phase interaction energies are calculated with various composite DFT, semi-empirical quantum mechanical (SQM), and force field (FF) methods and are evaluated using accurate DFT reference values. Of the employed QM methods, PBEh-3c proves to be the most robust for large systems with a relative mean absolute deviation (relMAD) of 8.5% with respect to the reference interaction energies. r2SCAN-3c yields an even smaller relMAD, at least for the subset of complexes for which the calculation could be converged, but is less robust for systems with smaller HOMO–LUMO gaps. The inclusion of Fock-exchange is therefore important for the description of very large non-covalent interaction (NCI) complexes in the gas phase. GFN2-xTB was found to be the best performer of the SQM methods with an excellent result of only 11.1% deviation. From the assessed force fields, GFN-FF and GAFF achieve the best accuracy. Considering their low computational costs, both can be recommended for routine calculations of very large NCI complexes, with GFN-FF being clearly superior in terms of general applicability. Hence, GFN-FF may be routinely applied in supramolecular synthesis planning.

1 Introduction

2 The LNCI16 Benchmark Set

3 Computational Details

4 Generation of Reference Values

5 Results and Discussion

6 Conclusions

Preserving the Integrity of Empirical Force Fields

Generalized force fields (FFs) act as extensions to biomolecular FFs to provide a wide coverage of organic molecules. However, their precise application to an arbitrary molecule presents a separate challenge. We show that MATCH assigns different atom types and bonded and nonbonded parameters than CGenFF, and the AM1-BCC charge model, commonly used with GAFF/GAFF2, does not exactly reproduce the performance of the RESP charge model. The results indicate the need for caution when employing FFs to ensure their integrity with respect to their implementation and validation.

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.

Nature of hydrogen-bond-enhanced halogen bonding viewed through electron density changes

Elucidating the mechanism of how we can achieve fine tuning of intermolecular interaction strength will be helpful for designing functionally important molecules. In the present study, a theoretical analysis is conducted, by examining the electron density changes, for two halogen-bonding iodinated systems whose halogen-bond strengths have been considered to be enhanced by the presence of a hydrogen-bond donating group (termed hydrogen-bond-enhanced halogen bonding). It is shown that, contrary to the expectation obtained from the enhancement of electrostatic potential along the line extended from the C-I bond, the anisotropy of electron distribution on the iodine atom remains nearly the same. This means that the hydrogen bond and halogen bond contribute almost independently and additively to the enhancement of electrostatic potential, indicating the nature of this enhancement and, in a more general sense, the relationship between the strength and the extent of directionality of halogen bonding.

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.

Simulating chalcogen bonding using molecular mechanics: a pseudoatom approach to model ebselen

The organoselenium compound ebselen has recently been investigated as a treatment for COVID-19; however, efforts to model ebselen in silico have been hampered by the lack of an efficient and accurate method to assess its binding to biological macromolecules. We present here a Generalized Amber Force Field modification which incorporates classical parameters for the selenium atom in ebselen, as well as a positively charged pseudoatom to simulate the σ-hole, a quantum mechanical phenomenon that dominates the chemistry of ebselen. Our approach is justified using an energy decomposition analysis of a number of density functional theory–optimized structures, which shows that the σ-hole interaction is primarily electrostatic in origin. Finally, our model is verified by conducting molecular dynamics simulations on a number of simple complexes, as well as the clinically relevant enzyme SOD1 (superoxide dismutase), which is known to bind to ebselen.

Graphical Abstract

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.

Hidden Halogen-Bonding Ability of Fluorine Manifesting in the Hydrogen-Bond Configurations of Hydrogen Fluoride

Elucidating how the intermolecular interactions of a covalently bonded fluorine atom are similar to and different from those of the other halogen atoms will be helpful for a better unified understanding of them. In the present study, the case of hydrogen fluoride is theoretically studied from this viewpoint by using the techniques of electron density analysis, molecular dynamics of liquid, and others. It is shown that the extra-point model, which locates an additional charge site on the line extended from (not within) the covalent bond and has been adopted for halogen-bonding systems as a key to the generation of proper stability and directionality, works well also in this case. A significantly bent hydrogen-bond configuration, which is characteristic of the intermolecular interactions of hydrogen fluoride, is reasonably well reproduced, meaning that it is a manifestation of the latent halogen-bonding ability, which is hidden by the strongly electronegative nature.

Topology Automated Force-Field Interactions (TAFFI): A Framework for Developing Transferable Force Fields

Force-field development has undergone a revolution in the past decade with the proliferation of quantum chemistry based parametrizations and the introduction of machine learning approximations of the atomistic potential energy surface. Nevertheless, transferable force fields with broad coverage of organic chemical space remain necessary for applications in materials and chemical discovery where throughput, consistency, and computational cost are paramount. Here, we introduce a force-field development framework called Topology Automated Force-Field Interactions (TAFFI) for developing transferable force fields of varying complexity against an extensible database of quantum chemistry calculations. TAFFI formalizes the concept of atom typing and makes it the basis for generating systematic training data that maintains a one-to-one correspondence with force-field terms. This feature makes TAFFI arbitrarily extensible to new chemistries while maintaining internal consistency and transferability. As a demonstration of TAFFI, we have developed a fixed-charge force-field, TAFFI-gen, from scratch that includes coverage for common organic functional groups that is comparable to established transferable force fields. The performance of TAFFI-gen was benchmarked against OPLS and GAFF for reproducing several experimental properties of 87 organic liquids. The consistent performance of these force fields, despite their distinct origins, validates the TAFFI framework while also providing evidence of the representability limitations of fixed-charge force fields.

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.

A Reduced Generalized Force Field for Biological Halogen Bonds

The halogen bond (or X-bond) is a noncovalent interaction that is increasingly recognized as an important design tool for engineering protein-ligand interactions and controlling the structures of proteins and nucleic acids. In the past decade, there have been significant efforts to characterize the structure-energy relationships of this interaction in macromolecules. Progress in the computational modeling of X-bonds in biological molecules, however, has lagged behind these experimental studies, with most molecular mechanics/dynamics-based simulation methods not properly treating the properties of the X-bond. We had previously derived a force field for biological X-bonds (ffBXB) based on a set of potential energy functions that describe the anisotropic electrostatic and shape properties of halogens participating in X-bonds. Although fairly accurate for reproducing the energies within biomolecular systems, including X-bonds engineered into a DNA junction, the ffBXB with its seven variable parameters was considered to be too unwieldy for general applications. In the current study, we have generalized the ffBXB by reducing the number of variables to just one for each halogen type and show that this remaining electrostatic variable can be estimated for any new halogenated molecule through a standard restricted electrostatic potential calculation of atomic charges. In addition, we have generalized the ffBXB for both nucleic acids and proteins. As a proof of principle, we have parameterized this reduced and more general ffBXB against the AMBER force field. The resulting parameter set was shown to accurately recapitulate the quantum mechanical landscape and experimental interaction energies of X-bonds incorporated into DNA junction and T4 lysozyme model systems. Thus, this reduced and generalized ffBXB is more readily adaptable for incorporation into classical molecular mechanics/dynamics algorithms, including those commonly used to design inhibitors against therapeutic targets in medicinal chemistry and materials in biomolecular engineering.

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.

A low cost, high accuracy method for halogen bonding complexes

The ONIOM scheme M052X/[Def2TZVP+Def2TZVPD.ECP(I)]:AM1 is shown to represent halogen bond (XB) geometries nearly as well as DFT while being more than two orders of magnitude faster in systems containing >40 atoms. This finding is shown to hold for 40 XB donors, which cover most known backbones, and for a range of neutral and anionic Lewis bases. Complexation free energies can be accurately computed using these geometries and a single-point energy calculation at the DFT level. This approach circumvents the unfavorable scaling of computing time associated with modeling large systems involving halogen bonding.

Dissecting the electric quadrupolar and polarization effects operating in halogen bonding through electron density analysis with a focus on bromine

The form of the electron density change (or difference) is usable as a kind of fingerprint of the electronic structural origin or mechanism that gives rise to intermolecular interactions. Here, this method is applied to halogen-bonding brominated systems to dissect the electric quadrupolar effect (arising from the anisotropic distribution of the valence electrons and intrinsic to the s²p_x ²p_y ²p_z electronic configuration) and the polarization effect (induced by a partial negative charge of the halogen-bond accepting atom). It is shown that a suitable location of the "extra point" for placing a partial positive charge to represent the former is crucial and is clearly found from the electron density difference from the spherically isotropic Br^- ion, while the latter consists of the dipolar polarization of the Br atom and the delocalized polarization of the whole molecule. A practical way for application to molecular dynamics simulations, etc., to represent these two factors is discussed.

Automated parameterization of quantum-mechanically derived force-fields including explicit sigma holes: A pathway to energetic and structural features of halogen bonds in gas and condensed phase

In classical molecular dynamics, general purpose atomistic force-fields (FFs) often deliver inaccurate results when dealing with halogen bonds (XBs), notwithstanding their crucial role in many fields of science, ranging from material design to drug development. Given the large dimensions of the systems of interest, it would be therefore desirable to increase the FF accuracy maintaining the simplicity of the standard Lennard-Jones (LJ) plus point charge description to avoid an excessive computational cost. A simple yet effective strategy consists in introducing a number of virtual sites able to mimic the so-called "explicit σ-hole." In this work, we present an automated FF parameterization strategy based on a global optimization of both LJ and charge parameters with respect to accurate quantum mechanical data, purposely computed for the system under investigation. As a test case, we report on two homologue series, characterized either by weak or strong XBs, namely, the di-halogenated methanes and the mono-, di-, and tri-substituted acetonitriles, taking into consideration Cl, Br, and I substituents. The resulting quantum mechanically derived FFs are validated for each compound in the gas and in the condensed phase by comparing them to general purpose and specific FFs without virtual sites and to highly accurate reference quantum mechanical data. The results strongly support the adoption of the specific FFs with virtual sites, which overcome the other investigated models in representing both gas phase energetics and the structural patterns of the liquid phase structure related to the presence of XBs.

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.

Spectrometric and computational studies of the binding of HIV-1 integrase inhibitors to viral DNA extremities

Three integrase strand transfer inhibitors are in intensive clinical use, raltegravir (RAL), elvitegravir (EVG) and dolutegravir (DTG). The onset of integrase resistance mutations limits their therapeutic efficiency. As put forth earlier, the drug affinity for the intasome could be improved by targeting preferentially the retroviral nucleobases, which are little, if at all, mutation-prone. We report experimental results of anisotropy fluorescence titrations of viral DNA by these three drugs. These show the DTG > EVG > RAL ranking of their inhibitory activities of the intasome to correspond to that of their free energies of binding, ∆Gs, to retroviral DNA, and that such a ranking is only governed by the binding enthalpies, ∆H, the entropy undergoing marginal variations. We sought whether this ranking might be reproduced through quantum chemistry (QC) Density Functional Theory calculations of intermolecular interaction energies between simplified models consisting of sole halobenzene ring and the highly conserved retroviral nucleobases G4 and C16. These calculations showed that binding of EVG has a small preference over DTG, while RAL ranked third. This indicates that additional interactions of the diketoacid parts of the drugs with DNA could be necessary to further enable preferential binding of DTG. The corresponding ∆Etotvalues computed with a polarizable molecular mechanics/dynamics procedure, Sum of Interactions Between Fragments Ab initio computed (SIBFA), showed good correlations with this ∆E(QC) ranking. These validations are an important step toward the use of polarizable molecular dynamics simulations on DTG or EVG derivatives in their complexes with the complete intasome, an application now motivated and enabled by the advent of currently developed and improved massively parallel software.

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.

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.

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.

Amino Acid Hot Spots of Halogen Bonding: A Combined Theoretical and Experimental Case Study of the 5-HT7 Receptor

A computational approach combining a structure-activity relationship library of halogenated and the corresponding unsubstituted ligands (called XSAR) with QM-based molecular docking and binding free energy calculations was used to search for amino acids frequently targeted by halogen bonding (hot spots) in a 5-HT₇R as a case study. The procedure identified two sets of hot spots, extracellular (D2.65, T2.64, and E7.35) and transmembrane (C3.36, T5.39, and S5.42), which were further verified by a synthesized library of halogenated arylsulfonamide derivatives of (aryloxy)ethylpiperidines. It was found that a halogen bond formed between T5.39 and a bromine atom at 3-position of the aryloxy fragment caused the most remarkable, 35-fold increase in binding affinity for 5-HT₇R when compared to the nonhalogenated analog. The proposed paradigm of halogen bonding hot spots was additionally verified on D₄ dopamine receptor showing that it can be used in rational drug design/optimization for any protein target.

Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations

Accurate in silico prediction of protein-ligand binding affinities has been a primary objective of structure-based drug design for decades due to the putative value it would bring to the drug discovery process. However, computational methods have historically failed to deliver value in real-world drug discovery applications due to a variety of scientific, technical, and practical challenges. Recently, a family of approaches commonly referred to as relative binding free energy (RBFE) calculations, which rely on physics-based molecular simulations and statistical mechanics, have shown promise in reliably generating accurate predictions in the context of drug discovery projects. This advance arises from accumulating developments in the underlying scientific methods (decades of research on force fields and sampling algorithms) coupled with vast increases in computational resources (graphics processing units and cloud infrastructures). Mounting evidence from retrospective validation studies, blind challenge predictions, and prospective applications suggests that RBFE simulations can now predict the affinity differences for congeneric ligands with sufficient accuracy and throughput to deliver considerable value in hit-to-lead and lead optimization efforts. Here, we present an overview of current RBFE implementations, highlighting recent advances and remaining challenges, along with examples that emphasize practical considerations for obtaining reliable RBFE results. We focus specifically on relative binding free energies because the calculations are less computationally intensive than absolute binding free energy (ABFE) calculations and map directly onto the hit-to-lead and lead optimization processes, where the prediction of relative binding energies between a reference molecule and new ideas (virtual molecules) can be used to prioritize molecules for synthesis. We describe the critical aspects of running RBFE calculations, from both theoretical and applied perspectives, using a combination of retrospective literature examples and prospective studies from drug discovery projects. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative binding free energy simulations, with a focus on real-world drug discovery applications. We offer guidelines for improving the accuracy of RBFE simulations, especially for challenging cases, and emphasize unresolved issues that could be improved by further research in the field.

Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches

Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.

Optimal affinity ranking for automated virtual screening validated in prospective D3R grand challenges

The goal of virtual screening is to generate a substantially reduced and enriched subset of compounds from a large virtual chemistry space. Critical in these efforts are methods to properly rank the binding affinity of compounds. Prospective evaluations of ranking strategies in the D3R grand challenges show that for targets with deep pockets the best correlations (Spearman ρ ~ 0.5) were obtained by our submissions that docked compounds to the holo-receptors with the most chemically similar ligand. On the other hand, for targets with open pockets using multiple receptor structures is not a good strategy. Instead, docking to a single optimal receptor led to the best correlations (Spearman ρ ~ 0.5), and overall performs better than any other method. Yet, choosing a suboptimal receptor for crossdocking can significantly undermine the affinity rankings. Our submissions that evaluated the free energy of congeneric compounds were also among the best in the community experiment. Error bars of around 1 kcal/mol are still too large to significantly improve the overall rankings. Collectively, our top of the line predictions show that automated virtual screening with rigid receptors perform better than flexible docking and other more complex methods.

Looking Back, Looking Forward at Halogen Bonding in Drug Discovery

Halogen bonding has emerged at the forefront of advances in improving ligand: receptor interactions. In particular the newfound ability of this extant non-covalent-bonding phenomena has revolutionized computational approaches to drug discovery while simultaneously reenergizing synthetic approaches to the field. Here we survey, via examples of classical applications involving halogen atoms in pharmaceutical compounds and their biological hosts, the unique advantages that halogen atoms offer as both Lewis acids and Lewis bases.

Improved Description of Sulfur Charge Anisotropy in OPLS Force Fields: Model Development and Parameterization

The atomic point-charge model used in most molecular mechanics force fields does not represent well the electronic anisotropy that is featured in many directional noncovalent interactions. Sulfur participates in several types of such interactions with its lone pairs and σ-holes. The current study develops a new model, via the addition of off-atom charged sites, for a variety of sulfur compounds in the OPLS-AA and OPLS/CM5 force fields to address the lack of charge anisotropy. Parameter optimization is carried out to reproduce liquid-state properties, torsional and noncovalent energetics from reliable quantum mechanical calculations, and electrostatic potentials. Significant improvements are obtained for computed free energies of hydration, reducing the mean unsigned errors from ca. 1.4 to 0.4-0.7 kcal/mol. Enhanced accuracy in directionality and energetics is also obtained for molecular complexes with sulfur-containing hydrogen and halogen bonds. Moreover, the new model reproduces the unusual conformational preferences of sulfur-containing compounds with 1,4-intramolecular chalcogen bonds. Transferability of the new force field parameters to cysteine and methionine is verified via molecular dynamic simulations of blocked dipeptides. The study demonstrates the effectiveness of using off-atom charge sites to address electronic anisotropy, and provides a parametrization methodology that can be applied to other systems.

Relationships between hydrogen bonds and halogen bonds in biological systems

The recent recognition that halogen bonding (XB) plays important roles in the recognition and assembly of biological molecules has led to new approaches in medicinal chemistry and biomolecular engineering. When designing XBs into strategies for rational drug design or into a biomolecule to affect its structure and function, we must consider the relationship between this interaction and the more ubiquitous hydrogen bond (HB). In this review, we explore these relationships by asking whether and how XBs can replace, compete against or behave independently of HBs in various biological systems. The complex relationships between the two interactions inform us of the challenges we face in fully utilizing XBs to control the affinity and recognition of inhibitors against their therapeutic targets, and to control the structure and function of proteins, nucleic acids and other biomolecular scaffolds.

The evaluation of QM/MM-driven molecular docking combined with MM/GBSA calculations as a halogen-bond scoring strategy

The combination of quantum mechanics/molecular mechanics-driven (QM/MM) molecular docking with binding free-energy calculations was successfully used to reproduce the X-ray geometries of protein–ligand complexes with halogen bonding. The procedure involves quantum-polarized ligand docking (QPLD) to obtain the QM-derived ligand atomic charges in the protein environment at the B3PW91/cc-pVTZ level and the MM/GBSA (generalized-Born/surface area) algorithm to calculate the binding free energies of resultant complexes. The performance was validated using a set of 106 X-ray complexes and compared with the Glide and AutoDock VinaXB scoring functions in terms of RMSD and the reconstruction of halogen-bond geometry (distance and σ-hole angle). The results revealed that docking and scoring using the QPLD–GBSA procedure outperformed the remaining scoring functions in the majority of instances. Additionally, a comparison of the orientation of the top ranked binding poses calculated using the fixed atomic charges of ligands obtained from force-field parameterization and by QM calculations in the protein environment provides strong evidence that the use of QM-derived charges is significant.

Exploiting directional long range secondary forces for regulating electrostatics-dominated noncovalent interactions

It has been well established that long range secondary electrostatic interactions (SEIs) have a significant effect on the stability of supramolecular complexes. However, general rules for exploiting SEIs in the rational design of diverse supramolecular complexes have been difficult to obtain. In this work, we outline a quantum chemical approach for understanding the strength of electrostatic interactions. This approach is seen to provide excellent correlation between the electrostatic force and the binding energy between two partners in hydrogen-bonded complexes, as well as that between two ions in ion-pair complexes. Furthermore, we illustrate how the understanding of the binding allows for the rational design of new complexes where the association constant between the two partners can be increased or decreased, as desired, by several orders of magnitude. Hence, the current work showcases a general, simple and powerful method of understanding and exploiting long range secondary electrostatic interactions.