The newly synthesized coumarin derivative with dopamine, 3-(1-((3,4-dihydroxyphenethyl)amino)ethylidene)-chroman-2,4-dione, was completely structurally characterized by X-ray crystallography. It was shown that several types of hydrogen bonds are present, which additionally stabilize the structure. The compound was tested in vitro against different cell lines, healthy human keratinocyte HaCaT, cervical squamous cell carcinoma SiHa, breast carcinoma MCF7, and hepatocellular carcinoma HepG2. Compared to control, the new derivate showed a stronger effect on both healthy and carcinoma cell lines, with the most prominent effect on the breast carcinoma MCF7 cell line. The molecular docking study, obtained for ten different conformations of the new compound, showed its inhibitory nature against CDKS protein. Lower inhibition constant, relative to one of 4-OH-coumarine, proved stronger and more numerous interactions with CDKS protein. These interactions were carefully examined for both parent molecule and derivative and explained from a structural point of view.
Ministarstvo Prosvete, Nauke i Tehnološkog RazvojaOI174028OI172016OI173051OI172015OI1720401. Introduction
Cell structures can be significantly influenced and damaged by sustained oxidative stress which is considered to be a major cause in the pathogenesis of many, if not all, diseases. Cancer is one of the most prominent causes of death in the modern world. Cancer initiation and progression phases have been closely related to oxidative stress which increases somatic mutations, neoplastic transformation, and generally genome instability [1, 2]. In the past several years, there have been an increased number of articles that connect oxidative stress with the development of various cancers, although the actual links are a matter of dispute [1, 2].
Nowadays, there is an increased interest in new drug discovery, to adequately promote the survival rate against cancer. Presently existing anticancer drugs are insufficient due to many side effects, nonselectivity, resistance, etc. Many of the naturally occurring molecules possessing antioxidant activity exhibit a synergistic effect when combined with other naturally occurring, or synthetic, molecules. This can lead to the creation of new types of a biological rationale for the treatment of different cancer types or their use as adjuvants with conventional therapeutic regimens. The present paper presents the design of a new coumarin/catecholamine derivative [1, 2].
Coumarins and their derivatives are widely present family of molecules in nature. They can be accumulated in fruits [3], vegetables [4, 5], trees [6], seeds [7], and vines. Coumarins have important biological activities, some of which include regulation of growth, control of respiration [8], defense against herbivores and microorganisms [9], and hormonal and signaling role [10]. The common structural elements are fused pyrone and benzene rings with a carboxyl group on the first ring. The structural diversity of coumarins allows different pharmacological properties: antibacterial, antifungal, antioxidant, and cytotoxic [11–13]. The synthetic coumarins also have great potential as drugs for the treatment of neurodegenerative [14, 15] and microbial [16, 17] diseases, HIV [18, 19], and cancer [20, 21]. They also have been used as regulators of reactive radical species [22, 23]. Several reactive derivatives of coumarin with various molecules were already synthesized, and their biological reactivity explained experimentally and theoretically [24–26].
Dopamine, on the other side, belongs to the important group of catecholamines, which function as hormones and neurotransmitters in the peripheral endocrine and central nervous system [27–30]. The other functions include control of movement, reward processing, attention, working memory, emotions, sleeping, and dreaming. Dopamine and other catecholamines, as well as their metabolites, are gaining significant attention, due to the hypothesis that the origin of some, if not all, neurodegenerative diseases lies in oxidative stress and changes at the molecular level [31–34]. More than 20% of oxygen in the body is used in the brain; therefore, a considerable amount of the present fatty acids is exposed to reactive oxygen and nitrogen species [35–38]. The low permeability of the blood-brain barrier limits the use of the external antioxidants, which increases the importance of molecules with good radical scavenging activities that are already present in the body [37]. The catechol structure, as well as various side groups, has been proven, in different experimental and theoretical studies, to enhance the antiradical potency of molecules [39–43]. Besides being good antiradical agents, catecholamines and their analogs also exhibited as potential antitumor agents [44–46]. The structure-activity studies show that catechol moiety is also the essential part of these molecules influencing their antitumor activity [45, 47, 48]. In the case of dopamine and its analogs, it has been concluded that the aminoethyl group could be replaced by methyl or aminomethyl groups without a noticeable change in activity towards P-388 leukemia cells [49]. One of the assumptions is that the formation of reactive quinone is responsible for the cytotoxic and genotoxic activity of dopamine [50].
This article presents data on the reactivity of novel coumarin/dopamine derivative towards various tumor and healthy cells lines, with special emphasis on the synergistic effect of these two molecules. The synthesis is carried out assuming that there is no formation of a quinone, which is harmful to healthy cells, and that the catechol unit is preserved and combined with the coumarin. The crystallographic structure is obtained for the synthesized derivative and explained in details. The molecular docking study is performed in order to gain better insight into the mechanism of the biological activity of new compound and possible differences in coumarin-protein and new derivative-protein interactions. The protein selected for molecular docking is cyclin-dependent kinases (CDKs), for both 4-hydroxycoumarin, as a parent molecule, and synthesized derivative. CDKs plays an important role in the control of cell division [51], and its deregulation can lead to the development of human diseases including cancer [52, 53]. Thus, the suppression of its activity could be a promising molecular route for anticancer therapy [54].
2. Materials and Methods2.1. Substances
Dopamine hydrochloride, 4-hydroxycoumarin (4-OH-coumarin), methanol, toluene, and 96% ethanol were purchased from Sigma-Aldrich. The starting compound, 3-acetyl-4-hydroxy-coumarin, was obtained as explained in [55].
2.2. Spectral Analysis
The vibrational spectra were recorded on the Perkin-Elmer Spectrum One FT-IT spectrometer in the range between 4000 and 400 cm-1. The KBr pellet technique was used. The NMR spectra (1H and 13C) were recorded on a Varian Gemini 200 spectrometer in CDCl3 as a solvent and with TMS as an internal standard. The elemental microanalysis for C, H, and N was performed on the Vario EL III C, H, and N Elemental Analyzer. The mass spectra were obtained on a 5973 Mass spectrometer (Agilent, Santa Clara, CA) with MS quadruple, temperature 150°C, and a mass scan range of 40–600 amu at 70 eV.
2.3. Cell Culture
In this study, the following cell lines were used: healthy human keratinocyte HaCaT (AddexBio T0020001), cervical squamous cell carcinoma SiHa (ATCC®, HTB-35™), breast carcinoma MCF7 (ATCC® HTB-22™), and hepatocellular carcinoma HepG2 (ATCC® HB-8065™). All cell lines were grown in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum (FBS), 4500 mg/L glucose, and 1x antimycotic/antibiotic (all from Invitrogen™, USA). The cells were maintained at 37°C in 5% CO2.
2.4. Cell Viability Assay
In order to assess whether selected compounds influence the viability of cell, the MTS Cell Proliferation Assay was applied, a colorimetric method for sensitive quantification of viable cells. The method is based on the reduction of MTS tetrazolium compound by viable cells generating a colored formazan product. This conversion is thought to be carried out by NAD(P)H-dependent dehydrogenase enzymes in metabolically active cells. The formazan dye produced by viable cells can be quantified by measuring the absorbance at 490–500 nm.
Cells (5 × 103 for Siha and MCF7 and 1 × 104 for HaCaT and HepG2/per well) were seeded in 96-well plates a day before treatment and then treated with various concentrations of 4-OH-coumarin or derivative 3 (100, 300, and 500 μM) for 48 hours. After 48 h, the effect of these treatments was monitored on cell’s viability using MTS Cell Proliferation Assay (Promega CellTiter 96® AQueous One Solution Cell Proliferation Assay) and colorimetric quantification was done using a plate reader (Plate Reader Infinite 200 pro, Tecan).
2.5. Statistical Analyses
Statistical analyses were performed with SPSS statistical software (version 20). The data represent the means±SEM from at least three independent experiments. Statistical analyses were performed by Student’s t-test, and p value ≤ 0.05 was considered significant.
2.6. Theoretical Methods
Software package AutoDock 4.0 was used for molecular docking simulations [56]. The three-dimensional crystal structure of CDK protein was procured from the Protein Data Bank (PDB ID: 1KE9). Preparation of protein for docking simulations was done by removing the cocrystallized ligand, water molecules, and cofactors in the Discovery Studio 4.0. The calculations of Kollman charges and adding of the polar hydrogen were performed using the AutoDockTools (ADT) graphical user interface. The optimization of coumarin and coumarin derivate structures was performed at B3LYP-D3BJ/6-311++G(d,p) level of the theory [57–61] in the gas phase. For this purpose, Gaussian 09 Program package [62] was used. The structure of coumarin derivate was taken from crystallographic data. This level of theory was proven to reproduce well the crystallographic structure and experimentally obtained spectra for similar compounds [24, 25, 26]. The docking simulations were completed independently. During these simulations, the protein retained a rigid structure in the ADT, while both docked compounds possessed flexibility. All bonds of the compounds were treated as rotatable. For calculations of partial charges, the Geistenger method was selected. The Lamarckian Genetic Algorithm (LGA) method was used in all calculations for protein-ligand flexible docking simulation [63]. The grid boxes with grid center 4.0 × 50.3 × 11.5 of CDK protein covered all of the protein binding sites and enabled free motion of ligands. The interactions between CDKs and the corresponding ligands were estimated and analyzed, for several most stable conformations, through the positions of the amino acids and the type of interaction that is achieved.
2.7. General Procedure for the Synthesis of 3-(1-((3,4-Dihydroxyphenethyl)amino)ethylidene)-chroman-2,4-dione
The synthesis of new coumarin derivative is presented in Scheme 1. The new compound (designated as compound 3 in Scheme 1) was obtained by refluxing and mixing of the reaction mixture (3-acetyl-4-hydroxicumarine (0.0014 mol; 0.3 g), dopamine hydrochloride (0.0014 mol; 0.23 g), and equimolar amounts of triethyl amine (0.0014 mol; 0.15 g) in 50 ml methanol for 3 h. The progress of reactions was monitored by TLC (toluene: acetone = 7: 3). When the reaction was completed, the mixture was cooled to room temperature. The obtained white crystals were filtered, air-dried, and recrystallized from methanol.
Synthesis of 3-(1-((3,4-dihydroxyphenethyl)amino)ethylidene)-chroman-2,4-dione.
Yield, 0.321 g (64.45 %), Anal. Calcd. for C19H17O5N (Mr = 339.32) %: C, 67.25; H, 5.05; N, 4.12. Found: C, 67.00; H, 5.10; N, 4.16.
IR (KBr, cm-1): 3305 (NH and OH), 2921, 2857 (CH), 1668 (C=O), 1605, 1569, 1531, (C=C), 1118 (C–O).
2.8. X-ray Data Collection and Structure Refinement
A summary of X-ray diffraction experiment and structure refinement for 3·MeOH is given in Table 1. The data collection was performed on an Oxford Diffraction Xcalibur Gemini ultra-diffractometer equipped with an AtlasS2 CCD detector using CuKα radiation. CrysAlis PRO 1.171.39.35c [64] was used for data collection, cell refinement, data reduction, and absorption correction. The structure was solved by SHELXT [65] and subsequent Fourier syntheses using SHELXL [66], implemented in WinGX program suit [67]. Anisotropic displacement parameters were refined for all nonhydrogen atoms. The hydrogen atoms bonded to nitrogen, oxygen atoms were found in the Fourier maps and refined freely, and aromatic and aliphatic carbon-bonded hydrogen atoms were placed in the calculated positions and refined riding on their parent C atoms with corresponding C–H distances and Uiso(H) = 1.2 or 1.5 Ueq(C), respectively. The analysis of bond distances and angles was performed using SHELXL and PLATON [68]. DIAMOND [69] was used for molecular graphics.
Crystal data and structure refinement of 3·MeOH.
Compound
3·MeOH
Empirical formula
C20H21NO6
Formula weight
371.38
Temperature
120(2) K
Wavelength
1.54184 Å
Crystal system
Triclinic
Space group
P-1
Unit cell dimensions
a = 7.9449(3) Å
α = 82.538(3)°
b = 8.8742(3) Å
β = 82.296(3)°
c = 13.1920(5) Å
γ = 70.039(3)°
Volume
862.79(6) Å3
Z; density (calculated)
2; 1.430 g·cm-3
Absorption coefficient
0.883 mm-1
F(000)
392
Crystal shape, color
Prism, white
Crystal size
0.311 × 0.201 × 0.072 mm3
θ range for data collection
3.395–67.361°
Index ranges
-9 ≤ h ≤ 6, -10 ≤ k ≤ 10, -15 ≤ l ≤ 15
Reflections collected/independent
8977/3086 [R(int) = 0.0226]
Absorption correction
Analytical
Max. and min. transmission
0.938 and 0.841
Data/restraints/parameters
3086/0/262
Goodness-of-fit on F2
1.034
Final R indices [I > 2σ(I)]
R1 = 0.0342, wR2 = 0.0876
R indices (all data)
R1 = 0.0398, wR2 = 0.0924
Largest diff. peak and hole
0.226; -0.255 e.Å-3
3. Results and Discussion3.1. Chemistry
The synthesis of coumarin derivative is presented in Scheme 1. The structure of synthesized compound 3 was determined by means of elemental, spectral (IR, 1H NMR, and 13C NMR), and X-ray structural analysis (in the form of 3·MeOH).
The formation of compound 3 was confirmed by IR spectra, with the presence of bands positioned at 3305 cm-1 assigned to NH and OH group vibrations. Also, stretching vibrations corresponding to the C=O and C–O groups from 2,4-dioxochroman moiety were identified at 1668 and 1118 cm-1, respectively.
In 1H NMR spectrum, the singlet positioned at 2.55 ppm was assigned to protons on C2′ atom. The resulting signals of aromatic protons of the 2,4-dioxochroman part were in the range from 7.25 to 7.91 ppm. Aromatic protons belonging to the phenyl group of dopamine part were detected in the range from 6.5 to 6.67 ppm. The protons of the methylene C2″–H and C1″–H groups were identified as triplets at 2.80 and 3.6 ppm. Signals of protons of the phenyl OH and enamine NH group were identified as broadened singlets at 8.81 and 13.66 ppm.
The 13C NMR spectra of compound 3 indicated the presence of aromatic carbon atoms of the dopamine part in the range from 115.87 to 145.46 ppm. The signals of carbon atoms of coumarine moiety were identified in the range from 96.18 to 179.60 ppm. The carbons of lactone (C–2) and ketone (C–4) showed resonances at 162.13 and 179.60 ppm. Signals at 18.35, 34.11, and 45.71 ppm were assigned to C2′, C1″, and C2″ carbons, respectively.
3.2. Crystallographic Structure
X-ray structure analysis revealed that 3 crystallizes in the triclinic P-1 space group. Its molecular structure is formed by the molecule of 3, consisting of bicyclic coumarine and 3,4-dihydroxyphenethyl fragments joined by aminoethylidine chain, and by a solvated molecule of methanol tied with the mentioned molecule by a hydrogen bond (Figure 1). Within the molecule, the dihedral angle between the planes of coumarin fragment and the 3,4-dihydroxyphenyl ring is 61.46(3)°.
Molecular structure with an atom numbering scheme of 3·MeOH. Displacement ellipsoids are drawn at 50% probability; hydrogen bonds are shown as red dashed lines.
A typical structural feature of that type of compounds is an intramolecular N − H···O hydrogen bond which forms a six-membered ring with S(6) graph-set motif [70] and thus the molecule occurs in a ketoamine tautomeric form. Considering the O3=C4−C3=C1′−N1−H1N1-conjugated bond ring system created, owing to the intramolecular hydrogen bond formation, the equalization of the C3−C4 and C3=C1′ (both 1.442(2) Å) bond lengths (Table S1) is observed although the bonds are formally single and double, respectively. This can be explained by the π-electron delocalization within the system, and thus, we may conclude that the above hydrogen bond can be classified as a resonance-assisted hydrogen bond [71]. Interestingly, in the related 3-(1-((m-toluidine)amino)ethylidene)-chroman-2,4-dione and in 3-(1-(2-hydroxyethylamino)ethylidene)-chroman-2,4-dione [72] compounds, C3−C4 bond (1.434(2) and 1.430(3) Å, respectively) is even slightly shorter than C3=C1′ bond (1.436(2) and 1.437(3) Å, respectively). The elongation of C4=O3 (1.248(2) Å) bond, which is markedly longer than C2=O2 (1.223(2) Å) bond not involved in a strong hydrogen bond, and shortening of C1′−N1 (1.309(2) Å) bond in 3 can also be observed. Very similar bond lengths were observed in related compounds, like 3-(1-((o-toluidine)amino)ethylidene)-chroman-2,4-dione [72], 3-(1-(phenylamino)ethylidene)-chroman-2,4-dione [24], 3-(1-(3-hydroxypropylamino)propylidene)chroman-2,4-dione [73], 2-(1-(2,4-dioxochroman-3-ylidene)ethylamino)-3-methylbutanoate [74], 3-[(1-benzylamino)ethylidene]-2H-chromene-2,4(3H)-dione [75], or 3-[1-((2-hydroxyphenyl)amino)ethylidene]-2H-chromene-2,4(3H)-dione compound [76]. All other bond lengths and angles (Table S1) in the molecule of 3 are within normal ranges [77].
Except for above discussed N1−H1N1···O3 intramolecular hydrogen bond, due to which the exocyclic C3=C1′ double bond has an E geometry, the molecules of 3 are stabilized in the solid state by intermolecular O−H···O and C−H···O hydrogen bonds. Further stabilization of the solid state structure occurs when molecule of methanol occupies empty space between molecules of 3 and is tied with the molecules by a pair of O−H···O hydrogen bonds (Table 2). Due to these bonds, the molecules of 3 and methanol are tied to form chains parallel with the [011] direction (Figure S2). These chains are further connected by π-π interactions between pyran-2,4-dione (py) and phenyl (ph) rings of coumarin moieties in adjacent chains into a 2D structure parallel with (011) (Figure 2). These π-π interactions are characterized by Cgpy···Cgpyiii and Cgpy···Cgphiii centroid-centroid distances of 3.6479(1) and 3.7569(1) Å, respectively (iii=1–x,1–y,1–z).
Molecular packing of 3·MeOH showing π-π interactions (black dashed lines) between neighboring chains formed by hydrogen bonds (red dashed lines). Hydrogen atoms not involved in hydrogen bonds are omitted for clarity.
3.3. Antitumor Activity
The effect of both 4-OH coumarin and its derivate 3 on selected cell lines was monitored and compared to the effect of DMSO that was used as a control (vehicle). As presented in Figure 3, 4-OH-coumarin exhibited a mild effect on cell’s viability in all tested cell types with the maximal effect on breast carcinoma cell line MCF7, inducing 20% reduction in viability of these cells 48 h after treatment at 500 μM concentration (p=0,023). On the other hand, its derivate 3 had, in comparison with the control sample, a significantly stronger effect on both healthy and carcinoma cell lines with the most prominent effect on the MCF7 cell line. In particular, 100 μM concentration of derivative 3 led to the reduction of cell’s viability to approximately 75% in HaCat cells (p=0,023), 72% in SiHa cells (p=0,008), 49% in MCF7 cells (p=0,008), and 62% in HepG2 cells (p=0,005). The reduction effect increased with the concentration of the derivative 3 and had the greatest effect on the MCF7 cells, with a reduction percentage of 43% at the concentration of 300 μM (p=0.001). The statistical significance is also observed when the effect of 3 is compared with the effect of 4-OH–coumarin (Figure 3, p values are indicated as asterisks). Taken together, the higher cytotoxic effect of derivate 3 is observed, both for healthy and for carcinoma cells, compared to 4-OH-coumarin, where the effect on the carcinoma cells was somewhat stronger.
Cell viability assay (MTS assay) on HaCaT, SiHa, MCF7, and HepG2 cells performed 48 h after treatment with either 4-OH-coumarin or compound 3. Relative cell viability of cells treated with selected compounds was calculated as a percentage of DMSO-treated cell viability that was set as 100%. Data are presented as the means±S.E.M. (standard error mean) of at least three independent experiments performed in triplicate for each concentration. Mean values of relative cell viability were compared with Student’s t-test, and p values are presented as ∗p≤0.05, ∗∗p≤0.01, and ∗∗∗p≤0.001. Each color corresponds to a bar presented on the histogram.)
It is clearly demonstrated that coumarin derivate 3 exhibits a cytotoxic activity against all analyzed cell types in vitro. The observed effect is somewhat stronger against carcinoma cell lines compared to healthy keratinocytes, with the most prominent effect on breast carcinoma cell line MCF7. Antitumor activity of natural and synthetic coumarin derivatives has been previously reported by various authors including high cytotoxicity against ovarian cancer cells, lung carcinoma cells, pancreatic carcinoma, hepatocarcinoma, and breast and colon carcinoma [78–82]. Depending on their structures, coumarins can act on various tumor cells by different mechanisms such as inhibition of the telomerase and protein kinase activities, downregulation of oncogene expression or induction of the caspase-9-mediated apoptosis, and suppression of cancer cell proliferation [82]. A further work is needed to elucidate the mechanism of action for this newly synthesized derivate 3.
3.4. Molecular Docking
The molecular docking studies were performed for the evaluation of the inhibitory nature of examined compounds against CDKS protein. These simulations gave the predicted protein-ligand binding energies and identified the potential ligand binding sites. The structure of the newly synthesized compound was optimized at the B3LYP-D3BJ/6-311++G(d,p) level of theory, based on the crystallographic positions of atoms (Tables S3 and S4). Ten different conformations were analyzed for both investigated compounds. Tables S5 and S6 give values of the estimated free energy of binding and inhibition constant values (Ki), the distance between respective active sites of ligand and amino acids, and pairwise interaction energies (Ei), as well as types of interactions for the investigated models with the lowest docked conformation energies. The most stable conformations are presented in Figures 4 and 5.
Docking positions of 4-OH-coumarin.
Docking positions of coumarin derivate.
The parent compound contains several polar groups, namely, a hydroxy group in position 4, a carbonyl group, and an oxygen atom in the pyrone ring. Since the rest of the structure makes the benzene ring, therefore, the most of interactions, presented in Table S5, include π-alkyl and π-σ hydrophobic interactions with leucine, valine, alanine, arginine, and phenylalanine in various positions [83–85]. These interactions are characterized by low pairwise interaction energy and large atomic distances (≥2.5 Å). Due to a large number of weak noncovalent interactions, all of the models obtained for 4-OH-coumarin have very low binding energy and high value of inhibition constant, which is in accordance with the obtained experimental result.
It is also evident that two types of hydrogen bonds are formed in given models. The first type is the conventional hydrogen bond. Only two bonds of this type, with significant values of pairwise interaction energy, have atomic distances lower than 2 Å. These bonds are formed upon the interaction between CDKs and LEU83 and GLU257. It should be pointed out that in models 3 and 4 there are two more bonds which deserve attention, THR97 and GLU195, respectively. It should be pointed out that the hydroxy group of 4-OH-coumarine behaves as the hydrogen atom donor in the interaction with THR97 (model 2), GLU195 (model 3), and GLU257 (model 9). There are more interactions in which the parent molecule behaves as the hydrogen atom acceptor. Other conventional hydrogen bonds are weak, due to high atomic distances and low pairwise interaction energy. These interactions are established with positively charged amino acids, lysine, histidine, and arginine. The second type of hydrogen bonds is the carbon-hydrogen bond. Predicted values of the interaction energy for this bond type are very low, and the atomic distances much larger (≥3 Å). These bonds are formed between donating groups of HIS268 and ARG260 and hydrogen-acceptor groups of the CDKs.
Table S6 presents molecular docking results for compound 3. When the structures of two investigated compounds are compared, it is notable that only compound 3 possesses a catechol moiety and longer aliphatic chain with a nitrogen atom. This increases the number of possible interactions with amino acids. If results from Tables S5 and S6 are compared, it is obvious that the relative abundance of hydrogen bonds increases with respect to other noncovalent interactions. There are also still hydrophobic π-σ (ILE10 in model 4 and PRO292 in model 8) and hydrophobic π→alkyl (in models 2, 3, 8, and 9) interactions. All these interactions have very low pairwise interaction energies. When these interactions are formed, atomic distances fall into the very wide range of values (3.5-4.40 Å). In spite the fact that these interactions are very weak, they additionally stabilize the structures.
When hydrogen bonds are concerned, there are also two types of bonds, the conventional and carbon-hydrogen bonds. Compound 3 behaves as the hydrogen atom donor with ASP in positions 145 and 86 (model 4), THR97 (model 6), VAL293 (model 8), ALA282, and HIS293 (model 2). It is important to point out that the number of possible hydrogen atom accepting amino acids increases, due to the fact that compound 3 has additional polar groups. The strength of formed conventional hydrogen bonds and interatomic distances are in the same range as in 4-OH-coumarine. Therefore, it can be concluded that the reactivity of coumarin part is conserved and that the additional groups contribute to the increasing number of interactions.
It should be emphasized that both values, the binding energy and inhibition constant, of the compound 3 are lower in comparison with the parent molecule. Moreover, it should be pointed out that the value of inhibition constant for compound 3 is more than six times lower than the one of 4-OH-coumarine. This indicates that compound 3 interacts better with CDKS protein. Obviously, the additional hydroxy groups and an aromatic ring, with an aliphatic chain containing a nitrogen atom, allow new interactions which inhibit the reactivity of the protein. This could be one of the reasons for the increased reactivity of compound 3 towards tumor cells.
4. Conclusion
The coumarin and dopamine derivative, 3-(1-((3,4-dihydroxyphenethyl)amino)-ethylidene)-chroman-2,4-dione, was synthesized under mild conditions. The new compound was analyzed by NMR, IR, microanalysis, and X-ray crystallography. The X-ray analysis showed that the similar structural motifs are present, as with other previously obtained derivatives with aminophenols. Several types of hydrogen bonds, both intramolecular and intermolecular, stabilize the structure within the crystal. The molecule is not planar, but there is the dihedral angle between the planes of coumarin fragment and 3,4-dihydroxyphenyl ring of 61.46(3)°. There are also π→π and stacking interactions within the crystal structure.
The antitumor activity was investigated against healthy and tumor cell lines both for 4-OH-coumarin and its derivative. Compared to 4-OH coumarin, new derivative showed a significantly stronger effect on both healthy and carcinoma cell lines. When treated with 100 μM solution of compound 3, the reduction of cell’s viability was approximately 75% in HaCat cells, 72% in SiHa cells, 49% in MCF7 cells, and 62% in HepG2. The most prominent effect was observed on the breast carcinoma MCF7 cell line.
The molecular docking study was performed in order to better understand the difference in binding between the two investigated molecules and CDK protein. Different interactions are possible due to the presence of polar groups in the coumarin structure. The hydrogen bonds are the strongest interactions observed, in which 4-OH-coumarin can act as the hydrogen atom acceptor and hydrogen atom donor. The most numerous are π→alkyl interactions with various amino acids. The strength of interactions, given as the pairwise interaction energy, is preserved in coumarin derivative. New interactions are established as a result of the presence of additional polar groups: catechol moiety and alkyl chain. The obtained result implies that a number of interactions determine the activity towards investigated protein.
In the end, it can be concluded that the presented results are promising. Future experiments, which would include new cell lines and new coumarin-neurotransmitter derivatives, in addition to all the above levels of testing, will be also supplemented with classical MD calculations in order to better explain the mechanisms of action.
Data Availability
The X-ray structure analysis and molecular docking data used to support the findings of this study are included within the supplementary information file.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
This research was funded by the Ministry of Education, Science and Technological Development of the Republic of Serbia through Grant nos. OI172040, OI172015, OI173051, OI172016, and OI174028. The instruments used for crystallographic measurements are from the ASTRA lab established within the Operation Program Prague Competitiveness (project CZ.2.16/3.1.00/24510).
Supplementary Materials
The Supplementary information file contains additional information on crystallographic data of the newly synthesized molecule, along with selected bond lengths and angles, as well as results of docking analysis with CDKs (type of interaction, atom distances, bonding energy, estimated free energy of binding, and inhibition constants).
HalliwellB.Oxidative stress and cancer: have we moved forward?2007401111110.1042/BJ200611312-s2.0-3384633165017150040ReuterS.GuptaS. C.ChaturvediM. M.AggarwalB. B.Oxidative stress, inflammation, and cancer: how are they linked?201049111603161610.1016/J.FREERADBIOMED.2010.09.0062-s2.0-7804937098720840865MeadJ. A. R.SmithJ. N.WilliamsR. T.Studies in detoxication. 72. The metabolism of coumarin and of o-coumaric acid1958681677410.1042/bj06800672-s2.0-000000215513522576RosselliS.MaggioA. M.FaraoneN.SpadaroV.Morris-NatschkeS. L.BastowK. F.LeeK.-H.BrunoM.The cytotoxic properties of natural coumarins isolated from roots of Ferulago campestris (Apiaceae) and of synthetic ester derivatives of aegelinol20094121701170620120111Atta-ur-RahmanShabbirM.Ziauddin SultaniS.JabbarA.Iqbal ChoudharyM.Cinnamates and coumarins from the leaves of Murraya paniculata199744468368510.1016/S0031-9422(96)00617-62-s2.0-0031080931FullerR. W.BokeschH. R.GustafsonK. R.McKeeT. C.CardellinaJ. H.IIMcMahonJ. B.CraggG. M.SoejartoD. D.BoydM. R.HIV-inhibitory coumarins from latex of the tropical rainforest tree Calophyllum teysmannii var. inophylloide19944161961196410.1016/S0960-894X(01)80543-62-s2.0-0028169099SpinoC.DodierM.SotheeswaranS.Anti-HIV coumarins from calophyllum seed oil19988243475347810.1016/S0960-894X(98)00628-32-s2.0-00325348979934455HatanoT.YasuharaT.FukudaT.NoroT.OkudaT.Phenolic constituents of licorice. II. Structures of licopyranocoumarin, licoarylcoumarin and glisoflavone, and inhibitory effects of licorice phenolics on xanthine oxidase198937113005300910.1248/cpb.37.30052-s2.0-00248319772632045HarborneJ. B.WaltonN. J.BrownD. E.Classes and functions of secondary products from plants1999LondonImperial College Press10.1142/9789812817273_0001BourgaudF.HehnA.LarbatR.DoerperS.GontierE.KellnerS.MaternU.Biosynthesis of coumarins in plants: a major pathway still to be unravelled for cytochrome P450 enzymes200652-329330810.1007/s11101-006-9040-22-s2.0-33845445814OjalaT.RemesS.HaansuuP.VuorelaH.HiltunenR.HaahtelaK.VuorelaP.Antimicrobial activity of some coumarin containing herbal plants growing in Finland2000731-229930510.1016/S0378-8741(00)00279-82-s2.0-0033797710CottigliF.LoyG.GarauD.FlorisC.CausM.PompeiR.BonsignoreL.Antimicrobial evaluation of coumarins and flavonoids from the stems of Daphne gnidium L20018430230510.1078/0944-7113-00036BasileA.SorboS.SpadaroV.BrunoM.MaggioA.FaraoneN.RosselliS.Antimicrobial and antioxidant activities of coumarins from the roots of Ferulago campestris (Apiaceae)200914393995210.3390/molecules140309392-s2.0-63449137987PatilP. O.BariS. B.FirkeS. D.DeshmukhP. K.DondaS. T.PatilD. A.A comprehensive review on synthesis and designing aspects of coumarin derivatives as monoamine oxidase inhibitors for depression and Alzheimer’s disease20132192434245010.1016/J.BMC.2013.02.0172-s2.0-84876107554ZhouX.WangX.-B.WangT.KongL.-Y.Design, synthesis, and acetylcholinesterase inhibitory activity of novel coumarin analogues200816178011802110.1016/J.BMC.2008.07.0682-s2.0-5024923788918701305PatilS. A.UnkiS. N.BadamiP. S.In vitro antibacterial, antifungal, and DNA cleavage studies of coumarin Schiff bases and their metal complexes: synthesis and spectral characterization201221124017402710.1007/s00044-011-9932-62-s2.0-84055211691RefatM. S.El-DeenI. M.AnwerZ. M.El-GholS.Bivalent transition metal complexes of coumarin-3-yl thiosemicarbazone derivatives: spectroscopic, antibacterial activity and thermogravimetric studies20099201-314916210.1016/J.MOLSTRUC.2008.10.0592-s2.0-59749103154BhavsarD.TrivediJ.ParekhS.SavantM.ThakrarS.BavishiA.RadadiyaA.ValaH.LunagariyaJ.ParmarM.PareshL.LoddoR.ShahA.Synthesis and in vitro anti-HIV activity of N-1,3-benzo[d]thiazol-2-yl-2-(2-oxo-2H-chromen-4-yl)acetamide derivatives using MTT method201121113443344610.1016/J.BMCL.2011.03.1052-s2.0-7995609113421515046KostovaI.Coumarins as inhibitors of HIV reverse transcriptase20064334736310.2174/1570162067777093932-s2.0-3374594118716842086ChenY.LiuH.-R.LiuH.-S.ChengM.XiaP.QianK.WuP.-C.LaiC.-Y.XiaY.YangZ.-Y.Morris-NatschkeS. L.LeeK.-H.Antitumor agents 292. Design, synthesis and pharmacological study of S- and O-substituted 7-mercapto- or hydroxy-coumarins and chromones as potent cytotoxic agents201249748510.1016/J.EJMECH.2011.12.0252-s2.0-8486278293822265685ThakurA.SinglaR.JaitakV.Coumarins as anticancer agents: a review on synthetic strategies, mechanism of action and SAR studies201510147649510.1016/J.EJMECH.2015.07.0102-s2.0-8493721751226188907KontogiorgisC.Hadjipavlou-LitinaD.Biological evaluation of several coumarin derivatives designed as possible anti-inflammatory/antioxidant agents2003181636910.1080/14756360310000692912-s2.0-0037320194KontogiorgisC. A.Hadjipavlou-LitinaD. J.Synthesis and antiinflammatory activity of coumarin derivatives200548206400640810.1021/JM05801492-s2.0-2644446124616190766AvdovićE. H.MilenkovićD.Dimitrić MarkovićJ. M.ĐorovićJ.VukovićN.VukićM. D.JevtićV. V.TrifunovićS. R.PotočňákI.MarkovićZ.Synthesis, spectroscopic characterization (FT-IR, FT-Raman, and NMR), quantum chemical studies and molecular docking of 3-(1-(phenylamino)ethylidene)-chroman-2,4-dione2018195314010.1016/J.SAA.2018.01.0232-s2.0-85041455013AvdovićE. H.MilenkovićD.Dimitrić-MarkovićJ. M.VukovićN.TrifunovićS. R.MarkovićZ.Structural, spectral and NBO analysis of 3-(1-(3-hydroxypropylamino)ethylidene)chroman-2,4-dione20171147697510.1016/J.MOLSTRUC.2017.06.0942-s2.0-85021114149MilenkovićD.AvdovićE. H.DimićD.BajinZ.RistićB.VukovićN.TrifunovićS. R.MarkovićZ. S.Reactivity of the coumarine derivative towards cartilage proteins: combined NBO, QTAIM, and molecular docking study2018149115916610.1007/s00706-017-2051-42-s2.0-85035102673HornykiewiczO.Dopamine miracle: from brain homogenate to dopamine replacement200217350150810.1002/mds.101152-s2.0-003658989712112197BermanR. M.NarasimhanM.MillerH. L.AnandA.CappielloA.OrenD. A.HeningerG. R.CharneyD. S.Transient depressive relapse induced by catecholamine depletion: potential phenotypic vulnerability marker?199956539540310.1001/archpsyc.56.5.3952-s2.0-003293185610232292KravitzH. M.SabelliH. C.FawcettJ.Dietary supplements of phenylalanine and other amino acid precursors of brain neuroamines in the treatment of depressive disorders1984841 Supplement1191236490415PaulingL.Orthomolecular psychiatry. Varying the concentrations of substances normally present in the human body may control mental disease19681603825265271564125310.1126/science.160.3825.2652-s2.0-0014422809LyrasL.CairnsN. J.JennerA.JennerP.HalliwellB.An assessment of oxidative damage to proteins, lipids, and DNA in brain from patients with Alzheimer’s disease1997685206120699109533SayreL. M.SmithM. A.PerryG.Chemistry and biochemistry of oxidative stress in neurodegenerative disease20018772173810.2174/09298670133729222-s2.0-003500343911375746ToshniwalP. K.ZarlingE. J.Evidence for increased lipid peroxidation in multiple sclerosis1992172205207153883410.1007/BF009668012-s2.0-0026500136EmeritJ.EdeasM.BricaireF.Neurodegenerative diseases and oxidative stress2004581394610.1016/j.biopha.2003.11.0042-s2.0-164248350914739060LemkeT. L.WilliamsD. A.FoyeW. O.ZitoS. W.2013Wolters Kluwer Health/Lippincott Williams & WilkinsSmythiesJ.Redox aspects of signaling by catecholamines and their metabolites20002357558310.1089/152308600501923322-s2.0-0033744516Gilgun-SherkiY.MelamedE.OffenD.Oxidative stress induced-neurodegenerative diseases: the need for antioxidants that penetrate the blood brain barrier200140895997510.1016/S0028-3908(01)00019-32-s2.0-003498636111406187AndersenJ. K.Oxidative stress in neurodegeneration: cause or consequence?2004107S18S2510.1038/nrn14342-s2.0-3142514196ShimizuT.NakanishiY.NakaharaM.WadaN.Moro-OkaY.HiranoT.KonishiT.MatsugoS.Structure effect on antioxidant activity of catecholamines toward singlet oxygen and other reactive oxygen species in vitro201047318119010.3164/jcbn.09-1122-s2.0-7804945325921103026SârbuC.CasoniD.Comprehensive evaluation of biogenic amines and related drugs’ antiradical activity using reactive 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical201311567968810.2478/s11532-013-0210-y2-s2.0-84875709539KładnaA.BerczyńskiP.KrukI.MichalskaT.Aboul-EneinH. Y.Superoxide anion radical scavenging property of catecholamines201328445045510.1002/bio.24752-s2.0-8488137705023319391DimićD.MilenkovićD.Dimitrić MarkovićJ.MarkovićZ.Antiradical activity of catecholamines and metabolites of dopamine: theoretical and experimental study20171920129701298010.1039/C7CP01716B2-s2.0-85024497737DimićD. S.MilenkovićD. A.MarkovićJ. M. D.MarkovićZ. S.Thermodynamic and kinetic analysis of the reaction between biological catecholamines and chlorinated methylperoxy radicals201811691166117810.1080/00268976.2017.14149672-s2.0-85038611550YamafujiK.MurakamiH.ShinozukaM.Antitumour activity of dopa, dopamine, noradrenalin or adrenalin and their reaction with nucleic acids197073319520310.1007/BF005243982-s2.0-0014711036WickM. M.Levodopa and dopamine analogs: melanin precursors as antitumor agents in experimental human and murine leukemia1979636991997466656DasguptaP. S.LahiriT.Antitumor effect of i.p. dopamine in mice bearing Ehrlich ascites carcinoma1987113436336810.1007/BF003977202-s2.0-0023185580WickM. M.Levodopa and dopamine analogs: dihydroxy and trihydroxybenzylamines as novel quinol antitumor agents in experimental leukemia in vivo1981659-108618677273019WickM. M.Levodopa and dopamine analogs as DNA polymerase inhibitors and antitumor agents in human melanoma1980405141414186768447DriscollJ. S.Catecholamine analogs as potential antitumor agents197968121519152110.1002/jps.26006812142-s2.0-0018641928StokesA. H.HastingsT. G.VranaK. E.Cytotoxic and genotoxic potential of dopamine199955665966510.1002/(SICI)1097-4547(19990315)55:6<659::AID-JNR1>3.0.CO;2-CAbdel LatifN. A.BatranR. Z.KhedrM. A.AbdallaM. M.3-Substituted-4-hydroxycoumarin as a new scaffold with potent CDK inhibition and promising anticancer effect: synthesis, molecular modeling and QSAR studies20166711612910.1016/J.BIOORG.2016.06.0052-s2.0-8497663491227372186MorganD. O.Cyclin-dependent kinases: engines, clocks, and microprocessors199713126129110.1146/annurev.cellbio.13.1.2612-s2.0-0031466305MeijerL.LeclercS.LeostM.Properties and potential applications of chemical inhibitors of cyclin-dependent kinases1999822-327928410.1016/S0163-7258(98)00057-62-s2.0-003303675810454205MorsyS. A.FarahatA. A.NasrM. N. A.TantawyA. S.Synthesis, molecular modeling and anticancer activity of new coumarin containing compounds201725687388310.1016/J.JSPS.2017.02.0032-s2.0-8501305947428951673SukdolakS.SolujićS.ManojlovićN.VukovićN.KrstićL. J.Hantzsch reaction of 3-(2-bromoacetyl)-4-hydroxy-chromen-2-one. Synthesis of 3-(thiazol-4-yl)-4-hydroxy coumarines200441459359610.1002/jhet.55704104182-s2.0-4344634288MorrisG. M.HueyR.LindstromW.SannerM. F.BelewR. K.GoodsellD. S.OlsonA. J.AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility200930162785279110.1002/jcc.212562-s2.0-70349932423DunningT. H.JrGaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen19899021007102310.1063/1.4561532-s2.0-33746614482BeckeA. D.JohnsonE. R.A density-functional model of the dispersion interaction200512315, article 15410110.1063/1.20652672-s2.0-2834443962816252936BeckeA. D.Density-functional thermochemistry. III. The role of exact exchange19939875648565210.1063/1.4649132-s2.0-0000189651BeckeA. D.Density-functional exchange-energy approximation with correct asymptotic behavior19883863098310010.1103/PhysRevA.38.30982-s2.0-4243553426GrimmeS.EhrlichS.GoerigkL.Effect of the damping function in dispersion corrected density functional theory20113271456146510.1002/jcc.217592-s2.0-7995294355921370243FrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.ScalmaniG.BaroneV.MennucciB.PeterssonG. A.NakatsujiH.CaricatoM.LiX.HratchianH. P.IzmaylovA. F.BloinoJ.ZhengG.SonnenbergJ. L.HadaM.EharaM.ToyotaK.FukudaR.HasegawaJ.IshidaM.NakajimaT.HondaY.KitaoO.NakaiH.VrevenT.MontgomeryJ.AJ.PeraltaJ. E.OgliaroF.BearparkM.HeydJ. J.BrothersE.KudinK. N.StaroverovV. N.KobayashiR.NormandJ.RaghavachariK.RendellA.BurantJ. C.IyengarS. S.TomasiJ.CossiM.RegaN.MillamJ. M.KleneM.KnoxJ. E.CrossJ. B.BakkenV.AdamoC.JaramilloJ.GompertsR.StratmannR. E.YazyevO.AustinA. J.CammiR.PomelliC.OchterskiJ. W.MartinR. L.MorokumaK.ZakrzewskiV. G.VothG. A.SalvadorP.DannenbergJ. J.DapprichS.DanielsA. D.FarkasÖ.ForesmanJ. B.OrtizJ. V.CioslowskiJ.FoxD. J.2009Gaussian Inc., Wallingfordhttp://www.gaussian.comFuhrmannJ.RurainskiA.LenhofH. P.NeumannD.A new Lamarckian genetic algorithm for flexible ligand-receptor docking20103191911191810.1002/jcc.214782-s2.0-77953189302Rigaku Oxford Diffraction2017Yarnton, EnglandRigaku Oxford DiffractionSheldrickG. M.SHELXT – integrated space-group and crystal-structure determination20157113810.1107/S20532733140263702-s2.0-84920370689SheldrickG. M.Crystal structure refinement with SHELXL20157113810.1107/S20532296140242182-s2.0-8493739880925567568FarrugiaL. J.WinGX suite for small-molecule single-crystal crystallography199932483783810.1107/S00218898990060202-s2.0-0141452964SpekA. L.Structure validation in chemical crystallography200965214815510.1107/S090744490804362X2-s2.0-5884916185719171970BrandenburgK.2012Bonn, GermanyCrystal Impact GbRBernsteinJ.DavisR. E.ShimoniL.ChangN. L.Patterns in hydrogen bonding: functionality and graph set analysis in crystals199534151555157310.1002/anie.1995155512-s2.0-33748502022GilliG.BellucciF.FerrettiV.BertolasiV.Evidence for resonance-assisted hydrogen bonding from crystal-structure correlations on the enol form of the .beta.-diketone fragment198911131023102810.1021/ja00185a0352-s2.0-10844256940MałeckaM.GrabowskiS. J.BudziszE.Crystal and molecular structures of 3-[1-(2-hydroxyethylamino)-ethylidene]-chroman-2,4-dione and 2-methoxy-3-[1-(benzylamino)-ethylidene]-2,3-dihydro-2,4-dioxo-2λ5-benzo[e][1,2]oxaphosphinane and DFT study of intramolecular H-bonds of related compounds20042971-323524410.1016/J.CHEMPHYS.2003.10.0292-s2.0-0346500436AvdovićE. H.StojkovićD. L. J.JevtićV. V.KosićM.RistićB.Harhaji-TrajkovićL.VukićM.VukovićN.MarkovićZ. S.PotočňákI.TrifunovićS. R.Synthesis, characterization and cytotoxicity of a new palladium(II) complex with a coumarin-derived ligand 3-(1-(3-hydroxypropylamino)ethylidene)chroman-2,4-dione. Crystal structure of the 3-(1-(3-hydroxypropylamino)ethylidene)-chroman-2,4-dione201746618819610.1016/J.ICA.2017.06.0152-s2.0-85020718781StojkovićD. L.JevtićV. V.VukovićN.VukićM.PotočňákI.ZelenI. R.ZarićM. M.MišićM. M.BaskićD.KaluđerovićG. N.TrifunovićS. R.Crystal and molecular structure of a new palladium(II) complex with a coumarin-valine derivate201758355055710.1134/S00224766170301792-s2.0-85021192899MałeckaM.BudziszE.3-[(1-Benzylamino)ethylidene]-2H-chromene-2,4(3H)-dione20066211o5058o506010.1107/S16005368060415112-s2.0-33750884626BrahmiaA.Ben AyedT.Ben HassenR.3-[1-(2-Hydroxyanilino)ethylidene]-3H-chromen-2,4-dione2013698o1296o129610.1107/S160053681301934X2-s2.0-8488147006424109374GroomC. R.BrunoI. J.LightfootM. P.WardS. C.The Cambridge Structural Database201672217117910.1107/S20525206160039542-s2.0-8496354827727048719SinghR. K.LangeT. S.KimK. K.BrardL.A coumarin derivative (RKS262) inhibits cell-cycle progression, causes pro-apoptotic signaling and cytotoxicity in ovarian cancer cells2011291637210.1007/s10637-009-9335-42-s2.0-7915147032019865799KhaghanzadehN.MojtahediZ.RamezaniM.ErfaniN.GhaderiA.Umbelliprenin is cytotoxic against QU-DB large cell lung cancer cell line but anti-proliferative against A549 adenocarcinoma cells20122016910.1186/2008-2231-20-692-s2.0-84870512215MusaM. A.JosephM. Y.LatinwoL. M.BadisaV.CooperwoodJ. S.In vitro evaluation of 3-arylcoumarin derivatives in A549 cell line201535265365925667442NasrT.BondockS.YounsM.Anticancer activity of new coumarin substituted hydrazide–hydrazone derivatives20147653954810.1016/J.EJMECH.2014.02.0262-s2.0-8489754414224607878AminK. M.EissaA. A. M.Abou-SeriS. M.AwadallahF. M.HassanG. S.Synthesis and biological evaluation of novel coumarin–pyrazoline hybrids endowed with phenylsulfonyl moiety as antitumor agents20136018719810.1016/J.EJMECH.2012.12.0042-s2.0-8487171311023291120NinkovićD. B.Vojislavljević-VasilevD. Z.MedakovićV. B.HallM. B.BrothersE. N.ZarićS. D.Aliphatic–aromatic stacking interactions in cyclohexane–benzene are stronger than aromatic–aromatic interaction in the benzene dimer20161837257912579510.1039/C6CP03734H2-s2.0-8498856888127711596NinkovićD. B.AndrićJ. M.MalkovS. N.ZarićS. D.What are the preferred horizontal displacements of aromatic–aromatic interactions in proteins? Comparison with the calculated benzene–benzene potential energy surface20141623111731117710.1039/C3CP54474E2-s2.0-8490126747224805772KamalA. A. M.PetreraL.EberhardJ.HartmannR. W.Structure–functionality relationship and pharmacological profiles of Pseudomonas aeruginosa alkylquinolone quorum sensing modulators201715214620463010.1039/C7OB00263G2-s2.0-85021776931