The genus Euphorbia is the largest among the plant family Euphorbiaceae, comprising about 2000 known species consisting of annuals to trees. They all contain latex and have unique flower structures. About 70 species have been reported in Iran that 17 of them are endemic (1) . E. erythradenia Bioss. is one of Iranian endemic spurges (2) which few phytochemical studies have been reported for this plant. Some species of this genus are used in folk medicines to cure skin diseases, gonorrhea, migraines, intestinal parasites, and warts (1) which most of them have been proved by modern phytochemical and pharmacological studies. In Iran, they are also used as a purgative.
The presence of several kinds of secondary metabolites including triterpenes, steroids, flavonoids, phenolics and aromatic compounds in Euphorbia species have been reported up until now (1, 3) , but what highlights the features of Euphorbia's secondary metabolites, is its specific terpenoids profile. Euphorbia's terpenoids, including diterpenes and triterpenes, have shown interesting various biological activities which are attracting medicinal researchers increasingly as a source of new potential drugs as well as lead compounds.
In this paper, as a part of our efforts to find new sources of bioactive terpenoids, one of Iranian Euphorbia species, E. erythradenia arial parts were phytochemically investigated and five known triterpene compounds extracted and their chemical structure were elucidated ( Figure 1 ). Considering the huge demand for finding new drugs for viral infections, specially acquired human immunodeficiency virus (HIV), anti HIV activity of the extracted triterpenoids plus the four previously reported ingenoid diterpenes from E. erythradenia (Figure 2 ), were evaluated, using single cycle replicable HIV virions. In order to find the important structural features of anti HIV-1 active compounds, ligand-protein binding studies were carried out using AutoDock software.
Column chromatography (CC): silica gel 63-200 μm; TLC: silica gel GF 254 plates (20×20 cm, 0.5 mm, Merck); detection by spraying with cerium sulphate in 10% aq. H 2 SO 4 and heating. Preparative HPLC: Agilent 1100 Series with a normal phase column (250×20 mm i.d.) packed with 5µm silica (YMC Co., Ltd., Kyoto, Japan) that was equipped with 515 Waters pump and UV-Vis detector. All solvents used were of HPLC grade (Caledon Laboratory Chemicals Ltd. Georgetown, Ont. Canada). NMR: Bruker AV-400 ( 1 H) and AV-100 ( 13 C), δ in ppm rel. to Me4Si and J in Hz; HR-ESI-MS: Waters Q-TOF Micro YA019 mass spectrometer in m/z and EI-MS spectra: Varian MAT 112 or MAT 312 spectrometers. four previously reported ingenoid diterpenes from E. erythradenia (Figure 2 ), were evaluated, using single cycle replicable HIV virions. In order to find the important structural features of anti HIV-1 active compounds, ligand-protein binding studies were carried out using AutoDock software. Cycloart-25-ene-3β,24-diol(1b); Cycloart-23-ene-3β,25-diol (2); 24, ; Oleanolic acid (4). 
Aerial flowering parts of E. erythradenia f. (Euphorbiaceae) were collected in September 2010 from populations growing in Gharbalbiz in the neighborhood of Mehriz city, Yazd province (I.R.Iran). Plant material was identified by Ali Mirhosseini, plant taxonomist (Agricultural and natural Resources Research Center of Yazd Province) and a voucher specimen (nos. 1947) was deposited in the Herbarium Center of Yazd (HCY).
The air-dried powdered plant (2.8 kg) was macerated for four days with methanol (9 L×3), at room temperature. Filtration and under vacuum evaporation resulted in a green gum (433.2 g, 15.47% of dried powder weight), that was suspended in 1 liter methanol 80% and fractionated with hexane (1 lit × 3) which led to 51.3 g (1.83% of dried powder weight) hexane fraction. Methanolic suspension was further fractionated three times with 1 liter of chloroform that resulted in 37.7g chloroform fraction (1.35% of dried powder weight). Chloroform fraction was subjected to silica gel column, using hexane/chloroform (30→100) and chloroform/methanol (0→20), to render several fractions: F [1] [2] [3] [4] . Subsequently, the richest fraction in triterpenoids, F 2 (based on TLC visualization after cerium sulfate solvent spraying and heating), was chromatographed on silica gel C.C. (hexane/ethyl acetate, 0→60) to afford five subfractions; TF 2a -TF 2e . Subfractions were injected into preparative HPLC using hexane/ethyl acetate gradient system at detection wavelength of 270 nm. Fractions were inspected and combined using TLC and heating after acidic cerium sulfate solvent spraying which resulted into pure compounds: TF 2b1 (1a & 1b, 15 mg), TF 2c2 (2, 18 mg), TF 2d5 (3, 13mg) . Moreover F 3 was subjected to column of silica gel and washed with hexane/ethyl acetate (0→20). Then TF 3.101-103 subfractions were combined and recrystallized with cold chloroform which resulted to compound TF 3.101 (4, 26 mg). Table 1 ; 13 C-NMR data (100 MHz, CDCl 3 ) see Table 2 . Table 1 ; 13 C-NMR data (100 MHz, CDCl 3 ) see Table 2 . Isolated diterpenes from E. erythradenia. 13-acetyl-3-butanoyl-5-angeloyl-20-nicotinyl-1,2,6,7tetrahydroingenol (5); 7,13-diacetyl-5-angeloyl-20-nicotinyl-3-propionyl-1,2,6,7-tetrahydroingenol (6); 5,13-diacetyl-7-benzoyl-20-nicotinyl-3-propionyl-1,2,6,7-tetrahydroingenol (7); 13-acetyl-5-angeloyl-11nicotinyl-3-butanoyl-1,2,6,7-tetrahydroingenol (8).
Column chromatography (CC): silica gel 63-200 μm; TLC: silica gel GF254 plates (20×20 cm, 0.5 mm, Merck); detection by spraying with cerium sulphate in 10% aq. H2SO4 and heating.
Preparative HPLC: Agilent 1100 Series with a normal phase column (250×20 mm i.d.) packed with 5µm silica (YMC Co., Ltd., Kyoto, Japan) that was equipped with 515 Waters pump and UV-Vis detector. All solvents used were of HPLC grade (Caledon Laboratory Chemicals Ltd.
Georgetown, Ont. Canada). NMR: Bruker AV-400 ( 1 H) and AV-100 ( 13 C), δ in ppm rel. to Isolated diterpenes from E. erythradenia. 13-acetyl-3-butanoyl-5-angeloyl-20-nicotinyl-1,2,6,7-tetrahydroingenol (5); 7,13-diacetyl-5-angeloyl-20-nicotinyl-3-propionyl-1,2,6,7tetrahydroingenol (6); 5,13-diacetyl-7-benzoyl-20-nicotinyl-3propionyl-1,2,6,7-tetrahydroingenol (7); 13-acetyl-5-angeloyl-11-nicotinyl-3-butanoyl-1,2,6,7-tetrahydroingenol (8). Table 1 ; 13 C-NMR data (100 MHz, CDCl 3 ) see Table 2 .
Cycloart-3β,24,25-triol ( Table 1 ; 13 C-NMR data (100 MHz, CDCl 3 ) see Table 2 . Table 1 ; 13 C-NMR (100 MHz, CDCl 3 ) see Table  2 . HREI-MS m/z 457.3623 (Calcd. for C 30 H 48 O 3 + H + ).
Single Cycle Replicable HIV-1 virions (SCR) were previously constructed and their characterization ascertained the production of infective HIV-1 virions with the capability of one cycle of replication (4, 5). Briefly, these single cycle replicable HIV-1 (SCR HIV-1) virions were produced by deleting a 2-kb fragment within the Pol region of the HIV-1 genome from the pNL4-3 strain. Pseudotyped SCR HIV-1 virions were constructed by co-transfection of HEK 293T cells with pmzNL4-3 (containing deleted genome), psPAX2, and pMD2G plasmids obtained from Addgene (www.addgene.org).
The pmzNL4-3 plasmid encodes the HIV-1 fulllength RNA, with packaging ability containing the aforementioned deletion in the Pol region; the psPAX2 plasmid encodes HIV Gag and Gag-Pro-Pol polyproteins, besides all the HIV-1 accessory proteins; and the pMD2G plasmid encodes the vesicular stomatitis virus surface glycoprotein (VSVG), which is required for assembly and the budding process of virus. These pseudotyped virions have the ability of infecting a broad spectrum of cell, even without the CD4 receptor (including Hela cells).
Hela cells, which were used as target cells in this experiment, were seeded at a density of 6 × 10 4 cells per well in 96-well plates. Each well was infected with 600 ng of p24 single cycle replicable HIV-1 (SCR HIV-1) virions. After 24 h of virus adsorption, cells were washed 3 times with pre-warmed DMEM to remove free virus particles. Cells were then incubated for 48 hours in a total volume of 200 μL per well of fresh medium containing various concentrations of compounds 1-8. Nevirapine (a HIV-1/2 RT inhibitor) was used as positive control. After 48 h, the p24 antigen (Ag) assay was performed on the supernatants by using a quantitative p24 ELISA method (HIV p24 ELISA, BioMerieux, France), according to the manufacturer's protocol. The IC 50 of compounds 1-8 was calculated according to the method described by Cheng et al.(6) . The therapeutic index (TI) was evaluated as the ratio of CC 50 to IC 50 .
The cellular toxicity of compounds 1-8 in Hela cells was assessed using a cell proliferation XTT kit (Roche Diagnostics, Germany), as described previously (7). Briefly, cells were plated in duplicate in 96-well plates in the presence or absence of various concentrations of compounds 1-8. After incubation at 37°C with 5% CO2 for 3 days, 50 μL of prepared XTT mixture was added to each well. The cells were incubated for an additional 4 hours to allow the production of XTT formazan. Absorbance was measured using an ELISA plate reader (BioTek ELx800) at a test wavelength of 450 nm and a reference wavelength of 690 nm. Percent inhibition was calculated using the following formula: Inhibition (%) = [100 -(A t /A s )] × 100, where A s is the absorbance of the solvent and A t , of the test sample, respectively. The cytotoxic concentration that resulted in a reduction of the number of viable cells by 50% (CC 50 ) was calculated from dose response curves.
Docking studies using AutoDock software The high resolution crystal structure of HIV protease (PDB code: 1AJX) complexed with Aha001 was retrieved from PDB Protein Data Bank and its ligand was deleted from the active site. Compound 7 was constructed on HyperChem 8.0.3 version and minimized using AMBER force field molecular mechanics and PolakRiebiere algorithm with RMS gradient = 0.01 Kcal/mol.
The receptor was kept rigid, and ligand was allowed to be flexible. Polar hydrogens and Kollman united atom partial charges were added to the individual protein atoms. A docking grid box was built with 28, 28 and 28 points in 12.79, 23.30 and 5.85 directions in the catalytic site of protein and the number of generations and maximum number of energy evaluations was placed on 100 and 2,700,000, respectively. Docking results were clustered with a root mean square deviation (RMSD) of 0.5 Å and evaluated by Pymol software.
Compound 1a & 1b showed a pair of doublets in the up-field area (δ H 0.32, J = 4.0 Hz and 0.55, J = 4.0 Hz) that develops hypothesis of cycloartane backbone of them. It should be mentioned that a pair of doublets in the up-field area in 1 H-NMR spectrum is a wellknown characteristic of cyclopropane ring of cycloartanes between natural products. Based on this hypothesis comparison of 13 C-NMR spectral data of cycloartane structures and Compound 1a & 1b 13 C-NMR was performed which proved our hypothesis.
Without pentacyclic structure, based on DEPT and broadband (BB) decoupled 13 C-NMR, remaining carbons included three methyl groups (δ C 18.3, 17.2 and 17.6) , four methylene groups that two of them are involved in olefinic bonds (δ C 31.5, 31.6 &110.9, 111.4), three methane groups which two of them are oxygenated (δ C 35.9, 76.3 and 76.7), and a couple of quaternary carbons which are involved in double bonds (δ C 147.5 and 147.7). The duplicity of most of the signals of the side chain in 13 C-NMR spectra (δ C 17. 2 and 17.6, 31.5 and 31.6, 110.9 and 111.4, 76.3 and 76.7, 147.5 and 147 .7) and intensity of them in comparison with the remaining peaks with identical multiplicity and cycloartane structure as a basic skeleton, indicated a mixture of two epimers at one of side chain's carbon. In this way, taking one peak of each two close peaks, plus the remaining carbons of side chain, two methyl groups (2 CH 3 ), two methylene groups (CH 2 and =CH 2 ), two methine groups (CH and CHO), and an olefinic quaternary carbon were involved in side chain which totally were in accordance of
It is crystal clear that the quaternary olefinic carbon is linked to an olefinic methylene group and therefore, the 1 H-NMR signal for this methylene group has appeared as a doublet in the downfield region (δ H 4.88, J = 1.2 Hz and 4.93, J = 1.2 Hz). A triplet peak at δ H 4.02 (J = 6.4 Hz) with the integration equals to one hydrogen, demonstrated neighboring oxygenated methine to an aliphatic CH 2 in one side and quaternary carbon in another side. A singlet at δ H 1.72 cleared that one of methyl groups is connected to quaternary carbon while another methyl linked to a methine because of its doublet multiplicity (δ H 0.88, J = 5.6). Finally, considering these data and literature (8) compound 1a was determined as (24R)-cycloart-25-ene-3β,24-diol which was mixed with its C-24 epimer compound 1b; (24S)-cycloart-25-ene-3β,24-diol.
Cycloartane skeleton of compound 2 was confirmed by a pair of doublets in the up-field area (δ H 0.26, J = 4.0 Hz and 0.49, J = 4.0 Hz) and comparison of 13 C-NMR spectrum between compound 2 and 1a. In addition to the carbons which are involved in cycloartane skeleton, 13 C-NMR (BB and DEPT) showed three methyl groups, a methylene group and three methine group that two of them are olefinic (δ C 124. Broadband decoupled 13 C-NMR and DEPT spectra, showed the presence of an aliphatic backbone which has a double bond consisting of a methine group and a quaternary carbon (δ C 122.6 and 143.6 respectively). Furthermore, a carbonyl group (δ C 183.2) was identified that along with the loss of m/z = 45 at EI-MS indicated the presence of a carboxylic acid functional group. The presence of a carboxylic group and a double bond, taking together with seven degrees of unsaturation, supported pentacyclic structure of compound 4. Olefinic proton was clear with a triplet (δ H 5.28, J = 3.4 Hz) and seven angular methyl singlets at δ H 0.75, 0.77, 0.90, 0.91, 0.93, 0.99 and 1.13 were also detected. In addition, two signals at δ H 3.22 (1H, dd, J = 11.0, 4.6 Hz) and 2.82 (1H, dd, J = 13.8, 4.2 Hz) represented an oxygenated and allylic protons, respectively. These data suggested oleanene skeleton in agreement with EI-MS m/z 248 and 207, which are diagnostically important peaks of an olean-12-ene structure that has a carboxylic functional group at C-17, as a result of retro Diels-Alder cleavage of ring C (11). According to the previously reports for oleanolic acid (12, 13), the hydroxyl group on position 3 has a beta configuration and therefore the carbinolic hydrogen (H-3) has an alpha configuration. The 1 H-NMR signal for this hydrogen has appeared as a doublet of doublets (1H, J = 11.0, 4.6 Hz) at δ H 3.22. Eventually, the foregoing characteristics as well as published papers (12, 13), determined that compound 4 is oleanolic acid.
Anti-human immunodeficiency virus-1 (HIV-1) activity of compounds 1-4 and previously 
(4)
reported ingenoid diterpenes from E. erythradenia was evaluated ( Table 3 ). The four ingenoid diterpenes had been identified previously and their pro-apoptotic effect has been established (14). Based on literatures reports, it is known that oleanolic acid (15, 16) and cycloartanes (17) partially and diterpenes especially are cytotoxic which our results confirm that. But compound 7 showed significant anti HIV activity (IC 50 = 0.008 μM, CC 50 = 3.264 μM, TI = 380.64) and higher therapeutic index than nevirapine, a non-nucleoside reverse transcriptase inhibitor (NNRTI) used to treat HIV-1 infection and AIDS. Ligand protein binding studies between compound 7 and HIV-1 protease (pdb 1AJX) demonstrated three hydrogen bond between Asp 25A, Asp 28A and Gly 49A of protease active site and compound 7 with a distance less than 3 A° that can be responsible for the observed anti HIV-1 activity (Figure 3 ).