Structure-based discovery of potent WD repeat domain 5 inhibitors that demonstrate efficacy and safety in preclinical animal models.

WD repeat domain 5 (WDR5) is a core scaffolding component of many multiprotein complexes that perform a variety of critical chromatin-centric processes in the nucleus. WDR5 is a component of the mixed lineage leukemia MLL/SET complex and localizes MYC to chromatin at tumor-critical target genes. As a part of these complexes, WDR5 plays a role in sustaining oncogenesis in a variety of human cancers that are often associated with poor prognoses. Thus, WDR5 has been recognized as an attractive therapeutic target for treating both solid and hematological tumors. Previously, small-molecule inhibitors of the WDR5-interaction (WIN) site and WDR5 degraders have demonstrated robust in vitro cellular efficacy in cancer cell lines and established the therapeutic potential of WDR5. However, these agents have not demonstrated significant in vivo efficacy at pharmacologically relevant doses by oral administration in animal disease models. We have discovered WDR5 WIN-site inhibitors that feature bicyclic heteroaryl P7 units through structure-based design and address the limitations of our previous series of small-molecule inhibitors. Importantly, our lead compounds exhibit enhanced on-target potency, excellent oral pharmacokinetic (PK) profiles, and potent dose-dependent in vivo efficacy in a mouse MV4:11 subcutaneous xenograft model by oral dosing. Furthermore, these in vivo probes show excellent tolerability under a repeated high-dose regimen in rodents to demonstrate the safety of the WDR5 WIN-site inhibition mechanism. Collectively, our results provide strong support for WDR5 WIN-site inhibitors to be utilized as potential anticancer therapeutics.

Structure-based design, synthesis, and biological characterization of indolylarylsulfone derivatives as novel human immunodeficiency virus type 1 inhibitors with potent antiviral activities and favorable drug-like profiles.

In the current antiretroviral landscape, continuous efforts are still needed to search for novel chemotypes of human immunodeficiency virus type 1 (HIV-1) inhibitors with improved drug resistance profiles and favorable drug-like properties. Herein, we report the design, synthesis, biological characterization, and druggability evaluation of a class of non-nucleoside reverse transcriptase inhibitors. Guided by the available crystallographic information, a series of novel indolylarylsulfone derivatives were rationally discovered via the substituent decorating strategy to fully explore the chemical space of the entrance channel. Among them, compound 11h bearing the cyano-substituted benzyl moiety proved to be the most effective inhibitor against HIV-1 wild-type and mutant strains (EC50 = 0.0039-0.338 µM), being far more potent than or comparable to etravirine and doravirine. Besides, 11h did not exhibit cytotoxicity at the maximum test concentration. Meanwhile, the binding target of 11h was further confirmed to be reverse transcriptase (IC50 = 0.055 µM). Preliminary structure-activity relationship were discussed to guide further optimization work. Molecular docking and dynamics simulation studies were investigated in detail to rationalize the biological evaluation results. Further drug-likeness assessment indicated that 11h possessed excellent physicochemical properties. Moreover, no apparent hERG blockade liability and cytochrome P450 inhibition were observed for 11h. Notably, 11h was characterized by favorable in vitro metabolic stability with moderate clearance rates and long half-lives in human plasma and liver microsomes. Overall, 11h holds great promise as an ideal Anti-HIV-1 lead compound due to its potent antiviral efficacy, low toxicity, and favorable drug-like profiles.

Structure-based design and synthesis of sulfonylureas as novel NLRP3 inhibitors for Alzheimer's disease.

Alzheimer's disease (AD) remains an incurable neurodegenerative condition that poses a threat to humanity. Immune signaling in the brain, particularly the NLR family pyrin domain containing 3 (NLRP3), is currently targeted for AD treatment. Based on the crystal structure of the NACHT domain of NLRP3 and its renowned inhibitor MCC950, we designed and synthesized nineteen sulfonylurea compounds and evaluated their capacity to inhibit caspase-1 and interleukin-1ß (IL-1ß). Of these, nine were selected for measuring their IC50 for caspase-1 and cytotoxicity analysis. Finally, three compounds were chosen to assess their inhibitory effect on IL-1ß in mice. The results showed that compound 5m had a superior ability to reduce IL-1ß levels in the brain compared to MCC950 at a lower dosing concentration, indicating that 5m has the potential to penetrate the blood-brain barrier (BBB) and inhibit inflammation both in vitro and in vivo. Docking studies of compound 5m on NLRP3 revealed a binding mode similar to MCC950. These findings suggest that compound 5m holds promise as an NLRP3 inhibitor for AD treatment.

FitScore: a fast machine learning-based score for 3D virtual screening enrichment.

Enhancing virtual screening enrichment has become an urgent problem in computational chemistry, driven by increasingly large databases of commercially available compounds, without a commensurate drop in in vitro screening costs. Docking these large databases is possible with cloud-scale computing. However, rapid docking necessitates compromises in scoring, often leading to poor enrichment and an abundance of false positives in docking results. This work describes a new scoring function composed of two parts - a knowledge-based component that predicts the probability of a particular atom type being in a particular receptor environment, and a tunable weight matrix that converts the probability predictions into a dimensionless score suitable for virtual screening enrichment. This score, the FitScore, represents the compatibility between the ligand and the binding site and is capable of a high degree of enrichment across standardized docking test sets.

Structure-based design and optimization of a new class of small molecule inhibitors targeting the P-stalk binding pocket of ricin.

Ricin, a category-B agent for bioterrorism, and Shiga toxins (Stxs), which cause food poisoning bind to the ribosomal P-stalk to depurinate the sarcin/ricin loop. No effective therapy exists for ricin or Stx intoxication. Ribosome binding sites of the toxins have not been targeted by small molecules. We previously identified CC10501, which inhibits toxin activity by binding the P-stalk pocket of ricin toxin A subunit (RTA) remote from the catalytic site. Here, we developed a fluorescence polarization assay and identified a new class of compounds, which bind P-stalk pocket of RTA with higher affinity and inhibit catalytic activity with submicromolar potency. A lead compound, RU-NT-206, bound P-stalk pocket of RTA with similar affinity as a five-fold larger P-stalk peptide and protected cells against ricin and Stx2 holotoxins for the first time. These results validate the P-stalk binding site of RTA as a critical target for allosteric inhibition of the active site.

A Phenotypic Approach to the Discovery of Potent G-Quadruplex Targeted Drugs.

G-quadruplex (G4) sequences, which can fold into higher-order G4 structures, are abundant in the human genome and are over-represented in the promoter regions of many genes involved in human cancer initiation, progression, and metastasis. They are plausible targets for G4-binding small molecules, which would, in the case of promoter G4s, result in the transcriptional downregulation of these genes. However, structural information is currently available on only a very small number of G4s and their ligand complexes. This limitation, coupled with the currently restricted information on the G4-containing genes involved in most complex human cancers, has led to the development of a phenotypic-led approach to G4 ligand drug discovery. This approach was illustrated by the discovery of several generations of tri- and tetra-substituted naphthalene diimide (ND) ligands that were found to show potent growth inhibition in pancreatic cancer cell lines and are active in in vivo models for this hard-to-treat disease. The cycles of discovery have culminated in a highly potent tetra-substituted ND derivative, QN-302, which is currently being evaluated in a Phase 1 clinical trial. The major genes whose expression has been down-regulated by QN-302 are presented here: all contain G4 propensity and have been found to be up-regulated in human pancreatic cancer. Some of these genes are also upregulated in other human cancers, supporting the hypothesis that QN-302 is a pan-G4 drug of potential utility beyond pancreatic cancer.

Structure-Based Design of Bicyclic Helical Peptides That Target the Oncogene ß-Catenin.

The inhibition of intracellular protein-protein interactions is challenging, in particular, when involved interfaces lack pronounced cavities. The transcriptional co-activator protein and oncogene ß­catenin is a prime example of such a challenging target. Despite extensive targeting efforts, available high-affinity binders comprise only large molecular weight Inhibitors. This hampers the further development of therapeutically useful compounds. Herein, we report the design of a considerably smaller peptidomimetic scaffold derived from the α-helical ß­catenin-binding motif of Axin. Sequence maturation and bicyclization provided a stitched peptide with an unprecedented crosslink architecture. The binding mode and site were confirmed by a crystal structure. Further derivatization yielded a ß-catenin inhibitor with single-digit micromolar activity in a cell-based assay. This study sheds a light on how to design helix mimetics with reduced molecular weight thereby improving their biological activity.

Structural-guided design to improve the catalytic performance of aldo-keto reductase KdAKR.

Aldo-keto reductases (AKRs) are important biocatalysts that can be used to synthesize chiral pharmaceutical alcohols. In this study, the catalytic activity and stereoselectivity of a NADPH-dependent AKR from Kluyveromyces dobzhanskii (KdAKR) toward t-butyl 6-chloro (5S)-hydroxy-3-oxohexanoate ((5S)-CHOH) were improved by mutating its residues in the loop regions around the substrate-binding pocket. And the thermostability of KdAKR was improved by a consensus sequence method targeted on the flexible regions. The best mutant M6 (Y28A/L58I/I63L/G223P/Y296W/W297H) exhibited a 67-fold higher catalytic efficiency compared to the wild-type (WT) KdAKR, and improved R-selectivity toward (5S)-CHOH (dep value from 47.6% to >99.5%). Moreover, M6 exhibited a 6.3-fold increase in half-life (t1/2 ) at 40°C compared to WT. Under the optimal conditions, M6 completely converted 200 g/L (5S)-CHOH to diastereomeric pure t-butyl 6-chloro-(3R, 5S)-dihydroxyhexanoate ((3R, 5S)-CDHH) within 8.0 h, with a space-time yield of 300.7 g/L/day. Our results deepen the understandings of the structure-function relationship of AKRs, providing a certain guidance for the modification of other AKRs.

Evaluating the druggability of TrmD, a potential antibacterial target, through design and microbiological profiling of a series of potent TrmD inhibitors.

The post-transcriptional modifier tRNA-(N1G37) methyltransferase (TrmD) has been proposed to be essential for growth in many Gram-negative and Gram-positive pathogens, however previously reported inhibitors show only weak antibacterial activity. In this work, optimisation of fragment hits resulted in compounds with low nanomolar TrmD inhibition incorporating features designed to enhance bacterial permeability and covering a range of physicochemical space. The resulting lack of significant antibacterial activity suggests that whilst TrmD is highly ligandable, its essentiality and druggability are called into question.

Peptide-based covalent inhibitors of protein-protein interactions.

Protein-protein interactions (PPI) are involved in all cellular processes and many represent attractive therapeutic targets. However, the frequently rather flat and large interaction areas render the identification of small molecular PPI inhibitors very challenging. As an alternative, peptide interaction motifs derived from a PPI interface can serve as starting points for the development of inhibitors. However, certain proteins remain challenging targets when applying inhibitors with a competitive mode of action. For that reason, peptide-based ligands with an irreversible binding mode have gained attention in recent years. This review summarizes examples of covalent inhibitors that employ peptidic binders and have been tested in a biological context.

Design and synthesis of 1H-pyrazolo[3,4-d]pyrimidine derivatives as hematopoietic progenitor kinase 1 (HPK1) inhibitors.

Despite immune checkpoint inhibitors' tremendous success in the treatment of tumors, the moderate response rate limits their widespread use. Hematopoietic progenitor kinase 1 (HPK1) is served as an essential negative regulator of T-cell receptor, which has been identified as a promising target for enhancing antitumor immunity. However, the development of a selective HPK1 inhibitor is still challenging. Herein, we reported a novel series of 1H-pyrazolo[3,4-d]pyrimidine derivatives as HPK1 inhibitors by structure-based rational design. The optimal compound 10n significantly inhibited HPK1 with an IC50 value of 29.0 nM and the phosphorylation of SLP76 at a concentration as low as 0.1 µM. Furthermore, compound 10n exhibited good selectivity over a panel of 25 kinases, including GLK from the same MAP4K family. Together, the current study provided a novel, potent, and selective HPK1 inhibitor, acting as a lead compound for the future development of cancer immunotherapy.

Homology modeling of Forkhead box protein C2: identification of potential inhibitors using ligand and structure-based virtual screening.

Overexpression of Forkhead box protein C2 (FOXC2) has been associated with different types of carcinomas. FOXC2 plays an important role in the initiation and maintenance of the epithelial-mesenchymal transition (EMT) process, which is essential for the development of higher-grade tumors with an enhanced ability for metastasis. Thus, FOXC2 has become a therapeutic target for the development of anticancer drugs. MC-1-F2, the only identified experimental inhibitor of FOXC2, interacts with the full length of FOXC2. However, only the DNA-binding domain (DBD) of FOXC2 has resolved crystal structure. In this work, a three-dimensional (3D) structure of the full-length FOXC2 using homology modeling was developed and used for structure-based drug design (SBDD). The quality of this 3D model of the full-length FOXC2 was evaluated using MolProbity, ERRAT, and ProSA modules. Molecular dynamics (MD) simulation was also carried out to verify its stability. Ligand-based drug design (LBDD) was carried out to identify similar analogues for MC-1-F2 against 15 million compounds from ChEMBL and ZINC databases. 792 molecules were retrieved from this similarity search. De novo SBDD was performed against the full-length 3D structure of FOXC2 through homology modeling to identify novel inhibitors. The combination of LBDD and SBDD helped in gaining a better insight into the binding of MC-1-F2 and its analogues against the full length of the FOXC2. The binding free energy of the top hits was further investigated using MD simulations and MM/GBSA calculations to result in eight promising hits as lead compounds targeting FOXC2.

Computational 3D Modeling-Based Identification of Inhibitors Targeting Cysteine Covalent Bond Catalysts for JAK3 and CYP3A4 Enzymes in the Treatment of Rheumatoid Arthritis.

This work aimed to find new inhibitors of the CYP3A4 and JAK3 enzymes, which are significant players in autoimmune diseases such as rheumatoid arthritis. Advanced computer-aided drug design techniques, such as pharmacophore and 3D-QSAR modeling, were used. Two strong 3D-QSAR models were created, and their predictive power was validated by the strong correlation (R2 values > 80%) between the predicted and experimental activity. With an ROC value of 0.9, a pharmacophore model grounded in the DHRRR hypothesis likewise demonstrated strong predictive ability. Eight possible inhibitors were found, and six new inhibitors were designed in silico using these computational models. The pharmacokinetic and safety characteristics of these candidates were thoroughly assessed. The possible interactions between the inhibitors and the target enzymes were made clear via molecular docking. Furthermore, MM/GBSA computations and molecular dynamics simulations offered insightful information about the stability of the binding between inhibitors and CYP3A4 or JAK3. Through the integration of various computational approaches, this study successfully identified potential inhibitor candidates for additional investigation and efficiently screened compounds. The findings contribute to our knowledge of enzyme-inhibitor interactions and may help us create more effective treatments for autoimmune conditions like rheumatoid arthritis.

A molecular evolution algorithm for ligand design in DOCK.

As a complement to virtual screening, de novo design of small molecules is an alternative approach for identifying potential drug candidates. Here, we present a new 3D genetic algorithm to evolve molecules through breeding, mutation, fitness pressure, and selection. The method, termed DOCK_GA, builds upon and leverages powerful sampling, scoring, and searching routines previously implemented into DOCK6. Three primary experiments were used during development: Single-molecule evolution evaluated three selection methods (elitism, tournament, and roulette), in four clinically relevant systems, in terms of mutation type and crossover success, chemical properties, ensemble diversity, and fitness convergence, among others. Large scale benchmarking assessed performance across 651 different protein-ligand systems. Ensemble-based evolution demonstrated using multiple inhibitors simultaneously to seed growth in a SARS-CoV-2 target. Key takeaways include: (1) The algorithm is robust as demonstrated by the successful evolution of molecules across a large diverse dataset. (2) Users have flexibility with regards to parent input, selection method, fitness function, and molecular descriptors. (3) The program is straightforward to run and only requires a single executable and input file at run-time. (4) The elitism selection method yields more tightly clustered molecules in terms of 2D/3D similarity, with more favorable fitness, followed by tournament and roulette.

Targeting non-alcoholic fatty liver disease: Design, X-ray co-crystal structure and synthesis of 'first-in-kind' inhibitors of serine/threonine kinase25.

We describe the synthesis of a series of 3-t-butyl 5-aminopyrazole p-substituted arylamides as inhibitors of serine-threonine25 (STK25), an enzyme implicated in the progression of non-alcoholic fatty liver disease (NAFLD). Appending a p-N-pyrrolidinosulphonamide group to the arylamide group led to a 'first-in kind' inhibitor with IC50 = 228 nM. A co-crystal structure with STK 25 revealed productive interactions which were also reproduced using molecular docking. A new series of triazolo dihydro oxazine carboxamides of 3-t-butyl 5-aminopyrazole was not active against STK25.

Structure-based optimization of hydroxylactam as potent, cell-active inhibitors of lactate dehydrogenase.

Structure-based design was utilized to optimize 6,6-diaryl substituted dihydropyrone and hydroxylactam to obtain inhibitors of lactate dehydrogenase (LDH) with low nanomolar biochemical and single-digit micromolar cellular potencies. Surprisingly the replacement of a phenyl with a pyridyl moiety in the chemical structure revealed a new binding mode for the inhibitors with subtle conformational change of the LDHA active site. This led to the identification of a potent, cell-active hydroxylactam inhibitor exhibiting an in vivo pharmacokinetic profile suitable for mouse tumor xenograft study.

Targeting the interaction of ß-catenin and TCF/LEF transcription factors to inhibit oncogenic Wnt signaling.

The Wnt/ß-catenin signaling pathway is crucially involved in embryonic development, stem cell maintenance and tissue renewal. Hyperactivation of this pathway is associated with the development and progression of various types of cancers. The transcriptional coactivator ß-catenin represents a pivotal component of the pathway and its interaction with transcription factors of the TCF/LEF family is central to pathway activation. Inhibition of this crucial protein-protein interaction via direct targeting of ß-catenin is considered a promising strategy for the inactivation of oncogenic Wnt signaling. This review summarizes advances in the development of Wnt antagonists that have been shown to directly bind ß-catenin.

Computational design, synthesis and biological evaluation of PDE5 inhibitors based on N2,N4-diaminoquinazoline and N2,N6-diaminopurine scaffolds.

We report the synthesis, and characterization of twenty-nine new inhibitors of PDE5. Structure-based design was employed to modify to our previously reported 2,4-diaminoquinazoline series. Modification include scaffold hopping to 2,6-diaminopurine core as well as incorporation of ionizable groups to improve both activity and solubility. The prospective binding mode of the compounds was determined using 3D ligand-based similarity methods to inhibitors of known binding mode, combined with a PDE5 docking and molecular dynamics based-protocol, each of which pointed to the same binding mode. Chemical modifications were then designed to both increase potency and solubility as well as validate the binding mode prediction. Compounds containing a quinazoline core displayed IC50s ranging from 0.10 to 9.39 µM while those consisting of a purine scaffold ranging from 0.29 to 43.16 µM. We identified 25 with a PDE5 IC50 of 0.15 µM, and much improved solubility (1.77 mg/mL) over the starting lead. Furthermore, it was found that the predicted binding mode was consistent with the observed SAR validating our computationally driven approach.

Structure-Based design of [(2-Hydroxyethoxy)methyl]-6-(phenylthio)-thymine derivatives as nonnucleoside HIV-1 reverse transcriptase Inhibitors: From HEPTs to Sulfinyl-substituted HEPTs.

The [(2-hydroxyethoxy)methyl]-6-(phenylthio)thymine (HEPT) analogs were reported to be a kind of promising lead compounds as nonnucleoside HIV-1 reverse transcriptase inhibitors. In this work, a series of novel sulfinyl-substituted analogs were designed by structure-based design strategy with the purpose of improving the activity of HEPT, followed by evaluating their anti-HIV-1 activity in MT-4 cells. Most of the final compounds had moderate to strong activity against wild-type HIV-1 strain (IIIB) with EC50 values in the range of 0.21-1.91 µM, which were around 4 âˆ¼ 32-fold better than the reference compound HEPT. Some of them showed higher sensitivity toward clinically relevant mutant L100I and E138K viruses than NVP. Selected compounds were further evaluated for their activity against wild-type reverse transcriptase (RT), and most of them exhibited nanomolar activity, suggesting a good correlation with the cell-based activity. The compounds 11h, 11l, and 11ab displayed the best anti-HIV-1 activity against wild-type HIV-1 strain (EC50 = 0.280, 0.209, and 0.290 µM) and nanomolar activity against mutant strains (L100I and E138K), superior to HEPT and NVP. Molecular modeling studies were also performed to elucidate the biological activity, providing a structural insight for follow-up research on HEPT optimization.

Smaller, Stronger, More Stable: Peptide Variants of a SARS-CoV-2 Neutralizing Miniprotein.

Based on the structure of a de novo designed miniprotein (LCB1) in complex with the receptor binding domain (RBD) of the SARS-CoV-2 spike protein, we have generated and characterized truncated peptide variants of LCB1, which present only two of the three LCB1 helices, and which fully retained the virus neutralizing potency against different SARS-CoV-2 variants of concern (VOC). This antiviral activity was even 10-fold stronger for a cyclic variant of the two-helix peptides, as compared to the full-length peptide. Furthermore, the proteolytic stability of the cyclic peptide was substantially improved, rendering it a better potential candidate for SARS-CoV-2 therapy. In a more mechanistic approach, the peptides also served as tools to dissect the role of individual mutations in the RBD for the susceptibility of the resulting virus variants to neutralization by the peptides. As the peptides reported here were generated through chemical synthesis, rather than recombinant protein expression, they are amenable to further chemical modification, including the incorporation of a wide range of non-proteinogenic amino acids, with the aim to further stabilize the peptides against proteolytic degradation, as well as to improve the strength, as well the breadth, of their virus neutralizing capacity.