Saif Ullah1,2, Julie Pelletier3, Jean Sévigny3,4, Jamshed Iqbal1,2. 1. Centre for Advanced Drug Research, COMSATS University Islamabad, Abbottabad Campus, Abbottabad22060, Pakistan. 2. Department of Pharmacy, COMSATS University Islamabad, Abbottabad Campus, Abbottabad22060, Pakistan. 3. Centre de Recherche Du CHU de Québec-Université Laval, QuébecG1V 4G2, QC, Canada. 4. Département de Microbiologie-infectiologie et D'immunologie, Faculté de Médecine, Université Laval, QuébecG1V 0A6, QC, Canada.
Abstract
Aberrant level of ectonucleotide pyrophosphatase/phosphodiesterase-1 and -3 is linked with numerous disorders, for instance, diabetes, cancer, osteoarthritis, chondrocalcinosis, and allergic reactions. These disorders may be cured or minimized by blocking the activity of ENPP1 and ENPP3 isozymes. In this study, arylamide sulphonates were synthesized, characterized, and evaluated for their capability to affect the activity of isozymes ENPP1 and ENPP3. Among the selective inhibitors of ENPP1, compounds 4f and 4q exhibited sub-micromolar IC50 values of 0.28 ± 0.08 and 0.37 ± 0.03 μM, respectively, followed by 7a, with IC50 equal to 0.81 ± 0.05 μM, whereas out of the selective inhibitors of isozyme ENPP3, 4t and 7d preferably lessened the activity to half of the maximal inhibitory concentration of 0.15 ± 0.04 and 0.16 ± 0.01 μM alternatively. In addition, many structures including 4c, 4g, 4k, 4l, 4n, 4o, 4r, 4s, 7b, 7c, and 7e inhibited the activity of both isozymes to a significant level. Enzyme kinetic study of compound 4j revealed an uncompetitive mode of inhibition of ENPP1 isozyme, while 7e competitively blocked the activity of ENPP3. Cell viability analysis revealed the compound 4o as a cytotoxic agent against MCF7 (human breast cancer cell line) with a percentage inhibition of 63.2 ± 2.51%, whereas compounds 4c, 4d, 4n, and 7d decreased the HeLa cell viability (human cervical cancer cell line) to more than 50%. The tested compounds were non-cytotoxic against HEK293 (a human embryonic kidney cell line). Molecular docking analysis of selected inhibitors of both isozymes produced optimistic interactions with the influential amino acids, such as Leu290, Lys295, Tyr340, Asp376, His380, and Pro323 of ENPP1, whereas residues Asn226, His329, Leu239, Tyr289, Pro272, Tyr320, and Ala205 of ENPP3 crystallographic structure formed interactions with the potent inhibitors.
Aberrant level of ectonucleotide pyrophosphatase/phosphodiesterase-1 and -3 is linked with numerous disorders, for instance, diabetes, cancer, osteoarthritis, chondrocalcinosis, and allergic reactions. These disorders may be cured or minimized by blocking the activity of ENPP1 and ENPP3 isozymes. In this study, arylamide sulphonates were synthesized, characterized, and evaluated for their capability to affect the activity of isozymes ENPP1 and ENPP3. Among the selective inhibitors of ENPP1, compounds 4f and 4q exhibited sub-micromolar IC50 values of 0.28 ± 0.08 and 0.37 ± 0.03 μM, respectively, followed by 7a, with IC50 equal to 0.81 ± 0.05 μM, whereas out of the selective inhibitors of isozyme ENPP3, 4t and 7d preferably lessened the activity to half of the maximal inhibitory concentration of 0.15 ± 0.04 and 0.16 ± 0.01 μM alternatively. In addition, many structures including 4c, 4g, 4k, 4l, 4n, 4o, 4r, 4s, 7b, 7c, and 7e inhibited the activity of both isozymes to a significant level. Enzyme kinetic study of compound 4j revealed an uncompetitive mode of inhibition of ENPP1 isozyme, while 7e competitively blocked the activity of ENPP3. Cell viability analysis revealed the compound 4o as a cytotoxic agent against MCF7 (human breast cancer cell line) with a percentage inhibition of 63.2 ± 2.51%, whereas compounds 4c, 4d, 4n, and 7d decreased the HeLa cell viability (human cervical cancer cell line) to more than 50%. The tested compounds were non-cytotoxic against HEK293 (a human embryonic kidney cell line). Molecular docking analysis of selected inhibitors of both isozymes produced optimistic interactions with the influential amino acids, such as Leu290, Lys295, Tyr340, Asp376, His380, and Pro323 of ENPP1, whereas residues Asn226, His329, Leu239, Tyr289, Pro272, Tyr320, and Ala205 of ENPP3 crystallographic structure formed interactions with the potent inhibitors.
The members of the ENPP
family are conserved and ubiquitous eukaryotic
enzymes with an extracellular active site and exist as membrane-bounded
glycoproteins.[1] The ENPP family contains
seven fundamental isozymes (ENPP1–7), which are numbered as
per their cloning order.[2] ENPP1, ENPP2,
ENPP3, and ENPP4 can hydrolyze a variety of nucleotides together with
dinucleoside polyphosphates, nucleotide sugars, and cyclic (di-)nucleotides
and release corresponding nucleoside 5′-monophosphates; for
example, ATP is hydrolyzed into AMP and pyrophosphate (PPi).[3] However, among the nucleotide hydrolyzing potential
of isozymes (ENPP1, ENPP2, ENPP3, and ENPP4), the nucleotide metabolizing
activity of ENPP2 is weak. ENPP2 prefers phospholipids as a substrate,
whereas ENPP6 and ENPP7 are alkaline sphingomyelinases; hence, instead
of nucleotidases, the subtypes (ENPP2, ENPP6, and ENPP7) may be identified
as phospholipases.[4−8]ENPP1 was initially found as the surface sign for B-lymphocytes
that secrete antibodies; therefore, it was named plasma-cell differentiation
antigen-1 or PC-1.[9] ENPP1 exists in several
tissues such as in the bone (osteoblast) and cartilage (chondrocytes),
where it performs an essential part in the mineralization process.[10] ENPP1 is the transmembrane glycoprotein of homodimeric
type II, accompanied by an N-terminal transmembrane
domain, a catalytic domain, a C-terminal nuclease-like domain, and
two somatomedin-B-like domains.[11] The subcellular
location of ENPP1 is determined by the transmembrane domain, which
is essential for the dimerization of monomers by multiple disulfide
bonds. It is also noteworthy that the SMB2 domain of the ENPP1 is
assumed to be the residue for the insulin receptor interaction. The
ENPP1 catalytic domain amino acid residues make 24–60% likeness
among the various isozymes of other human ENPPs (ENPP2–7) and
similarity to the alkaline phosphatase family (APs).[12] Like the APs, two Zn2+ cations are strongly
linked by six preserved Asp/His residues at the active site. In addition,
the nuclease-like domain is connected to the catalytic domain via
a linker region called the “lasso loop”.[13] If any mutation occurs in this linker region,
the catalytic activity of the isozyme is eliminated. The nuclease-like
domain does not itself reveal catalytic activity, but it is required
to translocate ENPPs from the endoplasmic reticulum to the Golgi apparatus,
necessary for the correct folding of the ENPPs. Moreover, this domain
includes a perceived “EF hand” Ca2+ binding
motif, which is significant for the catalytic activity of ENPP1. Though
PPi is essential for the prevention of ectopic mineralization, its
over-supply contributes to the mineral aggregation of calcium pyrophosphate
dihydrate (CPPD) in joints. This situation is usually associated with
age-related osteoarthritis, known as chondrocalcinosis. The role of
ENPP1 was also found for insulin receptor signaling. The over-expression
of ENPP1 was found to be linked with defective insulin-stimulated
autophosphorylation in patients with type 2 diabetes. Moreover, ENPP1
was also identified and associated with tumor gradation in human astrocytic
brain tumors.[14]The isozyme ENPP3
(a glycoprotein) is associated with the plasma
membrane of the cell. ENPP3 (CD203c) is expressed on surfaces of various
body organs such as epithelial and mucosal, particularly on mast cells
and basophils.[15−17] ENPP3 is actively involved in fluid homeostasis,
modulation of bile formation, and cerebral spinal fluid secretion.
ENPP3 is actively involved against the synchronization of nucleotide
sugar glycosylation of brain-specific proteins and involves intracellular
as well as extracellular polyphosphate hydrolysis. The molecular structure
of ENPP3 contains parallel β-sheets of eight trapped catalytic
phosphodiesterase (PDE) domain in the molecular structure of ENPP3
that hold five active N-glycosylation sites and eight
α-helices on each side, representing a major component of alkaline
phosphatase superfamily.[12,18] This PDE domain is
linked to the nuclease-like region (NUC) through the linker L2, two
somatomedin B (SMB1 and SMB2 domains) with the L1 link attached to
the NUC.[19] An elevated proportion of ENPP3
on the cell surface was detected in comparison to other inflammatory
mediators, which stimulated basophils with antigen-bound IgE. ATP
emissions are thus impaired by ENPP3 and served as a sign of recognition
that basophils of patients have allergic responsiveness.[20]Several structures with sulphonate moiety
are reported to possess
therapeutic importance such as suramin (Figure ) with a complex structure and poly sulphonate
groups, being widely used as a positive control of choice in different
biological studies.[21] In the previous study,
we have reported pyrrolo[2,3-b]pyridine derivatives
such as (N-(pyrrolo[2,3-b]pyridine-1-carbonyl)benzenesulfonamide),
which was found as a selective inhibitor of ENPP1 with an IC50 value of 4.50 ± 0.16 μM and carbohydrazide-based derivative,
which selectively blocked the activity of ENPP3 to half of the maximal
value of 0.15 μM.[22,23] Various previously
reported structures with sulphonate scaffolds including raloxifene
sulphonates, benzofuran, and benzothiophene sulphonates exhibited
significant enzyme inhibitory activity to sub-micromolar levels, for
example, raloxifene sulphonate derivative 3-(4-(2-(piperidin-1-yl)ethoxy)benzoyl)-2-(4-(tosyloxy)phenyl)benzo[b]thiophen-6-yl-4-methylbenzenesulfonate tabulated selective
inhibition of isozyme ENPP1 to an IC50 value of 0.45 μM,
and benzothiophene derivative (4-(benzo[b]thiophen-5-yl)phenyl
cyclohexanesulfonate) was endowed with IC50 = 0.12 μM
against ENPP1 and 1.89 μM toward ENPP3.[24,25] We have also reported sulphonate derivatives mainly based on non-aromatic
or saturated cyclic hydrocarbons. These structures exhibited considerable
inhibitory activity but limited preferable selectivity among the three
isozymes ENPP1, ENPP2, and ENPP3.[26]
Figure 1
Previously
reported structures of ENPP inhibitors and the target
compounds.
Previously
reported structures of ENPP inhibitors and the target
compounds.Structure–activity relationship of the designed
compounds.The selected structural determinants in the target
compounds are
deemed unique scaffolds due to their importance in medicinal chemistry
(Figure ). The compounds
with carboxamide moiety have been reported as antimycobacterial agents
and antibacterial drugs.[27−31] The compounds bearing the sulphonate and sulphamoyl group have shown
anticancer, anti-hyperglycemic, and antimicrotubule properties.[32−35] Similarly, 4-aminophenol has been reported as a pharmacophore of
the drug “Paracetamol”.[36] Considering the importance of the structural scaffolds of the reported
structures such as the carboxamide moiety, sulphonate group, and 4-aminophenol,
we moved toward the synthesis of arylamide-based sulphonates. The
combination of various structural components referred to the compound
as multiple choices of selectivity toward isozymes ENPP1 and ENPP3.
Figure 2
Structure–activity relationship of the designed
compounds.
Results and Discussion
Chemistry
The synthetic pathways
to design arylamide sulphonate derivatives 4a–t and 7a–f are illustrated in Schemes and 2, respectively. The required sulphonate derivatives (4a–t) were synthesized by a 2-step synthetic pathway. In the first step,
the base-catalyzed nucleophilic-electrophilic substitution reaction
of 4-aminophenol (1) with 4-methyl benzoyl chloride (2a) or benzoyl chloride (2b) led to the formation
of structure 3a or 3b, respectively. In
the second step, sulphonate derivatives 4a–l were
produced by treating 3a with substituted benzene sulfonyl
chloride; similarly, sulphonate derivatives from 4m to 4t were obtained by reaction of 3b with substituted benzene
sulphonate chlorides under the same conditions. The second series
of sulphonate derivatives (7a–f) were synthesized
by a reaction of 1 with 9-methyl-9H-fluorene-9-carbonyl
chloride (5) to yield the phenolic intermediate 6, which was further used with various sulfonyl chlorides
to get the final products.
Scheme 1
Synthesis of Sulphonate Derivatives (4a–t)
Scheme 2
Synthesis of Sulphonate Derivatives (7a–f)
Enzyme Inhibition Assay
The enzymatic
inhibitory potential of the synthesized molecules was determined against
the hydrolytic activity of ENPP1 and ENPP3 isozymes. A number of compounds
manifested selective and significant inhibition of the activity of
both isozymes. The results of the enzymatic evaluation of the compounds 4a–t and 7a–f are depicted in Tables and 2, respectively.
Table 1
ENPP1 and ENPP3 Isozyme Inhibition
at 100 μM Compound Concentration, Presented in IC50 ± SEM (μM) Values or Percentage Inhibition ±SD
IC50 ± SEM (μM)a or % inhibition ± SDb
codes
R1
R2
ENPP1
ENPP3
4a
CH3
4-Me(C6H4)
26.2 ± 2.11b
17.0 ± 3.24b
4b
CH3
phenyl
45.2 ± 1.53b
29.8 ± 2.11a
4c
CH3
4-I(C6H4)
14.41 ± 0.54a
4.23 ± 0.38a
4d
CH3
4-OCH3(C6H4)
0.18 ± 0.01a
45.7 ± 2.35b
4e
CH3
4-Cl(C6H4)
44.8 ± 1.52b
0.82 ± 0.12a
4f
CH3
4-n-propyl(C6H4)
0.28 ± 0.08a
8.4 ± 3.11b
4g
CH3
4-OCF3(C6H4)
0.45 ± 0.07a
0.19 ± 0.02a
4h
CH3
biphenyl
1.73 ± 0.03a
27.1 ± 4.54b
4i
CH3
CH2(C6H4)
0.13 ± 0.01a
46.7 ± 1.54b
4j
CH3
4-n-butyl(C6H4)
0.23 ± 0.03a
42.1 ± 3.43b
4k
CH3
4-F(C6H4)
7.61 ± 0.02a
0.37 ± 0.08a
4l
CH3
2-F(C6H4)
0.15 ± 0.01a
3.48 ± 0.34a
4m
H
4-OCH3(C6H4)
0.82 ± 0.05a
22.6 ± 2.12b
4n
H
CH2(C6H4)
3.96 ± 0.12a
4.12 ± 0.05a
4o
H
4-n-propyl(C6H4)
1.71 ± 0.04a
2.12 ± 0.08a
4p
H
biphenyl
1.51 ± 0.02a
33.5 ± 1.92b
4q
H
8-quinolinyl
0.37 ± 0.03a
19.5 ± 3.53b
4r
H
4-I(C6H4)
2.01 ± 0.34a
0.53 ± 0.19a
4s
H
4-F(C6H4)
5.76 ± 0.57a
2.95 ± 0.25a
4t
H
4-OCF3(C6H4)
41.4 ± 3.24b
0.15 ± 0.04a
Suramin
7.80 ± 0.09a
0.89 ± 0.16a
The results are obtained by means
of a triplicate assay; IC50 values are calculated for compounds
showing inhibition of enzyme activity to more than 50%; SEM = standard
error mean.
Percentage inhibition
of the compounds
exhibiting <50% inhibition; SD = standard deviation.
Table 2
ENPP1 and ENPP3 Isozyme Inhibition
at 100 μM Compound Concentration, Presented in IC50 ± SEM (μM) Values or %Age Inhibition ±SD
IC50 ± SEM (μM)a or % inhibition ± SDb
codes
R
ENPP1
ENPP3
7a
2-F(C6H4)
0.81 ± 0.05a
28.2 ± 3.22b
7b
CF3
0.90 ± 0.16a
4.16 ± 0.22a
7c
4-F(C6H4)
0.50 ± 0.03a
6.69 ± 1.89a
7d
methyl
9.2 ± 4.13b
0.16 ± 0.01a
7e
4-CH3(C6H4)
18.87 ± 0.54a
0.38 ± 0.13a
7f
phenyl
8.45 ± 0.25a
42.1 ± 2.44b
The results are obtained by means
of a triplicate assay; IC50 values are calculated for compounds
showing inhibition of enzyme activity to more than 50%; SEM = standard
error mean.
Percentage inhibition
of the compounds
exhibiting <50% inhibition; SD = standard deviation.
The results are obtained by means
of a triplicate assay; IC50 values are calculated for compounds
showing inhibition of enzyme activity to more than 50%; SEM = standard
error mean.Percentage inhibition
of the compounds
exhibiting <50% inhibition; SD = standard deviation.The results are obtained by means
of a triplicate assay; IC50 values are calculated for compounds
showing inhibition of enzyme activity to more than 50%; SEM = standard
error mean.Percentage inhibition
of the compounds
exhibiting <50% inhibition; SD = standard deviation.
Study of Structure–Activity Relationship
of Compounds Measured through Enzyme Inhibition Assay
The
effect of different structural variations at position “R1, R2” for compounds 4a–t and 7a–f was studied in terms of the structure–activity
relationship. Comparison of results of compounds 4a, 4b, 4f, and 4j by the structures
indexed interesting results. Compound 4a with 4-methyl
phenyl linkage inhibited both isozymes’ activity to a lesser
extent, merely 26.2% in the case of ENPP1 and only 17.0% for ENPP3.
These results were slightly modified when the p-methyl
group was removed with the phenyl ring only (4b). Again
structure 4b was devoid of ENPP1 inhibitory potential
but blocked the activity of ENPP3 to an IC50 value of 29.76
± 2.13 μM. Modification of phenyl ring substitution to
4-methoxy phenyl group (4d) inversed the inhibitory potential
more selectively to ENPP1 (IC50 = 0.18 ± 0.01 μM)
and blocked the ENPP3 isozyme to only 45.7%.In compliance with
the modification at phenyl ring with methyl linkage, elongation of
this linkage to three carbons atoms (4f) and four carbon
atoms (4j) significantly and selectively turned the activity
toward ENPP1, that is, compound 4f with a propyl linkage
blocked the ENPP1 to IC50 = 0.28 ± 0.08 μM,
whereas 4j with a butyl linker inhibited the isozyme
ENPP1 with an IC50 value of 0.23 ± 0.03 μM.
These three structures with hydrocarbon linkers (4d, 4f, 4j) poorly blocked the hydrolytic activity
of ENPP3 to very low or nearly 50% percent inhibition values such
as 45.7, 8.4, and 42.1%, respectively. The outcomes of 4d, 4f, and 4j suggested that the methoxy
group or long hydrocarbon side chain with the phenyl ring was favorable
to ENPP1 inhibition. A study of SAR for compounds 4b and 4h showed selective but opposite behavior toward isozymes
ENPP1 and ENPP3. The molecule with a single phenyl ring attachment
at position “R2” possessed the selective
potential to ENPP3, whereas the biphenyl attachment at the same position
turned the selectivity more specifically to isozyme ENPP1. The structures
possessing halogens with phenyl moiety produced variable but promising
outcomes. Comparison of SAR for 4c, 4e,
and 4k substituted with 4-iodophenyl, 4-chlorophenyl,
and 4-fluorophenyl at the “R1” position showed
almost non-selective behavior without significant inhibitory difference
among the two isozymes, except inhibitor 4e, which exhibited
comparatively more affinity for ENPP3 with IC50 equal to
0.82 ± 0.12 μM. The relative study of SAR of 4k and 4l revealed that the presence of fluorine atom
at the para position was much favorable for ENPP3, whereas fluorine
atom at the meta position shifted the potency toward ENPP1.The presence of the methoxy group on the phenyl ring rendered compound 4m toward the selectivity of ENPP1 (IC50 = 0.82
± 0.05 μM). The same selectivity pattern was observed in
the case of 4p with biphenyl attachment at the R2 position, reducing the activity of ENPP1 to half of the maximal
inhibitory concentration (IC50) of 1.51 ± 0.02 μM
and only 33.5% inhibition to ENPP3 isozyme. This pattern of selectivity
suggested the −OCH3 or phenyl group linkage at the
phenyl ring, which was favorable for ENPP1 inhibition. Elongation
of the side chain to three carbon atoms (4o) abolished
the selectivity pattern for either of the isozymes. A careful comparative
study of structures 4f and 4o highlighted
the point that the n-propyl linkage is much beneficial
toward the selective inhibition of ENPP1 activity in the presence
of methyl group on the other side of the structure. The presence of
8-quinolinyl moiety produced the most potent and selective structure
(4q) that blocked the ENPP1 activity to IC50 value equal to 0.37 ± 0.03 μM, suggesting the heterocyclic
attachment favorable for the inhibition of ENPP1 isozyme activity.
The halogenated structures 4r–t revealed compound 4t possessing −OCF3 as a selective inhibitor
for ENPP3 (0.15 ± 0.04 μM), contributing to the effect
of −OCF3 toward inhibition of ENPP3 isozyme activity. Moreover,
the methoxy group substitution was found preferable for only ENPP1
(4d, 4m), and the results of 4h and 4p demonstrated selective behavior toward ENPP1,
suggesting the importance of the biphenyl group for ENPP1 inhibition
(Table ).In
the next series of compounds (7a–f), aromatic
and aliphatic sulfonyl chlorides were used to produce complex but
promising structures. Interestingly, among the six compounds (7a–f), two compounds selectively blocked the activity
of ENPP1, 7a with 2-fluorophenyl substitution (IC50 = 0.81 ± 0.05 μM) and 7f (IC50 = 8.45 ± 0.25 μM) possessing a single phenyl
ring at position “R”, and the molecule 7d carrying a methyl group preferably inhibited the isozyme ENPP3 illustrating
an IC50 value of 0.16 ± 0.01 μM. The trifluoromethane
substitution produced structure 7b, which presented dual
inhibition of both isozymes (ENPP1, IC50 = 0.90 ±
0.16 μM; ENPP3, IC50 = 4.16 ± 0.22 μM),
whereas a single fluorine atom as a 4-fluorophenyl attachment preferably
blocked the ENPP1 activity over ENPP3.
Pattern of Enzyme Inhibition of Compounds
4j and 7e
The pattern of enzyme inhibition for a single active
compound against each isozyme was assessed by constructing the Lineweaver–Burk
graphs (Figure ).
The mode of inhibition for the comparatively selective and potent
inhibitor of ENPP1 (4j) and ENPP3 (7e) was
determined by performing enzyme kinetic studies. Compound 4j blocked ENPP1 uncompetitively, whereas 7e competitively
inhibited ENPP3 activity.
Figure 3
Presentation of Lineweaver–Burk plots
of enzyme kinetics
of 4j (left) and 7e (right). The experiments
were performed in triplicates.
Presentation of Lineweaver–Burk plots
of enzyme kinetics
of 4j (left) and 7e (right). The experiments
were performed in triplicates.
Cell Cytotoxicity Analysis (MTT Assay)
The MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide]
assay was executed to find out the cytotoxic effect of the compounds
on human breast cancer cell line (MCF7), human cervical cancer cell
line (HeLa), and human embryonic kidney cell line (HEK293). The compounds
were found cytotoxic against cancer cell lines, but none of the molecules
were toxic on normal cell lines. The results of the cytotoxicity assay
are tabulated in Table . The cytotoxic potential of the compounds was evaluated on the three
cancerous cell lines: (1) MCF7, a human breast cancer cell line; (2)
HeLa, a human cervical cancer cell line; and (3) HEK293, a human embryonic
kidney cell line. Neither of the tested compounds manifested a significant
effect on MCF7 cancer cell viability except for 4o with
63.2% inhibition checked at 100 μM compound concentration. The
molecules 4c, 4d, 4n, and 7d decreased the viability of HeLa cancer cells to above 50%.
None of the compounds affected the viability of normal cell lines
to a significant level.
Table 3
Cytotoxicity of 100 μM Compound
on MCF7, HeLa, and HEK293 Cell Lines
codes
MCF7 (% inhibition ±SD)a
HeLa cell
line (% inhibition ±SD)a
HEK293 (% inhibition ±SD)a
4a
1.1 ± 0.12
16.4 ± 0.92
5.3 ± 0.81
4b
7.4 ± 0.04
44.9 ± 1.81
6.4 ± 0.22
4c
33.2 ± 1.32
57.8 ± 2.72
15.2 ± 0.83
4d
31.5 ± 2.43
53.4 ± 1.93
13.1 ± 0.73
4e
24.2 ± 1.73
32.4 ± 3.82
8.3 ± 0.35
4f
25.7 ± 0.92
49.2 ± 2.12
4.2 ± 0.12
4g
11.7 ± 0.84
17.7 ± 0.72
3.2 ± 0.11
4h
9.1 ± 0.71
48.2 ± 1.43
4.8 ± 0.24
4i
5.7 ± 0.06
38.5 ± 1.14
11.2 ± 0.92
4j
13.3 ± 0.73
19.1 ± 1.11
27.6 ± 2.55
4k
11.2 ± 0.52
7.7 ± 0.73
2.1 ± 1.22
4l
0.5 ± 0.03
14.3 ± 3.52
17.4 ± 0.83
4m
4.1 ± 0.01
7.1 ± 0.81
2.3 ± 1.12
4n
7.3 ± 0.03
51.3 ± 0.82
5.8 ± 0.23
4o
63.2 ± 2.51
38.7 ± 2.43
12.6 ± 0.74
4p
8.4 ± 0.75
5.3 ± 0.07
4.3 ± 0.15
4q
11.2 ± 1.04
24.2 ± 1.44
4.6 ± 0.32
4r
10.2 ± 0.86
12.6 ± 0.73
11.3 ± 0.71
4s
2.1 ± 0.08
16.5 ± 1.12
20.6 ± 1.14
4t
29.6 ± 1.22
49.7 ± 2.53
8.6 ± 0.14
7a
12.7 ± 0.62
10.7 ± 0.72
14.6 ± 1.12
7b
14.1 ± 0.63
32.4 ± 2.63
24.5 ± 1.16
7c
12.4 ± 1.22
23.5 ± 1.64
9.2 ± 0.42
7d
0.2 ± 0.06
56.3 ± 1.17
5.5 ± 0.36
7e
0.3 ± 0.15
33.8 ± 1.56
12.4 ± 1.42
7f
5.6 ± 1.11
28.1 ± 1.33
6.2 ± 1.17
cisplatin
N.D
84.3 ± 2.72
65.5 ± 1.48
doxorubicin
88.8 ± 2.74
N.D
N.D
Cytotoxicity of compounds at 100
μM, calculated by taking a mean of the triplicate assays; %
inhibition ± standard deviation; ND = not determined.
Cytotoxicity of compounds at 100
μM, calculated by taking a mean of the triplicate assays; %
inhibition ± standard deviation; ND = not determined.
Molecular Docking Studies
The putative
mode of structural component’s interaction with the amino acid
residues of ENPP1 and ENPP3 isozymes was found by applying molecular
docking analysis. The in silico findings produce reliable outcomes
following the in vitro results. The 2D and 3D modes of interaction
are shown in Figures –8.
Figure 4
Illustration of 3D plausible binding modes for
selective inhibitor 4f with ENPP1 isozyme.
Figure 8
Presentation of 3D binding interaction of dual inhibitor 4g; (a) 3D binding modes with ENPP1 and (b) 3D binding modes
with an ENPP3 crystallographic structure.
Illustration of 3D plausible binding modes for
selective inhibitor 4f with ENPP1 isozyme.Illustration of 3D binding interactions of selective inhibitor 7a with residues of ENPP1 isozyme.Presentation of 3D binding modes of inhibitor 4t with
ENPP3 isozyme residues.Illustration of 3D binding modes of inhibitor 7d with
ENPP3 isozyme residues.Presentation of 3D binding interaction of dual inhibitor 4g; (a) 3D binding modes with ENPP1 and (b) 3D binding modes
with an ENPP3 crystallographic structure.
Molecular Docking Studies at the Crystal
Structure of ENPP1 Isozyme
The most active synthesized compounds
were consigned to molecular docking analysis to elaborate the inhibitory
results against the ENPP1 isozyme. The reported crystallographic structure
of human ENPP1, docked with reference ligand “N-{4-[(7-methoxyquinolin-4-yl)oxy]phenyl}sulfuric diamide”
(TZV), was downloaded from Protein Data Bank (PDB) with PDB ID = 6wew.[37] The bound ligand TZV was found in contact with several
amino acid residues of the catalytic site; such as the quinoline nucleus
was stacked between Tyr340 and Phe257 and through π–π
T-shaped interaction with Pro323 and Leu290. Moreover, TZV showed
H-bond interaction with Lys338, Lys295, Asp376, and His380. A π-alkyl
and π-anion connection was found with Leu290, Pro323, and Asp218.[37] The most selective and relatively potent inhibitor
of ENPP1 4f showed H-bonding linkage between the oxygen
atom of sulphonate group and amino acid residue Thr256. The inhibitor 4f rested itself through an array of π-linkages, that
is, π–π stacked, π–π T-shaped,
amide-π stacked, and π-alkyl bonds with residues His380,
Tyr371, Tyr340, Trp322, and Leu290. Moreover, 4f was
coordinated to zinc metal ion, Asp376, through π-anion and π-cation
bonds (Figure ).The docking studies of another selective inhibitor of isozyme ENPP1
(7a) manifested a wide range of occupancy of amino acids,
attributing to its selectivity. Residues Lys255, Asn277, and Thr256
coordinated through hydrogen bonds, predominately with sulphonate
group of 7a. Similar to 4f, compound 7a stabilized itself via π-connections such as a π-anion
link to Asp376, π–π stacked connection Tyr340,
π-alkyl interaction to residues Leu290, and Tyr340. Moreover,
Zn2+ chelation connected to the oxygen atom of sulphonate
moiety was observed as depicted in Figure .
Figure 5
Illustration of 3D binding interactions of selective inhibitor 7a with residues of ENPP1 isozyme.
Molecular Docking Studies at the Crystal
Structure of ENPP3 Isozyme
The selective inhibitor of ENPP3 4t stabilized itself through an array of several amino acid
residues. H-bond interactions were observed between oxygen atoms of
sulphonate group and residues Ala205, Asn226, and Lys204 and hydrogen
of amide moiety and residue Asn482. Several π-interactions contributed
to the stability and activity of 4t in the active site.
These linkages include π-anion connection with Asp167: π-alkyl
stick bridging Ala205; Leu239: π–π stacked connection
with His329: π-sulfur attachment with His483: π-anion
linkage with His483; and Asp167: π-lone pair link among the
lone pair of asparagine residue (Asn477) and phenyl ring of phenyl
acetamide moiety. A π-cation linkage was noticed between the
zinc ion and sulphonate group oxygen. More noticeably, an unfavorable
bond was observed among Asp325 and oxygen atoms of sulphonate moiety,
contributing to the high affinity and stability of the molecule 4t (Figure ).
Figure 6
Presentation of 3D binding modes of inhibitor 4t with
ENPP3 isozyme residues.
The molecular docking study of another selective structure
for ENPP3 (7d) revealed promising interactions with amino
acid residues. The structure 7d balanced itself through
hydrogen bonding of sulphonate group oxygens with Asn226 and His329.
The three aromatic rings and the methyl group connected to the heterocyclic
scaffold formed an array of π–π stacked and π-alkyl
linkages with residues Tyr289, Pro272, Leu239, Tyr320, and Ala205.
In addition, one π-anion interaction was found between the phenyl
ring of phenyl methanesulfonate attachment and Asp325 as shown in Figure .
Figure 7
Illustration of 3D binding modes of inhibitor 7d with
ENPP3 isozyme residues.
Molecular Docking Studies at the Crystal
Structure of ENPP1 and ENPP3 Isozymes
The observation of
the plausible binding mode of the relatively potent but dual inhibitor 4g gave insights into the binary nature or non-selective behavior
of the structure. In the case of ENPP1, inhibitor 4g formed
an array of π-connections; Tyr371 (π–π stacked,
π-sigma bond), Tyr340 (π–π T shaped), Trp322
(π–π stacked), and Pro323 and Phe321 (π-alkyl
linkage). Similar to inhibitors 4f and 7a, compound 4g was also chelated with zinc metal ions. 4g blocked the catalytic activity of ENPP3 to a low micromolar
level, tabulated as the second active inhibitor to ENPP3 (0.19 ±
0.02 μM). The binding modes gave the results in obedience to
experimental data. As a result of docking into the ENPP3 crystalline
structure, 4g was linked via four H-bonds: two with Lys204
and Asn226 through the oxygen of sulphonate scaffold and two with
Asn482 and His329. The whole structure was surrounded by Pi-linkages of several amino acid residues such as π–π
bonds with residues Tyr289, Tyr320, Phe206, and His329 and π-alkyl
and alkyl bonds with Leu239, Ala205, and Pro272. In addition, the
important connections were observed with His483 (π-lone pair
bond) and Asp167 (π-anion bond). Zinc ion chelation was also
observed through the oxygen atom of the sulphonate group (Figure ).
Conclusions
A new series of compounds
possessing specifically sulphonate group
was synthesized, characterized, and examined for the affinity and
preference to block the activity of ENPP1 and ENPP3 isozymes of the
ectonucleotidase pyrophosphatase/phosphodiesterase (ENPP) family.
Among the series, 4a–t compounds 4d, 4f, 4h, 4i, 4j, 4m, 4p, and 4q selectively
blocked the activity of ENPP1 to less or high preferentially, whereas
in the group of structures 7a–f, structure 7a and 7f preferably inhibited the activity of
ENPP1. Among both series, compounds 4f, 4h, 4m, 4p, 4q, and 7a were found as more potent and particular inhibitors of ENPP1 with
IC50 values of 0.28 ± 0.08, 1.73 ± 0.03, 0.82
± 0.05, 1.51 ± 0.02, 0.37 ± 0.03, and 0.81 ± 0.05
μM, alternatively. Similarly, compounds 4b, 4e, 4t, and 7d showed discriminatory
selective behavior for ENPP3 with most active structures 4t (0.15 ± 0.04 μM) and 7d (IC50 = 0.16 ± 0.01 μM). In the same way, compounds 4c, 4g, 4k, 4l, 4n, 4o, 4r, 4s, 7b, 7c, and 7e exhibited dual inhibitory
behavior to ENPP1 and ENPP3, particularly 4g as a more
potent structure that tabulated IC50 = 0.45 ± 0.07
μM against ENPP1 and IC50 = 0.19 ± 0.02 μM
toward ENPP3. The cancer cell line (MCF7, HeLa) and normal cell line
(HEK293) viability assay revealed that compounds 4c, 4d, 4n, and 7d were cytotoxic against
HeLa cancer cell lines. The structure 4o significantly
decreased the MCF7 cancer cell line viability to 63.2%; however, none
of the molecules showed cytotoxicity toward normal cell lines and
placed the compounds within a safety profile. The plausible pattern
of interactions of selective and non-selective structures with salient
amino acids comprising both proteins ENPP1 and ENPP3 complied with
in vitro results. In summation, we found several selective and non-cytotoxic
inhibitors for both isozymes ENPP1 and ENPP3. These compounds may
be further evaluated as therapeutic agents for the ailments associated
with ENPP1 and ENPP3 isozymes.
Experimental Section
Reagents, Solvents, and Apparatus
The chemicals and solvents used in this study were obtained from
commercial sources. Chemicals were used as such, whereas solvents
were dried or distilled where required. 1H NMR and 13C NMR spectra were obtained by a Bruker spectrometer (AV
250, AV 300, and AV 500) in solvents CDCL3 and DMSO. FTIR spectra
were detected by an ATR apparatus. The high resonance mass spectra
(HRMS) ESI was detected on a device Finnigan MAT 95 XP with an HP-5
capillary column using helium as carrier gas (Thermo Electron Corporation).
The melting point of the compounds was obtained with the use of the
Büchi apparatus. Column chromatography was carried out by using
silica gel 60 A with a mesh size of 60–200. Analytical thin-layer
chromatography was performed on silica gel plates (0.20 mm, 60 A),
and spots were detected through a UV absorbance lamp 254 nm/366 nm.
General Procedure of Synthesis
Synthesis of N-(4-Hydroxyphenyl)-4-methylbenzamide
(3a) and N-(4-Hydroxyphenyl)benzamide
(3b)
The synthesis of compounds 3a and 3b was carried out by treating 4-aminophenol (1) (4.12 mmol; 0.450 g) with 4-methyl benzoyl chloride (2a) (2.50 mmol; 0.386 g) or benzoyl chloride (2b) (2.50 mmol; 0.350 g; 0.423 mL), respectively. The reactants were
combined with stirring at 0 °C in acetone (50 mL) in the presence
of anhydrous potassium carbonate (3.62 mmol; 0.500 g). After that,
the reaction mixture was stirred for 4 h at room temperature and filtered.
The filtrate proceeded to dryness, the residue was dissolved in ethyl
acetate (40 mL), and extraction was carried out by using dilute HCl
(30 mL). The organic layer was separated and washing was performed
with a saline solution (2 × 20 mL). The separated organic layer
was dried of water molecules by using anhydrous sodium sulfate (Na2SO4) and filtered. The organic solvent was evaporated
to dryness, and obtained compounds were used for the next step.
Synthesis of Arylamide Sulphonates (4a–l, 4m–t)
The targeted
compounds 4a–l were synthesized by the reaction
of 3a (0.485 mmol; 0.110 g) with various benzene sulfonyl
chloride derivatives (0.985 mmol), and compounds from 4m–t were obtained by treating 3b (0.485 mmol; 0.103 g)
with appropriate benzene sulfonyl derivatives. The separate reactions
were performed in dry tetrahydrofuran (THF) (12 mL), and triethylamine
(2.10 mmol) was added dropwise to the reaction mixture at 0 °C.
Each reaction mixture was stirred at room temperature till completion
and quenched up to three times with water and ethyl acetate in equal
volumes (3 × 20 mL). The organic layer was dried with anhydrous
Na2SO4 followed by two times washing with a
saline solution (2 × 30 mL), filtered, and subjected to dryness
on the rotary evaporator. The crude residue was purified by recrystallization
with absolute ethanol.
Synthesis of N-(4-Hydroxyphenyl)-9-methyl-9H-fluorene-9-carboxamide (6)
Compound 6 was synthesized by the reaction of 4-aminophenol (1) (4.12 mmol; 0.450 g) with 9-methyl-9H-fluorene-9-carbonyl
chloride (5) (2.50 mmol; 0.605 g) in acetone (50 mL)
and base anhydrous potassium carbonate (3.62 mmol; 0.500 g). The reaction
proceeded as per conditions followed for the synthesis of 3a and 3b.
Synthesis of Arylamide Sulphonates 7a–F
The desired sulphonate derivatives 7a–f were produced by the reaction of N-(4-hydroxyphenyl)-9-methyl-9H-fluorene-9-carboxamide (6) (0.485 mmol; 0.152 g) with alkyl or aryl sulfonyl chloride derivatives
(0.980 mmol) in dry tetrahydrofuran (THF) (12 mL). Triethylamine (2.10
mmol) was added dropwise to the stirring mixture at 0 °C. The
reaction proceeded as per the protocol used for the synthesis of compounds 4a–l and 4m–t.
Transfection
was done in COS-7 cells using plasmids expressing human ENPP1 (GenBank
accession no. NM_006208) or ENPP3 (GenBank accession no. NM_005021) and
“Lipofectamine” as transfection reagent.[38,39] In a 10 cm culture dish plate, COS-7 cells (70–80% confluent
cells) were incubated with DMEM/F-12 (lacking FBS) including plasmid
DNA (6 μg) and 24 μL of Lipofectamine reagent and incubated
for 5 h at 37 °C; then an equal volume of DMEM/F-12 supplemented
with 20% FBS was added to stop the process of transfection. The transfected
cells were isolated from the media after 48–72 h.[40]
Preparation of Protein Aliquots
The transfected cells were washed with cold Tris-saline buffer (4
°C) and collected in the presence of a harvesting buffer that
was composed of 95 mM NaCl and 0.1 mM phenylmethylsulfonyl fluoride
(PMSF), 45 mM Trisbase, pH 7.5. The scrapped cells were subjected
to washing through centrifugation twice at 300g for
5 min at 4 °C.[30] The cells were sonicated
after re-suspension in the harvesting buffer with aprotinin (10 μg/mL)
and centrifuged at 850g at 4 °C for 5 min. The
supernatant was collected, and aliquots were preserved at −80
°C after adding glycerol (7.5%). For estimation of protein, Bradford
microplate assay was applied using bovine serum albumin as ref (41).
Enzyme Inhibition Assay
The effect
of synthesized molecules on the activity of isozymes ENPP1 and ENPP3
was evaluated by using the previously reported method with minor changes.[22,38] Initial screening of the compounds was performed at 100 μM
concentration per well in triplicate in the buffer: 5 mM MgCl2, 0.1 mM ZnCl2, and 50 mM Tris–HCl, adjusted
to pH 9.5–9.6. Briefly, 30 ng ENPP1 and 32 ng ENPP3 proteins
per well were added followed by the addition of 100 μM synthesized
compounds dissolved in 10% dimethyl sulfoxide (DMSO). The reaction
was initiated by adding a substrate (thymidine 5′-monophosphate
para-nitrophenyl ester) 300 μM per well and allowed to incubate
for 35 min at 37 °C. The reading was taken with the help of a
microplate reader (BioTek FLx800, Instruments, Inc., USA) at a wavelength
of 405 nm. The compounds presenting more than 50% inhibition of either
isozyme were subjected to serial dilutions, and the IC50 values were calculated through the non-linear regression analysis
curve fitting program PRISM 5.0 (GraphPad, San Diego, CA, USA).
Enzyme Kinetic Study
The enzyme
kinetic study was performed by using a serial concentration of substrate
0, 1.25, 2.5, 5.0, 7.5, and 10 mM. The compound concentrations for 4j were 0.0, 0.5, 1.0, and 1.5 μM, and for inhibitor 4e, the used concentrations were 0.0, 0.3, 0.6, and 1.2 μM.
The assay plate was prepared by adding the concentrations of enzymes
(ENPP1 = 30 ng, ENPP3 = 32 ng) per well to the assay buffer. After
adding the serial concentrations of compounds and substrates, the
assay mixture was incubated at 37 °C, and the reading was noted
for up to 35 min with 5 min intervals by using a microplate reader
(BioTek FLx800, Instruments, Inc., USA). The data were analyzed by
PRISM 5.0 (Graph Pad, San Diego, CA, USA), and Lineweaver–Burk
graphs were plotted.
Cytotoxicity Assay (MTT Assay)
The MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide]
assay was applied to check the cytotoxic effect of the compounds with
minor modifications in the already reported method.[42] The cells were seeded in 96-well plates at a quantity of
about 4×104 cells per well and incubated for 24 h
in a CO2 incubator. After incubation for 24 h, serum-free
media were added to replace the existing media of cell culture, and
the cells were treated with test compounds along with positive control
at 100 μM concentrations per well. Afterward, the incubation
was given for the next 24 h; MTT reagent (0.5 mg/mL) was added and
further incubated at 37 °C for 4 h. Dimethylsulfoxide was added
to solubilize the violet crystals and a stopping reagent composed
of isopropanol (50%) and sodium dodecyl sulfate (10%) was added. The
absorbance was measured at 570 nm with the help of a microplate reader
(FLUOstar Omega Microplate Reader, Germany), and results were compiled.
Molecular Docking Analysis
For
the in silico observations, the available ligand-bounded human crystallographic
structures of ENPP1 (PDB ID = 6wew) and ENPP3 (PDB ID = 6C02) were downloaded
from PDB.[37,43] The minimization of energy and molecule
database of the docked molecules was generated with the help of the
Molecular Operating Environment.[44] The
protein preparation such as energy minimization, selection of pocket
site, removal of water molecules, protonation, and docking studies
was performed by using the default system of LeadIT (BioSolveIT GmbH,
Germany).[45] Among the selected number of
poses for each compound, poses that possessed low binding energy and
considerable affinity with the amino acids residues were selected
for further visualization by using Discovery Studio Visualizer DS.[46]
Authors: Mohammad H Semreen; Mohammed I El-Gamal; Saif Ullah; Saquib Jalil; Sumera Zaib; Hanan S Anbar; Joanna Lecka; Jean Sévigny; Jamshed Iqbal Journal: Bioorg Med Chem Date: 2019-04-26 Impact factor: 3.641
Authors: Saif Ullah; Mohammed I El-Gamal; Randa El-Gamal; Julie Pelletier; Jean Sévigny; Mahmoud K Shehata; Hanan S Anbar; Jamshed Iqbal Journal: Eur J Med Chem Date: 2021-03-10 Impact factor: 6.514
Authors: William N Addison; Fereshteh Azari; Esben S Sørensen; Mari T Kaartinen; Marc D McKee Journal: J Biol Chem Date: 2007-03-23 Impact factor: 5.157
Authors: Anne Zuse; Peter Schmidt; Silke Baasner; Konrad J Böhm; Klaus Müller; Matthias Gerlach; Eckhard G Günther; Eberhard Unger; Helge Prinz Journal: J Med Chem Date: 2007-10-31 Impact factor: 7.446