Mezna Saleh Altowyan1, Saied M Soliman2, Matti Haukka3, Nora Hamad Al-Shaalan1, Aminah A Alkharboush1, Assem Barakat4. 1. Department of Chemistry, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia. 2. Department of Chemistry, Faculty of Science, Alexandria University, P.O. Box 426, Ibrahimia, Alexandria 21321, Egypt. 3. Department of Chemistry, University of Jyväskylä, P.O. Box 35, Jyväskylä FI-40014 Finland. 4. Department of Chemistry, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia.
Abstract
A new series of spirooxindoles based on ethylene derivatives having furan aryl moiety are reported. The new hybrids were achieved via [3 + 2] cycloaddition reaction as an economic one-step efficient approach. The final constructed spirooxindoles have four contiguous asymmetric carbon centers. The structure of 3a is exclusively confirmed using X-ray single crystal diffraction. The supramolecular structure of 3a is controlled by O···H, H···H, and C···C intermolecular contacts. It includes layered molecules interconnected weak C-H···O (2.675 Å), H···H (2.269 Å), and relatively short Cl···Br interhalogen interactions [3.4500(11)Å]. Using Hirshfeld analysis, the percentages of these intermolecular contacts are 10.6, 25.7, 6.4, and 6.2%, respectively. The spirooxindoles along with ethylene derivatives having furan aryl moiety were assessed against breast (MCF7) and liver (HepG2) cancer cell lines. The results indicated that the new chalcone 3b showed excellent activity in both cell lines (MCF7 and HepG2) with IC50 = 4.1 ± 0.10 μM/mL (MCF7) and 3.5 ± 0.07 μM/mL (HepG2) compared to staurosporine with 4.3 and 2.92 folds. Spirooxindoles 6d (IC50 = 4.3 ± 0.18 μM/mL), 6f (IC50 = 10.3 ± 0.40 μM/mL), 6i (IC50 = 10.7 ± 0.38 μM/mL), and 6j (IC50 = 4.7 ± 0.18 μM/mL) exhibited potential activity against breast adenocarcinoma, while compounds 6d (IC50 = 6.9 ± 0.23 μM/mL) and 6f (IC50 = 3.5 ± 0.11 μM/mL) were the most active hybrids against human liver cancer cell line (HepG2) compared to staurosporine [IC50 = 17.8 ± 0.50 μM/mL (MCF7) and 10.3 ± 0.23 μM/mL (HepG2)]. Molecular docking study exhibited the virtual mechanism of binding of compound 3b as a dual inhibitor of EGFR/CDK-2 proteins, and this may highlight the molecular targets for its cytotoxic activity.
A new series of spirooxindoles based on ethylene derivatives having furan aryl moiety are reported. The new hybrids were achieved via [3 + 2] cycloaddition reaction as an economic one-step efficient approach. The final constructed spirooxindoles have four contiguous asymmetric carbon centers. The structure of 3a is exclusively confirmed using X-ray single crystal diffraction. The supramolecular structure of 3a is controlled by O···H, H···H, and C···C intermolecular contacts. It includes layered molecules interconnected weak C-H···O (2.675 Å), H···H (2.269 Å), and relatively short Cl···Br interhalogen interactions [3.4500(11)Å]. Using Hirshfeld analysis, the percentages of these intermolecular contacts are 10.6, 25.7, 6.4, and 6.2%, respectively. The spirooxindoles along with ethylene derivatives having furan aryl moiety were assessed against breast (MCF7) and liver (HepG2) cancer cell lines. The results indicated that the new chalcone 3b showed excellent activity in both cell lines (MCF7 and HepG2) with IC50 = 4.1 ± 0.10 μM/mL (MCF7) and 3.5 ± 0.07 μM/mL (HepG2) compared to staurosporine with 4.3 and 2.92 folds. Spirooxindoles 6d (IC50 = 4.3 ± 0.18 μM/mL), 6f (IC50 = 10.3 ± 0.40 μM/mL), 6i (IC50 = 10.7 ± 0.38 μM/mL), and 6j (IC50 = 4.7 ± 0.18 μM/mL) exhibited potential activity against breast adenocarcinoma, while compounds 6d (IC50 = 6.9 ± 0.23 μM/mL) and 6f (IC50 = 3.5 ± 0.11 μM/mL) were the most active hybrids against human liver cancer cell line (HepG2) compared to staurosporine [IC50 = 17.8 ± 0.50 μM/mL (MCF7) and 10.3 ± 0.23 μM/mL (HepG2)]. Molecular docking study exhibited the virtual mechanism of binding of compound 3b as a dual inhibitor of EGFR/CDK-2 proteins, and this may highlight the molecular targets for its cytotoxic activity.
In the last decade, spirooxindole
scaffold has been recognized
as a potential pharmacophore in drug discovery, particularly for cancer
research development.[1−3] Spirooxindole has a unique rigid structural architecture
with the diversity of pharmaceutical activities which made this nucleus
a privileged structure in the new drug discovery. The exploration
and discovery of novel drugs for anticancer therapy with lower toxicity
and high selectivity is always an area of intensive research. Many
examples having spirooxindoles have been reported so far for cancer
therapy research and exhibited durable regression for cancer treatment
with oral administration in the preclinical advanced stages; as an
example MI-888 is a lead compound having spiro-pyrrolidinyl oxindole
for p53-MDM2 protein–protein interaction inhibitors.[4] MI-773 and MI-219 are other two representative
examples having spirooxindole-pyrrolidine derivatives which exhibited
cytotoxicity and interfered with the proteasomal degradation of p53.[5,6]Several structural pharmacophores have been coadministrated
with
the spirooxindole privileged structures inspired by research for structural
complexities with their diverse bioactivities.[7,8] Due
to the urgent need to discover a new cancer agent with more targeting
to cancer cells and less harmful for the normal tissue, chemists have
synthesized many spirooxindoles for this purpose. In particular, spirooxindoles
containing furan scaffold, which are called oxa-spirooxindoles, widely
exist in many natural and biologically active molecules.[8−15] The designed, synthesized, and reported spirooxindoles show anticancer
activity and can serve as an anti-tumor agent, CB2 receptor agonist,
antagonists of progesterone receptors, antagonists of progesterone
receptors, Nav1.7 blocker (XEN907), and selective cyclooxygenase COX-1
with TNF-α and IL-6 Inhibitors.[16−19]Spirooxindole core has
been continuously attracting the attention
of researchers and has become a dynamic area of research due to its
outstanding pharmacological properties. Barakat et al. extensively
studied this spirooxindole scaffold recently and have reported so
many examples so far focusing on the drug discovery research.[20−28] In this library, Barakat has reported the synthesis of a new class
of new spiro-heterocycles coadministrated with different pharmacophores
such thiochromene, benzofuran, benzothiophene, cyclohexanone, pyrrole,
and rhodanine scaffolds by the three-component [3 + 2] cycloaddition
reaction in a regio- and stereo-selective fashion. All synthesized
compounds were subjected to anticancer activity against a variety
of cancer cell lines such as PC3, HeLa, MCF-7, MDA-MB231, and so forth,
and many of them showed high efficacy against the tested cell lines.Based on these findings and in continuation of our research program
toward the synthesis of multifunctionalized spirooxindole drug skeleton
for drug research development, we report here the new spirooxindole
system appending the furan structural pharmacophore. The molecular
features of the new spirooxindole derivative compound were elucidated
based on X-ray diffraction of a single crystal, Hirshfeld analysis,
and atoms-in-molecules calculations. Also, the spirooxindoles along
with ethylene derivatives having the aryl-furan moiety were assessed
against two cancer cell lines including breast adenocarcinoma (MCF7)
and liver cancer cell line (HepG2) (Figure ).
Figure 1
Some biologically active spirooxindole-based
pharmacophores.
Some biologically active spirooxindole-based
pharmacophores.
Results and Discussion
Chemistry
New spirooxindoles having
aryl-furan motif were designed and synthesized according to Scheme . The ethylene derivatives 3a–c having the aryl-furan scaffold was mandatory as
a dipolarphile for the [3 + 2] cycloaddition reaction approach and
were prepared from the acetophenones with the aryl-furan carbaldehydes
in basic condition to afford the corresponding chalcones in precipitated
form in a high chemical yield. The ethylene derivative 3a was successfully obtained in a single crystalline form by slow diffusion/evaporation
in DCM/EtOH, and the crystal was suitable for X-ray diffraction analysis.
The required materials for the [3 + 2] cycloaddition reaction were
ethylene derivative having aryl-furan motif 3a–c, four amino acids 5a–d, and two substituted
isatins 4a,b, which achieved the desired spirooxindoles 6a–j. Ten examples were successfully synthesized in
stereoselective and high chemical yield up to 94%. A set of trials
were carried out to get any of those final compounds in a crystalline
form and only one compound (6b) was provided as crystal;
it was not good enough to provide the optimum X-ray data for publication
quality but at least gave us the main skeleton of the final compound.
Other spectrophotometric tools were employed to prove the chemical
structure. The plausible mechanism is depicted in Scheme based on the previous reported
literature.[29−33]
Scheme 1
Synthetic Route for the Ethylene Derivatives 3a–c Engrafted Aryl-Furan Ring
Scheme 2
A Plausible Mechanism for the 32CA Reaction of Azomethine Ylide to
Ethylene Derivative 3a–c to Afford the Spirooxindole
Analogues 6a–j
Crystal Structure Description of 3a
The X-ray structure of 3a is shown in Figure which is found in
good agreement with its spectral characterizations. Compound 3a crystallized in the monoclinic crystal system and the P21/c space group with z = 4 and one molecular unit as asymmetric formula (Table ). The unit cell parameters
are a = 19.4870(9) Å, b = 13.8512(6)
Å, c = 5.8412(3) Å, β = 92.720(6)°,
and V = 1574.87(13) Å3. The list
of bond distances and angles is given in Table . The structure comprised three aromatic
planar ring systems. These rings are abbreviated A, B, and C for simplicity (Figure ). It is worth to note that the three rings
are twisted with one another to different extents. The mean plane
of ring A is found twisted with respect to ring B by 12.95°. The corresponding values for ring C is 8.8° with respect to the mean plane of ring B.
Figure 2
Thermal ellipsoids at 30% probability level showing atom numbering
of 3a.
Table 1
Crystal
Data
3a
CCDC
2165902
empirical formula
C19H12BrClO2
fw
387.65
temp (K)
120(2)
λ (Å)
0.71073
cryst syst.
monoclinic
space group
P21/c
a (Å)
19.4870(9)
b (Å)
13.8512(6)
c (Å)
5.8412(3)
β (deg)
92.720(6)
V (Å3)
1574.87(13)
Z
4
ρcalc (Mg/m3)
1.635
μ (Mo Kα) (mm–1)
2.786
no. reflns.
14253
unique reflns.
3887
completeness to θ = 25.242°
99.8%
GOOF (F2)
1.065
Rint
0.0558
R1a (I ≥ 2σ)
0.0547
wR2b (I ≥ 2σ)
0.1134
R1 =
Σ||Fo| – |Fc||/Σ|Fo|.
wR2 =
{Σ[w(Fo2 – Fc2)2]/Σ[w(Fo2)2]}1/2.
Table 2
Selected Bond Lengths [Å] and
Angles [deg] for 3a
atoms
distance
atoms
distance
Br(1)–C(1)
1.889(4)
O(2)–C(13)
1.368(4)
Cl(2)–C(17)
1.732(4)
O(2)–C(10)
1.380(4)
O(1)–C(7)
1.234(4)
Thermal ellipsoids at 30% probability level showing atom numbering
of 3a.R1 =
Σ||Fo| – |Fc||/Σ|Fo|.wR2 =
{Σ[w(Fo2 – Fc2)2]/Σ[w(Fo2)2]}1/2.The supramolecular
structure of 3a is controlled by
different types of intermolecular contacts (Figure ). The most important contacts are emphasized
by different colors in the packing scheme shown in Figure . H···O (magenta),
H···H (purple), C···C (orange), and
Br···Cl (turquoise) are the main contacts in this crystal
structure. The presence of short Br···Cl (3.450 Å)
and C···C (C8vC15; 3.323 Å and C9···C13;
3.399 Å) contacts revealed the presence of significant interhalogen
and π–π stacking interactions, respectively. The
rest of the interactions are O···H and H···H
and are depicted in Table .
Figure 3
Most important intermolecular contacts in the crystal structure
of 3a.
Figure 4
Packing of the molecular
units via H···O (magenta),
H···H (purple), C···C (orange), and
Br···Cl (turquoise) in 3a.
Table 3
Different Contacts and Their Distances
(Å)
contact
distance
symm. code
Br1···Cl2
3.4500(11)
–1 + x, y, –1 + z
C8···C15
3.323(4)
1 – x, 1 – y, 1 – z
C9···C13
3.399(4)
1 – x, 1 – y, 1 – z
O1···H16
2.675
1 – x, 1 – y, 2 – z
H3···H16
2.269
1 – x, 1 – y, 2 – z
Most important intermolecular contacts in the crystal structure
of 3a.Packing of the molecular
units via H···O (magenta),
H···H (purple), C···C (orange), and
Br···Cl (turquoise) in 3a.
Hirshfeld Surface Analysis
Crystal
structure stability is governed by intermolecular interactions among
different fragments in the crystal. Hirshfeld analysis is a simple
tool for decomposition of different intermolecular contacts in the
crystal. Different Hirshfeld surfaces including dnorm, shape index, and curvedness maps are shown in Figure . The dnorm map contains a number of red spots which represent
the regions at which the important and short distance contacts occurred.
Summary of all contacts and their percentages based on Hirshfeld calculations
is presented in Figure .
Figure 5
Hirshfeld surfaces of 3a.
Figure 6
Intermolecular
contacts and their percentages in 3a.
Hirshfeld surfaces of 3a.Intermolecular
contacts and their percentages in 3a.Analysis of these interactions using fingerprint plot and dnorm maps is given in Figure . The major contacts are H···H,
H···C, and O···H interactions, and their
percentages are 25.7, 27.6, and 10.6%, respectively. It is worth noting
that the red spots shown in the dnorm map
are related to the Br···Cl (6.2%), C···C
(6.4%), and H···H contacts. Hence, the hydrogen–hydrogen
interactions are not only the most common but also considered strong.
The H3···H16 contact (2.06 Å) is the shortest
of these interactions. Also, the presence of short Br1···Cl2
(3.45 Å) revealed the presence of interhalogen interactions,
whereas the short C13···C9 (3.399 Å) and C16···C9
(3.323 Å) contacts revealed the importance of the π–π
stacking interactions. The latter is further confirmed by the presence
of red/blue triangles in the shape index map (Figure ). The O1···H16 contact appeared
as a blue region in the dnorm map. The
corresponding interaction distance is found to be 2.674 Å based
on Hirshfeld analysis. This value is slightly greater than the vdWs
radii sum (2.61 Å) of the O and H atoms. Other contacts shown
in this figure are considered weak interactions and have less contribution
in the molecular packing of 3a.
Figure 7
Decomposed dnorm maps (A) and fingerprint
plots (B) of short contacts in 3a.
Decomposed dnorm maps (A) and fingerprint
plots (B) of short contacts in 3a.
Biological Activity of the Synthesized Compounds
The ethylene derivatives having the aryl-furan moiety and the spirooxindoles
were assessed against two cancer cell lines including breast and liver
carcinoma by the MTT assay. Interestingly, the new chalcone 3b discovered is the most active member between the synthesized
compounds against both cell lines, with IC50 = 4.1 ±
0.10 and 3.5 ± 0.07 μM/mL for MCF7 and HepG2, respectively,
and more potent than the standard drug used as staurosporine [IC50 = 17.8 ± 0.50 μM/mL (MCF7) and 10.3 ± 0.23
μM/mL (HepG2)]. The other two chalcones 3a and 3c show moderate activity (Table ).
Table 4
Cytotoxicity Results
of the Ethylene
Derivatives Having the Aryl-Furan Moiety and the Spirooxindoles against
Breast and Liver Carcinoma
IC50 values are expressed
as mean ± SD of three independent trials.
IC50 values are expressed
as mean ± SD of three independent trials.Although compounds 3a and 3b are two
isomers with different positions for chlorine substitutions (o/p-substitution),
they exhibited different cytotoxicity against the two tested cell
lines. This may be due to the relevance of chlorine substitution to
interact with the corresponding amino acids in the effective target
binding site, which will affect the stability of drug–target
complexes and the ability of the complex to have a biological response.For the synthesized spirooxindoles, for the breast cancer cell
line (MCF7), spirooxindole having the bulky fused with the pyrrolidine
ring 6d exhibited the most active hybrid between these
series with IC50 = 4.3 ± 0.18 μM/mL more potent
with 4 folds and the standard staurosporine (IC50 = 17.8
± 0.50 μM/mL). Next in the reactivity towards cytotoxicity
was compound spirooxindole 6j which instead of the Br-atom
on the benzoyl ring in compound 6d replaced by more electron
withdrawing group, the reactivity slightly decreased which the IC50 equal to 4.7 ± 0.18 μM/mL. Compounds 6f (IC50 = 10.3 ± 0.40 μM/mL) and 6i (IC50 = 10.7 ± 0.38 μM/mL) still provided
better cytotoxicity compared to staurosporine. The remaining spirooxindoles
provided moderate toxicity against breast adenocarcinoma. Spirooxindoles 6f (IC50 = 3.5 ± 0.11 μM/mL) and 6d (IC50 = 6.9 ± 0.23 μM/mL) were the
most active hybrids against human liver cancer cell line (HepG2) compared
to staurosporine [IC50 = 10.3 ± 0.23 μM/mL (HepG2)].
The remaining spirooxindoles provided moderate toxicity, and the IC50 ranged from 11.8 ± 0.37 to 91.4 ± 2.86 μM/mL.
Molecular Docking Study
As cyclic-dependent
kinase (CDK-2) and epidermal growth factor receptor (EGFR) are key
proteins that mediate and trigger the proliferation of cancer cells,
their inhibition is an interesting target for apoptosis induction
upon treatment of a chemotherapeutic drug.[34−36] A molecular
docking study was performed to highlight the virtual mechanism of
binding toward EGFR and CDK-2 proteins. As seen in Figure , compound 3b was
docked inside the EGFR binding site with the binding energy of −19.63
kcal/mol, and it formed one H-bind interaction with Met 769 as the
H-bond acceptor. Additionally, compound 3b was docked
inside the CDK-2 active site with the binding energy of −18.6
Kcal/mol, and it formed one H-bond with Lys 89. Hence, compound 3b had a good binding affinity toward EGFR and CDK-2 proteins,
and these targets may be the effective target for its cytotoxic activity.
Figure 8
Binding
disposition and ligand–receptor interactions of
compound 3b inside the EGFR (A) and CDK-2 (B) proteins.
Three-dimensional images were generated by Chimera-UCSF. Cocrystallized
ligand (yellow-colored) and docked compound (cyan-colored).
Binding
disposition and ligand–receptor interactions of
compound 3b inside the EGFR (A) and CDK-2 (B) proteins.
Three-dimensional images were generated by Chimera-UCSF. Cocrystallized
ligand (yellow-colored) and docked compound (cyan-colored).
Conclusions
We concluded
that a set of spirooxindoles having aryl-furan moiety
have been successfully synthesized and assessed against two cancer
cell lines. The results provided promising data for breast adenocarcinoma
and could be utilized as a lead compound for further development.
Structural analysis of the newly synthesized compound was performed
with the aid of X-ray single crystal structure and Hirshfeld calculations.
The reported X-ray structure of 3a agreed very well with
its spectral characterizations. Its supramolecular structure is controlled
by many intermolecular contacts such as O···H, H···H,
and C···C intermolecular contacts as well as the Br···Cl
interhalogen interactions. All intermolecular contacts occurring in
the crystal are quantitatively analyzed based on Hirshfeld calculations.
The percentages of these contacts are 10.6, 25.7, 6.4, and 6.2%. The
biological activity concluded that the new chalcone 3b showed excellent activity in both MCF7 and HepG2 with IC50 = 4.1 ± 0.1 μM/mL (MCF7) and 3.5 ± 0.07 μM/mL
(HepG2) compared to staurosporine with 4.3 and 2.92 folds. The synthesized
spirooxindole compounds 6d (IC50 = 4.3 ±
0.18 μM/mL), 6f (IC50 = 10.3 ±
0.40 μM/mL), 6i (IC50 = 10.7 ±
0.38 μM/mL), and 6j (IC50 = 4.7 ±
0.18 μM/mL) show high antiproliferative activity in
vitro against breast adenocarcinoma. On the other hand, compounds 6d (IC50 = 6.9 ± 0.23 μM/mL) and 6f (IC50 = 3.5 ± 0.11 μM/mL) were the
most active hybrids against human liver cancer cell line (HepG2) compared
to staurosporine [IC50 = 17.8 ± 0.50 μM/mL (MCF7)
and 10.3 ± 0.23 μM/mL (HepG2)]. These compounds could be
useful for further cancer research development.
Materials
and Methods
All technical instruments and chemicals used
in this study are
provided in the Supporting Information.
General Procedure for the Synthesis of the
Chalcones (Ethylene Derivatives 3a–c)
The synthesis of ethylene derivatives 3a–c was
performed according to the reported procedure by mixing equimolar
of the acetophenone derivatives 1a,b [1-(4-bromophenyl)ethan-1-one 1a; 1-(4-nitrophenyl)ethan-1-one 1b] with the
corresponding aldehydes 2a,b [5-(4-chlorophenyl)furan-2-carbaldehyde 2a; 5-(2-chlorophenyl)furan-2-carbaldehyde 2b] in the presence of basic condition (NaOH, 2 equiv) in ethanolic
solution to give the corresponding ethylene derivatives in a precipitated
form which were used for the next step without any further purification.
General
Method for the Synthesis of the Spirooxindoles
Scaffold 6a–j
An equimolar (0.5 mmol)
of the ethylene derivatives 3a–c with isatin derivatives 4a,b and amino acids 5a–d in methanol
was heated under reflux for 5 h to provide the final products 6a–j as solid materials in almost quantitative yield
upon slow evaporation overnight. In the case of compound 6b, the crystal quality was not good enough to solve the X-ray structure
in a suitable form enough for publication quality but at least provided
a reasonable structure for the spiro-compound.
The
crystal of 3a was immersed in cryo-oil, mounted in a
loop, and measured at a temperature of 120 K. The X-ray diffraction
data were collected on a Rigaku Oxford Diffraction Supernova diffractometer
using Mo Kα radiation. The CrysAlisPro[37] software package was used for cell refinement and data reduction.
An analytical absorption correction (CrysAlisPro[37]) was applied to the intensities before structure solution.
The structure was solved by intrinsic phasing (SHELXT[38]) method. Structural refinement was carried out using SHELXL[39] software with the SHELXLE[40] graphical user interface. Hydrogen atoms were positioned
geometrically and constrained to ride on their parent atoms, with
C–H = 0.95 Å and Uiso = 1.2Ueq (parent atom).
Hirshfeld
Surface Analysis
Crystal
Explorer 17.5 program[41] was used to perform
the Hirshfeld topology analysis.
Molecular
Docking Study
The investigated
compounds were docked against the protein structures of EGFR (PDB
= 1M17) and
CDK-2 (PDB = 2A4L) using AutoDock Vina software following routine work.[42−45] Vina was used to improve protein and ligand structures and favor
them energetically. Proteins and compound structures were prepared
and optimized using Maestro. Then, the binding sites inside proteins
were determined using grid-box dimensions around the cocrystallized
ligands. Binding activities interpreted molecular docking results
in terms of binding energy and ligand–receptor interactions.
The visualization was then done with Chimera.
Authors: Gehad Lotfy; Mohamed M Said; El Sayed H El Ashry; El Sayed H El Tamany; Abdullah Al-Dhfyan; Yasmine M Abdel Aziz; Assem Barakat Journal: Bioorg Med Chem Date: 2017-01-16 Impact factor: 3.641
Authors: Heba S A ElZahabi; Mohamed S Nafie; Dina Osman; Nehal H Elghazawy; Dalia H Soliman; Abdelghany Ali H El-Helby; Reem K Arafa Journal: Eur J Med Chem Date: 2021-06-05 Impact factor: 6.514
Authors: Yasmine M Abdel Aziz; Gehad Lotfy; Mohamed M Said; El Sayed H El Ashry; El Sayed H El Tamany; Saied M Soliman; Marwa M Abu-Serie; Mohamed Teleb; Sammer Yousuf; Alexander Dömling; Luis R Domingo; Assem Barakat Journal: Front Chem Date: 2021-12-14 Impact factor: 5.221
Authors: Assem Barakat; Matti Haukka; Saied M Soliman; M Ali; Abdullah Mohammed Al-Majid; Ayman El-Faham; Luis R Domingo Journal: Molecules Date: 2021-11-30 Impact factor: 4.411