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.

High-confidence placement of low-occupancy fragments into electron density using the anomalous signal of sulfur and halogen atoms

Fragment-based drug design using X-ray crystallography is a powerful technique to enable the development of new lead compounds, or probe molecules, against biological targets. This study addresses the need to determine fragment binding orientations for low-occupancy fragments with incomplete electron density, an essential step before further development of the molecule. Halogen atoms play multiple roles in drug discovery due to their unique combination of electronegativity, steric effects and hydrophobic properties. Fragments incorporating halogen atoms serve as promising starting points in hit-to-lead development as they often establish halogen bonds with target proteins, potentially enhancing binding affinity and selectivity, as well as counteracting drug resistance. Here, the aim was to unambiguously identify the binding orientations of fragment hits for SARS-CoV-2 nonstructural protein 1 (nsp1) which contain a combination of sulfur and/or chlorine, bromine and iodine substituents. The binding orientations of carefully selected nsp1 analogue hits were focused on by employing their anomalous scattering combined with Pan-Dataset Density Analysis (PanDDA). Anomalous difference Fourier maps derived from the diffraction data collected at both standard and long-wavelength X-rays were compared. The discrepancies observed in the maps of iodine-containing fragments collected at different energies were attributed to site-specific radiation-damage stemming from the strong X-ray absorption of I atoms, which is likely to cause cleavage of the C-I bond. A reliable and effective data-collection strategy to unambiguously determine the binding orientations of low-occupancy fragments containing sulfur and/or halogen atoms while mitigating radiation damage is presented.

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.

New Hybrid Tetrahydropyrrolo[3,2,1-ij]quinolin-1-ylidene-2-thioxothiazolidin-4-ones as New Inhibitors of Factor Xa and Factor XIa: Design, Synthesis, and In Silico and Experimental Evaluation

Despite extensive research in the field of thrombotic diseases, the prevention of blood clots remains an important area of study. Therefore, the development of new anticoagulant drugs with better therapeutic profiles and fewer side effects to combat thrombus formation is still needed. Herein, we report the synthesis and evaluation of novel pyrroloquinolinedione-based rhodanine derivatives, which were chosen from 24 developed derivatives by docking as potential molecules to inhibit the clotting factors Xa and XIa. For the synthesis of new hybrid derivatives of pyrrolo[3,2,1-ij]quinoline-2-one, we used a convenient structural modification of the tetrahydroquinoline fragment by varying the substituents in positions 2, 4, and 6. In addition, the design of target molecules was achieved by alkylating the amino group of the rhodanine fragment with propargyl bromide or by replacing the rhodanine fragment with 2-thioxoimidazolidin-4-one. The in vitro testing showed that eight derivatives are capable of inhibiting both coagulation factors, two compounds are selective inhibitors of factor Xa, and two compounds are selective inhibitors of factor XIa. Overall, these data indicate the potential anticoagulant activity of these molecules through the inhibition of the coagulation factors Xa and XIa.

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.

Comprehensive analysis of commercial fragment libraries

Screening of fragment libraries is a valuable approach to the drug discovery process. The quality of the library is one of the keys to success, and more particularly the design or choice of a library has to meet the specificities of the research program. In this study, we made an inventory of the commercial fragment libraries and we established a methodology which allows any library to be positioned in relation to the complete offer currently on the market, by addressing the following questions: does this chemical library look like another chemical library? What is the coverage of the current chemical space by this chemical library? What are the characteristic structural features of the fragments of this chemical library? We based our analysis on 2D and 3D chemical descriptors, framework class generation and the generative topographic map. We identified 59 270 scaffolds, which can be searched in a dedicated web site (https://gtmfrag.drugdesign.unistra.fr) and developed a model which accounts for fragment diversity while being easy to interpret (download at 10.5281/zenodo.5534434).

A Biological Take on Halogen Bonding and Other Non-Classical Non-Covalent Interactions

Classical hydrogen bonds have, for many decades, been the dominant non-covalent interaction in the toolbox that chemists and chemical engineers have used to design and control the structures of compounds and molecular assemblies as novel materials. Recently, a set of non-classical non-covalent (NC-NC) interactions have emerged that exploit the properties of the Group IV, V, VI, and VII elements of the periodic table (the tetrel, pnictogen, chalcogen, and halogen bonds, respectively). Our research group has been characterizing the prevalence, geometric constraints, and structure-function relationship specifically of the halogen bond in biological systems. We have been particularly interested in exploiting the biological halogen bonds (or BXBs) to control the structures, stabilities, and activities of biomolecules, including the DNA Holliday junction and enzymes. In this review, we first provide a set of criteria for how to determine whether BXBs or any other NC-NC interactions would have biological relevance. We then navigate the trail of studies that had led us from an initial, very biological question to our current point in the journey to establish BXBs as a tool for biomolecular engineering. Finally, we close with a perspective on future directions for this line of research.

Combined MD/QTAIM techniques to evaluate ligand-receptor interactions. Scope and limitations

In medicinal chemistry, it is extremely important to evaluate, as accurately as possible, the molecular interactions involved in the formation of different ligand-receptor (L-R) complexes. Evaluating the different molecular interactions by quantum mechanics calculations is not a simple task, since formation of an L-R complex is a dynamic process. In this case, the use of combined techniques of molecular dynamics (MD) and quantum calculations is one the best possible approaches. In this work we report a comparative study using combined MD and QTAIM (Quantum Theory of Atoms In Molecules) calculations for five biological systems with different levels of structural complexity. We have studied Acetylcholinesterase (AChE), D2 Dopamine Receptor (D2DR), beta Secretase (BACE1), Dihydrofolate Reductase (DHFR) and Sphingosine Kinase 1 (SphK1). In these molecular targets, we have analyzed different ligands with diverse structural characteristics. The inhibitory activities of most of them have been previously measured in our laboratory. Our results indicate that QTAIM calculations can be extremely useful for in silico studies. It is possible to obtain very accurate information about the strength of the molecular interactions that stabilize the formation of the different L-R complexes. Better correlations can be obtained between theoretical and experimental data by using QTAIM calculations, allowing us to discriminate among ligands with similar affinities. QTAIM analysis gives fairly accurate information for weak interactions which are not well described by MD simulations. QTAIM study also allowed us to evaluate and determine which parts of the ligand need to be modified in order to increase its interactions with the molecular target. In this study we have discussed the importance of combined MD/QTAIM calculations for this type of simulations, showing their scopes and limitations.

Computational study of the substituent effect of halogenated fused-ring heteroaromatics on halogen bonding

Halogen bonding (XB) has been applied in many fields from crystal engineering to medicinal chemistry. Compared with the well-studied XB of simple halogenated aromatics, little research has been done on the XB of halogenated fused-ring heteroaromatics, a prevalent substructure in organic compounds. With 1H-pyrrolo[3,2-b]pyridines (PPs) as examples of novel fused-ring heteroaromatics with hydrogen bond donor and acceptor and XB donor, the XB formed by the halogenated heteroaromatics was explored in this study. With 4 different substituents, viz., -CH₃, -NH₂, -F, and -CONH₂, at different positions, 339 derivatives of brominated PP (Br-PP) were designed for calculating their electrostatic potential of the σ-hole of the halogen atom (V_S,max) and binding energy with ammonia as XB acceptor (E_int) at M06-2X/6-311++G(d,p) level by PCM model in dichloromethane. The calculated V_S,max values ranging from -1.3 to 35.1 kcal/mol and the calculated E_int ranging from -0.82 to -2.37 kcal/mol demonstrated that the XB is complicated and highly tunable. Noticeably, the electron-withdrawing substituents, especially at ortho-position, do not always increase the values of V_S,max, while the electron-donating substituents do not always decrease V_S,max. Similar results were observed from the calculation on 339 iodinated PPs at M06-2X/6-311++G(d,p) level. The complexity of the XB formed by the halogenated fused ring heteroaromatics indicated a great potential of tuning its strength by different substituents at different positions and revealed a necessity of quantum chemistry calculation for predicting the XB.Graphical abstract.

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.

Embracing the Diversity of Halogen Bonding Motifs in Fragment-Based Drug Discovery-Construction of a Diversity-Optimized Halogen-Enriched Fragment Library

Halogen bonds have recently gained attention in life sciences and drug discovery. However, it can be difficult to harness their full potential, when newly introducing them into an established hit or lead structure by molecular design. A possible solution to overcome this problem is the use of halogen-enriched fragment libraries (HEFLibs), which consist of chemical probes that provide the opportunity to identify halogen bonds as one of the main features of the binding mode. Initially, we have suggested the HEFLibs concept when constructing a focused library for finding p53 mutant stabilizers. Herein, we broaden and extent this concept aiming for a general HEFLib comprising a huge diversity of binding motifs and, thus, increasing the applicability to various targets. Using the construction principle of feature trees, we represent each halogenated fragment by treating all simple to complex substituents as modifiers of the central (hetero)arylhalide. This approach allows us to focus on the proximal binding interface around the halogen bond and, thus, its integration into a network of interactions based on the fragment's binding motif. As a first illustrative example, we generated a library of 198 fragments that unifies a two-fold strategy: Besides achieving a diversity-optimized basis of the library, we have extended this "core" by structurally similar "satellite compounds" that exhibit quite different halogen bonding interfaces. Tuning effects, i.e., increasing the magnitude of the σ-hole, can have an essential influence on the strength of the halogen bond. We were able to implement this key feature into the diversity selection, based on the rapid and efficient prediction of the highest positive electrostatic potential on the electron isodensity surface, representing the σ-hole, by VmaxPred.

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.

A systematic analysis of atomic protein-ligand interactions in the PDB

As the protein databank (PDB) recently passed the cap of 123 456 structures, it stands more than ever as an important resource not only to analyze structural features of specific biological systems, but also to study the prevalence of structural patterns observed in a large body of unrelated structures, that may reflect rules governing protein folding or molecular recognition. Here, we compiled a list of 11 016 unique structures of small-molecule ligands bound to proteins - 6444 of which have experimental binding affinity - representing 750 873 protein-ligand atomic interactions, and analyzed the frequency, geometry and impact of each interaction type. We find that hydrophobic interactions are generally enriched in high-efficiency ligands, but polar interactions are over-represented in fragment inhibitors. While most observations extracted from the PDB will be familiar to seasoned medicinal chemists, less expected findings, such as the high number of C-H···O hydrogen bonds or the relatively frequent amide-π stacking between the backbone amide of proteins and aromatic rings of ligands, uncover underused ligand design strategies.

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.

Identification of a Small Molecule Inhibitor of RAD52 by Structure-Based Selection

It has been reported that inhibition of RAD52 either by specific shRNA or a small peptide aptamer induced synthetic lethality in tumor cell lines carrying BRCA1 and BRCA2 inactivating mutations. Molecular docking was used to screen two chemical libraries: 1) 1,217 FDA approved drugs, and 2) 139,735 drug-like compounds to identify candidates for interacting with DNA binding domain of human RAD52. Thirty six lead candidate compounds were identified that were predicted to interfere with RAD52 –DNA binding. Further biological testing confirmed that 9 of 36 candidate compounds were able to inhibit the binding of RAD52 to single-stranded DNA in vitro. Based on molecular binding combined with functional assays, we propose a model in which the active compounds bind to a critical “hotspot” in RAD52 DNA binding domain 1. In addition, one of the 9 active compounds, adenosine 5’-monophosphate (A5MP), and also its mimic 5-aminoimidazole-4-carboxamide ribonucleotide (AICAR) 5’ phosphate (ZMP) inhibited RAD52 activity in vivo and exerted synthetic lethality against BRCA1 and BRCA2–mutated carcinomas. These data suggest that active, inhibitory RAD52 binding compounds could be further refined for efficacy and safety to develop drugs inducing synthetic lethality in tumors displaying deficiencies in BRCA1/2-mediated homologous recombination.

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.

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.

Halogen bonding in the crystal structures of 1,2-diiodo alkenes

The crystal structures of 1,2-diiodoolefins are governed by extensive halogen bonding involving I⋯I, I⋯O and I⋯C interactions. These interactions are discussed considering nine new crystal structures determined in our laboratory and several additional crystal structures from the literature that complement our data.

Inhibitors of c-Jun N-terminal kinases: an update

The c-Jun N-terminal kinases (JNKs) are serine/threonine kinases implicated in the pathogenesis of various diseases. Recent advances in the development of novel inhibitors of JNKs will be reviewed. Significant progress in the design of JNK inhibitors displaying selectivity versus other kinases has been achieved within the past 4 years. However, the development of isoform selective JNK inhibitors is still an open task.

Thermodynamic parameters for binding of some halogenated inhibitors of human protein kinase CK2

The interaction of human CK2α with a series of tetrabromobenzotriazole (TBBt) and tetrabromobenzimidazole (TBBz) analogs, in which one of the bromine atoms proximal to the triazole/imidazole ring is replaced by a methyl group, was studied by biochemical (IC50) and biophysical methods (thermal stability of protein-ligand complex monitored by DSC and fluorescence). Two newly synthesized tri-bromo derivatives display inhibitory activity comparable to that of the reference compounds, TBBt and TBBz, respectively. DSC analysis of the stability of protein-ligand complexes shows that the heat of ligand binding (Hbind) is driven by intermolecular electrostatic interactions involving the triazole/imidazole ring, as indicated by a strong correlation between Hbind and ligand pKa. Screening, based on fluorescence-monitored thermal unfolding of protein-ligand complexes, gave comparable results, clearly identifying ligands that most strongly bind to the protein. Overall results, additionally supported by molecular modeling, confirm that a balance of hydrophobic and electrostatic interactions contribute predominantly, relative to possible intermolecular halogen bonding, in binding of the ligands to the CK2α ATP-binding site.

Force Field Model of Periodic Trends in Biomolecular Halogen Bonds

The study of the noncovalent interaction now defined as a halogen bond (X-bond) has become one of the fastest growing areas in experimental and theoretical chemistry--its applications as a design tool are highly extensive. The significance of the interaction in biology has only recently been recognized, but has now become important in medicinal chemistry. We had previously derived a set of empirical potential energy functions to model the structure-energy relationships for bromines in biomolecular X-bonds (BXBs). Here, we have extended this force field for BXBs (ffBXB) to the halogens (Cl, Br, and I) that are commonly seen to form stable X-bonds. The ffBXB calculated energies show a remarkable one-to-one linear relationship to explicit BXB energies determined from an experimental DNA junction system, thereby validating the approach and the model. The resulting parameters allow us to interpret the stabilizing effects of BXBs in terms of well-defined physical properties of the halogen atoms, including their size, shape, and charge, showing periodic trends that are predictable along the Group VII column of elements. Consequently, we have established the ffBXB as an accurate computational tool that can be applied, for example, for the design of new therapeutic compounds against clinically important targets and new biomolecular-based materials.