Targeting Histidine Side Chains in Molecular Design through Nitrogen–Halogen Bonds

Synthesis of Amino Acids Bearing Halodifluoromethyl Moieties and Their Application to p53-Derived Peptides Binding to Mdm2/Mdm4

Introduction

Therapeutic peptides are a significant class of drugs in the treatment of a wide range of diseases. To enhance their properties, such as stability or binding affinity, they are usually chemically modified. This includes, among other techniques, cyclization of the peptide chain by bridging, modifications to the backbone, and incorporation of unnatural amino acids. One approach previously established, is the use of halogenated aromatic amino acids. In principle, they are thereby enabled to form halogen bonds (XB). In this study, we focus on the -R-CF₂X moiety (R = O, NHCO; X = Cl, Br) as an uncommon halogen bond donor. These groups enable more spatial variability in protein-protein interactions. The chosen approach via Fmoc-protected building blocks allows for the incorporation of these modified amino acids in peptides using solid-phase peptide synthesis.

Results and Discussion

Using a competitive fluorescence polarization assay to monitor binding to Mdm4, we demonstrate that a p53-derived peptide with Lys24Nle(εNHCOCF₂X) exhibits an improved inhibition constant K_i compared to the unmodified peptide. Decreasing K_i values observed with the increasing XB capacity of the halogen atoms (F ≪ Cl < Br) indicates the formation of a halogen bond. By reducing the side chain length of Nle(εNHCOCF₂X) to Abu(γNHCOCF₂X) as control experiments and through quantum mechanical calculations, we suggest that the observed affinity enhancement is related to halogen bond-induced intramolecular stabilization of the α-helical binding mode of the peptide or a direct interaction with His54 in human Mdm4.

Revisiting a challenging p53 binding site: a diversity-optimized HEFLib reveals diverse binding modes in T-p53C-Y220C

The cellular tumor antigen p53 is a key component in cell cycle control. The mutation Y220C heavily destabilizes the protein thermally but yields a druggable crevice. We have screened the diversity-optimized halogen-enriched fragment library against T-p53C-Y220C with STD-NMR and DSF to identify hits, which we validated by 1H,15N-HSQC NMR. We could identify four hits binding in the Y220C cleft, one hit binding covalently and four hits binding to an uncharacterized binding site. Compound 1151 could be crystallized showing a flip of C220 and thus opening subsite 3. Additionally, 4482 was identified to alkylate cysteines. Data shows that the diversity-optimized HEFLib leads to multiple diverse hits. The identified scaffolds can be used to further optimize interactions with T-p53C-Y220C and increase thermal stability.

A Theoretical Study of the Halogen Bond between Heteronuclear Halogen and Benzene

Halogen bonds play an important role in many fields, such as biological systems, drug design and crystal engineering. In this work, the structural characteristics of the halogen bond between heteronuclear halogen XD (ClF, BrCl, IBr, ICl, BrF and IF) and benzene were studied using density functional theory. The structures of the complexes between heteronuclear halogen and benzene have Cs symmetry. The interaction energies of the complexes between heteronuclear halogen XD (ClF, BrCl, IBr, ICl, BrF and IF) and benzene range from -27.80 to -37.18 kJ/mol, increasing with the increases in the polarity between the atoms of X and D, and are proportional to the angles of a between the Z axis and the covalent bond of heteronuclear halogen. The electron density (ρ) and corresponding Laplacian (∇²ρ) values indicate that the interaction of the heteronuclear halogen and benzene is a typical long-range weak interaction similar to a hydrogen bond. Independent gradient model analysis suggests that the van der Waals is the main interaction between the complexes of heteronuclear halogen and benzene. Symmetry-adapted perturbation theory analysis suggests that the electrostatic interaction is the dominant part in the complexes of C₆H₆⋯ClF, C₆H₆⋯ICl, C₆H₆⋯BrF and C₆H₆⋯IF, and the dispersion interaction is the main part in the complexes of C₆H₆⋯BrCl, C₆H₆⋯IBr.

Quadruple Target Evaluation of Diversity-Optimized Halogen-Enriched Fragments (HEFLibs) Reveals Substantial Ligand Efficiency for AP2-Associated Protein Kinase 1 (AAK1)

Fragment-based drug discovery is one of the most utilized approaches for the identification of novel weakly binding ligands, by efficiently covering a wide chemical space with rather few compounds and by allowing more diverse binding modes to be found. This approach has led to various clinical candidates and approved drugs. Halogen bonding, on the other hand, has gained traction in molecular design and lead optimization, but could offer additional benefits in early drug discovery. Screening halogen-enriched fragments (HEFLibs) could alleviate problems associated with the late introduction of such a highly geometry dependent interaction. Usually, the binding mode is then already dominated by other strong interactions. Due to the fewer competing interactions in fragments, the halogen bond should more often act as an anchor point for the binding mode. Previously, we proposed a fragment library with a focus on diverse binding modes that involve halogens for gaining initial affinity and selectivity. Herein, we demonstrate the applicability of these HEFLibs with a small set of diverse enzymes: the histone-lysine N-methyltransferase DOT1L, the indoleamine 2,3-dioxygenase 1 (IDO1), the AP2-associated protein kinase 1 (AAK1), and the calcium/calmodulin-dependent protein kinase type 1G (CAMK1G). We were able to identify various binding fragments via STD-NMR. Using ITC to verify these initial hits, we determined affinities for many of these fragments. The best binding fragments exhibit affinities in the one-digit micromolar range and ligand efficiencies up to 0.83 for AAK1. A small set of analogs was used to study structure-affinity relationships and hereby analyze the specific importance of each polar interaction. This data clearly suggests that the halogen bond is the most important interaction of fragment 9595 with AAK1.

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.

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.

Exploring the Role of N6-Substituents in Potent Dual Acting 5'-C-Ethyltetrazolyladenosine Derivatives: Synthesis, Binding, Functional Assays, and Antinociceptive Effects in Mice ∇

Structural determinants of affinity of N⁶-substituted-5'-C-(ethyltetrazol-2-yl)adenosine and 2-chloroadenosine derivatives at adenosine receptor (AR) subtypes were studied with binding and molecular modeling. Small N⁶-cycloalkyl and 3-halobenzyl groups furnished potent dual acting A₁AR agonists and A₃AR antagonists. 4 was the most potent dual acting human (h) A₁AR agonist (K_i = 0.45 nM) and A₃AR antagonist (K_i = 0.31 nM) and highly selective versus A_2A; 11 and 26 were most potent at both h and rat (r) A₃AR. All N⁶-substituted-5'-C-(ethyltetrazol-2-yl)adenosine derivatives proved to be antagonists at hA₃AR but agonists at the rA₃AR. Analgesia of 11, 22, and 26 was evaluated in the mouse formalin test (A₃AR antagonist blocked and A₃AR agonist strongly potentiated). N⁶-Methyl-5'-C-(ethyltetrazol-2-yl)adenosine (22) was most potent, inhibiting both phases, as observed combining A₁AR and A₃AR agonists. This study demonstrated for the first time the advantages of a single molecule activating two AR pathways both leading to benefit in this acute pain model.

Targeting the protein backbone with aryl halides: systematic comparison of halogen bonding and π⋯π interactions using N-methylacetamide

Systematic plane scans reveal a seamless transition from σ-hole interactions with the carbonyl oxygen to interactions with the amide π-electrons at increasing distances.

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

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

Experimental and Theoretical Evaluation of the Ethynyl Moiety as a Halogen Bioisostere

Bioisosteric replacements are widely used in medicinal chemistry to improve physicochemical and ADME properties of molecules while retaining or improving affinity. Here, using the p53 cancer mutant Y220C as a test case, we investigate both computationally and experimentally whether an ethynyl moiety is a suitable bioisostere to replace iodine in ligands that form halogen bonds with the protein backbone. This bioisosteric transformation is synthetically feasible via Sonogashira cross-coupling. In our test case of a particularly strong halogen bond, replacement of the iodine with an ethynyl group resulted in a 13-fold affinity loss. High-resolution crystal structures of the two analogues in complex with the p53-Y220C mutant enabled us to correlate the different affinities with particular features of the binding site and subtle changes in ligand binding mode. In addition, using QM calculations and analyzing the PDB, we provide general guidelines for identifying cases where such a transformation is likely to improve ligand recognition.

Evaluating the potential of halogen bonding in molecular design: automated scaffold decoration using the new scoring function XBScore

We present a QM-derived empirical scoring function for the interaction between aromatic halogenated ligands and the carbonyl oxygen atom of the protein backbone. Applying this scoring function, we developed an algorithm that evaluates the potential of protein-bound ligands to form favorable halogen-bonding contacts upon scaffold decoration with chlorine, bromine, or iodine. Full recovery of all existing halogen bonds in the PDB involving the protein backbone was achieved with our protocol. Interestingly, the potential for introducing halogen bonds through scaffold decoration of unsubstituted aromatic carbon atoms appears to easily match the number of previously known halogen bonds. Our approach can thus be used as a blueprint for integration of halogen bonding into general empirical scoring functions, which at present ignore this interaction. Most importantly, we were able to identify a substantial number of protein-ligand complexes where the benefits and challenges of introducing a halogen bond by molecular design can be studied experimentally.

Thermodynamic and structural characterization of halogen bonding in protein-ligand interactions: a case study of PDE5 and its inhibitors

The significance of halogen bonding in protein-ligand interactions has been recognized recently. We present here the first comprehensive thermodynamic and structural characterization of halogen bonding in PDE5-inhibitor interactions. ITC studies reveal that binding strength of the halogen bonding between chlorine, bromine, and iodine of inhibitor and the protein is -1.57, -3.09, and -5.59 kJ/mol, respectively. The halogens interact with the designed residue Y612 and an unexpected buried water molecule.

Halogen-enriched fragment libraries as chemical probes for harnessing halogen bonding in fragment-based lead discovery

Halogen bonding has recently experienced a renaissance, gaining increased recognition as a useful molecular interaction in the life sciences. Halogen bonds are favorable, fairly directional interactions between an electropositive region on the halogen (the σ-hole) and a number of different nucleophilic interaction partners. Some aspects of halogen bonding are not yet understood well enough to take full advantage of its potential in drug discovery. We describe and present the concept of halogen-enriched fragment libraries. These libraries consist of unique chemical probes, facilitating the identification of favorable halogen bonds by sharing the advantages of classical fragment-based screening. Besides providing insights into the nature and applicability of halogen bonding, halogen-enriched fragment libraries provide smart starting points for hit-to-lead evolution.