Design, synthesis and SAR of thienopyridines as potent CHK1 inhibitors

Insight into crystal structures and identification of potential styrylthieno[2,3-b]pyridine-2-carboxamidederivatives against COVID-19 Mpro through structure-guided modeling and simulation approach

Anti-SARS-CoV-2 drugs are urgently needed to prevent the pandemic and for immunization. Their protease inhibitor treatment for COVID-19 has been used in clinical trials. In Calu-3 and THP1 cells, 3CL SARS-CoV-2 Mpro protease is required for viral expression, replication, and the activation of the cytokines IL-1, IL-6, and TNF-. The Mpro structure was chosen for this investigation because of its activity as a chymotrypsin-like enzyme and the presence of a cysteine-containing catalytic domain. Thienopyridine derivatives increase the release of nitric oxide from coronary endothelial cells, which is an important cell signaling molecule with antibacterial activity against bacteria, protozoa, and some viruses. Using DFT calculations, global descriptors are computed from HOMO-LUMO orbitals; the molecular reactivity sites are analyzed from an electrostatic potential map. NLO properties are calculated, and topological analysis is also part of the QTAIM studies. Both compounds 1 and 2 were designed from the precursor molecule pyrimidine and exhibited binding energies (-14.6708 kcal/mol and -16.4521 kcal/mol). The binding mechanisms of molecule 1 towards SARS-COV-2 3CL Mpro exhibited strong hydrogen bonding as well as Vdw interaction. In contrast, derivative 2 was bound to the active site protein's active studied that several residues and positions, including (His41, Cys44, Asp48, Met49, Pro52, Tyr54, Phe140, Leu141, Ser144, His163, Ser144, Cys145, His164, Met165, Glu166, Leu167, Asp187, Gln189, Thr190, and GLn192) are critical for the maintenance of inhibitors inside the active pocket. Molecular docking and 100 ns MD simulation analysis revealed that Both compounds 1 and 2 with higher binding affinity and stability toward the SARS-COV-2 3CL Mpro protein. Binding free energy calculations and other MD parameters support the finding.Communicated by Ramaswamy H. Sarma.

The pyridazine heterocycle in molecular recognition and drug discovery

The pyridazine ring is endowed with unique physicochemical properties, characterized by weak basicity, a high dipole moment that subtends π-π stacking interactions and robust, dual hydrogen-bonding capacity that can be of importance in drug-target interactions. These properties contribute to unique applications in molecular recognition while the inherent polarity, low cytochrome P450 inhibitory effects and potential to reduce interaction of a molecule with the cardiac hERG potassium channel add additional value in drug discovery and development. The recent approvals of the gonadotropin-releasing hormone receptor antagonist relugolix (24) and the allosteric tyrosine kinase 2 inhibitor deucravacitinib (25) represent the first examples of FDA-approved drugs that incorporate a pyridazine ring. In this review, the properties of the pyridazine ring are summarized in comparison to the other azines and its potential in drug discovery is illustrated through vignettes that explore applications that take advantage of the inherent physicochemical properties as an approach to solving challenges associated with candidate optimization.

Graphical Abstract

Design, synthesis and evaluation of novel thienopyridazine derivatives as Chk1/2 inhibitors

In order to search for novel checkpoint kinase 1/2 (Chk1) inhibitors, we have designed and synthesized a series of new compounds incorporating thienopyridazine core. Bioevaluation showed that compounds 10j, 10i, 13e and 10o exhibited relatively good inhibitory activity. Notably, compound 10o displayed high selectivity against a panel of kinases and inhibited Chk1/2 signaling pathway stimulated by DNA damage drugs in cellular level. Molecular docking of 10o to the ATP-binding site of Chk1 kinase domain indicated the existence of polar interactions between 10o and the ATP-ribose-binding residues of Chk1. In mouse HT-29 xenografts, a synergistic effect was observed. Co-treatment by CPT-11 and 10o significantly diminished the tumor volume, indicating the great potential of 10o as a candidate of Chk1/2 inhibitor.

Three-Dimensional Analysis of Binding Sites for Predicting Binding Affinities in Drug Design

Understanding the interaction between drug molecules and proteins is one of the main challenges in drug design. Several tools have been developed recently to decrease the complexity of the process. Artificial intelligence and machine learning methods offer promising results in predicting the binding affinities. It becomes possible to do accurate predictions by using the known protein-ligand interactions. In this study, the electrostatic potential values extracted from 3-dimensional grid cubes of the drug-protein binding sites are used for predicting binding affinities of related complexes. A new algorithm with a dynamic feature selection method was implemented, which is derived from Compressed Images For Affinity Prediction (CIFAP) study, to predict binding affinities of Checkpoint Kinase 1 and Caspase 3 inhibitors.

Applications of amide isosteres in medicinal chemistry

Isosteric replacement of amide groups is a classic practice in medicinal chemistry. This digest highlights the applications of most commonly employed amide isosteres in drug design aiming at improving potency and selectivity, optimizing physicochemical and pharmacokinetic properties, eliminating or modifying toxicophores, as well as providing novel intellectual property of lead compounds.

Acoustic Droplet Ejection Enabled Automated Reaction Scouting

Miniaturization and acceleration of synthetic chemistry are critically important for rapid property optimization in pharmaceutical, agrochemical, and materials research and development. However, in most laboratories organic synthesis is still performed on a slow, sequential, and material-consuming scale and not validated for multiple substrate combinations. Herein, we introduce fast and touchless acoustic droplet ejection (ADE) technology into small-molecule chemistry to transfer building blocks by nL droplets and to scout a newly designed isoquinoline synthesis. With each compound in a discrete well, 384 random derivatives were synthesized in an automated fashion, and their quality was monitored by SFC-MS and TLC-UV-MS analysis. We exemplify a pipeline of fast and efficient nmol scouting to mmol- and mol-scale synthesis for the discovery of a useful novel reaction with great scope.

Regiospecific alkyl addition of (hetero)arene-fused thiophenes enabled by a visible-light-mediated photocatalytic desulfuration approach

A unique photoredox desulfuration approach enabling the regiospecific alkyl addition of (hetero)arene-fused thiophenes is presented.

Thieno[2,3-b]pyridine derivatives: a new class of antiviral drugs against Mayaro virus

Mayaro virus (MAYV) is an arthropod-borne virus and a member of the family Togaviridae, genus Alphavirus. Its infection leads to an acute illness accompanied by long-lasting arthralgia. To date, there are no antiviral drugs or vaccines against infection with MAYV and resources for the prevention or treatment of other alphaviruses are very limited. MAYV has served as a model to study the antiviral potential of several substances on alphavirus replication. In this work we evaluated the antiviral effect of seven new derivatives of thieno[2,3-b]pyridine against MAYV replication in a mammalian cell line. All derivatives were able to reduce viral production effectively at concentrations that were non-toxic for Vero cells. Molecular modeling assays predicted low toxicity risk and good oral bioavailability of the substances in humans. One of the molecules, selected for further study, demonstrated a strong anti-MAYV effect at early stages of replication, as it protected pre-treated cells and also during the late stages, affecting virus morphogenesis. This study is the first to demonstrate the antiviral effect of thienopyridine derivatives on MAYV replication in vitro, suggesting the potential application of these substances as antiviral molecules against alphaviruses. Additional in vivo research will be needed to expand the putative therapeutic applications.

Ru-Catalysed synthesis of fused heterocycle-pyridinones and -pyrones

Heterocycle-pyridinones and heterocycle-pyranones have been prepared by Ru-catalysed oxidative coupling of N-unprotected primary heterocycle-amides and heterocycle-carboxylic acids with internal alkynes.

OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins

The parametrization and validation of the OPLS3 force field for small molecules and proteins are reported. Enhancements with respect to the previous version (OPLS2.1) include the addition of off-atom charge sites to represent halogen bonding and aryl nitrogen lone pairs as well as a complete refit of peptide dihedral parameters to better model the native structure of proteins. To adequately cover medicinal chemical space, OPLS3 employs over an order of magnitude more reference data and associated parameter types relative to other commonly used small molecule force fields (e.g., MMFF and OPLS_2005). As a consequence, OPLS3 achieves a high level of accuracy across performance benchmarks that assess small molecule conformational propensities and solvation. The newly fitted peptide dihedrals lead to significant improvements in the representation of secondary structure elements in simulated peptides and native structure stability over a number of proteins. Together, the improvements made to both the small molecule and protein force field lead to a high level of accuracy in predicting protein-ligand binding measured over a wide range of targets and ligands (less than 1 kcal/mol RMS error) representing a 30% improvement over earlier variants of the OPLS force field.

DockBench: An Integrated Informatic Platform Bridging the Gap between the Robust Validation of Docking Protocols and Virtual Screening Simulations

Virtual screening (VS) is a computational methodology that streamlines the drug discovery process by reducing costs and required resources through the in silico identification of potential drug candidates. Structure-based VS (SBVS) exploits knowledge about the three-dimensional (3D) structure of protein targets and uses the docking methodology as search engine for novel hits. The success of a SBVS campaign strongly depends upon the accuracy of the docking protocol used to select the candidates from large chemical libraries. The identification of suitable protocols is therefore a crucial step in the setup of SBVS experiments. Carrying out extensive benchmark studies, however, is usually a tangled task that requires users' proficiency in handling different file formats and philosophies at the basis of the plethora of existing software packages. We present here DockBench 1.0, a platform available free of charge that eases the pipeline by automating the entire procedure, from docking benchmark to VS setups. In its current implementation, DockBench 1.0 handles seven docking software packages and offers the possibility to test up to seventeen different protocols. The main features of our platform are presented here and the results of the benchmark study of human Checkpoint kinase 1 (hChk1) are discussed as validation test.

Computational studies on the binding mechanism between triazolone inhibitors and Chk1 by molecular docking and molecular dynamics

Chk1, a serine/threonine protein kinase that participates in transducing DNA damage signals, is an attractive target due to its involvement in tumor initiation and progression. As a novel Chk1 inhibitor, the triazolone's bioactivity mechanism is not clear. In this study, we carried out an integrated computational study that combines molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations to identify the key factors necessary for the bioactivities. With the aim of discerning the structural features that affect the inhibitory activity of triazolones, MK-8776, a Chk1 inhibitor that reached the clinical stage, was also used as a reference for simulations. A comparative analysis of the triazolone inhibitors at the molecular level offers valuable insight into the structural and energetic properties. A general feature is that all the studied inhibitors bind in the pocket characterized by residues Leu14, Val22, Ala35, Glu84, Tyr85, Cys86, and Leu136 of Chk1. Moreover, introducing hydrophobic groups into triazolone inhibitors is favorable for binding to Chk1, which is corroborated by residue Leu136 with a relatively large difference in the contribution between MK-8776 and five triazolones to the total binding free energies. A hydrogen bond between the polar hydrogen atoms at R1 and Cys86 can facilitate proper placement of the inhibitor in the binding pocket of Chk1 that favors binding. However, the introduction of hydrophilic groups into the R2 position diminishes binding affinity. The information provided by this research is of benefit for further rational design of novel promising inhibitors of Chk1.

Docking-based 3D-QSAR study of pyridyl aminothiazole derivatives as checkpoint kinase 1 inhibitors

Checkpoint kinase 1 (Chk1) is a promising target for the design of novel anticancer agents. In the present work, molecular docking simulations and three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were performed on pyridyl aminothiazole derivatives as Chk1 inhibitors. AutoDock was used to determine the probable binding conformations of all the compounds inside the active site of Chk1. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models were developed based on the docking conformations and alignments. The CoMFA model produced statistically significant results with a cross-validated correlation coefficient (q2) of 0.608 and a coefficient of determination (r2) of 0.972. The reliable CoMSIA model with q2 of 0.662 and r2 of 0.970 was obtained from the combination of steric, electrostatic and hydrogen bond acceptor fields. The predictive power of the models were assessed using an external test set of 14 compounds and showed reasonable external predictabilities (r(2)pred) of 0.668 and 0.641 for CoMFA and CoMSIA models, respectively. The models were further evaluated by leave-ten-out cross-validation, bootstrapping and progressive scrambling analyses. The study provides valuable information about the key structural elements that are required in the rational design of potential drug candidates of this class of Chk1 inhibitors.

Structure-based design, discovery and development of checkpoint kinase inhibitors as potential anticancer therapies

INTRODUCTION

Checkpoint kinase (CHK) inhibitors offer the promise of enhancing the effectiveness of widely prescribed cancer chemotherapies and radiotherapy by inhibiting the DNA damage response, as well as the potential for single agent efficacy.

AREAS COVERED

This article surveys structural insights into the checkpoint kinases CHK1 and CHK2 that have been exploited to enhance the selectivity and potency of small molecule inhibitors. Furthermore, the authors review the use of mechanistic cellular assays to guide the optimisation of inhibitors. Finally, the authors discuss the status of the current clinical candidates and emerging new clinical contexts for CHK1 and CHK2 inhibitors, including the prospects for single agent efficacy.

EXPERT OPINION

Protein-bound water molecules play key roles in structural features that can be targeted to gain high selectivity for either enzyme. The results of early phase clinical trials of checkpoint inhibitors have been mixed, but significant progress has been made in testing the combination of CHK1 inhibitors with genotoxic chemotherapy. Second-generation CHK1 inhibitors are likely to benefit from increased selectivity and oral bioavailability. While the optimum therapeutic context for CHK2 inhibition remains unclear, the emergence of single agent preclinical efficacy for CHK1 inhibitors in specific tumour types exhibiting constitutive replication stress represents exciting progress in exploring the therapeutic potential of these agents.

Discovery of checkpoint kinase inhibitor (S)-5-(3-fluorophenyl)-N-(piperidin-3-yl)-3-ureidothiophene-2-carboxamide (AZD7762) by structure-based design and optimization of thiophenecarboxamide ureas

Checkpoint kinases CHK1 and CHK2 are activated in response to DNA damage that results in cell cycle arrest, allowing sufficient time for DNA repair. Agents that lead to abrogation of such checkpoints have potential to increase the efficacy of such compounds as chemo- and radiotherapies. Thiophenecarboxamide ureas (TCUs) were identified as inhibitors of CHK1 by high throughput screening. A structure-based approach is described using crystal structures of JNK1 and CHK1 in complex with 1 and 2 and of the CHK1-3b complex. The ribose binding pocket of CHK1 was targeted to generate inhibitors with excellent cellular potency and selectivity over CDK1and IKKβ, key features lacking from the initial compounds. Optimization of 3b resulted in the identification of a regioisomeric 3-TCU lead 12a. Optimization of 12a led to the discovery of the clinical candidate 4 (AZD7762), which strongly potentiates the efficacy of a variety of DNA-damaging agents in preclinical models.