Isolation of active constituents from plant material

Two plants named; A. annua L. (Voucher sample number AR001, Figures S1a and S1b) and V. agnus-castus L (Voucher sample number AR005, Figures S1c and S1d) were used in this study. Aerial parts of A. annua were collected from a cultivated area in Cairo, while V. agnus-castus fruits were collected from Giza zoo. Plant materials were identified by Mrs. Teresa Labib, Taxonomist at the El-Orman Botanical Garden, Giza. To complement the authentication of plants, samples were also compared with the Kew garden’s (https://www.kew.org/) reference samples revealing matched descriptions. A voucher sample of each plant was kept at the Pharmacognosy Department, Faculty of pharmacy, Alexandria University.

Fig. 1 2D ligand interaction diagram for docking poses of (a) androstenedione, (b) β-sitosterol and its (c) 3D interaction diagram in the active site of aromatase crystalline structure (3EQM). 

Dried aerial parts of A. annua (3.5 kg) were extracted with 70% ethanol, filtered and concentrated to obtain an 850 g residue. The hydroalcoholic extract (ethanol: water, 1:1 v/v) was partitioned with organic solvents of increasing polarities (starting with light petroleum then methylene chloride (CH₂Cl₂) followed by ethyl acetate and finally n-butanol) yielding 400, 62.8, 136.7 and 48 g, respectively. Twenty grams from the methylene chloride fraction was repeatedly chromatographed using normal silica gel columns (CH₂Cl₂ 100%, then gradual increase in the polarity by addition of methanol (MeOH) to 100%) to yield subfractions 21 (eluted from 0.5% MeOH in CH₂Cl₂) , 24 (eluted from 1% MeOH in CH₂Cl₂ ) and 40 (eluted from 20% MeOH in CH₂Cl₂)‎‎‎. Compound (1) was generated from subfraction 21 by crystallization and further purification on Sephadex LH20 column chromatography. Chromatographic purification of subfraction 24 (using light petroleum and CH₂Cl₂) followed by crystallization and purification on Sephadex LH20 yielded compounds (2) and (3), respectively. Compound (4) was obtained from subsequent normal phase chromatographic separation of subfraction 40 using MeOH and CH₂Cl₂ then crystallization. Similarly, fresh V. agnus-castus fruits (0.5 kg) were dried, ground and extracted with 70% ethanol. The extract was filtered and evaporated under reduced pressure. The residue obtained (21.6 gm) was suspended in ethanol: water 1:1 v/v then fractionated using the same solvents as for A. annua yielding 11.5, 6.9, 2.8 and 1.4 g, respectively. The CH₂Cl₂ fraction was subjected to silica gel column chromatography to afford 48 fractions, from which, selected subfraction 4 (eluted from 100 % CH₂Cl₂ ) was furtherly subjected to normal phase column chromatographic separation with n-hexane and ethyl acetate with increasing polarity, then fractionally crystallized to afford compounds (5) and (6). Crystallization of subfractions 13 (eluted from 2 % MeOH in CH₂Cl₂) and 33 (eluted from 15 % MeOH in CH₂Cl₂) using light petroleum and methanol yielded compounds (7) and (8), respectively.

Reagents, substrates and equipment

Potassium phosphate buffer, pH 7.4, stored at room temperature. Insect cell expressed baculovirus expressed aromatase + cDNA-expressed NADPH cytochrome P450 reductase 0.5 nmol/0.5 ml or 5 mg/ml (BD Supersomes™ enzymes, BD Biosciences), stored at -80^oC. Fluorometric substrate: Dibenzylfluorescein (Santa Cruz Biotechnology, Headquarters Canada). Prepared as 0.4 μM, stored at -20^oC. NADPH regenerating System: solution A, 5 ml (Corning® Gentest™). This solution contains 20 mg (26 mM) nicotinamide adenine dinucleotide phosphate (NADP+) and 20 mg (66 mM) glucose-6-phoshate, stored at -20⁰C. NADPH regenerating System, solution B, (Corning^® Gentest™) containing glucose-6-phosphate dehydrogenase (G6PDH) solution (40 U) in 1 ml of 5 mM sodium citrate, stored at -20^oC. Assay stop solution: 2 N sodium hydroxide, Stored at room temperature. Albumin as control protein was prepared at a concentration 0.38 mg/ml. Chrysin standard (Sigma Aldrich Co.) used as positive control for aromatase inhibitory activity.[¹]

Aromatase inhibition of the isolated compounds was evaluated by measuring the decrease in fluorescence intensity due to fluorescein, the hydrolysis product of dibenzylfluorescein (DBF), following the procedure previously described by Stresser [¹¹] using a 96 well plate. Briefly, the test sample (50 μl) was preincubated with 100 μl NADPH regenerating system (containing 1.5 ml cofactor solution A and 0.3 ml of cofactor solution B then completed to 15 ml with 0.5 M potassium phosphate, pH 7.4 buffer) for 30 min at 37 °C. The final concentration of each compound in the first well was 3.3 μM, which was used for preparing six further serial dilutions. After then, 100 μl of the enzyme and substrate mixture (containing 10 μl DBF, 40 μl of aromatase, and 0.38 mg/ml albumin completed to 10 ml with 50 mM potassium phosphate, pH 7.4) were added to incubate for an additional 30 min at 37 °C. Finally, the fluorescence of reaction solvent was measured using PerkinElmer HTS 7000 Bio Assay Reader with excitation at 485 nm and emission at 535 nm. Chrysin was tested as positive control.

Rigid and induced fit docking studies

The isolated compounds were subjected to rigid and induced fit docking using GLIDE and IFD (Prime module) programs of Schrodinger^® suite software.

Ligands were prepared using LigPrep module. They were optimized and energetically minimiz‎‎ed through OPLS 3 force field algorithm. Ionization states were predicted at pH 7 +/- 2 and Epik tool was selected to generate tautomers. Ligands were then desalted and specified chiralities were retained. Structures of unidentified stereochemistry are automatically subjected to all possible conformational changes for selection of best pose. The optimized structures were saved as maestro format (.mae) to retain 3D conformers. The crystal structure of human aromatase (3EQM) was first extracted from PDB and prepared by the aid of protein preparation wizard. The protein was first preprocessed where bond orders and hydrogens are assigned, zero order bonds to metals and disulphide bonds were created. Further, all the water molecules were deleted beyond 5 Å from the active site and iron states in heme at PH 7+/- 3 were generated where the one having the lowest state penalty was selected. Once the protein structures were reviewed and modified, Hydrogen bonds were assigned using PROPKA at PH =7 followed by energy minimization using OPLS 3 force field until the relative mean standard deviation (RMSD) of the minimized structure compared to the crystal structure reached 0.30 Ǻ. In the IFD protocol included refinement step where ligands were initially docked into the rigid protein using softened-potential Glide docking. The receptor grid was defined as an enclosing box at the centroid of workspace ligand (androstenedione) using an automated-generated box size (20 Å) with van der Waals radii scaling of 0.5 for both protein and ligands. Protein residues within a distance of 5.0 Å from any ligand poses were kept free to move, and the side chains were restrained, thus taking into account the protein flexibility. A maximum number of 20 poses were retained for each ligand. The protein backbone, side-chains as well as the ligand were energetically minimized using Prime program and the refined 20 ligand pose-complexes were scored by Prime energy. Ligands were re-docked into protein structures within 30 kcal/mol of the lowest energy structure obtained from the refinement step. The binding energy was computed as Induced fit score (IFD Score) for each ligand pose. The ligand-protein complexes scoring the lowest binding energy expressed by IFD Score, were subjected to deep interaction analyses using Maestro Ligand interaction 2-D and 3D diagrams.[⁷]

In-silico prediction of pharmacokinetic ADME&T properties

The drug-likeness behaviour of the isolated eight compounds from A. annua and V. agnus-castus was investigated through the evaluation of their pharmacokinetic (ADME&T‎) profiles using QikProp module embedded in Schrodinger^® Suite 2017 B software package. [¹³]

The compounds' structures were compiled from Pubchem database as sdf files. Structures were then imported and prepared by LigPrep module using OPLS3 force field algorithm. Herein, optimiz‎‎ation through energy minimiz‎‎ation and determination of ionization states at the specified PH (7.0 +/-3) were conducted. Compounds were saved as .mae files then imported to the qikprop wizard for ADME&T‎ prediction processed in normal mode. Physiochemical and pharmaceutical related descriptors were calculated for each compound where to assess its pharmacokinetic properties. The acceptable range for each descriptor was determined based on the values calculated for 95% of known drugs.

Aromatase inhibitory potential of the major constituents of A. annua and V. agnus-castus

Hydro-alcoholic extracts of each plant were fractionated using light petroleum, CH₂Cl₂, EtOAc and n-BuOH fractions. Further chromatographic purification of the CH₂Cl₂ fractions resulted in the isolation of eight compounds; artemisinin (1) [¹⁴], myricetin-3,7,4'-trimethyl ether (2) [¹⁵], scopoletin (3) [¹⁶] and domesticoside (4) [¹⁷] from A. annua L. and penduletin-4'-methyl ether (5) [¹⁸], casticin (6) [¹⁹], β-sitosterol (7) and β-sitosterol-3-O-β-D-glucopyranoside (8) from V. agnus-castus L.. The structures of isolated compounds were determined based on studying and analyzing their physical and chemical (Table S1), spectral data (Figures S2 to S7, Tables S2 to S5). as well as reviewing literature Structure elucidation of beta-sitosterol and the glycosidic derivative was estimated by two methods: Co-chromatography using reference standards and solvent systems; Toluene: ethyl acetate (7:3 v/v) for beta-sitosterol and Methylene chloride: methanol (9:1 v/v) for the glycosylated derivative. 2D- chromatographic runs were also performed for further identity confirmation. In addition, samples and reference standards solutions were analyzed by HPLC coupled with ELSD detector giving the same retention time values. Of note, myricetin-3,7,4’- trimethyl ether was reported for the first time from genus Artemisia.

Table 1 Aromatase inhibition of the tested compounds and positive control expressed as IC₅₀ values along with their rigid and flexible XP GScores. 

Name	IC50 (μM)	Rigid	Flexible
		XP GScore (Kcal/mol)	XP GScore (Kcal/mol)
β -sitosterol	0.13	-------	-9.61
Myricetin-3,7,4'-trimethyl ether	0.25	--------	-11.17
Domesticoside	0.39	-9.073	-10.58
Casticin	0.79	-4.912	-10.24
Penduletin -4'- methyl ether	1.00	-------	-9.47
Artemisinin	1.44	-4.858	-6.65
β -sitosterol-3- O- β -D-glucopyranoside	10.12	-------	---------
Scopoletin	42.57	-6.128	-7.95
Chrysin (5,7 dihydroxyflavone) (positive control)	1.25	-4.47	-8.55

The isolated compounds were tested for their aromatase inhibitory activity expressed as IC₅₀ values via constructing dose response curves (Figures S8 and S9). Results in Table 1 revealed that β-sitosterol was the most active compound possessing the lowest IC₅₀ values of 0.13 μM. On the contrary, scopoletin was considered the least active compound with an IC₅₀ value of 42.57 μM. Results also revealed that the isolated tested flavones exhibited lower IC₅₀ values compared to the positive control; chrysin (a known natural aromatase inhibitor) thus, could be considered as relatively potent aromatase inhibitors or drug leads. It was worthy to mention that this is the first report to test these isolated compounds for aromatase inhibition.

Molecular docking simulation studies

From Table 1, it was noted that the XP GScores obtained from IFD for artemisinin, scopoletin, domesticoside and casticin were more favorable (more negative) than their corresponding values in rigid docking. It was also interesting to observe that myricetin-3,7,4'-trimethyl ether, penduletin -4'-methyl ether and β-sitosterol exhibited good GLIDE score values using IFD even though they achieved no scores in rigid docking. This was probably attributed to the considered protein flexibility that allowed for better interactions and higher degrees of binding to the target protein.[²⁰] Interestingly, β-sitosterol-3-O-β-D-glucopyranoside showed no XP GScore in docking due to its non-sterically possible poses generated as adopted by GLIDE hierarchical sampling methodology.[²¹] IFD results were consistent with those obtained from the in-vitro assay where compounds exhibiting more negative docking scores were more likely to be more potent aromatase inhibitors as indicated by their lower IC₅₀ values.

As shown in Fig. 1, the binding mode of the natural substrate androstenedione (co-crystallised ligand) inside the enzyme's active site demonstrated the presence of Hem600 metal, hydrogen bond interactions with Met374 alongside hydrophobic interactions with some backbone residues as Ile133, Phe134, Phe221, Trp224, Ile305, Ala306, Val370, Leu372, Val373, Met374 and Leu477. Additionally, one positive charge interaction with Arg115, another for negative charge interaction with Asp309 and two polar interactions with the amino acid residues Thr310 and Ser478 were also observed. These interactions were in close agreement to those previously reported in literature.[⁷]

Insight analysis of the compounds' 2D interaction diagrams (Figures S10 and S11) indicated that their binding sites were in close proximity to those of androstenedione. This reflected that binding was successfully attained in the exact active site or pocket. Comparable interaction patterns to androstenedione could be observed for β-sitosterol [⁷] (Fig. 1), however, β -sitosterol engaged in extra nine hydrophobic interactions with Ile132, Phe148, Met303, Leu152, Pro308, Ala307, Phe203, Leu304 and Met311 beside one extra negative charged interaction with Glu302 and a polar interaction with Ser199. Based on both in-vitro and docking results, β -sitosterol exhibited the highest aromatase inhibitory activity with an IC₅₀ value of 0.13 μM and a GScore of -9.6 kcal/mol. These findings revealed that this compound may serve as a drug lead for designing potent aromatase inhibitors.

The ranking of the flavones from the view of both in-silico (IFD) and in-vitro studies was established as: myricetin-3,7,4'-trimethyl ether > casticin > penduletin -4'-methyl ether > chrysin. To relate between these data, understanding of the detailed binding modes of these flavones to their target receptor was performed seeking for a clue to their biological activities. Fig.(S10) illustrated the interaction patterns of the isolated flavones compared to both chrysin and vorozole (Fig. S11e) as known potent non-steroidal aromatase inhibitors where more hydrophobic interactions with backbone amino acid residues were observed. Comparative analysis of the isolated flavones’ interaction diagrams revealed the presence of two polar interactions with Asp309 and Glu302 in compounds myricetin-3,7,4'-trimethyl ether and casticin compared to only one with Asp309 in penduletin-4'-methyl ether. IFD conformations suggested that myricetin-3,7,4'-trimethyl ether had favorable hydrophobic interactions with some amino acid residues as well as three hydrogen bond interactions with Met374, Leu477 and Asp309, however, casticin possessed similar hydrogen bond interactions with more hydrophobic interactions and an extra metal co-ordination bond with the prosthetic Hem600. This activity order was corroborated by their IC₅₀ values (Table 1).

Numerous interactions of domesticoside with aromatase enzyme (Figures S11a and S11b) together with the GScore value (-10.58 Kcal/mol) verified its high anti-aromatase activity (IC₅₀ = 0.39 μM) in relation with the other compounds. 2D interaction diagrams of the weak inhibitory compounds; artemisinin and scopoletin are also presented in Figures S11c and S11d, respectively.

Pharmacokinetic ADME&T prediction

The ADME&T‎ properties of the active compounds were evaluated using QikProp which predicts some physiochemical properties. From the results tabulated (Table S6), it was observed that all compounds properties lied within the acceptable range for 95% of known drugs except for β-sitosterol having 6 stars hence, lying outside the acceptable range.

Prediction of oral absorption

The oral absorption process of a certain compound depends on many factors such as its solubility, permeability alongside its interactions with metabolizing enzymes and transporters in the gut wall.

QPlog S and CIQP log S were computed and used for determination of compounds' solubility. All Compounds except β-sitosterol and its β-D- glucoside were found to have acceptable water solubility values suggesting good solubility in the gut, higher contact with its walls and so higher absorption.[²²] Sub cut-off values of the above two compounds may be probably attributed to their large molecular weights. Likewise, all compounds showed variable degrees of permeability through the gut wall when calculating their QPPCaco (penetration across gut/blood barrier) values. For instance, β-sitosterol exhibited great QPPCaco value (3381.617) when compared to that of domesticoside (66.86).

Overall oral absorption of drugs is assessed by the following descriptors: qualitative human oral absorption, % human oral absorption and Jorgensen's “Rule of Three” (ro3).[²³] From all the tested compounds, only β-sitosterol and its β-D- glucoside violate one of these rules thus, might have some concerns regarding their oral bioavailability.

Prediction of plasma-protein binding (distribution)

The drug efficacy is directly related to its binding affinity to plasma proteins (ex: albumin). The more the drug binds to plasma proteins, the less free drug will be available to diffuse to the target organ. In this study, the predicted plasma-protein binding has been estimated by the QPlog KHSA; the prediction of binding to human serum albumin.[²⁴] With the exception of compound β-sitosterol, all compounds had values that were within the standard range of drug-like compounds. This might be attributed to the hydrophobic nature of the compound which contributed to its strong binding affinity to plasma proteins.[²³]

Prediction of first pass metabolism

QikProp also predicted the number of metabolic reactions the compound is subjected to when orally adminstered. That gave an indication of the accessibility of a compound to the target organ upon reaching the blood stream.[²⁵] The predicted values showed variations within the compounds, from 1 [as in Artemisinin] to 6 [as in myricetin 3,7,4' methyl ether, domesticoside, casticin and β-sitosterol-3- O-β-D- glucopyranoside, even though, all values were satisfactory.

Prediction of oral bioavailability

Oral bioavailability is defined as the part of drug dose that reaches circulation when orally administered. [²³] The oral bioavailability of a compound is mainly based on the process of its absorption from gut wall as well as its degree of metabolism in the liver (first pass metabolism). [²⁵] QikProp predicted some physiochemical properties such as molecular weight (mol_MW), solubility (QPlogS and CIQPlogS), number of hydrogen bonds donated to (donorHB) and accepted (accptHB) from water molecules in medium. These four molecular properties constitute the Lipinski rule of 5 which measures the drug-likeness or oral bioavailability of compounds i.e; a compound with a pre-determined biological activity would be considered orally active if it does not score more than 1 in Lipinski's rule. This rule gave indication for the compound pharmacokinetics in living organisms.[⁸] Results showed that all compounds were great drug like candidates showing acceptable values according to Lipinski rule with the exception of β-sitosterol-β-D- glucoside that violate 2 rules raising concerns regarding its bioavailability. Moreover, the optimum values of some descriptors such as polar surface area (PSA) and number of rotatable bonds (#rotor) had positive impact on the drug oral bioavailability.[²⁶] All compounds' values fell within the recommended range indicating high possibility for finding orally bioavailable candidates within these tested compounds.

Predicition of excretion profile

Ideal drugs are characterised by being safe, i.e; rapidly eliminated from the body without causing harm. Excretion process depends on the molecular weight and the partition coefficient (LogP) of the compound thus, small and hydrophilic compounds are more rapidly cleared through kidneys than large, hydrophobic ones.[²³] Accordingly, all compounds showed satisfactory LogP values, except for QPlogPo/w and QPlogPw values of compound β-sitosterol which might be due to its high degree of hydrophobicity.

Prediction of drug related cardiotoxicity

Clinical studies revealed that potassium channels denoted by HERG K⁺ - that control heart beating and electrical activity - contributed to the cardiotoxicity caused by many therapeutic agents as they block these channels leading to serious arrhythmia. Hence, predicting the IC₅₀ values of HERG K⁺ blockers using the parameter QPlogHERG might help in evaluating their in-silico cardiotoxicities.[¹³] Results showed that compounds are devoid of cardiotoxic side effects while, a special warning alarm should be directed to compounds domesticoside, penduletin -4'-methyl ether and β-sitosterol-3- O- β-D-glucopyranoside whose values were slightly below the cut off value (-5).

From the above calculated pharmacokinetic parameters, it was concluded that all compounds' properties - except for some of β-sitosterol - lied within the acceptable range for human use, thereby considered as drug leads for potent aromatase inhibitors.