Amit K Keshari1, Ashok K Singh1, Vinit Raj1, Amit Rai1, Prakruti Trivedi2, Balaram Ghosh2, Umesh Kumar3, Atul Rawat3, Dinesh Kumar3, Sudipta Saha1. 1. Department of Pharmaceutical Sciences, Babasaheb Bhimrao Ambedkar University, Vidya Vihar, Raebareli Road, Lucknow, Uttar Pradesh. 2. Department of Pharmacy, Birla Institute of Technology & Science Pilani, Hyderabad Campus, Hyderabad, Telangana State. 3. Centre of Biomedical Research (CBMR), Sanjay Gandhi Post-Graduate Institute of Medical Sciences Campus, Raebareli Road, Lucknow, Uttar Pradesh, India.
Abstract
In our efforts to address the rising incidence of hepatocellular carcinoma (HCC), we have made a commitment to the synthesis of novel molecules to combat Hep-G2 cells. A facile and highly efficient one-pot, multicomponent reaction has been successfully devised utilizing a p-toluenesulfonic acid (p-TSA)-catalyzed domino Knoevenagel/Michael/intramolecular cyclization approach for the synthesis of novel 5H-benzo[h]thiazolo[2,3-b]quinazoline and indeno[1,2-d] thiazolo[3,2-a]pyrimidine analogs bearing a bridgehead nitrogen atom. This domino protocol constructed one new ring by the concomitant formation of multiple bonds (C-C, C-N, and C=N) involving multiple steps without the use of any metal catalysts in one-pot, with all reactants effi-ciently exploited. All the newly synthesized compounds were authenticated by means of Fourier transform infrared spectroscopy, liquid chromatography-mass spectrometry, proton nuclear magnetic resonance spectroscopy, and carbon-13 nuclear magnetic resonance spectroscopy, together with elemental analysis, and their antitumor activity was evaluated in vitro on a Hep-G2 human cancer cell line by sulforhodamine B assay. Computational molecular modeling studies were carried out on cancer-related targets, including interleukin-2, interleukin-6, Caspase-3, and Caspase-8. Two compounds (4A and 6A) showed growth inhibitory activity comparable to the positive control Adriamycin, with growth inhibition of 50% <10 μg/mL. The results of the comprehensive structure-activity relationship study confirmed the assumption that two or more electronegative groups on the phenyl ring attached to the thiazolo[2,3-b]quinazoline system showed the optimum effect. The in silico simulations suggested crucial hydrogen bond and π-π stacking interactions, with a good ADMET (absorption, distribution, metabolism, excretion, and toxicity) profile and molecular dynamics, in order to explore the molecular targets of HCC which were in complete agreement with the in vitro findings. Considering their significant anticancer activity, 4A and 6A are potential drug candidates for the management of HCC.
In our efforts to address the rising incidence of hepatocellular carcinoma (HCC), we have made a commitment to the synthesis of novel molecules to combat Hep-G2 cells. A facile and highly efficient one-pot, multicomponent reaction has been successfully devised utilizing a p-toluenesulfonic acid (p-TSA)-catalyzed domino Knoevenagel/Michael/intramolecular cyclization approach for the synthesis of novel 5H-benzo[h]thiazolo[2,3-b]quinazoline and indeno[1,2-d] thiazolo[3,2-a]pyrimidine analogs bearing a bridgehead nitrogen atom. This domino protocol constructed one new ring by the concomitant formation of multiple bonds (C-C, C-N, and C=N) involving multiple steps without the use of any metal catalysts in one-pot, with all reactants effi-ciently exploited. All the newly synthesized compounds were authenticated by means of Fourier transform infrared spectroscopy, liquid chromatography-mass spectrometry, proton nuclear magnetic resonance spectroscopy, and carbon-13 nuclear magnetic resonance spectroscopy, together with elemental analysis, and their antitumor activity was evaluated in vitro on a Hep-G2humancancer cell line by sulforhodamine B assay. Computational molecular modeling studies were carried out on cancer-related targets, including interleukin-2, interleukin-6, Caspase-3, and Caspase-8. Two compounds (4A and 6A) showed growth inhibitory activity comparable to the positive control Adriamycin, with growth inhibition of 50% <10 μg/mL. The results of the comprehensive structure-activity relationship study confirmed the assumption that two or more electronegative groups on the phenyl ring attached to the thiazolo[2,3-b]quinazoline system showed the optimum effect. The in silico simulations suggested crucial hydrogen bond and π-π stacking interactions, with a good ADMET (absorption, distribution, metabolism, excretion, and toxicity) profile and molecular dynamics, in order to explore the molecular targets of HCC which were in complete agreement with the in vitro findings. Considering their significant anticancer activity, 4A and 6A are potential drug candidates for the management of HCC.
One in six human deaths globally is due to multifaceted disease, such as cancer. Liver cancer accounted for 788,000 deaths in 2015 and is the second most common cause of cancer-related deaths, worldwide. The most prevalent primary liver cancer with poor prognosis is hepatocellular carcinoma (HCC), which has potentially lethal human malignacy worldwide, particularly in Asia and Africa.1,2 In men, HCC is the fifth most commonly diagnosed cancer and the second leading cause of cancer deaths; in women, it is the seventh most frequently diagnosed cancer, and worldwide, overall, it is the sixth commonest cause of cancer-related deaths.3 In spite of continuous efforts to develop novel therapeutic strategies to treat HCC, it remains a challenge. As a consequence, therefore, increased attention has been paid, especially in the field of liver cancer therapy, to the discovery and development of safe and novel anticancer agents, together with improved cytotoxicity toward cancerous cells.4Domino reactions involving simultaneous formation of C–C and C–heteroatom multiple bonds in a single flask are a means to achieving economical methods for the manufacture of especially unusual, fused heterocyclic, medicinally privileged scaffolds without any separation of intermediates throughout the process. Moreover, domino reactions, being attractively appealing, are frequently connected with savings in terms of energy and reaction periods, are highly convergent, have fewer environmental impacts, and show atom economy in a single-step process starting from multiple reactants.5,6Conventional multistep chemical reactions suffer from various limitations compared with the one-pot domino technique, as they require a large number of synthetic operations, including isolation and purification of the products of each individual step. The multistep synthetic approach, therefore, has led to synthetic inefficiency, is time-consuming, or contains side reactions generating large amounts of waste.7A key challenge in recent drug development processes is to design a rapid, versatile, and efficient synthesis that provides target molecules containing structural complexity and diversity with a choice of fascinating biologic activities.8Thiazole-based heterocycles are central to modern chemical synthesis due to their synthetic feasibility and their incorporation into various types of therapeutically useful agents and are of paramount interest in the development of important pharmacophores in the drug discovery endeavor.9,10 Functionalized quinazoline and fused quinazolines have long been of increasing interest in the field of synthetic organic and medicinal chemistry on account of the wide range of biologic activities,11–14 including anticancer,15 anti-inflammatory,16 antituberculosis,17 anticonvulsant,18 antimalarial,19 antihypertensive,20 antidiabetic,21 and so on. Furthermore, suitably substituted thiazolo[2,3-b]quinazolines also serve as versatile building blocks in synthetic chemistry.Nowadays, bridgehead (ring junction) nitrogen atom-containing thiazolo[3,2-a]pyrimidines command much attention as privileged scaffolds comprising a vital class of heterocyclic structures possessing exciting and varied phar-macologic activities,22–24 such as: being potential inhibitors of cyclin-dependent kinase,25 CDC25B phosphatase enzymes,26 IspF proteins,27 and YycG histidine kinase,28 antibacterial,29 analgesic, and anti-inflammatory activities,30 acting as a corticotropin-releasing factor receptor antagonist,31 hypolipidemic activity,32 antiviral including anti-HIV effects, anticancer,33,34 cardiotonic, and inotropic effects,35 acting as potent A2A adenosine receptor inverse agonists with antinociceptive activity,36 and so on.To the best of our knowledge, among the earlier reported methods devoted to the synthesis of especially 5H-benzo[h] thiazolo[2,3-b]quinazolines,37,38 only a few are described, and metal-free multicomponent domino reactions (MDRs)39,40 are still rare. In this regard, most of the reported methods have suffered from various drawbacks, for example, harsh reaction conditions, multistep synthetic routes, costly reagents/catalysts, prolonged reaction periods, tedious workups after each step, and poor availability of starting materials. There is no route for the direct construction of the indeno[1,2-d] thiazolo[3,2-a]pyrimidine moiety through domino reactions, and a literature survey revealed only one report to be available, which includes the synthesis of 5-phenylindeno[1,2-d] thiazolo[3,2-a]pyrimidine-3,6(2H,5H)-dione and 5-phenyl-2,3-dihydro indeno[1,2-d]thiazolo[3,2-a]pyrimidin-6(5H)-one in two steps.41In order to arrive at a highly proficient and convergent synthetic strategy for the construction of these two vital structural elements, while at the same time avoiding several limitations of earlier reports, and in our venture toward the development of a modern synthetic approach, we first report herein an operationally simple and straightforward metal-free, one-pot MDR for obtaining a panel of novel 5H-benzo[h]thiazolo[2,3-b]quinazoline (1A–15A) and indeno[1,2-d]thiazolo[3,2-a]pyrimidine (1B–15B) analogs using the reactions of highly substituted α-tetralone or α-indanone with some aromatic aldehydes and distinctive 2-aminothiazoles in ethanol (EtOH) in the presence of catalytic amounts of p-toluenesulfonic acid (p-TSA; 20 mol%).7,42,43 Readily available and cheap starting materials, together with an environmentally benign and mild acidic catalyst, were employed to achieve these diversely decorated skeletons in impressive yields. They have also been synthesized by conventional two-step reactions and compared with the preferred MDRs described above.23,44–46 Additionally, all of the newly synthesized compounds were evaluated for their antitumor activity in vitro on a Hep-G2humancancer cell line by sulforhodamine B (SRB) assay and structure–activity relationship prediction on the basis of in vitro findings. Furthermore, computational molecular modeling studies with ADMET (absorption, distribution, metabolism, excretion, and toxicity) profiling and molecular dynamic (MD) simulation were carried out on cancer-related targets, including interleukin (IL)-2, IL-6, Caspase-3, and Caspase-8 receptor sites, to identify the potential modes of action of the named compounds.
Materials and methods
General information
All of the reactions described below were carried out with commercially available chemicals of reagent grade that were used as received without further purification unless otherwise noted. Reagents were purchased from Merck and Sigma-Aldrich chemical companies. Melting points (m.p.) were determined on an m.p. apparatus and are uncorrected. Compounds were named following the International Union of Pure and Applied Chemistry rules as applied by Advanced Chemistry Development/ChemSketch. Elemental analyses were performed with a Euro EA 3000 Elemental Analyzer for C, H, N, and the results were within ±0.4% of the theoretical values. Infrared (IR) spectra were recorded on a Thermo Scientific Nicolet 6700 Fourier transform IR spectroscopy spectrometer. Nuclear magnetic resonance (NMR) spectra were obtained on a Bruker 800 MHz NMR spectrometer (1H 800 MHz, 13C 200 MHz) (Bruker, Rheinstetten, Germany) and processed in TopSpin 2.1. Chemical shifts are expressed in parts per million (ppm) downfield from tetramethylsilane as the internal standard. Coupling constants are expressed in Hz. All of the compounds were analyzed for mass data using a liquid chromatography–mass spectrometry (LCMS)-2020 mass spectrometer (Schimadzu, Tokyo, Japan). All the compounds were dissolved in 1:1 (v/v) mixtures of acetonitrile:methanol, and 10 μL of the resulting solution was injected to acquire the data. The analysis was performed in electrospray ionization mode using a capillary column at a flow rate of 0.2 mL/min with a 50% water/methanol (1:1) mixture for 120 s. Data analysis was performed using lab solution LCMS data processing software. Reaction progress as well as the purity of the compounds were evaluated with thin layer chromatography (TLC) plates by using ethylacetate:n-hexane (3:7) as the eluent, and the developed chromatogram was visualized under ultraviolet light and iodine vapors.
General experimental procedures for characterization of the synthesized compounds
One-pot efficient synthesis of substituted 5H-benzo[h]thiazolo[2,3-b]quinazolines (1A–15A)
A mixture of substituted tetralone (1 mmol), appropriate aromatic aldehydes (1 mmol), and distinctive 2-aminothiazoles (1 mmol) in EtOH (5.0 mL) in the presence of 20 mol% p-TSA was heated under reflux for 3–4 h. The reaction mixture was poured into ice-cold water. A solid product was obtained, which was filtered, washed thoroughly with distilled water, and recrystallized from EtOH. Pure crystals were obtained (Scheme 1). The progress of the reaction was monitored by TLC on precoated silica gel-G plates using 30% ethylacetate:n-hexane as the solvent system. TLC revealed just a single spot, which proved the presence of a single product.
Scheme 1
One-pot efficient synthetic route to the titled compound (1A–15A).
Two-step synthesis of substituted 5H-benzo[h] thiazolo[2,3-b]quinazolines (1A–15A)
These compounds were synthesized as follows.
Step I
A mixture of substituted tetralones (1 mmol), appropriate aldehydes (1 mmol), piperidine (0.1 mmol), and glacial acetic acid (0.1 mmol) in EtOH (5.0 mL) was heated under reflux for 10–12 h. The reaction mixture was cooled to room temperature for a few minutes and then poured into ice-cold water. A solid product, substituted benzylidene-tetralones, was obtained, which was filtered, washed thoroughly with distilled water, and dried. The dry residue was recrystallized from EtOH. The progress of the reaction was monitored by TLC on precoated silica gel-G plates using 40% ethylacetate:n-hexane as the solvent system.
Step II
A solution of distinctive 2-aminothiazoles (1 mmol) and appropriate benzylidene-tetralones in glacial acetic acid (2.0 mL) was heated under reflux for 18–20 h. The solvent was evaporated under vacuum and then the residue obtained was dissolved in chloroform, washed with water, and the organic layer separated, dried, and evaporated. The resulting solid obtained was recrystallized from EtOH to yield the target compounds (1A–15A; Scheme 2).
Scheme 2
Two-step synthetic route to the titled compound (1A–15A).
One-pot efficient synthesis of substituted indeno[1,2-d]thiazolo[3,2-a]pyrimidines (1B–15B)
A mixture of substituted α-indanone (1 mmol), appropriate aromatic aldehydes (1 mmol), and distinctive 2-aminothiazoles (1 mmol) in EtOH (5.0 mL) in the presence of 20 mol% p-TSA was heated under reflux for 3–4 h. The reaction mixture was then poured into ice-cold water. A solid product was obtained, which was filtered, washed thoroughly with distilled water, and recrystallized from EtOH. Pure crystals were obtained (Scheme 3). The progress of the reaction was monitored by TLC on precoated silica gel-G plates using 30% ethylacetate:n-hexane as the solvent system. TLC revealed just a single spot, which proved the presence of a single product.
Scheme 3
One-pot efficient synthetic route to the titled compound (1B–15B).
Two-step synthesis of substituted indeno[1,2-d] thiazolo[3,2-a]pyrimidines (1B–15B)
These compounds were synthesized as follows.A mixture of substituted indanones (1 mmol), appropriate aldehydes (1 mmol), piperidine (0.1 mmol), and glacial acetic acid (0.1 mmol) in EtOH (5.0 mL) was heated under reflux for 10–12 h. The reaction mixture was cooled to room temperature for a few minutes and then poured into ice-cold water. A solid product, substituted benzylidene-indanones, was obtained, which was filtered, washed thoroughly with distilled water, and dried. The dry residue was recrystallized from EtOH. The progress of the reaction was monitored by TLC on precoated silica gel-G plates using 40% ethylacetate:n-hexane as the solvent system.A solution of distinctive 2-aminothiazoles (1 mmol) and appropriate benzylidene-indanones in glacial acetic acid (2.0 mL) was heated under reflux for 18–20 h. The solvent was evaporated under vacuum and then the residue obtained was dissolved in chloroform, washed with water, and the organic layer separated, dried, and evaporated. The resulting solid obtained was recrystallized from EtOH to yield the target compounds (1B–15B; Scheme 4).
Scheme 4
Two-step synthetic route to the titled compound (1B–15B).
The Hep-G2 cells derived from humanliver carcinoma were purchased from the National Cell Repository NCCS, Pune, India and later cultured in the lab of ACTREC, Tata Memorial Centre, Navi Mumbai, India for in vitro testing purposes. As per the instructions, the cells were grown in Roswell Park Memorial Institute medium 1640 containing 10% fetal bovine serum. Then, 100 μL cells containing media were seeded into 96-well plates at a density of 5×103 cells/well and incubated at 37°C in a humidified atmosphere containing 5% CO2 for 24 h prior to addition of the experimental drugs. Initially, all the experimental drugs were solubilized in DMSO at 100 mg/mL and further diluted to 1.0 mg/mL using water and stored frozen prior to use. Aliquots of 100 μL (containing 90 μL of medium) of various dilutions were added to appropriate wells, resulting in the required drug concentrations of 10, 20, 40, and 80 μg/mL, and were maintained at standard conditions for 48 h. For each of the experiments, a well-known anticancer drug Adriamycin (ADR) was used as a positive control at concentrations of 10, 20, 40, and 80 μg/mL. The assay was terminated by the gentle addition of 50 μL of cold 30% (w/v) trichloroacetic acid and the mixture was incubated for 60 min at 4°C. The supernatant was discarded and the plates were washed five times using tap water and then air dried. Next, SRB solution (50 μL) at 0.4% (w/v) in 1% acetic acid was added to each well and incubated for 20 min at room temperature. After staining, the unbound dye was recovered and the residual dye was removed by washing five times with 1% acetic acid and then air dried. Subsequently, the bound stain was eluted with 10 mM trizma base and the absorbance was examined at a wavelength of 540 nm on an enzyme-linked immunosorbent assay plate reader.The growth inhibition of 50% (GI50) was calculated using the formula [(Ti-Tz)/(C-Tz)] ×100%, where Tz represents time zero growth, C is the control growth, and Ti is the test growth at four concentration levels in the presence of the drug.47,48
Molecular docking
The structures of the ligands were designed by Chem Draw Ultra 12.0, and the geometry was optimized six times through GaussView 5.0 (trial version). In addition, the National Centre for Biotechnology Information49 and the Protein Data Bank50 were utilized as chemical sources to obtain the established four homological cancer protein targets, that is, IL-2 (1Z92), IL-6 (1IL6), Caspase-3 (1QX3), and Caspase-8 (1IBC). DS Visualizer software51 was used to remove the co-crystal with the assigned target. Next, the active site was identified using the CASTp database.52 For a particular system, validation of the docking protocol was performed using re-docking studies, where we found excellent agreement between the localization of ligands upon docking, and from the crystal structure of the assigned targets. We demonstrated the reliability and quality of the docking method in reproducing experimentally observed binding-mode-assigned targets. In addition, the docking analyses of the test set were carried out using Autodock4.1,53 together with the Lamarkian genetic algorithm for automated flexible ligand docking, and the binding energy was estimated as negative kcal/mol. In addition, probable H-bonds and π-bonds were also assessed.54−56
Prediction of ADME properties
The ADME and drug-like properties of selected ligands were predicted by employing MedChem Designer and QikProp tools. Chemical structures were optimized with LigPrep. Additionally, the percentage absorbance (% ABS) and Lipinski’s violation were evaluated for this study.57
MD simulation
The nature of the inhibitor used on the active site domain of IL-6 was investigated by means of MD simulation studies. The energies of dock configurations were minimized to eliminate the unfavorable atomic contacts as starting conformations for dynamic simulation in Elmar Krieger MD simulation tools (trial version).58 An AMBER03 force field was assigned to execute a real-time MD simulation.59 The complex was solvated through an HOH model at density =0.997 g/L inside the 10Å simulation cell boundary and then adjusted to the physiologic pH at 7.4. Moreover, physiologic NaCl solution with 0.9% mass fraction Na+ and Cl− ion concentration was used to maintain and neutralize the simulation cell boundary. Then, the MD simulation was run for 3,000 ps at a temperature of 298 K and 1 bar pressure to obtain snapshot (sim) trajectories. Finally, the sim trajectories were analyzed and the resulting data were plotted using Sigma Plot 11.0 tools.
Results and discussion
Design
These two core structural motifs are formed by the fusion of three biodynamic privileged heterosystems in such a way that one nitrogen atom occupies a bridge head position, therefore being common to both the heterocyclic rings, that is, the thiazole and the pyrimidine rings, and possessing unique structural diversity. Thiazolo[2,3-b]quinazoline and thiazolo[3,2-a]pyrimidine, as heterocyclic skeletons, exhibit a diverse range of pharmacologic activities. In view of the great medicinal importance associated with these remarkable structures, we are interested in studying MDRs with the aim of evolving an efficient and highly convergent synthetic procedure for the library of densely functionalized 5H-benzo[h]thiazolo[2,3-b]quinazoline (1A−15A) and indeno[1,2-d]thiazolo[3,2-a]pyrimidine (1B−15B) analogs in a single operation.We initiated our study of the proposed method with 10 mol% p-TSA as the acidic catalyst (Scheme 5) by investigating the conditions for the reaction comprising 5-methoxy-2,3-dihydro-1H-inden-1-one (1 mmol), 4-bromobenzaldehyde (1 mmol), 4-methylthiazol-2-amine (1 mmol), and EtOH (5.0 mL) at 60°C to yield the corresponding product (4B, 51% yield) as a model in 4.5 h (entry 1). With this finding in hand, when the acidic catalyst was changed from 10% to 20% and the temperature from 60°C to 80°C, the product yield increased to 69% after 3 h reaction time (entry 4). We examined the reaction in different organic solvents, including EtOH, methanol (MeOH), ethylacetate (EtOAc), acetonitrile (CH3CN), and toluene (PhCH3) under varied temperature conditions (Table 1). The proposed reaction is optimized for the amount of acidic catalyst required, as well as for temperature, to improve the yields of the desired molecules (entries 1–8), and the optimum amount was found to be 20 mol% p-TSA in EtOH at 80°C. The best results were obtained for 20 mol% p-TSA in 3 h at 80°C (entry 4, 69% yield).
In order to explore the scope of these conditions, we embarked on the synthesis of a panel of novel 5H-benzo[h] thiazolo[2,3-b]quinazoline (1A–15A) and indeno[1,2-d] thiazolo[3,2-a]pyrimidine (1B–15B) analogs bearing a bridgehead nitrogen atom from a range of substrates by the advancement of facile and efficient MDRs of highly substituted α-tetralone (a) or α-indanone (a′) with some aromatic aldehydes (b) and distinctive 2-aminothiazoles (c) in EtOH in the presence of 20 mol% p-TSA. The α-tetralone, 6-methoxy-1-tetralone, 7-methoxy-1-tetralone or α-indanone, 5-methoxy-1-indanone, and the substituted benzaldehydes with 4-bromo, trimethoxy, 3-hydroxy-4-methoxy (vanillin), 4-chloro were reacted with distinctive 2-aminothiazole, 4-methylthiazol-2-amine, 5-methylthiazol-2-amine to afford a library of proposed target molecules, that is, (1A–15A) and (1B–15B), respectively.The one-pot, three-component reaction went smoothly in EtOH in the presence of p-TSA and gave the targeted compounds (1A–15A) and (1B–15B) in impressive yields up to 86%. The various substitutions of synthesized derivatives are given in Table 2A and B. All the newly synthesized compounds were authenticated on the basis of Fourier transform infrared spectroscopy, LCMS, and 1H and 13C NMR spectral and elemental analyses.
Table 2
Various substitutions and in vitro cytotoxicity data of synthesized derivatives against human hepatoma (Hep-G2) cell lines: A (1A–15A) and B (1B–15B)
Comp code
X1
X2
R1
R2
R3
Y1
Y2
*GI50 (µg/mL)
LC50 (µg/mL)
TGI (µg/mL)
A
1A
–H
–H
–H
–Br
–H
–H
–H
>80
NE
NE
2A
–OCH3
–H
–H
–Br
–H
–H
–CH3
NE
NE
NE
3A
–H
–OCH3
–OCH3
–OCH3
–OCH3
–H
–H
NE
NE
NE
4A
–H
–OCH3
–OCH3
–OH
–H
–H
–H
<10
>80
>80
5A
–OCH3
–H
–H
–Br
–H
–H
–H
>80
>80
>80
6A
–OCH3
–H
–OCH3
–OCH3
–OCH3
–H
–H
<10
NE
NE
7A
–H
–H
–H
–OH
–H
–H
–H
76.1
NE
NE
8A
–H
–OCH3
–H
–Cl
–H
–CH
–H
>80
NE
NE
9A
–OCH3
–H
–H
–Cl
–H
–CH
–H
>80
NE
NE
10A
–H
–H
–H
–Cl
–H
–CH
–H
>80
NE
NE
11A
–H
–OCH3
–H
–Cl
–H
–H
–H
>80
NE
NE
12A
–H
–H
–H
–Cl
–H
–H
–H
>80
NE
NE
13A
–H
–OCH3
–H
–Cl
–H
–H
–CH
>80
NE
NE
14A
–OCH3
–H
–H
–Cl
–H
–H
–CH
>80
NE
NE
15A
–H
–H
–H
–Cl
–H
–H
–CH
.80
NE
NE
Comp code
X
R1
R1
R2
R3
Y1
Y2
*GI50 (µg/mL)
LC50 (µg/mL)
TGI (µg/mL)
B
1B
–H
–H
–H
–Br
–H
–H
–H
>80
NE
NE
2B
–H
–Cl
–H
–H
–H
–H
–H
>80
NE
NE
3B
–H
–H
–H
–Br
–H
–H
–CH3
>80
NE
NE
4B
–OCH3
–H
–H
–Br
–H
–CH3
–H
>80
NE
NE
5B
–H
–H
–H
–Br
–H
–CH3
–H
>80
NE
NE
6B
–OCH3
–H
–H
–Br
–H
–H
–H
>80
NE
NE
7B
–OCH
–H
–H
–Br
–H
–H
–CH3
>80
NE
NE
8B
–H
–H
–OCH3
–OCH3
–OCH3
–H
–H
>80
NE
NE
9B
–OCH3
–H
–OCH3
–OCH3
–OCH3
–H
–H
>80
NE
NE
10B
–H
–H
–H
–Cl
–H
–CH3
–H
>80
NE
NE
11B
–OCH3
–H
–H
–Cl
–H
–CH3
–H
>80
NE
NE
12B
–H
–H
–H
–Cl
–H
–H
–H
>80
NE
NE
13B
–OCH3
–H
–H
–Cl
–H
–H
–H
>80
NE
NE
14B
–H
–H
–H
–Cl
–H
–H
–CH3
>80
NE
NE
15B
–OCH3
–H
–H
–Cl
–H
–H
–CH3
>80
NE
NE
ADR
<10
NE
<10
Notes:
GI50 value of ≤10 μg/mL (or 1 μM) is considered to demonstrate activity in case of pure compounds (synthetic compound). Bold values indicate the active compounds.
Abbreviations: ADR, Adriamycin, positive control compound; GI50, concentration of drug causing 50% inhibition of cell growth; LC50, concentration of drug causing 50% cell kill; NE, “not effective” even at the concentration >80 μg/mL; TGI, concentration of drug causing total inhibition of cell growth.
Desired compounds have also been synthesized by a multistep synthesis involving the reaction of distinctive 2-aminothiazoles with the appropriate benzylidene-tetralones/indanones achieved by the reaction of highly substituted α-tetralone or α-indanone with different aromatic aldehydes, piperidine and glacial acetic acid in EtOH. This reaction proceeds in two steps with a comparatively longer reaction time (18–20 h) and giving moderate yields (40%–50%).
Plausible mechanism
These reactions presumably proceed through a Knoevenagel condensation between substituted α-tetralone or α-indanone and some appropriate aromatic aldehydes in the first step to construct α,β-unsaturated ketones, respectively, which undergo a Michael-type addition approach with the nucleophilic endocyclic nitrogen of the distinctive 2-aminothiazole under the maintained reaction conditions. Then, successive intramolecular cyclization occurred with the loss of a water molecule to give 5H-benzo[h] thiazolo[2,3-b]quinazolines (1A−15A) and indeno[1,2-d] thiazolo[3,2-a]pyrimidines (1B−15B). In this setting, the domino approach and the reaction sequence of Knoevenagel condensation/Michael-type addition/intramolecular cyclization were done in a single step in a one-pot procedure in EtOH. The plausible mechanism for the reaction is delineated in Scheme 6.
Scheme 6
Plausible scenario to account for the formation of (1A–15A) and same for the preparation of (1B–15B).
A library of 30 5H-benzo[h]thiazolo[2,3-b]quinazoline and indeno[1,2-d]thiazolo[3,2-a]pyrimidine derivatives were screened against Hep-G2 cells. The inhibition activities (GI50) of the synthetic compounds 1A–15A, 1B–15B, and ADR on Hep-G2 cells are summarized in Table 2A and B. The effects of treatment with the most active members (4A and 6A) of this study and with ADR on Hep-G2 are shown by the plotted growth curve (Figure 1) and the microscopic images (Figure 2).
Figure 1
Growth Curve of (1A–15A and 1B–15B): Human Hepatoma Cell Line (Hep-G2).
Abbreviation: ADR, Adriamycin.
Figure 2
Microscopic pictures showing the effect of treatments with the active compounds (A) 4A and (B) 6A and the reference drug (C) ADR on Hep-G2 human liver cancer cell line.
Note: Scale bar 50 μm.
Abbreviation: ADR, Adriamycin.
Close examination of the activity data (Table 2A and B) revealed that the presence of –OH and –OCH3 groups in the phenyl ring, preferably at the R1, R2, and R3 positions, in conjugation with the introduction of the –OCH3 group in the tetralone ring at X1 or X2 is crucial for inhibitory activity. Substitution of 3-methoxy-4-hydroxy on the phenyl ring and 6-methoxy on the tetralone ring led to compound 4A, while substitutions of 3,4,5-trimethoxy on the phenyl ring and 7-methoxy on the tetralone ring led to compound 6A with better cytotoxicity (GI50 <10 μg/mL). In contrast, substitutions with –Br, –Cl, and –OH groups on the phenyl ring at the R2 position were ineffective against Hep-G2 cells. The methyl substitution on the thiazole ring did not contribute to producing inhibitory activity. Replacement of the tetralone ring with the indanone ring system, that is, 1B–15B, reduced the anticancer activity against Hep-G2humanliver cancer cell lines.The growth curve of the in vitro findings suggested that the percentage growth inhibition values of the potent compounds 4A and 6A was ≤50% at 10 μg/mL concentration, but they did not move toward a negative value. Therefore, it might be expected in future that both compounds could lead to the death of cancerous cells while minimizing that of normal cells.
Molecular docking studies
Computational analysis was performed using various well-known liver cancer targets, specifically IL-2, IL-6, Caspase-3, and Caspase-8, through Autodoc 4.1, together with Lamarkian genetic algorithm for automated docking of flexible ligands. Molecular docking poses of both the potent compounds 4A and 6A with the related targets IL-2, IL-6, Caspase-3, and Caspase-8 are presented in Figure 3 which indicates the interactions of amino acids with ligands, H-bonds, and π-bonds as well as the bond lengths. The docking affinities (kcal/mol), numbers of H-bonds and π-bonds, and the interactions of amino acids of the potent members with them are given in Table 3. Compounds 4A and 6A exhibited strong affinity (binding interaction energy ranged from −6.1 to −7.87 kcal/mol), while others showed good binding affinity with preferred molecular targets. Compound 4A displayed excellent affinity with IL-6 (−7.87 kcal/mol, 2H and 4π-bonds), IL-2 (−7.68 kcal/mol, 5π-bonds), Caspase-3 (−6.69 kcal/mol, 1H and 7π-bonds), and Caspase-8 (−6.48 kcal/mol, 1π-bond). A similar pattern was observed for compound 6A with IL-2 (−7.63 kcal/mol, 1H and 1π-bonds), IL-6 (−7.43 kcal/mol, 2H and 4π-bonds), and Caspase-8 (−6.76 kcal/mol, 1H and 5π-bond), but with less affinity for Caspase-3 (−6.14 kcal/mol, 1H and 6π-bonds). Moreover, both the active members of the present study have strong affinity binding to IL-2, Caspase-3, Caspase-8, and the binding energies, particularly for the IL-6 receptor site, are predominantly high (from −7.43 to −7.87 kcal/mol). Accordingly, it might be assumed that the promising cytotoxic properties of these compounds, which were indicated by the in vitro antitumor activity on Hep-G2 cells, might be better mediated through an IL-6-dependent mechanism (Table 3).
Figure 3
Docking images of active compounds (A, C, E, G) 4A and (B, D, F, H) 6A for caspase-3, caspase-8, IL-2 and IL-6.
Abbreviation: IL, interleukin.
Table 3
Docking affinity of active compounds with assigned anticancer receptors
Ligands
Receptors
Binding affinity (kcal/mol)
Amino acids involved in interaction
H-bonds
π-bonds
4A
IL-2
−7.680
ARG A 38 THR A 41 ASP B 6 ALA B 13 THR B 14 PHE B 15 LYS B 16 TYR B 20 GLU B 106 THR B 115 GLU B 116 ARG B 117 ILE B 118 TYR B 119 PHE B 121
0
5
IL-6
−7.871
ASN A 62 LEU A 63 PRO A 66 LYS A 67 MET A 68 ALA A 69 PHE A 75 LEU A 166 ARG A 169 SER A 170 GLU A 173 PHE A 174 SER A 177
2
4
Caspase-3
−6.695
THR A 62 SER A 65 TYR A 204 TRP A 206 ARG A 207 SER A 209 PHE A 250 SER A 251 PHE A 256 HOH A 645 HOH A 684 HOH A 696 HOH A 708 HOH A 733 HOH A 736
1
7
Caspase-8
−6.480
SER A 65 TYR A 204 TRP A 206 ARG A 207 ASN A 208 SER A 209 TRP A 214 GLU A 248 SER A 249 PHE A 250 SER A 251
0
1
6A
IL-2
−7.630
ARG A 38 THR A 41 ASP B 6 GLU B 9 PHE B 15 LYS B 16 ALA B 17 TYR B 20 GLU B 113 THR B 115 GLU B 116 ARG B 117 ILE B 118 TYR B 119 PHE B 121
1
1
IL-6
−7.438
GLU A 43 THR A 44 LYS A 47 SER A 48 MET A 50 ARG A 105 PHE A 106 GLU A 107 SER A 108 GLN A 157 TRP A 158 ASP A 161 THR A 164
0
3
Caspase-3
−6.142
SER A 65 TYR A 204 TRP A 206 ARG A 207 ASN A 208 SER A 209 TRP A 214 SER A 249 PHE A 250 SER A 251 PHE A 256 HOH A 645 HOH A 665 HOH A 684 HOH A 696 HOH A 736
1
6
Caspase-8
−6.760
SER A 63 ARG A 64 SER A 65 TYR A 204 SER A 205 TRP A 206 ARG A 207 ASN A 208 SER A 209 TRP A 214 GLU A 248 SER A 249 PHE A 250 SER A 251
1
5
Abbreviation: IL, interleukin.
A study of pharmacokinetic parameters was carried out utilizing QikProp version 4.5 tools to predict the ADME properties of both series (1A–15A) and (1B–15B) and the ranges for the calculated properties of all members, along with their average values are summarized in Table 4. In addition, we also calculated % ABS, number of H-bond acceptors (n–OH), number of H–bond donors (n-OHNH), octanol/water partition coefficients (QPlogPo/w), and Lipinski’s violation. Interestingly, it was found that the % ABS obtained for all members was 100% and the QPlogPo/w prediction was found to be within the accepted range of −2.0 to 6.5. Moreover, all members followed the violated Lipinski parameters.
Table 4
Pharmacokinetic parameters important for oral bioavailability and protein binding parameters of synthesized compounds
MD simulation was performed on the active inhibitor 4A with IL-6 to explore the binding poses in depth. Compound 4A displayed high binding affinity (low docking energy) and stable complex for IL-6 receptor. Therefore, we decided to investigate the influence of compound 4A on the active site of IL-6 receptor. The stability of the system under simulation was assessed using the root-mean-square deviation (RMSD) of atomic positions of the backbone atoms, relative to the starting structures. The binding energy versus time plot indicates the change in conformation and, therefore, the binding energy, with respect to time. As shown in Figures 4 and 5, the RMSD, the binding energy, and the potential energy of IL-6 with compound 4A, including the complex, were computed using MD trajectory frames. Furthermore, we monitored the structural stability of the backbone structure throughout the process through a graphic profile. We observed a dramatic fluctuation in RMSD at 1,000 ps (time), whereas no significant fluctuations were observed after 1,500 ps (time), and it achieved an almost steady state, indicating the stability of the backbone structure with the ligand at about 1,500 ps time in the MD simulation. The binding energy and the potential energy of the complex were estimated as a function of time, which indicated that the potential energy (kJ/mol) did not show further fluctuations after 100 ps, whereas the average complex binding energy was observed to be about −0.9 kg/mol. The fluctuations in the residue of the backbone structure are illustrated in Figure 4. In summary, we examined these data, where we found the structural stability of compound 4A with the active site domain of IL-6 receptor.
Figure 4
Docking complex of 4A with the IL-6 receptor.
Notes: Structural conformational changes (A) before MD simulation and (B) after MD simulation: backbone of the active site domain complex, which shows the contraction of ligand with amino acids residue.
Abbreviations: IL, interleukin; MD, molecular dynamic.
Figure 5
The stability profile of ligand and protein complex under the MD simulation.
Notes: (A) Average RMSD versus time graph that shows the convergence of simulated structure toward an equilibrium state with respect to the reference structure (initial structure). (B) Potential energy of complex versus time graph that shows the stability of ligand and protein complex and (C) binding energy of complex versus time graph that also shows the stability of ligand and protein complex.
We have successfully devised a highly proficient and operationally simple metal-free, one-pot MDR for obtaining a series of novel, hitherto unreported, 5H-benzo[h]thiazolo[2,3-b] quinazoline (1A–15A) and indeno[1,2-d]thiazolo[3,2-a] pyrimidine (1B–15B) analogs displaying potent anticancer activity against Hep-G2 cells as an alternative approach for the treatment of HCC. Considering this library of novel compounds, we concluded that 5H-benzo[h]thiazolo[2,3-b] quinazoline, along with a substituted phenyl ring (hydrophobic side chain) establishes an important pharmacophoric structure, and the R1, R2, and R3 positions of the phenyl ring, as well as the X1 and X2 positions of the tetralone ring system, are the key reactive sites that could be modified with various groups to elicit greater antitumorigenic potential.We have also demonstrated that substitutions with more electronegative groups (–OH, –OCH3) on the hydrophobic side chain directly linked to thiazolo[2,3-b]quinazoline led to active members 4A and 6A, eliciting enhanced anti-tumorigenic activity, with GI50 <10 μg/mL, which was confirmed by docking analyses. Additionally, 3-methoxy-4- hydroxyphenyl–substituted 5H-benzo[h]thiazolo[2,3-b] quinazoline led to 4A, which displayed excellent antitum-origenic activity among a library of 30 novel synthesized compounds with values of GI50 <10 μg/mL. Various computational approaches demonstrated effective oral absorption and protein binding. These compounds, therefore, might be stable in some form of pharmaceutical dosage. Additionally, MD simulation supported our hypothesis regarding the stability of compound 4A with IL-6 protein during the simulation run. The results from the MD simulation run of the active inhibitor showed very less fluctuation with the active site domain of IL-6 and it achieved an almost steady state. We concluded that the compounds would be bound stably to IL-6. Further studies on the biologic activities of the potent compounds are underway in our laboratory.This approach, using a one-pot, multicomponent reaction sequence involving a domino Knoevenagel condensation/Michael-type addition followed by an intramolecular cyclization, where the desired molecules are obtained in a one-flask domino manner with atom and step economy in impressive yields (up to 86%) from readily available and low-priced starting materials, is a resource-effective and desirable route. Additionally, the target compounds were also synthesized by conventional two-step reactions, which were compared with this novel approach on the basis of obtained yields (40%–50%). Finally, we conclude that this newly developed method is both more useful and more profitable in each and every aspect, compared with the conventional route. Due to the importance of these two core structural motifs, that is, 5H-benzo[h]thiazolo[2,3-b]quinazoline and indeno[1,2-d] thiazolo[3,2-a]pyrimidine, especially in the areas of pharmaceutical and medicinal chemistry, we suggest that the protocols that we have outlined should open up new avenues of investigation, with enormous implications for achieving diversity in chemical synthesis.
Authors: Xiangfeng Chen; Youguo Du; Huanliang Sun; Feidong Wang; Lingsheng Kong; Min Sun Journal: Bioorg Med Chem Lett Date: 2013-12-25 Impact factor: 2.823
Authors: P Skehan; R Storeng; D Scudiero; A Monks; J McMahon; D Vistica; J T Warren; H Bokesch; S Kenney; M R Boyd Journal: J Natl Cancer Inst Date: 1990-07-04 Impact factor: 13.506
Authors: Pranesh Kumar; Ashok K Singh; Vinit Raj; Amit Rai; Amit K Keshari; Dinesh Kumar; Biswanath Maity; Anand Prakash; Sabyasachi Maiti; Sudipta Saha Journal: Int J Nanomedicine Date: 2018-02-16
Authors: Amit K Keshari; Ashok K Singh; Umesh Kumar; Vinit Raj; Amit Rai; Pranesh Kumar; Dinesh Kumar; Biswanath Maity; Sneha Nath; Anand Prakash; Sudipta Saha Journal: Drug Des Devel Ther Date: 2017-10-13 Impact factor: 4.162