Bilal Ahmad Khan1, Muhammad Ashfaq2, Shabbir Muhammad3, Khurram Shahzad Munawar4,5, Muhammad Nawaz Tahir2, Abdullah G Al-Sehemi3, Saleh S Alarfaji3. 1. Department of Chemistry, University of Azad Jammu and Kashmir, Muzaffarabad 13100, Azad Jammu and Kashmir, Pakistan. 2. Department of Physics, University of Sargodha, Sargodha, Punjab 40100, Pakistan. 3. Department of Chemistry, College of Science, King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia. 4. Institute of Chemistry, University of Sargodha, Sargodha 40100, Pakistan. 5. Department of Chemistry, University of Mianwali, Mianwali 42200, Pakistan.
Abstract
Ethyl 4-(4-fluorophenylamino)-2,6-bis(4-(trifluoromethyl)phenyl)-1-(4-fluoro-phenyl)-1,2,5,6-tetrahydropyridine-3-carboxylate (FTEAA) has been synthesized efficiently in an iodine-catalyzed five-component reaction of 4-fluoroaniline, 4-trifluoromethyl benzaldehyde, and ethyl acetoacetate in methanol at 55 °C for 12 h. Various spectro-analytical techniques such as 1H and 13C NMR and Fourier-transform infrared spectroscopy have validated the structure of FTEAA. Further confirmation of the structure of FTEAA has been established on the basis of single-crystal X-ray diffraction analysis. The supramolecular assembly of FTEAA in terms of strong and comparatively weak noncovalent interactions is fully investigated by Hirshfeld surface analysis, the interaction energy between pairs of molecules, and energy frameworks. The void analysis is conducted to explore the strength and stability of the crystal structure. Furthermore, molecular docking analysis was computationally performed to see the potential intermolecular interactions between the selected proteins and FTEAA. The binding interaction energies are found to be -8.8 and -9.6 kcal/mol for the proteins MAO-B (PDB ID: 2V5Z) and MAO-A (PDB ID: 2Z5X), respectively. These reasonably good binding energies (more negative values) indicate the efficient associations between the FTEAA and target proteins. The proteins and FTEAA were also analyzed for intermolecular interactions. FTEAA and proteins interact in a variety of ways, like conventional hydrogen bonds, carbon-hydrogen bonds, alkyl, π-alkyl, and halide interactions.
Ethyl 4-(4-fluorophenylamino)-2,6-bis(4-(trifluoromethyl)phenyl)-1-(4-fluoro-phenyl)-1,2,5,6-tetrahydropyridine-3-carboxylate (FTEAA) has been synthesized efficiently in an iodine-catalyzed five-component reaction of 4-fluoroaniline, 4-trifluoromethyl benzaldehyde, and ethyl acetoacetate in methanol at 55 °C for 12 h. Various spectro-analytical techniques such as 1H and 13C NMR and Fourier-transform infrared spectroscopy have validated the structure of FTEAA. Further confirmation of the structure of FTEAA has been established on the basis of single-crystal X-ray diffraction analysis. The supramolecular assembly of FTEAA in terms of strong and comparatively weak noncovalent interactions is fully investigated by Hirshfeld surface analysis, the interaction energy between pairs of molecules, and energy frameworks. The void analysis is conducted to explore the strength and stability of the crystal structure. Furthermore, molecular docking analysis was computationally performed to see the potential intermolecular interactions between the selected proteins and FTEAA. The binding interaction energies are found to be -8.8 and -9.6 kcal/mol for the proteins MAO-B (PDB ID: 2V5Z) and MAO-A (PDB ID: 2Z5X), respectively. These reasonably good binding energies (more negative values) indicate the efficient associations between the FTEAA and target proteins. The proteins and FTEAA were also analyzed for intermolecular interactions. FTEAA and proteins interact in a variety of ways, like conventional hydrogen bonds, carbon-hydrogen bonds, alkyl, π-alkyl, and halide interactions.
The piperidine scaffold,
the most important structural feature
of many pharmaceuticals and found in the top-selling small molecule
pharmaceutical, is ubiquitous in many natural products and medicinal
products. Moreover, the piperidine moiety behaves as the main chromophore
in several commercial drugs, including Aricept, alogliptin, desloratadine,
raloxifene, and Ritalin (Figure ).[1,2]
Figure 1
Commercial drugs bearing the piperidine
moiety.
Commercial drugs bearing the piperidine
moiety.In addition, piperidines have widespread pharmacological
applications
as antibacterial, antifungal, anticancer, anti-tuberculosis, elastase
inhibitors, antiobesity, antianalgesic, acetylcholine serine protease
inhibitor, antihypertensive, antidepressant, neuroprotective, and
antipsychotics.[3,4] Due to the significant relevance
of piperidine derivatives belonging to the pharmacologically active
scaffold in different therapeutic drugs, their synthesis has been
extensively studied in the literature.[5,6]Dysthymia
and nervous ailments such as Parkinson’s and Alzheimer’s
diseases arise because of an over-release of monoamine oxidases, which
diminishes the levels of monoaminergic transmitters and, hence, the
communication in the central nervous system. The communication and
intracellular levels of monoamine oxidases can be maintained by administering
monoamine oxidase inhibitors.[7] Among monoamine
oxidases, monoamine oxidase A has the potential to catalyze the deamination
of adrenaline and noradrenaline as well as serotonin; however, monoamine
oxidase B catalyzes the disintegration of phenyl ethyl amine (a central
nervous stimulant) and phenyl methyl amine.[8] Additionally, monoamine oxidases A and B can break down dopamine,
tryptamine, and tyramine.[9] Monoamine oxidases
are administrated by the use of cordyline,[10] moclobemide (MAO-A),[11] and selegiline
(MOA-B).[12]Owing to their tendency
to break down monoamine neurotransmitters,
monoamine oxidases are effective protein targets for drug development.
The molecules disclosing inhibition potential in opposition to MAO-A
and B could serve as potential drugs to cure depression and Parkinson’s
and Alzheimer’s diseases.[13]Recently, cis- and trans-isomers
of 1-propargyl-4-styrylpiperidine have been reported as potent inhibitors
of MAO-A and MAO-B having IC50 values of 0.7261 + 0.0269
and 0.3422 + 0.0224 μM, respectively[14] (Figure ).
Figure 2
(a) Piperidine-based
geometric isomers selectively inhibit monoamine
oxidase A and B. (b) Tetrahydropyridine-based inhibitors of monoamine
oxidase A and B.
(a) Piperidine-based
geometric isomers selectively inhibit monoamine
oxidase A and B. (b) Tetrahydropyridine-based inhibitors of monoamine
oxidase A and B.Tetrahydropyridines, structures with one double
bond in the piperidine
nucleus, possess good monoamine oxidase substrate properties. 1-Methyl-4-phenyl-1,2,3,6-tetrahydropyridine
(MPTP) and its analogues are inhibitors of monoamine oxidase B.[15] Installation of various substituents on the
tetrahydropyridine has shown an increase in activity toward monoamine
oxidase B inhibition (Figure b).[16] The introduction of fluorine
and trifluoromethyl groups onto the core structure of MPTP provided
better enzyme substrates for monoamine oxidase A and B[17] as the fluorine and trifluoromethyl groups are
known to impart enhanced metabolic stability, lipophilicity, and dipole
moment to drug-like molecules.[18,19]The current study
involved the synthesis and characterization of
a densely substituted tetrahydropyridine-based derivative, ethyl 4-(4-fluorophenylamino)-2,6-bis(4-(trifluoromethyl)phenyl)-1-(4-fluoro-phenyl)-1,2,5,6-tetrahydropyridine-3-carboxylate (FTEAA). Furthermore, the synthesized derivative is evaluated
for potential dual inhibition of MAO-A and B proteins. Several important
structural aspects are determined through the X-ray diffraction (XRD)
technique for the single crystal structure of FTEAA.
As intermolecular interactions play a vital role in crystal formation,
the supramolecular assembly of FTEAA in terms of strong
as well as, comparatively, weak noncovalent interactions is fully
investigated by Hirshfeld surface (HS) analysis. Furthermore, docking
studies to examine the binding interactions of FTEAA within
the active pockets of MAO-A and B proteins supplement the current
experimental findings regarding inhibition activity.
Results and Discussion
Single Crystal Analysis
Single crystal
XRD is utilized for crystal structure determination. The Cambridge
structure database search ratifies the originality of the crystal
structure of FTEAA. The experimental details are listed
in Table , while selected
bond lengths and bond angles are given in Table .
Table 1
Single-crystal Analysis of FTEAA
crystal data
FTEAA
CCDC
2173544
chemical formula
C34H26F8N2O2
Mr
646.57
crystal system, space group
triclinic, P1̅
temperature
(K)
150 (2)
a, b, c (Å)
9.9381 (6), 11.7975
(7),
14.0204 (9)
α, β, γ (deg)
106.603 (5), 105.040 (6),
99.032 (5)
V (Å3)
1473.11 (17)
Z
2
density (calculated)g/cm–3
1.458
F(000)
664
radiation type
Cu Kα
wavelength (λ)
1.54184
μ (mm–1)
1.088
crystal size (mm)
0.28 × 0.24 × 0.17
data collection
diffractometer
Bruker APEXII CCD
absorption correction
multiscan (SADABS; Bruker,
2007)
no.
of measured, independent,
and observed [I > 2σ(I)] reflections
9288, 4663, 3772
Rint
0.034
θ range for data collection (deg)
3.474–62.769
index
ranges
–11 ≤ h ≤ 11, –13 ≤ k ≤ 12, –15 ≤ l ≤ 16
data
refinement
R[F2 > 2σ(F2)], wR(F2), S
0.042, 0.116, 1.04
no. of reflections
4663
no. of parameters
426
H-atom treatment
H-atom parameters constrained
Δρmax, Δρmin (e Å–3)
0.32, −0.27
Table 2
Selected Bond Lengths (Å) and
Bond Angles (deg) of FTEAA
selected
bond lengths
selected bond angles
N1–C1
1.470 (2)
C1–N1–C29
121.33 (15)
N1–C5
1.453 (2)
C5–N1–C29
119.98 (14)
N1–C29
1.394 (2)
C1–N1–C5
118.59 (14)
N2–C3
1.354 (3)
C3–N2–C16
125.50 (14)
N2–C16
1.431 (3)
C2–C13–O1
124.40 (18)
O1–C13
1.231 (2)
C2–C13–O2
113.77 (15)
O2–C13
1.343 (2)
O1–C13–O2
121.81 (16)
O2–C14
1.451 (2)
C13–O2–C14
116.17 (14)
F4–C19
1.354 (3)
O2–C14–C15
107.23 (16)
F8–C32
1.365 (2)
C18–C19–F4
118.9 (2)
In FTEAA (Figure a and Table ), the 1,2,3,6-tetrahydropyridine group A (C1–C5/N1)
is puckered with a total puckering amplitude of 0.681 (2) Å,
θ = 84.69 (17)° and φ = 248.83(17) Å. The C-atom
(C4) deviation from the mean plane of group A is greater than that
of other atoms in the mean plane, with a value of 0.4187 (1) Å.
The conformation of the tetrahydropyridine ring is “a flatted
boat” conformation (Figure b), which is substantiated by the puckering of C1 and
C4 atoms to the extent of 0.5010 (3) and 0.6573 (3) Å, respectively,
from the least square plane defined by (C2/C3/C5/N1) atoms. The flattened
boat conformation is also confirmed by the torsional angles associated
with the tetrahydropyridine ring that are C(5)–N(1)–C(1)–C(2)
= 37.1 (2)°, N(1)–C(1)–C(2)–C(3) = −46.8
(2)°, C(1)–C(2)–C(3)–C(4) = 2.7 (2)°,
C(2)–C(3)–C(4)–C(5) = 49.0 (2)°, C(1)–N(1)–C(5)–C(4)
= 12.6 (2)°, and C(3)–C(4)–C(5)–N(1) = −55.94
(19)°. The phenyl ring B (C6–C11), ethyl formate group
C (C13–C15/O1/O2), 4-fluoroaniline group D (C16–C21/N2/F4),
phenyl ring E (C22–C27), and fluorophenyl ring F (C29–C34/F8)
are planar with root mean square deviations of 0.0087, 0.0126, 0.0171,
0.0039, and 0.0076 Å, respectively. The mean plane of group A
makes the dihedral angles of 82.32 (6), 32.79 (9), 72.39 (6), 82.42
(6), and 28.04 (9)° with ring B and groups D, E, and F, respectively.
The dihedral angles infer that group A is nearly perpendicular to
ring B and group E. One of the trifluoro groups (F1A–F3A) is
disordered over two sets of locations with an occupancy ratio of 0.508(6):
0.492(6). The mean plane of the major part (F1A–F3A) is orientated
at a dihedral angle of 11.89 (3)°, relative to the mean plane
of the minor part (F1B–F3B). The anisotropic displacement parameter
(ADP) of the major part of each disordered atom is made equal to the
ADP of the minor part. The DFIX and DANG restraints are used for stabilizing
disordered parts in refinement. The molecular configuration is stabilized
by intramolecular N–H···O bonding to form an
S (6) loop. The molecules are connected and arranged in the shape
of dimers via C–H···F bonding (Figure ). As a result of dimerization,
an R22(8) loop
is formed. The dimers are then interconnected via C–H···O
and N–H···F bonding. An R22(14) loop is generated by C–H···O
bonding. The crystal packing is further stabilized by intra and as
well as intermolecular C–H···π interactions
(Table ). Intermolecular
C–H···π interactions connect the molecular
units in the form of a dimer. One such dimer is shown in Figure .
Figure 3
(a) ORTEP diagram of FTEAA drawn at a probability
level of 50%. Hydrogen atoms are visible in the form of small circles
of arbitrary radii. Only the main portion of the disordered trifluoro
moiety is displayed for simplicity. (b) Conformation of the tetrahydropyridine
ring in FTEAA with all atoms other than the tetrahydropyridine
ring was omitted for clarity.
Figure 4
Packing diagram of FTEAA. Only selected hydrogen
atoms
and the main part of the disordered trifluoro group are shown for
simplicity.
Figure 5
Graphical illustration of C–H···π
interactions
in FTEAA.
(a) ORTEP diagram of FTEAA drawn at a probability
level of 50%. Hydrogen atoms are visible in the form of small circles
of arbitrary radii. Only the main portion of the disordered trifluoro
moiety is displayed for simplicity. (b) Conformation of the tetrahydropyridine
ring in FTEAA with all atoms other than the tetrahydropyridine
ring was omitted for clarity.Packing diagram of FTEAA. Only selected hydrogen
atoms
and the main part of the disordered trifluoro group are shown for
simplicity.Graphical illustration of C–H···π
interactions
in FTEAA.The Cambridge structural database (CSD 9) provides
nine crystal
structures that are closely related to the crystal structure of FTEAA. The reference codes of these nine structures are COJCIP,[20] GATCUB,[21] KAGSAQ,[22] KAGSEU,[17] LETBET,[23] MEJTAZ,[24] MEWTOZ,[25] MONVAO,[26] and WIHCOH.[27] The S (6) intramolecular H-bonding loop was
formed through intramolecular N–H···O bonding
between the amino group and carbonyl O-atom in all the nine crystal
structures of the CSD. The same H-bonded loop is formed in FTEAA by the intramolecular H-bonding between the amino group and carbonyl
O-atom. The tetrahydropyridine ring adopted a boat conformation, with
some deviation in all the nine crystal structures of the CSD except
in the crystal structures with reference code COJCIP and MEJTAZ. The
authors of the crystal structure with reference code MEJTAZ did not
describe the conformation of the tetrahydropyridine ring as well as
intra- and intermolecular H-bonding. The tetrahydropyridine ring adopted
a twisted chair conformation in the crystal structure with reference
code COJCIP. Just like in FTEAA, the intermolecular C–H···O
and C–H···F H-bonding interactions are present
in the crystal structures with reference codes COJCIP, LETBET, MEWTOZ,
and WIHCOH, whereas C–H···F bonding is absent
in KAGSEU and C–H···O bonding is absent in KAGSAQ. Figure S1 shows the simulated powder pattern
of FTEAA. The sharp peaks indicated high crystallinity,
and the most intense peaks in the simulated powder XRD pattern indicated
the preferred growth direction.
HS Exploration
HS analysis is achieved
on Crystal Explorer version 21.5.[28] It
is used for the exploration of the noncovalent interactions in the
crystal packing of single crystals.[29] The
properties of the single crystal depend on the packing of molecules
and noncovalent interactions that interlinked the molecules with each
other. The HS concept was developed when scientists were trying to
divide the electronic density of crystals into molecular fragments
for integration purposes.[30,31] HS can be plotted by
using various properties like dnorm, shape index, curvedness, and
so forth. The hydrogen bonding information can be extracted by drawing
a graph of the HS over dnorm, where dnorm is the normalized distance.[32,33] This surface contained three colors namely red, blue, and white.
The red and blue regions show the short interatomic contacts and long
interatomic contacts, respectively.[34,35] The white
regions stand for the contacts having an interatomic distance equal
to the sum of van der Waals radii. The short contacts or the contacts
that are suitable for the hydrogen bonding interactions are depicted
by the red spots on the HS (Figure a). The green dotted lines in Figure a represent H-bonding interactions. π···π
stacking interactions can be visualized by HS that is plotted over
the shape index.[36,37] The existence of consecutive
red and blue triangle areas on the HS around aromatic rings indicates
their involvement in π···π stacking interactions.[38] The presence of such regions on the HS indicates
the existence of the π···π stacking interactions
in FTEAA (Figure b). However, these π···π stacking
interactions are very weak as the intercentroid distance for this
interaction ranges from 4.5501 (13) to 5.9959 (12) Å.
Figure 6
HS plotted
over (a) dnorm in the range of −0.2348 to 1.5584
au; (b) shape index in the range of −1 to 1 au.
HS plotted
over (a) dnorm in the range of −0.2348 to 1.5584
au; (b) shape index in the range of −1 to 1 au.2D fingerprint plots express the involvement of
the interatomic
contacts in the packing of the crystal.[39,40] di and de represents the distances from the HS to the nearest
atom both inside and outside of it, respectively. Figure a is the 2D plot for all kinds
of probable interactions in FTEAA. The fluorine atoms
play an important role in the packing of the crystal, as the most
significant contact for FTEAA is H–F or F–H
(34.1%), as shown in Figure b. Figure f represents HS associated with H–F or F–H contact,
which clearly shows the involvement of the fluorine atoms in the crystal
packing interactions. The interacting region is shown by the red and
blue regions of the surface. The other important contact is H–H,
whose contribution to crystal packing is very close to that of H–F
contacts (Figure c).
The blue regions around H-atoms in Figure g show short H–H contacts. As expected,
while looking at the hydrogen bond geometry table (Table ), C–H contacts must
have a noteworthy involvement (15%) in the packing of the crystal.
The most important interatomic contact that involves O-atoms is O–H,
with a percentage contribution of 4.1% (Figure e). The contacts with the relatively smaller
contribution to the packing of the crystal are manifested in Figure S2.
Figure 7
(a) 2D plots for overall interaction.
(b–i) 2D plots of
important contacts and their HS.
Table 3
Hydrogen-bond Geometry (Å, deg)a
D–H···A
D–H
H···A
D···A
∠(D–H···A)°
N2–H2···O1
0.88
2.11
2.683 (2)
122
N2–H2···F5i
0.88
2.34
3.067 (2)
140
C11–H11···O1ii
0.95
2.44
3.345 (2)
159
C24–H24···F7iii
0.95
2.51
3.356 (2)
148
C15–H15C···F7iv
0.98
2.62
3.570 (3)
164
C31–H31···F8v
0.95
2.48
3.412 (2)
168
C–H···π
C–H
H···π
C··· π
∠(C–H··· π)°
C21–H21···Cg1
0.95
2.95
3.508 (2)
119
C15–H15B···Cg2vi
0.98
2.73
3.624 (2)
152
Symmetry codes: (i) −x + 2, −y + 2, −z + 1; (ii) −x + 2, −y + 1, −z + 1; (iii) −x + 1, −y + 2, −z +
1; (iv) x, y–1, z; (v) −x, −y + 1,
−z; (vi) 1 – x, 1
– y, 1 – z. Cg1 and
Cg2 are the centroid of phenyl rings (C22–C27) and (C29–C34),
respectively.
(a) 2D plots for overall interaction.
(b–i) 2D plots of
important contacts and their HS.Symmetry codes: (i) −x + 2, −y + 2, −z + 1; (ii) −x + 2, −y + 1, −z + 1; (iii) −x + 1, −y + 2, −z +
1; (iv) x, y–1, z; (v) −x, −y + 1,
−z; (vi) 1 – x, 1
– y, 1 – z. Cg1 and
Cg2 are the centroid of phenyl rings (C22–C27) and (C29–C34),
respectively.The enrichment ratio for the pair of chemical moieties
in the crystal
provides the tendency of that pair of chemical species to make the
crystal packing interactions. The enrichment ratio is obtained by
dividing the proportion of the actual contact in the crystal packing
with the proportion of random contact calculated theoretically. Table S1 shows the enrichment ratio of the pair
of chemical moieties in FTEAA. The O–H contact
has the maximum tendency to generate crystal packing interactions
in FTEAA, with an enrichment ratio of 1.43. The other
contacts of higher propensity are C–H, F–H, and F–C
with enrichment ratios of 1.19, 1.05, and 1.04, respectively.To explore the interactions between pairs of molecules, interaction
energy is calculated by creating a bunch of molecules of size 3.8
Å around the reference molecule (Figure a). The B3LYP/6-31G(d,p) level electron density
model is used for the calculations.[41] Fifteen
symmetry-related molecules are present around the reference molecule,
and the coordinates of atoms that are involved in the interaction
energy calculation are given in Table S2. The total energy between the molecular pair is the sum of four
energies known as electrostatic, polarization, dispersion, and repulsion.
The negative energy value designates that the energy is attractive,
whereas the positive energy value indicates repulsive energy. The
total energy is greater for the molecule pairs for which the distance
between the molecular center is small. Figure b shows that for most of the molecular pairs,
the dispersion energy term dominates over other energy terms. Electrostatic
energy is mostly attractive but, in some cases, it may be repulsive.
For example, in the present case, it is repulsive for the molecular
pairs, with the distance between the molecular centers being 10.04,
13.21, and 16.36 Å.
Figure 8
(a) Chemical diagram of the color-coded interaction
mapping inside
3.8 Å of the reference molecule; (b) interaction energy results.
(a) Chemical diagram of the color-coded interaction
mapping inside
3.8 Å of the reference molecule; (b) interaction energy results.Energy frameworks for a single crystal are important
to realize
the topology of the single crystal.[42] Interaction
energies are used as input to form energy frameworks. The size or
width of the cylinders is directly proportional to the strength of
the interaction energy type. The width of the cylinders is smaller
in the case of electrostatic energy (Figure a) as compared to the width of the cylinders
in the case of dispersive energy (Figure b), which indicates that the contribution
of dispersive energy is greater in defining the total energy as compared
to the contribution of electrostatic energy. Figure c is the energy frame for the total energy,
which is obtained by adding up the contributions of electrostatic
energy and dispersive energy.
Figure 9
Energy framework for (a) coulomb energy, (b)
dispersion energy,
and (c) total energy.
Energy framework for (a) coulomb energy, (b)
dispersion energy,
and (c) total energy.A void analysis is important to understand the
response of the
single crystals to the applied stress. It is related to the mechanical
properties of single crystals. The voids in the crystal packing are
calculated by using the Hartree–Fock theory.[43,44] The electron density of all the atoms present in the molecule is
added up while assuming that the atoms are spherically symmetrical
(Figure ). The volume
of the void is determined to be 216.02 Å3. As voids
occupy a very small amount of space (14.6%) in the crystal packing,
this suggests that the molecules are tightly packed and that the crystal
packing has no big cavities.
Figure 10
Graphical views of voids in the crystal packing
of FTEAA along (a) c-axis and (b) b-axis.
Graphical views of voids in the crystal packing
of FTEAA along (a) c-axis and (b) b-axis.
1H NMR and 13C NMR Spectroscopy
In the 1H NMR spectrum, three protons of the methyl
group were detected as a triplet at δ = 1.47 ppm. Two methylene
protons (OCH) were observed
as a multiplet at δ = 4.51–4.32 ppm. Two methylene protons
present in the piperidine ring were observed as two doublets of doublets
representing one proton each at δ = 2.65 and 2.81 ppm. A multiplet
for N–CH–CH2 of the ring was
found at δ = 5.09–5.19 ppm. Another N–CH– was found as a singlet at δ = 6.39 ppm. Protons present
in the trifluoro methylated aromatic ring were observed downfield
as compared to the fluoro-containing aromatic ring. All the sixteen
aromatic protons were observed between δ = 7.60–6.72
ppm. The proton present in the −NH was found as
a broad singlet at δ = 10.19 ppm.Similarly, the 13C NMR spectrum of FTEAA has shown carbon peaks
in total agreement with the carbon atoms present in FTEAA. The carbon atom of the methyl group was found at δ = 14.43
ppm, methylene attached to methyl was found at δ = 33.26 ppm,
the carbonyl carbon was observed at δ = 167.54 ppm, and the
alkenic carbon attached to carbonyl was found at δ = 97.11 ppm,
whereas the second alkenic carbon attached to nitrogen atom was observed
at δ = 137.08 ppm. All the aromatic carbon atoms were observed
in the range of δ = 155.39–112.51 ppm. Trifluoro methylated
carbon atom was found as a quartet at 125.01 ppm.
Monoamine Oxidase Inhibition
FTEAA was evaluated for its inhibitory potential against MAO-A
and MAO-B enzymes. For this purpose, clorgyline and deprenyl were
utilized as the standard drugs for MAO-A and MAO-B, respectively. FTEAA showed excellent inhibition of both MAOs in the lower
micromolar range. The inhibitory concentration (IC50) values
of the tested compound and standard inhibitors are shown in Table .
Table 4
Monoamine Oxidase A and B Inhibition
IC50 ± SEM (μM)/%age inhibition
code
MAO-A
MAO-B
FTEAA
0.52 ± 0.03
1.02 ± 0.11
clorgyline
0.0045 ± 0.0003
61.35 ± 1.13
deprenyl
67.25 ± 1.02
0.0196 ± 0.001
Computational Procedure
Binding Energy and Intermolecular Analysis
The FTEAA was analyzed after the docking that was
carried out with two of the selected proteins, 2V5Z and 2Z5X. The docking studies
gave rise to nine different conformational modes in which FTEAA could fit into the binding pocket of proteins. However, for the
analysis of resulting interactions of FTEAA, its binding
mode with the lowest root mean square deviation value was preferred
over others. The binding energies were −8.8 and −9.6
kcal/mol for proteins 2V5Z and 2Z5X, respectively (see Table ). These reasonably good binding energies (more negative values)
indicate the effective docking process and efficient associations
between the FTEAA and proteins. The protein and FTEAA were then analyzed for intermolecular interactions.
There are diverse interactions between FTEAA and proteins,
including conventional hydrogen bonds; carbon–hydrogen bonds;
and alkyl, π-alkyl, and halide interactions. Visualizing these
interactions is essential and tells us the extent of association between
the macromolecules (proteins) and FTEAA. The docking
interaction visualization results of the FTEAA with the 2V5Z protein have been
verified and shown in Figure . It is obvious that favorable hydrogen bonding has been developed
with ARG36 and ARG38 residues at 2.24 and 2.2 Å distances, respectively.
Furthermore, the halide interactions can be seen with GLU390 and ASP37
residues. The alkyl and π-alkyl interactions can be observed
with PRO234, ARG36, and ARG38 amino acid residues.
Table 5
Binding Energies (kcal/mol) between FTEAA and Selected Proteins, the Number of H-bonds, and Hydrophobic
and Electrostatic Interactions for Docked Proteins
y
binding energy
number and types of interactions
proteins
kcal/mol
no. of H-bonds
hydrophobic
interactions
electrostatic
Interactions
2V5Z
–8.8
ARG36 (2.24 Å), ARG38 (2.20 Å)
3
1
2Z5X
–9.6
ARG76 (2.30 Å), GLN474 (2.73 Å)
1
1
Figure 11
Interactions of the FTEAA with the 2V5Z protein; (a) overview
of FTEAA adjustment in the binding pocket of protein
with interacting residues, where green lines show hydrogen bond interactions
between FTEAA and amino acid residues of 2V5Z; (b) 3-D diagram
of the focused view having FTEAA and protein residues
around it; (c) clear view having various residues and FTEAA interactions in a 2-dimensional illustration.
Interactions of the FTEAA with the 2V5Z protein; (a) overview
of FTEAA adjustment in the binding pocket of protein
with interacting residues, where green lines show hydrogen bond interactions
between FTEAA and amino acid residues of 2V5Z; (b) 3-D diagram
of the focused view having FTEAA and protein residues
around it; (c) clear view having various residues and FTEAA interactions in a 2-dimensional illustration.The docking interaction visualization results of the FTEAA with the 2Z5X protein have been demonstrated in Figure . The docking developed desirable
hydrogen
bonding between the FTEAA and ARG76 with a 2.3 Å
interaction distance. Another hydrogen bond emerged with GLN474 at
the 2.73 Å interaction length. Other considerable associations
include halide (fluorine) interactions with several residues, including
GLU474, LYS199, VAL473, and GLU450. The π-alkyl associations
were observed with LYS440 and TRP472. Carbon–hydrogen bonding
was also observed with certain residues.
Figure 12
Interactions of FTEAA with the 2Z5X protein; (a) overview
of FTEAA adjustment in the binding pocket of the protein
with interacting residues, where green lines show hydrogen bond interactions
between FTEAA and amino acid residues of 2Z5X; (b) 3-D diagram
of the focused view having FTEAA and protein residues
around it and interaction distances; (c) 2-D diagram of the clear
view has various residues and FTEAA interactions.
Interactions of FTEAA with the 2Z5X protein; (a) overview
of FTEAA adjustment in the binding pocket of the protein
with interacting residues, where green lines show hydrogen bond interactions
between FTEAA and amino acid residues of 2Z5X; (b) 3-D diagram
of the focused view having FTEAA and protein residues
around it and interaction distances; (c) 2-D diagram of the clear
view has various residues and FTEAA interactions.
ADMET and Drug Likeness
Usually,
the in silico computationally designed drugs undergo analysis for
their pharmacokinetic properties before being experimented on in vivo
or practicing them clinically. These drugs should follow four out
of five of Lipinski’s rules, which include that the molecular
weight of the drug must not exceed 500 g/mol, hydrogen bond donors
must be under five, hydrogen bond acceptors must essentially be less
than ten, its lipophilicity (log P) should not cross
the digit five, and the number of rotatable bonds must be less than
10. Table represents
the potential of FTEAA to be used as an inhibitor drug.
It may be seen that FTEAA truly follows the rule of five,
and its properties come under the prescribed range. This indicates
the physiochemical favorability of this drug to inhibit the proteins
efficiently and control the disease spread.
Table 6
Predicted Pharmacokinetic Analysis
of FTEAA
compound
intestinal
absorption (%)
total clearance (log/mL/min/kg)
max.
tolerable
dose (log mg/kg/day)
RO5 rule
FTEAA
95.664
0.812
0.493
yes
The toxicological traits and pharmacokinetics of the FTEAA are determined by evaluating ADMET analysis. This analysis
is done
immediately after the docking process so that inhibitors with poor
excretion and absorption properties and hepatotoxicity may be rejected
early. In the FTEAA, it may be seen that intestinal absorption
is high, at about 95.66% (see Table ). The total clearance is 0.446 log/mL/min/kg for the
entitled compound. The tolerable dose is prescribed to be 0.513 log
mg/kg/day in the case of both the proteins.
Conclusions
Efficient synthesis of
highly functionalized piperidine FTEAA was accomplished
by following a single-pot five-component procedure.
The structure of FTEAA was validated with NMR (1H and 13C), Fourier-transform infrared (FT-IR) spectroscopy,
and single-crystal XRD analysis. The supramolecular environment is
elaborated by HS analysis, from which it is inferred that the H–F/F–H
contact is a noteworthy donor to the crystal packing, and the O–H
contact has a higher propensity to make the crystal packing interactions
relative to other contacts. The void analysis confirmed the strength
of the crystal packing. The supramolecular environment is further
explored by the interaction energy between the molecular pairs and
by the construction of the energy frameworks. FTEAA is
found to be an effective inhibitor of both monoamine oxidase A and
B with IC50 values of 0.52 ± 0.03 and 1.02 ±
0.11 μM. Additionally, the in silico molecular docking analysis
also indicated that there is a reasonably good association between
the FTEAA ligand and selected MAO proteins. The binding
interaction energies were calculated through molecular docking methods.
The entitled ligand shows binding interaction energies of −8.8
and −9.6 kcal/mol for proteins 2V5Z and 2Z5X, respectively. The negative binding energies
were further supported by the study and visualization of intermolecular
interactions among entitled ligand and proteins. The ADMET and drug-likeness
studies indicated the potential of FTEAA to be used as
a monoamine oxidase inhibitor.
Experimental Section
Material and Techniques
All the chemicals,
reagents, and solvents were acquired from Sigma-Aldrich and utilized
without further purification. The melting point of FTEAA was measured by using the Gallenkamp melting point apparatus (MP-D).
NMR spectra (1H and 13C) were measured on a
Bruker AV-400 spectrometer (400 MHz). The FT-IR spectrum was recorded
in the ATR (attenuated total reflectance) mode by using a Shimadzu
FT-IR spectrophotometer.
Preparation of Ethyl 4-(4-fluorophenylamino)-2,6-bis(4-(trifluoromethyl)phenyl)-1-(4-fluorophenyl)-1,2,5,6-tetrahydropyridine-3-carboxylate
To prepare compound FTEAA, a 5 mL methanolic solution
of 4-fluoroaniline (2 mmol, 186 mg), ethyl acetoacetate (1 mmol, 130
mg), and molecular iodine (10 mol %, 25 mg) was stirred at 55 °C
for 25 min and 4-trifluoromethyl benzaldehyde (2 mmol, 348 mg) was
added. After 12 h, the contents were cooled down to room temperature
and subsequent precipitates were filtered, washed with absolute ethanol,
dried, and recrystallized from 10% ethyl acetate/n-hexane to obtain FTEAA.Colorless solid; yield
73%; mp 218–220 °C; solubility: acetone, chloroform, dichloromethane,
ethyl acetate and ethanol; 1H NMR (400 MHz, CDCl3): δ ppm 1.47 (t, J = 7.15 Hz, 3H), 2.65 (dd, J = 15.31, 2.76 Hz, 1H), 2.81 (dd, J =
15.31, 5.77 Hz, 1H), 4.28–4.41 (m, 1H), 4.41–4.57 (m,
1H), 5.09–5.19 (m, 1H), 6.26–6.32 (m, 2H), 6.33–6.38
(m, 2H) 6.39 (s, 1H), 6.72–6.87 (m, 4 H), 7.20–7.31
(m, 2H), 7.43 (d, J = 8.28 Hz, 2H), 7.52–7.60
(m, 4H) 10.19 (br s, 1H) (Figure S3). 13C NMR (101 MHz, CDCl3): δ ppm 14.43, 33.26,
54.65, 57.35, 59.66, 76.36, 76.68, 97.11, 112.51, 116.67, 122.44,
122.56, 124.93, 124.97, 125.01, 125.28, 125.31, 125.40, 125.82, 126.36,
126.54, 128.28, 128.66, 128.84, 129.44, 137.08, 145.92, 146.27, 146.28,
147.73, 155.39, 167.54 (Figure S4). FT-IR
ν(cm–1): 3242 (N–H), 3080 (ArCH), 2981
(CH2), 1649 (C=O), 1591 (C=C), 1115 (C–O),
764 (C=C).The one-step five-component synthetic route
was implemented to
access the desired piperidine FTEAA as shown in Scheme . The reported procedure
was followed to prepare the desired compound.[45,46] The purity was established based on thin-layer chromatography. The
structure of FTEAA was confirmed on the basis of 1H and 13C NMR, FT-IR, and single-crystal XRD analysis.
For
the single crystal analysis of FTEAA by X-ray, a Bruker
Kappa Apex-II CCD diffractometer is used along with APEX-II software
for the collection of data. The diffractometer contains an X-ray source
that produces Cu Kα radiations, and data is collected at a low
temperature of 150 K for the best possible results. SHELXT-2014[47] and SHELXL 2019/2[48] are used for the solution and refinement of the crystal structure,
respectively. All the nonhydrogen atoms are assigned ADP in refinement,
whereas all the H-atoms are assigned isotropic displacement parameters
and refined with the help of the riding model. Initially, one of the
trifluoro groups has unusual thermal ellipsoids. The problem is identified
and resolved by solving the disorder in the trifluoro group by using
DFIX, DANG, and EADP restraints. ORTEP-3,[49] PLATON,[50] and Mercury[51] were used for the graphical illustration of the single-crystal
XRD data.
Monoamine Oxidase Activity
The assay
was conducted in white 96-well plates as previously described.[52] Before use, the enzymes MAO-A and MAO-B were
irreversibly blocked at 25 °C for 15–20 min with clorgyline
and deprenyl, respectively. The assay volume for MAO-A was 200 μL,
containing 145 μL buffer, 10 μL of test compound (100
μM end concentration), and 20 μL (26 μg end concentration)
of the freshly prepared enzyme. The assay volume for MAO-B was 100
μL, containing 45 μL of buffer, 10 μL of test compound
(100 μM end concentration), and 20 μL (5 μg end
concentration) of the freshly prepared enzyme. The reaction mixture
was incubated at 37 °C for 10 min separately. After 15 min, 10
μL of the substrate (0.3 mM final concentration) was added,
and at the end, 10 μL of Amplex red (fluorogenic substrate)
was added to each well. The change in fluorescence was observed by
using a FLUOstar Omega fluorescence plate reader (BMG Labtech GmbH,
Ortenberg, Germany). The species with more than 50% inhibition of
either MAO-A or MAO-B were further investigated and IC50 values were calculated with the help of the PRISM 5.0 (GraphPad,
San Diego, California, USA) program.
Computational Methodology
Protein Preparation
The 3-dimensional
structures of MAO-B (PDB ID: 2V5Z, resolution 1.60 Å) and MAO-A (PDB ID: 2Z5X) with a resolution
of 2.20 Å were attained from the Protein Bank database.[53] The MGL Tools were used to prepare the proteins
for molecular docking analysis.[54] Co-crystalized
ligands and water molecules were removed from macromolecules and polar
hydrogens, and Kollman Charges were added.[55]
Ligand Preparation
The experimental
structure of the FTEAA was fully optimized by using density
functional theory methods with B3LYP, along with a 6-311G** basis
set.[56] The FTEAA was prepared
in the PDBQT format in the MGL tool, and torsions were set by converting
rotatable bonds into nonrotatable ones so that the FTEAA could not undergo structural variation or bond modification during
the docking process. Some further details for molecular docking methodology
can be seen from our previous work with similar types of calculations.[57,58]
Molecular Docking
The prepared proteins
and FTEAA were subjected to molecular docking using Autodock
Vina 4.2. A configuration file was prepared, and the grid centers
and grid dimensions were set according to the prereported literature.
However, the grid dimensions were also confirmed by locating the residues
interacting with co-crystallized inhibitors.[59] After docking, the 3-D and 2-D interactions of FTEAA with amino acid residues were visualized by the Biovia Discovery
Studio visualizer (2020).
ADMET Analysis
The physicochemical
properties of the FTEAA were estimated through ADMET
analysis, which includes adsorption, desorption, metabolism, excretion,
and toxicity. This analysis was performed using a reliable online
tool called pkCSM.[60] ADMET properties indicate
the drug-like properties of the FTEAA.
Authors: Peter R Spackman; Michael J Turner; Joshua J McKinnon; Stephen K Wolff; Daniel J Grimwood; Dylan Jayatilaka; Mark A Spackman Journal: J Appl Crystallogr Date: 2021-04-27 Impact factor: 3.304
Authors: Akbar Ali; Muhammad Khalid; Muhammad Fayyaz Ur Rehman; Sadia Haq; Arif Ali; Muhammad Nawaz Tahir; Muhammad Ashfaq; Faiz Rasool; Ataualpa Albert Carmo Braga Journal: ACS Omega Date: 2020-06-18