Synthesis, Anticancer Activity, and Molecular Modeling of New Halogenated Spiro[pyrrolidine-thiazolo-oxindoles] Derivatives

: A one-pot, single-step, and an atom-economical process towards the synthesis of highly functionalized spirooxindoles analogues was e ﬃ ciently conducted to produce a satisfactory chemical yields (70–93%) with excellent relative diastereo-, and regio-selectivity. An in vitro antiproliferative assay was carried out on di ﬀ erent cancer cell lines to evaluate the biological activity of the synthesized tetrahydro-1’ H -spiro[indoline-3,5’-pyrrolo[1,2- c ]thiazol]-2-one 5a–n. The prepared hybrids were then tested in vitro for their antiproliferative e ﬀ ects against three cancer cell lines, namely, HepG2 (liver cancer), MCF-7 (breast cancer), and HCT-116 (colon cancer). The spirooxindole analogue 5g exhibited a broad activity against HepG2, MCF-7, and HCT-116 cell lines of liver, breast, and colorectal cancers when compared to cis platin. Modeling studies including shape similarity, lipophilicity scores, and physicochemical parameters were calculated. The results of this study indicated that spirooxindole analogue 5g retained a good physiochemical parameters with acceptable lipophilicity scores.


Introduction
The design of highly complex spiro-heterocycles with multifunctional and potential pharmaceutical efficacy has attracted considerable attention from synthetic and medicinal chemists [1]. One of the most privileged aza-heterocyclic scaffolds is spiro[pyrrolidine-oxindole] [2], which is present in natural products and useful as a building block for the synthesis of significant biologically active compounds. This class of aza-heterocyclic compounds has gained great interest, owing to several reports of its pharmaceutical potency, including anticancer [3], antitumor [4], 5-HT3 receptor antagonist [5], acetylcholinesterase-inhibitory [6], antibacterial [7], antibiotic [8], and MDM2-p53 inhibitor [9] effects; selective cyclooxygenase COX-1 with TNF-α and IL-6 inhibitors [10]; and potential hypoglycemic dual inhibitory activity against α-amylase and α-glucosidase [11] (Figure 1). To date, prolonged efforts have been exerted to expand divergent complexity and to develop efficient synthetic routes for these valuable privileged aza-heterocyclic scaffolds, which would remarkably enhance their bioactivity [1,12]. In particular, [3+2] cycloaddition is one of the most efficient synthetic approaches to produce these valuable scaffolds with stereoselective method and high yield [13]. To extend our previous research, we explored the effect of halogen substitution on the isatin ring.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 2 of 10 antagonist [5], acetylcholinesterase-inhibitory [6], antibacterial [7], antibiotic [8], and MDM2-p53 inhibitor [9] effects; selective cyclooxygenase COX-1 with TNF-α and IL-6 inhibitors [10]; and potential hypoglycemic dual inhibitory activity against α-amylase and α-glucosidase [11] (Figure 1). To date, prolonged efforts have been exerted to expand divergent complexity and to develop efficient synthetic routes for these valuable privileged aza-heterocyclic scaffolds, which would remarkably enhance their bioactivity [1,12]. In particular, [3+2] cycloaddition is one of the most efficient synthetic approaches to produce these valuable scaffolds with stereoselective method and high yield [13]. To extend our previous research, we explored the effect of halogen substitution on the isatin ring. Our previous studies [9] revealed that the presence of dihalide substitution on acyl moiety substantially increased the anticancer activity of the resulting product(s). Moreover, it was reported [9] that chlorinated indole moiety retained better activity, as illustrated in Figure 1. (a). Subsequently, this study was designed to introduce two bromo atoms on the indole ring, presumably to enhance Our previous studies [9] revealed that the presence of dihalide substitution on acyl moiety substantially increased the anticancer activity of the resulting product(s). Moreover, it was reported [9] that chlorinated indole moiety retained better activity, as illustrated in Figure 1a.
According to Scheme 2 and based on our previous study [9], the reaction proceeds via one pot reaction, in which initially 5,7-dibromoisatin 2 reacted with thioproline 3 affording the azomethine ylide after the removal of carbon dioxide from the intermediate. Subsequently, the azomethine ylide reacted with the bis-benzylidine 1a-n to provide the target compounds in a regioselective and diastereoselective manner. The reaction proceeded via path A regio-selectively to afford the regioisomer products 5a-n, while the second regio-isomers 5a-n' did not occur (path B). There are possible diastereoselective products that could be formed, but in this case only diastereoselective compounds 5a-n occurred via the path C not D.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 4 of 10 According to Scheme 2 and based on our previous study [9], the reaction proceeds via one pot reaction, in which initially 5,7-dibromoisatin 2 reacted with thioproline 3 affording the azomethine ylide after the removal of carbon dioxide from the intermediate. Subsequently, the azomethine ylide reacted with the bis-benzylidine 1a-n to provide the target compounds in a regioselective and diastereoselective manner. The reaction proceeded via path A regio-selectively to afford the regioisomer products 5a-n, while the second regio-isomers 5a-n' did not occur (path B). There are possible diastereoselective products that could be formed, but in this case only diastereoselective compounds 5a-n occurred via the path C not D.

Scheme 2.
Plausible reaction mechanism of the synthesized compounds.

Biological Activity
The compounds were subjected to an initial evaluation for potential cytotoxic activity against different cancer cell lines, namely, HepG2, MCF-7, and HCT-116 cells, at 50 µM. Cell viability was measured using MTT assay. Among the screened 14 compounds, 3 (5h, 5i, and 5j) did not show any cytotoxic activity against HepG2 cells. The concentration of the active compounds that killed 50% of the cells (IC50) was evaluated against HepG2 cells. Compound 5g (IC50 = 5.00 ± 0.66 µM) was the most potent active compound, showing more potent activity than that of the standard chemotherapeutic drug cisplatin (IC50 = 9.00 ± 0.76 µM) ( Table 1). Moderate anticancer activity against HepG2 cells was observed for compounds 5a and 5m (IC50 = 10.00 ± 0.47 and 17.00 ± 0.68 µM, respectively).

Biological Activity
The compounds were subjected to an initial evaluation for potential cytotoxic activity against different cancer cell lines, namely, HepG2, MCF-7, and HCT-116 cells, at 50 µM. Cell viability was measured using MTT assay. Among the screened 14 compounds, 3 (5h, 5i, and 5j) did not show any cytotoxic activity against HepG2 cells. The concentration of the active compounds that killed 50% of the cells (IC 50 ) was evaluated against HepG2 cells. Compound 5g (IC 50 = 5.00 ± 0.66 µM) was the most potent active compound, showing more potent activity than that of the standard chemotherapeutic drug cisplatin (IC 50 = 9.00 ± 0.76 µM) ( Table 1). Moderate anticancer activity against HepG2 cells was observed for compounds 5a and 5m (IC 50 = 10.00 ± 0.47 and 17.00 ± 0.68 µM, respectively).

Effect of the Dibromo on the Anticancer Activity
The structure-activity relationship between the previously reported spirooxindole analogues 4b, 4c, 4d, 4f, and 4i-n [9] and the diboromo-substituted spiroxindoles 5b, 5c, 5d, 5f, and 5i-n is described. In fact, the IC50 values of Table 1 clearly show that the replacement of the H atoms of the previously reported compounds 4b, 4c, 4d, 4f, and 4i-n with that of its analogues with the Br produced a significant decrease in the inhibitory growth effect on the HEPG2 cell line. On the other hand, compounds 5b, 5c, 5f, 5k, and 5m (dibromo-substituted) showed better activity against HCT-116 cells than their dibromo-unsubstituted indole counterparts. Compounds 5d, 5l, and 5n showed less activity than the compounds 4d, 4l, and 4n, respectively. Compounds 5i and 5j were not active and compounds 4i and 4j presented some activity (Table 1).

Shape Alignment by Rapid Overlay Chemical Structure (ROCS) Analysis
Shape and electrostatic potential are two fundamental molecular descriptors for computational drug discovery, because in protein ligand binding, the shape of a ligand has to conform in large degree to the shape of a protein binding site. The electrostatic potentials presented in the binding site have to complement the electrostatic potential of the protein. Accordingly, it is very important to model and understand protein ligand bindings correctly. The 3D shape structure exhibits good neighborhood behavior, in which high similarity in shape reflects high similarity in biology. Shape similarity can have different applications, such as virtual screening, lead-hopping, molecular alignment, pose generation, and predictions.
ROCS is a tool used in shape similarity studies. ROCS requires a query, which is an active molecule with some biological activities in at least one 3D conformation. It also requires a database of the molecules of the compounds of interest. Consistent with these standards, our compounds (database set) showed similarity to standard compound BI-0225 ( Figure 2). Compound 5g showed high similarity to BI-0225 in terms of its oxindole moiety and oxoindole ring.

Effect of the Dibromo on the Anticancer Activity
The structure-activity relationship between the previously reported spirooxindole analogues 4b, 4c, 4d, 4f, and 4i-n [9] and the diboromo-substituted spiroxindoles 5b, 5c, 5d, 5f, and 5i-n is described. In fact, the IC50 values of Table 1 clearly show that the replacement of the H atoms of the previously reported compounds 4b, 4c, 4d, 4f, and 4i-n with that of its analogues with the Br produced a significant decrease in the inhibitory growth effect on the HEPG2 cell line. On the other hand, compounds 5b, 5c, 5f, 5k, and 5m (dibromo-substituted) showed better activity against HCT-116 cells than their dibromo-unsubstituted indole counterparts. Compounds 5d, 5l, and 5n showed less activity than the compounds 4d, 4l, and 4n, respectively. Compounds 5i and 5j were not active and compounds 4i and 4j presented some activity (Table 1).

Shape Alignment by Rapid Overlay Chemical Structure (ROCS) Analysis
Shape and electrostatic potential are two fundamental molecular descriptors for computational drug discovery, because in protein ligand binding, the shape of a ligand has to conform in large degree to the shape of a protein binding site. The electrostatic potentials presented in the binding site have to complement the electrostatic potential of the protein. Accordingly, it is very important to model and understand protein ligand bindings correctly. The 3D shape structure exhibits good neighborhood behavior, in which high similarity in shape reflects high similarity in biology. Shape similarity can have different applications, such as virtual screening, lead-hopping, molecular alignment, pose generation, and predictions.
ROCS is a tool used in shape similarity studies. ROCS requires a query, which is an active molecule with some biological activities in at least one 3D conformation. It also requires a database of the molecules of the compounds of interest. Consistent with these standards, our compounds (database set) showed similarity to standard compound BI-0225 ( Figure 2). Compound 5g showed high similarity to BI-0225 in terms of its oxindole moiety and oxoindole ring.

Effect of the Dibromo on the Anticancer Activity
The structure-activity relationship between the previously reported spirooxindole analogues 4b, 4c, 4d, 4f, and 4i-n [9] and the diboromo-substituted spiroxindoles 5b, 5c, 5d, 5f, and 5i-n is described. In fact, the IC50 values of Table 1 clearly show that the replacement of the H atoms of the previously reported compounds 4b, 4c, 4d, 4f, and 4i-n with that of its analogues with the Br produced a significant decrease in the inhibitory growth effect on the HEPG2 cell line. On the other hand, compounds 5b, 5c, 5f, 5k, and 5m (dibromo-substituted) showed better activity against HCT-116 cells than their dibromo-unsubstituted indole counterparts. Compounds 5d, 5l, and 5n showed less activity than the compounds 4d, 4l, and 4n, respectively. Compounds 5i and 5j were not active and compounds 4i and 4j presented some activity (Table 1).

Shape Alignment by Rapid Overlay Chemical Structure (ROCS) Analysis
Shape and electrostatic potential are two fundamental molecular descriptors for computational drug discovery, because in protein ligand binding, the shape of a ligand has to conform in large degree to the shape of a protein binding site. The electrostatic potentials presented in the binding site have to complement the electrostatic potential of the protein. Accordingly, it is very important to model and understand protein ligand bindings correctly. The 3D shape structure exhibits good neighborhood behavior, in which high similarity in shape reflects high similarity in biology. Shape similarity can have different applications, such as virtual screening, lead-hopping, molecular alignment, pose generation, and predictions.
ROCS is a tool used in shape similarity studies. ROCS requires a query, which is an active molecule with some biological activities in at least one 3D conformation. It also requires a database of the molecules of the compounds of interest. Consistent with these standards, our compounds (database set) showed similarity to standard compound BI-0225 ( Figure 2). Compound 5g showed high similarity to BI-0225 in terms of its oxindole moiety and oxoindole ring.

Effect of the Dibromo on the Anticancer Activity
The structure-activity relationship between the previously reported spirooxindole analogues 4b, 4c, 4d, 4f, and 4i-n [9] and the diboromo-substituted spiroxindoles 5b, 5c, 5d, 5f, and 5i-n is described. In fact, the IC50 values of Table 1 clearly show that the replacement of the H atoms of the previously reported compounds 4b, 4c, 4d, 4f, and 4i-n with that of its analogues with the Br produced a significant decrease in the inhibitory growth effect on the HEPG2 cell line. On the other hand, compounds 5b, 5c, 5f, 5k, and 5m (dibromo-substituted) showed better activity against HCT-116 cells than their dibromo-unsubstituted indole counterparts. Compounds 5d, 5l, and 5n showed less activity than the compounds 4d, 4l, and 4n, respectively. Compounds 5i and 5j were not active and compounds 4i and 4j presented some activity (Table 1).

Shape Alignment by Rapid Overlay Chemical Structure (ROCS) Analysis
Shape and electrostatic potential are two fundamental molecular descriptors for computational drug discovery, because in protein ligand binding, the shape of a ligand has to conform in large degree to the shape of a protein binding site. The electrostatic potentials presented in the binding site have to complement the electrostatic potential of the protein. Accordingly, it is very important to model and understand protein ligand bindings correctly. The 3D shape structure exhibits good neighborhood behavior, in which high similarity in shape reflects high similarity in biology. Shape similarity can have different applications, such as virtual screening, lead-hopping, molecular alignment, pose generation, and predictions.
ROCS is a tool used in shape similarity studies. ROCS requires a query, which is an active molecule with some biological activities in at least one 3D conformation. It also requires a database of the molecules of the compounds of interest. Consistent with these standards, our compounds (database set) showed similarity to standard compound BI-0225 ( Figure 2). Compound 5g showed high similarity to BI-0225 in terms of its oxindole moiety and oxoindole ring.

Effect of the Dibromo on the Anticancer Activity
The structure-activity relationship between the previously reported spirooxindole analogues 4b, 4c, 4d, 4f, and 4i-n [9] and the diboromo-substituted spiroxindoles 5b, 5c, 5d, 5f, and 5i-n is described. In fact, the IC 50 values of Table 1 clearly show that the replacement of the H atoms of the previously reported compounds 4b, 4c, 4d, 4f, and 4i-n with that of its analogues with the Br produced a significant decrease in the inhibitory growth effect on the HEPG2 cell line. On the other hand, compounds 5b, 5c, 5f, 5k, and 5m (dibromo-substituted) showed better activity against HCT-116 cells than their dibromo-unsubstituted indole counterparts. Compounds 5d, 5l, and 5n showed less activity than the compounds 4d, 4l, and 4n, respectively. Compounds 5i and 5j were not active and compounds 4i and 4j presented some activity (Table 1).

Shape Alignment by Rapid Overlay Chemical Structure (ROCS) Analysis
Shape and electrostatic potential are two fundamental molecular descriptors for computational drug discovery, because in protein ligand binding, the shape of a ligand has to conform in large degree to the shape of a protein binding site. The electrostatic potentials presented in the binding site have to complement the electrostatic potential of the protein. Accordingly, it is very important to model and understand protein ligand bindings correctly. The 3D shape structure exhibits good neighborhood behavior, in which high similarity in shape reflects high similarity in biology. Shape similarity can have different applications, such as virtual screening, lead-hopping, molecular alignment, pose generation, and predictions.
ROCS is a tool used in shape similarity studies. ROCS requires a query, which is an active molecule with some biological activities in at least one 3D conformation. It also requires a database of the molecules of the compounds of interest. Consistent with these standards, our compounds (database set) showed similarity to standard compound BI-0225 (Figure 2). Compound 5g showed high similarity to BI-0225 in terms of its oxindole moiety and oxoindole ring.

Predicted Pharmacokinetics and Pharmacodynamics Parameters
Absorption, distribution, metabolism, excretion, and toxicity (ADMET) prediction for drug candidates is mandatory in the drug design process, as these parameters contribute to determining the failure of approximately 60% of all drugs in the development and approval phases. It is wellknown that ADMET prediction is performed at the last stage of the drug development process with high cost and effort. At present, ADMET is determined at the beginning of drug discovery stages in order to eliminate molecules with poor ADMET properties from the drug discovery pipeline with the aim to save research costs. In this regard, computational tools were used to predict ADMET properties in this study [14].
The Caco-2 cell, percentage of human intestinal absorption (HIA), and skin permeability models have all been suggested as reliable in vitro models to estimate oral drug absorption and transdermal delivery [15]. Drug penetration to the blood brain barrier (BBB) provides insight into drugs that act on the central nervous system and on plasma protein binding (PPB). Compared to the other compounds, 5g showed the lowest BBB penetration value (0.017) and a low value in the Caco-2 cell model (18.80). All compounds showed high PPB and HIA values, as well as very low skin permeability values in the range of −1.80 to −2.79 (Table 2).

Ligand Efficiency (LE) and Lipophilic Efficiency (LipE)
In the current study, for optimization assessment, LE was calculated [18]. The parameter LE has a crucial role in "lead optimization for drug-like candidate" properties [19]. Compounds with the

Predicted Pharmacokinetics and Pharmacodynamics Parameters
Absorption, distribution, metabolism, excretion, and toxicity (ADMET) prediction for drug candidates is mandatory in the drug design process, as these parameters contribute to determining the failure of approximately 60% of all drugs in the development and approval phases. It is well-known that ADMET prediction is performed at the last stage of the drug development process with high cost and effort. At present, ADMET is determined at the beginning of drug discovery stages in order to eliminate molecules with poor ADMET properties from the drug discovery pipeline with the aim to save research costs. In this regard, computational tools were used to predict ADMET properties in this study [14].
The Caco-2 cell, percentage of human intestinal absorption (HIA), and skin permeability models have all been suggested as reliable in vitro models to estimate oral drug absorption and transdermal delivery [15]. Drug penetration to the blood brain barrier (BBB) provides insight into drugs that act on the central nervous system and on plasma protein binding (PPB). Compared to the other compounds, 5g showed the lowest BBB penetration value (0.017) and a low value in the Caco-2 cell model (18.80). All compounds showed high PPB and HIA values, as well as very low skin permeability values in the range of −1.80 to −2.79 (Table 2).

Ligand Efficiency (LE) and Lipophilic Efficiency (LipE)
In the current study, for optimization assessment, LE was calculated [18]. The parameter LE has a crucial role in "lead optimization for drug-like candidate" properties [19]. Compounds with the highest activity were selected for evaluation against sensitive cancer cell lines (breast and colon cancer cells).
The compounds had an LE value in the range of 0.19-0.26 except for compound 5n (Table 3). All compounds exhibited higher LE values in breast cancer cells than in colon cancer cells, especially compounds 5c, 5e, and 5l (LE = 0.26), all of which were structural isomers.
The recommended LE value should be in the range of 0.3. The acceptable LE value should be higher than 0.3.

Lipophilic Efficiency (LipE) or Ligand Lipophilic Efficiency (LEE)
Lip E or LLE is an avenue to determine compound affinity with respect to its lipophilicity. Nowadays, the lipophilic efficiency (LipE) index (LEE), which includes lipophilicity and potency, is becoming more and more popular in drug design. It allows for the normalization of observed potency with changes in the lipophilicity, and it is considered an effective and practical tool for keeping lipophilicity under control to avoid any "molecular obesity".
LipE or LLE is calculated as the difference between the potency and log P as illustrated in the following equation: Lip E = pIC50 -cLog P According to data revealed in Table 3, compound 5g showed best value in comparison to other derivatives between both cell lines. highest activity were selected for evaluation against sensitive cancer cell lines (breast and colon cancer cells). LE was calculated using the following equation [20]: LE = (pIC50 × 1.37)/NHA IC50 = half-maximal inhibitory concentration (in terms of molar concentration); NHA = nonhydrogen atom.
The compounds had an LE value in the range of 0.19-0.26 except for compound 5n (Table 3). All compounds exhibited higher LE values in breast cancer cells than in colon cancer cells, especially compounds 5c, 5e, and 5l (LE = 0.26), all of which were structural isomers.
The recommended LE value should be in the range of 0.3. The acceptable LE value should be higher than 0.3.

Lipophilic Efficiency (LipE) or Ligand Lipophilic Efficiency (LEE)
Lip E or LLE is an avenue to determine compound affinity with respect to its lipophilicity. Nowadays, the lipophilic efficiency (LipE) index (LEE), which includes lipophilicity and potency, is becoming more and more popular in drug design. It allows for the normalization of observed potency with changes in the lipophilicity, and it is considered an effective and practical tool for keeping lipophilicity under control to avoid any "molecular obesity".
LipE or LLE is calculated as the difference between the potency and log P as illustrated in the following equation: Lip E = pIC50 -cLog P According to data revealed in Table 3, compound 5g showed best value in comparison to other derivatives between both cell lines. highest activity were selected for evaluation against sensitive cancer cell lines (breast and colon cancer cells). LE was calculated using the following equation [20]: LE = (pIC50 × 1.37)/NHA IC50 = half-maximal inhibitory concentration (in terms of molar concentration); NHA = nonhydrogen atom.
The compounds had an LE value in the range of 0.19-0.26 except for compound 5n (Table 3). All compounds exhibited higher LE values in breast cancer cells than in colon cancer cells, especially compounds 5c, 5e, and 5l (LE = 0.26), all of which were structural isomers.
The recommended LE value should be in the range of 0.3. The acceptable LE value should be higher than 0.3.

Lipophilic Efficiency (LipE) or Ligand Lipophilic Efficiency (LEE)
Lip E or LLE is an avenue to determine compound affinity with respect to its lipophilicity. Nowadays, the lipophilic efficiency (LipE) index (LEE), which includes lipophilicity and potency, is becoming more and more popular in drug design. It allows for the normalization of observed potency with changes in the lipophilicity, and it is considered an effective and practical tool for keeping lipophilicity under control to avoid any "molecular obesity".
LipE or LLE is calculated as the difference between the potency and log P as illustrated in the following equation: Lip E = pIC50 -cLog P According to data revealed in Table 3, compound 5g showed best value in comparison to other derivatives between both cell lines. highest activity were selected for evaluation against sensitive cancer cell lines (breast and colon cancer cells). LE was calculated using the following equation [20]: LE = (pIC50 × 1.37)/NHA IC50 = half-maximal inhibitory concentration (in terms of molar concentration); NHA = nonhydrogen atom.
The compounds had an LE value in the range of 0.19-0.26 except for compound 5n (Table 3). All compounds exhibited higher LE values in breast cancer cells than in colon cancer cells, especially compounds 5c, 5e, and 5l (LE = 0.26), all of which were structural isomers.
The recommended LE value should be in the range of 0.3. The acceptable LE value should be higher than 0.3.

Lipophilic Efficiency (LipE) or Ligand Lipophilic Efficiency (LEE)
Lip E or LLE is an avenue to determine compound affinity with respect to its lipophilicity. Nowadays, the lipophilic efficiency (LipE) index (LEE), which includes lipophilicity and potency, is becoming more and more popular in drug design. It allows for the normalization of observed potency with changes in the lipophilicity, and it is considered an effective and practical tool for keeping lipophilicity under control to avoid any "molecular obesity".
LipE or LLE is calculated as the difference between the potency and log P as illustrated in the following equation: Lip E = pIC50 -cLog P According to data revealed in Table 3, compound 5g showed best value in comparison to other derivatives between both cell lines. highest activity were selected for evaluation against sensitive cancer cell lines (breast and colon cancer cells). LE was calculated using the following equation [20]: LE = (pIC50 × 1.37)/NHA IC50 = half-maximal inhibitory concentration (in terms of molar concentration); NHA = nonhydrogen atom.
The compounds had an LE value in the range of 0.19-0.26 except for compound 5n (Table 3). All compounds exhibited higher LE values in breast cancer cells than in colon cancer cells, especially compounds 5c, 5e, and 5l (LE = 0.26), all of which were structural isomers.
The recommended LE value should be in the range of 0.3. The acceptable LE value should be higher than 0.3.

Lipophilic Efficiency (LipE) or Ligand Lipophilic Efficiency (LEE)
Lip E or LLE is an avenue to determine compound affinity with respect to its lipophilicity. Nowadays, the lipophilic efficiency (LipE) index (LEE), which includes lipophilicity and potency, is becoming more and more popular in drug design. It allows for the normalization of observed potency with changes in the lipophilicity, and it is considered an effective and practical tool for keeping lipophilicity under control to avoid any "molecular obesity".
LipE or LLE is calculated as the difference between the potency and log P as illustrated in the following equation: Lip E = pIC50 -cLog P According to data revealed in Table 3, compound 5g showed best value in comparison to other derivatives between both cell lines. highest activity were selected for evaluation against sensitive cancer cell lines (breast and colon cancer cells). LE was calculated using the following equation [20]: LE = (pIC50 × 1.37)/NHA IC50 = half-maximal inhibitory concentration (in terms of molar concentration); NHA = nonhydrogen atom.
The compounds had an LE value in the range of 0.19-0.26 except for compound 5n (Table 3). All compounds exhibited higher LE values in breast cancer cells than in colon cancer cells, especially compounds 5c, 5e, and 5l (LE = 0.26), all of which were structural isomers.
The recommended LE value should be in the range of 0.3. The acceptable LE value should be higher than 0.3.

Lipophilic Efficiency (LipE) or Ligand Lipophilic Efficiency (LEE)
Lip E or LLE is an avenue to determine compound affinity with respect to its lipophilicity. Nowadays, the lipophilic efficiency (LipE) index (LEE), which includes lipophilicity and potency, is becoming more and more popular in drug design. It allows for the normalization of observed potency with changes in the lipophilicity, and it is considered an effective and practical tool for keeping lipophilicity under control to avoid any "molecular obesity".
LipE or LLE is calculated as the difference between the potency and log P as illustrated in the following equation: Lip E = pIC50 -cLog P According to data revealed in Table 3, compound 5g showed best value in comparison to other derivatives between both cell lines.