Development of CDK4/6 Inhibitors: A Five Years Update

The inhibition of cyclin dependent kinases 4 and 6 plays a role in aromatase inhibitor resistant metastatic breast cancer. Three dual CDK4/6 inhibitors have been approved for the breast cancer treatment that, in combination with the endocrine therapy, dramatically improved the survival outcomes both in first and later line settings. The developments of the last five years in the search for new selective CDK4/6 inhibitors with increased selectivity, treatment efficacy, and reduced adverse effects are reviewed, considering the small-molecule inhibitors and proteolysis-targeting chimeras (PROTACs) approaches, mainly pointing at structure-activity relationships, selectivity against different kinases and antiproliferative activity.     

Pharmacophore and docking-based virtual screening approach for the design of new dual inhibitors of Janus kinase 1 and Janus kinase 2

Janus kinase 1 and 2, non-receptor protein tyrosine kinases, are implicated in various cancerous diseases. Involvement of these two enzymes in the pathways that stimulate cell proliferation in cancerous conditions makes them potential therapeutic targets for designing new dual-targeted agents for the treatment of cancer. In the present study, two separate pharmacophore models were developed and the best models for JAK1 (AAADH.25) and JAK2 (ADRR.92) were selected on the basis of their external predictive ability. Both models were subjected to a systematic virtual screening (VS) protocol using a PHASE database of 1.5 million molecules. The hits retrieved in VS were investigated for ADME properties to avoid selection of molecules with a poor pharmacokinetic profile. The molecules considered to be within the range of acceptable limits of ADME properties were further employed for docking simulations with JAK1 and JAK2 proteins to explore the final hits that possess structural features of both pharmacophore models as well as display essential interactions with both of them. Thus, the new molecules obtained in this way should show inhibitory activity against JAK1 and JAK2 and may serve as novel therapeutic agents for the treatment of cancerous disease conditions.     

Molecular modeling study of checkpoint kinase 1 inhibitors by multiple docking strategies and prime/MM-GBSA calculation

Developing chemicals that inhibit checkpoint kinase 1 (Chk1) is a promising adjuvant therapeutic to improve the efficacy and selectivity of DNA-targeting agents. Reliable prediction of binding-free energy and binding affinity of Chk1 inhibitors can provide a guide for rational drug design. In this study, multiple docking strategies and Prime/Molecular Mechanics Generalized Born Surface Area (Prime/MM-GBSA) calculation were applied to predict the binding mode and free energy for a series of benzoisoquinolinones as Chk1 inhibitors. Reliable docking results were obtained using induced-fit docking and quantum mechanics/molecular mechanics (QM/MM) docking, which showed superior performance on both ligand binding pose and docking score accuracy to the rigid-receptor docking. Then, the Prime/MM-GBSA method based on the docking complex was used to predict the binding-free energy. The combined use of QM/MM docking and Prime/MM-GBSA method could give a high correlation between the predicted binding-free energy and experimentally determined pIC(50) . The molecular docking combined with Prime/MM-GBSA simulation can not only be used to rapidly and accurately predict the binding-free energy of novel Chk1 inhibitors but also provide a novel strategy for lead discovery and optimization targeting Chk1.     

Identification of novel Plasmodium falciparum PI4KB inhibitors as potential anti-malarial drugs: Homology modeling,molecular docking and molecular dynamics simulations

The current study was set to discover selective Plasmodium falciparum phosphatidylinositol-4-OH kinase type III beta (pfPI4KB) inhibitors as potential antimalarial agents using combined structure-based and ligand-based drug discovery approach. A comparative model of pfPI4KB was first constructed and validated using molecular docking techniques. Performance of Autodock4.2 and Vina4 software in predicting the inhibitor-PI4KB binding mode and energy was assessed based on two Test Sets: Test Set I contained five ligands with resolved crystal structures with PI4KB, while Test Set II considered eleven compounds with known IC50 value towards PI4KB. The outperformance of Autodock as compared to Vina was reported, giving a correlation coefficient (R²) value of 0.87 and 0.90 for Test Set I and Test Set II, respectively. Pharmacophore-based screening was then conducted to identify drug-like molecules from ZINC database with physicochemical similarity to two potent pfPI4KB inhibitors –namely cpa and cpb. For each query inhibitor, the best 1000 hits in terms of TanimotoCombo scores were selected and subjected to molecular docking and molecular dynamics (MD) calculations. Binding energy was then estimated using molecular mechanics–generalized Born surface area (MM-GBSA) approach over 50 ns MD simulations of the inhibitor-pfPI4KB complexes. According to the calculated MM-GBSA binding energies, ZINC78988474 and ZINC20564116 were identified as potent pfPI4KB inhibitors with binding energies better than those of cpa and cpb, with ΔG_binding ≥ −34.56 kcal/mol. The inhibitor-pfPI4KB interaction and stability were examined over 50 ns MD simulation; as well the selectivity of the identified inhibitors towards pfPI4KB over PI4KB was reported.     

Synthesis and biological evaluation of substituted 2-anilino-7H-pyrrolopyrimidines as PDK1 inhibitors

An efficient and scalable route for a series of novel substituted 2-anilino-7H-pyrrolopyrimidine compounds as potential inhibitors of PDK1, an important regulator of the PI3K/Akt pathway that is dysregulated in many cancers, was developed and is described. The synthetic strategy was designed around Suzuki and Buchwald–Hartwig cross-couplings of a boronate fragment and various customised anilines sequentially with 2,4-dichloro-7-tosyl-7H-pyrrolopyrimidine. All fragments were constructed separately and cross-coupled to provide access to a range of novel compounds. Biological evaluation of these was undertaken, with modest inhibition observed.     

In silico design novel (5-imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine derivatives as inhibitors for glycogen synthase kinase 3 based on 3D-QSAR,molecular docking and molecular dynamics simulation

Glycogen Synthase Kinase 3 (GSK-3) is a member of cellular kinase with various functions, such as glucose regulation, cellular differentiation, neuronal function and cell apoptosis. It has been proved as an important therapeutic target in type 2 diabetes mellitus and Alzheimer's disease. To better understand their structure–activity relationships and mechanism of action, an integrated computational study, including three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, and molecular dynamics (MD), was performed on 79 (5-Imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine GSK-3 inhibitors. In this paper, we constructed 3D-QSAR using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) method. The results showed that the CoMFA model (q ² = 0.743, r² = 0.980) and the CoMSIA model (q² = 0.813, r² = 0.976) had stable and reliable predictive ability. The electrostatic and H-bond donor fields play important roles in the models. The contour maps of the model visually showed the relationship between the activity of compounds and their three-dimensional structure. Molecular docking was used to identify the key amino acid residues at the active site of GSK-3 and explore its binding mode with ligands. Based on 3D-QSAR models, contour maps and the binding feature between GSK-3 and inhibitor, we designed 10 novel compounds with good potential activity and ADME/T profile. Molecular dynamics simulation results validated that Ile62, Val70 and Lys85 located in the active site play a key role for GSK-3 complexed with inhibitors. These results might provide important information for designing GSK-3 inhibitors with high activity.     

The design and synthesis of novel 3-[2-indol-1-yl-ethyl]-1H-indole derivatives as selective inhibitors of CDK4

We present the design, synthesis and biological activity of novel 3-[2-indol-1-yl-ethyl]-1H-indole selective inhibitors of CDK4.     

Design of selective TACE inhibitors using molecular docking studies: Synthesis and preliminary evaluation of anti-inflammatory and TACE inhibitory activity

Tumor necrosis factor-α (TNF-α) converting enzyme (TACE) has been considered one of the principal therapeutic targets for the treatment of TNF-dependent pathologies. Several TACE inhibitors have been reported, but none of them has been successfully passed to phase II clinical trials. In the present work, we attempted to design highly selective new non-hydroxamate sulfonamide TACE inhibitors. The docking study was performed on one of the crystal structures of TACE, selected based on its resolution and R value, to tackle the flexibility issue of the active site. The results allowed us to distinguish the analogues with a higher binding affinity toward the active site of TACE and to identify the substituent of analogues needed for binding with the surrounding site of the enzyme. Finally the analogues were docked on crystal structures of six different matrix metalloproteinases (MMPs) for a selectivity study of TACE over MMPs. Some of these analogues were synthesized and subjected to preliminary testing for in vivo anti-inflammatory activity and TACE inhibitory activity.     

Discovery of potential inhibitors targeting the kinase domain of polynucleotide kinase/phosphatase (PNKP): Homology modeling,virtual screening based on multiple conformations,and molecular dynamics simulation

In recent years, the level of interest has been increased in developing the DNA-repair inhibitors, to enhance the cytotoxic effects in the treatment of cancers. Polynucleotide kinase/phosphatase (PNKP) is a critical human DNA repair enzyme that repairs DNA strand breaks by catalyzing the restoration of 5’-phosphate and 3’-hydroxyl termini that are required for subsequent processing by DNA ligases and polymerases. PNKP is the only protein that repairs the 3′-hydroxyl group and 5′-phosphate group, which depicts PNKP as a potential therapeutic target. Besides, PNKP is the only DNA-repair enzyme that contains the 5′-kinase activity, therefore, targeting this kinase domain would motivate the development of novel PNKP-specific inhibitors. However, there are neither crystal structures of human PNKP nor the kinase inhibitors reported so far. Thus, in this present study, a sequential molecular docking-based virtual screening with multiple PNKP conformations integrating homology modeling, molecular dynamics simulation, and binding free energy calculation was developed to discover novel PNKP kinase inhibitors, and the top-scored molecule was finally submitted to molecular dynamics simulation to reveal the binding mechanism between the inhibitor and PNKP. Taken together, the current study could provide some guidance for the molecular docking based-virtual screening of novel PNKP kinase inhibitors.     

Tin(IV) complex of 3,4-dihydroxybenzaldehyde 4-ethylthiosemicarbazone: Synthesis,spectral studies and molecular modeling

The tin(IV) complex with 3,4-dihydroxybenzaldehyde 4-ethylthiosemicarbazone, [Sn(BETSC)Cl₂] has been synthesised and characterized using elemental analysis, Fourier transform infrared spectroscopy (FT-IR), UV-vis and NMR spectra. Theoretical calculations have also been performed by B3LYP method using LANL2DZ basis set. The assignments of bands observed in FT-IR spectrum for H₂BETSC and its Sn(IV) complex have been made using DFT method. Some significant differences in vibrational structures between the H₂BETSC and [Sn(BETSC)Cl₂] have been observed and discussed. Molecular docking was performed for the ligand H₂BETSC against dihydrofolate reductase (DHFR).     

Virtual screening of B-Raf kinase inhibitors: A combination of pharmacophore modelling,molecular docking, 3D-QSAR model and binding free energy calculation studies

B-Raf kinase has been identified as an important target in recent cancer treatment. In order to discover structurally diverse and novel B-Raf inhibitors (BRIs), a virtual screening of BRIs against ZINC database was performed by using a combination of pharmacophore modelling, molecular docking, 3D-QSAR model and binding free energy (ΔG_bind) calculation studies in this work. After the virtual screening, six promising hit compounds were obtained, which were then tested for inhibitory activities of A375 cell lines. In the result, five hit compounds show good biological activities (IC₅₀ < 50 μM). The present method of virtual screening can be applied to find structurally diverse inhibitors, and the obtained five structurally diverse compounds are expected to develop novel BRIs.     

A comparative QSAR analysis and molecular docking studies of phenyl piperidine derivatives as potent dual NK1R antagonists/serotonin transporter (SERT) inhibitors

Depression is a critical mood disorder that affects millions of patients. Available therapeutic antidepressant agents are associated with several undesirable side effects. Recently, it has been shown that Neurokinin 1 receptor (NK₁R) antagonists can potentiate the antidepressant effects of serotonin-selective reuptake inhibitors (SSRIs). In this study, a series of phenyl piperidine derivatives as potent dual NK₁R antagonists/serotonin transporter (SERT) inhibitors were applied to quantitative structure–activity relationship (QSAR) analysis. A collection of chemometrics methods such as multiple linear regression (MLR), factor analysis–based multiple linear regression (FA-MLR), principal component regression (PCR), and partial least squared combined with genetic algorithm for variable selection (GA-PLS) were applied to make relations between structural characteristics and NK₁R antagonism/SERT inhibitory of these compounds. The best multiple linear regression equation was obtained from GA-PLS and MLR for NK₁R and SERT, respectively. Based on the resulted model, an in silico-screening study was also conducted and new potent lead compounds based on new structural patterns were designed for both targets. Molecular docking studies of these compounds on both targets were also conducted and encouraging results were acquired. There was a good correlation between QSAR and docking results. The results obtained from validated docking studies indicate that the important amino acids inside the active site of the cavity that are responsible for essential interactions are Glu33, Asp395 and Arg26 for SERT and Ala30, Lys7, Asp31, Phe5 and Tyr82 for NK₁R receptors.     

Design,synthesis, crystal structures and biological evaluation of some 1,3-thiazolidin-4-ones as dual CDK2/EGFR potent inhibitors with potential apoptotic antiproliferative effects

A series of novel thiazolidine-4-one derivatives was synthesized by reacting 1,4-disubstituted hydrazine carbothioamides with diethyl azodicarboxylate. The structures were confirmed by spectroscopic data as well as single-crystal X-ray analyses. The antiproliferative activity of the synthesized compounds was investigated against four human cancer cell lines using an MTT assay. Compounds 5d, 5e, and 5f revealed the most potent antiproliferative activity with GI₅₀ values ranging from 0.70 µM to 1.20 µM, compared to doxorubicin GI₅₀ value = 1.10 µM. Compounds 5d, 5e, and 5f were further investigated for their inhibitory activities against CDK2 and EGFR as potential targets for their molecular mechanism. Compounds 5e and 5f have showed potent inhibitory activity to CDK2 enzyme with IC₅₀ values of 18 and 14 nM, which is more potent than the reference dinaciclib (IC₅₀ = 20 nM). Moreover, compounds 5e and 5f were the most potent EGFR inhibitors, with IC₅₀ values of 93 and 87 nM, respectively, compared to the reference erlotinib (IC₅₀ = 70 nM). In addition, the most potent derivatives were tested for their apoptotic activity against caspases 3, 8, and 9, and the results showed that compounds 5d, 5e, and 5f revealed a greater increase in active caspases 3,8 and 9 than doxorubicin. Also, compounds 5d, 5e, and 5f elevated cytochrome C levels in the MCF-7 human breast cancer cell line by about 15.5, 15.8, and 16.5 times, respectively. Finally, a molecular docking study was performed to investigate the binding sites of these compounds within the active sites of CDK2 and EGFR targets, and the results confirmed that the most potent CDK2 and EGFR inhibitor 5h also have showed the highest docking score.     

A merged molecular docking,ADME-T and dynamics approaches towards the genus of Arisaema as herpes simplex virus type 1 and type 2 inhibitors

An attempt toward screening of phytoconstituents (Arisaema genus) against herpes viruses (HSV-1 and HSV-2) was carried out using in silico approaches. Human HSV-1 and HSV-2 are accountable for cold sores genital herpes, respectively. Two drug targets, namely thymidine kinase (TK; PDB: 2ki5) serine protease (PDB: 1at3) were selected for HSV-1 and HSV-2. Initially, molecular docking tool was employed to screened apex hits phytoconstituents against herpes infections. ADME-T studies of top ranked were also further highlighted to achieve their effectiveness. Following, molecular dynamics studies were also examined to further optimize the stability of ligands. Glide scores and binding interactions of phytoconstituents were compared with Acyclovir, the main drug used in treatment of HSV, the screened top hits exhibited more glide scores and better binding for both HSV-1 and HSV-2 receptors. Additionally, ADME-T showed an ideal range for top hits while molecular dynamics results also illustrated stability of models. Ultimately, the whole efforts reveal to top three most promising hits for HSV-1 (39, 21, 19) and HSV-2 (20, 51, 19) receptors which can be explored further in wet lab experiments as promising agents against HSV infections.     

Synthesis of substituted 4-(1H-indol-6-yl)-1H-indazoles as potential PDK1 inhibitors

The development of a preparative route to a series of novel 4-(1H-indol-6-yl)-1H-indazole compounds as potential PDK1 inhibitors is described. The synthetic strategy centres on the late-stage Suzuki cross-coupling of N-unprotected indazole and indole fragments. The use of a monoligated palladium catalyst system was found to be highly beneficial in the cross-coupling reaction. The indazole and indole fragments were constructed by diazotisation/cyclisation and S_NAr/reductive cyclisation sequences, respectively.     

The selectivity and promiscuity of brain-neuroregenerative inhibitors between ROCK1 and ROCK2 isoforms: An integration of SB-QSSR modelling,QM/MM analysis and in vitro kinase assay

The Rho-associated kinases (ROCKs) have long been recognized as an attractive therapeutic target for various neurological diseases; selective inhibition of ROCK1 and ROCK2 isoforms would result in distinct biological effects on neurogenesis, neuroplasticity and neuroregeneration after brain surgery and traumatic brain injury. However, the discovery and design of isoform-selective inhibitors remain a great challenge due to the high conservation and similarity between the kinase domains of ROCK1 and ROCK2. Here, a structure-based quantitative structure–selectivity relationship (SB-QSSR) approach was used to correlate experimentally measured selectivity with the difference in inhibitor binding to the two kinase isoforms. The resulting regression models were examined rigorously through both internal cross-validation and external blind validation; a nonlinear predictor was found to have high fitting stability and strong generalization ability, which was then employed to perform virtual screening against a structurally diverse, drug-like compound library. Consequently, five and seven hits were identified as promising candidates of 1-o-2 and 2-o-1 selective inhibitors, respectively, from which seven purchasable compounds were tested in vitro using a standard kinase assay protocol to determine their inhibitory activity against and selectivity between ROCK1 and ROCK2. The structural basis, energetic property and biological implication underlying inhibitor selectivity and promiscuity were also investigated systematically using a hybrid quantum mechanics/molecular mechanics (QM/MM) scheme.     

Comparison of azacyclic urea A-98881 as HIV-1 protease inhibitor with cage dimeric N-benzyl 4-(4-methoxyphenyl)-1,4- dihydropyridine as representative of a novel class of HIV-1 protease inhibitors: A molecular modeling study

The functional groups of cage dimeric N-alkyl substituted 3,5-bis(hydroxymethyl)-4-(4-methoxyphenyl)-1,4-dihydropyridines are similar to those of cyclic and azacyclic ureas that are potent inhibitors of HIV-1 protease of the dihydroxyethylene- and hydroxyethylene type, respectively. In the following study the conformity of common functional groups is investigated concerning their orientation in space as well as in the enzyme HIV-1 protease. Starting from X-ray crystal data of the centrosymmetric cage dimeric N-benzyl derivative with ester groups, the derivative with hydroxymethylene groups was built and a systematic conformational search was performed for the conformationally important torsion angles considering electrostatic and van der Waals interactions. From the huge number of conformations those comprising centrosymmetrical and C2-symmetrical energy minima were selected and minimized. The three remaining conformers were fitted to the azacyclic urea A-98881 selected from the HIV-1 protease enzyme- inhibitor complex using the centroids of the corresponding aromatic residues and additionally by the field fit option of the Advanced CoMFA module of SYBYL. Interestingly, the energetically most favourable one, which, additionally, possesses C2-symmetry like the active site cavity of HIV-1 protease, showed the best fit. Comparing the electrostatic potential (EP) of the latter with the EP of A-98881 the aromatic residues show excellent accordance. Slight differences in the extent of the EP were found in the areas of the hydroxymethylene groups of the cage dimer and the single hydroxy group as well as the urea carbonyl group of A- 98881, respectively. In order to compare the binding possibilities to the enzyme HIV-1 protease for the cage dimer and A-98881, their interaction fields with certain probes (CH3 for alkyl, NHamide, and carbonyl, O– of COO–), representing the decisive functional groups of the active site, have been calculated using GRID and projected into the enzyme placing the structures according to the position of A-98881 in the enzyme- inhibitor complex. The strongest calculated fields of the O– probe were found near Asp 25 for both structures. Another respective conformity consists in the overlap of the fields for the NHamide probe near Ile 50 and 50 for the investigated cage dimer and A-98881.