Bijo Mathew1, Jong Min Oh2, Mohamed A Abdelgawad3, Ahmed Khames4, Mohammed M Ghoneim5, Sunil Kumar1, Lekshmi R Nath6, Sachithra Thazhathuveedu Sudevan1, Della Grace Thomas Parambi3, Clement Agoni7, Mahmoud E S Soliman7, Hoon Kim2. 1. Department of Pharmaceutical Chemistry, Amrita School of Pharmacy, Amrita Vishwa Vidyapeetham, AIMS Health Sciences Campus, Kochi 682041, India. 2. Department of Pharmacy, and Research Institute of Life Pharmaceutical Sciences, Sunchon National University, Suncheon 57922, Republic of Korea. 3. Department of Pharmaceutical Chemistry, College of Pharmacy, Jouf University, Sakaka, Al Jouf 72341, Saudi Arabia. 4. Department of Pharmaceutics and Industrial Pharmacy, College of Pharmacy, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia. 5. Department of Pharmacy Practice, Faculty of Pharmacy, AlMaarefa University, Ad Diriyah 13713, Saudi Arabia. 6. Department of Pharmacognosy, Amrita School of Pharmacy, Amrita Vishwa Vidyapeetham, AIMS Health Sciences Campus, Kochi 682041, India. 7. Molecular Bio-Computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Westville Campus, Durban 4001, South Africa.
Abstract
Fifteen multiconjugated dienones (MK1-MK15) were synthesized and evaluated to determine their inhibitory activities against monoamine oxidases (MAOs) A and B. All derivatives were found to be potent and highly selective MAO-B inhibitors. Compound MK6, with an IC50 value of 2.82 nM, most effectively inhibited MAO-B, like MK12 (IC50 = 3.22 nM), followed by MK5, MK13, and MK14 (IC50 = 4.02, 4.24, and 4.89 nM, respectively). The selectivity index values of MK6 and MK12 for MAO-B over MAO-A were 7361.5 and 1780.5, respectively. Compounds MK6 and MK12 were competitive reversible inhibitors of MAO-B, with K i values of 1.10 ± 0.20 and 3.0 ± 0.27 nM, respectively. Cytotoxic studies showed that MK5, MK6, MK12, and MK14 exhibited low toxicities on Vero cells, with IC50 values of 218.4, 149.1, 99.96, and 162.3 μg/mL, respectively, which were much higher than those for their effective nanomolar-level concentrations. Also, MK5, MK6, MK12, and MK14 decreased cell damage in H2O2-induced cells via a significant scavenging effect of reactive oxygen species. Molecular modeling was performed to rationalize the potential inhibitory activities of MK5, MK6, MK12, and MK14 toward MAO-B and their possible binding mechanisms, showing high-affinity binding pocket interactions and conformation perturbations of the compounds with MAO-B, which were interpreted as the conformational dynamics of MAO-B. This study concluded that all the compounds tested were more potent MAO-B inhibitors than the reference drugs, and leading compounds could be further explored for their effectiveness in various kinds of neurodegenerative disorders.
Fifteen multiconjugated dienones (MK1-MK15) were synthesized and evaluated to determine their inhibitory activities against monoamine oxidases (MAOs) A and B. All derivatives were found to be potent and highly selective MAO-B inhibitors. Compound MK6, with an IC50 value of 2.82 nM, most effectively inhibited MAO-B, like MK12 (IC50 = 3.22 nM), followed by MK5, MK13, and MK14 (IC50 = 4.02, 4.24, and 4.89 nM, respectively). The selectivity index values of MK6 and MK12 for MAO-B over MAO-A were 7361.5 and 1780.5, respectively. Compounds MK6 and MK12 were competitive reversible inhibitors of MAO-B, with K i values of 1.10 ± 0.20 and 3.0 ± 0.27 nM, respectively. Cytotoxic studies showed that MK5, MK6, MK12, and MK14 exhibited low toxicities on Vero cells, with IC50 values of 218.4, 149.1, 99.96, and 162.3 μg/mL, respectively, which were much higher than those for their effective nanomolar-level concentrations. Also, MK5, MK6, MK12, and MK14 decreased cell damage in H2O2-induced cells via a significant scavenging effect of reactive oxygen species. Molecular modeling was performed to rationalize the potential inhibitory activities of MK5, MK6, MK12, and MK14 toward MAO-B and their possible binding mechanisms, showing high-affinity binding pocket interactions and conformation perturbations of the compounds with MAO-B, which were interpreted as the conformational dynamics of MAO-B. This study concluded that all the compounds tested were more potent MAO-B inhibitors than the reference drugs, and leading compounds could be further explored for their effectiveness in various kinds of neurodegenerative disorders.
Monoamine oxidases
(MAOs) are the prime metabolizing enzymes of
various biogenic amines via oxidative deamination.[1] The alteration of biogenic amine concentrations in the
brain by MAO directly correlates with several neurological disorders,
such as Parkinson’s disease (PD) and Alzheimer’s disease
(AD).[2] This oxidative degradation generates
toxic byproducts, such as hydrogen peroxide, reactive oxygen species
(ROS), and ammonia, which can trigger oxidative stress with mitochondrial
dysfunction in neural cells.[3,4] High levels of MAO-B
have been observed in the substantia nigra of PD patients, and progressively
reversible and highly selective MAO-B inhibitors have proved efficient
for relieving the symptoms of PD patients.[5]Chalcones are simple organic compounds with enone-based linkers
between phenyl and hetero nuclei.[6] Numerous
structural manipulations have been applied to the chalcone scaffold
as a selective MAO-B inhibitor, changing the hetero nucleus by placing
various electron-withdrawing and/or electron-donating groups on the
two aromatic/heteroaromatic rings.[7] These
studies have demonstrated that factors such as the length, electron
delocalization, and hydrophobicity of rings around linkers play a
crucial role in the development of MAO inhibitors.[8−15]In 2013, Desideri et al. reported that an extended conjugation
in the chalcone framework could exhibit remarkable MAO-B inhibition;
compounds (2E,4E)-5-(4-chlorophenyl)-1-(2-hydroxy-4-methoxyphenyl)penta-2,4-dien-1-one
and (2E,4E)-5-(4-chlorophenyl)-1-(2,4-dihydroxyphenyl)penta-2,4-dien-1-one
were the most potent human MAO-B inhibitors, with IC50 values
of 4.51 and 11.35 nM, respectively.[16] A
previous research report from our laboratory showed that extended
conjugations in chalcone linkers exhibited MAO-B inhibitory activity;
an unsaturation extension on the furan-based chalcone (F1) improved the MAO-B inhibitory activity (Ki = 0.0041 μM) and, to a greater extent, the selectivity
index (SI = 172.4).[17] The extended conjugation
was also appraised in 1-[4-(morpholin-4-yl)phenyl]-5-phenylpenta-2,4-dien-1-one
(MO10) as a potent selective MAO-B inhibitor (IC50 = 0.044 μM; Ki = 0.0080
± 0.003 μM), with an SI value of 366.13.[18] We hypothesized that the electronic feature of the carbonyl
group in a linker could be enhanced by two carbon–carbon double
bond units. A recent study provided evidence that the presence of
halogens on various MAO-B inhibitors significantly impacts the energetic
stability of the inhibitor-binding cavities of these enzymes.[19]Evidence from the recent design of MAO-B
inhibitors showed the
importance of halogens and extended conjugation in chalcone scaffolds.
We obtained the framework of the title compound, which involved (1)
extended conjugation of the three carbon linkers by adding an olefinic
linkage to improve the electrophilic nature of the spacers, (2) addition
of halogens at the para position of the chalcone ring A, and (3) introduction
of electron-donating methoxyl and electron-withdrawing nitro groups
at the para position of ring B of conjugated dienones (Figure ).
Figure 1
Design strategy for conjugated
dienones as potent MAO-B inhibitors.
Design strategy for conjugated
dienones as potent MAO-B inhibitors.On the other hand, AD is associated with a decrease in neurotransmitters,
specifically acetylcholine (ACh), and with an increase in acetylcholinesterase
(AChE) and/or butyrylcholinesterase (BChE).[20] In addition, β-secretase (β-site amyloid precursor protein-cleaving
enzyme 1, BACE1) has been considered a target for AD treatment because
it induces AD through the production of amyloid-β peptides.[21] Recently, multifunctional agents targeting MAO-A,
MAO-B, AChE, BChE, and BACE1 have been studied for the effective treatment
of AD.[22,23]The current study synthesized a series
of conjugated dienones (MK1–MK15) and investigated
their in vitro MAO-A and
MAO-B inhibitory profiles, including AChE, BChE, and BACE1 enzymes.
The lead molecules were further subjected to kinetics, reversibility
studies, assessment of cytotoxicity on normal cell lines, ROS assay,
and molecular dynamics (MD) simulation.
Results and Discussion
Chemistry
The multiconjugated dienones were synthesized
by a pyrrolidine-catalyzed reaction between various substituted cinnamaldehyde
derivatives and halogenated acetophenones (Scheme ). All final derivatives were characterized
using 1H NMR, 13C NMR, and mass spectrometry
(see the Supporting Information).
Scheme 1
Synthesis
of Multiconjugated Ketones (MK1–MK15): (a) Pyrrolidine
and (b) Ethanol
Biochemistry
MAO Inhibition
Studies
All compounds more effectively
inhibited MAO-B than MAO-A and had strong inhibitory activities against
MAO-B, with residual activities of <50% at 1 μM (Table ). In general, (2E,4E)-1,5-diphenylpenta-2,4-dien-1-one
derivatives (MK1–MK15) showed strong inhibitory
activities against MAO-B. Compound MK6 most potently
inhibited MAO-B with an IC50 value of 2.82 nM, followed
by MK12 (IC50 = 3.22 nM) (Table ). The −Br atom at the
para position of MK6 (a parent of the second subseries)
increased the MAO-B inhibitory activity compared to the −Cl
atom at the para position of MK1 (IC50 = 11.17
nM—a parent of the first subseries). In the first subseries
containing the −Cl atom at the para position, an NO2 atom at the para position of MK5 (IC50 =
4.02 nM) increased the MAO-B inhibitory activity compared to the parent MK1. In the second subseries containing the −Br atom
at the para position, all derivatives substituted by other groups
decreased the inhibitory activity against MAO-B. In the third subseries
containing the −F atom at the para position, the methoxy group
of MK12 had the most effective inhibitory activity against
MAO-B (IC50 = 3.22 nM), followed by the nitro group of MK13 (IC50 = 4.24 nM). However, when the −F
atom of MK13 was replaced by the −H atom of MK15 (IC50 = 12.4 nM), the inhibitory activity
decreased. Moreover, the methoxy group and −F atom of MK12 effectively inhibited both MAO-B and MAO-A (Table ). MK6 was selective for MAO-B, with an SI value of 7361.5 over MAO-A (Table ). Multitarget analyses
showed that all compounds weakly inhibited AChE, BChE, and β-secretase
(BACE1) at 10 μM (Table ).
Table 1
Inhibition of MAO-A, MAO-B, AChE,
BChE, and BACE1 by the MK seriesa
residual
activity (%)
IC50 (μM, nM)
MAO-A
MAO-B
AChE
BChE
BACE1
MAO-A
MAO-B
compounds
10 μM
1 μM
10 μM
10 μM
10 μM
(μM)
(nM)
SIb
MK1
83.0 ± 2.71
2.17 ± 1.49
84.0 ± 0.23
94.1 ± 4.72
77.1 ± 0.0087
16.7 ± 0.40
11.17 ± 1.65
1498.2
MK2
97.8 ± 0.16
14.3 ± 0.43
82.4 ± 0.98
95.4 ± 2.86
81.4 ± 0.97
23.5 ± 0.86
101.15 ± 1.63
232.4
MK3
73.6 ± 2.83
5.81 ± 0.27
81.4 ± 1.06
83.2 ± 0.46
84.8 ± 0.53
16.9 ± 0.50
13.07 ± 0.81
1290.2
MK4
77.0 ± 6.20
1.61 ± 2.28
90.7 ± 0.57
94.9 ± 7.16
72.5 ± 0.90
15.6 ± 0.63
9.67 ± 0.47
1609.8
MK5
72.9 ± 0.44
5.81 ± 0.27
74.1 ± 1.42
71.6 ± 1.05
73.3 ± 2.51
19.2 ± 0.94
4.02 ± 0.13
4783.0
MK6
84.0 ± 1.22
–1.69 ± 1.02
74.7 ± 0.66
97.5 ± 3.58
74.1 ± 1.48
20.6 ± 1.02
2.82 ± 0.39
7361.5
MK7
85.4 ± 1.99
8.90 ± 0.15
75.8 ± 0.73
99.2 ± 1.09
84.5 ± 0.91
21.7 ± 2.01
15.50 ± 0.06
1400.0
MK8
69.6 ± 2.03
7.85 ± 0.22
86.1 ± 5.10
80.7 ± 7.70
80.9 ± 1.67
24.1 ± 0.038
16.69 ± 0.62
1443.8
MK9
74.9 ± 1.29
2.60 ± 0.57
57.7 ± 0.00
87.4 ± 6.95
65.0 ± 1.61
22.5 ± 0.047
13.48 ± 1.06
1671.7
MK10
69.1 ± 1.32
3.40 ± 0.14
84.0 ± 0.73
87.9 ± 1.86
69.4 ± 1.03
19.5 ± 0.59
40.62 ± 0.87
481.0
MK11
76.9 ± 3.02
3.45 ± 0.21
77.8 ± 2.19
89.2 ± 6.58
71.4 ± 1.03
27.5 ± 0.82
19.75 ± 0.53
1394.7
MK12
40.9 ± 1.24
–2.71 ± 0.73
93.3 ± 0.73
94.6 ± 1.11
81.4 ± 1.03
5.70 ± 0.72
3.22 ± 0.04
1780.5
MK13
71.5 ± 8.06
1.95 ± 0.067
91.5 ± 0.80
97.8 ± 1.55
76.7 ± 0.64
23.7 ± 0.17
4.24 ± 0.13
5635.2
MK14
61.4 ± 6.27
–0.95 ± 1.35
84.7 ± 5.62
89.0 ± 1.55
77.3 ± 1.23
18.9 ± 1.02
4.89 ± 0.17
3871.5
MK15
54.4 ± 5.37
–3.52 ± 0.41
68.8 ± 0.80
80.2 ± 1.55
89.5 ± 0.92
12.7 ± 0.57
12.40 ± 7.92
1025.8
toloxatone
1.08 ± 0.025
lazabemide
110.00 ± 16.00
clorgyline
0.0070 ± 0.00070
pargyline
140.0 ± 5.90
Results are the means ± standard
errors from duplicate or triplicate experiments.
SI values are expressed for MAO-B
compared with that for MAO-A. For tacrine (a reference compound for
AChE and BChE), IC50 was confirmed by values of 270.0 ±
19.0 and 60.0 ± 2.2 nM, respectively. For donepezil (a reference
compound for AChE and BChE), IC50 was confirmed by values
of 9.5 ± 1.9 and 180.0 ± 3.8 nM, respectively. For quercetin
(a reference compound for BACE1), IC50 was confirmed by
the value of 13.4 ± 0.035 μM. For BACE inhibitor IV (a
reference compound for BACE1), IC50 was confirmed by the
value of 0.44 ± 0.064 μM.
Results are the means ± standard
errors from duplicate or triplicate experiments.SI values are expressed for MAO-B
compared with that for MAO-A. For tacrine (a reference compound for
AChE and BChE), IC50 was confirmed by values of 270.0 ±
19.0 and 60.0 ± 2.2 nM, respectively. For donepezil (a reference
compound for AChE and BChE), IC50 was confirmed by values
of 9.5 ± 1.9 and 180.0 ± 3.8 nM, respectively. For quercetin
(a reference compound for BACE1), IC50 was confirmed by
the value of 13.4 ± 0.035 μM. For BACE inhibitor IV (a
reference compound for BACE1), IC50 was confirmed by the
value of 0.44 ± 0.064 μM.Interestingly, all tested compounds showed significant
MAO-B inhibitory
activities compared to the reference drugs. The lead compounds MK6 and MK12 had potent MAO-B inhibitory activities
(IC50 = 2.82 ± 0.39 and 3.22 ± 0.04 nM, respectively),
which were 39 and 34 times more potent, respectively, than that of
the reference reversible MAO-B drug lazabemide. These lead molecules
also showed 50 and 44 times higher inhibitory activities, respectively,
than the reference irreversible MAO-B inhibitor pargyline.
Kinetic
Study
Based on kinetic studies of MK6 and MK12 for MAO-B, Lineweaver–Burk plots showed
that the lines for MK6 and MK12 met at a
point on the y-axis (Figure A,C), and their secondary plots had Ki values of 1.10 ± 0.20 and 3.00 ±
0.27 nM, respectively (Figure B,D). These results suggest that MK6 and MK12 are competitive inhibitors that bind at the active site
of MAO-B.
Figure 2
Lineweaver–Burk plots for MAO-B inhibition by MK6 and MK12 (A,C) and their respective secondary plots
(B,D) of slopes vs inhibitor concentrations.
Figure 3
Recoveries
of MAO-B inhibition by MK6 and MK12 using
dialysis experiments.
Lineweaver–Burk plots for MAO-B inhibition by MK6 and MK12 (A,C) and their respective secondary plots
(B,D) of slopes vs inhibitor concentrations.Recoveries
of MAO-B inhibition by MK6 and MK12 using
dialysis experiments.
Reversibility Studies
In the experiments, the concentration
of MK6 or MK12 was 6.0 nM and that of lazabemide
(a reference reversible inhibitor) and pargyline (a reference irreversible
inhibitor) was 0.22 and 0.28 μM, respectively. The relative
activities for undialyzed (AU) and dialyzed
(AD) samples were compared to determine
their reversibility patterns. The inhibition of MAO-B by MK6 and MK12 was recovered from 34.9% (AU) to 79.2% (AD) and from
33.4 to 74.7%, respectively (Figure ). These recovery values were similar to those of lazabemide,
a reversible reference inhibitor against MAO-B (i.e., from 34.9 to
83.2%), and could be distinguished from pargyline, an irreversible
reference inhibitor against MAO-B (i.e., from 35.9 to 35.9%). These
results indicated that MK6 and MK12 were
reversible inhibitors of MAO-B.
Cytotoxicity Studies of
Vero Cells
We evaluated the
biological safety of the effective compounds, such as MK5, MK6, MK12, and MK14, on
a normal epithelial cell line from the kidney of an African green
monkey (Vero cells) using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium
bromide method. The Vero cells were treated with different concentrations
(1–500 μg/mL) of the tested compounds for 24 h, and the
relative cell viability was calculated at 570 nm using an ELISA microplate
reader. The results showed that the compounds exhibited a percentage
decrease in cell viability in a concentration-dependent manner (Figures a–7a). The IC50 values of MK5, MK6, MK12, and MK14 were calculated as 218.4, 149.1, 99.96, and
162.3 μg/mL, respectively (Figures b–7b), from
a dose–response curve plotted using the GraphPad Prism 6.0
software, whereas the EC50 values of the respective compounds
were at the nanomolar level. The biological safety of the compounds
was found to be in the order of MK5 > MK14 > MK6 > MK12. Additionally, the
effect
of the compounds on cellular morphology was analyzed using phase-contrast
microscopy. The cell membrane integrity and reduction in cell numbers
associated with cellular viability were demonstrated morphologically
in the Vero cells. The results indicated that MK5, MK6, and MK14 exhibited no signs of toxicity
at a concentration of 100 μg/mL (Figures c, 5c, and 7c), whereas MK12 (Figure c) was the most toxic of the
four tested compounds. The cells exposed to higher concentrations
(300 μg/mL) showed the marked morphological alterations typically
associated with cytotoxicity, such as a marked reduction in cellular
density, cellular shrinkage, and blebbing. This study revealed that MK5, MK6, MK12, and MK14 were biologically safe compounds with IC50 values ≥100
μg/mL, and the effective concentration of the compounds was
at a level of <10 nM.
Figure 4
Effect of MK5 on the cell viability
of Vero cells:
(a) cell viability >75% up to 300 μg/mL; (b) representation
of a dose–response curve with an IC50 value of 218.4
μg/mL (696.1 μM); (c) morphological studies of Vero cells
with different concentrations under a phase-contrast microscope, exposed
for 24 h. The control value was 100%, and the data were presented
as the means ± SEs from three independent experiments.
Figure 7
Effect of MK14 on the cell viability of Vero
cells:
(a) cell viability >80% at 60 μg/mL; (b) representation of
a
dose–response curve with an IC50 value of 162.3
μg/mL (614.0 μM); (c) morphological studies of Vero cells
with different concentrations under a phase-contrast microscope, exposed
for 24 h. The control value was 100%, and the data were presented
as the means ± SEs from three independent experiments.
Figure 5
Effect of MK6 on the cell viability of Vero
cells:
(a) cell viability > 70% at 80 μg/mL; (b) representation
of
a dose–response curve with an IC50 value of 149.1
μg/mL (476.1 μM); (c) morphological studies of Vero cells
with different concentrations under a phase-contrast microscope, exposed
for 24 h. The control value was taken as 100%, and the data were presented
as the means ± SEs from three independent experiments.
Figure 6
Effect of MK12 on the cell viability of Vero
cells:
(a) cell viability >80% at 40 μg/mL; (b) representation of
a
dose–response curve with an IC50 value of 99.96
μg/mL (354.1 μM); (c) morphological studies of Vero cells
with different concentrations under a phase-contrast microscope, exposed
for 24 h. The control value was 100%, and the data were presented
as the means ± SEs from three independent experiments.
Effect of MK5 on the cell viability
of Vero cells:
(a) cell viability >75% up to 300 μg/mL; (b) representation
of a dose–response curve with an IC50 value of 218.4
μg/mL (696.1 μM); (c) morphological studies of Vero cells
with different concentrations under a phase-contrast microscope, exposed
for 24 h. The control value was 100%, and the data were presented
as the means ± SEs from three independent experiments.Effect of MK6 on the cell viability of Vero
cells:
(a) cell viability > 70% at 80 μg/mL; (b) representation
of
a dose–response curve with an IC50 value of 149.1
μg/mL (476.1 μM); (c) morphological studies of Vero cells
with different concentrations under a phase-contrast microscope, exposed
for 24 h. The control value was taken as 100%, and the data were presented
as the means ± SEs from three independent experiments.Effect of MK12 on the cell viability of Vero
cells:
(a) cell viability >80% at 40 μg/mL; (b) representation of
a
dose–response curve with an IC50 value of 99.96
μg/mL (354.1 μM); (c) morphological studies of Vero cells
with different concentrations under a phase-contrast microscope, exposed
for 24 h. The control value was 100%, and the data were presented
as the means ± SEs from three independent experiments.Effect of MK14 on the cell viability of Vero
cells:
(a) cell viability >80% at 60 μg/mL; (b) representation of
a
dose–response curve with an IC50 value of 162.3
μg/mL (614.0 μM); (c) morphological studies of Vero cells
with different concentrations under a phase-contrast microscope, exposed
for 24 h. The control value was 100%, and the data were presented
as the means ± SEs from three independent experiments.
ROS Assay
The effect of MK5, MK6, MK12, and MK14 on intracellular
ROS scavenging
was tested using Vero cells. The cells were exposed to H2O2 for excess ROS generation for 10 min, and the respective
drugs were treated as described previously.[24] The images were taken under a fluorescence microscope. Following
drug treatment, the intracellular ROS generation in H2O2-treated Vero cells decreased, and calculation of the intensity
of fluorescence after 24 h of drug treatment indicated balanced pro-oxidant
and antioxidant levels in the cell system (Figure a–d). The effective compounds MK5, MK6, MK12, and MK14 exhibited a significant ROS scavenging effect compared with H2O2-treated control cells. Excess generation of
ROS has been reported to induce oxidative stress in the brain, thereby
leading to neuronal damage in neurodegenerative diseases. Severe nerve
damage can be controlled by balancing ROS generation and scavenging
by antioxidants.[25] This study concluded
that the effective compounds efficiently controlled the ROS produced
by H2O2 treatment, and they were biologically
safe. The compounds may be considered future therapeutics if their
efficacy is further confirmed by preclinical trials.
Figure 8
Effects of MK5, MK6, MK12, and MK14 on ROS
levels induced by H2O2. Vero cells were pretreated
with respective drugs for 24
h. These cells were exposed to 100 μg/mL H2O2 for 10 min, and ROS production was evaluated.
Effects of MK5, MK6, MK12, and MK14 on ROS
levels induced by H2O2. Vero cells were pretreated
with respective drugs for 24
h. These cells were exposed to 100 μg/mL H2O2 for 10 min, and ROS production was evaluated.
Computational Studies
Computational Analysis Based on the MM/PBSA
Method
The binding modes of the lead molecules MK5, MK6, MK12, and MK14 were
established using
the MM/PBSA method. Binding
free energy provides insights into the binding affinity of a compound
with its target, and it is an important parameter for hit-to-lead
and lead optimization in drug discovery.[26] Binding affinity estimations of the compounds would therefore provide
insights into the molecular basis of their activity against MAO-B.
In this report, we employed the MM/PBSA approach, which has been widely
used to estimate binding free energies due to its reliability and
cheaper cost than experimental methods.[27] To calculate binding free energy, snapshots taken between 100 and
200 ns were used to ensure that all the simulated models had equilibrated.
The MM/PBSA method considers several energy contributions, including
van der Waals (vdW) and electrostatic interactions, polar solvation
energy, and nonpolar solvent-accessible surface area energy. Entropy
energy contributions were not considered because this research investigated
each compound’s binding only to MAO-B but with different binding
modes; hence, entropic energy contributions would have produced minimal
differences between binding modes. The degree of binding affinity
reflected the strength of interactions between the compounds and MAO-B
and therefore indicated their inhibitory potential. The MM/PBSA calculation
results, presented in Table , showed estimated total binding free energies of −40.86,
−43.47, −41.69, and −38.92 kcal/mol for MK5, MK6, MK12, and MK14, respectively. Of the compounds, MK14 exhibited the
highest binding free energy, while MK12 (−38.92
kcal/mol) had the lowest binding free energy. Overall, all the compounds
had favorable binding affinities with MAO-B, characterized by the
high energy contributions of vdW and electrostatic interactions with
binding site residues, which supported our experimental findings,
with IC50 values for MK5, MK6, MK12, and MK14 of 0.0040, 0.0028, 0.0032,
and 0.0049 μM, respectively. These values corresponded with
high binding affinities and reflected the favorable functional strength
of each compound as a potential drug.
Table 2
MM/PBSA
Binding Free Energy Profiles
of MK5, MK6, MK12, and MK14
compounds
ΔEvdw (kcal/mol)
ΔEele (kcal/mol)
ΔGgas (kcal/mol)
ΔGsol (kcal/mol)
ΔGbind (kcal/mol)
MK5
–45.41 ± 0.04
–7.97 ± 0.04
–53.38 ± 0.05
12.52 ± 0.03
–40.86 ± 0.04
MK6
–46.77 ± 0.05
–10.70 ± 0.06
–57.47 ± 0.07
14.00 ± 0.04
–43.47 ± 0.07
MK12
–50.44 ± 0.04
–4.59 ± 0.04
–55.03 ± 0.05
13.35 ± 0.03
–41.69 ± 0.05
MK14
–43.16 ± 0.04
–9.76 ± 0.05
–52.92 ± 0.05
13.99 ± 0.03
–38.92 ± 0.04
Binding Site Energetics
That Characterize MK5, MK6, MK12, and
MK14 Binding
Having estimated the binding free energies of MK5, MK6, MK12, and MK14 toward MAO-B, we proceeded to explore the energetics of each binding
site by quantifying the energy contribution of each binding site residue
using the per-residue energy decomposition component of the MM/PBSA
approach in AMBER 18.[28] Decomposition of
the energetics of the binding site residues allowed us to identify
residues that were crucial to the binding of each residue while providing
a molecular perspective on the possible binding mechanism of each
compound. Residues that contributed total energies ≤−1
kcal/mol were considered crucial to the binding of the corresponding
compound and could inform future drug design processes for novel MAO-B
inhibitors. As shown in Figure A, the major interactions constituting the binding of MK5 included LEU171 (−1.44 kcal/mol), ILE198 (−1.12
kcal/mol), ILE199 (−1.47 kcal/mol), TYR398 (−1.97 kcal/mol),
and TYR435 (−1.78 kcal/mol). These residues were shown to correspondingly
engage in high-affinity interactions with MK5, as shown
in Figure A1. The major residues involved in the binding mechanism of MK6 included LEU171 (−1.51 kcal/mol), TYR188 (−2.07 kcal/mol),
ILE199 (−1.85 kcal/mol), TYR398 (−2.15 kcal/mol), and
TYR435 (−2.51 kcal/mol), as evidenced by high-affinity interactions,
such as conventional hydrogen bonds, π–π stacked–T-shaped
interactions, and π-alkyl–alkyl interactions, as shown
in Figure B,B1. These high-affinity interactions culminated in the highest
binding free energy of MK6 compared to the other compounds.
Figure 9
Per-residue
energy contribution plot and corresponding ligand interaction
profile of MAO-B binding pocket residues in a complex of MK5 (A) and MK6 (B). The ligand interaction profile highlighted
the accompanying interactions and revealed the contributed binding
free energies.
Per-residue
energy contribution plot and corresponding ligand interaction
profile of MAO-B binding pocket residues in a complex of MK5 (A) and MK6 (B). The ligand interaction profile highlighted
the accompanying interactions and revealed the contributed binding
free energies.Based on the per-residue energy
decomposition of the MK12–MAO-B complex, the major
binding site residues that were
identified as crucial in the binding process included TYR60 (−1.35
kcal/mol), PHE343 (−1.62 kcal/mol), TYR398 (−3.05 kcal/mol),
and TYR435 (−1.34 kcal/mol), as shown in Figure A. The MK12 binding
was also characterized by notable interactions, such as π–π
stacked–T-shaped interactions, π–sulfur interaction,
and halogen interaction, as shown in Figure A1. Likewise, the crucial residues
that contributed to the binding of MK14 included CYS172
(−1.59 kcal/mol), TYR188 (−1.53 kcal/mol), ILE198 (−1.33
kcal/mol), GLN206 (−1.21 kcal/mol), TYR398 (−2.24 kcal/mol),
and TYR435 (−1.27 kcal/mol), as shown in Figure B. These residues mediated
high-affinity interactions, such as π–π stacked–T-shaped
interactions, π-alkyl–alkyl interactions, and conventional
hydrogen bonds, as shown in Figure B1. Overall, MK5, MK6, MK12, and MK14 bind favorably to MAO-B,
as shown by high-affinity interactions with specific residues, and
a pool of vdW interactions, which stabilized each compound within
the MAO-B binding pocket to facilitate the inhibitory activity.
Figure 10
Per-residue
energy contribution plot and corresponding ligand interaction
profile of MAO-B binding pocket residues in a complex of MK12 (A) and MK14 (B). The ligand interaction profile highlighted
the accompanying interactions and revealed the contributed binding
free energies.
Per-residue
energy contribution plot and corresponding ligand interaction
profile of MAO-B binding pocket residues in a complex of MK12 (A) and MK14 (B). The ligand interaction profile highlighted
the accompanying interactions and revealed the contributed binding
free energies.
Structural and Conformational
Implications of the Binding of MK5, MK6, MK12, and MK14
The therapeutic binding
of chemical compounds with biological
targets is usually associated with various conformational and structural
changes that interfere with the normal functions of the biological
targets.[29] The MD simulation performed
in this research allowed for a nanosecond assessment of the structural
changes associated with the MAO-B binding of each compound. The crucial
parameters assessed to provide insights into the structural changes
included enzyme structure stability, enzyme structure flexibility,
and enzyme folding/unfolding dynamics.[30,31] These were
computed by estimating the C-α root-mean-square deviation (RMSD),[32] C-α root-mean-square fluctuation (RMSF),[33] and solvent accessible surface area (SASA)[34] of the trajectories generated after the 200
ns MD simulation.The C-α RMSD measures atomistic deviations
and reflects the stability and convergence of simulated models. As
observed in Figure , all the simulated models converged after about 100 ns, leading
to ensuing deviations being attributed to the presence or absence
of a bound ligand. The initial increase in deviations from the start
of the simulation up to about 100 ns resulted from initial atomic
expansions. As shown in the table and Figure A, although the binding of MK5, MK6, MK12, and MK14 generally
increased the RMSD of the C-α atoms of MAO-B, as evidenced by
the relatively higher average RMSD of the inhibitor-bound systems,
the structure of MAO-B remained generally stable over the simulation
period, with an average RMSD below 2 Å. As shown in Table , MK5, MK6, MK12, and MK14 exhibited
average RMSDs of 1.71, 1.89, 1.56, and 1.84 Å, respectively,
and the unbound MAO-B also had an average RMSD of 1.56 Å. This
suggested that the binding of the compounds with MAO-B was characterized
by an increase in the stability of MAO-B—a feature that could
favor binding interaction dynamics.[35] The
RMSF (a parameter that predicted the residue flexibility of individual
residues of MAO-B in the presence or absence of compounds) was also
calculated. As observed from the RMSF plots in Figure and Table , individual residues in the unbound MAO-B exhibited
relatively lower average RMSFs, suggesting lower residue flexibility.
Comparatively, individual residues in the inhibitor-bound MAO-B, except
for the MK14-bound system, exhibited relatively higher
average RMSFs than the unbound MAO-B, suggesting that the binding
of MK5, MK6, and MK12 induced
residue flexibility, which could interfere with the function of MAO-B,
leading to the observed inhibitory activity and high binding affinity. MK14, on the other hand, decreased the flexibility of individual
residues, as evidenced by an average RMSF of 7.80 Å. The decreased
residue flexibility could have impeded crucial binding interactions
resulting in the low binding affinity calculated for MK14.
Figure 11
Structural and conformational analysis. (A) Comparative RMSD plots
of the inhibitor-bound MAO-B and the unbound MAO-B; (B) comparative
RMSF plots showing per-residue fluctuations across the 200 ns simulation
period for the inhibitor-bound MAO-B and the unbound MAO-B; (C) comparative
SASA plots for the inhibitor-bound MAO-B and the unbound MAO-B across
the 200 ns MD simulation period.
Table 3
Average RMSD, RMSF, and SASA Estimations
for Simulated Models over 200 ns
parameter
APO
MK5
MK6
MK12
MK14
RMSD (Å)
1.56
1.71
1.89
1.56
1.84
RMSF (Å)
9.69
10.71
14.06
11.62
7.80
SASA (Å2)
21,423.24
20,308.47
20,970.32
20,973.83
20,457.29
Structural and conformational analysis. (A) Comparative RMSD plots
of the inhibitor-bound MAO-B and the unbound MAO-B; (B) comparative
RMSF plots showing per-residue fluctuations across the 200 ns simulation
period for the inhibitor-bound MAO-B and the unbound MAO-B; (C) comparative
SASA plots for the inhibitor-bound MAO-B and the unbound MAO-B across
the 200 ns MD simulation period.Furthermore, using SASA calculations, we investigated the impact
of inhibitor binding on the folding and unfolding of MAO-B, considering
the importance of this phenomenon for enzyme functioning.[36] Unfolding/folding of the enzyme structure could
impede or expose individual residues to solvent surfaces, thereby
interfering with binding interactions and enzyme functioning. As observed
in Figure and Table , the unbound MAO-B
had a relatively higher average SASA of 21423.24 Å2 compared with all the bound conformations of MAO-B. This suggested
that, following binding of MK5, MK6, MK12, and MK14, the individual residues of MAO-B
underwent structural rearrangement consistent with structural folding,
the burial of hydrophobic residues, and a reduction in exposure to
solvent surfaces, which tended to influence the functioning of MAO-B.
Of the compounds, the MK5-bound complex had the highest
folding, with an average SASA of 20308.47 Å2, whereas
the least folding occurred in the MK12 complex, with
an average SASA of 20973.83 Å2. Overall, the similarity
in the binding dynamics of the compounds suggested similarity in the
structural mechanisms of inhibition characterized by distortion of
the conformational dynamics of MAO-B.
Conclusions
In this study, 15 halogen-bearing multiconjugated dienones were
synthesized and evaluated for their human MAOs, ChEs, and BACE1 inhibition.
Surprisingly, all derivatives showed a potent selective MAO-B inhibitory
activity in the nanomolar range compared to the reference drugs. MK6 had the most potent inhibitory activity against MAO-B,
with an IC50 value of 2.82 nM, followed by MK12 (IC50 = 3.22 nM). Kinetic and reversibility studies showed
that MK6 and MK12 were competitive and reversible
inhibitors of MAO-B. These compounds exhibited no distinct signs of
toxicity on normal Vero cells in in vitro toxicity studies. Additionally,
pro-oxidant and antioxidant levels were retained by MK6 and MK12. The MD studies provided novel insights into
the binding modes of the inhibitor-binding cavity of MAO-B. Therefore,
this study suggests that MK6 and MK12 have
therapeutic potential for the treatment of various neurodegenerative
disorders, such as AD and PD.
Materials and Methods
Synthesis
Acetophenone/halogenated
acetophenones (0.01
M) were added to a stirred solution of the respective cinnamaldehyde
derivatives (0.01 M) in 20 mL of ethanol using a micropipette. Pyrrolidine
(0.01 M) was added to the mixture immediately, which soon changed
to a brown or orange color. Overnight stirring resulted in multiconjugated
ketones after the addition of ice cubes, which were filtered under
suction, washed thoroughly with water, and then dried in a desiccator
overnight.[37]
MAO inhibitory activities
were assayed
by recombinant MAO-A and MAO-B using kynuramine (0.06 mM) and benzylamine
(0.3 mM) as substrates.[38] Toloxatone and
clorgyline were used as reference compounds for MAO-A, and lazabemide
and pargyline were used for MAO-B. The Km of benzylamine
for MAO-B was 0.17–0.18 mM.[39] For
multitarget analysis, AChE, BChE, and BACE1 inhibitory activities
were tested as described previously.[40]
Enzyme Inhibition and Kinetic Studies
The inhibitory
activities of MK1–MK15 against MAOs
were first screened at 10 μM. For the compounds showing <50%
residual activities, we determined the IC50 values of the
compounds. The SI values of MAO-B were expressed by calculating IC50 (MAO-A)/IC50 (MAO-B). Enzyme kinetics were determined
for compounds MK6 and MK12 with MAO-B at
five different substrate concentrations. The inhibition patterns were
analyzed using Lineweaver–Burk plots and their secondary plots
for three inhibitor concentrations.[41−43]
Inhibition Reversibility
of MK6 and MK12
The dialysis method
was used for the reversibility test
of MAO-B inhibition by MK6 or MK12 after
preincubation with the enzyme for 30 min at ∼2 × IC50 (i.e., 6.0 nM), as previously described.[44,45] For reference compounds, MAO-B was preincubated with lazabemide
(a reference reversible MAO-B inhibitor) or pargyline (a reference
irreversible MAO-B inhibitor) at 0.22 and 0.28 μM, respectively.
Reversibility patterns were determined by comparing the activities
of dialyzed (AD) and undialyzed (AU) samples.
Cytotoxicity and ROS Assays
The cytotoxicities and
ROS quenching abilities of the lead compounds were evaluated as previously
described.[46,47]
Computational Methodology
Detailed procedures for enzyme
refining, ligand preparation, molecular docking, and dynamic simulations
are described in the Supporting Information.