Isolation of dihydrobenzofuran derivatives from ethnomedicinal species Polygonum barbatum as anticancer compounds

Background Ethnomedicinally, the family Polygonaceae is famous for the management of cancer. Various species of this family have been reported with anticancer potentials. This study was designed to isolate anticancer compounds from ethnomedicinally important species Polygonum barbatum. Methods The column chromatography was used for the isolation of compounds from the solvent fraction of P. barbatum. The characterization of isolated compounds was performed by various spectroscopic techniques like UV, IR, mass spectrometry and 1D-2D NMR spectroscopy. Keeping in view the ethnomedicinal importance of the family, genus and species of P. barbatum, the isolated compounds (1–3) were screened for anticancer potentials against oral cancer (CAL-27) and lungs cancer (NCI H460) cell lines using MTT assay. Active compound was further investigated for apoptosis by using morphological changes and flow cytometry analysis. In vivo anti-angiogenic study of the isolated compounds was also carried using chorioallantoic membrane assay. Docking studies were carried out to explore the mechanism of anticancer activity. Results Three dihydrobenzofuran derivatives (1–3) have been isolated from the ethyl acetate fraction of P. barbatum. The structures of isolated compounds were elucidated as methyl (2S,3S)-2-(3,4-dimethoxyphenyl)-4-((E)-3-ethoxy-3-oxoprop-1-en-1-yl)-7-methoxy-2,3-dihydrobenzo-furan-3-carboxylate (1), (E)-3-((2S,3S)-2-(3,4-dimethoxyphenyl)-7-methoxy-3-(methoxy carbonyl)-2,3-dihydrobenzofuran-4-yl)acrylic acid (2) and (2S,3S)-4-((E)-2-carboxyvinyl)-2-(3,4-dimethoxyphenyl)-7-hydroxy-2,3-dihydrobenzofuran-3-carboxylic acid (3). The compound 1 was found to be more potent with IC50 of 48.52 ± 0.95 and 53.24 ± 1.49 against oral cancer cells as compared to standard drug (IC50 = 97.76 ± 3.44 μM). Both compound also inhibited lung cancer cells but at higher concentrations. Morphological and flow cytometry analysis further confirms that compound 1 induces apoptosis after 24 to 48 h treatment. In antiangiogenesis assay, compounds 1, 2 and 3 exhibited IC50 values of 8.2 ± 1.1, 13.4 ± 1.1 and 57.7 ± 0.3 μM respectively. The docking studies revealed that the compounds under study have the potential to target the DNA and thymidylate synthase (TS). Conclusion Based on its overwhelming potency against the tested cell lines and in angiogenesis assay, compound 1 can be further evaluated mechanistically and can be developed as anticancer drug candidate. Electronic supplementary material The online version of this article (10.1186/s40659-018-0209-0) contains supplementary material, which is available to authorized users.


Introduction
Polygonum barbatum is a perennial herb, mostly found in shady and marshy places near river banks, ponds and other aquatic areas in different countries of South-East Asia [25,26]. The species belongs to the family Polygonaceae represented by 19 genera and 103 species in Pakistan out of total 55 genera and approximately 1200 species worldwide [21]. Traditionally, various species of Polygonaceae family are used for the management of cardiovascular diseases, hypertension and cancer [14,29]. Based on the ethnomedicinal background, some species have been verified scientifically for the management of cancer [1,2,4,17]. The leaf extract of P. barbatum has been used in Chinese traditional medicine for the treatment of ulcer, and root extract as astringent and as carminative agent [10,20,30].
Anti-inflammatory, antinociceptive and diuretic activities of P. barbatum has been previously reported [10,25]. Similarly, P. barbatum also possess potential analgesic, brine shrimps, cholinergic, as well as spasmolytic potentials [10,11,24]. In indigenous system of medicine, it has been used as potential anticancer agent against oral cancer [18]. P. barbatum has also been reported to have antimicrobial activity, suggesting it to be a source of potential antibiotic agents [18].
Additionally, the therapeutic and biological potentials of some species of the genus Polygonum have also been previously reported [5-8, 19, 21, 22]. Phytochemical investigations of various species of this genus showed the presence of secondary metabolites, particularly anthraquinones, flavonoids, as well as phenylpropanoid [9,16]. Furthermore, secondary metabolites such as acetophenone, viscozulenic acid and sitosterone have been previously isolated from P. barbatum [26]. Similarly, Farooq et al. reported sesquiterpenes from P. barbatum possessing potential anti-proliferative activity and inhibitory effect on cancer cell migration [13].
Cancer, a second leading cause of death, resulted in 8.8 million deaths in 2015 globally [32]. In Pakistan, the three deadliestof them are breast cancer, oral cancer and lungs cancer with death rate of 26.76, 8.71 and 5.53 per 100,000 population, respectively [33]. Despite of the current advances in therapeutic approaches, there is an immense need of new potential compounds that can be used as safe anticancer drugs. Keeping in view the ethnomedicinal importance of Polygonaceae and potential therapeutic efficiency of P. barbatum, we were able to isolate three dihydrobenzofuran derivatives from its ethyl acetate fraction. To study the anticancer effect, we targeted oral squamous cell carcinoma and large cell lung carcinoma. Furthermore, we also assessed the potential of these isolated dihydrobenzofuran derivatives on normal vascular development in chick embryos using a multifaceted and low-cost chicken chorioallantoic membrane (CAM) assay.

Methods
The Polygonum barbatum (whole plant) was collected from different areas of Khyber Pakhtunkhwa in October 2015 and specimen having voucher # 66130 has been deposited in herbarium at Department of Botany, postgraduate college Abbottabad, Pakistan.
The shade dried plant (5.4 kg) was ground into fine powder, extracted with methanol and filtered twice. The crude extract (245 g) was obtained from filtrate by using vacuum rotary evaporator and partitioned into n-hexane (85 g), ethyl acetate (48 g) and n-butanol (94 g) fraction.
Further investigation and purification of compounds was done through column chromatography of ethyl acetate fraction selected based on TLC.

Instrumentation
The NMR spectra were recorded by using Bruker AMX-400 spectrometer ( 1 H-NMR at 400 MHz and 13 C-NMR at 100 MHz). The values of chemical shift and coupling constant were recorded as ppm and Hertz (Hz) respectively. The Varian MAT-312 spectrometer was used to record HR-EI-MS and EI-MS spectra. The E. Merck 230-400 mesh silica gel was used to pack columns while pre-coated silica gel plates were employed for TLC analysis. In addition, Ceric sulphate solution in 10% sulphuric acid (H 2 SO 4 ) has been used for detection of UV active compounds.

Anticancer activity
The isolated compounds (1-3) were investigated for anticancer activity against human oral squamous cell carcinoma cell line (CAL 27) and human large cell lung carcinoma cell line (NCI H460) as described recently [13]. The cell lines were purchased from ATCC and maintained in Dulbecco's Modified Eagle's Medium (DMEM) supplemented with 10% Fetal Bovine Serum (FBS), 100 µg/mL streptomycin and 100 units/mL penicillin at 37 °C in a humidified 5% CO 2 incubator. The CAL 27 and NCI H460 cells were seeded in 96 well plates with the density of 1.5 × 10 4 and 1 × 10 4 cells per well, respectively and incubated for 24 h. The cells were then treated with different concentrations of compound 1, 2 and 3 for 48 h. 1% DMSO was used as vehicle control and empty wells were used as negative control. The standard anticancer drugs namely 5-Fluorouracil and Cisplatin were used against CAL 27 and NCI H460, respectively. Quantification of metabolically active cells were determined using 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide (MTT) dye which reduces to insoluble Formazan crystals by the action of mitochondrial succinate dehydrogenase. The soluble MTT dye was aspirated off and insoluble crystals were dissolved using DMSO. The optical density of each well was determined by Spec-traMax microplate reader at 570 nm. IC 50 values were calculated using software EZ-Fit Enzyme Kinetics by Perrella Scientific. For morphological changes cells treated at IC 50 and IC 70 were observed under phase contrast microscopy (Nikon, Japan) at 24 and 48 h respectively.

YO-PRO-1 apoptotic assay
Oral cancer (CAL-27) cells (5 × 10 5 ) were seeded in a 6-well plate and incubated for 24 h in presence of 5% CO 2 . After incubation compound 1 was added at IC 50 and IC 70 concentration and further incubated for 48 h. Following incubation cells were harvested by using trypsin (0.05%) and pallet out after washing three times at 400×g for 5 min to remove media and trypsin. Cell pallets were suspended in 1 mL PBS and stained with Vybrant Apoptosis Assay kit # 4 (Invitrogen, USA) containing YO-PRO-1 and Propidium Iodide (PI) dye according to manufacturer protocol. Briefly, cells were incubated with 1 μL YO-PRO-1 stock solution (Component A) and 1 μL PI stock solution (Component B) for 30 min on ice. Later, cells were analyzed on FACSCalibur (Becton-Dickinson, USA). YO-PRO-1 and PI were excited at 488 nm, and fluorescence was measured at 530 and 620 nm, respectively. A total of 10,000 events were acquired from each sample. The percentages of live, apoptotic and dead cells were determined using CellquestPro software.

Chorioallantoic membrane (CAM) assay
Chorioallantoic membrane (CAM) assay was performed according to our previously reported procedure [15]. The experiments were carried out as per the approval of the ethical committee, Department of Pharmacy, University of Malakand, Pakistan according to the animals Bye-Laws 2008 (Scientific Procedure Issue-1). The fertilized domestic chicken eggs were purchased from poultry form Chakdara, Pakistan. The fertilized eggs were incubated for 5-6 days at 37 °C with slowly shaking for at least two to three times a day. After the incubation period, the seven-day old eggs were examined under the flash light for detection of the embryo head. After that a small hole was bored at the narrow end of the eggs and 0.5-1 mL of albumin was sucked with the help of eighteen-gauge hypodermic needle so that yolk sacs fell down from the shell membrane. The shell around the embryo air sac was removed through forceps and the shell membrane at the base of air sac was peel away. On 8th day, a thermanox cover slip was carefully placed on the surface of CAM loaded with different samples for each concentration (10 µL) and were placed in incubator. After 24 to 48 h, the numbers of vessels on the surface of CAM were observed. The numbers of blood vessels in each egg was counted by figuring out the actual number of blood vessels in 1 cm 2 area of the CAM adjacent to the sample application site. Vessels radially converging in the direction of the center were counted under a microscope. At least twenty eggs were used for each sample dose. The % of increase and inhibition were calculated using formula; Median inhibitory concentrations (IC 50 ) for anti-angiogenic assay were calculated via regression analysis of dose response curve by adding all the concentration and percent activity.

Docking studies
Docking studies were carried out using Molecular Operating Environment [27] software. Crystal structures of DNA and Human Thymidylate Synthase (TS) was retrieved from Protein Data Bank with code numbers 3LPV and 1JU6 respectively. Validation of docking procedure was carried out by using re-dock method. Ligand and DNA/enzyme preparation, active site determination % inhibition = No. of blood vessels (Norml Saline) − No. of blood vessles (Test Sample) No.of blood vessels (Norml Saline) ×100 and docking procedure was carried out by using our previously reported procedures [3,15]. The lowest energy minimized pose was used for further analysis. Ligandinteraction module of MOE was used to calculate the 2D ligand-enzyme interactions. The view of the docking results and analysis of their surface with graphical representations were done using MOE and discovery studio visualizer [12].

Isolation and Characterization of benzofuran derivatives
The ethyl acetate fraction of P. barbatum upon column chromatography gave three benzofuran derivatives (1-3).
Using various spectroscopic analyses and published literature, the structures of isolated compounds were determined.   Table S1 in Additional file 1).
The 13 C-NMR along with DEPT spectrum showed signals for 24 carbon atoms comprising of one methyl, one methylene, nine methine, seven quaternary and four methoxy groups. In addition, two carbonyl carbon atoms were observed at δ C 170.1 (C-9) and δ C 174.6 (C-9′) while the olefinic carbons at position H-7 and H-8 resonated at δ C 143.7 and 120.7 respectively (as shown in Fig. 1). The only downfield methylene appeared at δ C 66.2, while the methyl group resonated at δ C 19.4 ppm. The two methine carbons of dihydrobenzofuran centered at δ C 82.4 for C-7′ and δ C 58.1 for C-8′. The chemical shift values of aromatic carbon and methoxy groups have been shown in Additional file 1: Table S2.
The HMBC spectrum revealed strong interaction of protons at H-5 and H-6 with C-4 and C-1 while protons of methoxy group substituent at this position (C-4) also showed strong correlation with C-4. The olefinic proton at position H-7 showed strong HMBC correlation with C-1 and C-8 while other olefinic proton at H-8 showed HMBC correlation with C-1, C-7 and C-9 (carbonyl carbon of ethyl ester moiety). The placement of methyl ester as substituent on dihydrobenzofuran ring was confirmed by HMBC correlation of methoxy protons were connected with carbonyl carbon at C-9′ as evident from Fig. 2. This placement was further supported by H-7′ and H-8′ protons which showed strong correlations with C-9′ carbon. The position of two methoxy group at position 3′ and 4′ was also confirmed by their HMBC correlation as depicted in Fig. 2. The aromatic protons at position 2′ and 6′ showed strong correlation with C-1′ and C-7′ confirming the placement of this aromatic ring at dihydrobenzofuran ring.
The relative configuration of aliphatic protons at position H-7′ and H-8′ was suggested trans to each other by Nuclear Overhauser Effect (NOE) as well as the coupling constant value (i.e., 5.2 Hz) suggested them to be at Trans position to each other. All the spectral data of 1 H-NMR, 13 C-NMR, HR-EI-MS, HMBC, and COSY along with comparison of previous literature   [34]. The trans orientation of H-7′ and H-8′ was further supported by NOE experiment. When H-7′ was irradiated, H-8′ did not show any increase in intensity of H-8′ confirmed the trans orientation. In addition, singlets for methoxy groups were appeared at δ H 3.60 (3H, s, 4-OMe), 3.62 (3H, s, 3′-OMe) and 3.65 (3H, s, 4′-OMe), while the methoxylgroup of ester moiety appeared at δ H 3.53 as singlet in 1 H-NMR spectrum.
The 13 C-NMR spectra of compound 2 were quite identical to compound 1 that showed the presence of 22 carbon atoms including one methoxyl group of ester moiety, nine methine, seven quaternary and three methoxy groupson aromatic rings along with two carbonyl carbon atoms. The chemical shift value of olefinic carbon atoms were δ C 142.4 (C-7) and 119.6 (C-8) while two carbonyl carbon resonated at δ C 172.1 (C-9) and 174.1 (C-9′). The carbon atoms of five membered ring of dihydrobenzofuran resonated at δ C 84.1 (C-7′) and 59.2 (C-8′) while chemical shift values of aromatic carbon atoms has been shown in Additional file 1: Table S2.

Compound 3
Compound 3 was isolated as yellow gummy solid, that showed molecular ion peak of 386.997 in HR-EI-MS suggesting a molecular formula as C 20  Similarly, 13 C-NMR data was also in close agreement with those of compound (1-2). The 13 C-NMR spectra suggested the presence of two carboxylic acid groups in compound 3, which appeared at δ C 173.1 (C-9) and δ C 175.6 (C-9′). It was also evident that the fused benzene ring with furan ring was substituted with hydroxyl group at position 4 with chemical shift value of δ C 146.6.
The 1 H-NMR and 13 C-NMR spectra of compound 3 were quite similar to compounds 1 and 2. The differences observed were presence of carboxylic acid instead of ester substituent at position-8 while at position-4 a hydroxyl (-OH) group was observed as substituent rather than methoxy group. The difference was retrieved from UV, IR and NMR data which was further supported by mass spectrum. The chemical shift values of both 1 H-NMR and 13 C-NMR spectra have been given in Additional file 1: Tables S1 and S2. Similarly, HMBC correlation was also found similar to compound 1 and 2 as shown in Fig. 2. All spectra data and comparison of literature suggested the structure of compound 3 as (2S,3S)-4-((E)-2-carboxyvinyl)-2-(3,4-dimethoxyphenyl)-7-hydroxy-2,3-dihydrobenzofuran-3-carboxylic acid.

In vitro anticancer activity
As widely spread natural compounds, benzofuran derivatives have been extensively studied for their biological activities [23]. Previously, a series of synthetic dihydrobenzofuran derivatives showed in vitro anticancer potential [23]. Most of these compounds have been reported to have cytotoxic potential against many human cancer cell lines including leukemia (HL-60-TB and K-562), breast cancer (MDA-MB-435. MDA-N, BT-549 and MCF-7), human cervical cancer (HeLa) and human lung carcinoma (A549) cell lines [23,28]. Based on the previous reports, we screened the anticancer potential of the three isolated dihydrobenzofuran derivative from P. barbatum against oral (CAL-27) and lung cancer (NCI-H460) cells. Lung and oral cancer was selected by  Table 1). Whereas, these two compounds also showed anticancer activities against lungs cancer cells (NCI-H460) at higher doses as compared to control drug (Cisplatin; Table 1). Compound 1 was selected for further apoptosis analysis by observing morphological changes and FACS analysis considering its IC 50 which was almost half the IC 50 of standard drugs. After treatment with compound 1 at IC 50 concentration cell morphology started to change after 24 h and with the passage of time cell detach from the surface with loss in cell shape suggesting apoptosis as compared to control. On the other hand, at higher concentration, i.e. IC 70 cells were undergone apoptosis after 24 h (Fig. 3). Flow cytometry analysis (FACS) also correspond with the morphological data that after 48 h of treatment significant apoptosis (either early or late) and necrosis was observed (Fig. 4).

Antiangiogenic studies
Chicken chorioallantoic membrane (CAM) is a low cost in vivo tumor biology model and is used widely for the tumor angiogenesis studies. Tumor growth depends on blood supply which in turn depends on the number of vessels. Hence, the inhibition of neovascularization is one of the targets. We recently reported dihydropyrimidine-based compounds that have potential to inhibit of new blood vessels formation [15]. Moreover, we have also employed CAM assay for the antitumor and antiangiogenic potentials of crude extracts of some medicinal plants. In current research, we also studied the potential of our isolated dihydrobenzofuran derivatives on normal vascular development in chick embryos using CAM model. Normal saline was used as negative control. Dexamethasone and methotrexate was used as positive control [15].
The antiangiogenic activities of the isolated compounds 1, 2 and 3 were based on the number of blood vessels formed in the chorioallantoic membrane (CAM) of chick embryos. Results of anti-angiogenic assay of isolated compounds are shown in Table 2. Additional file 1:   Table S3 showed the average number of blood vessels of each concentration in the CAM of chick embryo using different five concentration of each compound. Chick embryos treated with distilled water yielded the highest mean of 30.4 blood vessels formed. This indicates that water, being the negative control, is the least antiangiogenic among the five treatments. The compound 3 (IC 50 = 57.7 ± 0.3 μM, Table 2), yielded moderate number of blood vessels with average mean of about 10.8. The compound 2 (IC 50 = 13.4 ± 1.1 μM) yielded an average of 8.6 blood vessels formed. While, the compound 1 (IC 50 = 8.2 ± 1.1 μM) yielded a mean of only 5.0 blood vessels which were most active against antiangiogenic activity as given in Fig. 5. The IC 50 value was calculated via regression analysis of dose response curve by adding all the concentrations and percent activities ( Figure S1 in Additional file 1).

Docking studies
In cancer research, cancer cell lines have emerged as an alternative tool for the study of biological mechanisms. However, to hypothesize a biomolecular target that help us to explore the mechanism is a difficult task. The standard drugs used in current research are cisplatin and 5-fluorouracil. Cisplatin is a chemotherapeutic agent that binds with DNA to form a bifunctional 1,2-intrastand