Design and synthesis new indole-based aromatase/iNOS inhibitors with apoptotic antiproliferative activity

The present study details the design, synthesis, and bio-evaluation of indoles 3–16 as dual inhibitors of aromatase and inducible nitric oxide synthase (iNOS)with antiproliferative activity. The study evaluates the antiproliferative efficacy of 3–16 against various cancer cell lines, highlighting hybrids 12 and 16 for their exceptional activity with GI50 values of 25 nM and 28 nM, respectively. The inhibitory effects of the most active hybrids 5, 7, 12, and 16, on both aromatase and iNOS were evaluated. Compounds 12 and 16 were investigated for their apoptotic potential activity, and the results showed that the studied compounds enhance apoptosis by activating caspase-3, 8, and Bax and down-regulating the anti-apoptotic Bcl-2. Molecular docking studies are intricately discussed to confirm most active hybrids’ 12- and 16-binding interactions with the aromatase active site. Additionally, our novel study discussed the ADME characteristics of derivatives 8–16, highlighting their potential as therapeutic agents with reduced toxicity.


Introduction
Aromatase, a member of the cytochrome P450 family (CYP19), is required for estrogen biosynthesis.CYP19 catalyzes a series of three hydroxylation reactions to transform C19 androgens (androstenedione and testosterone) into aromatic C18 estrogens (estradiol and estrone) (Graham Lorence et al., 1995).It is well established that elevated estrogen levels, both before and after menopause, promote hormone-dependent breast cancer (HDBC) and the metastasis of cancer cells to other parts of the body in women (Zhao et al., 2016).Aromatase inhibition reduces estrogen levels, which proves essential for controlling hormone-sensitive breast cancer.Aromatase is an important enzyme in the production of estrogen, especially after menopause.Aromatase inhibitors have been shown to be effective chemopreventive medicines in HDBC (Riemsma et al., 2010;Maghraby et al., 2023).
Indole is a valuable scaffold in medicinal chemistry because it is the building block for many important natural substances, such as serotonin, tryptophan, and tryptamine (Kaushik et al., 2013).Even though indole has a variety of therapeutic applications, the indole nucleus received increased attention in the AI class after the approval of Zindoxifene (compound VI, Figure 2), a strong antiestrogen (Fürstner et al., 1994).Given the importance of indole heterocycles in medicinal chemistry, Prior et al. synthesized and tested a series of novel 2-aryl indoles for aromatase inhibition (Prior et al., 2017).With an IC 50 value of 1.61 μM, compound VII (Figure 2) containing a nitrile (CN) group at indole's C-3 position was the most potent.The CN group lost its potency twice when it was moved to the indole's C-5 position (IC 50 = 3.34 μM).According to SAR, analogs with EWG substituted at position C-3 rather than C-5 exhibited excellent aromatase inhibition activity.
Pingaew et al. synthesized a series of bis-and tris-indoles with sulfonamides substituted with EWG (electron withdrawing group) and EDG (electron donating group) to develop effective AI (Pingaew et al., 2018).The bis-indole derivative (compound VIII, Figure 2) was discovered to be the most potent derivative, with an IC 50 value of 0.7 µM.The biological results showed that activity decreased when EWG was used instead of the methoxy group.Furthermore, it was discovered that bis-indole derivatives were more effective than tris-indoles.
Amoroso and colleagues developed indole aryl sulfonamides as effective AIs (Fantacuzzi et al., 2020).The most potent AI was derivative IX (Figure 2), which had an indole ring at the C-5 position and an IC 50 value of 0.16 µM.The potency is lost when aryl Structure of some nonsteroidal aromatase inhibitors (I-V).
sulfonamide is directly substituted with the 3-/6-position of indole.According to biological findings, the attachment position of the indole ring plays an important role in the aromatase inhibition activity of indole aryl sulfonamides.
On the other hand, nitric oxide (NO) is a small biomolecule that has various effects on tumor biology (Pervin et al., 2007).In general, depending on the tumor microenvironment and concentration, NO can influence cell proliferation, migration, and apoptosis (Burke et al., 2013;Shami et al., 1998).Endothelial (eNOS) and neuronal (nNOS) are two constitutive NOS isoforms; eNOS is crucial for maintaining vascular homeostasis, and its uncoupling is linked to the emergence of cardiovascular disorders (Pourbagher-Shahri et al., 2021).Conversely, the nNOS is primarily linked to learning and memory, smooth muscle relaxation, and synaptic plasticity (Wang et al., 1999).Aside from the constitutive NOS isoforms, an inducible NOS isoform (iNOS) is expressed in response to pro-inflammatory stimuli and is thus vital in the immune system.However, iNOS can be overexpressed in a variety of pathological conditions that result in uncontrolled NO production, including inflammatory bowel disease, rheumatoid arthritis, and cancer; thus, iNOS inhibition could be a therapeutic strategy for these diseases (Cinelli et al., 2020;Maccallini et al., 2020).In particular, iNOS overexpression is thought to be a predictor of poor outcomes in TNBC (Triple Negative Breast Cancer) patients, as it is associated with decreased relapse-free survival (Glynn et al., 2010), and iNOS can modulate many potential therapeutic targets for TNBC (Walsh et al., 2016).TNBC proliferation and migration were suppressed in vitro and in vivo by iNOS inhibition using L-NAME (compound X) and 1400 W (compound XI) (phase 1/2 clinical trials, Figure 3) (Granados-Principal et al., 2015).Furthermore, iNOS regulates tumor progression in an inflammatory environment by modulating the EGFR/MAPK pathway (Garrido et al., 2017).As a result, safer iNOS inhibitors with improved pharmacokinetic properties are desired to obtain new valuable therapeutic tools for treating breast cancer.
Additionally, apoptosis (programmed cell death) is a crucial control mechanism that triggers cell death when DNA damage surpasses the capability of the repair processes (Denmeade and Isaacs, 1996).Apoptosis, which is a normal factor of development, helps to control cell number and proliferation.Defects in apoptotic signalling contribute to several human disorders, including cancer.These flaws allow tumor cells to live longer than expected, reducing their reliance on exogenous survival factors and shielding them from oxidative stress and hypoxia, resulting in tumor growth and expansion (Cotter, 2009).These flaws allow for the accumulation of genetic changes that promote angiogenesis, disrupt cell proliferation, interfere with differentiation, and increase invasiveness during tumor formation.Therefore, restoring normal apoptotic equilibrium is a reliable practice for cancer treatment (Plati et al., 2008).
Frontiers in Chemistry frontiersin.org04 Al-Wahaibi et al. 10.3389/fchem.2024.1432920pyrazino [1,2-a]indol-1(2H)-ones (5-7, 9, 10-12, and 14-16).The newly synthesized compounds will be tested in vitro against a panel of four cancer cell lines as antiproliferative agents.All the newly synthesized compounds were then tested for aromatase inhibitory activity.The most potent aromatase inhibitors will be tested for their ability to inhibit iNOS.Moreover, the apoptotic potential of the most active derivatives was evaluated.Finally, docking analysis and ADME experiments were conducted for the most potent derivatives.
2 Results and discussion

Chemistry
Scheme 1 depicts the synthesis of the key intermediates and target compounds 3-7.Compound 1 was reacted with L-threonine benzyl ester 2 in the presence of BOP and DIPEA as coupling reagents, yielding compound 3 as a white solid after chromatography on silica gel with EtOAc/hexane (2:3) as an eluent.The oxidation of compound 3 with Dess-Martin periodinane reagent produced a crude product 4 that was used in the following step without further purification.Refluxing 4 in toluene with PTSA yields crude product purified by flash chromatography using EtOAc and n-hexane to afford 5 in 82% yield as a white solid.The structure of 5 was validated using 1 H NMR, 13 C NMR, and HRMS spectroscopy.The 1 H NMR spectrum of 5 revealed the pyrazine NH characteristic signal in the form of a singlet of one proton at δ 8.41 ppm corresponding to NH and a singlet signal of benzyl methylene protons at δ 5.36 ppm (CH 2 Ph).The spectrum also revealed a singlet signal of 3H at δ 3.19 ppm, corresponding to the CH 3 group, as well as protons from the ethyl group in the form of a quartet signal (2H) at δ 3.30 ppm and a triplet signal (3H) at δ 1.30 ppm.The 13 C NMR spectrum of 5 revealed two signals at δ 162 and δ 156 ppm for (ester C=O) and (NH-C=O), respectively, and aromatic carbons.The 5 structure was also validated by HRESI-MS with the [M + H] + ion at m/z 393.1013.Compound 5 was hydrolysed with aqueous NaOH to yield carboxylic acid 6.Carboxamide 7 was synthesized by coupling carboxylic acid 6 with benzylamine in the presence of DIPEA in DCM using BOP as the coupling reagent.
Compounds 6 and 7's structures were confirmed using 1 H NMR,  The synthesis of compounds 8-12 is outlined in Scheme 2. Compound 8 was synthesized using compound 1 and ethyl phenylalanine (compound 1a) as described in the general procedure for the synthesis of compound 3 (Scheme 1), yielding 8 as an oily product after purification by flash chromatography on silica gel short column using EtOA and n-hexane.Compound 9 is obtained by cyclizing 8 in toluene with PTSA.It is then purified using flash chromatography with EtOAc and n-hexane to yield 9 as a white solid in 75% yield.Various spectroscopic analysis methods were used to confirm the structure of compound 9.The 1 H NMR spectrum 9 revealed a pyrazine NH characteristic signal at δ 8.31 ppm and a methine proton multiplet signal at δ 4.81-4.79ppm (NHCH).The spectrum also revealed two doublet signals of one proton each at δ 3.29 and 3.10 ppm, corresponding to the CH 2 Ph group and ethyl and aromatic protons.Scheme 3 depicts the synthesis of compounds 13-16.Compound 13 was synthesized from compound 1 and methyl-2amino-2-phenylacetate (compound 1b) using the typical process for compound 8 synthesis (Scheme 2), producing 13 as a solid product.By cyclizing 13 in toluene with PTSA, compound 14 is obtained as a white solid.Compound 15 was obtained in 89% yield by reducing 14 in aqueous ethanol with NaBH 4 , followed by dehydration with PTSA in toluene to get 16.Various spectroscopic methods of analyses were employed to confirm the structures of compounds 13-16 (Supplementary Material).

Cell viability assay
This test investigates the effects of the new compounds 3-16 on normal cell lines to determine their safety.The normal human mammary gland epithelial (MCF-10 A) cell line was used to test the viability of the tested compounds.After 4 days of incubation on MCF-10 A cells within 50 µM of each investigated compound, cell viability was determined using the MTT assay (Gomaa et al., 2020;Marzouk et al., 2020).Table 1 outcomes demonstrate that none of the compounds examined were cytotoxic, and all compounds showed more than 89% cell viability at 50 µM.
Generally, the investigated compounds 3-16 had significant antiproliferative activity with GI 50 values ranging from 24 nM to 72 nM versus the four cancer cell lines evaluated, compared to the standard Erlotinib, which had a GI 50 value of 33 nM.Compounds 5,7,8,12,13,and 16 were the most potent six derivatives, with GI 50 values ranging from 24 nM to 38 nM, making compounds 5, 7, 12, and 16 (GI 50 = 32, 29, 24, and 26, respectively) more potent than Erlotinib (GI 50 = 33 nM), and even more potent than Erlotinib against breast (MCF-7) cancer cell lines.Compound 12 (R = 3-benzyl) was the most potent of the newly synthesized derivatives 3-16, with a GI 50 value of 24 nM, 1.4-fold higher than the reference Erlotinib (GI 50 = 33 nM).Based on the results, the type of substitution detected on the nitrogen atom at position two of the pyrazino [1,2-a]indol-1(2H)-one moiety appears to be crucial for activity.Compound 11 (R = 3-benzyl), which has the same backbone as compound 12 but with a free NH group at position 2, displayed a GI 50 of 41 nM (1.7-fold less potent than 12), demonstrating that the N-methyl group at position two is important for the antiproliferative activity.
Moreover, the 3-position substitution of the pyrazino [1,2-a] indol-1(2H)-one moiety is critical for activity.Compounds 7 (R = 3benzylcarboxamide) and 16 (R = 3-phenyl) had GI 50 values of 29 nM and 26 nM, respectively, indicating they were at least 1.4fold more effective than compound 11.These data indicate that the third position substitution significantly affects the antiproliferative activity of these compounds, with the highest activity shown in compounds with a benzyl substitution, followed by carboxamide and phenyl substitutions.
Another area of interest is the difference in efficacy between open-ring compounds and their closed-ring counterparts.Openring compounds 8 and 13 showed promising antiproliferative effects with GI 50 values of 38 nM and 35 nM, respectively, being equipotent to the reference Erlotinib (GI 50 = 33 nM), whereas compounds 3 and 4 showed weak activity with GI 50 values of 59 nM and 62 nM, respectively.Based on these findings, it is possible to conclude that the kind of functional groups and the type of substituents dictate the nature of these compounds' activity.Finally, Compounds 6, 10, and 15 exhibited low potency, as seen by their GI 50 values of 66, 54, and 72 nM, respectively.This suggests that having a high level of hydrophilicity does not confer any advantage for antiproliferative activity.

Assay for aromatase inhibitory action
Using ketoconazole and letrozole as reference medications, compounds 5, 7, 8, 12, 13, and 16 were examined for their ability to inhibit aromatase, a possible target for their antiproliferative effect (Maghraby et al., 2023).Table 2 displays the results as IC 50 values.The examined compounds 5, 7, 8, 12, 13, and 16 displayed significant aromatase inhibitory effect with IC 50 values ranging from 10 nM to 33 nM, in comparison to the reference drugs ketoconazole (IC 50 > 100 nM) and letrozole (IC 50 = 2 nM).The assay's results align with the findings of the antiproliferative assay, which demonstrated that compounds 7, 12, and 16, the most potent antiproliferative derivatives, were also the most potent aromatase inhibitors, with IC 50 values of 15 nM, 10 nM, and 12 nM, respectively.
Once again, compound 12 was the most potent aromatase inhibitor with an IC 50 value of 10 ± 1 nM, at least 10-fold more potent than the standard Ketoconazole but 5-fold less potent than the Letrozole.Compound 16 was the second most active aromatase inhibitor, with an IC 50 value of 12 ± 1 nM, like compound 12.Finally, compounds 5, 8, and 13 demonstrated good inhibitory action with IC 50 values of 20 ± 2, 33 ± 3, and 28 ± 2 nM, respectively.
These findings indicate that antiproliferative and anti-aromatase activity are linked and that any parameter that affects one has the same effect on the other.

Nitric oxide synthase inhibition
The L-citrulline assay method with fluorometric detection was used to test compounds 5, 7, 12, and 16 (most potent anti-aromatase derivatives) for their ability to inhibit iNOS (Carrión et al., 2023).The investigated compounds were tested at 1 µM against the iNOS and 10 µM against the eNOS to determine their isoform selectivity.Because of its critical role in cardiovasculature, inhibiting the constitutive enzyme eNOS should be avoided.Table 3 shows the results, which were expressed as enzyme percent inhibition compared to 1400 W as the positive control (100% inhibition at 10 µM).Results revealed that compounds 7, 12, and 16 were effective iNOS inhibitors, with enzyme percent inhibition values of 73, 87, and 81, respectively.Compound 7 had moderate inhibitory action against eNOS isoform, whereas compounds 12 and 16 were completely inactive against eNOS.Moreover, compound 5 was active against eNOS with a percent % inhibition value of 49%, even though it maintained a good inhibition of iNOS (62%), Table 3.
As a result, compounds 12 and 16 were chosen for IC 50 evaluation (Table 4).Interestingly, derivative 12 had an IC 50 value against iNOS of 0.075 µM, compared to reference 1,400 W (0.082 µM), as well as excellent isoform selectivity against eNOS (>670).Additionally, despite being less potent than reference 1400W, compound 16 demonstrated its efficacy against iNOS (IC 50 = 0.098 µM) while being inactive against eNOS.Based on their inhibition potency, derivatives 12 and 16 were thought to be possible aromatase inhibitors with iNOS inhibitory action.

Apoptosis assays
Compounds 12 and 16, the most potent derivatives in all in vitro investigations, were examined for their ability to initiate the apoptosis cascade and demonstrate proapoptotic activity.

Caspase-3 and caspase-8 expression levels
When cells get specific signal instructions, they undergo apoptosis, which results in several significant changes.Caspases, considered the main workers in apoptosis, are activated early in the process, cleaving critical cellular components such as nuclear proteins like DNA repair enzymes or structural proteins found in the cytoskeleton that are required for optimal cellular function.Caspases can activate DNases, which break down nuclear DNA (Chauhan et al., 2001;Kroemer and Reed, 2000).Compounds 12 and 16 were tested as caspase-3/8 activators against the human epithelial (A-594) cancer cell line (Youssif et al., 2019), with the results presented in Table 5.
In comparison to the untreated control cells, compound 12 treatment at its IC 50 concentration increased the expression levels of active caspase-3 and caspase-8 in A-594 cells by 9 and 29 folds, respectively (Table 5; Figure 5).Compound 16 treatment causes induction of caspase-3 and caspase-8 levels to be 8 and 25 times greater than those of untreated control cells.In all cases, compounds 12 and 16 were more effective caspase-3 and 8 activators than the control Staurosporine.

Proapoptotic BAX and anti-apoptotic Bcl-2 expression levels
Tumor cells can develop resistance to apoptosis through mutation or downregulation of proapoptotic proteins (e.g., Bax and Bak) or over-expression of anti-apoptotic proteins (e.g., Bcl-2 and Bcl-xL), which promote cell survival and are considered a hallmark in more than half of all cancers.This makes cancer cells more resistant to apoptotic stimuli and conventional cytotoxic anticancer drugs (Miyashita et al., 1994).
In the current study, treating lung (A-549) cancer cell lines with compounds 12 and 16 at IC 50 values increased pro-apoptotic Bax expression levels by 41 and 35 folds, respectively, while decreasing anti-apoptotic Bcl-2 expression levels by approximately 8 and 6 folds, Table 5 and Figure 5.As a result, Compounds 12 and 16 significantly increased the Bax/Bcl-2 ratio compared to the untreated control cells.

Docking study into aromatase
This research conducted a computational docking analysis to explore the binding interactions of the most active derivatives 7, 12, and 16 with the human placental aromatase cytochrome P450 (CYP19A1).Aromatase, a crucial enzyme in estrogen biosynthesis, belongs to the cytochrome P450 family and features a heme group at its catalytic centre, essential for its enzymatic activity.
The active site is located within a deep hydrophobic pocket, surrounded by key residues such as Met374, Arg115, Ala306, Val370, and Val373, which are critical for substrate binding and catalysis.These residues not only stabilize the binding of substrates and inhibitors but also play a significant role in the enzyme's selectivity and function.The careful arrangement of these residues around the heme group ensures precise catalytic activity, making aromatase a pivotal target in breast cancer therapy.The Discovery Studio software was used for this study, combined with the crystallographic structure of human placental aromatase cytochrome P450, which was co-crystallized with testosterone (PDB: 5JKW) (Ibrahim et al., 2020;Shaykoon et al., 2020;Bhat et al., 2022).
The molecular systems underwent energy minimization using the OPLS-AA force field, a method essential for ensuring conformational stability (Bernacki et al., 2005).Before starting the docking experiment, careful preparation of the protein structure was carried out to enhance accuracy.This included protonation steps, which significantly improved the reliability of the following docking analysis.The effectiveness of the docking process was validated by re-docking the co-crystallized testosterone into the active site of the aromatase enzyme, which yielded an S score of −8.43 kcal/mol.The S score is a measure used to assess the binding strength between compounds and their target protein, with lower scores indicating stronger binding affinities.Notably, this result was marked by a crucial hydrogen bond interaction between the hydroxyl group of testosterone and the amino acids Met374 and Arg115.Additionally, there were multiple interactions with Ala306, Val370, and Val373, highlighting their importance in ligand stabilization (Figure 6).
To further illustrate the binding interactions, cartoon and surface representations of the aromatase binding pocket are provided, highlighting the spatial arrangement of these critical residues and their interaction with the docked ligands.These visual models help in understanding how the structural features of the protein contribute to its interaction with different ligands (Figure 6).
The analysis of docking scores revealed that hybrid 12, which had the highest in vitro aromatase activity, showed a score of −7.86 kcal/mol.In contrast, hybrid 7, with the lowest in vitro aromatase activity, registered a score of −6.95 kcal/mol.Additionally, the docking score for compound 16, when tested against aromatase, was determined to be −7.55 kcal/mol.
Regarding the interactions between the compounds and the aromatase enzyme, compound 12 demonstrated significant interactions.Hydrogen bonds are fundamental in stabilizing the ligand within the active site of the enzyme, and their strength is largely dependent on both the bond distance and the bond angle.Caspase 3, 8, Bax, and Bcl-2 expression levels of compounds 12 and 16.
Frontiers in Chemistry frontiersin.org The carbonyl oxygen in its pyrazine ring acted as a hydrogen bond acceptor, forming a bond with crucial Met374 residue.The hydrogen bond distance between the carbonyl oxygen of compound 12 and the hydrogen of the Met374 side chain was measured at 2.50 Å.The bond angle formed was calculated to be approximately 167 °, indicating a nearly linear and therefore strong hydrogen bond, which is optimal for interaction stability.Additionally, the aromatic rings of compound 12 were involved in several hydrophobic interactions with key residues in the binding pocket.Notably, these interactions included pi-pi T-shaped interactions with Phe134 and alkyl/pi-alkyl interactions with Val370, Leu372, and Leu477.These non-covalent interactions collectively enhance the stability and binding affinity of compound 12 within the aromatase active site, Figure 7.The combination of these different interactions played a key role in the strong binding affinity of compound 12 within the active site of the aromatase enzyme.
Compound 16, like compound 12, showed significant interactions with the aromatase enzyme.In compound 16, the carbonyl oxygen of the pyrazine moiety acted as a hydrogen bond acceptor with the Arg115 residue.This hydrogen bond plays a crucial role in anchoring compound 16 within the active site.The hydrogen bond distance between the carbonyl oxygen of compound 16 and the hydrogen of the Arg115 side chain was measured at 2.66 Å.The bond angle formed was approximately 170 °, which indicates a strong and nearly optimal hydrogen bond for stabilizing the compound within the active site.Additionally, its aliphatic methyl group was involved in a pialkyl bond interaction with the Phe134 residue, as illustrated in Figure 8.These interactions will likely play a crucial role in stabilizing the compound within the binding pocket of the aromatase enzyme.
Similarly, pyrazoline oxygen of compound 7 establishes a hydrogen bond with the crucial residue Arg115 when it interacts with the aromatase enzyme, Figure 9.This hydrogen bond is significant as it represents the primary interaction stabilizing compound 7 within the active site.The hydrogen bond distance between the carbonyl oxygen of compound 7 and the hydrogen of the Arg115 side chain was measured at 2.76 Å.The bond angle was approximately 160 °, which supports a moderately strong hydrogen bond for stabilizing the compound within the active site.However, this is the only significant interaction observed for compound 7.The length of its amide bridge may not be optimal for a snug fit in the aromatase active site.Compound 7 also lacks the extensive hydrophobic interactions in compounds 12 and 16.
In conclusion, these results offer a basic insight into the potential inhibitory effects of the examined compounds and establish a foundation for further experimental validation and optimization in drug discovery.The detailed elucidation of these molecular interactions is crucial in directing the rational design of new therapeutic agents, particularly those targeting cancerrelated pathways.

ADMET studies
This supplementary section presents a comprehensive overview of the ADME (Absorption, Distribution, Metabolism, Excretion) studies for the newly synthesized compounds, specifically compounds 8-16.These studies are essential to evaluate the Cartoon and surface representations of the aromatase binding pocket docked with testosterone (PDB: 5JKW).The model shows the hydrogen bond interaction between the hydroxyl group of testosterones and the amino acids Met374 and Arg115, with additional interactions with Ala306, Val370, and Val373.Frontiers in Chemistry frontiersin.org12 compounds' pharmacokinetics and potential therapeutic efficacy (Al-Wahaibi et al., 2023d).The absorption studies show that all compounds have good absorption levels but vary in their Polar Surface Area (PSA), influencing their solubility and permeability.In terms of distribution, all compounds exhibit high plasma protein binding (>90%), which is crucial for their distribution in the bloodstream and free concentration.Additionally, the Blood-Brain Barrier (BBB) permeability ranges from very high to medium, indicating potential central nervous system activity for some compounds.The metabolism analysis reveals that several compounds (10-16) are inhibitors of the CYP2D6 enzyme, which could lead to significant drug-drug interactions and affect drug metabolism.The excretion studies indicate that most compounds have low solubility levels, potentially impacting their excretion rates.The Lipophilicity (AlogP98) values, ranging from moderate to high, suggest implications for membrane permeability and possible bioaccumulation Table 6.Overall, these ADME profiles offer valuable insights into the pharmacokinetic behaviors of the   Frontiers in Chemistry frontiersin.orgcompounds.The high plasma protein binding and variation in BBB permeability necessitate carefully evaluating their therapeutic window and CNS effects (Figure 10).The inhibition of CYP2D6 by several compounds highlights the need to assess potential drug interactions further.The solubility and lipophilicity profiles underscore the importance of optimization in future drug formulation strategies.This data is instrumental in guiding these compounds' clinical development and safety assessments, paving the way for their potential therapeutic application.The predicted ADME study of hybrids 8-16.Comprehensive profiles of absorption, distribution, metabolism, and excretion properties are depicted, including plasma protein binding, blood-brain barrier permeability, CYP2D6 inhibition, solubility, and lipophilicity (AlogP98).

Structure activity relationship (SAR)
In the present work, we report the design and synthesis of new indole-based derivatives able to inhibit the aromatase enzyme and/or iNOS isozyme.Compounds 12 and 16 were found to be the most intriguing since they are strong aromatase inhibitors with IC 50 values of 10 and 12 nM, respectively.Moreover, they are potent and selective iNOS inhibitors (IC 50 = 0.075 and 0.098 µM, respectively), and they were potential antiproliferative agents against the MCF-7 breast cancer cell line (IC 50 = 25 and 28 nM, respectively).Notably, at 50 µM concentrations, 12 and 16 demonstrated no effect on nontumor cells MCF-10A, indicating the promise tumor-cell selectivity of these derivatives.The molecular docking studies have successfully demonstrated the specific binding interactions of compounds 12 and 16, with the aromatase active site.This detailed analysis is pivotal in understanding their mechanism of action in breast cancer therapy.Furthermore, the comprehensive assessment of the absorption, distribution, metabolism, and excretion (ADME) properties of these hybrids underscores their potential as therapeutic agents.These insights, particularly their reduced toxicity profiles, are instrumental in their candidacy for effective breast cancer treatments.Based on the findings, 12 and 16 are promising prospects for the development of new agents for cancer therapy.Moreover, additional investigation into the mechanism of action, in vivo carcinogenic animal models, and lead optimization remains ongoing in our lab.

Synthesis of benzyl
To a solution of compound 3 (0.34 g, 0.82 mmol) in DCM (30 mL), Dess-Martin periodinane (0.45 g, 1.00 mmol) was added with stirring at 0 °C for 30 min, then the temperature was warmed to rt for 5 h.A saturated solution of Na 2 S 2 O 3 (10 mL) and NaHCO 3 (10 mL) was added to the reaction mixture by stirring at rt for 15 min.The organic layer was separated, dried over MgSO 4, and evaporated under reduced pressure to yield a crude product 4, which was used for the next step without further purification.

Assay of cell viability effect
The normal human mammary gland epithelial (MCF-10A) cell line was used to test the viability of the tested compounds.After 4 days of incubation on MCF-10A cells within 50 µM of each investigated compound, cell viability was determined using the MTT assay (Gomaa et al., 2020;Marzouk et al., 2020).Refer to Supplementary Appendix SA1 (supp.File) for more details.

Antiproliferative assay
The MTT assay was used to investigate the antiproliferative activity of 3-16 versus four human cancer cell lines using erlotinib as a control (Youssif et al., 2022;Mahmoud et al., 2022).See Supplementary Appendix SA1 for more details.Compounds 5,7,8,12,13,and 16 were tested for their capacity to inhibit aromatase, a potential target for their antiproliferative activity, using ketoconazole and letrozole as reference drugs (Maghraby et al., 2023).See Supplementary Appendix SA1 for more information.

Assay of nitric oxide synthase inhibition
The L-citrulline assay method with fluorometric detection was used to test compounds 5, 7, 12, and 16 (most potent antiaromatase derivatives) for their ability to inhibit iNOS (Carrión et al., 2023).Refer to Supplementary Appendix SA1 for more details.(Pawar and Rohane, 2021).We utilized the Protein Preparation Wizard to prepare the selected proteins for docking analysis.Following protein preparation, we meticulously mapped the ligands onto a three-dimensional model and conducted energy minimization using LigPrep.We generated a receptor grid tailored to the selected binding site to optimize potential binding interactions using the Receptor Grid Generation Tool.Subsequently, the Glide tool was employed to comprehensively assess both docking scores and the diverse binding modes exhibited by the ligands.This rigorous approach facilitated a thorough exploration of the binding affinities and interaction patterns between the ligands and the selected proteins, enhancing our understanding of their potential as inhibitors in breast cancer therapy.

In silico ADMET analysis
In our study, ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) studies were conducted using BIOVIA I Discovery Studio 2021 (Mohammed Ali et al., 2023).The chemical structures of all compounds were imported, and ADMET descriptors were predicted using integrated models.These models encompassed assessments based on Lipinski's Rule of Five and evaluations of absorption, distribution, metabolism, excretion, and toxicity.The resulting data were meticulously analysed to determine the drug-likeness and safety profiles of the compounds under investigation.
Validation, Visualization, Writing-original draft, Writing-review and editing.SB: Visualization, Writing-review and editing.

Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article.This work was funded by Princess Nourah Bint Abdulrahman University Researchers Supporting Project Number (PNURSP2024R3), Princess Nourah Bint Abdulrahman University, Riyadh, Saudi Arabia.The authors also acknowledge support from the KIT-Publication Fund of the Karlsruhe Institute of Technology.

FIGURE 3
FIGURE 3Structures of compounds X and XI.
10.74 ppm of pyrazine NH and the other at δ 9.05 ppm of NHCH 2 Ph, the characteristic signals of ethyl group in the form of a quartet at δ 3.25 ppm (2H) and triplet at δ 1.25 ppm (3H), and a singlet signal of CH 3 group at δ 2.73 ppm (3H).Additionally, the spectrum revealed the presence of aromatic proton-specific signals.HRESI-MS m/z of 7 calcd for [M-H] -C 22 H 19 ClN 3 O 2 : 392.1171, found: 392.1168.
The 13 C NMR spectrum of 9 revealed two signals for (CHC = O) and (NH-C = O), respectively, at δ 165.56 and δ 157.82 ppm and aromatic carbons.Compound 10 was obtained in quantitative yield by reducing 9 with NaBH 4 in aqueous ethanol.The 1 H NMR spectrum of 10 revealed the presence of a doublet signal at δ 7.69 (1H) corresponding to CHOH and a singlet signal at δ 4.89 of the hydroxyl group (CHOH).Moreover, the 13 C NMR revealed the presence of only one carbonyl carbon signal at δ 161.50 of (NH-C=O).HRESI-MS m/z of 10 calcd for [M + H] + C 20 H 20 ClN 2 O 2 : 355.1208, found: 355.1210 validated compound 10.Compound 12 was synthesized from 10 in two steps: dehydration with PTSA in toluene to yield 11, then alkylation with MeI in DMF in the presence of NaI to afford 12 in 88% yield as a white solid.The compound 12 1 H NMR spectrum revealed a singlet signal at δ 3.86 (2H, PhCH 2 ), while the CH 3 group signal overlapped with that of the ethyl group signal as a multiplet signal at δ 3.26-3.19(5H, CH 3 , CH 2 CH 3 ) and a triplet signal at δ 1.25 (3H, CH 2 CH 3 ).HRESI-MS with the [M + H] + ion at m/z 351.1257 also validated the 12 structure.

FIGURE 7
FIGURE 7 Docking representation models of compound 12 within the binding site of aromatase (PDB: 5JKW) (A) 2D-docked model of compound 12 (B) 3Ddocked model of compound 12.The carbonyl oxygen of pyrazine forms a hydrogen bond with Met374, and aromatic rings engage in hydrophobic interactions with Val370, Leu372, and Leu477.

FIGURE 8
FIGURE 8Docking representation models of compound 16 within the binding site of aromatase (PDB: 5JKW) (A) 2D-docked model of compound 16 (B) 3D-docked model of compound 16.The carbonyl oxygen of pyrazine forms a hydrogen bond with Arg115, and the aliphatic methyl group interacts with Phe134.

FIGURE 9
FIGURE 9 Docking representation models of compound 7 within the binding site of aromatase (PDB: 5JKW) (A) 2D-docked model of compound 7 (B) 3Ddocked model of compound 7.The carbonyl oxygen of pyrazine establishes a hydrogen bond with Arg115.
1. Open-ring (Scaffold A) compounds, namely, compounds 8 and 13, exhibited significant antiproliferative effects, comparable to the reference Erlotinib.However, compounds 3 and 4 displayed limited activity, indicating that the nature of functional groups and Scaffold B (closed-ring) compounds, the antiproliferative effects is determined by the following factors: a) The carbonyl oxygen in pyrazine ring (position-1) is essential for activity.It acted as a hydrogen bond acceptor, forming a bond with crucial amino acid Arg115.b) The type of substitution detected on the nitrogen atom at position two of the pyrazino [1,2-a]indol-1(2H)-one moiety appears to be crucial for activity demonstrating that the N-methyl group (R 1 = Me) at position two is important for the antiproliferative activity.c) The third position substitution significantly affects the antiproliferative activity of these compounds, with the highest activity shown in compounds with a benzyl substitution (R 2 = CH 2 Ph), followed by carboxamide (R 2 = CONHCH 2 PH) and phenyl (R 2 = Ph) substitutions d) A high level of hydrophilicity (R 2 = COOH and/or R 3 = OH) is not favoured for antiproliferative action.

TABLE 1
Cell viability assay and antiproliferative action of compounds 3-16.