In Silico Targeting of Fascin Protein for Cancer Therapy: Benchmarking, Virtual Screening and Molecular Dynamics Approaches.

Fascin is an actin-bundling protein overexpressed in various invasive metastatic carcinomas through promoting cell migration and invasion. Therefore, blocking Fascin binding sites is considered a vital target for antimetastatic drugs. This inspired us to find new Fascin binding site blockers. First, we built an active compound set by collecting reported small molecules binding to Fascin's binding site 2. Consequently, a high-quality decoys set was generated employing DEKOIS 2.0 protocol to be applied in conducting the benchmarking analysis against the selected Fascin structures. Four docking programs, MOE, AutoDock Vina, VinaXB, and PLANTS were evaluated in the benchmarking study. All tools indicated better-than-random performance reflected by their pROC-AUC values against the Fascin crystal structure (PDB: ID 6I18). Interestingly, PLANTS exhibited the best screening performance and recognized potent actives at early enrichment. Accordingly, PLANTS was utilized in the prospective virtual screening effort for repurposing FDA-approved drugs (DrugBank database) and natural products (NANPDB). Further assessment via molecular dynamics simulations for 100 ns endorsed Remdesivir (DrugBank) and NANPDB3 (NANPDB) as potential binders to Fascin binding site 2. In conclusion, this study delivers a model for implementing a customized DEKOIS 2.0 benchmark set to enhance the VS success rate against new potential targets for cancer therapies.

Discovery of Ircinianin Lactones B and C-Two New Cyclic Sesterterpenes from the Marine Sponge Ircinia wistarii.

Two new ircinianin-type sesterterpenoids, ircinianin lactone B and ircinianin lactone C (7 and 8), together with five known entities from the ircinianin compound family (1, 3-6) were isolated from the marine sponge Ircinia wistarii. Ircinianin lactones B and C (7 and 8) represent new ircinianin terpenoids with a modified oxidation pattern. Despite their labile nature, the structures could be established using a combination of spectroscopic data, including HRESIMS and 1D/2D NMR techniques, as well as computational chemistry and quantum-mechanical calculations. In a broad screening approach for biological activity, the class-defining compound ircinianin (1) showed moderate antiprotozoal activity against Plasmodium falciparum (IC50 25.4 µM) and Leishmania donovani (IC50 16.6 µM).

Predicting the Magnitude of σ-Holes Using VmaxPred, a Fast and Efficient Tool Supporting the Application of Halogen Bonds in Drug Discovery.

Halogen bonding as a modern molecular interaction has received increasing attention not only in materials sciences but also in biological systems and drug discovery. Thus, there is a growing demand for fast, efficient, and easily applicable tailor-made tools supporting the use of halogen bonds in molecular design and medicinal chemistry. The potential strength of a halogen bond is dependent on several properties of the σ-hole donor, e.g., a (hetero)aryl halide, and the σ-hole acceptor, a nucleophile with n or π electron density. Besides the influence of the interaction geometry and the type of acceptor, significant tuning effects on the magnitude of the σ-hole can be observed, caused by different (hetero)aromatic scaffolds and their substitution patterns. The most positive electrostatic potential on the isodensity surface ( Vmax), representing the σ-hole, has been widely used as the standard descriptor for the magnitude of the σ-hole and the strength of the halogen bond. Calculation of Vmax using quantum-mechanical methods at a reasonable level of theory is time-consuming and thus not applicable for larger numbers of compounds in drug discovery projects. Herein we present a tool for the prediction of this descriptor based on a machine-learned model with a speedup of 5 to 6 orders of magnitude relative to MP2 quantum-mechanical calculations. According to the test set, the squared correlation coefficient is greater than 0.94.

Scaffold Effects on Halogen Bonding Strength.

Halogen bonds have become increasingly popular interactions in molecular design and drug discovery. One of the key features is the strong dependence of the size and magnitude of the halogen's σ-hole on the chemical environment of the ligand. The term σ-hole refers to a region of lower electronic density opposite to a covalent bond, e.g., the C-X bond. It is typically (but not always) associated with a positive electrostatic potential in close proximity to the extension of the covalent bond. Herein, we use a variety of 30 nitrogen-bearing heterocycles, halogenated systematically by chlorine, bromine, or iodine, yielding 468 different ligands that are used to exemplify scaffold effects on halogen bonding strength. As a template interaction partner, we have chosen N-methylacetamide representing the ubiquitously present protein backbone. Adduct formation energies were obtained at a MP2/TZVPP level of theory. We used the local maximum of the electrostatic potential on the molecular surface in close proximity to the σ-hole, V S,max, as a descriptor for the magnitude of the positive electrostatic potential characterizing the tuning of the σ-hole. Free optimization of the complexes gave reasonable correlations with V S,max but was found to be of limited use because considerable numbers of chlorinated and brominated ligands lost their halogen bond or showed significant secondary interactions. Thus, starting from a close to optimal geometry of the halogen bond, we used distance scans to obtain the best adduct formation energy for each complex. This approach provided superior results for all complexes exhibiting correlations with R2 > 0.96 for each individual halogen. We evaluated the dependence of V S,max from the molecular surface onto which the positive electrostatic potential is projected, altering the isodensity values from 0.001 au to 0.050 au. Interestingly, the best overall fit using a third-order polynomial function (R2 = 0.99, RMSE = 0.562 kJ/mol) with rather smooth transitions between all halogens was obtained for V S,max calculated from an isodensity surface at 0.014 au.

Fluorescence polarization-based competition binding assay for c-Jun N-terminal kinases 1 and 2.

In order to evaluate the isoform selectivity of novel inhibitors within the c-Jun N-terminal kinase (JNK) family, a fluorescence polarization-based competition binding assay, previously developed for JNK3, was extended to the other isoforms JNK1 and JNK2. The assay is based on the displacement of a versatile fluorescent pyridinylimidazole-based probe and was validated by testing the precursor of the probe as well as standard JNK inhibitors.

Fluorescence polarization-based assays for detecting compounds binding to inactive c-Jun N-terminal kinase 3 and p38α mitogen-activated protein kinase.

Two fluorescein-labeled pyridinylimidazoles were synthesized and evaluated as probes for the binding affinity determination of potential kinase inhibitors to the c-Jun N-terminal kinase 3 (JNK3) and p38α mitogen-activated protein kinase (MAPK). Fluorescence polarization (FP)-based competition binding assays were developed for both enzymes using 1-(3',6'-dihydroxy-3-oxo-3H-spiro[isobenzofuran-1,9'-xanthen]-5-yl)-3-(4-((4-(4-(4-fluorophenyl)-2-(methylthio)-1H-imidazol-5-yl)pyridin-2-yl)amino)phenyl)thiourea (5) as an FP probe (JNK3: Kd = 3.0 nM; p38α MAPK: Kd = 5.7 nM). The validation of the assays with known inhibitors of JNK3 and p38α MAPK revealed that both FP assays correlate very well with inhibition data received by the activity assays. This, in addition to the viability of both FP-based binding assays for the high-throughput screening procedure, makes the assays suitable as inexpensive prescreening protocols for JNK3 and p38α MAPK inhibitors.

Using Surface Scans for the Evaluation of Halogen Bonds toward the Side Chains of Aspartate, Asparagine, Glutamate, and Glutamine.

Using halogen-specific Connolly type molecular surfaces, we herein invented a new type of surface-based interaction analysis employed for the study of halogen bonding toward model systems of biologically relevant carboxylates (ASP/GLU) and carboxamides (ASN/GLN). Database mining and statistical assessment of the PDB revealed that such interactions are widely underrepresented at the moment. We observed important distance-dependent adaptions of the binding modes of halobenzenes from a preferential oxygen-directed to a bifurcated interaction geometry of the carboxylate. In addition, halogen···π contacts perpendicular to the nitrogen atom of the carboxamide become increasingly important for the lighter halogens. Our analysis on a MP2/TZVPP level of theory is backed by CCSD(T)/CBS reference calculations. To put the vast interaction energies into perspective, we also performed COSMO-RS calculations of the solvation free energy. Facilitating the visualization of our results mapped onto any binding site of choice, we aim to inspire more design studies showcasing these underrepresented interactions.

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

We target the gatekeeper MET146 of c-Jun N-terminal kinase 3 (JNK3) to exemplify the applicability of X···S halogen bonds in molecular design using computational, synthetic, structural and biophysical techniques. In a designed series of aminopyrimidine-based inhibitors, we unexpectedly encounter a plateau of affinity. Compared to their QM-calculated interaction energies, particularly bromine and iodine fail to reach the full potential according to the size of their σ-hole. Instead, mutation of the gatekeeper residue into leucine, alanine, or threonine reveals that the heavier halides can significantly influence selectivity in the human kinome. Thus, we demonstrate that, although the choice of halogen may not always increase affinity, it can still be relevant for inducing selectivity. Determining the crystal structure of the iodine derivative in complex with JNK3 (4X21) reveals an unusual bivalent halogen/chalcogen bond donated by the ligand and the back-pocket residue MET115. Incipient repulsion from the too short halogen bond increases the flexibility of Cε of MET146, whereas the rest of the residue fails to adapt being fixed by the chalcogen bond. This effect can be useful to induce selectivity, as the necessary combination of methionine residues only occurs in 9.3% of human kinases, while methionine is the predominant gatekeeper (39%).

Machine learning estimates of natural product conformational energies.

Machine learning has been used for estimation of potential energy surfaces to speed up molecular dynamics simulations of small systems. We demonstrate that this approach is feasible for significantly larger, structurally complex molecules, taking the natural product Archazolid A, a potent inhibitor of vacuolar-type ATPase, from the myxobacterium Archangium gephyra as an example. Our model estimates energies of new conformations by exploiting information from previous calculations via Gaussian process regression. Predictive variance is used to assess whether a conformation is in the interpolation region, allowing a controlled trade-off between prediction accuracy and computational speed-up. For energies of relaxed conformations at the density functional level of theory (implicit solvent, DFT/BLYP-disp3/def2-TZVP), mean absolute errors of less than 1 kcal/mol were achieved. The study demonstrates that predictive machine learning models can be developed for structurally complex, pharmaceutically relevant compounds, potentially enabling considerable speed-ups in simulations of larger molecular structures.

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.

pROC-Chemotype Plots Enhance the Interpretability of Benchmarking Results in Structure-Based Virtual Screening.

Recently, we have reported a systematic comparison of molecular preparation protocols (using MOE or Maestro) in combination with two docking tools (GOLD or Glide), employing our DEKOIS 2.0 benchmark sets. Herein, we demonstrate how comparable settings of data preparation protocols can affect the profile and AUC of pROC curves based on variations in chemotype enrichment. We show how the recognition of different classes of chemotypes can affect the docking performance, particularly in the early enrichment, and monitor changes in this recognition behavior based on score normalization and rescoring strategies. For this, we have developed "pROC-Chemotype", which is an automated protocol that matches and visualizes ligand chemotype information together with potency classes in the pROC profiles obtained by docking. This tool enhances the understanding of the influence of chemotype recognition in early enrichment, but also reveals trends of impaired recognition of chemotype classes at the end of the score-ordered rank. Identifying such issues helps to devise score-normalization strategies to overcome this potential bias in an intuitive manner. Furthermore, strong perturbations in chemotype ranking between different methods can help to identify the underlying reasons (e.g., changes in the protonation/tautomerization state). It also assists in the selection of appropriate scoring functions that are capable to retrieve more potent and diverse hits. In summary, we demonstrate how this new tool can be utilized to identify and highlight chemotype-specific behavior, e.g., in dataset preparation. This can help to overcome some chemistry-related bias in virtual screening campaigns. pROC-Chemotype is made freely available at www.dekois.com.

Kinetic mechanism of p53 oncogenic mutant aggregation and its inhibition.

Aggregation of destabilized mutants of the tumor suppressor p53 is a major route for its loss of activity. In order to assay drugs that inhibit aggregation of p53, we established the basic kinetics of aggregation of its core domain, using the mutant Y220C that has a mutation-induced, druggable cavity. Aggregation monitored by light scattering followed lag kinetics. Electron microscopy revealed the formation of small aggregates that subsequently grew to larger amorphous aggregates. The kinetics of aggregation produced surprising results: progress curves followed either by the binding of Thioflavin T or the fluorescence of the protein at 340 nm fitted well to simple two-step sequential first-order lag kinetics with rate constants k(1) and k(2) that were independent of protein concentration, and not to classical nucleation-growth. We suggest a mechanism of first-order formation of an aggregation competent state as being rate determining followed by rapid polymerization with the higher order kinetics. By measuring the inhibition kinetics of k(1) and k(2), we resolved that the process with the higher rate constant followed that of the lower. Further, there was only partial inhibition of k(1) and k(2), which showed two parallel pathways of aggregation, one via a state that requires unfolding of the protein and the other of partial unfolding with the ligand still bound. Inhibition kinetics of ligands provides a useful tool for probing an aggregation mechanism.

Lithocholic acid is an endogenous inhibitor of MDM4 and MDM2.

The proteins MDM2 and MDM4 are key negative regulators of the tumor suppressor protein p53, which are frequently upregulated in cancer cells. They inhibit the transactivation activity of p53 by binding separately or in concert to its transactivation domain. MDM2 is also a ubiquitin ligase that leads to the degradation of p53. Accordingly, MDM2 and MDM4 are important targets for drugs to inhibit their binding to p53. We found from in silico screening and confirmed by experiment that lithocholic acid (LCA) binds to the p53 binding sites of both MDM2 and MDM4 with a fivefold preference for MDM4. LCA is an endogenous steroidal bile acid, variously reported to have both carcinogenic and apoptotic activities. The comparison of LCA effects on apoptosis in HCT116 p53(+/+) vs. p53(-/-) cells shows a predominantly p53-mediated induction of caspase-3/7. The dissociation constants are in the µM region, but only modest inhibition of binding of MDM2 and MDM4 is required to negate their upregulation because they have to compete with transcriptional coactivator p300 for binding to p53. Binding was weakened by structural changes in LCA, and so it may be a natural ligand of MDM2 and MDM4, raising the possibility that MDM proteins may be sensors for specific steroids.

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.

Targeting histidine side chains in molecular design through nitrogen-halogen bonds.

Halogen bonds are directional noncovalent interactions that can be used to target electron donors in a protein binding site. In this study, we employ quantum chemical calculations to explore halogen···nitrogen contacts involving histidine side chains. We characterize the energetics on the MP2 level of theory using SCS-MP2 and CCSD(T)/CBS as reference calculations and elucidate their energy profile in suboptimal geometries. We derive simple rules allowing medicinal chemists and chemical biologists to easily determine preferred areas of interaction in a binding site and exploit them for scaffold decoration and design. Our work shows that nitrogen-halogen bonds are valuable interactions that are this far underexploited in patent applications, lead structure, and clinical candidate selection. We highlight their potential to increase binding affinities and suggest that they can significantly contribute to inducing and tuning subtype selectivities.

Evaluation and optimization of virtual screening workflows with DEKOIS 2.0--a public library of challenging docking benchmark sets.

The application of molecular benchmarking sets helps to assess the actual performance of virtual screening (VS) workflows. To improve the efficiency of structure-based VS approaches, the selection and optimization of various parameters can be guided by benchmarking. With the DEKOIS 2.0 library, we aim to further extend and complement the collection of publicly available decoy sets. Based on BindingDB bioactivity data, we provide 81 new and structurally diverse benchmark sets for a wide variety of different target classes. To ensure a meaningful selection of ligands, we address several issues that can be found in bioactivity data. We have improved our previously introduced DEKOIS methodology with enhanced physicochemical matching, now including the consideration of molecular charges, as well as a more sophisticated elimination of latent actives in the decoy set (LADS). We evaluate the docking performance of Glide, GOLD, and AutoDock Vina with our data sets and highlight existing challenges for VS tools. All DEKOIS 2.0 benchmark sets will be made accessible at http://www.dekois.com.

Trisubstituted and tetrasubstituted pyrazolines as a novel class of cell-growth inhibitors in tumor cells with wild type p53.

Derivatives with scaffolds of 1,3,5-tri-substituted pyrazoline and 1,3,4,5-tetra-substituted pyrazoline were synthesized and tested for their inhibitory effects versus the p53(+/+) HCT116 and p53(-/-) H1299 human tumor cell lines. Several compounds were active against the two cell lines displaying IC50 values in the low micromolar range with a clearly more pronounced effect on the p53(+/+) HCT116 cells. The compound class shows excellent developability due to the modular synthesis, allowing independent optimization of all three to four key substituents to improve the properties of the molecules.

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-CF2X 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(εNHCOCF2X) exhibits an improved inhibition constant Ki compared to the unmodified peptide. Decreasing Ki 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(εNHCOCF2X) to Abu(γNHCOCF2X) 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.

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.

Halogen-enriched fragment libraries as leads for drug rescue of mutant p53.

The destabilizing p53 cancer mutation Y220C creates a druggable surface crevice. We developed a strategy exploiting halogen bonding for lead discovery to stabilize the mutant with small molecules. We designed halogen-enriched fragment libraries (HEFLibs) as starting points to complement classical approaches. From screening of HEFLibs and subsequent structure-guided design, we developed substituted 2-(aminomethyl)-4-ethynyl-6-iodophenols as p53-Y220C stabilizers. Crystal structures of their complexes highlight two key features: (i) a central scaffold with a robust binding mode anchored by halogen bonding of an iodine with a main-chain carbonyl and (ii) an acetylene linker, enabling the targeting of an additional subsite in the crevice. The best binders showed induction of apoptosis in a human cancer cell line with homozygous Y220C mutation. Our structural and biophysical data suggest a more widespread applicability of HEFLibs in drug discovery.