Discovery of novel wee1 inhibitors via structure-based virtual screening and biological evaluation

Wee1 plays a critical role in the arrest of G2/M cell cycle for DNA repair before entering mitosis. Many cancer cells have been identified as overexpression of Wee1. In this research, pharmacophore modeling, molecular docking and molecular dynamics simulation approaches were constructed to identify novel potential Wee1 inhibitors. A compound 8 was found to have a novel skeleton against Wee1 with an IC₅₀ value of 22.32 µM and a K_i value of 13.11 µM. Kinetic assays were employed to evaluate the compound 8 as a competitive inhibitor. Compound 8 was tested against A-549 tumor cell lines with IC₅₀ value of 17.8 µM. To investigate the intermolecular interaction of Wee1 and compound 8, further molecular dynamics simulations were performed. It indicates that the binding mode of compound 8 and reference ligand is similar. The active core scaffold of compound 8 could represent a promising lead compound for studying Wee1 and be used for further structural optimization to design more potent Wee1 inhibitors.     

Design,synthesis, and biological evaluation of diarylpyrazole derivatives as antitumor agents targeting microtubules

A new series of novel diarylpyrazole derivatives as microtubule destabilizers were synthesized and evaluated for the anti-proliferative activities. Anti-proliferative assays were performed on the human cervix adenocarcinoma cell line (HeLa) and human gastric adenocarcinoma cell line (SGC-7901), and the compound 9s containing indole ring showed great anti-proliferative activity against HeLa cells with IC₅₀ value of 1.9 ± 0.11 μM. Further biological studies showed that 9s was able to inhibit tubulin polymerization, disrupt the cytoskeleton, block the cell cycle in the G2/M phase, and induce cell apoptosis in a concentration-dependent manner. In addition, the results of molecular docking studies showed that compound 9s could bind tightly to the colchicine binding site of tubulin through hydrogen bonding interaction. These preliminary results recommend that compound 9s is likely to be a microtubule destabilizer that deserves further investigation.     

Investigation of the binding mode of a novel cruzain inhibitor by docking,molecular dynamics,ab initio and MM/PBSA calculations

Chagas disease remains a major health problem in South America, and throughout the world. The two drugs clinically available for its treatment have limited efficacy and cause serious adverse effects. Cruzain is an established therapeutic target of Trypanosoma cruzi, the protozoan that causes Chagas disease. Our group recently identified a competitive cruzain inhibitor (compound 1) with an IC₅₀?=?15 µM that is also more synthetically accessible than the previously reported lead, compound 2. Prior studies, however, did not propose a binding mode for compound 1, hindering understanding of the structure–activity relationship and optimization. Here, the cruzain binding mode of compound 1 was investigated using docking, molecular dynamics (MD) simulations with ab initio derived parameters, ab initio calculations, and MM/PBSA. Two ligand protonation states and four binding poses were evaluated. A careful ligand parameterization method was employed to derive more physically meaningful parameters than those obtained by automated tools. The poses of unprotonated 1 were unstable in MD, showing large conformational changes and diffusing away from the binding site, whereas the protonated form showed higher stability and interaction with negatively charged residues Asp161 and Cys25. MM/PBSA also suggested that these two residues contribute favorably to binding of compound 1. By combining results from MD, ab initio calculations, and MM/PBSA, a binding mode of 1 is proposed. The results also provide insights for further optimization of 1, an interesting lead compound for the development of new cruzain inhibitors.     

Molecular modeling studies to discover novel mIDH2 inhibitors with high selectivity for the primary and secondary mutants

Mutant isocitrate dehydrogenase 2 (mIDH2) is an emerging target for the treatment of cancer. AG-221 is the first mIDH2 inhibitor approved by the FDA for acute myeloid leukemia treatment, but its acquired resistance has recently been observed, necessitating the development of new inhibitor. In this study, a multi-step virtual screening protocol was employed for the analysis of a large database of compounds to identify potential mIDH2 inhibitors. To this end, we firstly utilized molecular dynamics (MD) simulations and binding free energy calculations to elucidate the key factors affecting ligand binding and drug resistance. Based on these findings, the receptor-ligand interaction-based pharmacophore (IBP) model and hierarchical docking-based virtual screening were sequentially carried out to assess 212,736 compounds from the Specs database. The resulting hits were finally ranked by PAINS filter and ADME prediction and the top compounds were obtained. Among them, six molecules were identified as mIDH2 putative inhibitors with high selectivity by interacting with the capping residue Asp312. Furthermore, subsequent docking and MD experiments demonstrated that compound V2 might have potential inhibitory activity against the AG-221-resistant mutants, thereby making it a promising lead for the development of novel mIDH2 inhibitors.     

Multicomponent access to novel proline/cyclized cysteine tethered monastrol conjugates as potential anticancer agents

The versatility of multicomponent Biginelli’s reaction is exploited in the development of proline and cyclized cysteine tethered conjugates of monastrol, a kinesin Eg5 inhibitor. Ten new conjugates are synthesized focusing on structural replacement of the ester moiety (C-5 position) of the monastrol backbone with amino acid based amide moieties. On cytotoxic evaluation, conjugate 24 has shown promising in vitro cytotoxic activity against leukemia. Molecular docking studies revealed that the conjugates 19 and 24 exhibit better interaction at kinesin Eg5 receptor compared to monastrol. Moreover, computational calculations and predictions of important molecular properties suggest that these new amino acid based conjugates could be further improved to provide potential anticancer agents.     

Indole-derived chalcones as anti-dermatophyte agents: In vitro evaluation and in silico study

A series of indole-derived methoxylated chalcones were described as anti-dermatophyte agents. The in vitro antifungal susceptibility testing against different dermatophytes revealed that most of compounds had potent activity against the dermatophyte strains. In particular, the 4-ethoxy derivative 4d with MIC values of 0.25−2 μg/ml was the most potent compound against Trichophyton interdigitale, Trichophyton veruccosum and Microsporum fulvum. Moreover, the 4-butoxy analog 4i displaying MIC values in the range of 1−16 μg/ml had the highest inhibitory activity against Trichophyton mentagrophytes, Microsporum canis, and Arthroderma benhamiae. To predict whether the synthesized compounds interact with tubulin binding site of dermatophytes, the 3D-structure of target protein was modeled by homology modeling and then used for molecular docking and molecular dynamics (MD) simulation studies. Docking simulation revealed that the promising compound 4d can properly bind with tubulin. The molecular dynamics analysis showed that interactions of compound 4d with the active site of target protein have binding stability throughout MD simulation. The results of this study could utilize in the design of more effective antifungal drugs with tubulin inhibition mechanism against keratinophilic fungi.     

Fluorogenic Probes/Inhibitors of β‐Lactamase and their Applications in Drug‐Resistant Bacteria

β‐Lactam antibiotics are generally perceived as one of the greatest inventions of the 20^th century, and these small molecular compounds have saved millions of lives. However, upon clinical application of antibiotics, the β‐lactamase secreted by pathogenic bacteria can lead to the gradual development of drug resistance. β‐Lactamase is a hydrolase that can efficiently hydrolyze and destroy β‐lactam antibiotics. It develops and spreads rapidly in pathogens, and the drug‐resistant bacteria pose a severe threat to human health and development. As a result, detecting and inhibiting the activities of β‐lactamase are of great value for the rational use of antibiotics and the treatment of infectious diseases. At present, many specific detection methods and inhibitors of β‐lactamase have been developed and applied in clinical practice. In this Minireview, we describe the resistance mechanism of bacteria producing β‐lactamase and further summarize the fluorogenic probes, inhibitors of β‐lactamase, and their applications in the treatment of infectious diseases. It may be valuable to design fluorogenic probes with improved selectivity, sensitivity, and effectiveness to further identify the inhibitors for β‐lactamases and eventually overcome bacterial resistance.     

Structure-based design and in-parallel synthesis of inhibitors of AmpC β-lactamase

Background: Group I β-lactamases are a major cause of antibiotic resistance to β-lactams such as penicillins and cephalosporins. These enzymes are only modestly affected by classic β-lactam-based inhibitors, such as clavulanic acid. Conversely, small arylboronic acids inhibit these enzymes at sub-micromolar concentrations. Structural studies suggest these inhibitors bind to a well-defined cleft in the group I β-lactamase AmpC; this cleft binds the ubiquitous R1 side chain of β-lactams. Intriguingly, much of this cleft is left unoccupied by the small arylboronic acids.Results: To investigate if larger boronic acids might take advantage of this cleft, structure-guided in-parallel synthesis was used to explore new inhibitors of AmpC. Twenty-eight derivatives of the lead compound, 3-aminophenylboronic acid, led to an inhibitor with 80-fold better binding (2; K_i 83 nM). Molecular docking suggested orientations for this compound in the R1 cleft. Based on the docking results, 12 derivatives of 2 were synthesized, leading to inhibitors with K_i values of 60 nM and with improved solubility. Several of these inhibitors reversed the resistance of nosocomial Gram-positive bacteria, though they showed little activity against Gram-negative bacteria. The X-ray crystal structure of compound 2 in complex with AmpC was subsequently determined to 2.1 Å resolution. The placement of the proximal two-thirds of the inhibitor in the experimental structure corresponds with the docked structure, but a bond rotation leads to a distinctly different placement of the distal part of the inhibitor. In the experimental structure, the inhibitor interacts with conserved residues in the R1 cleft whose role in recognition has not been previously explored.Conclusions: Combining structure-based design with in-parallel synthesis allowed for the rapid exploration of inhibitor functionality in the R1 cleft of AmpC. The resulting inhibitors differ considerably from β-lactams but nevertheless inhibit the enzyme well. The crystal structure of 2 (K_i 83 nM) in complex with AmpC may guide exploration of a highly conserved, largely unexplored cleft, providing a template for further design against AmpC β-lactamase.     

Design,synthesis, biological evaluation and molecular dynamics of LAR inhibitors

In this study, firstly, the pharmacophore model was established based on LAR inhibitors. ZINC database and drug-like database were screened by Hypo-1-LAR model, and the embryonic compound ZINC71414996 was obtained. Based on this compound, we designed 9 compounds. Secondly, the synthetic route of the compound was determined by consulting Reaxys and Scifinder databases, and 9 compounds (1a-1i) were synthesized by nucleophilic substitution, Suzuki reaction and so on. Meanwhile, their structures were confirmed by ¹H NMR and ¹³C NMR. Thirdly, the Enzymatic assays was carried out, the biological evaluation of compounds 1a-1i led to the identification of a novel PTP-LAR inhibitor 1c, which displayed an IC₅₀ value of 4.8 μM. At last, molecular dynamics simulation showed that compounds could interact strongly with the key amino acids LYS1350, LYS1352, ARG1354, TYR1355, LYS1433, ASP1435, TRP1488, ASP1490, VAL1493, SER1523, ARG1528, ARG1561, GLN1570, LYS1681, thereby inhibiting the protein activity. This study constructed the pharmacophore model of LAR protein, designed small-molecule inhibitors, conducted compound synthesis and enzyme activity screening, so as to provide a basis for searching for drug-capable lead compounds.     

Discovery of novel 1,2,4-triazine-chalcone hybrids as anti-gastric cancer agents via an axis of ROS-ERK-DR5 in vitro and in vivo

In this work, a series of novel 1,2,4-triazine-chalcone hybrids were designed through the molecular hybridization strategy, synthesized by two step chlorinations and further aldol condensation and evaluated their antiproliferative activity against MGC-803, HCT-116, PC-3, EC-109 and A549 cells. Compound 9l displayed significant antiproliferative activity against MGC-803, HCT-116, PC-3, EC-109 and A549 cell lines with IC₅₀ values of 0.41, 0.43, 0.61, 0.78 and 0.52 μM, respectively. Subsequent mechanistic investigations suggested that compound 9l induced the generation of ROS and inhibited the activation of the ERK pathway. Compound 9l induced extrinsic cell apoptosis by up-regulating DR5 dependent on the generation of ROS, while up-regulation of DR5 caused by compound 9l relied on the inhibition of ERK. Thus, compound 9l inhibited the gastric cancer cells via an axis of ROS-ERK-DR5 in vitro. Compound 9l also showed potent activity on cell proliferation inhibition, and was effective in suppressing the growth of MGC-803 xenograft tumor in nude mice without obvious toxicity. Therefore, compound 9l is to be reported as anti-gastric cancer agent in vitro and in vivo via an axis of ROS-ERK-DR5.     

Structure elucidation,DNA binding,DFT, molecular docking and cytotoxic activity studies on novel single crystal (E)-1-(2-fluorobenzylidene)thiosemicarbazide

Compound 3 {(E)-1-(2-fluorobenzylidene)thiosemicarbazide} – a new Schiff base of thiosemicarbazide has been synthesized, characterized and reported for crystal structure. Planer side chain in the crystal structure was observed co-planer with aromatic ring plane and molecules were connected into centrosymmetric dimmers via intermolecular hydrogen bonding. DFT geometry optimization and the relevant quantum parameters indicated unstable and reactive nature of compound 3. Experimental and theoretical findings for DNA binding by UV–visible, cyclic voltammetry and molecular docking studies showed consistency in kinetic (K_b) and thermodynamic (ΔG) parameters and that compound 3 significantly interacted with DNA via intercalation. Viscometric analysis further comprehended intercalation as possible binding mode of the compound with DNA and non-denaturing of DNA in the presence of 10% aqueous DMSO. Docked parameters further assured the drug like characteristics of the investigated compound as fit in Lipinski’s criteria. Dose dependant cytotoxic activity of compound 3 against human Huh-7 cell line indicated its anti-cancer potential at 100?µg/ml concentration.     

Het(aryl)isatin to het(aryl)aminoindoline scaffold hopping: A route to selective inhibitors of matrix metalloproteinases

Matrix metalloproteinases (MMPs) are a large family of zinc-dependent endoproteases known to exert multiple regulatory roles in tumor progression. A variety of chemical classes have been explored for targeting individual MMP isoforms. In the present study, we further developed our isatin based scaffold BB0223107 capable of binding to and inactivating MMP-2 in a zinc-independent manner (Agamennone et al., 2016). Forty four new compounds were synthesized based on the modified BB0223107. All compounds were tested in enzyme inhibition assays against MMP-2, ?8 and ?13. SAR studies demonstrated that 5-het(aryl)-3-aminoindolin-2-ones (37–39) were active toward MMP-2 and MMP-13. The most potent compounds 33 and 37 displayed an IC₅₀ of 3 µM against MMP-13 and showed a negligible activity toward MMP-8; almost all new compounds were inactive toward MMP-8. Replacement of the isatin ring with a biaryl system (compound 33) did not decrease the potency against MMP-13 but reduced the selectivity. Structure-based computational studies were carried out to rationalize the inhibitory activity data. The analysis of binding geometries confirmed that all fragments occupied the S1′ site in the three enzymes while no ligand was able to bind the catalytic zinc ion. To the best of our knowledge, this is the first example of 3-aminoindolin-2-one-based MMP inhibitors that, based on the computer modeling study, do not coordinate the zinc ion. Thus, the het(aryl)-3-aminoindolin-2-one derivatives emerge as a drug-like and promising chemotype that, along with the hetaryl variations, represents an alternative and thrifty tool for chemical space exploration aimed at MMP inhibitor design.     

Structure-based virtual screening of influenza virus RNA polymerase inhibitors from natural compounds: Molecular dynamics simulation and MM-GBSA calculation

The resistances of matrix protein 2 (M2) protein inhibitors and neuraminidase inhibitors for influenza virus have attracted much attention and there is an urgent need for new drug. The antiviral drugs that selectively act on RNA polymerase are less prone to resistance and possess fewer side effects on the patient. Therefore, there is increased interest in screening compounds that can inhibit influenza virus RNA polymerase. Three natural compounds were found by using molecular docking-based virtual screening, which could bind tightly within the polymerase acidic protein-polymerase basic protein 1 (PA-PB1) subunit of influenza virus polymerase. Firstly, their drug likeness properties were evaluated, which showed that the hepatotoxicity values of all the three compounds indicating they had less or no hepatotoxicity, and did not have the plasma protein biding (PPB) ability, the three compounds needed to be modified in some aspects, like bulky molecular size. The stability of the complexes of PA-hits was validated through molecular dynamics (MD) simulation, revealing compound 2 could form more stable complex with PA subunit. The torsional conformations of each rotatable bond of the ligands in PA subunit were also monitored, to investigate variation in the ligand properties during the simulation, compound 3 had fewer rotatable bonds, indicating that the molecule had stronger rigidity. The bar charts of protein–ligand contacts and contacts over the course of trajectory showed that four key residues (Glu623, Lys643, Asn703 and Trp706) of PA subunit that participated in hydrogen-bond, water bridge and hydrophobic interactions with the hit compounds. Finally, the binding free energy and contributed energies were calculated by using MM-GBSA method. Out of the three compounds, compound 1 showed the lowest total binding free energy. Among all the interactions, the contribution of the covalent binding and the van der Waals energy were more than other items, compound 1 formed more stable hydrogen bonds with the residues of PA subunit binding pocket. This study smoothed the path for the development of novel lead compounds with improved binding properties, high drug likeness, and low toxicity to humans for the treatment of influenza, which provided a good basis for further research on novel and effective influenza virus PA-PB1 interaction inhibitors.     

Evaluation of quercetin as a potential β-lactamase CTX-M-15 inhibitor via the molecular docking,dynamics simulations,and MMGBSA

Antimicrobial resistance (AMR) threatens millions of people around the world and has been declared a global risk by the World Economic Forum. One of the important AMR mechanisms in Enterobacteriaceae is the production of extended-spectrum β-lactamases. The most common ESBL, CTX-M β-lactamases, is spread to the world by CTX-M-15 and CTX-M-14. Sulbactam, clavulanic acid, and tazobactam are first-generation β-lactamase inhibitors and avibactam is a new non-β-lactam β-lactamase inhibitor. We studied that avibactam, sulbactam, clavulanic acid, tazobactam, and quercetin natural flavonoids were docked to target protein CTXM-15. Subsequently, the complexes were simulated using the molecular dynamics simulations method during 100 ns for determining the final binding positions of ligands. Clavulanic acid left CTX-M-15 and other ligands remained in the binding site after the simulation. The estimated binding energies were calculated during 100 ns simulation by the MMGBSA-MMPBSA method. The estimated free binding energies of avibactam, sulbactam, quercetin, tazobactam, and clavulanic acid were sorted as –33.61 kcal/mol, –16.04 kcal/mol, –14 kcal/mol, –12.68 kcal/mol, and –2.95 kcal/mol. As a result of both final binding positions and free binding energy calculations, Quercetin may be evaluated an alternative candidate and a more potent β-lactamases inhibitor for new antimicrobial combinations to CTX-M-15. The results obtained in silico studies are predicted to be a preliminary study for in vitro studies for quercetin and similar bioactive natural compounds. These studies are notable for the discovery of natural compounds that can be used in the treatment of infections caused by β-lactamase-producing pathogens.     

Extra precision docking,free energy calculation and molecular dynamics studies on glutamic acid derivatives as MurD inhibitors

The binding modes of well known MurD inhibitors have been studied using molecular docking and molecular dynamics (MD) simulations. The docking results of inhibitors 1-30 revealed similar mode of interaction with Escherichia coli-MurD. Further, residues Thr36, Arg37, His183, Lys319, Lys348, Thr321, Ser415 and Phe422 are found to be important for inhibitors and E. coli-MurD interactions. Our docking procedure precisely predicted crystallographic bound inhibitor 7 as evident from root mean square deviation (0.96 Å). In addition inhibitors 2 and 3 have been successfully cross-docked within the MurD active site, which was pre-organized for the inhibitor 7. Induced fit best docked poses of 2, 3, 7 and 15/2Y1O complexes were subjected to 10 ns MD simulations to determine the stability of the predicted binding conformations. Induce fit derived docked complexes were found to be in a state of near equilibrium as evident by the low root mean square deviations between the starting complex structure and the energy minimized final average MD complex structures. The results of molecular docking and MD simulations described in this study will be useful for the development of new MurD inhibitors with high potency.     

Synthesis,anticancer activity,toxicity evaluation and molecular docking studies of novel phenylaminopyrimidine—(thio)urea hybrids as potential kinase inhibitors

Thirty-two novel urea/thiourea compounds as potential kinase inhibitor were designed, synthesized and evaluated for their cytotoxic activity on breast (MCF7), colon (HCT116) and liver (Huh7) cancer cell lines. Compounds 10, 19 and 30 possessing anticancer activity with IC₅₀ values of 0.9, 0.8 and 1.6 μM respectively on Huh7 cells were selected for further studies. These hit compounds were tested against liver carcinoma panel. Real time cell electronic sensing assay was used to evaluate the effects of the compounds 10, 19 and 30 on the growth pattern of liver cancer cells. Apoptotic cell death and cell cycle analysis upon treatment of liver carcinoma cells with hit compounds were determined. A significant apoptotic cell death was detected upon treatment of Huh7 and Mahlavu cells with compound 30 after 48 h of treatment. Additionally, compound 10 caused cell cycle arrest at G0/G1 phase. Mutagenicity of hit compounds was evaluated. Assertively, these compounds were not found to be mutagenic on Salmonella typhimurium strains TA98 and TA100. To understand the binding modes of the synthesized compounds, molecular docking studies were performed using the crystal data of VEGFR and Src-kinase enzymes in correlation with anticancer activities.     

Multicomponent synthesis,cytotoxicity, and computational studies of novel imidazopyridazine-based N-phenylbenzamides

A one-pot multicomponent synthesis and application of new imidazopyridazine based N-phenylbenzamides is described. An atom-economical method involving dimethyl phthalate, substituted anilines, and pyridazine-4,5-diamine provided the desired compounds in 120–150 min with 80–85% yield. The reaction was catalyzed with phosphoric acid, and glycerol was used as a safer, greener solvent. Anticancer evaluation against selected cancer cell lines revealed that compound 4e was the most active from the series and exhibited IC₅₀ values below 9.1 µM. Compounds 4h and 4d also displayed good and comparable IC₅₀ values (10.2–12.1 µM). Molecular docking and molecular dynamic studies showed that compound 4e exhibit good binding affinity and stable complex formation with ABL1-kinase protein, respectively. Additional computational predictions such as ADME and drug-likeness demonstrated the potential of the new benzamides as leads for further development.     

Enantioselective synthesis of new oxazolidinylthiazolidines as enzyme inhibitors

The synthesis of new oxazolidinylthiazolidines bicycles, oxygen analogues of bisthiazolidines, also known as metallo-β-lactamase inhibitors is described. The reaction of β-aminoalcohols and 2,5-dihydroxy-1,4-dithiane led to oxazolidinylthiazolidines and/or dithioazabicycles as the main products. The distribution pattern depends mainly on the aminoalcohol substituents. In a one-pot reaction, four new bonds are formed in good yields and with high atom efficiency. When the oxazolidinylthiazolidines are formed, two stereogenic centres are generated with high enantiospecificity. The reaction mechanism is discussed based on crystallographic data and interconversion studies. Two oxazolidinylthiazolidines were evaluated as inhibitors of the potent lactamase NDM-1 and compound 4f displayed competitive inhibition with K_i = 1.6 ± 0.6 μM.