1272
M. Wos et al. / Biomedicine & Pharmacotherapy 93 (2017) 1269–1276
Table 3
The IC50 and SI of 1, 4 derivatives and Cisplatin on cancer cell lines (A549, HepG2, MCF-7) and normal cell line (BJ) after 96-h incubation. The IC50 values were expressed as
mean values ꢁ SD from three separate experiments. IC50 ꢀ inhibitory concentration that caused 50% inhibition of cell proliferation; SI ꢀ the ratio between IC50 on normal cell
line and IC50 on cancer cell line.
Comp.
IC50 ꢁ SD [
m
g/mL]
SI
BJ
A549
HepG2
MCF-7
BJ/
BJ/
BJ/
A549
HepG2
MCF-7
1
4
2.87 ꢁ 1.10
3.09 ꢁ 1.19
1.19 ꢁ 0.56
4.96 ꢁ 1.96
6.80 ꢁ 1.25
2.38 ꢁ 1.16
3.64 ꢁ 1.18
11.88ꢁ 3.91
2.82 ꢁ 1.79
2.53 ꢁ 1.21
3.62 ꢁ 1.07
2.97 ꢁ 1.33
0.58
0.45
0.5
0.79
0.26
0.42
1.13
0.85
0.40
Cisplatin
2.6. Molecular modeling
2.6.1. Homology modeling
compounds (1–9) were prepared in reaction of carboxylic acid
hydrazide with isothiocyanate using methanol as solvent
a
(Scheme 1). All the synthesized compounds were characterized
by spectroscopic methods. 1H NMR spectra of thiosemicarbazide
derivatives exhibited proton signals typical of the NH group in the
range of 7.98–11.93 ppm. Also, in the IR spectra the signals from the
Homology modeling was applied to construct molecular model
a-glucosidase from Saccharomyces cerevisiae (UniProt entry
of
P53341) in complex with glucose. Sequence alignment was
performed using Muscle software [46]. Homology modeling was
performed using Modeller v. 9.17 [47]. X-ray structure of
group C .
O and C = S are observed in the range 1600–1800 cmꢀ1
¼
The above observations suggest that the compounds exist mainly
in the thione form. This is consistent with the literature data [9,11].
a
-glucosidase from Geobacillus sp. HTA-462 (PDB ID: 2ZE0 [48])
was used as a template (sequence identity 40%). The glucose
moiety was transferred to the enzyme catalytic site from the
sucrose mutase-sucrose complex (PDB ID: 2PWE [49]). A popula-
tion of 100 homology models was generated and subsequently
assessed by the Modeller objective function and Discrete Opti-
mized Protein Energy profiles [50]. The best model was subjected
to quality assessments using the Schrödinger software module for
Ramachandran plots.
For derivative 4 as a model compound X-ray structure
determination was carried out It was performed in order to
confirm the synthesis pathway, assumed molecular structure of 4
and identification of its tautomeric form in the crystalline phase.
This compound crystallizes with the methanol solvent molecule
lying on special position with twofold rotation symmetry. The
structure and conformation of the molecule 4 in the crystal is
shown in Fig. 1. In the crystalline state, the molecule exists in N1-
amino/N3-amino/N4-amino/S2-thione/O5-keto tautomeric form,
as evidenced by the C2–S2 and C5–O5 bond length of 1.213(3) and
1.684(2) Å, respectively, typical for the carbonyl and thione group
[55]. Moreover, the positions of the amino-H atom are localized at
the difference electron-density map in the vicinity of N1, N3 and
N4 atoms. The carbonylthiosemicarbazide part of the molecule
adopt trans-cis-gauche-trans conformation with the C11–N1–C2–
N3, N1–C2–N3–N4, C2–N3–N4–C5 and N3–N4–C5–C51 of 174.2(2),
ꢀ8.6(3), 103.9(2) and 171.73(18)ꢃ, respectively. The torsion angle
O5–C5–C51–C56 of ꢀ20.6(4)ꢃ shows that the carbonyl group is
slightly distorted from coplanarity with the benzene ring, similarly
as disordered nitro group with the torsion angles C55–C54–N6–
O61 of 6.0(5)ꢃ, C55–C54–N6–O62A of ꢀ152.8(4)ꢃ and C55–C54–
N6–O62 B of 152.9(7)ꢃ. The thione C2 = S2 group has trans and cis
conformation with respect to N3–N4 and N1–C11 bonds,
respectively; the respective torsion angles are 172.38(16) and
2.6.2. Compound preparation
The investigated compounds were modeled using the LigPrep
protocol from the Schrödinger Suite [51]. In order to sample
different protonation states of ligands in physiological pH, Epik
module was used with default settings [52].
2.6.3. Molecular docking
Molecular docking was performed using Glide from the
Schrödinger suite of software [53]. The grid file was generated
indicating glucose as a reference ligand. Molecular docking was
performed using the SP (standard precision) protocol of Glide. 100
poses were generated for each ligand. PyMol v. 0.99 [54] was used
for visualization of results.
3. Results and discussion
3.1. Chemistry
The starting reagents for the synthesis of thiosemicarbazides 1–
9 were 2-, 3-, 4-pyridyl carboxylic acid hydrazide or 4-nitrophenyl
carboxylic acid hydrazide. These compounds were obtained by the
reaction of appropriate ester with hydrazine hydrate in anhydrous
ethanol by the procedure described earlier [31]. The title
Table 4
Probability to inhibit (Pa) and not to inhibit (Pi)
software PASS.
a-glucosidase calculated using
Compound
Pa
Pi
1
2
3
4
5
6
7
8
9
0.919
0.901
0.915
0.886
0.777
0.927
0.926
0.915
0.900
0.009
0.012
0.010
0.014
0.034
0.009
0.009
0.010
0.012
Fig. 3. The
a-glucosidase inhibitory activity and estimation of IC50 of nitro
derivatives and reference compound acarbose. The results were presented as mean
values ꢁ SD from three separate experiments. *Statistically significance obtained at
P < 0.05 compared to the control (unpaired Students t-Test).