Next Article in Journal
Risk Assessment by Presurgical Tractography Using Navigated TMS Maps in Patients with Highly Motor- or Language-Eloquent Brain Tumors
Next Article in Special Issue
Therapeutic Approach to Low-Grade Serous Ovarian Carcinoma: State of Art and Perspectives of Clinical Research
Previous Article in Journal
Breath Analysis for Early Detection of Malignant Pleural Mesothelioma: Volatile Organic Compounds (VOCs) Determination and Possible Biochemical Pathways
Previous Article in Special Issue
A Novel Role for NUAK1 in Promoting Ovarian Cancer Metastasis through Regulation of Fibronectin Production in Spheroids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Azasteroid Alkylators as Dual Inhibitors of AKT and ERK Signaling for the Treatment of Ovarian Carcinoma

by
Panagiotis Dalezis
1,
Eleni Geromichalou
1,
Aikaterini Polonifi
1,
Sofia Sagredou
1,
Nikolaos Nikoleousakos
1,
Michael Nikolaou
1,
Vasiliki Sarli
2,
Mihalis I. Panayiotidis
3,4,* and
Dimitrios T. Trafalis
1,*
1
Laboratory of Pharmacology, Medical School, National & Kapodistrian University of Athens, 11527 Athens, Greece
2
Department of Chemistry, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
3
Department of Electron Microscopy & Molecular Pathology, The Cyprus Institute of Neurology & Genetics, 2371 Nicosia, Cyprus
4
The Cyprus School of Molecular Medicine, The Cyprus Institute of Neurology & Genetics, 1683 Nicosia, Cyprus
*
Authors to whom correspondence should be addressed.
Cancers 2020, 12(5), 1263; https://doi.org/10.3390/cancers12051263
Submission received: 12 March 2020 / Revised: 10 April 2020 / Accepted: 24 April 2020 / Published: 16 May 2020
(This article belongs to the Special Issue Preclinical and Clinical Advances in Ovarian Cancer)

Abstract

:
(1) Background: Previous findings show that lactam steroidal alkylating esters display improved therapeutic efficacy with reduced toxicity. The aim of this study was to evaluate the anticancer activity of two newly synthesized aza-steroid alkylators (ENGA-L06E and ENGA-L08E) against human ovarian carcinoma cells, and consequently, the dual inhibition of RAS/PI3K/AKT and RAS/RAF/MEK/ERK signaling pathways, both of which are closely associated with ovarian cancer; (2) Methods: The in vitro cytostatic and cytotoxic effects of ENGA-L06E and ENGA-L08E were evaluated in a panel of five human ovarian cancer cell lines, as well as in in vivo studies. ENGA-L06E and ENGA-L08E, in addition to another two aniline-mustard alkylators, POPAM and melphalan (L-PAM), were utilized in order to determine the acute toxicity and antitumor efficacy on two human ovarian xenograft models. Also, in silico studies were performed in order to investigate the dual inhibition of ENGA-L06E and ENGA-L08E on RAS/PI3K/AKT and RAS/RAF/MEK/ERK signaling pathways; (3) Results: Both, in vitro and in vivo studies demonstrated that ENGA-L06E and ENGA-L08E were significantly more effective with a lower toxicity profile in comparison to POPAM and L-PAM alkylators. Moreover, in silico studies demonstrated that the two new aza-steroid alkylators could act as efficient inhibitors of the phosphorylation of AKT and ERK1/2 molecules; and (4) Conclusions: Both ENGA-L06E and ENGA-L08E demonstrated high anticancer activity through the inhibition of the PI3K-AKT and KRAS-ERK signaling pathways against human ovarian carcinoma, and thus constituting strong evidence towards further clinical development.

1. Introduction

Ovarian carcinoma (OC) is the most lethal gynecological malignancy due to drug resistance in current therapeutic approaches. There are several histological types of OC that are characterized by distinct molecular genetic alterations associated with different histopathological features. On the other hand, the understanding of the culpable molecular pathways is necessary for the development of novel therapeutics. Nowadays, the therapeutic strategy depends more on tumor stage and grade than on the histological type. The World Health Organization (WHO) classifies ovarian neoplasms according to their histological differentiation, namely epithelial, sex cord-stromal, and germ cell neoplasms, where epithelial ovarian tumors constitute the most common form of ovarian carcinomas [1,2]. In addition, non-epithelial ovarian cancers represent 10–15% of all ovarian cancers and include diverse tumors of germ cell origin, sex cord-stromal cell origin, as well as extremely rare types of ovarian cancers, such as small-cell carcinomas and sarcomas [3]. In platinum-sensitive OC patients or those who recur ≥6 months from the initial diagnosis of the disease, treatment has traditionally consisted of platinum-based chemotherapy. Unfortunately, there is a considerable number of patients who show resistance to platinum-based chemotherapy and so current research efforts are focused on the discovery of new therapeutic options. The majority of advanced-stage OC patients, even those who respond well to frontline therapy, will ultimately recur and succumb to their disease [1,2]. Consequently, new therapeutic options are required targeting the discovery of new drug compounds. To this end, the PI3K/AKT/mTOR signaling pathway is one of the most thoroughly investigated ones as it is implicated in several cellular cascades, including growth, motility, survival, proliferation, protein synthesis, autophagy, transcription, as well as angiogenesis. Its deregulation has been shown in several tumors, including ovarian cancer. In this setting, targeted inhibitors have been introduced, while clinical trials have assessed their safety in ovarian cancer patients. Meanwhile, the ongoing phase I and II studies are evaluating the oncologic outcome of these compounds in combination with conventional chemotherapy [4,5].
The discovery of inhibitor molecules of the mammalian rapamycin (mTOR), poly-ADP-ribose polymerase (PARP), MAPK/ERK, and components of the epidermal growth factor receptor (EGFR) pathway are on the frontline of clinical research. In ovarian cancer, important hormones using the extracellular signal-regulated kinase (ERK) pathway (e.g., gonadotropins and gonadotropin-releasing hormones) activate the KRAS-BRAF-MEK axis, and consecutively, ERK which can then activate transcription factors, including myc or elk-1. On the other hand, the mitogen-activated protein kinase (MAPK) pathway develops oncogenic potential through permanent signaling activity as a result of the Kirsten rat sarcoma viral oncogene (KRAS) or v-raf murine sarcoma viral oncogene homolog B (BRAF) mutations, thus inducing cell proliferation. In addition, activation of ERK can activate downstream protein kinases or transcription factors that enhance tumor development [2]. Overall, the MAPK/ERK pathway can contribute to therapy-induced tumor-growth suppression. For instance, ERK-mediated mitogen-activated protein kinase (MAPK) phosphatase-1 (MKP-1) expression is involved in cisplatin resistance, and thus, its targeting could potentially overcome cisplatin resistance in patients with ovarian cancer [6]. Two newly synthesized homo-aza-steroidal alkylating esters of the p-bis-(2-chloroethyl) aminophenoxy-propionic acid (POPAM) were designed in an attempt to improve the therapeutic efficacy in tumors resistant to classical and non-classical alkylating drugs. In general, the lactam steroid alkylators combine two active compounds in a single molecule by incorporating modified steroids, that bear a lactam (–NH-CO-) chemical moiety in one or more steroid rings, and thus functioning as cytotoxic anticancer agents. These hybrid aza-steroid alkylators produce more potent antitumor effects, depending on the alkylating moiety, as well as on the modified steroidal carrier bearing the lactam group, thus increasing anticancer activity significantly and improving the therapeutic ratio while reducing overall systemic toxicity. Finally, these molecules produce a multi-targeted antitumor effect by inhibiting various cellular cascades important for the survival of cancer cells [7,8,9,10,11,12].
In the present study, we have evaluated the cytostatic and cytotoxic effect of two new aza-steroid alkylators (ENGA-L06E and ENGA-L08E) against aniline-mustard alkylators POPAM and melphalan (L-PAM) in a panel of five human ovarian cancer cell lines, as well as on two human ovarian tumor xenograft models. More specifically, we have investigated their inhibitory effects on the RAS/PI3K/AKT and RAS/RAF/MEK/ERK signaling pathways and whether their dual inhibition can potentiate antitumor efficacy by utilizing in silico computational methods.

2. Results

We have tested the new hybrid homo-aza-steroidal alkylators against five human ovarian adenocarcinoma cell lines, in vitro, in comparison with the azasteroid esters’ alkylating moiety POPAM. The two conjugates, ENG-L06E and ENGA-L08E, displayed significantly higher anticancer activity, cytostatic, and cytotoxic effects on the tested cell lines than the non-steroidal alkylator POPAM. ENGA-L06E demonstrated rather more potent, cytostatic, and cytotoxic effects than the ENGA-L08E. Furthermore, TP53 mutations, in all human ovarian cancer cell lines, didn’t affect the cytostatic or cytotoxic effects of ENG-L06E and ENGA-L08E. Similarly, KRAS (in OVCAR-5), MEK1 (in OVCAR-3), PIK3CA (in SK-OV-3) mutations, as well as the mismatch repair (MMR) status didn’t invoke the cytostatic or cytotoxic effects of ENG-L06E and ENGA-L08E significantly. In addition, the cytostatic and cytotoxic impact of ENG-L06E and ENGA-L08E was shown to be independent of the steroid receptor expression in all cell lines (Table 1). Both ENG-L06E and ENGA-L08E were more effective against UWB1.289 cells, by means of exerting significantly higher cytostatic and cytotoxic activity, which are known to be Breast cancer type 1 (BRCA1)-null (mutated, lack wild-type BRCA1 gene), in comparison to the other tested cells bearing either the BRCA1 wild type (w.t) or where at BRAC1 is restored (UWB1.289 + BRCA1 cells). Finally, both ENG-L06E and ENGA-L08E were fairly effective against OVCAR-5 and SK-OV-3 cells (both of which exert low sensitivity to the alkylator cisplatin) as well as OVCAR-3 cells, which are (i) resistant to the aniline mustard alkylator melphalan, (ii) not sensitive to cisplatin and neither to Adriamycin [13,14,15,16] (Table 1).
The analytical results of the inhibitory effects of ENGA-L06E and ENGA-L08E on ERK1/2 and AKT1/2 phosphorylation, in all studied cancer cell lines, and the corresponding kinetics according to the Michaelis Menten model are illustrated in Figure 1 and Figure 2, as well as in Table A1 and Table A2. According to the results, both compounds prevented the phosphorylation of ERK1/2 and AKT in all cell lines tested, effectively. The results on pERK1/2 and pAKT inhibitory effect induced by aza-steroidal alkylators ENGA-L06E and ENGA-L08E demonstrated similar patterns of kinetics (Vmax, Km). The highest effect on pERK1/2 inhibition resulted at 60 min and intervened in a dose and time-dependent relevance. In all the tested cell lines, taking together the kinetic parameters and IC50s (the concentrations of tested compounds that induce 50% inhibition of phosphorylation), the ENGA-L06E was shown to be more effective than ENGA-L08E in inhibiting ERK1/2 phosphorylation (Figure 1, Table A1). However, both compounds provoked higher inhibition of pERK1/2 in UWB1.289 cells by about 85% at 60 min. On the other hand, the inhibition of AKT protein phosphorylation caused by ENGA-L08E was somewhat more potent than the inhibitory effect of ENGA-L06E, in all tested cell lines. The highest effect on pAKT inhibition resulted at 60 min and took place in a dose and time-dependent manner, similar to that of ERK1/2 phosphorylation (Figure 2 and Table A2). Overall, ENGA-L08E exhibited a higher inhibitory effect on pAKT in OVCAR-3 cells while ENGA-L06E in SK-OV-3 and UWB1.289 cells (Figure 2). The inhibitory activity of both compounds on ERK1/2 and AKT phosphorylation was not affected significantly by mutations on TP53, PIK3CA, microsatellite instability (MSI) status, as well as steroid receptor expression and sensitivity to cytotoxic drugs. However, KRAS mutation in OVCAR-5 cells may affect AKT phosphorylation, while MEK1 mutations, in OVCAR-3 cells, may derogate the inhibitory effect on ERK1/2 phosphorylation (Table A1 and Table A2). Interestingly, both ENG-L06E and ENGA-L08E show an optimal inhibitory effect, on both ERK1/2 and AKT phosphorylation, in BRCA1-null UWB1.289 cells (Table A1 and Table A2).
The inhibitory effect of ENGA-L06E and ENGA-L08E on the phosphorylation status of ERK1/2 and AKT, on both ovarian cancer cell lines, is presented in Figure 3. More specifically, the inhibition curves demonstrate t levels of percent expression of ERK1/2/GAPDH and AKT/GAPDH, which reflect more accurately the inhibition of pERK1/2 and pAKT in treated cells. In more detail, the phosphorylation of ERK1/2 and AKT was inhibited under the treatment of both compounds at 25 and 50 μΜ. As the inhibition curves depict, the two tested hybrid molecules yielded approximately 50–70% inhibition at 25 μΜ whilst phosphorylation of targeted proteins was almost totally inhibited (>85%) at higher concentrations (50 μΜ). Moreover, the inhibition of phosphorylation of ERK1/2 and AKT starts in 1 h while it gradually increases in 2 h after which it leads to a plateau (Figure 3). Overall, ENGA-L08E induced a higher inhibition of AKT phosphorylation than ENGA-L06E, whereas ENGA-L06E demonstrated a higher inhibition on ERK1/2 phosphorylation than ENGA-L08E (in both cell lines (Figure 3). Finally, preliminary experimental data (not shown) have shown that the tested azasteroid alkylators inhibit phosphorylation of AKT and ERK in the tumor xenograft tissue in vivo. More detailed quantitative studies on the inhibition of AKT and ERK phosphorylation and the effects on the expression of other predictive biomarkers on tumor cells, such as PD-L1, DNA mismatch repair (MMR) system, in correlation with the treatment schedules, time, and dose administration of the azasteroid alkylators in vivo, are currently under investigation.
ENGA-L06E, ENGA-L08E, L-PAM, and POPAM were evaluated for their in vivo acute toxicity, and such values are illustrated inTable 2. The conjugation of the alkylating agent with the aza-steroidal moiety resulted in significantly reduced acute toxicity effect and an increase of the LD10 and LD50 values. The LD50 and LD10 values induced by ENGA-L06E and ENGA-L08E were significantly decreased compared to POPAM and L-PAM. Moreover, ENGA-L08E had a significantly higher acute toxicity effect than ENGA-L06E (Table 2).
In addition, ENGA-L06E and ENGA-L08E were administrated i.p. at the LD10 dose and caused significant reduction in the size of SK-OV-3 (Figure 4) and OVCAR-5 (Figure 5) human xenograft tumor in mice. More specifically, 23 and 44 days after implantation, as well as 10 (d-10) and 16 days (d-16) after drug administration, the sizes of the primary tumors were 928 ± 16.4 mm3 (Figure 4) and 836 ± 67.4 mm3 (Figure 5) for the control group compared to 390 ± 21.4 mm3 and 478 ± 18.3 mm3 (Figure 4), as well as 113 ± 14 mm3 and 220 ± 19.1 mm3 (Figure 5) for the groups receiving ENGA-L08E and ENGA-L06E, respectively.
Furthermore, treatment with ENGA-L08E and ENGA-L06E significantly increased the lifespan of SK-OV-3 as well as of OVCAR-5-bearing mice (Table A5 and Table A6). On the other hand, POPAM and L-PAM showed no significant antitumor activity on SK-OV-3 (Table A3) and OVCAR-5-bearing mice (Table A4). Finally, the tumor response and anticancer activity assessment of ENGA-L06E, ENGA-L08E, POPAM, and L-PAM, according to clinical criteria, demonstrated a significant increase of overall survival (OS), progression-free survival (PFS), and 100% response rate (RR) for ENGA-L06E and ENGA-L08E, respectively. On the other hand, treatment with POPAM and L-PAM showed no significant accrual of OS, PFS, and 100% RR in either SK-OV-3 or OVCAR-5-bearing mice (Table A5).
The crystallographic structures of the proteins were used during the induced fit docking (IFD) simulations of hybrid steroid organic molecules ENGA-L06E and ENGA-L08E. Their crystal structures were prepared using the Protein Preparation Wizard [17], missing loops and side chains were prepared using Prime, ionized at pH 7.4 using PROPKA [18], and minimized using the OPLS3 force field [19]. The structures were designed using Maestro 10.2 and ionized at the target pH 7.4 using the Ligprep module [20]. Docking simulations were performed in order to elucidate the mechanism of action of the studied compounds at the molecular level, providing an interpretation of their biologic activity as anticancer drugs. The first step was to evaluate the active site for ERK1 and ERK2, a site that is targeted by the same co-crystallized ligand 38Z0, for both proteins (Figure 6A). A grid was generated from the ERK1/2 co-crystallized inhibitor, and the 38Z0 was docked, resulting respectively, in satisfactory overlapping with the crystallographic data. Hydrogen atoms of ligands were omitted for clarity (Figure 6B). ENGA-L06E was found to be less active in both ERK1 and ERK2 than the co-crystalized ligand 38Z0, whereas ENGA-L08E was more active than the ENGA-L06E. In addition, ENGA-L08E bound stronger than the co-crystalized ligand 38Z0 to ERK2 protein. It was also noticeable that in the ERK1 protein, the molecule ENGA-L08E bound efficiently compared to 38Z0 (Table A6).
Stereochemistry is an important physicochemical property capable of explaining the differences in binding capacity. To this end, ENGA-L08E (with the hot pink color of C atoms) enters the cavity with the same direction as the 38Z0, in contrast to ENGA-L06E that approaches the binding site from the alkylating agent moiety while the steroidal moiety of lactam testosterone remains out of the cavity. Generally, both ENGA-L06E and ENGA-L08E (consisting of a modified steroidal skeleton) came from precursor steroids, testosterone and estrone, respectively, conjugated with the same alkylating agent. The former, due to its aromatic ring and the carboxyl group, which functions as a linker, increase the binding capacity of ENGA-L06E and ENGA-L08E via bonding to residues Lys71 with Halogen and pi-cation bond, and Asp184 with aromatic H-bond. Therefore, the binding efficiency of ENGA-L08E is influenced by the presence of estrone’s aromatic ring in contrast to ENGA-L06E, in which the precursor steroid was testosterone (Figure 7A,B). On the other hand, we evaluated the active site of AKT1 and AKT2 (in the grey grid), which are targeted by two different co-crystallized ligands, 0R4 and X39, respectively. In doing so, we overlapped the two proteins, and AKT1 (the dark green color) was found to be slightly bigger than the AKT2 (the yellow color), while there was a clear difference in the location of the active site where the co-crystallized ligands were located (Figure 8A). In addition to the docking results for both proteins (Table A7), we decided to emphasize on AKT1 protein’s cavity (Figure 8B) given that both ENGA-L06E and ENGA-L08E seem to be inactive for AKT2 (Table A7). The specificity of the active center (Figure 9A), in contrast to the active center of ERK1/2 protein, lies in the fact that it crosses the protein, and in order to inhibit it, the tested molecules should fit into the cavity the same way the co-crystallized ligand 0R4 does (Figure 8B; Table A7).
After the evaluation of the active site for AKT1, we generated a specific grid for AKT1/2 in a similar way as described above. Subsequently, the hybrid molecules ENGA-L06E and ENGA-L08E were docked in the identified cavity and compared to the co-crystallized ligand OR4. The most active molecule appeared to be ENGA-L08E, whereas ENGA-L06E was also considered as comparable to the co-crystallized ligand (Table A7). A closer view of the docking solution revealed an extended hydrophobic region, as defined by residues Asn54, Gln59, Lys268 for ENGA-L08E and Gln203, Ser205, Thr211, Lys268 for ENGA-L06E (Figure 9C). Although ENGA-L06E further appears to make more contacts with the residues in the active site of the protein, at the same time, it presents slightly higher energy than ENGA-L08E (Table A7) because of the lack of hydrogen bonds when compared to ENGA-L08E (Figure 9C). The superimposition of the three molecules: 0R4 (in green color), ENGA-L06E (in purple color) and ENGA-L08E (in the hot pink color) revealed that the alkylating moiety of ENGA-L08E, and by extension the whole molecule, is precisely adjusted above the co-crystallized ligand in contrast to the ENGA-L06E molecule (Figure 9B). More specifically, in Figure 10A,B, it is clearly represented how the ENGA-L08E molecule is adapted into the AKT1 protein’s cavity. The specificity of this particular active center (Figure 9A), in contrast to the active center of ERK1/2 protein, lies in that it crosses the protein and in order to inhibit it, the tested molecules should fit into the cavity the same way the co-crystallized ligand 0R4 fits (Figure 8B; Table A7). Subsequently, the hybrid molecules ENGA-L06E and ENGA-L08E were docked in the identified cavity and compared to the co-crystallized ligand. The most active molecule appeared to be ENGA-L08E, while ENGA-L06E was also considered comparable to the co-crystalized ligand (Table A7). A closer view of the docking solution in Figure 9C reveals an extended hydrophobic region, as defined by residues Asn54, Gln59, Lys268 for ENGA-L08E, and residues Gln203, Ser205, Thr211, Lys268 for ENGA-L06E. Although ENGA-L06E appears to make more contacts with the residues in the active site of the protein, at the same time, it presents slightly higher energy than ENGA-L08E (Table A7) probably because of the lack of hydrogen bonds (Figure 9C). Complementary to the above data, the superimposition of the three molecules 0R4 (in green color), ENGA-L06E (in purple color) and ENGA-L08E (in the hot pink color) revealed that the alkylating moiety of ENGA-L08E, and by extension the whole molecule, is precisely adjusted above the co-crystalized ligand, in contrast to the ENGA-L06E molecule (Figure 9B). More specifically, from the backside of the protein, it is clearly represented howtheENGA-L08E molecule is adapted into the AKT1 protein’s cavity (Figure 10B).

3. Discussion

The new homo-aza-steroidal alkylating esters of POPAM were designed in the context of improving therapeutic efficacy against tumors resistant to alkylating drugs [7,9,21]. Lactam steroid alkylators (LSAs) combine two active compounds in a single molecule, incorporating modified steroids that bear the lactam (–NH-CO-) chemical moiety in one or more steroid rings, with the capacity to act as cytotoxic anticancer agents. These hybrid molecules were aimed to increase anticancer activity, reduce systemic toxicity, and improve the overall therapeutic ratio. Furthermore, LSAs induce a multi-targeted antitumor effect by inhibiting the function of critical cellular pathways important for cancer cell survival cancer cells [7,8,9,10,11,12,22,23,24,25].
Of particular importance, RAS/RAF/MEK/ERK and PI3K/PTEN/Akt/mTOR are two signaling pathways capable of regulating a wide array of cellular functions, including proliferation, differentiation, survival, metabolism, and migration. Various studies indicate that AKT and ERK are significantly implicated in cancer cell survival and that their synergism promotes cancer progression and metastasis, thus leading to poor therapeutic outcomes and tumor recurrence. Indeed, RAS/RAF/MEK/ERK and PI3K/PTEN/Akt/mTOR are important pathways. Their activation occurs in approximately 70% of all ovarian cancers. Furthermore, PI3K catalytic subunit alpha (PIK3CA) gene amplification has also been associated with genomic instability, p53 mutation, and a lack of response to chemotherapy. Evidence from ovarian cancer cell lines and animal models revealed that activation of the PI3K/AKT pathway might lead to chemotherapy resistance. Suppression of apoptosis contributed to platinum and taxane resistance. The chemotherapy-induced apoptosis has been restored by specific PI3K inhibitors in vitro and in vivo [26]. Consequently, targeting the regulation of downstream functions by AKT and ERK pathways is of great importance and therapeutic relevance in cancer [27,28]. Even though the downregulation of phosphorylated ERK or AKT may induce cell apoptosis or cell cycle arrest, the crosstalk between these two signaling pathways may counteract the effect of single inhibition. For instance, it has been shown that inhibition of the MEK/ERK pathway enhances AKT phosphorylation while ERK activation attenuates AKT phosphorylation. Therefore, inhibition of either pathway alone may not contribute significantly to the context of an anticancer therapy, whereas inhibition of both pathways may interrupt the transduction of survival signal(s), and consequently, achieve a significant therapeutic effect [29,30]. Previous research has demonstrated that the constitutive activation of AKT in various ovarian cancer cell lines (e.g., SK-OV-3, OVCAR-4, and OVCAR-5) promotes cell survival and resistance to cisplatin. Thus, Akt is a potential therapeutic target as its inhibition might increase cytotoxic sensitivity selectively in ovarian cancer that has a dependence on AKT kinase activity for survival [30]. On the other hand, therapies targeted towards MAPK/ERK components have various response rates in solid tumors, including ovarian cancer. The mechanisms responsible for the acquired chemoresistance in ovarian tumors are many, including efflux pumps, repair mechanisms, signaling survival pathways (e.g., PI3K/AKT, MAPK, estrogen signaling), etc. Specifically, PARP1 inhibition causes a loss of ERK2 stimulation by decreasing the activity of critical pro-angiogenic factors, including the vascular endothelial growth factor (VEGF) and the hypoxia-inducible factor (HIF) [17]. Thus, many natural and synthetic agents capable of suppressing the MAPK/ERK signaling pathway have been tested in order to overcome chemoresistance, and thus a series of targeted inhibitors have emerged [6,31]. Previous studies have shown that LSAs induce an anticancer response through the combination of the direct induction of cellular DNA damage with the inhibition of PARP enzyme(s), and consequently, the inhibition of DNA repair [12,25]. To these ends, inhibitors of nuclear PARP enzymes (e.g., PARP-1) have significantly improved clinical outcomes in ovarian cancer, especially in patients with BRCA1/2 gene mutations or additional homologous recombination DNA repair pathway deficiencies. Cancers with deficiency of DNA repair (by homologous recombination) are showing greater sensitivity to PARP inhibitors. Beyond BRCA1/2 gene mutations, the loss, perturbation, or overloading of function of various proteins involved in homologous recombination [e.g., ATM, CHEK2, BARD1, BRIP1, MRE11, RAD50, NBS1, RAD51C, RAD51D, RAD52, PALB2, and DNA-dependent protein kinase (DNA-PK)] have been suggested to be synthetically lethal with PARP inhibition. In that case, synthetic lethality could be potentially optimized when PARP inhibitors are combined with DNA-damaging agents, or with molecularly targeted agents capable of impairing DNA repair, such as the case of azasteroid alkylators [32,33].
The response to PARP inhibitors and platinum-based chemotherapy can be predicted by ADP-ribosylation (ADPRylation), which is also catalyzed by PARP enzymes, independently of the BRCA1/2 or homologous recombination deficiency (HRD) status [34]. In general, PARP-1 is activated by direct interaction with pERK2, which then amplifies ERK-signals, resulting in enhanced pERK2-catalyzed phosphorylation of target transcription factors and increased gene expression. Consequently, the pharmacologic inhibition of PARP can cause loss of ERK2 stimulation, thus blocking the ERK signaling network. Moreover, selective ERK inhibition induces tumor regression in MEK inhibitor-resistant models, and PARP inhibition could be similarly effective in MEK or ERK-resistant tumors. In fact, MEK, ERK, and PARP inhibitors, as well as their combinations, have recently entered in early phase clinical trials [35,36]. On the other hand, given that RAS mutant tumors are usually resistant to PARP inhibitors and other anticancer drugs, the inhibition of MEK and ERK can reverse PARP inhibitor resistance in KRAS mutant tumors, including ovarian cancer. Furthermore, the effects of treatment with combinations of PARP and MEK inhibitors are independent of BRCA1/2 and p53 mutation status, suggesting a generalized synergistic activity [37].
In general, the PI3K/AKT signaling pathway plays a critical role in the therapeutic management of ovarian cancer. More specifically, inhibition of PARP1 activity can lead to the upregulation of the PI3K/AKT pathway, thus promoting resistance to PARP inhibitors (PARPi). For instance, pharmacologic inhibition of PARP-1 promotes PI3K/AKT activity and mTOR signaling resulting in decreased cell death as well as PARPi response. Activation of the AKT survival pathway may counteract the cytotoxic effects of PARP inhibition, suggesting that the AKT pathway inhibition may enhance PARP inhibition in the context of antitumor therapy. As a result, combination(s) of small molecules targeting PARP, as well as the PI3K/AKT pathway, would potentially provide a greater clinical benefit [36,38]. On the other hand, the Ras-ERK and PI3K-mTOR signaling pathways intersect in order to regulate each other and to also co-regulate downstream functions. Both pathways can negatively regulate each other’s activity, and cross-inhibition is often emerging. Moreover, when one pathway is chemically blocked, the other pathway can effectively be activated. Also, the Ras-ERK pathway cross-activates PI3K-AKT-mTORC1 signaling as long as Ras-GTP can directly bind and allosterically activate Phosphoinositide 3-kinase (PI3K). Finally, both pathways are often activated by genetic alterations in the upstream signaling molecules, which could further result in increased transcription of various genes promoting cell growth and malignant transformation. So, the rationale of dual inhibition of both of these pathways in suppressing ovarian cancer appears to be promising [39]. In fact, the extent and importance of this crosstalk and its significance in cancer therapeutics is becoming evident since AKT and ERK signaling are both aberrantly activated in a wide range of other human cancers [40]. Nevertheless, the utilization of dual-targeted inhibitors for PI3K/AKT/mTOR and RAS/ERK signaling pathways is a potentially effective therapeutic approach against ovarian cancer [41].
On the other hand, steroid hormones promote epithelial ovarian cancer growth, and their receptor expression is associated with disease outcome. Estrogen (ER-α/-β), and progesterone receptors PR(-A/-B) are frequently expressed in ovarian cancer with an important variability relating to histological subtype, grade, and stage. Hormone therapy is frequently used, but the significance of activity by specific agents or tumor characteristics is unknown. Although the activity of endocrine therapy in advanced ovarian cancer is considerable, limited clinical trials have shown modest response rates to endocrine therapy. However, a small subset of patients that respond very well was identified with few side effects [42,43,44]. The steroid components of ENGA-L06E and ENGA-L08E are testosterone and estrone derivatives, respectively. It has been reported that the steroidal moiety of homo-aza-steroid alkylators is binding to estrogen receptors (ER-α), in breast cancer cells with positive ER expression, and consequently was shown to be more sensitive to their cytostatic and cytotoxic effects [25,45]. However, both ENGA-L06E and ENGA-L08E compounds were shown significant anticancer activity against all human ovarian carcinoma cell lines tested and regardless of the steroid receptor expression.

4. Materials and Methods

4.1. Synthesis of Chemical Compounds

All modified hybrid steroidal esters and comparative alkylators were synthesized with methods that have been previously described [23]. More specifically, we have tested the two previously synthesized hybrid homo-azasteroid alkylators (e.g., ENGA-L06E and ENGA-L08E) along with the proper reference alkylator aniline mustards (e.g., POPAM and L-PAM). The lactam-steroidal carrier of ENGA-L06E is a derivative of testosterone while ENGA-L06E of estrone. The chemical structures of all studied compounds are shown in Figure 11.

4.2. Cell Lines

Five human ovarian cancer cell lines (UWB1.289, UWB1.289+BRCA1, SK-OV-3, OVCAR-3, and OVCAR-5) were treated with the newly synthesized compounds at concentrations of 1–100 μΜ in order to test their cytostatic and cytotoxic activity. All cell lines were obtained from the American Type Culture Collection (ATCC), USA except OVCAR-5, which was kindly donated from BRFFA, Academy of Athens, Greece, and were grown in different culture media (according to the manufacturer’s recommendations), 37 °C and 5% CO2. Finally, all specific molecular characteristics and mutated genes’ profiles, in all five human ovarian cancer cell lines, concerning chemotherapeutic drug tolerance, are shown in Table A8. The cell lines have been tested with standards procedures for microbial contamination, including mycoplasma.

4.3. Determination of Cell Viability by the MTT Assay

Cells were plated in 96-well plate at a density of 1 × 104 cells/mL per well and maintained for 72 h at 37 °C in a 5% CO2 incubator and grown as monolayers. After 24 h, cells were treated with 0.1–100 μmol/L of the tested compounds for 48 h. The viability of cultured cells was estimated by utilizing the (3-(4,5-imethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) metabolic assay, as described previously [23,24]. Briefly, MTT (Sigma, St Louis, Missouri, USA) was dissolved in PBS in a concentration of 5 mg/mL, filter sterilized, and stored at 4 °C. Then, 50 μL of stock solution was added to each culture and incubated for 3 h at 37 °C. Formazan crystals were solubilized by DMSO (100 μL). The absorbance of the converted dye was measured at a wavelength of 540 nm on ELISA reader (VersaMax Microplate Reader Orleans, USA). The mean concentrations of each drug that generated 50% or total (100%) growth inhibition (GI50 and TGI, respectively), as well as the drug concentrations that produced cytotoxicity against 50% of the cultured cells [half maximal inhibitory concentration (IC50)], were calculated using the linear regression method [46]. Using seven absorbance measurements [time 24 h (Ct24), control growth 72 h (Ct72), and test growth in the presence of the drug at five concentration levels (Tt72x)], the percentage of growth was calculated at each level of the drug concentrations. The percentage growth inhibition was calculated according to National Cancer Institute (NCI) as: [(Tt72x) − (Ct24)/(Ct72) − (Ct24)] × 100 for concentrations for which Tt72x > Ct24 and [(Tt72x) − (Ct24)/Ct24] × 100 for concentrations for which Tt72x < Ct24. The GI50 was calculated from [(Tt72x) − (Ct24)/(Ct72) − (Ct24)] × 100 = 50 while the TGI from [(Tt72x) − (Ct24)/(Ct72) − (Ct24)] × 100 = 0, and IC50 from [(Tt72x) − (Ct24)/Ct24] × 100 = 50. All experiments were carried out in triplicate.

4.4. Determination of ERK 1/2 and AKT 1/2/3 Phosphorylation by ELISA

Phosphorylation of ERK 1/2 (Thr202/Tyr204, Thr185/Tyr187) and AKT 1/2/3 (Ser473) was determined by the Activation InstantOne™ ELISA assay kit (Invitrogen, eBioscience, (San Diego, CA, USA) at 450 nm and according to manufacturer’s recommendations. All cell lines were cultured in 12-well plates, seeding 105cells per well. After 24 h, cells were treated with ENGA-L06E and ENGA-L08E, at concentrations of 10, 20, and 60 μΜ, and time exposures of 0, 15, 30, and 60 min. Cells were washed in PBS, extracted in RIPA buffer, and the protein content was quantified by the Bradford assay. Then, 50 μg/mL of protein concentration was added in a 96-well plate, in triplicates, followed by the addition of 50 μL/well of mixed antibodies AKT/ERK 1,2 and pAKT/ERK 1,2 for 1 h. Samples were washed three times, and then 100 μL/well of detection reagent was added and incubated for 30 min. Finally, the reaction was terminated by the addition of 100 μL/well of STOP solution, after which the absorption of all samples was measured at 450 nm.

4.5. Determination of ERK 1/2 and AKT 1/2/3 Phosphorylation by Western Immunoblotting

In order to confirm that ENGA-L06E and ENGA-L08E inhibit the phosphorylation of ERK1/2 and AKT 1/2/3, western immunoblotting was utilized [46]. Cells were treated with the tested compounds for 0, 0.25, 0.5, 1, 3, 6, 12, and 24 h, and samples were prepared in the presence of ice-cold RIPA buffer and protease inhibitor cocktail mixture. The proteinaceous extract was vortexed, centrifuged at 14,000× g, for 10 min, at 4 °C and the total protein content was quantitated by the BCA assay. Then, an equal quantity (25 μg) from each protein sample was loaded into SDS-PAGE electrophoresis followed by transfer to a nitrocellulose membrane with further incubation at 5% skim milk in Tween 20-Tris buffered saline (TBST). After blocking, the membrane was incubated overnight with primary antibodies [ERK1/2, pERK1/2, pAKT, AKT, and Glyceraldehyde 3-phosphate dehydrogenase (GAPDH] at specified concentrations. Finally, the membrane was washed and then incubated with the HRP-conjugated rabbit anti-mouse secondary antibody. Signals were detected using an automated image acquisition system, while all acquired images were further normalized to GAPDH and semi-quantified by the ImageJ software (NIH, Bethesda, Washington D.C., USA).

4.6. Determination of In Vivo Toxicity and Antitumor Efficacy

CB17 severe combined immunodeficiency (SCID) male and female mice were obtained from the National Center for Scientific Research (NCSR) Demokritos, Institute of Biology, Athens, Greece. They were fed water and an irradiated standard rodent diet ad libitum and housed under pathogen-free conditions. All animals were sacrificed by cervical disruption when shown signs of cachexia, suffering, or the tumor volume reached 2000 mm3. For intraperitoneal (i.p.) treatment, stock solutions of the four tested compounds were prepared immediately before use. They were suspended in corn oil, in the desired concentration, following initial dissolution in 10% dimethylsulfoxide (DMSO) without a record of any observable toxic effects. Briefly, acute toxicity induced by the tested compounds was determined after a single i.p. injection (into groups of ten CB17 SCID mice, which were then observed for 30 days) and as described previously [11,47]. The LD10 and LD50 (lethal dose for 10% and 50% of animals, respectively) were determined after graphical estimation (30-day curves) where the percentage of deaths due to the toxicity of each dose was shown in the ordinate while the administered doses were indicated on the abscissa. Two animal models were established based on either SK-OV-3 or OVCAR-5 being implanted subcutaneously (s.c.) into CB17 SCID mice. Overall, at Day 0 (D0), 4 × 106 cells were implanted s.c., and mice were randomized into either test (n = 5) or control (n = 7) groups while their tumor volumes were in the 250–300 mm3 range. After randomization, all animals were treated i.p. with either the tested compounds (at a single LD10 dose; test group) or saline buffer (at a single 0.2 mL dose; control group). The first day that each animal model was treated with each of the tested compounds was determined as Day 1 (d-1). Finally, all stock solutions for each of the tested compounds were prepared immediately before use. The subcutaneously tumors were measured, and the animals were weighed twice a week starting on the first day of treatment. Tumors were measured using an analog caliper. Tumor volume (mm3) was calculated according to the formula D × d2/2, where D and d refer to the larger and smaller perpendicular dimensions collected at each measurement, respectively. The diameters were determined from two orthogonal measurements on each tumor. The mean tumor volume was determined by average, all tumor volumes in each group [48]. On the other hand, for the evaluation of antitumor efficacy, the following parameters were considered: (i) weekly mean tumor area change and tumor inhibition (TI) were determined by utilizing the formula: TI(%) = [1-(TWT − TWZ)/(TWC − TWZ)] × 100 [where TWT is the tumor area (mm3) in treated animals at the time of evaluation; TWZ is the tumor area (mm3) at the time of initiation of treatment (day-1), and TWC is the tumor volume (mm3) in untreated animals (controls) at the time of evaluation]; (ii) percentage of increase in median lifespan [(T/C)%] of treated (T) over control (C) animals. All mice were sacrificed when the tumor volume had reached 2000mm3 (endpoint) [47,49]. Furthermore, three other clinical criteria were also considered as part of our final determination(s): (iii) time of progression-free survival (PFS) which represents the time, in days, that tumor volume presented growth >20% (after treatment with each test compound); (iv) overall survival (OS), which represents the median time, in days [where the treated and untreated animals lived from the day of tumor cells inoculation (D0)]; (v) response rate (RR) which represents the percentage of dose-related response to the tested compounds, as it is interpreted by tumor volume regression (>50%).
The protocol was approved by the ethical committee of the National & Kapodistrian University of Athens and was conducted according to the European Directive 86/609/EEC guidelines for the care and use of laboratory animals. The study received a permit from the Veterinary Directorate of the Prefecture of Athens (Approval #: 2023/2017) according to the Greek legislation conforming to the 2010/53/EU Council Directive.

4.7. Elucidation of Underlined Molecular Mechanism(s) by in Silico Approaches

Crystal structures of AKT1, AKT2, ERK1, and ERK2 were available in the RCSB Protein Data Bank of Brookhaven National Laboratory from which we have chosen the ones with the highest resolution (Å) (the higher the resolution, the lower the Å value). The chosen PDB carries a co-crystalized ligand with a carbon chain quite large to resemble as much as possible with our long chain steroid hybrid molecules in order to fit properly at the proteins’ binding site cavity. These structures were superimposed, showing no conformational variability for the backbone and for most sidechains. Waters were deleted with Maestro, the graphical interface of Schrödinger software (http://www.schrodinger.com), prior to docking. Molecular docking studies were carried out for the best-fitted compounds to the model, while the final selection criteria were compounds docking scores and the presence of crucial interactions for binding to the studied proteins [50]. Finally, the resulting poses were examined manually, and the most promising ones were used again in docking with Induced Fit Docking protocol 2013-2 [51,52], Glide, version 5.9, Prime, version 3.2 (Schrödinger, LLC, New York, NY) [53,54]. To rationally design the steroidal conjugates to be able to bind to the binding sites of ERK1/2 and AKT1/2, we used for the induced fit docking (IFD) process the different available crystal structures shown in Table A9. All of our newly synthesized compounds are sketched and converted into three-dimensional MOL2 files using Schrödinger Release 2015-2, Maestro Version 10.5, and minimized using Schrödinger Release 2015-2, LigPrep 3.4 [54], and the OPLS3 [20] force field. In all cases, the protonation state for the compounds was adjusted at physiological pH using LigPrep (Schrödinger, http://www.schrodinger.com).

4.8. Statistical Analysis

All experiments were repeated in triplicate. Data were presented as mean ± standard deviation (SD). Statistical differences between experimental groups were compared using the Student’s t-test (two-tailed paired). Statistical significance was considered at p < 0.05. Microsoft Excell (Microsoft Hellas, Athens, Greece) was used to analyze the data.

5. Conclusions

Two new homo-aza-steroidal alkylating ester compounds (e.g., ENGA-L06E and ENGA-L08E) were tested in an in vitro model of human ovarian carcinoma (consisting of five cell lines namely UWB1.289, UWB1.289+BRCA1, SK-OV-3, OVCAR-5, and OVCAR-3) as well as an in vivo model consisting of SK-OV-3 and OVCAR-5 tumor xenografts. Both compounds were found to exert very high anticancer activity by inhibiting both the PI3K-AKT and KRAS-ERK signaling pathways. Detailed in silico computational studies were also in accordance and confirmed results of the in vitro studies. As has been previously reported, both of these compounds are capable of inhibiting PARP activity, and consequently, overcome resistance of ovarian cancer cells to alkylating agents, including cisplatin. Moreover, mutations on TP53, KRAS, MEK1 and PIK3CA, microsatellite instability (MSI) status, and steroid receptor expression did not affect the anticancer effect of both compounds. Thus, lactam steroidal alkylators (e.g., ENGA-L06E and ENGA-L08E) provide a novel combination of pharmacological properties exerting significant antitumor activity against human ovarian carcinoma independently of the RAS-ERK and PI3K-AKT mutation, MMR and steroid receptor status, as well as, existing tumor cell resistance to alkylating agents. Conclusively, these data support evidence for further development of lactam steroid alkylators in clinical studies on the treatment of epithelial ovarian carcinoma.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/5/1263/s1, Figure S1: Densitometry readings / intensity ratios of expression and quantification data of Figure 3, Figure S2: Whole western blots, with molecular weight markers (M.W.M.), presenting the expression levels of pERK 1,2 and pAKT in SK-OV-3 and OVCAR-3 ovarian cancer cell lines after treatment with hybrid aza-steroidal alkylators ENGA-L06E and ENGA-L08E at a concentration of 25 μM in 1 and 3 h intervals.

Author Contributions

Conceptualization, P.D. and D.T.T.; methodology, P.D., V.S. and D.T.T.; software, V.S.; formal analysis, P.D., V.S., M.I.P. and D.T.T.; investigation, P.D., E.G., A.P., S.S., N.N., V.S. and D.T.T.; resources, M.I.P. and D.T.T.; data curation, P.D., E.G., A.P., S.S., N.N., V.S. and D.T.T.; writing-original draft preparation, P.D., E.G., A.P., S.S., N.N., V.S., M.I.P. and D.T.T.; writing-review and editing, M.I.P., M.N. and D.T.T.; supervision, D.T.T.; project administration, D.T.T.; funding acquisition, M.I.P. and D.T.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was covered by the international patent W02017001439 and funded by GALENICA S.A. and ENERGON BIOTECH S.A., Greece as well as, in part, by the Hellenic Society of Medical Oncology (HESMO). The APC was funded by The Cyprus Institute of Neurology & Genetics, Nicosia, Cyprus.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The results of inhibitory effects of ENGA-L06E and ENGA-L08E on ERK1,2 phosphorylation, in all studied cancer cell lines, as well as the corresponding kinetics according to Michaelis Menten model (significance level p < 0.01).
Table A1. The results of inhibitory effects of ENGA-L06E and ENGA-L08E on ERK1,2 phosphorylation, in all studied cancer cell lines, as well as the corresponding kinetics according to Michaelis Menten model (significance level p < 0.01).
pERK 1/2 INHIBITION Model Michaelis Menten/Equation: y = Vmax ∗ x/(Km + x)
ENGA-L06EENGA-L08E
OVCAR-5 cells15 min30 min60 min15 min30 min60 min
Reduc. Chi-Sqr3.1261.9500.8673.4357.3543.869
Adj. R-Square0.9950.9980.9990.9940.9900.995
Vmax ± SE68.95 ± 3.2174.99 ± 2.4780.67 ± 1.5263.85 ± 3.6267.94 ± 4.5469.19 ± 2.95
Km ± SE7.58 ± 1.337.06 ± 0.915.64 ± 0.488.97 ± 1.756.09 ± 1.744.17 ± 0.97
IC50 (μM)20.0214.129.1832.3816.9610.86
OVCAR-3 cells
Reduc. Chi-Sqr4.2441.0352.7250.0811.6000.941
Adj. R-Square0.9910.9990.9980.99980.9980.999
Vmax ± SE53.57 ± 3.4180.62 ± 1.9893.91 ± 3.6156.77 ± 0.4961.29 ± 1.8863.79 ± 1.33
Km ± SE5.91 ± 1.648.87 ± 0.7511.21 ± 1.326.61 ± 0.234.01 ± 0.692.63 ± 0.42
IC50 (μM)4914.4712.7648.8517.779.51
SK-OV-3 cells
Reduc. Chi-Sqr2.0740.8950.9061.4490.1260.086
Adj. R-Square0.9980.9990.9990.9960.99980.9999
Vmax ± SE79.14 ± 2.6282.95 ± 1.6984.71 ± 1.6045.92 ± 2.0771.49 ± 0.5574.32 ± 0.43
Km ± SE7.59 ± 0.957.23 ± 0.576.10 ± 0.496.52 ± 1.204.83 ± 0.193.72 ± 0.13
IC50 (μM)13.0110.978.78>5011.257.67
UWB1.289 cells
Reduc. Chi-Sqr2.1190.4043.3601.3220.4042.372
Adj. R-Square0.9980.99970.9980.99870.99970.99844
Vmax ± SE76.21 ± 2.1184.08 ± 0.8889.63 ± 2.5373.97 ± 1.6484.08 ± 0.8888.83 ± 2.13
Km ± SE3.57 ± 0.602.81 ± 0.212.80 ± 0.573.36 ± 0.472.81 ± 0.212.84 ± 0.49
IC50 (μM)6.814.123.537.024.253.67
UWB1.289+BRCA1 cells
Reduc. Chi-Sqr0.4650.4830.1010.1811.9081.626
Adj. R-Square0.9990.9990.99980.99970.9980.998
Vmax ± SE59.43 ± 1.1072.79 ± 1.1474.35 ± 0.4956.97 ± 0.6166.56 ± 1.9971.72 ± 1.86
Km ± SE5.41 ± 0.465.67 ± 0.404.52 ± 0.153.49 ± 0.233.53 ± 0.653.72 ± 0.57
IC50 (μM)28.6812.439.2925.0210.658.56
Table A2. The results of inhibitory effects of ENGA-L06E and ENGA-L08E on AKT1,2 phosphorylation, in all studied cancer cell lines, as well as the corresponding kinetics according to Michaelis Menten model (significance level p < 0.01).
Table A2. The results of inhibitory effects of ENGA-L06E and ENGA-L08E on AKT1,2 phosphorylation, in all studied cancer cell lines, as well as the corresponding kinetics according to Michaelis Menten model (significance level p < 0.01).
pAKT 1/2 INHIBITION Model Michaelis Menten/Equation: y = Vmax ∗ x/(Km + x)
ENGA-L06EENGA-L08E
OVCAR-5 cells15 min30 min60 min15 min30 min60 min
Reduc. Chi-Sqr0.4290.3212.1981.3760.0861.047
Adj. R-Square0.9990.9990.9960.9960.99980.998
Vmax ± SE57.76 ± 3.4959.24 ± 0.4962.49 ± 0.2560.45 ± 3.6561.17 ± 0.6261.47 ± 2.03
Km ± SE12.77 ± 2.127.71 ± 0.237.37 ± 0.1114.15 ± 2.247.64 ± 0.286.48 ± 0.85
IC50 (μM)>5049.8429.51>5034.1928.26
OVCAR-3 cells
Reduc. Chi-Sqr11.2073.44512.9780.2330.2503.633
Adj. R-Square0.9830.9940.9780.9990.99970.996
Vmax ± SE57.71 ± 5.5160.76 ± 2.9762.20 ± 5.1265.97 ± 1.6075.64 ± 0.9978.39 ± 4.18
Km ± SE4.49 ± 2.225.25 ± 1.204.47 ± 1.9115.34 ± 0.946.48 ± 0.348.06 ± 1.51
IC50 (μM)29.1324.4218.31>5012.7514.19
SK-OV-3 cells
Reduc. Chi-Sqr0.0140.10640.0301.8482.8672.637
Adj. R-Square0.99990.99980.99990.9960.9950.996
Vmax ± SE79.23 ± 0.3986.10 ± 1.1387.27 ± 0.5979.22 ± 0.8981.78 ± 2.5382.20 ± 1.68
Km ± SE15.16 ± 0.1916.13 ± 0.5215.92 ± 0.2717.31 ± 0.4615.85 ± 1.2114.03 ± 0.76
IC50 (μM)25.9422.3421.3729.6224.9421.78
UWB1.289 cells
Reduc. Chi-Sqr0.7026.1083.4830.6700.4093.296
Adj. R-Square0.9990.9920.9960.9990.99950.997
Vmax ± SE62.12 ± 1.4164.94 ± 4.2271.12 ± 3.3463.02 ± 1.3368.23 ± 0.9970.28 ± 2.75
Km ± SE4.11 ± 0.504.31 ± 1.454.98 ± 1.103.67 ± 0.453.00 ± 0.292.68 ± 0.77
IC50 (μM)16.9314.4111.7714.078.246.60
UWB1.289+BRCA1 cells
Reduc. Chi-Sqr4.1303.1084.3140.5874.0698.907
Adj. R-Square0.9910.9940.9930.9990.9930.986
Vmax ± SE59.78 ± 4.6757.10 ± 3.3863.12 ± 4.1957.59 ± 1.4260.90 ± 3.8263.84 ± 5.49
Km ± SE8.81 ± 2.306.00 ± 1.496.76 ± 1.755.49 ± 0.605.80 ± 1.555.37 ± 2.07
IC50 (μM)44.4240.2625.7836.2326.6219.42
Table A3. In vivo antitumor effects of compounds ENGA-L06E, ENGA-L08E, L-PAM and POPAM on SKOV-3 human ovarian adenocarcinoma xenograft bearing mice.
Table A3. In vivo antitumor effects of compounds ENGA-L06E, ENGA-L08E, L-PAM and POPAM on SKOV-3 human ovarian adenocarcinoma xenograft bearing mice.
CompoundsTreatment
Schedule
(Days)
Dose
(mg/kg)
MST a
(Days)
TI% bT/C% cCures d
ENGA-L08Ed-175 (i.p.)76.5 ± 3.9
p < 0.01
76.4187.92/5
ENGA-L06Ed-1120 (i.p.)79.6 ± 4.0
p < 0.01
67.6195.62/5
L-PAMd-114 (i.p.)48.5 ± 2.310.5119.70/5
POPAMd-115 (i.p.)47.7 ± 1.512.0117.20/5
Controlsd-1saline40.7 ± 3.0--100.00/7
Tumor cells were inoculated at day 0; a MST: mean survival time (days), animals were sacrificed when the tumor volume had reached 2000mm3 (endpoint); b TI: % Tumor Inhibition; c T/C: mean survival time (%) of drug-treated animals (T) versus saline-treated controls (C); d Cures: Defined as mice alive for 120 days after tumor inoculation (D0).
Table A4. In vivo antitumor effects of compounds ENGA-L06E, ENGA-L08E, L-PAM and POPAM on OVCAR-5 human ovarian adenocarcinoma xenograft bearing mice.
Table A4. In vivo antitumor effects of compounds ENGA-L06E, ENGA-L08E, L-PAM and POPAM on OVCAR-5 human ovarian adenocarcinoma xenograft bearing mice.
CompoundsTreatment
Schedule
(Days)
Dose
(mg/kg)
MST a
(Days)
TI% bT/C% cCures d
ENGA-L08d-175 (i.p.)94.8 ± 6.3
p < 0.01
131.0187.83/5
ENGA-L06Ed-1120 (i.p.)93.2 ± 2.4
p < 0.01
117.0184.80/5
L-PAMd-114 (i.p.)61.0 ± 4.324.5120.80/5
POPAMd-114 (i.p.)59.7 ± 3.36.22118.30/5
Controlsd-1saline50.5 ± 5.8--100.00/7
Tumor cells were inoculated at day 0; a MST: mean survival time (days), animals were sacrificed when the tumor volume had reached 2000 mm3 (endpoint); b TI: % Tumor Inhibition; c T/C: mean survival time (%) of drug-treated animals (T) versus saline-treated controls (C); d Cures: Defined as mice alive for 120 days after tumor inoculation (D0).
Table A5. Evaluation of the in vivo antitumor effect with the clinical criteria: Median Overall Survival (OS), Progression Free Survival (PFS) and Response Rate (RR).
Table A5. Evaluation of the in vivo antitumor effect with the clinical criteria: Median Overall Survival (OS), Progression Free Survival (PFS) and Response Rate (RR).
Clinical Research Criteria
Median Overall Survival
(OS) (Days)
Progression Free Survival
(PFS) (Days)
Response Rate (%) (RR)
DrugsSK-OV-3OVCAR-5SK-OV-3OVCAR-5SK-OV-3OVCAR-5
ENGA-L061129811 27100100
ENGA-L0810512713 34100100
L-PAM49600000
POPAM49570000
Control50
(>2000 mm3)
400000
Table A6. Docking scores and Glide e-model in kcal/ mol of molecules ENGA-L06E and ENGA-L08E compared to the co-crystal ligand 38Z0 in proteins ERK1 and ERK2.
Table A6. Docking scores and Glide e-model in kcal/ mol of molecules ENGA-L06E and ENGA-L08E compared to the co-crystal ligand 38Z0 in proteins ERK1 and ERK2.
ERK1 (pdb ID: 4QTB)ERK2 (pdb ID: 4QTA)
Docking Score
(kcal mol−1)
Glide e-Model
(kcal mol−1)
Docking Score
(kcal mol−1)
Glide e-Model
(kcal mol−1)
Co-crystalized ligand 38Z0−10.372−129.250−8.244−104.505
ENGA-L06E−5.155−60.642−6.848−71.318
ENGA-L08E−7.156−81.061−8.982−79.097
Table A7. Docking scores and Glide e-model in kcal/ mol of molecules ENGA-L06E and ENGA-L08E compared to the co-crystal ligand in proteins AKT-1 and AKT-2.
Table A7. Docking scores and Glide e-model in kcal/ mol of molecules ENGA-L06E and ENGA-L08E compared to the co-crystal ligand in proteins AKT-1 and AKT-2.
AKT1 (pdb ID: 4EJN)AKT2 (pdb ID: 2X39)
Docking Score
(kcal mol−1)
Glide e-Model
(kcal mol−1)
Docking Score
(kcal mol−1)
Glide e-Model
(kcal mol−1)
Co-crystalized ligand −12,188−131,104−10,158−103,940
ENGA-L06E−7.060−84.653−4.273−59.558
ENGA-L08E−9.446−86.525−4.017−57.739
Table A8. Human ovarian cancer cell lines to be used for the in vitro analysis of the compounds biological activity.
Table A8. Human ovarian cancer cell lines to be used for the in vitro analysis of the compounds biological activity.
Cancer TypeHuman Cell Line DesignationMutated OncogenesSpecial Characteristics
Ovarian AdenocarcinomaSK-OV-3 TP53 (p.S90fs)
PIK3CA (p.H1047R)
Tumor Necrosis Factor; Diphtheria Toxin; Cisplatin and Adriamycin not sensitive; Estrogen receptor α (ERα) positive /Progesterone receptor negative MSI-High *
Epithelial Ovarian AdenocarcinomaOVCAR-3MEK1 (p.G159R; p.R160K; p.A158) TP53 (p.R248Q)Androgen/Estrogen/Progesterone receptor positive;
Cisplatin Adriamycin not sensitive; Melphalan resistant; MSS **
Epithelial Ovarian AdenocarcinomaOVCAR-5KRAS (p.G12V)Cisplatin not sensitive; Estrogen receptor α (ERα) negative/GPER positive; MSS **
Ovarian CarcinomaUWB1.289p53 (625delAG)
BRCA1-null
(2594delC)
p53, cytokeratin 7 (CK-7) positive; calretinin positive; Wilms’ tumor protein (WT) positive; Estrogen/Progesterone receptor negative; MSS **
Ovarian CarcinomaUWB1.289+BRCA1p53
(625delAG)
Estrogen/Progesterone receptor negative
* MSI-high: High Microsatellite Instability, ** MSS: Microsatellite stable.
Table A9. Crystallographic structures of ERK1/2 and AKT1/2 obtained from Brookhaven RCSB PDB (Protein Data Bank).
Table A9. Crystallographic structures of ERK1/2 and AKT1/2 obtained from Brookhaven RCSB PDB (Protein Data Bank).
ProteinPDB ID’sYearResolution (Å)Co-Crystallized Inhibitor
ERK14QTB20141.4038Z
ERK24QTA20141.4538Z
AKT14EJN20122.190R4
AKT22X3920101.93X39

References

  1. Armbruster, S.; Coleman, R.L.; Rauh-Hain, J.A. Management and Treatment of Recurrent Epithelial Ovarian Cancer. Hematol. Oncol. Clin. North. Am. 2018, 32, 965–982. [Google Scholar] [CrossRef]
  2. Smolle, E.; Taucher, V.; Pichler, M.; Petru, E.; Lax, S.; Haybaeck, J. Targeting signaling pathways in epithelial ovarian cancer. Int. J. Mol. Sci. 2013, 14, 9536–9555. [Google Scholar] [CrossRef] [Green Version]
  3. Boussios, S.; Zarkavelis, G.; Seraj, E.; Zerdes, I.; Tatsi, K.; Pentheroudakis, G. Non-epithelial Ovarian Cancer: Elucidating Uncommon Gynaecological Malignancies. Anticancer Res. 2016, 36, 5031–5042. [Google Scholar] [CrossRef] [Green Version]
  4. Gasparri, M.L.; Bardhi, E.; Ruscito, I.; Papadia, A.; Farooqi, A.A.; Marchetti, C.; Bogani, G.; Ceccacci, I.; Mueller, M.D.; Benedetti, P.P. PI3K/AKT/mTOR Pathway in Ovarian Cancer Treatment: Are We on the Right Track? Geburtshilfe Frauenheilkd 2017, 77, 1095–1103. [Google Scholar] [CrossRef] [Green Version]
  5. Cheaib, B.; Auguste, A.; Leary, A. The PI3K/Akt/mTOR pathway in ovarian cancer: Therapeutic opportunities and challenges. Chin. J. Cancer 2015, 34, 4–16. [Google Scholar] [CrossRef]
  6. Wang, J.; Zhou, J.Y.; Wu, G.S. ERK-dependent MKP-1-mediated cisplatin resistance in human ovarian cancer cells. Cancer Res. 2007, 67, 11933–11941. [Google Scholar] [CrossRef] [Green Version]
  7. Camoutsis, C.; Trafalis, D.T.P. An overview on the antileukemic potential of D-homo-aza-and respective 17β-acetamido-steroidal alkylating esters. Invest. New Drugs 2003, 21, 47–54. [Google Scholar] [CrossRef]
  8. Camoutsis, C.; Trafalis, D.; Pairas, G.; Papageorgiou, A. On the formation of 4-[N, N-bis (2-chloroethyl) amino] phenyl acetic acid esters of hecogenin and aza-homo-hecogenin and their antileukemic activity. Il Farmaco 2005, 60, 826–829. [Google Scholar] [CrossRef]
  9. Trafalis, D.T.P.; Camoutsis, C.; Papageorgiou, A. Research on the anti-tumour effect of steroid lactam alkylator (NSC-294859) in comparison with conventional chemotherapeutics in malignant melanoma. Melanoma Res. 2005, 15, 273–281. [Google Scholar] [CrossRef]
  10. Trafalis, D.T.; Camoutsis, C.; Dalezis, P.; Papageorgiou, A.; Kontos, M.; Karamanakos, P.; Giannakos, G.; Athanassiou, A.E. Antitumor effect of a- and d- lactam androgen nitrogen mustards on non-small cell lung carcinoma. J. BUON 2004, 9, 275–282. [Google Scholar]
  11. Trafalis, D.T. Hybrid aza-steroid alkylators in the treatment of colon cancer. Cancer Lett. 2006, 243, 202–210. [Google Scholar] [CrossRef]
  12. Trafalis, D.T.; Polonifi, A.; Dalezis, P.; Nikoleousakos, N.; Katsamakas, S.; Sarli, V. Targeting on poly (ADP-ribose) polymerase activity with DNA-damaging hybrid lactam-steroid alkylators in wild-type and BRCA1-mutated ovarian cancer cells. Chem. Biol. Drug Des. 2017, 90, 854–866. [Google Scholar] [CrossRef]
  13. Hamilton, T.C.; Young, R.C.; McKoy, W.M.; Grotzinger, K.R.; Green, J.A.; Chu, E.W.; Whang-Peng, J.; Rogan, A.M.; Green, W.R.; Ozols, R.F. Characterization of a Human Ovarian Carcinoma Cell Line (NIH:OVCAR-3) with Androgen and Estrogen Receptors. Cancer Res. 1983, 43, 5379–5389. [Google Scholar]
  14. Yan, X.D.; Li, M.; Yuan, Y.; Mao, N.; Pan, L.Y. Biological comparison of ovarian cancer resistant cell lines to cisplatin and Taxol by two different administrations. Oncol. Rep. 2007, 17, 1163–1169. [Google Scholar] [CrossRef] [Green Version]
  15. Godwin, A.K.; Meister, A.; O’Dwyer, P.J.; Huang, C.S.; Hamilton, T.C.; Anderson, M.E. High resistance to cisplatin in human ovarian cancer cell lines is associated with marked increase of glutathione synthesis. Proc. Natl. Acad. Sci. USA 1992, 89, 3070–3074. [Google Scholar] [CrossRef] [Green Version]
  16. Surnar, B.; Kolishetti, N.; Basu, U.; Ahmad, A.; Goka, E.; Marples, B.; Kolb, D.; Lippman, M.E.; Dhar, S. Reduction of Cisplatin-Induced Ototoxicity without Compromising Its Antitumor Activity. Biochemistry 2018, 57, 6500–6513. [Google Scholar] [CrossRef]
  17. Boussios, S.; Karathanasi, A.; Cooke, D.; Neille, C.; Sadauskaite, A.; Moschetta, M.; Zakynthinakis-Kyriakou, N.; Pavlidis, N. PARP Inhibitors in Ovarian Cancer: The Route to “Ithaca”. Diagnostics 2019, 9, 55. [Google Scholar] [CrossRef] [Green Version]
  18. Olsson, M.H.M.; Søndergaard, C.R.; Rostkowski, M.; Jensen, J.H. PROPKA3: Consistent treatment of internal and surface residues in empirical pKa predictions. J. Chem. Theory Comput. 2011, 7, 525–537. [Google Scholar] [CrossRef]
  19. Shivakumar, D.; Williams, J.; Wu, Y.; Damm, W.; Shelley, J.; Sherman, W. Prediction of Absolute Solvation Free Energies using Molecular Dynamics Free Energy Perturbation and the OPLS Force Field. J. Chem. Theory Comput. 2010, 6, 1509–1519. [Google Scholar] [CrossRef]
  20. Fosgerau, K.; Hoffmann, T. Peptide therapeutics: Current status and future directions. Drug Discov. Today 2015, 20, 122–128. [Google Scholar] [CrossRef] [Green Version]
  21. Trafalis, D.T.; Tsavdaridis, D.; Camoutsis, C.; Karayiani, V.; Mourelatos, D.; Dalezis, P.; Athanassiou, A.; Pangalis, G.A.; Papageorgiou, A. Preclinical studies on NSC290205 aza-steroid alkylator activity in combination with Adriamycin against lymphoid leukaemia. Br. J. Haematol. 2005, 128, 343–350. [Google Scholar] [CrossRef]
  22. Geromichalos, G.D.; Geromichalou, E.; Camoutsis, C.; Kontos, M.; Dalezis, P.; Papageorgiou, A.; Grivas, A.A.; Tsigris, C.; Trafalis, D.T. In silico/in vitro study of hybrid D-modified steroidal alkylator anticancer activity using uridine phosphorylase as target protein. Anticancer Res. 2011, 31, 831–842. [Google Scholar]
  23. Trafalis, D.; Geromichalou, E.; Dalezis, P.; Nikoleousakos, N.; Sarli, V. Synthesis and evaluation of new steroidal lactam conjugates with aniline mustards as potential antileukemic therapeutics. Steroids 2016, 115, 1–8. [Google Scholar] [CrossRef]
  24. Trafalis, D.; Dalezis, P.; Geromichalou, P.; Sagredou, S.; Sfakidou, E.; Voura, M.; Grammatikopoulou, M.; Gabriel, C.; Sarli, V. Discovery of steroidal lactam conjugates of POPAM-NH2 with potential anticancer activity. Fut. Med. Chem. 2019, 12, 19–35. [Google Scholar] [CrossRef]
  25. Trafalis, D.T.P.; Geromichalos, G.D.; Koukoulitsa, C.; Papageorgiou, A.; Karamanakos, P.; Camoutsis, C. Lactandrate: A D-homo-aza-androsterone alkylator in the treatment of breast cancer. Breast Cancer Res. Treat. 2006, 97, 17–23. [Google Scholar] [CrossRef]
  26. Boussios, S.; Karihtala, P.; Moscetta, M.; Karathanasi, A.; Sadauskaite, A.; Rassy, E.; Pavlidis, N. Combined Strategies with Poly (ADP-Ribose) Polymerase (PARP) Inhibitors for the Treatment of Ovarian Cancer: A Literature Review. Diagnostics 2019, 9, 87. [Google Scholar] [CrossRef] [Green Version]
  27. Ye, Q.; She, Q.B. Integration of AKT and ERK Signaling Pathways in Cancer: Biological and Therapeutic Implications. J. Pharmacol. Clin. Toxicol. 2013, 1, 1009. [Google Scholar]
  28. Ye, Q.; Cai, W.; Zheng, Y.; Evers, B.M.; She, Q.B. ERK and AKT signaling cooperate to translationally regulate survivin expression for metastatic progression of colorectal cancer. Oncogene 2014, 33, 1828–1839. [Google Scholar] [CrossRef] [Green Version]
  29. Hayashi, H.; Tsuchiya, Y.; Nakayama, K.; Satoh, T.; Nishida, E. Down-regulation of the PI3-kinase/Akt pathway by ERK MAP kinase in growth factor signaling. Genes Cells 2008, 13, 941–947. [Google Scholar] [CrossRef]
  30. Altomare, D.; Wang, H.; Skele, K.; Rienzo, A.; Klein-Szanto, A.; Godwin, A.; Testa, J. AKT and mTOR phosphorylation is frequently detected in ovarian cancer and can be targeted to disrupt ovarian tumor cell growth. Oncogene 2004, 23, 5853–5857. [Google Scholar] [CrossRef] [Green Version]
  31. Burotto, M.; Chiou, V.L.; Lee, J.M.; Kohn, E.C. The MAPK pathway across different malignancies: A new perspective. Cancer 2014, 120, 3446–3456. [Google Scholar] [CrossRef] [Green Version]
  32. Boussios, S.; Karihtala, P.; Moschetta, M.; Abson, C.; Karathanasi, A.; Zakynthinakis-Kyriakou, N.; Ryan, J.E.; Sheriff, M.; Rassy, E.; Pavlidis, N. Veliparib in ovarian cancer: A new synthetically lethal therapeutic approach. Invest. New Drugs 2020, 38, 181–193. [Google Scholar] [CrossRef]
  33. Liu, F.W.; Tewari, K.S. New Targeted Agents in Gynecologic Cancers: Synthetic Lethality, Homologous Recombination Deficiency and PARP Inhibitors. Curr. Treat. Options Oncol. 2016, 17, 12. [Google Scholar] [CrossRef] [Green Version]
  34. Conrad, L.B.; Lin, K.Y.; Nandu, T.; Gibson, B.A.; Lea, J.S.; Kraus, W.L. ADP-ribosylation Levels and Patterns Correlate with Gene Expression and Clinical Outcomes in Ovarian Cancers. Mol. Cancer Ther. 2020, 19, 282–291. [Google Scholar] [CrossRef] [Green Version]
  35. Cohen-Armon, M. PARP-1 activation in the ERK signaling pathway. Trends Pharmacol. Sci. 2007, 28, 556–560. [Google Scholar] [CrossRef]
  36. Weaver, A.N.; Yang, S.E. Beyond DNA repair: Additional functions of PARP-1 in cancer. Front. Oncol. 2013, 3, 290. [Google Scholar] [CrossRef] [Green Version]
  37. Sun, C.; Fang, Y.; Yin, J.; Chen, J.; Ju, Z.; Zhang, D.; Chen, X.; Vellano, C.P.; Jeong, K.J.; Ng, P.K.; et al. Rational combination therapy with PARP and MEK inhibitors capitalizes on therapeutic liabilities in RAS mutant cancers. Sci. Transl. Med. 2017, 9, 392. [Google Scholar] [CrossRef] [Green Version]
  38. Bitler, B.G.; Watson, Z.L.; Wheeler, L.J.; Behbakht, K. PARP inhibitors: Clinical utility and possibilities of overcoming resistance. Gynecol. Oncol. 2017, 147, 695–704. [Google Scholar] [CrossRef] [Green Version]
  39. Mendoza, M.C.; Er, E.E.; Blenis, J. The Ras-ERK and PI3K-mTOR Pathways: Cross-talk and Compensation. Trends Biochem. Sci. 2011, 36, 320–328. [Google Scholar] [CrossRef] [Green Version]
  40. Cao, Z.; Liao, Q.; Su, M.; Huang, K.; Jin, J.; Cao, D. AKT and ERK dual inhibitors: The way forward? Cancer Lett. 2019, 459, 30–40. [Google Scholar] [CrossRef]
  41. Sheppard, K.E.; Cullinane, C.; Hannan, K.M.; Wall, M.; Chan, J.; Barber, F.; Foo, J.; Cameron, D.; Neilsen, A.; Ng, P.; et al. Synergistic inhibition of ovarian cancer cell growth by combining selective PI3K/mTOR and RAS/ERK pathway inhibitors. Eur. J. Cancer 2013, 49, 3936–3944. [Google Scholar] [CrossRef]
  42. Lenhard, M.; Tereza, L.; Heublein, S.; Ditsch, N.; Himsl, I.; Mayr, D.; Friese, K.; Jeschke, U. Steroid hormone receptor expression in ovarian cancer: Progesterone receptor B as prognostic marker for patient survival. BMC Cancer 2012, 12, 553. [Google Scholar] [CrossRef] [Green Version]
  43. Paleari, L.; Gandini, S.; Provinciali, N.; Puntoni, M.; Colombo, N.; DeCensi, A. Clinical benefit and risk of death with endocrine therapy in ovarian cancer: A comprehensive review and meta-analysis. Gynecol. Oncol. 2017, 146, 504–513. [Google Scholar] [CrossRef] [Green Version]
  44. Modugno, F.; Laskey, R.; Smith, A.L.; Andersen, C.L.; Haluska, P.; Oesterreich, S. Hormone response in ovarian cancer: Time to reconsider as a clinical target? Endocr. Relat. Cancer 2012, 19, 255–279. [Google Scholar] [CrossRef] [Green Version]
  45. Kontos, M.; Nikolopoulou, M.; Trafalis, D.T.; Geromichalos, G.D.; Koukoulitsa, C.; Camoutsis, J.C.; Bastounis, E.; Karamanakos, P. The effect of an estrone D-lactam steroid ester derivative on breast cancer cells and its predicted binding interactions with the ligand binding domain of estrogen receptor-alpha. Oncol. Res. 2005, 16, 129–142. [Google Scholar] [CrossRef]
  46. Finlay, G.J.; Wilson, W.R.; Baguley, B.C. Comparison of in vitro activity of cytotoxic drugs towards human carcinoma and leukaemia cell lines. Eur. J. Cancer Clin. Oncol. 1986, 22, 655–662. [Google Scholar] [CrossRef]
  47. Goldim, A. Experimental evaluation of antitumor drugs in the USA and USSR and clinical correlations. Natl. Cancer Inst. Monogr. 1980, 55, 1–179. [Google Scholar]
  48. Tomayko, M.M.; Reynolds, C.P. Determination of subcutaneous tumor size in athymic (nude) mice. Cancer Chemother. Pharmacol. 1989, 24, 148–154. [Google Scholar] [CrossRef]
  49. Mahmood, T.; Yang, P.C. Western blot: Technique, theory, and trouble shooting. North. Am. J. Med. Sci. 2012, 4, 429–434. [Google Scholar]
  50. Geromichalos, G.D.; Alifieris, C.E.; Geromichalou, E.G.; Trafalis, D.T. Overview on the current status of virtual high-throughput screening and combinatorial chemistry approaches in multi-target anticancer drug discovery; Part, I. J. BUON 2016, 21, 764–779. [Google Scholar]
  51. Farid, R.; Day, T.; Friesner, R.A.; Pearlstein, R.A. New insights about HERG blockade obtained from protein modeling, potential energy mapping, and docking studies. Bioorg. Med. Chem. 2006, 14, 3160–3173. [Google Scholar] [CrossRef]
  52. Sherman, W.; Day, T.; Jacobson, M.P.; Friesner, R.A.; Farid, R. Novel procedure for modeling ligand/receptor induced fit effects. J. Med. Chem. 2006, 49, 534–553. [Google Scholar] [CrossRef]
  53. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein-Ligand Complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef] [Green Version]
  54. Halgren, T.A.; Murphy, R.B.; Friesner, R.A.; Beard, H.S.; Frye, L.L.; Pollard, W.T.; Banks, J.L. Glide: A New Approach for Rapid, Accurate Docking and Scoring. 2. Enrichment Factors in Database Screening. J. Med. Chem. 2004, 47, 1750–1759. [Google Scholar] [CrossRef]
Figure 1. Michaelis Menten model (y = Vmax ∗ x/(Km + x) inhibition curves of the ERK phosphorylation by ENGA-L06E and ENGA-L08E in SKOV-3, OVCAR-3, OVCAR-5, UWB1.289, UWB1.289+BRCA1 human ovarian cancer cell lines, treated with 10, 20, and 50 μM at different time conditions (15, 30, 60 min) (p < 0.01).
Figure 1. Michaelis Menten model (y = Vmax ∗ x/(Km + x) inhibition curves of the ERK phosphorylation by ENGA-L06E and ENGA-L08E in SKOV-3, OVCAR-3, OVCAR-5, UWB1.289, UWB1.289+BRCA1 human ovarian cancer cell lines, treated with 10, 20, and 50 μM at different time conditions (15, 30, 60 min) (p < 0.01).
Cancers 12 01263 g001
Figure 2. Michaelis Menten model (y = Vmax ∗ x/(Km + x) inhibition curves of AKT phosphorylation by ENGA-L06E and ENGA-L08E in SKOV-3, OVCAR-3, OVCAR-5, UWB1.289, UWB1.289+BRCA1 human ovarian cancer cell lines, treated with 10, 20, and 50 μM at different time conditions (15, 30, 60 min) (p < 0.01).
Figure 2. Michaelis Menten model (y = Vmax ∗ x/(Km + x) inhibition curves of AKT phosphorylation by ENGA-L06E and ENGA-L08E in SKOV-3, OVCAR-3, OVCAR-5, UWB1.289, UWB1.289+BRCA1 human ovarian cancer cell lines, treated with 10, 20, and 50 μM at different time conditions (15, 30, 60 min) (p < 0.01).
Cancers 12 01263 g002
Figure 3. Expression (western blots) and quantification (inhibition curves) levels of ERK1,2 and AKT phosphorylation. Levels of pERK 1,2 and pAKT in SK-OV-3 and OVCAR-3 cell lines after treatment with hybrid aza-steroidal alkylators ENGA-L06E and ENGA-L08E at a concentration of 25 μM in different time intervals. Densitometry readings / intensity ratios of expression and quantification data are presented in Figures S1A,B and S2.
Figure 3. Expression (western blots) and quantification (inhibition curves) levels of ERK1,2 and AKT phosphorylation. Levels of pERK 1,2 and pAKT in SK-OV-3 and OVCAR-3 cell lines after treatment with hybrid aza-steroidal alkylators ENGA-L06E and ENGA-L08E at a concentration of 25 μM in different time intervals. Densitometry readings / intensity ratios of expression and quantification data are presented in Figures S1A,B and S2.
Cancers 12 01263 g003
Figure 4. (a) In-vivo growth inhibition of the SKOV-3 human ovarian adenocarcinoma induced by ENGA-L06E, ENGA-L08E, L-PAM, and POPAM intraperitoneal administration; (b) The tumor volume of the OVCAR-5 adenocarcinoma for all studied compounds on day 9 after drug administration, when there was the highest decrease in tumor volume (induced by ENGA-L08E); (c) Representative tumor images of SKOV-3 human ovarian cancer xenograft removed from mice 14 days after treatment administration (day 1) of ENGA-L06, ENGA-L08, L-PAM, and POPAM drugs.
Figure 4. (a) In-vivo growth inhibition of the SKOV-3 human ovarian adenocarcinoma induced by ENGA-L06E, ENGA-L08E, L-PAM, and POPAM intraperitoneal administration; (b) The tumor volume of the OVCAR-5 adenocarcinoma for all studied compounds on day 9 after drug administration, when there was the highest decrease in tumor volume (induced by ENGA-L08E); (c) Representative tumor images of SKOV-3 human ovarian cancer xenograft removed from mice 14 days after treatment administration (day 1) of ENGA-L06, ENGA-L08, L-PAM, and POPAM drugs.
Cancers 12 01263 g004
Figure 5. (a) In-vivo growth inhibition of OVCAR-5 human ovarian adenocarcinoma induced by ENGA-L06E, ENGA-L08E, L-PAM, and POPAM intraperitoneal administration; (b) The tumor volume of the OVCAR-5 adenocarcinoma for all studied compounds, on day 16 after drug administration, when there was the highest decrease in tumor volume (induced by ENGA-L08E); (c) Representative tumor images of OVCAR-5 human ovarian cancer xenograft removed from mice 28 days after treatment administration (day 1) of ENGA-L06, ENGA-L08, L-PAM, and POPAM drugs.
Figure 5. (a) In-vivo growth inhibition of OVCAR-5 human ovarian adenocarcinoma induced by ENGA-L06E, ENGA-L08E, L-PAM, and POPAM intraperitoneal administration; (b) The tumor volume of the OVCAR-5 adenocarcinoma for all studied compounds, on day 16 after drug administration, when there was the highest decrease in tumor volume (induced by ENGA-L08E); (c) Representative tumor images of OVCAR-5 human ovarian cancer xenograft removed from mice 28 days after treatment administration (day 1) of ENGA-L06, ENGA-L08, L-PAM, and POPAM drugs.
Cancers 12 01263 g005
Figure 6. (A) Alignment of ERK1 (PDB ID 4QTB) (ribbons with silver) and ERK2 (PDB ID 4QTA) (ribbons with gold) with co-crystallized ligands 38Z0 (green color for ERK1 and orange color for ERK2) and (B) Solid surface representation of the crystal structure of ERK1. The crystallographic co-crystallized 38Z0 rendered in green C atoms for ERK1 and the co-crystallized 38Z0 rendered in orange C atoms for ERK2.
Figure 6. (A) Alignment of ERK1 (PDB ID 4QTB) (ribbons with silver) and ERK2 (PDB ID 4QTA) (ribbons with gold) with co-crystallized ligands 38Z0 (green color for ERK1 and orange color for ERK2) and (B) Solid surface representation of the crystal structure of ERK1. The crystallographic co-crystallized 38Z0 rendered in green C atoms for ERK1 and the co-crystallized 38Z0 rendered in orange C atoms for ERK2.
Cancers 12 01263 g006
Figure 7. (A) The active site cavity of ERK and (B) Superimposition of calculated molecules ENGA-L06E and ENGA-L08E given in purple and pink colored carbons, respectively, with labeling of the critical bonding residues.
Figure 7. (A) The active site cavity of ERK and (B) Superimposition of calculated molecules ENGA-L06E and ENGA-L08E given in purple and pink colored carbons, respectively, with labeling of the critical bonding residues.
Cancers 12 01263 g007
Figure 8. (A) Alignment of AKT1 with PDB ID 4EJN (ribbons with dark green) and AKT2 with PDB ID 2X39 (ribbons with yellow) with the co-crystallized inhibitor 0R4 (with green color) and X39 (with yellow color), respectively. (B) Surface of AKT1 and the binding site (grey) of the co-crystallized inhibitor 0R4.
Figure 8. (A) Alignment of AKT1 with PDB ID 4EJN (ribbons with dark green) and AKT2 with PDB ID 2X39 (ribbons with yellow) with the co-crystallized inhibitor 0R4 (with green color) and X39 (with yellow color), respectively. (B) Surface of AKT1 and the binding site (grey) of the co-crystallized inhibitor 0R4.
Cancers 12 01263 g008
Figure 9. (A) The active site cavity of AKT1 and (B) Superimposition of the calculated molecules ENGA-L06E and ENGA-L08E, given in purple and pink colored carbons, respectively, with the co-crystallized ligand of AKT1, 0R4 (with green color) and (C) The critical network bonding residues are labeled in yellow and purple dashed lines.
Figure 9. (A) The active site cavity of AKT1 and (B) Superimposition of the calculated molecules ENGA-L06E and ENGA-L08E, given in purple and pink colored carbons, respectively, with the co-crystallized ligand of AKT1, 0R4 (with green color) and (C) The critical network bonding residues are labeled in yellow and purple dashed lines.
Cancers 12 01263 g009
Figure 10. (A) Front side and (B) Backside of the active site cavity of AKT1 (PDB ID 4EJN) with hybrid molecules ENGA-L06E and ENGA-L08E docked and superimposition of calculated molecules ENGA-L06E and ENGA-L08E given in purple and pink colored carbons, respectively.
Figure 10. (A) Front side and (B) Backside of the active site cavity of AKT1 (PDB ID 4EJN) with hybrid molecules ENGA-L06E and ENGA-L08E docked and superimposition of calculated molecules ENGA-L06E and ENGA-L08E given in purple and pink colored carbons, respectively.
Cancers 12 01263 g010
Figure 11. Chemical structures of 2-amino-3-(4-(bis (2-chloroethyl) amino) phenyl) propanoic acid (L-phenylalanine mustard, L-PAM), 3-(4-(bis(2-chloroethyl) amino) phenoxy) propanoic acid (POPAM), 5a,7a-dimethyl-2-oxo-2,3,4,5,5a,5b,6,7,7a,8,9,10,10a,10b,11,12-hexa decahydrocyclopenta [5,6] naphtho [1,2-d] azepin-8-yl 3-(4-(bis(2-chloroethyl)amino) phenoxy) propanoate (ENGA-L06E) and 12a-methyl-2-oxo-1,2,3,4,4a,4b,5,6,10b,11,12,12a-dodecahydronaphtho[2,1-f]quinolin-8-yl 3-(4-(bis (2-chloroethyl)amino) phenoxy) propanoate (ENGA-L08E).
Figure 11. Chemical structures of 2-amino-3-(4-(bis (2-chloroethyl) amino) phenyl) propanoic acid (L-phenylalanine mustard, L-PAM), 3-(4-(bis(2-chloroethyl) amino) phenoxy) propanoic acid (POPAM), 5a,7a-dimethyl-2-oxo-2,3,4,5,5a,5b,6,7,7a,8,9,10,10a,10b,11,12-hexa decahydrocyclopenta [5,6] naphtho [1,2-d] azepin-8-yl 3-(4-(bis(2-chloroethyl)amino) phenoxy) propanoate (ENGA-L06E) and 12a-methyl-2-oxo-1,2,3,4,4a,4b,5,6,10b,11,12,12a-dodecahydronaphtho[2,1-f]quinolin-8-yl 3-(4-(bis (2-chloroethyl)amino) phenoxy) propanoate (ENGA-L08E).
Cancers 12 01263 g011
Table 1. In vitro growth inhibition/cytostatic [50% Growth Inhibition (GI50) and Total Growth Inhibition (TGI) μΜ] and cytotoxic effects (IC50 μM) of compounds POPAM, ENGA-L06E, and ENGA-L08E on OVCAR-3, OVCAR-5, SK-OV-3, UBW 1.289, and UWB1.289+BRCA1 cancer cell lines.
Table 1. In vitro growth inhibition/cytostatic [50% Growth Inhibition (GI50) and Total Growth Inhibition (TGI) μΜ] and cytotoxic effects (IC50 μM) of compounds POPAM, ENGA-L06E, and ENGA-L08E on OVCAR-3, OVCAR-5, SK-OV-3, UBW 1.289, and UWB1.289+BRCA1 cancer cell lines.
POPAMENGA-L06EENGA-L08E
GI50TGIIC50GI50TGIIC50GI50TGIIC50
OVCAR-356 ± 2.792 ± 7.5>10022 ± 0.932 ± 1.843 ± 2.930 ± 2.140 ± 2.756 ± 3.2
OVCAR-552 ± 2.387 ± 7.8>10025 ± 1.442 ± 2.865 ± 3.411 ± 0.660 ± 2.792 ± 8.4
SK-OV-350 ± 2.9>100>10018 ± 0.952 ± 3.376 ± 7.510 ± 0.576 ± 7.2>100
UWB1.28922 ± 1.986 ± 7.4>1003.5 ± 0.48 ± 0.820 ± 1.18 ± 0.920 ± 1.545 ± 3.3
UWB1.289+BRCA1>100>100>10012 ± 0.728 ± 1.940 ± 3.518 ± 1.076 ± 7.288 ± 8.1
Table 2. Acute toxicity of the compounds L-PAM, POPAM, ENGA-L06E, ENGA-L08E in CB17 SCID mice, where LD50 and LD10 are the lethal doses for 50% and 10% of the treated mice population.
Table 2. Acute toxicity of the compounds L-PAM, POPAM, ENGA-L06E, ENGA-L08E in CB17 SCID mice, where LD50 and LD10 are the lethal doses for 50% and 10% of the treated mice population.
CompoundsLD50 (mg/kg)LD10 (mg/kg)
L-PAM20.814
POPAM2015
ENGA-L06E150120
ENGA-L08E9575

Share and Cite

MDPI and ACS Style

Dalezis, P.; Geromichalou, E.; Polonifi, A.; Sagredou, S.; Nikoleousakos, N.; Nikolaou, M.; Sarli, V.; Panayiotidis, M.I.; Trafalis, D.T. Azasteroid Alkylators as Dual Inhibitors of AKT and ERK Signaling for the Treatment of Ovarian Carcinoma. Cancers 2020, 12, 1263. https://doi.org/10.3390/cancers12051263

AMA Style

Dalezis P, Geromichalou E, Polonifi A, Sagredou S, Nikoleousakos N, Nikolaou M, Sarli V, Panayiotidis MI, Trafalis DT. Azasteroid Alkylators as Dual Inhibitors of AKT and ERK Signaling for the Treatment of Ovarian Carcinoma. Cancers. 2020; 12(5):1263. https://doi.org/10.3390/cancers12051263

Chicago/Turabian Style

Dalezis, Panagiotis, Eleni Geromichalou, Aikaterini Polonifi, Sofia Sagredou, Nikolaos Nikoleousakos, Michael Nikolaou, Vasiliki Sarli, Mihalis I. Panayiotidis, and Dimitrios T. Trafalis. 2020. "Azasteroid Alkylators as Dual Inhibitors of AKT and ERK Signaling for the Treatment of Ovarian Carcinoma" Cancers 12, no. 5: 1263. https://doi.org/10.3390/cancers12051263

APA Style

Dalezis, P., Geromichalou, E., Polonifi, A., Sagredou, S., Nikoleousakos, N., Nikolaou, M., Sarli, V., Panayiotidis, M. I., & Trafalis, D. T. (2020). Azasteroid Alkylators as Dual Inhibitors of AKT and ERK Signaling for the Treatment of Ovarian Carcinoma. Cancers, 12(5), 1263. https://doi.org/10.3390/cancers12051263

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop