Abstract

Quinoline-based molecules are major constituents in natural products, active pharmacophores, and have excellent biological activities. Using 2H-thiopyrano[2,3-b]quinoline derivatives and CB1a protein (PDB ID: 2IGR), the molecular docking study has been revealed in this article. The study of in silico molecular docking analysis of such derivatives to determine the binding affinity, residual interaction, and hydrogen bonding of several 2H-thiopyrano[2,3-b]quinolines against CB1a is reported here. The current work demonstrated that 2H-thiopyrano[2,3-b]quinoline derivatives could be effective antitumor agents to produce potent anticancer medicines in the near future.

1. Introduction

Cancer is marked by the uncontrolled, rapid, and pathological growth of improperly altered cells and used in a broader sense. Despite considerable advancements in cancer therapy in recent times, cancer remains the world’s second biggest cause of death, trailing only cardiovascular (CVS) diseases [1]. Chemotherapeutic drug resistance is still a major issue in the fight against cancer; in addition, chemotherapy is further hampered by a lack of selectivity. Anticancer medications, in general, kill both healthy and cancerous cells, and they frequently have major side effects. Many efforts have been undertaken to find new chemotherapeutic medicines with low side effects and to develop safe and effective strategies to treat this disease [2]. Among all cancer types, cervical cancer has the world’s fourth most frequent malignancy among women, with the third highest fatality rate [3]. Cervical cancer alone accounts for about 12% of all cancers in women worldwide, according to a WHO (World Health Organization) report [4], and it is more common in developing nations. Until now, chemotherapeutic medication treatment for cervical cancer has been associated with a poor prognosis and a slew of negative side effects [5]. Therefore, scientists are working for the development of new safer chemicals having potent anticancer activity. For the identification of efficient cellular targets [6], early detection is required which improves the efficacy of the available chemotherapeutic approaches in the field of cancer research [7, 8].

A variety of substituted quinoline derivatives are widely known for their crucial role in the design of innovative pharmacological moieties for therapeutic use, as evidenced by the large number of commercially available medications that include this heterocycle [9]. Quinoline compounds are pharmacologically important and have piqued the interest of chemists and biologists alike due to their wide range of biological actions [10]. Quinoline entities are fundamental building blocks for many naturally occurring chemicals; particularly, bioactive quinoline heterocycles are found in alkaloids of quinoline derived from diverse plant families such as Fumariaceae, Papavaraceae, Rutaceae, and Berberidaceae [10, 11]. Many writers have reviewed the wide range of biological actions of quinoline compounds [12]. Insecticidal [13], antimicrobial [14], antimalarial [15], antiamoebic [16], analgesic [17], vasorelaxing [18], antidiabetic [19], antimycobacterial [20], anticancer [21], anti-inflammatory [22], antihypertensive [23], antiulcer [24], and anti-HIV (human immunodeficiency virus) [25]. Figure 1 shows some clinically used pharmacological based compounds having quinoline heterocyclic core.

The quinoline core is a significant structural unit in many naturally occurring compounds [2629], with interesting biological activity, and many conventional medications. Quinoline core having substitutions at different positions has shown substantial anticancer efficacy against a variety of targets, including topoisomerase I [30], tubulin [31], protein kinase [32], and so on. Some well-known quinolone-based anticancer medicines have been approved by the FDA (Food and Drug Administration) e.g., camptothecin, irinotecan, topotecan, and others. Furthermore, few such quinoline-based anticancer drugs are still in clinical testing phase such as bosutinib, lenvatinib, cabozantinib, and farnesyltransferase inhibitors (tipifarnib) [33].

Several well-known compounds having 2H-thiopyrano[2,3-b]quinolines can be synthesized in a highly efficient manner and create one C-S and one C-C bond with a high extraction economy, likely via a Michael addition, intramolecular aldol-dehydration domino process. It is common knowledge that combining pharmacophore units with different biological activities results in a new hybrid entity with higher biological processes and efficacy than the parent medications. Based on the previous discoveries, unique bioactive hybrid compounds based on quinoline having better anticancer efficacy have been developed nowadays.

In this article, we have reported the binding affinity and several interactions of compounds 14 (as shown in Figure 2) with the anticancer-based protein CB1a (CB1a, a novel anticancer peptide derived from natural antimicrobial peptide cecropin B: PDB ID: 2IGR). The values of binding affinity for 1–4 lies between −5.3 and6.1 Kcal/mol. The interaction also shows that compound 4 shows highest binding affinity (−6.1 Kcal/mol) as compared to other derivatives of thiopyrano[2,3-b]quinolines. For the docking study, we used several software packages such as AutoDock vina 4, discovery studio, and protein-ligand interaction profiler. These software packages help in the study and analysis of docking position, docking size, binding affinity, energy range as well as exhaustiveness. Several amino acids are attached with ligands (thiopyrano[2,3-b]quinoline), and its analogues show several interactions; some are as follows: ILE A-8, ILE A-18, LYS A-7, LYS A-8, LYS A-10, LYS A-11, LYS A-13, LYS A-16, LYS A-25, VAL A-14, LYS A-16, LYS A-26, PHE A-15, TRP A-12, TRP A-25, LYS A-10, and GLU A-9.

2. Materials and Methods

Synthesis of Thiopyrano[2,3-b]quinolines.

Meth Cohn et al. synthesized 2-chloro-3-carbaldehyde quinoline derivatives in high yield by using acetanilide using phosphoryl chloride solution under the presence of dimethyl formamide (DMF) at 70–80°C [34].

Srivastava and Singh et al. synthesized 3-formyl-quinoline-2-thione derivative by using 2-chloro-3-carbaldehydes quinoline derivative using sodium sulphide (Na2S) in the presence of dimethyl formamide (DMF) solvent at room temperature (Scheme 1) [35].

Singh et al. revealed the synthesis of thiopyrano[2,3-(b)]quinoline derivatives by using 3-formyl-quinoline-2-thione derivative and acrylonitrile in the presence of cheap base triethylamine (TEA) and dimethyl formamide (DMF) at room temperature [36]. It is a rapid and efficient one-pot reaction for the synthesis of 3-cyanothiopyrano[2,3-(b)]quinoline derivatives (Scheme 2). The whole reaction was synthesized very easy by the Domino Michael addition followed by cyclization.

Kumar et al. demonstrated the synthesis of 3-nitro-2-phenyl-2H-thiopyrano[2,3-b] quinoline derivatives using 2-mercaptoquinoline-3-carbaldehyde and substituted trans-β-nitrostyrenes in presence of triethylamine (TEA) at 100°C (Scheme 3) [37].

3. Molecular Docking Studies

Cecropins are a category of antimicrobial peptides abundantly available in Hyalophora cecropia’s immunological hemolymph. The N-terminal regions of native cecropins include basic residues, while the C-terminal parts have hydrophobic residues. Cecropin B (CB) has the most potent antibacterial properties and to examine the effects on cells and synthetic liposomes, CB derivatives cecropin B1 (CB1) and cecropin B3 (CB3) were developed. CB1 was made by replacing the C-terminal segment with CB’s N-terminal sequence, while CB3 was made by replacing the N-terminal segment with CB’s C-terminal sequence. CB and its equivalents have been shown in previous research to break membranes, and some of them can destroy cancer cells.

Preliminary SAR interpreting has been conducted using various methods such as molecular docking modelling and docking studies, where the interpretability is considerably better explicit. In this approach, pharmacophore techniques can provide a significant benefit and a clearer look at the structural elements that contribute to the structure activity relationship (SAR) [38]. In this study, pharmacophore generation was discovered using AutoDock vina 4 (The Scripps Research institute) [39] and discovery studio.

The scoring of ligands (binding affinity, ligand internal energy, and distance) is described in the result and discussion section. Using the default settings, the “protein-ligand interaction profiler” [40] calculated the binding affinity of ligands in stable ligand-protein complexes required for ligand binding with the receptor.

The molecular docking study shows that the value of dielectric constant is −0.1465 and the binding spacing is 0.375. The default setting of exhaustiveness is 8, and the RMSD values are calculated relative to the best mode and use only movable heavy atoms. Two variants of RMSD metrics are provided, rmsd/lb (RMSD lower bound) and rmsd/ub (RMSD upper bound). In addition to this, the energy of binding modes in the output is 4. The “protein-ligand interaction profiler” (PLIP) was used to calculate the interaction of protein and ligands in the stable ligand-protein complex required for ligand binding to the receptor.

Based on molecular docking study of four compounds having thiopyrano[2,3-b]quinoline core 14 against CB1a (PDB ID: 2IGR). The docking result is given in Table 1.

In general, the IC50 of CB1a, CB, melittin, magainin II, as well as the cancer chemotherapeutic agent were determined using cytotoxicity assays on cancer and noncancer cells. CB1a has a stronger cytotoxic activity (lower IC50 values) against leukaemia cells and stomach cancer than CB and magainin II. CB1a and CB are highly effective towards cancerous cells but not normal human cells. CB1a has an 8-fold and 2-3-fold higher cytotoxic activity than CB against AGS and leukaemia cell lines, respectively. As a result, molecular docking study of the synthesized compounds was performed in the current study to explore their binding pattern with a unique anticancer peptide derived from the natural antimicrobial peptide cecropin B (CB1).

Because of its excellent selectivity, CB1a can be developed as a powerful anticancer drug. CB1a > CB1 > CB > CB3 is the order of anticancer activity of CB analogues against tumor cell lines [41]. Therefore, for this study, CB1a, a novel anticancer peptide derived from natural antimicrobial peptide cecropin B, was used. The synthetic compounds (14) were found to have binding affinity from −5.3 to6.1 kcal/mol (Table 1) with the best result achieved using compound 4 (−6.1 kcal/mol). The hydrogen bond, residual interaction, and pi-pi interaction of the four compounds were summarized in (Table 1).

The compounds (1) show similar residual interactions with amino acid residues PHE-15, ILE A-8, ILE A-18, LYS A-7, LYS A-10, LYS A-11, VAL A-14, and TRP A-12 as shown in Figure 3. The in-silico interaction results match with the in vitro analysis of the synthesized compounds against the structure of CB1a, among which compound 1 shows low binding affinity as comparison to other derivatives with the value of −5.3 kcal/mol.

The compounds (2) show similar residual interactions with amino acid residues PHE A-15, TRP A-12, LYS A-16, and LYS A-26 as shown in Figure 4. The data also show that there is additional hydrogen bonding interaction with amino acid residue GLU-32. The in-silico interaction results match the in vitro analysis of the synthesized compounds against the structure of CB1a. The binding affinity of compound 2 is found as −5.5 kcal/mol.

The compounds (3) also show similar residual interactions with amino acid residues PHE A-15, ILE A-8, ILE A-18, LYS A-7, LYS A-10, LYS A-11, LYS A-13, VAL A-14, TRP A-12, TRP A-25, and GLU A-9 as shown in Figure 5. It shows additional hydrogen bonding with amino acid residue. The in-silico interaction results match the in vitro analysis of the synthesized compounds against the structure of CB1a, a novel anticancer peptide derived from natural antimicrobial peptide cecropin B among which compound 3 shows low binding affinity as comparison to other derivative and the value is −5.9 kcal/mol.

The compounds (4) show similar residual interactions with amino acid residues PHE A-15, ILE A-8, ILE A-18, LYS A-7, LYS A-11, VAL A-14, LYS A-16, LYS A-26, LYS A-10, and GLU A-9 as shown in Figure 6. The in-silico interaction of compound 4 results matches the in vitro analysis of the synthesized compounds against the structure of CB1a, which shows high binding affinity and good activity as compared to other derivatives since the value is −6.1 kcal/mol.

4. Conclusion

The current work demonstrates a useful study to find the possible effective antitumor agents having thiopyrano[2,3-b]quinoline derivatives, which could produce potent anticancer medicines in near future. To study the interaction of ligand and protein, several software packages were used in this work such as AutoDock vina 4, discovery studio, and protein-ligand interaction profiler (PLIP). The molecular docking of compounds having thiopyrano[2,3-b]quinoline core and its derivatives against CB1a (PDB ID: 2IGR) were studied in this article. The in silico molecular docking study of these compounds reveals that these frameworks have potential bioactivity having high binding affinity, adequate residual interaction, and hydrogen bonding interaction against the protein CB1a. The binding affinity value for thiopyrano[2,3-b]quinoline and its derivatives 14 is −5.3 to6.1 Kcal/mol. The study also revealed that several amino acids show interaction with thiopyrano[2,3-b]quinoline ligand and their analogues. Some of the identified amino acids are ILE-8, ILE-18, LYS-7, LYS-11, VAL-14, LYS-16, LYS-26, PHE-15, LYS-15, TRP-12, TRP-25, LYS-10, and GLU-9. In general, the reported core already exhibits antibacterial activity and can be further refined to serve as anticancer compounds in near future.

Abbreviations

WHO:World Health Organization
HIV:Human immunodeficiency virus
FDA:Food and Drug Administration
PDB:Protein data bank
PLIP:Protein-ligand interaction profiler
TEA:Triethylamine
DMF:Dimethyl formamide
SAR:Structure activity relationship
RMSD:Root mean square deviation.

Data Availability

The data are available at the Department of Applied Chemistry, Amity School of Applied Sciences, Amity University Madhya Pradesh, Gwalior, Madhya Pradesh-474 005, India.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Supplementary Materials

All the interaction details and necessary data are available at the supplementary data. (Supplementary Materials)