Tackling Halogenated Species with PBSA: Effect of Emulating the σ-Hole

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.

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.

MM/PB(GB)SA benchmarks on soluble proteins and membrane proteins

Predicting protein-ligand binding free energy rapidly and accurately remains a challenging question in modern drug discovery. Molecular mechanics/Poisson-Boltzmann (Generalized Born) surface area (MM/PB(GB)SA) has emerged as an essential tool for accelerating cost-efficient binding free energy calculation. This study presents benchmarks with three membrane-bound protein systems and six soluble protein systems. Different parameters were sampled for different benchmarks to explore the highest accuracy. These include ligand charges, protein force fields, extra points, GB models, nonpolar optimization methods, internal dielectric constants and membrane dielectric constants. Comparisons of accuracy were made between MM/PB(GB)SA, docking and free energy perturbation (FEP). The results reveal a competitive performance between MM/PB(GB)SA and FEP. In summary, MM/PB(GB)SA is a powerful approach to predict ligand binding free energy rapidly and accurately. Parameters of MM/PB(GB)SA calculations, such as the GB models and membrane dielectric constants, need to be optimized for different systems. This method can be served as a powerful tool for drug design.

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.

Molecular investigation of the dual inhibition mechanism for targeted P53 regulator MDM2/MDMX inhibitors

Inhibitors that competitively bind MDM2/MDMX can block the inhibition of P53 by MDM2/MDMX and restart its tumor-suppressive effect. Molecular studies targeting MDM2/MDMX inhibitors have always been a hot topic in anticancer drug design. Although numerous inhibitors have been designed previously against MDM2/MDMX, their dual inhibition efficacy has not been demonstrated, and few studies assessed the general causes affecting the dual inhibition of MDM2/MDMX by these inhibitors. Here, molecular dynamics simulations and alanine scanning combined with the interaction entropy method were employed to precisely investigate whether 16 inhibitors could dually inhibit MDM2/MDMX and the similarities and differences in the interaction modes. Thereby addressing the key residue sites affecting dual inhibition. Residues L54/M53, I61/60, M62/61, Y67/66, and V93/92 of MDM2/MDMX, which are in corresponding positions in both protein structures, provide significant conditions for these inhibitors to bind to MDM2/MDMX tightly. In addition, most of these inhibitors prefer to bind MDM2 than MDMX, and residues H96 and I99 in MDM2 are attractive targets for inhibitors, resulting in inhibitors binding to MDM2/MDMX with different affinity. These key residues should be considered in the development of dual inhibitors. For these 16 inhibitors, most have dual inhibitory potential for MDM2/MDMX based on the binding affinity of the complexes. Still, it is questionable whether they can exert excellent dual inhibition considering the assessment of the hot-spots. At least their binding affinity for MDMX is not superior to that for MDM2 due to the difference in energy of the van der Waals interactions at the key sites. Furthermore, based on the analysis of three representative inhibitors (TUZ/HRH and HRQ with different binding preferences for MDM2/MDMX), 3-chloropyridine in TUZ leads to the differential binding affinity between the inhibitor and MDM2/MDMX. It readily forms hydrophobic interactions with the surrounding residues H96 and I99. But this phenomenon does not occur in the TUZ-MDMX system, implying the critical role of residues H96/P95 and I99/L98. And the completely different binding mechanism of HRQ binding to MDM2/MDMX explains its inability to inhibit MDM2 well. Thus, we are cautious about its dual inhibitory ability. Besides, HRH is more prone to strong van der Waals interactions with MDM2 than MDMX whereas its 2-chlorofluorobenzene is detrimental to this. We hope that these findings will provide reliable molecular insights for the screening and optimization of targeting MDM2/MDMX dual inhibitors.

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.

Transmembrane Anion Transport Mediated by Halogen Bonds: Using Off-Center Charges

Synthetic anion transporters are promising therapeutic agents designed to emulate the specialized role of certain transmembrane proteins that maintain the ion concentration in cells. In the last few years, besides hydrogen bonds and ion pairs, halogen bonds have also been explored to promote the association between the synthetic molecule and the anion and their subsequent transport. This interaction is due to an anisotropic charge distribution on the halogen, and therefore, modeling halogen bonds is not a trivial task using classical force field methods that typically rely on point-charge models.Herein, a computational protocol capable of dealing with halogen bonds is presented. This protocol takes advantage of the addition of an off-center particle during the charge fitting procedure, and the resulting set of charges can be used along with the classical force field parameters from GAFF or GROMOS 54A7.

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.

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.

Halogen bond in separation science: A critical analysis across experimental and theoretical results

The halogen bond (XB) is a noncovalent interaction involving a halogen acting as electrophile and a Lewis base. In the last decades XB has found practical application in several fields. Nevertheless, despite the pivotal role of noncovalent interactions in separation science, investigations of XB in this field are still in their infancy, and so far a limited number of studies focusing on solid phase extraction, liquid-liquid microextraction, liquid-phase chromatography, and gas chromatography separation have been published. In addition, in the last few years, our groups have been systematically studying the potentiality of XB for HPLC enantioseparations. On this basis, in the present paper up-to-date results emerging from focused experiments and theoretical analyses performed by our laboratories are integrated with a descriptive presentation of XB features and the few studies published until now in separation science. Then, the aim of this article is to provide a comprehensive and critical discussion of the topic, and account for some still open issues in the application of XB to separate chemical mixtures.