Christina L Koumpoura1, Michel Nguyen1, Christian Bijani1, Laure Vendier1, Elena G Salina2, Silvia Buroni3, Giulia Degiacomi3, Sandrine Cojean4, Philippe M Loiseau4, Françoise Benoit-Vical1, Alfonso T García-Sosa5, Michel Baltas1. 1. Laboratoire de Chimie de Coordination du CNRS-UPR8241, Inserm ERL 1289 Team "New antiplasmodial molecules and pharmacological approaches", 205 route de Narbonne, BP 44099, Toulouse Cedex 31077, France. 2. Bach Institute of Biochemistry, Research Center of Biotechnology of the Russian Academy of Sciences, Moscow 119071, Russia. 3. Department of Biology and Biotechnology "Lazzaro Spallanzani", University of Pavia, Pavia 27100, Italy. 4. Antiparasite Chemotherapy, UMR 8076 CNRS BioCIS, Faculty of Pharmacy, University Paris-Saclay, Châtenay-Malabry 92290, France. 5. Department of Molecular Technology, Institute of Chemistry, University of Tartu, Ravila 14a, Tartu 50411, Estonia.
Abstract
The first effective synthetic approach to naphthofuroquinones via a reaction involving lawsone, various aldehydes, and three isocyanides under microwave irradiation afforded derivatives in moderate to good yields. In addition, for less-reactive aldehydes, two naphtho-enaminodione quinones were obtained for the first time, as result of condensation between lawsone and isocyanides. X-ray structure determination for 9 and 2D-NMR spectra of 28 confirmed the obtained structures. All compounds were evaluated for their anti-infectious activities against Plasmodium falciparum, Leishmania donovani, and Mycobacterium tuberculosis. Among the naphthofuroquinone series, 17 exhibited comparatively the best activity against P. falciparum (IC50 = 2.5 μM) and M. tuberculosis (MIC = 9 μM) with better (P. falciparum) or equivalent (M. tuberculosis) values to already-known naphthofuroquinone compounds. Among the two naphtho-enaminodione quinones, 28 exhibited a moderate activity against P. falciparum with a good selectivity index (SI > 36) while also a very high potency against L. donovani (IC50 = 3.5 μM and SI > 28), rendering it very competitive to the reference drug miltefosine. All compounds were studied through molecular modeling on their potential targets for P. falciparum, Pfbc1, and PfDHODH, where 17 showed the most favorable interactions.
The first effective synthetic approach to naphthofuroquinones via a reaction involving lawsone, various aldehydes, and three isocyanides under microwave irradiation afforded derivatives in moderate to good yields. In addition, for less-reactive aldehydes, two naphtho-enaminodione quinones were obtained for the first time, as result of condensation between lawsone and isocyanides. X-ray structure determination for 9 and 2D-NMR spectra of 28 confirmed the obtained structures. All compounds were evaluated for their anti-infectious activities against Plasmodium falciparum, Leishmania donovani, and Mycobacterium tuberculosis. Among the naphthofuroquinone series, 17 exhibited comparatively the best activity against P. falciparum (IC50 = 2.5 μM) and M. tuberculosis (MIC = 9 μM) with better (P. falciparum) or equivalent (M. tuberculosis) values to already-known naphthofuroquinone compounds. Among the two naphtho-enaminodione quinones, 28 exhibited a moderate activity against P. falciparum with a good selectivity index (SI > 36) while also a very high potency against L. donovani (IC50 = 3.5 μM and SI > 28), rendering it very competitive to the reference drug miltefosine. All compounds were studied through molecular modeling on their potential targets for P. falciparum, Pfbc1, and PfDHODH, where 17 showed the most favorable interactions.
Known
as dyes in the chemical industry, quinones are also an important
family of major impact as biologically active compounds.[1] The quinone fragment is present in many natural
products with important biological roles in humans, animals, and plants,[2] rendering it one of the most attractive structures
in medicinal chemistry. In this respect, various non-natural molecules
bearing the quinone scaffold have been prepared and evaluated for
anticancer, antifungal, antimalarial, and other biological activities.[3−7]Among various classes of quinones, naphthoquinones and especially
naphthofuroquinones have attracted increasing interest for their pharmacological
activities. Naphthofuroquinone derivatives can be obtained from natural
sources[8,9] but also via synthetic approaches, giving
rise to compounds with various structures and biological properties,
for example, cytotoxicity against tumor cells,[10] inhibition of tyrosine kinases receptors,[11] and antimalarial[12] or antidiabetic
activity[13] (Figure ).
Figure 1
Examples of naphthofuroquinones with interesting
biological activities.
The structure of the starting material 2-hydroxy-1,4-naphthoquinone
(1, lawsone) is given for comparison.
Examples of naphthofuroquinones with interesting
biological activities.
The structure of the starting material 2-hydroxy-1,4-naphthoquinone
(1, lawsone) is given for comparison.Due to their pharmacological activities, many synthetic efforts
have been dedicated to their preparation. Two main synthetic routes
have been adopted: (i) reaction of 2,3-dichloro-1,4-naphthoquinone
with 1,3-dicarbonyl reagents[14] and (ii)
strategies based on the [3 + 2] annulation of 2-hydroxy-1,4-naphthoquinone[15] with various reagents.[16−18] We can also
refer to two other synthetic approaches using either 3-iodolawsone[19] or thio-substituted 1,4-naphthoquinone.[20]Lawsone seems to be an important fragment
while searching for antimalarial
agents. Atovaquone (7) (Figure )[21] is reported
to be the leading drug interacting with the cytochrome bc1 complex, a specific mitochondrial target of P. falciparum.[22] Cytochrome bc1 and
dihydroorotate dehydrogenase (DHODH)[23−25] are the two main biological
targets of Plasmodium mitochondrion,
and they are also biochemically related. Atovaquone acts as a competitive
inhibitor at the Q0 site of the cytochrome bc1 complex. Unfortunately, strong P. falciparum resistance against atovaquone arises rapidly in the field.
Figure 2
Characteristic
examples of lawsone derivatives with antimalarial
activities against P. falciparum.
Characteristic
examples of lawsone derivatives with antimalarial
activities against P. falciparum.In this respect, efforts to design
and prepare naphthoquinones
and naphthofuroquinones based on lawsone scaffold have been pursued
worldwide:[26]Borgati et al.[12] reported the synthesis
of a series of naphthoquinones and naphthofuroquinones through functionalization
of lawsone either on the 2-hydroxy group affording alkoxy derivatives
or on the C-3 alkene carbon atom, followed by cyclization, thus leading
to naphthofuroquinones. Among all the synthesized compounds, naphthofuroquinone 6 exhibited the lowest IC50 value against the chloroquine-resistant P. falciparum W2 strain (Figure ). Based on docking studies, the authors
proposed that compounds of this series, especially 6,
could act as inhibitors of either cytochrome bc1 or
DHODH.Oramas-Royo et al.[27] reported
the synthesis
and antiplasmodial activities of a series of 1,2,3-triazole-naphthoquinone
conjugates. Compound 8 (Figure ) exhibited the second lowest IC50 value against the chloroquine-resistant P. falciparum W2 strain, while molecular docking on the potential target PfDHODH
showed very favorable interactions.In this context, we report
our efforts in constructing naphtho(furo)quinone
derivatives obtained via a three (or two)-component reaction of lawsone
with aldehydes and isocyanides. First, their syntheses will be described.
Then, their biological activities will be reported against three important
pathogens, P. falciparum, L. donovani, and M. tuberculosis, for which the problem of multiresistance to the current therapy
is pressing. In addition, we report docking studies of several drugs
on cytochrome bc1 and DHODH molecular targets.
Results
and Discussion
Chemistry
Multicomponent reactions
are an important
synthetic tool in medicinal chemistry.[28] They are one-pot reactions employing more than two starting compounds,
where the majority of atoms are incorporated into the final compound.[29] Nowadays, many basic multicomponent reactions
are named ones (Ugi, Passerini, Strecker, Biginelli, etc.), while
many others have been developed for the synthesis of heterocyclic
compounds such as 2-styryl quinolines, quinazolines,[30,31] and imidazoles.[32] Concerning the naphthofuran
system, synthesis of 2-amino-naphtho[2,3-b]furan-4,9-diones
has been reported by Teimouri et al.[33] as
one-pot three-component condensation among lawsone (1), an aldehyde, and an alkyl-isocyanide under reflux in toluene.
Later on, Jiménez-Alonso et al.[34] reported a methodology using ethylenediaminediacetic acid (EDDA)
as a catalyst, thus reducing the reflux time in toluene.The
mechanistic aspect of the synthesis of these adducts is considered
as follows: first, a Knoevenagel condensation between lawsone (1) and an aldehyde leads to the formation of a conjugated
enone; this step is followed by a [4 + 1] cycloaddition reaction of
the heterodiene moiety of the Knoevenagel adduct with an isocyanide.
This reaction affords an intermediate iminolactone that leads to the
final naphthofuroquinone derivative after rearrangement (Scheme ).
Scheme 1
Possible Mechanism
of the Multicomponent Domino Reaction of Lawsone
To optimize the reaction conditions (catalyst, solvent,
reaction
time, and activation method), the reaction between lawsone (1), 4-chlorobenzaldehyde, and tert-butyl
isocyanide (Scheme and Table ) was
first used as a model reaction. In fact, embelin has been recently
reported[35] to undergo this three-component
reaction under microwave irradiation. However, until now, there was
no equivalent study for lawsone (1). Microwave-assisted
organic synthesis has attracted great attention in the last decades
and has proven to be particularly useful in the construction of important
heterocyclic systems[36] such as benzothiazoles[37] and benzoxazoles.[38] Microwave-assisted multicomponent reactions are also part of this
progress, having been emerged as useful tools for the elaboration
of relevant heterocycles.[39]
Scheme 2
Model Domino
Reaction Used for the Optimization of Reaction Conditions
Table 1
Reaction Conditions and Yields for
the Model Domino Reactionc
methodology
reaction time
equivalents (lawsone/aldehyde: isocyanide)
catalyst
solvent
lawsone
conversiona
isolated
yieldb
reflux
120 min
1:1.2:1.2
EDDA
PhMe
90%
30%
irradiation (μw)
2 × 15 min
1:1:1
Et3N
DCE
12%
11%
irradiation (μw)
60 min
1:1:1
EDDA
DCE
60%
59%
irradiation
(μw)
2 × 60 min
1:1.5:1
EDDA
DCE
100%
66%
thermal activation
60 min
1:1.5:1
EDDA
DCE
100%
59%
thermal activation
60 min
1:1:1
Et3N
DCE
70%
57%
Calculated
based on the 1H NMR spectrum of the crude reaction mixture
at the end of the reaction.
Calculated based on the starting
amount of lawsone.
μw
stands for microwave
Calculated
based on the 1H NMR spectrum of the crude reaction mixture
at the end of the reaction.Calculated based on the starting
amount of lawsone.μw
stands for microwaveThe
reaction was first conducted in the presence of solvents (toluene
and 1,2-dichloroethane DCE) and catalysts (Et3N and EDDA),
favoring the Knoevenagel condensation. All reactions were conducted
under reflux for 2 h and monitored by thin-layer chromatography (TLC)
every 30 min. Our best result under these conventional conditions
was obtained when using EDDA in toluene: the obtained yield of the
final compound 9 was 30%. A monowave 50 thermic reactor
was also used. The reaction in DCE in the presence of 10% of EDDA
(1 h of heating at 182 °C, 15 bar) provided the target naphthofuroquinone
derivative 9 in 59% yield after PuriFlash purification.Then, we performed the reaction under microwave irradiation for
1 h (160 °C, 6 bars, 300 W). When operating in the presence of
1:1:1 mole equiv of reactants, after work up and purification, we
obtained the deep violet-colored naphthofuroquinone derivative 9 in 59% yield. Some traces of a yellow compound were also
isolated, without further identification until now. By increasing
the ratio of aldehyde to 1.5 equiv, the yield of 9 was
slightly increased up to 66%. The use of more equivalents of aldehyde
did not lead to better results. Conversely, the purification and isolation
of the naphthofuroquinone became problematic in the presence of excess
of aldehyde since the aldehyde was eluted in very near proximity to
the desired naphthofuroquinone.The microwave irradiation using
1:1:1 equiv of reactants and 10%
of EDDA as a catalyst was therefore considered as the optimal conditions.
Then, we performed this reaction with a series of aldehydes and three
different isocyanides, namely, tert-butyl, n-butyl, and cyclohexyl isocyanides. Table summarizes the obtained results.
Table 2
Structures and Isolated Yields of
Productsa
Calculated
based on the 1H NMR spectrum of the crude reaction mixture
at the end of the reaction.
Calculated based on the starting amount of lawsone.
Calculated
based on the 1H NMR spectrum of the crude reaction mixture
at the end of the reaction.
Calculated based on the starting amount of lawsone.When using p-chlorobenzaldehyde
(entry 1), the
domino adduct 9 was obtained in 59% yield with respect
to the starting amount of aldehyde (the conversion of aldehyde was
70%). A similar result was evidenced when using 3,4,5-trimethoxy benzaldehyde
(entry 3): the domino reaction afforded adduct 11 in
55% (with aldehyde conversion = 60%). Reaction with m-bromobenzaldehyde (entry 2) gave a more complex mixture, leading
after purification to adduct 10 in 30% yield. Despite
the fact that the aldehyde conversion was rather poor (30%), we managed
to isolate a second product from the reaction mixture. It was identified
as compound 28 (yield 20%), where the tert-butyl isocyanide has reacted directly with lawsone (entry 20). Noteworthily,
this compound is a mixture of cis- and trans-isomers in a 1:1 molar ratio. This point will be further discussed
below. Concerning the reaction with 3,4-dimethoxy benzaldehyde, the
aldehyde conversion was only 45% and a complex mixture was obtained
in which compound 28 (entry 20) was not detected. Purification
led to compound 12 in 25% yield (entry 4).Two
heteroaromatic bicyclic aldehydes were then chosen to carry
out the domino reaction. The reaction of N-methyl
indazole aldehyde in the presence of lawsone and tert-butyl isocyanide (ratio 1:1:1) in the presence of EDDA (10%) under
microwave irradiation for 2 h provided the domino compound 13 with 30% yield after purification (the aldehyde conversion was 30%),
indicating a high selectivity for this domino reaction (entry 5).
In addition, the same adduct 28 (cis- and trans-isomers) as above was isolated in 40%
yield. On the other hand, reaction of indole aldehyde (entry 6) afforded
a complex mixture. A poor yield of the domino adduct 14 (35%) with an 80% aldehyde conversion was obtained, while compound 28 (entry 20) was not detected. Finally, starting from 5-nitrofuraldehyde,
the domino compound 15 was obtained in poor yield (19%)
with 20% aldehyde conversion (entry 7).Piperidine-4-carbaldehydes
possessing either a Boc or a benzyl
protecting group on the nitrogen atom provided very similar results
(entries 8 and 9). In both cases, the corresponding naphthofuroquinones 16 and 17 were obtained in 52 and 53% yield,
respectively, after purification. Compound 28 (entry
20) was not detected. The conversion of the aldehyde was higher than
90%.The three-component reaction was also studied with commercially
available cyclohexyl and n-butyl isocyanides. All
reactions were conducted with five aldehydes representative of those
studied above: N-benzylpiperidine carbaldehyde, p-chlorobenzaldehyde, 3,4,5-trimethoxybenzaldehyde, N-methyl-indazole carbaldehyde, and 5-nitrofuraldehyde.In the reactions of cyclohexyl isocyanide (entries 10–14)
and using p-chlorobenzaldehyde, 3,4,5-trimethoxybenzaldehyde, N-methyl-indazole, or N-benzyl piperidine
carbaldehyde, the corresponding naphthofuroquinones 18–21 were obtained in 46, 42, 40, and 42% yields,
respectively (entries 10–12 and 14). Finally, the reactions
of cyclohexyl isocyanide with 5-nitrofuraldehyde provided a complex
mixture from which 22 was isolated in 12% yield (entry
13).Concerning reactions with n-butyl isocyanide
(entries
15–19), the best results in terms of yield and aldehyde conversion
were obtained with aromatic aldehydes (entries 15 and 16). Naphthofuroquinones 23 and 24 were obtained in 55 and 40% yield,
respectively, after purification. N-Benzyl piperidine
carbaldehyde led to a 30% yield of naphthofuroquinone 27 with 35% aldehyde conversion (entry 19), while 5-nitrofuraldehyde
resulted in only 10% of aldehyde conversion and afforded compound 26 (entry 18) in very poor yield in a complex mixture. Reaction
of n-butyl isocyanide with N-methyl-indazole
carbaldehyde afforded compound 29 with 50% yield (entry
21), while naphthofuroquinone 25 could be identified
in a complex mixture. Compound 25 was unstable in different
purification conditions.To summarize, a library of 19 novel
naphthofuroquinones was synthesized
under microwave irradiation in moderate yields (three-component reaction, Scheme ). In some cases,
where probably aldehydes are less reactive, a two-component condensation
between lawsone and tert-butyl or n-butyl isocyanide occurred, affording novel naphtho-enaminodione
quinones 28 and 29 (two-component reaction, Scheme ).
Scheme 3
Summary of All Types
of Condensation Products Obtained
X-ray and 2D-NMR Analysis of Compound 9
Compound 9 was recrystallized in methylene chloride.
The obtained single crystals, which appeared as purple platelets,
were analyzed by X-ray diffraction (Figure ). Mo-Kα radiation was used, and the
compound showed a P2/ space group (crystallographic Data
in the Supporting Information). Inspection
of the structure shows that the aromatic halogen ring exhibited a
dihedral angle of −44.2° with respect to the furan ring
(C6b/C1/C3/C2), while
the carbon atom of the tert-butyl group (C10) was almost coplanar with the furan ring. The observed dihedral
angle (O1/C2/N1/C10) had
a small value of 16.3°, thus indicating a potential useful space
for coordination between the oxygen atom of the carbonyl group (C9 = O2) and the oxygen atom of the furan ring (O1).
Figure 3
X-ray structure of compound 9.
X-ray structure of compound 9.Compound 9 was also analyzed by 2D-NMR at 298 K. All
the 1H and 13C signals were assigned based on
the chemical shifts, spin–spin coupling constants, splitting
patterns, and signal intensities and by using 1H–1H COSY45, 1H–13C HSQC, and 1H–13C HMBC experiments (see the Experimental Section).In addition, the
attribution of C4 and C9 of compound 9 was based on the
previous study of Borgati et al.,[40] where
the authors performed a complete, comparative 1H and 13C signal assignment of p-naphthoquinones
and ortho- and para-naphthofuroquinones.
In their case, the presence of a proton at
C3 allowed the differentiation/identification of C4 and C9 by HMBC
correlation with a Δd between C4 and C9 as
high as 8 ppm. Then, we assigned the chemical shifts of C9 and C4
carbon atoms at 169.6 and 182.0 ppm, respectively, which is consistent
with the reliable assignment of Borgati et al.The 1H and 13C chemical shifts of 9 are given
in Table . The 13C and HSQC spectra show 20 different carbon signals
for 9. From the 13C data, it was possible
to identify one cyclohexanedione (δC 169.58 and 181.95),
one phenyl chloride group (δC 129.2, 130.8, and 133.8),
six sp2-hybridized carbons (δC from 99.0 to 159.3),
and four sp3-hybridized carbons (δC 30.1 to 54.1).
Table 3
1H and13C NMR
Data Assignments of Compound 9 in CDCl3 at
298 Ka
1H and/or13C numbering
1H chemical shift, ppm
13C chemical shift, ppm
C-2
159.3
C-9a
144.2
C-9
169.6
C-8a
133.2
CH-8
8.19
126.1
CH-7
7.72
133.8
CH-6
7.63
132.4
CH-5
8.04
126.5
C-4a
133.2
C-4
182.0
C-3a
130.3
C-3
99.0
C-1
133.6
CH-2b
7.44
130.8
CH-3b
7.47
129.2
C-4b
133.8
CH-5b
7.47
129.2
CH-6b
7.44
130.8
NH
4.86
C-10
54.1
CH3-11,12,13
1.52
30.1
For the sake of clarity, the chemical
shifts are reported starting from the furan (C2), then the naphthoquinone
moiety (C9a to C3a), followed by the p-chlorophenyl
ring (C1 to C6b), and finally the N-tert-butyl group (HN1, C10 to C13).
For the sake of clarity, the chemical
shifts are reported starting from the furan (C2), then the naphthoquinone
moiety (C9a to C3a), followed by the p-chlorophenyl
ring (C1 to C6b), and finally the N-tert-butyl group (HN1, C10 to C13).
2D-NMR Analysis of Compound 28
A complete
2D-NMR analysis was also conducted for compounds 28 (28a/28a′) and 29 (29a/29a′), especially using HMBC correlations. Both compounds are present
in two stereoisomers in a 1:1 ratio.Compound 28a/28a′ is analyzed here (for compound 29a/29a′, see
the Supporting Information) at a temperature
of 298 K. All the 1H and 13C signals were assigned
on the basis of chemical shifts, spin–spin coupling constants,
splitting patterns, and signal intensities and by using 1H–1H COSY45, 1H–13C HSQC, and 1H–13C HMBC experiments
(see the Experimental Section). The 1H and 13C chemical shifts of compound 28a/28a′ are given in Table along with its structure indicating the evidenced correlations.
Table 4
1H and 13C NMR
Data Assignments of 28a/28a′ in CDCl3 at 298 Ka
1H and/or 13C numbering
1H chemical shift, ppm
13C chemical shift, ppm
28a
28a′
28a
28a′
C-7
108.6
109.3
C-8
176.8
177.5
C-9
181.7
182.0
C-12
132.3
132.3
CH-13
8.18
8.19
128.3
127.7
CH-14
7.72
7.72
133.3
133.4
CH-15
7.84
7.82
135.6
135.4
CH-16
8.16
8.27
126.5
127.3
C-11
134.5
136.3
C-10
184.5
181.9
CH-1
8.55
8.67
157.9
157.4
NH
12.07
12.22
C-3
56.5
56.5
CH3-4,5,6
1.52
1.52
29.6
29.6
The chemical shifts
are reported
starting by the ring (C7 to C10), followed by the alkene side chain.
The chemical shifts
are reported
starting by the ring (C7 to C10), followed by the alkene side chain.For 28a as for 29a, a cyclohexanetrione
scaffold was identified, with the 13C chemical shifts of
the three carbonyl functions at 176.8, 181.7, and 184.5 ppm for 28a and 177.5, 182.0, and 181.9 ppm for 28a’. These chemical shift values are not compatible with a putative
imino-enol form. In fact, an enol form of lawsone derivatives was
reported by Perez et al.[41] in 2007. These
authors reported the synthesis of 2-hydroxy-3-substituted naphthoquinones
obtained by the Heck reaction. The 13C NMR analyses of
their compounds exhibited chemical shifts for the 1,4-naphthoquinone
ring at 181–183 ppm, while the chemical shift of the carbon
atom bearing the hydroxyl group of the enol function was close to
160 ppm. In 2021, Olyaei et al.[42] reported
the synthesis of lawsone enaminone derivatives. They identified two
forms and attributed the alkene proton (1H NMR in DMSO)
of the enaminone at 9.2 and 9.3 ppm in DMSO. The authors did not report
any 13C data. These previous assignments unambiguously
demonstrate that compounds 28 and 29 are
the two cis/trans stereoisomers of an enaminone system
in a 1:1 ratio. For compound 28a/28a′ analyzed
here, we have also identified eight sp2-hybridized carbons δC from 108.6 to 157.9 ppm (for 28a) and from 109.3
to 157.4 ppm (for 28a′) and four sp3-hybridized
carbons δC from 29.6 to 56.5 for both forms.
Biological
Activities
The antiplasmodial activity of
the synthesized naphthoquinones was evaluated against the P. falciparum resistant strain F32-ART (Table ). For compounds showing
the best antimalarial activities, that is, the lowest IC50 values, their cytotoxicity was determined using Vero cells in order
to calculate their selectivity index (SI). The anti-infectious activities
of these compounds were also determined regarding two other pathogens,
namely, L. donovani and M. tuberculosis H37Rv. In addition, we have also
included in Table the evaluated activities of simple synthons (lawsone and iodolawsone)
along with the reference drugs atovaquone for P. falciparum, miltefosine for L. donovani, and
streptomycin for M. tuberculosis. We
have also included two antimalarial compounds from the literature:
MMV007571[43−45] and XCV[46] (Figure ), which are supposed to target
the cytochrome bc1 complex and/or the dihydroorotate
dehydrogenase (DHODH). The docking scores of all compounds regarding
these two mitochondrial targets were also evaluated, and the obtained
values are summarized in Table .
Table 5
Biological Evaluation of Naphthoquinones
Derivatives Reported Herein
compound
P. falciparum F32ARTIC50 (μM)
CC50a (μΜ)
against Vero cells
SI for Plasmodium parasites
L. donovanib IC50 (μM)
CC50c (μΜ)
against macrophages
SI for Leishmania parasites
M. tuberculosis H37Rv MIC (μΜ)
9
>10
>147
10
>10
>147
11
>10
>147
12
>10
>147
13
>10
>147
14
>10
42
15
>10
>147
16
>10
6.97
49.22
7.06
>147
17
2.5
86
34
9
18
11
50
4.5
9.41
54.21
5.76
39
19
6
>100
>16
>147
20
10
150
21
>10
20
22
3.6
36
10
17–34
23
>10
84
24
9
>100
>10
147
26
>10
84
27
4
24
6
18
28
1.4
>50
>36
3.50
>100
>28.5
62–124
29
1.9
>50
>26
6.33
46.48
7.34
>147
lawsone (1)
>10
6.27
44
7
>147
612
11.65
9.94d
iodo-lawsone
>50
>100
>2
13.2
56.23
4.3
>147
atovaquone (7)
0.001
6
6000
11.12
44
3.95
>147
streptomycin
0.25
miltefosine
1.46
22.59
15.47
Cytotoxicity was
evaluated against
Vero cells and expressed as CC50, and the corresponding
SI was relative to the ratio CC50/IC50 for P. falciparum.
The molecules were tested against
both LV9 Leishmania donovani axenic
amastigote forms and intramacrophage amastigote forms. The values
written on this table correspond to the activities on the second form.
Cytotoxicity was evaluated
against
macrophage RAW 264.7 cells and expressed as CC50. The SI
calculated corresponds to the ratio CC50/IC50 (intramacrophage amastigote forms).
SI evaluated on Hep GAA16 cells[12] IC50: inhibitory concentration 50%.
CC50: cytotoxic concentration 50%.
Figure 4
Structures of compounds also analyzed for comparison purposes.
Table 6
Computational Evaluation of Naphthoquinone
Derivatives as Mitochondrial Antimalarial Agents
compound
P. falciparum F32-ARTIC50 (μM)
P. falciparumbc1 docking score (kcal/mol)
P. falciparum DHODH docking score (kcal/mol)
9
>10
–9.7
–6.8
10
>10
–7.1
–7.7
11
>10
0.0
–7.0
12
>10
–9.2
–8.3
13
>10
0.0
–6.8
14
>10
–3.5
–6.9
15
>10
–3.5
–7.3
16
>10
–8.9
–9.4
17
2.5
–12.2
–11.0
18
11
–9.5
–8.2
19
6
–8.8
–9.0
20
10
–2.6
–6.3
21
>10
–10.0
–7.3
22
3.6
–6.4
–9.3
23
>10
–11.2
–11.1
24
9
–10.4
–9.2
26
>10
–10.5
–9.0
27
4
–10.8
–9.3
28
1.4
a = −7.2
a = −7.0
a′ = −7.0
a′ = −7.0
29
1.9
a = −7.5
a = −7.1
a′ = −7.3
a′ = −7.5
lawsone (1)
>10
–6.5
–6.6
6
11.65
–8.7
–8.7
iodo-lawsone
>50
–6.9
–7.1
atovaquone (7)
0.001
–12.1
–10.7
MMV007571
1
–8.6
–7.3
XCV
0.073
–11.2
Structures of compounds also analyzed for comparison purposes.Cytotoxicity was
evaluated against
Vero cells and expressed as CC50, and the corresponding
SI was relative to the ratio CC50/IC50 for P. falciparum.The molecules were tested against
both LV9 Leishmania donovani axenic
amastigote forms and intramacrophage amastigote forms. The values
written on this table correspond to the activities on the second form.Cytotoxicity was evaluated
against
macrophage RAW 264.7 cells and expressed as CC50. The SI
calculated corresponds to the ratio CC50/IC50 (intramacrophage amastigote forms).SI evaluated on Hep GAA16 cells[12] IC50: inhibitory concentration 50%.
CC50: cytotoxic concentration 50%.
Activities against P. falciparum
Depending on the isocyanide used for the synthesis, three
series can be distinguished bearing the tert-butyl, N-butyl, and N-cyclohexyl scaffolds. Concerning
the first series of compounds (9–17), all but one exhibited IC50 values above 10 μM
against P. falciparum. Conversely,
compound 17 had an IC50 value of 2.5 μM,
which, owing to its low cytotoxicity on Vero cells (86 μM),
confers to the molecule promising pharmacological characteristics.
Compound 17 possesses a (1-benzyl-piperidin-4-yl) substituent
at position 3 of the furan ring, while an analogous functionality
(4-chlorophenyl cyclohexyl) is present in atovaquone. Interestingly,
the N-Boc derivative 16 of compound 17 exhibits an IC50 value higher than 10 μM,
indicating that this protecting group results in loss of antiplasmodial
activity of the compound.Concerning the second naphthofuroquinones
series (18–22), bearing an N-cyclohexyl ring, four compounds out of five evaluated,
exhibited IC50 values in the range of 3.6–11 μΜ.
Compound 21, bearing a nitro-furan moiety at C3 of the
furan ring, exhibited an IC50 value higher than 10 μM,
similar to compounds 15 and 26 of the other
two families bearing the same substituent. Compound 22 exhibits the relatively better activity with an IC50 value
of 3.6 μM, compared to compounds 18, 19, and 20 (11, 6, and 10 μM, respectively). Nevertheless, 19 possessing a 3,4,5-trimethoxy phenyl group at C3 of the
furan ring has a higher SI value (>16 vs 34 for 17),
rendering this compound potentially the most interesting of this series,
from a pharmacological point of view.In the third series of
compounds (23–27) bearing the n-butyl group, compound 27 having the same
scaffold as compound 17 of
the previous series at C3 exhibits the best activity with an IC50 value of 4 μM. Nevertheless, its higher cytotoxicity
on Vero cells (24 μM) and consequently its weak SI (6) renders
compound 27 less interesting in comparison to 17. In this series, compound 24 having a 3,4,5-trimethoxy
phenyl group attached to C3 has an IC50 in the same range
as 27 but exhibits a much lower cytotoxicity value (>100
μM), rendering 25 potentially more interesting.Finally, concerning these series of naphthofuroquinones, it is
worth pointing out that many of our compounds are 2- to 5-fold more
active against P. falciparum than derivative 6 reported by Borgati et al.[12]Among the other evaluated compounds, while lawsone and iodolawsone
are weak to very weak inhibitors of P. falciparum in vitro, compounds 28 and 29 seem more
interesting. For both couples of compounds bearing either a N-tert-butyl or N-butyl
group, the IC50 values are 1.4 and 1.9 μM, respectively,
associated with a weak cytotoxicity (CC50 > 50 μM)
resulting in promising SIs (>36 and >26, respectively). Noteworthily,
it is the first time that enaminodione systems possessing also an
α,β-dione functionality were evaluated against P. falciparum. These compounds keep the 1,4-quinone
scaffold where the 2-hydroxy group of the atovaquone is modified to
a ketone group. The 2-hydroxy function of atovaquone plays a crucial
role in the interaction of the drug with the bc1 complex,
and this interaction is believed to be disturbed in atovaquone-resistant Plasmodium. Thus, modification at this position may
be an interesting starting point for further development of drugs
able to circumvent atovaquone resistance.
Activities against L. donovani
Naphthoquinones are well known
for their reported antileishmanial
activities. All our new compounds were first evaluated in vitro on
the axenic form of L. donovani, the
parasite responsible for visceral leishmaniasis in humans. Six of
them, namely, compounds 8, 16, 18, 28, 29, and iodolawsone, exhibited IC50 values lower than 10 μM, thus justifying an evaluation
on the L. donovani intramacrophage
amastigote model, which is closer to the pathological conditions.
Before carrying out this experiment, it was necessary to evaluate
the cytotoxicity of the compounds on the RAW 264.7 macrophage as the
host cell of the in vitro Leishmania model. All the
six compounds exhibited a reasonable cytotoxicity, with CC50 values higher than 40 μM (Table ). Compound 28 was less cytotoxic
with a CC50 value higher than 100 μM and consequently
a SI higher than 28. Then, the six compounds were eligible to be evaluated
on the intramacrophage amastigote model. Except iodolawsone (IC50 = 13.2 μM and SI = 4.3), the compounds exhibited IC50 values lower than 10 μM, and compounds 28 and 29 showed potent activities in close range of the
reference drug miltefosine. Notably, compound 28 showed
the best activity (IC50 = 3.5 μM) and the best SI
(> 28) that is better than that of miltefosine. Compared to lawsone,
the substitutions leading to 28 increased the antileishmanial
activity, while they significantly lowered the cytotoxicity. It should
be noted that 28 was also the most active and selective
compound against P. falciparum (IC50 = 1.4 μM and SI > 36).
Activities against M. tuberculosis H37Rv
Among all the compounds
evaluated, naphthofuroquinones 17, 22, and 27 exhibit interesting
activities against M. tuberculosis H37Rv
strain as well as against P. falciparum.Compound 17 (MIC = 9 μM) contains a benzyl-N-piperidine scaffold at position C3 of the furan ring and
a tert-butyl group attached to the nitrogen atom
at position C2. Compound 22 (MIC = 17–34 μM)
has the same benzyl-N-piperidine scaffold at C3,
while a N-cyclohexyl fragment is attached to C2.
Compound 27 (MIC = 18 μM) contains the same benzyl-N-piperidine scaffold at C3 and a N-butyl
function at C2. Naphthofuroquinone derivatives have been reported
by Ho-Yeon Song[47,48] to possess antituberculous activities.
In a recent article, the authors reported a series of benzonaphthofurandiones
obtained by direct intramolecular annulation of 2-chloro-3-phenoxynaphthoquinone
derivatives. The compound library was screened against the M. tuberculosis H37Rv strain and exhibited MIC values
in the range of 1.21–74.6 μM. BNF15, the most potent
of these compounds (Figure ), has an O–CH2–N–(CH3)2 substituent on the benzofuran fragment.
The authors reported that BNF15 was effective against all drug-sensitive
and drug-resistant M. tuberculosis isolates
tested and effectively killed intracellular M. tuberculosis and nontuberculous mycobacteria. The authors point out that (i)
the naphthofuran system could be a valuable fragment for compounds
targeting M. tuberculosis and (ii)
a fine tuning of the substitution pattern of the benzofuran frame
is crucial in order to enhance the activities by more than two log
units.
Figure 5
Compound BNF15, active against drug-resistant Mycobacterium
tuberculosis.
Compound BNF15, active against drug-resistant Mycobacterium
tuberculosis.
Computational Studies on Mitochondrial Targets of P. falciparum
In order to rationalize the
activities against P. falciparum, we
carried out molecular docking on the two targets expected to interact
with our compounds, that is, cytochrome bc1 (bc1)
and dihydroorotate dehydrogenase (DHODH). In fact, naphthoquinone
derivatives based on lawsone scaffold are known to interfere with
these two mitochondrial targets of P. falciparum.[26]The AlphaFold2-modeled structure
of the cytochrome enzyme (Q7HP03) was used to characterize the predicted
interactions between the compounds described and the P. falciparum target bc1.[49] All local quality measures of the predicted
protein structure were assessed as very good (Figure ), that is, per-residue confidence scores
(pLDDT) for the model were all very high near the binding site region
(central groove inside the α-helices, near Tyr 263) and had
low expected position error (dark green in Figure , in angstrom) for the predicted alignment
errors (PAE), that is, the values for expected position error at residue x when the predicted and true structures are aligned on
residue y.
Figure 6
Predicted per-residue confidence score (pLDDT)
for the AlphaFold2
model for Pf bc1. Binding site between the α-helices
shown in stick representation.
Figure 7
PAE for
residues in the AlphaFold2 model for Pf bc1.
Predicted per-residue confidence score (pLDDT)
for the AlphaFold2
model for Pf bc1. Binding site between the α-helices
shown in stick representation.PAE for
residues in the AlphaFold2 model for Pf bc1.For the DHODH target, the X-ray crystal structure
of PfDHODH (Protein Data Bank (PDB) 7l)[50] was
retrieved from the PDB having a resolution of 1.6 Å, adequate
for docking studies.[51] Hence, we docked
all naphthoquinone conjugates into the putative quinone binding tunnel
formed by the N-terminal domain to check if our compounds are susceptible
to interact with this target.The naphthofuroquinone derivatives
were docked into the binding
site of P. falciparumbc1 and also PfDHODH, and their predicted interactions
showed similar binding modes and interaction partners as compared
to the known inhibitors atovaquone, MMV007571, and XCV. Good docking
scores and binding poses were thus obtained (Table and Figure ). Compounds 17, 19, 22, 24, and 27 were predicted to
exhibit the best binding.
Figure 8
Binding modes and interactions for compounds
in the binding site
of Pf bc1.
Binding modes and interactions for compounds
in the binding site
of Pf bc1.The naphthoquinone ring was predicted to be able to bind in two
different binding modes: flipped and unflipped with respect to atovaquone,
and their interactions were also similar to those of another known
inhibitor, MMV007571. There is good enrichment of the top-ranked compounds
in Table with known
strong inhibitors appearing in the top-ranked compounds (in bold).Several compounds that were experimentally active had a stronger
predicted interaction with both Pf bc1 and DHODH
(Table ), such as 17, 27, and atovaquone.Among the enaminone
compounds 28 and 29, compound 29a’ (−7.45 kcal/mol) was predicted
to make good use of bridging water molecules in the binding site with
better predicted docking scores in the hydrated binding site, the
latter which may be advantageous for the binding energy and/or specificity
for some compounds (Figure ).[52,53] This is also the case for compound
XCV (−10.79 kcal/mol). Compound 11 binds mostly
through van der Waals contacts, both in bc1 and in
DHODH. Compound 22 probably acts through inhibition of P. falciparum DHODH as it has a stronger predicted
interaction with this enzyme than with P. falciparumbc1, as compared to controls. XCV in DHODH makes
strong use of interactions (hydrogen bonds and π–π
stacking) with His185, as well as hydrogen bonds to Arg285. Atovaquone
makes π–π stacking interactions with Phe171 in
DHODH and with Phe264 in bc1. MMV007571 also interacts
mainly through van der Waals contacts in bc1. Binding
to several P. falciparum targets at
the same time may lead to a compound possessing stronger inhibition
properties, as well as the potential to avoid resistance.
Figure 9
Binding modes
and interactions for compounds in the binding site
of P. falciparum DHODH. X-ray crystal
structure co-crystallized with inhibitor XCV in white, self-docked
in magenta (RMSD = 0.18 Å), and water molecules as spheres.
Binding modes
and interactions for compounds in the binding site
of P. falciparum DHODH. X-ray crystal
structure co-crystallized with inhibitor XCV in white, self-docked
in magenta (RMSD = 0.18 Å), and water molecules as spheres.The ADME properties for all compounds were predicted
with Swiss-ADME.[54] All compounds passed
Lipinski’s rule
of five for bioavailability, most of them also passing Ghose, Veber,
Egan, and Muegge filters or lead-likeness.[55] PAINS warnings[56] were present only for
the well-known quinone group.The good properties of predicted
binding, experimental inhibition,
interactions inside the hydrated and nonhydrated binding sites, multitarget
binding,[57] and favorable molecular properties
for further optimization of the compounds make them valuable features
to continue their development as possible leads for malaria therapeutics.
Conclusions
Based on lawsone, the first synthetic approach
of a three (or two)-component
reaction involving various aldehydes and three commercially available
isocyanides is reported. Under microwave irradiation, this approach
afforded a series of naphthoquinones in moderate to good yields. In
addition, with the less-reactive aldehydes, we obtained for the first
time two naphtho-enaminodiones quinones derived from the direct condensation
of lawsone with the corresponding isocyanide derivative.Among
the naphthofuroquinone series, compound 17 exhibited
the best activity against P. falciparum, which is more than 4 times better than the previous optimal reported
in the naphthofuroquinone series.[12] In
addition, compound 17 exhibited the best SI = 34 among
the synthesized naphthofuroquinones. Considering the molecular docking
of 17 on the potential targets Pfbc1
and PfDHODH, this compound showed very favorable interactions, thus
making it a valuable starting point for further development.Noteworthily, the two naphtho-enaminodione quinone compounds 28 and 29 also exhibited antimalarial activities
in the micromolar range. Compound 28 exhibited the best
activity against P. falciparum (IC50 = 1.4 μM) and the highest SI (> 36), while molecular
docking predicted favorable interactions with both targets. Then,
naphtho-enaminodiones quinones should be thus considered as interesting
compounds in terms of novelty in synthesis but also due to their biological
activities.Against M. tuberculosis, the same
compound 17 exhibited the best activity (MIC = 9 μM).
It is about 7-fold less active than the known naphthofuroquinone derivative
BNF15. A fine tuning of the substitution pattern should be important
in order to increase the activity of compound 17 against M. tuberculosis.Finally,
compounds 28 and 29 have potent
activities against L. donovani that
are close to that of the reference drug miltefosine. More importantly,
compound 28 showed the best activity (IC50 = 3.5 μM) and the best SI (> 28) that is higher than that
of miltefosine. Compound 28 is one of the rare compounds
which are active against both P. falciparum and L. donovani.In addition,
compound 28 is worth to be evaluated
in vivo against the L. donovani/BALB/c
mouse model.
Methods
Experimental Section
Our reagents and solvents were
purchased from Sigma-Aldrich, TCI, Alfa Aesar, and Fluorochem and
used as received without any further purification. Microwave irradiation
reactions were performed in a CEM Discover SP Microwave model 909150/SN:
DC 9208 apparatus. A monowave 50 thermic reactor was obtained from
Anton Paar. TLC was performed on silica gel 60 F254 plates
(Merck). The compounds and the reaction mixtures were visualized on
the TLC plates by irradiation with UV light. For flash column chromatography,
a PuriFlash XS520Plus system was used in combination with PF-30SIHP-JP-F0040
columns. When PE/A* is mentioned in the Experimental
Section, A* stands for the DCM/AcOEt (8:2) system of solvents. 1H and 13C NMR spectra for the reported compounds
were recorded on a Bruker Avance I 300 MHz (300 MHz for 1H and 75 MHz for 13C), on a Bruker Avance III Nanobay
400 MHz (400.0 MHz for 1H and 101 MHz for 13C), and on a Bruker Avance 600 MHz (600 MHz for 1H and
151 MHz for 13C) equipped with a 5 mm triple-resonance
inverse Z-gradient probe (TBI 1H, 31P, BB).
Chemical shifts (δ) and coupling constants (J) are expressed in ppm and Hz, respectively. The NMR experiments
were performed in CDCl3 or DMSO-d6 and referenced to the solvent signal. High-resolution mass
spectrometry (HRMS) analyses were carried out on an XevoG2QTof (Waters)
using electrospray ionization. Melting points were determined using
a Stuart SMP3 apparatus, and the obtained values are not corrected.
UV–visible spectra of selected compounds (CH2Cl2 solutions at 0.01 or 0.02 mg/mL) were recorded from 200 to
800 nm on Cary 3500 spectrophotometers (see the Supporting Information). Infrared spectra of selected compounds
were measured using the Perkin Elmer Frontier MIR/FIR spectrometer
and reported in the Experimental Section. The samples were used as such (powders).
Syntheses
Microwave
Irradiation Procedure
A suspension of lawsone
(0.4 mmol) in 1,2-dichloroethane (7 mL) was charged with the corresponding
aldehyde (0.4 mmol), isocyanide (0.4 mmol), and 10% mol of EDDA. The
reaction mixture was then stirred and irradiated (200 W, 7 bar) for
1–4 h at 160 °C until completion of the reaction (monitored
by TLC). The solvent was removed under vacuum evaporation, and the
obtained residue was diluted in DCM. The organic phase was washed
with 5% aqueous solution NaHCO3, water, and brine. The
organic layer was dried over Na2SO4, filtered,
and evaporated under reduced pressure. The obtained crude products
were then purified by FCC, affording the corresponding products in
moderate to good yields.
Conventional Procedure
A mixture
of lawsone (0.6 mmol)
in toluene (10 mL) was charged with the corresponding aldehyde (0.6
mmol), isocyanide (0.6 mmol), and 10% mol of EDDA. The obtained suspension
was stirred and set under reflux overnight until completion of the
reaction (monitored by TLC). The solvent was removed under vacuum
evaporation, and the obtained residue was diluted in DCM. The organic
phase was then washed with 5% aqueous solution NaHCO3,
water, and brine. The organic layer was dried over Na2SO4, filtered, and evaporated under reduced pressure. The obtained
crude products were then purified by FCC, affording the corresponding
products in moderate to good yields.
The compound was synthesized
by following the above-mentioned alternative procedure (0.5 mmol lawsone).
Reaction time 24 h. Aldehyde conversion 22%. The crude product was
purified by FCC with PE/A* (5:5) and then with PE/A* (4:6) to yield 129 mg (50%) of the product as a yellow solid
(mp 159 °C decomposition). R (PE/A 4:6) = 0.20. 1H NMR (400 MHz, CDCl3): δ 11.72 (s, 1H), 11.58 (s, 1H), 8.54 (d, J = 14.5 Hz, 1H), 8.42 (d, J = 14.5 Hz, 1H), 8.20
(d, J = 7.7 Hz, 1H), 8.15–8.08 (m, 1H), 7.75
(t, J = 7.5 Hz, 1H), 7.67 (t, J =
7.5 Hz, 1H), 3.58 (dq, J = 13.8, 6.7 Hz, 2H), 1.78–1.64
(m, 2H), 1.42 (dq, J = 14.7, 7.4 Hz, 2H), 0.95 (td, J = 7.2, 2.6 Hz, 3H). 13C NMR (75 MHz, CDCl3): δ 181.6 (C=O), 181.5 (C=O), 177.8 (C=O),
162.0 (CH), 161.5 (CH), 135.1 (CH), 134.9 (CH), 134.1 (C), 133.2 (CH),
133.1 (CH), 132.8 (C), 128.1 (CH), 127.7 (CH), 127.2 (CH), 126.7 (CH),
109.8 (C), 109.2 (C), 51.1 (CH2), 32.1 (CH2),
19.7 (CH2), 13.6 (CH3). HRMS calcd for C15H16NO3+ [M + H]+ = 258.1130; found, 258.1134.
X-ray Analysis of Compound 9
Data were
collected at low temperature (100 K) on a Bruker APEX II diffractometer
using a micro-focus-sealed X-ray tube, Mo-Kα radiation (λ
= 0.71073 Å), and equipped with an Oxford Cryosystems Cryostream
Cooler Device. The structures have been solved by Direct Methods using
a SHELXS-97[58] and refined by means of least-squares
procedures on F2 with the aid of the program
SHELXL2016[58] included in the software package
WinGX version 1.63.[59] The Atomic Scattering
Factors were taken from International tables for X-ray crystallography.[60] All hydrogen atoms were placed geometrically,
except for H1 carried by the N1 atom which was located by Fourier
difference maps. They were refined using an overlap model.All
nonhydrogen atoms were anisotropically refined, and in the last cycles
of refinement, a weighting scheme was used, where weights are calculated
from the following formula: w = 1/[σ2(Fo2) + (aP)2 + bP], where P = (Fo2 + 2Fc2)/3.Drawing of molecules were performed with
the program ORTEP32[61] with 30% probability
displacement ellipsoids
for nonhydrogen atoms.
2D NMR Analysis of Compounds 9, 28, and 29
In the present study, 1H and 13C NMR spectroscopies were used for the
characterization
of compounds 9, 28, and 29.
NMR samples were prepared by dissolving 10–20 mg of each compound
in 600 μL of CDCl3. All spectra were recorded on
a Bruker Avance 600 spectrometer equipped with a 5 mm triple-resonance
inverse Z-gradient probe (TBI 1H, 31P, and BB).
All chemical shifts for 1H and 13C are relative
to TMS using 1H (residual) or 13C chemical shifts
of the solvent as a secondary standard. Compound 9 was
also analyzed by 2D-NMR experiment at a temperature of 298 K. All
the 1H and 13C signals were assigned on the
basis of chemical shifts, spin–spin coupling constants, splitting
patterns, and signal intensities and by using 1H–1H COSY45, 1H–13C HSQC, and 1H–13C HMBC experiments. Gradient-enhanced 1H COSY45 was performed by including eight scans for per increment. 1H–13C correlation spectra using a gradient-enhanced
HSQC sequence (delay was optimized for 1JCH of
145 Hz) was obtained with 16 scans per increment. Gradient-enhanced
HMBC experiment was performed, allowing 62.5 ms for long-range coupling
evolution (64 scans were accumulated). Typically, 2048 t2 data points
were collected for 256 t1 increments.
Computational and Docking
Experiments
The AlphaFold2-predicted
structure of the protein Pf bc1 was downloaded from
the European Bioinformatics Institute website[62] and preprocessed with PrepWizard [Schrödinger LLC. 2021]
to add hydrogens, minimize the structure, and resolve the ionization
states and clashes. The binding site was specified in analogy to the
structure 4pd4 from the PDB.[63]The
X-ray crystal protein structure (7l01, resolution = 1.60 Å) for Pf DHOHD was downloaded from the PDB[63] and preprocessed with PrepWiz [Schrödinger LLC.
2021] both with and without water molecules. Co-factors (FMN—flavin
mononucleotide and ORO—orotic acid) were kept in the structure.
The co-crystallized ligand DSM782 (XCV, N-(1-(5-cyano-1H-pyrazol-3-yl)ethyl)-3-methyl-4-(1-(6-(trifluoromethyl)pyridin-3-yl)cyclopropyl)-1H-pyrrole-2-carboxamide) was also extracted and re-docked
(self-dock) for comparison and as a control.Chemical compounds
were imported as SMILES or drawn, energy-minimized
with Maestro [Schrödinger LLC. 2021], and then processed with
LigPrep [Schrödinger LLC. 2021] for assigning tautomers and
ionization states around pH 7 ± 2. The synthesized compounds
and known inhibitors were used as controls.Molecules were docked
into the binding site of bc1 using Glide
XP [Schrödinger LLC. 2021], including aromatic hydrogens as
donors, halogen acceptors, and other settings.[64,65] The input file is provided below:FORCEFIELD OPLS_2005GRID_CENTER −1.9900714545454548, 7.686090977272727, 9.326240250000001GRIDFILE glide-grid_AF-Q7_2.zipHBOND_ACCEP_HALO TrueHBOND_DONOR_AROMH TrueINNERBOX 12, 12, 12OUTERBOX 29.189783718905936,
29.189783718905936, 29.189783718905936RECEP_FILE glide-grid_AF-Q7_2.maegzMolecules were also docked into the binding site of DHODH, as defined
by the co-crystallized ligand inhibitor XCV (N-[(1R)-1-(5-cyano-1H-pyrazol-3-yl)ethyl]-3-methyl-4-{1-[6-(trifluoromethyl)pyridin-3-yl]cyclopropyl}-1H-pyrrole-2-carboxamide), in several runs using both the
crystallographic waters as well as without them.[52,53] The input file is given below:EPIK_PENALTIES FalseEXPANDED_SAMPLING TrueFORCEFIELD OPLS_2005GRIDFILE
glide-grid_4_nowats.zipHBOND_ACCEP_HALO TrueHBOND_DONOR_AROMH
TrueINCLUDE_INPUT_RINGS TrueLIGANDFILE todock.sdfPOSTDOCK_XP_DELE 0.5PRECISION XPWRITE_XP_DESC False
Antiplasmodial Activity and Cytotoxicity
The different
compounds were evaluated in vitro for their antiplasmodial activity
against the P. falciparum resistant
strain F32-ART, selected after 144 intermittent and increasing doses
of artemisinin, according to our published procedures.[66,67] Briefly, the antimalarial effect was first determined by Sybr Green
at two doses (1 and 10 μM), in triplicate, and on two independent
experiments. For compounds inducing around or more than 50% parasite
growth inhibition at 10 μM, a new chemosensitivity assay using
four doses, each one tested in triplicate, was then performed to determine
their exact IC50 values. For the best antimalarial compounds,
their cytotoxicity was obtained using Vero cells, leading to the calculation
of their SI as the ratio cytotoxicity/activity.
Antituberculous
Activity: MIC Determination by Resazurin Reduction
Microplate Assay
To determine the in vitro activity of the
compounds (MIC90) in M. tuberculosis H37Rv, the resazurin reduction microplate assay was performed as
previously described.[68] Briefly, serial
2-fold dilutions (starting from 64 μg/mL (for H37Rv)) of each
drug were prepared in 96-well black plates (Fluoronunc, Thermo Fisher,
Waltham, MA, USA) in 100 μL of Middlebrook 7H9 medium, without
the addition of Tween 80. Then, log-phase cultures were diluted (OD600 = 0.0005) and added in a 96-well black plate. Growth controls
containing no compound and sterile controls without inoculum were
also included. After 7 days of incubation at 37 °C, 10 μL
of resazurin (0.025% w/v) was added to each well, and bacterial viability
was assessed after a further 24 h of incubation using a Fluoroskan
Microplate Fluorometer (Thermo Fisher Scientific, Waltham, MA, USA;
excitation = 544 nm and emission = 590 nm). Bacterial viability was
calculated as a percentage of resazurin turnover in the absence of
compound. Streptomycin was used as a positive control. Results were
expressed as the average of at least three independent replicates.
Antileishmanial Activity and Cytotoxicity on RAW 264.7 Macrophages
Cell
Lines
The mouse monocyte/macrophage cell line
RAW 264.7 and L. donovani (MHOM/ET/67/HU3,
also called LV9) promastigotes and axenic amastigotes were maintained
according to the protocols described in Pomel et al., 2021.[69]
Evaluation of Compound Cytotoxicity on RAW
264.7 Macrophages
Cytotoxicity was evaluated on RAW 264.7
macrophages using the resazurin
method as detailed in Pomel et al., 2021.[69]
In Vitro Antileishmanial Evaluation on L. donovani Axenic Amastigotes
This evaluation was performed using
the SYBR Green method as previously described.[68] IC50 values were calculated using the ICEstimator
version 1.2 software (http://www.antimalarial-icestimator.net/runregression1.2.htm). Miltefosine was used as the reference drug.
In Vitro
Antileishmanial Evaluation on Intramacrophage Amastigotes
Determination of cytotoxicity, as presented above, was used to
select the highest drug concentrations that could be studied on the L. donovani intramacrophage amastigote model using
RAW 264.7 cells. Macrophages were infected with L.
donovani axenic amastigotes according to a ratio of
10 parasites per macrophage. In these conditions, the percentage of
infected macrophages was around 80%, and the mean number of amastigotes
per infected macrophage was 4 to 5 in the untreated controls. The
in vitro treatment was applied 24 h post-infection, and the
treatment duration was 48 h. The results of the effect of the
compounds are given as percentage reduction of parasite growth, measured
using the SYBR Green incorporation method. The activity of the compounds
is expressed as IC50, calculated using the ICEstimator
version 1.2 software (Pomel et al., 2021). Miltefosine was used as
the reference drug.
Authors: Birgit Viira; Anastasia Selyutina; Alfonso T García-Sosa; Maarit Karonen; Jari Sinkkonen; Andres Merits; Uko Maran Journal: Bioorg Med Chem Date: 2016-04-08 Impact factor: 3.641
Authors: Xiaoyi Deng; David Matthews; Pradipsinh K Rathod; Margaret A Phillips Journal: Acta Crystallogr F Struct Biol Commun Date: 2015-04-21 Impact factor: 1.056
Authors: Benjamin K Dickerman; Brendan Elsworth; Simon A Cobbold; Catherine Q Nie; Malcolm J McConville; Brendan S Crabb; Paul R Gilson Journal: Sci Rep Date: 2016-11-22 Impact factor: 4.379
Authors: Maged Henary; Carl Kananda; Laura Rotolo; Brian Savino; Eric A Owens; Giancarlo Cravotto Journal: RSC Adv Date: 2020-04-07 Impact factor: 3.361