Targeting the Gatekeeper MET146 of C-Jun N-Terminal Kinase 3 Induces a Bivalent Halogen/Chalcogen Bond

Evaluation of a Covalent Library of Diverse Warheads (CovLib) Binding to JNK3, USP7, or p53

Purpose

Over the last few years, covalent fragment-based drug discovery has gained significant importance. Thus, striving for more warhead diversity, we conceived a library consisting of 20 covalently reacting compounds. Our covalent fragment library (CovLib) contains four different warhead classes, including five α-cyanoacacrylamides/acrylates (CA), three epoxides (EO), four vinyl sulfones (VS), and eight electron-deficient heteroarenes with a leaving group (SNAr/SN).

Methods

After predicting the theoretical solubility of the fragments by LogP and LogS during the selection process, we determined their experimental solubility using a turbidimetric solubility assay. The reactivities of the different compounds were measured in a high-throughput 5,5'-dithiobis-(2-nitrobenzoic acid) DTNB assay, followed by a (glutathione) GSH stability assay. We employed the CovLib in a (differential scanning fluorimetry) DSF-based screening against different targets: c-Jun N-terminal kinase 3 (JNK3), ubiquitin-specific protease 7 (USP7), and the tumor suppressor p53. Finally, the covalent binding was confirmed by intact protein mass spectrometry (MS).

Results

In general, the purchased fragments turned out to be sufficiently soluble. Additionally, they covered a broad spectrum of reactivity. All investigated α-cyanoacrylamides/acrylates and all structurally confirmed epoxides turned out to be less reactive compounds, possibly due to steric hindrance and reversibility (for α-cyanoacrylamides/acrylates). The SNAr and vinyl sulfone fragments are either highly reactive or stable. DSF measurements with the different targets JNK3, USP7, and p53 identified reactive fragment hits causing a shift in the melting temperatures of the proteins. MS confirmed the covalent binding mode of all these fragments to USP7 and p53, while additionally identifying the SNAr-type electrophile SN002 as a mildly reactive covalent hit for p53.

Conclusion

The screening and target evaluation of the CovLib revealed first interesting hits. The highly cysteine-reactive fragments VS004, SN001, SN006, and SN007 covalently modify several target proteins and showed distinct shifts in the melting temperatures up to +5.1 °C and -9.1 °C.

Iodinated 4,4'-Bipyridines with Antiproliferative Activity Against Melanoma Cell Lines

In the last decade, biological processes involving halogen bond (HaB) as a leading interaction attracted great interest. However, although bound iodine atoms are considered powerful HaB donors, few iodinated new drugs were reported so far. Recently, iodinated 4,4'-bipyridines showed interesting properties as HaB donors in solution and in the solid state. In this paper, a study on the inhibition activity of seven halogenated 4,4'-bipyridines against malignant melanoma (MM) cell proliferation is described. Explorative dose/response proliferation assays were first performed with three 4,4'-bipyridines by using four MM cell lines and the normal BJ fibroblast cell line as control. Among them, the A375 MM cell line was the most sensitive, as determined by MTT assays, which was selected to evaluate the antiproliferative activity of all 4,4'-bipyridines. Significantly, the presence of an electrophilic iodine impacted the biological activity of the corresponding compounds. The 3,3',5,5'-tetrachloro-2-iodo-4,4'-bipyridine showed significant antiproliferation activity against the A375 cell line, and lower toxicity on BJ fibroblasts. Through in silico studies, the stereoelectronic features of possible sites determining the bioactivity were explored. These results pave the way for the utilization of iodinated 4,4'-bipyridines as templates to design new promising HaB-enabled inhibitors of MM cell proliferation.

Recent developments in molecular modeling tools and applications related to pharmaceutical and biomedical research

In modern pharmaceutical and biomedical research, molecular modeling represents a useful tool to explore processes and their mechanistic bases at the molecular level. Integrating experimental and virtual analysis is a fruitful approach to study ligand-receptor interaction in chemical, biochemical and biological environments. In these fields, molecular docking and molecular dynamics are considered privileged techniques for modeling (bio)macromolecules and related complexes. This review aims to present the current landscape of molecular modeling in pharmaceutical and biomedical research by examining selected representative applications published in the last years and highlighting current topics and trends of this field. Thus, a systematic compilation of all published literature has not been attempted herein. After a brief overview of the main theoretical and computational tools used to investigate mechanisms at molecular level, recent applications of molecular modeling in drug discovery, ligand binding and for studying protein conformation and function will be discussed. Furthermore, specific sections will be devoted to the application of molecular modeling for unravelling enantioselective mechanisms underlying the enantioseparation of chiral compounds of pharmaceutical and biomedical interest as well as for studying new forms of noncovalent interactivity identified in biochemical and biological environments. The general aim of this review is to provide the reader with a modern overview of the topic, highlighting advancements and outlooks as well as drawbacks and pitfalls still affecting the applicability of theoretical and computational methods in the field of pharmaceutical and biomedical research.

Characterization of 2,4-Dianilinopyrimidines Against Five P. falciparum Kinases PfARK1, PfARK3, PfNEK3, PfPK9, and PfPKB

Plasmodium kinases are increasingly recognized as potential novel antiplasmodial targets for the treatment of malaria, but only a small subset of these kinases have had structure-activity relationship (SAR) campaigns reported. Herein we report the discovery of CZC-54252 (1) as an inhibitor of five P. falciparum kinases PfARK1, PfARK3, PfNEK3, PfPK9, and PfPKB. 39 analogues were evaluated against all five kinases to establish SAR at three regions of the kinase active site. Nanomolar inhibitors of each kinase were discovered. We identified common and divergent SAR trends across all five kinases, highlighting substituents in each region that improve potency and selectivity for each kinase. Potent analogues were evaluated against the P. falciparum blood stage. Eight submicromolar inhibitors were discovered, of which 37 demonstrated potent antiplasmodial activity (EC50 = 0.16 μM). Our results provide an understanding of features needed to inhibit each individual kinase and lay groundwork for future optimization efforts toward novel antimalarials.

How many kinases are druggable? A review of our current understanding

There are over 500 human kinases ranging from very well-studied to almost completely ignored. Kinases are tractable and implicated in many diseases, making them ideal targets for medicinal chemistry campaigns, but is it possible to discover a drug for each individual kinase? For every human kinase, we gathered data on their citation count, availability of chemical probes, approved and investigational drugs, PDB structures, and biochemical and cellular assays. Analysis of these factors highlights which kinase groups have a wealth of information available, and which groups still have room for progress. The data suggest a disproportionate focus on the more well characterized kinases while much of the kinome remains comparatively understudied. It is noteworthy that tool compounds for understudied kinases have already been developed, and there is still untapped potential for further development in this chemical space. Finally, this review discusses many of the different strategies employed to generate selectivity between kinases. Given the large volume of information available and the progress made over the past 20 years when it comes to drugging kinases, we believe it is possible to develop a tool compound for every human kinase. We hope this review will prove to be both a useful resource as well as inspire the discovery of a tool for every kinase.

Principles and Applications of CF2X Moieties as Unconventional Halogen Bond Donors in Medicinal Chemistry, Chemical Biology, and Drug Discovery

As an orthogonal principle to the established (hetero)aryl halides, we herein highlight the usefulness of CF2X (X = Cl, Br, or I) moieties. Using tool compounds bearing CF2X moieties, we study their chemical/metabolic stability and their logP/solubility, as well as the role of XB in their small molecular crystal structures. Employing QM techniques, we analyze the observed interactions, provide insights into the conformational flexibilities and preferences in the potential interaction space. For their application in molecular design, we characterize their XB donor capacities and its interaction strength dependent on geometric parameters. Implementation of CF2X acetamides into our HEFLibs and biophysical evaluation (STD-NMR/ITC), followed by X-ray analysis, reveals a highly interesting binding mode for fragment 23 in JNK3, featuring an XB of CF2Br toward the P-loop, as well as chalcogen bonds. We suggest that underexplored chemical space combined with unconventional binding modes provides excellent opportunities for patentable chemotypes for therapeutic intervention.

Synthesis and Conformational Analysis of N-Aromatic Acetamides Bearing Thiophene: Effect of Intramolecular Chalcogen-Chalcogen Interaction on Amide Conformational Stability

The conformations of aromatic amides bearing an N-(2-thienyl) or N-(3-thienyl) group were investigated in solution and in the crystal state. NMR spectral data indicate that the conformational preferences of these amides in solution are dependent not only on the relative π-electron densities of the N-aromatic moieties, but also on the three-dimensional relationship between carbonyl oxygen and the N-aromatic moieties. A comparison of the conformational preferences of N-(2-thienyl)amides and N-(3-thienyl)amides revealed that the Z-conformers of N-(2-thienyl)acetamides are stabilized by 1,5-type intramolecular S···O═C interactions between amide carbonyl and thiophene sulfur. The crystal structures of these compounds were similar to the solution structures. The stabilization energy due to 1,5-type intramolecular S···O═C interaction in N-aryl-N-(2-thienyl)acetamides and N-methyl-N-(2-thienyl)acetamide was estimated to be ca. 0.74 and 0.93 kcal/mol, respectively.

Protecting-Group-Free Synthesis of Meridianin A-G and Derivatives and Its Antibiofilm Evaluation

Herein, a protecting-group-free protocol was developed to realize a time and step economy diversification of the Meridianin alkaloid. A broad range of substituents are tolerated to deliver the products in moderate to high yields, and the first synthesis of Meridianin B was achieved. The simplicity of this protocol enables the rapid construction of a Meridianin derivative library for antibiofilm evaluation. Preliminary results reveal that Meridianin derivatives were capable of inhibiting the Acinetobacter baumannii biofilm and lowering the antibiotic MIC synergistically.

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.

Identification of Novel 2,4,5-Trisubstituted Pyrimidines as Potent Dual Inhibitors of Plasmodial PfGSK3/PfPK6 with Activity against Blood Stage Parasites In Vitro

Essential plasmodial kinases PfGSK3 and PfPK6 are considered novel drug targets to combat rising resistance to traditional antimalarial therapy. Herein, we report the discovery of IKK16 as a dual PfGSK3/PfPK6 inhibitor active against blood stage Pf3D7 parasites. To establish structure–activity relationships for PfPK6 and PfGSK3, 52 analogues were synthesized and assessed for the inhibition of PfGSK3 and PfPK6, with potent inhibitors further assessed for activity against blood and liver stage parasites. This culminated in the discovery of dual PfGSK3/PfPK6 inhibitors 23d (PfGSK3/PfPK6 IC₅₀ = 172/11 nM) and 23e (PfGSK3/PfPK6 IC₅₀ = 97/8 nM) with antiplasmodial activity (23dPf3D7 EC₅₀ = 552 ± 37 nM and 23ePf3D7 EC₅₀ = 1400 ± 13 nM). However, both compounds exhibited significant promiscuity when tested in a panel of human kinase targets. Our results demonstrate that dual PfPK6/PfGSK3 inhibitors with antiplasmodial activity can be identified and can set the stage for further optimization efforts.

Potent GCN2 Inhibitor Capable of Reversing MDSC-Driven T Cell Suppression Demonstrates In Vivo Efficacy as a Single Agent and in Combination with Anti-Angiogenesis Therapy

General control nonderepressible 2 (GCN2) protein kinase is a cellular stress sensor within the tumor microenvironment (TME), whose signaling cascade has been proposed to contribute to immune escape in tumors. Herein, we report the discovery of cell-potent GCN2 inhibitors with excellent selectivity against its closely related Integrated Stress Response (ISR) family members heme-regulated inhibitor kinase (HRI), protein kinase R (PKR), and (PKR)-like endoplasmic reticulum kinase (PERK), as well as good kinome-wide selectivity and favorable PK. In mice, compound 39 engages GCN2 at levels ≥80% with an oral dose of 15 mg/kg BID. We also demonstrate the ability of compound 39 to alleviate MDSC-related T cell suppression and restore T cell proliferation, similar to the effect seen in MDSCs from GCN2 knockout mice. In the LL2 syngeneic mouse model, compound 39 demonstrates significant tumor growth inhibition (TGI) as a single agent. Furthermore, TGI mediated by anti-VEGFR was enhanced by treatment with compound 39 demonstrating the complementarity of these two mechanisms.

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.

BINANA 2: Characterizing Receptor/Ligand Interactions in Python and JavaScript

BINding ANAlyzer (BINANA) is an algorithm for identifying and characterizing receptor/ligand interactions and other factors that contribute to binding. We recently updated BINANA to make the algorithm more accessible to a broader audience. We have also ported the Python3 codebase to JavaScript, thus enabling BINANA analysis in the web browser. As proof of principle, we created a web-browser application so students and chemical-biology researchers can quickly visualize receptor/ligand complexes and their unique binding interactions.

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.

Inhibitors of the Hippo Pathway Kinases STK3/MST2 and STK4/MST1 Have Utility for the Treatment of Acute Myeloid Leukemia

Serine/threonine-protein kinases 3 and 4 (STK3 and STK4, respectively) are key components of the Hippo signaling pathway, which regulates cell proliferation and death and provides a potential therapeutic target for acute myeloid leukemia (AML). Herein, we report the structure-based design of a series of pyrrolopyrimidine derivatives as STK3 and STK4 inhibitors. In an initial screen, the compounds exhibited low nanomolar potency against both STK3 and STK4. Crystallization of compound 6 with STK4 revealed two-point hinge binding in the ATP-binding pocket. Further characterization and analysis demonstrated that compound 20 (SBP-3264) specifically inhibited the Hippo signaling pathway in cultured mammalian cells and possessed favorable pharmacokinetic and pharmacodynamic properties in mice. We show that genetic knockdown and pharmacological inhibition of STK3 and STK4 suppress the proliferation of AML cells in vitro. Thus, SBP-3264 is a valuable chemical probe for understanding the roles of STK3 and STK4 in AML and is a promising candidate for further advancement as a potential therapy.

Charge Assisted S/Se Chalcogen Bonds in SAM Riboswitches: A Combined PDB and ab Initio Study

In this study, we provide experimental (Protein Data Bank (PDB) inspection) and theoretical (RI-MP2/def2-TZVP level of theory) evidence of the involvement of charge assisted chalcogen bonding (ChB) interactions in the recognition and folding mechanisms of S-adenosylmethionine (SAM) riboswitches. Concretely, an initial PDB search revealed several examples where ChBs between S-adenosyl methionine (SAM)/adenosyl selenomethionine (EEM) molecules and uracil (U) bases belonging to RNA take place. While these interactions are usually described as a merely Coulombic attraction between the positively charged S/Se group and RNA, theoretical calculations indicated that the σ holes of S and Se are involved. Moreover, computational models shed light on the strength and directionality properties of the interaction, which was also further characterized from a charge-density perspective using Bader’s “Atoms in Molecules” (AIM) theory, Non-Covalent Interaction plot (NCIplot) visual index, and Natural Bonding Orbital (NBO) analyses. As far as our knowledge extends, this is the first time that ChBs in SAM–RNA complexes have been systematically analyzed, and we believe the results might be useful for scientists working in the field of RNA engineering and chemical biology as well as to increase the visibility of the interaction among the biological community.

Assessing reversible and irreversible binding effects of kinase covalent inhibitors through ADP-Glo assays

Typical enzymatic inhibition assays often demonstrate improved potency for kinase covalent inhibitors compared to reversible inhibitors. This can primarily be attributed to the irreversible mode of action and could affect the evaluations of the ATP-competitive nature of covalent inhibitors, hindering optimization of these compounds. Here, we describe a version of ADP-Glo assay, in which modification of inhibitor incubation time in the presence or absence of ATP enables a quick assessment of relative reversible and irreversible effects of kinase covalent inhibitors.

For complete details on the use and execution of this protocol, please refer to Schröder et al. (2020).

•

Simple enzymatic assays for quick assessment of kinase covalent inhibitors

•

Modification of incubation time allows probing ATP-competitive nature of inhibitors

•

The assays enable the evaluation of relative reversible and irreversible effects

Selenium chalcogen bonds are involved in protein-carbohydrate recognition: a combined PDB and theoretical study

In this manuscript the ability of selenium carbohydrates to undergo chalcogen bonding (ChB) interactions with protein residues has been studied at the RI-MP2/def2-TZVP level of theory. An inspection of the Protein Data Bank (PDB) revealed SeA (A = O, C and S) intermolecular contacts involving Se-pyranose ligands and ASP, TYR, SER and MET residues. Theoretical models were built to analyse the strength and directionality of the interaction together with "Atoms in Molecules" (AIM), Natural Bonding Orbital (NBO) and Non Covalent Interactions plot (NCIplot) analyses, which further assisted in the characterization of the ChBs described herein. We expect that the results from this study will be useful to expand the current knowledge regarding biological ChBs as well as to increase the visibility of the interaction among the carbohydrate chemistry community.

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.

Chalcogen bonds: Hierarchical ab initio benchmark and density functional theory performance study

We have performed a hierarchical ab initio benchmark and DFT performance study of D2 Ch•••A- chalcogen bonds (Ch = S, Se; D, A = F, Cl). The ab initio benchmark study is based on a series of ZORA-relativistic quantum chemical methods [HF, MP2, CCSD, CCSD(T)], and all-electron relativistically contracted variants of Karlsruhe basis sets (ZORA-def2-SVP, ZORA-def2-TZVPP, ZORA-def2-QZVPP) with and without diffuse functions. The highest-level ZORA-CCSD(T)/ma-ZORA-def2-QZVPP counterpoise-corrected complexation energies (ΔECPC ) are converged within 1.1-3.4 kcal mol-1 and 1.5-3.1 kcal mol-1 with respect to the method and basis set, respectively. Next, we used the ZORA-CCSD(T)/ma-ZORA-def2-QZVPP (ΔECPC ) as reference data for analyzing the performance of 13 different ZORA-relativistic DFT approaches in combination with the Slater-type QZ4P basis set. We find that the three-best performing functionals are M06-2X, B3LYP, and M06, with mean absolute errors (MAE) of 4.1, 4.2, and 4.3 kcal mol-1 , respectively. The MAE for BLYP-D3(BJ) and PBE amount to 8.5 and 9.3 kcal mol-1 , respectively.

A Quantitative Molecular Orbital Perspective of the Chalcogen Bond

We have quantum chemically analyzed the structure and stability of archetypal chalcogen-bonded model complexes D2 Ch⋅⋅⋅A- (Ch = O, S, Se, Te; D, A = F, Cl, Br) using relativistic density functional theory at ZORA-M06/QZ4P. Our purpose is twofold: (i) to compute accurate trends in chalcogen-bond strength based on a set of consistent data; and (ii) to rationalize these trends in terms of detailed analyses of the bonding mechanism based on quantitative Kohn-Sham molecular orbital (KS-MO) theory in combination with a canonical energy decomposition analysis (EDA). At odds with the commonly accepted view of chalcogen bonding as a predominantly electrostatic phenomenon, we find that chalcogen bonds, just as hydrogen and halogen bonds, have a significant covalent character stemming from strong HOMO-LUMO interactions. Besides providing significantly to the bond strength, these orbital interactions are also manifested by the structural distortions they induce as well as the associated charge transfer from A- to D2 Ch.

Reliable Comparison of Pnicogen, Chalcogen, and Halogen Bonds in Complexes of 6-OXF2-Fulvene (X = As, Sb, Se, Te, Be, I) With Three Electron Donors

The pnicogen, chalcogen, and halogen bonds between 6-OXF₂-fulvene (X = As, Sb, Se, Te, Br, and I) and three nitrogen-containing bases (FCN, HCN, and NH₃) are compared. For each nitrogen base, the halogen bond is strongest, followed by the pnicogen bond, and the chalcogen bond is weakest. For each type of bond, the binding increases in the FCN < HCN < NH₃ pattern. Both FCN and HCN engage in a bond with comparable strengths and the interaction energies of most bonds are < -6 kcal/mol. However, the strongest base NH₃ forms a much more stable complex, particularly for the halogen bond with the interaction energy going up to -18 kcal/mol. For the same type of interaction, its strength increases as the mass of the central X atom increases. These bonds are different in strength, but all of them are dominated by the electrostatic interaction, with the polarization contribution important for the stronger interaction. The presence of these bonds changes the geometries of 6-OXF₂-fulvene, particularly for the halogen bond formed by NH₃, where the F-X-F arrangement is almost vertical to the fulvene ring.

C-Jun N-terminal kinase inhibitors: Structural insight into kinase-inhibitor complexes

The activation of c-Jun N-terminal kinases (JNKs) plays an important role in physiological processes including neuronal function, immune activity, and development, and thus, JNKs have been a therapeutic target for various diseases such as neurodegenerative diseases, inflammation, and cancer. Efforts to develop JNK-specific inhibitors have been ongoing for several decades. In this process, the structures of JNK in complex with various inhibitors have contributed greatly to the design of novel compounds and to the elucidation of structure-activity relationships. Almost 100 JNK structures with various compounds have been determined. Here we summarize the information gained from these structures and classify the inhibitors into several groups based on the binding mode. These groups include inhibitors in the open conformation and closed conformation of the gatekeeper residue, non-ATP site binders, peptides, covalent inhibitors, and type II kinase inhibitors. Through this work, deep insight into the interaction of inhibitors with JNKs can be gained and this will be helpful for developing novel, potent, and selective inhibitors.

Double Chalcogen Bonds: Crystal Engineering Stratagems via Diffraction and Multinuclear Solid-State Magnetic Resonance Spectroscopy

Group 16 chalcogens potentially provide Lewis-acidic σ-holes, which are able to form attractive supramolecular interactions with electron rich partners through chalcogen bonds. Here, a multifaceted experimental and computational study of a large series of novel chalcogen-bonded cocrystals, prepared using the principles of crystal engineering, is presented. Single-crystal X-ray diffraction studies reveal that dicyanoselenadiazole and dicyanotelluradiazole derivatives work as promising supramolecular synthons with the ability to form double chalcogen bonds with a wide range of electron donors including halides and oxygen- and nitrogen-containing heterocycles. Extensive ⁷⁷ Se and ¹²⁵ Te solid-state nuclear magnetic resonance spectroscopic investigations of cocrystals establish correlations between the NMR parameters of selenium and tellurium and the local chalcogen bonding geometry. The relationships between the electronic environment of the chalcogen bond and the ⁷⁷ Se and ¹²⁵ Te chemical shift tensors were elucidated through a natural localized molecular orbital density functional theory analysis. This systematic study of chalcogen-bond-based crystal engineering lays the foundations for the preparation of the various multicomponent systems and establishes solid-state NMR protocols to detect these interactions in powdered materials.

Direct investigation of chalcogen bonds by multinuclear solid-state magnetic resonance and vibrational spectroscopy

We report a multifaceted experimental and computational study of three self-complementary chalcogen-bond donors as well as a series of seven chalcogen bonded cocrystals.

Theoretical Investigation of Cyano-Chalcogen Dimers and Their Importance in Molecular Recognition

In this manuscript the different noncovalent interactions established between (HYCN)₂ dimers (Y=S, Se and Te) have been studied at the MP2 and CCSD(T) level of theory. Several homodimers have been taken into account, highlighting the capacity of these compounds to act both as electron donor and acceptor. The main properties studied were geometries, binding energy (E_b ), and molecular electrostatic potential (MEP). Given the wide application of chalcogen bonds, and more specifically of cyano-chalcogen moieties in molecular recognition, natural bond orbital (NBO), "atoms-in-molecules" (AIM), and electron density shift (EDS) analysis were also used to analyse the different noncovalent interactions upon complexation. The presence of hydrogen, chalcogen and dipole-dipole interactions was confirmed and their implications on molecular recognition were analysed.

Replacement of Protein Binding-Site Waters Contributes to Favorable Halogen Bond Interactions

Halogen bond interaction between a protein electronegative atom and a ligand halogen atom is increasingly attracting attention in the field of structure-based drug design. Nevertheless, gaps in understanding make it desirable to better examine the role of forces governing the formation of favorable halogen bond interactions, and the development of effective and efficient computational approaches to "design in" favorable halogen bond interactions in lead optimization process are warranted. Here, we analyzed the binding-site water properties of crystal structures with characterized halogen bond interactions between ligand halogen atoms and protein backbone carbonyl groups and, thus, found that halogen atoms involved in halogen bond interactions frequently replace calculated binding-site waters upon ligand binding. Moreover, we observed that the preferential directionality of halogen bond interactions aligns well with the orientations of these replaced waters, and these replaced waters exhibited differential energetic characteristics as compared to waters that are displaced by halogen atoms that do not form halogen bond interactions. Our discovery that replacement of calculated binding-site waters contributes to the formation of favorable halogen bond interactions suggests a practical approach for rational drug design utilizing halogen bond interactions with protein backbone carbonyl groups.

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.

Designing Allele-Specific Inhibitors of Spastin, a Microtubule-Severing AAA Protein

The bump-hole approach is a powerful chemical biology strategy to specifically probe the functions of closely related proteins. However, for many protein families, such as the ATPases associated with diverse cellular activities (AAA), we lack structural data for inhibitor-protein complexes to design allele-specific chemical probes. Here we report the X-ray structure of a pyrazolylaminoquinazoline-based inhibitor bound to spastin, a microtubule-severing AAA protein, and characterize the residues involved in inhibitor binding. We show that an inhibitor analogue with a single-atom hydrogen-to-fluorine modification can selectively target a spastin allele with an engineered cysteine mutation in its active site. We also report an X-ray structure of the fluoro analogue bound to the spastin mutant. Furthermore, analyses of other mutant alleles suggest how the stereoelectronics of the fluorine-cysteine interaction, rather than sterics alone, contribute to the inhibitor-allele selectivity. This approach could be used to design allele-specific probes for studying cellular functions of spastin isoforms. Our data also suggest how tuning stereoelectronics can lead to specific inhibitor-allele pairs for the AAA superfamily.

Copper(i) complexes of functionalized sulfur-containing ligands: structural and theoretical insights into chalcogen bonding

Four new copper(i) thiocyanate complexes were studied using geometrical parameters and the lump–hole approach for justification of the strength and nature of chalcogen bonding.

Halogen bonding in differently charged complexes: basic profile, essential interaction terms and intrinsic σ-hole

Studies on halogen bonds (XB) between organohalogens and their acceptors in crystal structures revealed that the XB donor and acceptor could be differently charged, making it difficult to understand the nature of the interaction, especially the negatively charged donor's electrophilicity and positively charged acceptor's nucleophilicity. In this paper, 9 XB systems mimicking all possibly charged halogen bonding interactions were designed and explored computationally. The results revealed that all XBs could be stable, with binding energies after removing background interaction as strong as -1.2, -3.4, and -8.3 kcal mol-1 for Cl, Br, and I involved XBs respectively. Orbital and dispersion interactions are found to be always attractive while unidirectional intermolecular electron transfer from a XB acceptor to a XB donor occurs in all XB complexes. These observations could be attributed to the intrinsic σ-hole of the XB donor and the intrinsic electronic properties of the XB acceptor regardless of their charge states. Intramolecular charge redistribution inside both the donor and the acceptor is found to be system-dependent but always leads to a more stable XB. Accordingly, this study demonstrates that the orbital-based origin of halogen bonds could successfully interpret the complicated behaviour of differently charged XB complexes, while electrostatic interaction may dramatically change the overall bonding strength. The results should further promote the application of halogens in all related areas.

Chalcogen Bonding "2S-2N Squares" versus Competing Interactions: Exploring the Recognition Properties of Sulfur

Chalcogen bonding (CB) is the focus of increased attention for its applications in medicinal chemistry, materials science, and crystal engineering. However, the origin of sulfur's recognition properties remains controversial, and experimental evidence for supporting theories is still emerging. Here, a comprehensive evaluation of sulfur CB interactions is presented by investigating 2,1,3-benzothiadiazole X-ray crystallographic structures gathered from the Cambridge Structure Database (CSD), Protein Data Bank (PDB), and own laboratory findings. Through the systematic analysis of substituent effects on a subset library of over thirty benzothiadiazole derivatives, the competing interactions have been categorized into four main classes, namely 2S-2N CB square, halogen bonding (XB), S⋅⋅⋅S, and hydrogen-bonding (HB). A geometric model is employed to characterize the 2S-2N CB square motifs and discuss the role of electrostatic, dipole, and orbital contributions toward the interaction.

Switching Between Bicyclic and Linear Peptides - The Sulfhydryl-Specific Linker TPSMB Enables Reversible Cyclization of Peptides

Phage display-selected bicyclic peptides have already shown their great potential for the development as bioactive modulators of therapeutic targets. They can provide enhanced proteolytic stability and improved membrane permeability. Molecular design of new linker molecules has led to a variety of new synthetic approaches for the generation of chemically constrained cyclic peptides. This diversity can be useful for the development of novel peptide-based therapeutic, diagnostic, and scientific tools. Herein, we introduce 1,3,5-tris((pyridin-2-yldisulfanyl)methyl)benzene (TPSMB) as a planar, trivalent, sulfhydryl-specific linker that facilitates reversible cyclization and linearization via disulfide bond formation and cleavage of bicyclic peptides of the format CXnCXnC, where X is any proteinogenic amino acid except cysteine. The rapid and highly sulfhydryl-specific reaction of TPSMB under physiological conditions is demonstrated by selecting bicyclic peptide binders against c-Jun N-terminal kinase 3 (JNK3) as a model target. While model peptides remain stably cyclized for several hours in presence of typical blood levels of glutathione in vitro, high cytosolic concentrations of glutathione linearize these peptides completely within 1 h. We propose that reversible linkers can be useful tools for several technical applications where target affinity depends on the bicyclic structure of the peptide.

Exploring Halogen Bonds in 5-Hydroxytryptamine 2B Receptor-Ligand Interactions

Here, we predicted the potential halogen bonding interaction between compound 2 and the 5-hydroxytryptamine 2B (5-HT_2B) receptor and systematically assessed this interaction via structure-activity relationship analysis and molecular dynamics simulations. A physics-based computational protocol was then developed to further explore the opportunity of "designing in" halogen bonding interactions in structure-based ligand design for the 5-HT_2B receptor, which not only facilitated the identification of previously uncharacterized halogen bonds in known 5-HT_2B ligands but also enabled the rational design of halogen bonding interactions for the optimization of 5-HT_2B ligands. As a proof-of-concept, a series of halogen-substituted analogues of doxepin was synthesized and evaluated, which showed improved in vitro and in vivo potency.

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.

Halogen bonding for the design of inhibitors by targeting the S1 pocket of serine proteases

Halogen bonding (or X bonding) has attracted increasing interest due to its significant role in molecular recognition in biological systems. Trypsin-like serine proteases have many physiological and pathophysiological functions. There is therefore extensive interest in generating specific inhibitors for pharmacological intervention in their enzymatic activity. We study here if it is possible to use halogenated compounds as the P1 group to bind to the S1 specificity pocket of trypsin-like serine proteases to avoid the low bioavailability of the amidine or guanidine P1 group that is typically used in many inhibitors. We used 4-chlorobenzylamine (ClBA), 4-bromobenzylamine (BrBA) and 4-iodobenzylamine (IBA) as probes to test their binding modes to a trypsin-like serine protease, urokinase-type plasminogen activator (uPA), which has been recognized as a marker for breast cancer and an important target for inhibitor development. The results showed that these compounds inhibited uPA with stronger efficacies compared with their non-halogenated analogues. We also determined the high-resolution crystal structures of uPA in complex with BrBA and IBA, respectively. The structures revealed that BrBA bound to the S1 pocket of uPA via halogen bonds, but IBA did not make halogen bonds with uPA, demonstrating that the iodine may not be the best choice as a target moiety for serine proteases. These results advocate halogen bonding, especially bromine bonding, as an efficient strategy for the future design of novel inhibitors against trypsin-like serine proteases to provide strong potency and promote bioavailability.

Intermolecular and very strong intramolecular C–Se⋯O/N chalcogen bonds in nitrophenyl selenocyanate crystals

Different bonding strengths of C–Se⋯O/N chalcogen bonds involved in polymorphic o-NSC (1a/1b) and monomorphic p-NSC (2) result in different thermal properties.

Theoretical and experimental characterization of 1,4-N⋯S σ-hole intramolecular interactions in bioactive N-acylhydrazone derivatives

Sigma-hole (σ-hole) bonds are interactions that are gaining special attention in medicinal chemistry.

2-Alkylsulfanyl-4(5)-aryl-5(4)-heteroarylimidazoles: An Overview on Synthetic Strategies and Biological Activity

2-Alkylsulfanyl-4(5)-aryl-5(4)-heteroarylimidazoles represent an important class of ATP-competitive protein kinase inhibitors, offering the possibility of multiple interactions with different regions of the target enzyme. The necessity of exploring the effects of diverse chemical decorations around the imidazole core prompted the design of several synthetic routes aimed at achieving both efficiency and flexibility. Additionally, the optimization of established protocols and the extensive use of transition metal-catalyzed cross-coupling reactions have been broadening the spectrum of preparative methodologies within the last decade. This review summarizes the progress in the development of synthetic strategies leading to 2-alkylsulfanyl-4(5)-aryl-5(4)-heteroarylimidazoles and 1-alkyl-2-alkylsulfanyl-4(5)-aryl-5(4)-heteroarylimidazoles and offers a glance at the biological activities of this class of compounds.

Halogen Bonds Formed between Substituted Imidazoliums and N Bases of Varying N-Hybridization

Heterodimers are constructed containing imidazolium and its halogen-substituted derivatives as Lewis acid. N in its sp³, sp² and sp hybridizations is taken as the electron-donating base. The halogen bond is strengthened in the Cl < Br < I order, with the H-bond generally similar in magnitude to the Br-bond. Methyl substitution on the N electron donor enhances the binding energy. Very little perturbation arises if the imidazolium is attached to a phenyl ring. The energetics are not sensitive to the hybridization of the N atom. More regular patterns appear in the individual phenomena. Charge transfer diminishes uniformly on going from amine to imine to nitrile, a pattern that is echoed by the elongation of the C-Z (Z=H, Cl, Br, I) bond in the Lewis acid. These trends are also evident in the Atoms in Molecules topography of the electron density. Molecular electrostatic potentials are not entirely consistent with energetics. Although I of the Lewis acid engages in a stronger bond than does H, it is the potential of the latter which is much more positive. The minimum on the potential of the base is most negative for the nitrile even though acetonitrile does not form the strongest bonds. Placing the systems in dichloromethane solvent reduces the binding energies but leaves intact most of the trends observed in vacuo; the same can be said of ∆G in solution.

Discovering protein-ligand chalcogen bonding in the protein data bank using endocyclic sulfur-containing heterocycles as ligand search subsets

The chalcogen bond, the noncovalent, electrostatic attraction between covalently bonded atoms in group 16 and Lewis bases, is present in protein-ligand interactions based on X-ray structures deposited in the Protein Data Bank (PDB). Discovering protein-ligand chalcogen bonding in the PDB employed a strategy that focused on searching the database for protein complexes of five-membered, heterocyclic ligands containing endocyclic sulfur with endo electron-withdrawing groups (isothiazoles; thiazoles; 1,2,3-, 1,2.4-, 1,2,5-, 1,3,4-thiadiazoles) and thiophenes with exo electron-withdrawing groups, e.g., 2-chloro, 2-bromo, 2-amino, 2-alkylthio. Out of 930 ligands investigated, 33 or 3.5% have protein-ligand S---O interactions of which 31 are chalcogen bonds and two appear to be S---HO hydrogen bonds. The bond angles for some of the chalcogen bonds found in the PDB are less than 90°, and an electrostatic model is proposed to explain this phenomenon.

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.