R. Vieira Veloso et al.
Bioorganic Chemistry 109 (2021) 104709
3.6. Molecular docking studies
Table 2
Protein targets and parameters of binding after molecular docking with the 11m
molecule.
Protein
After identifying the activity and pharmacokinetics proprieties, 11m
was chosen for a virtual screening by molecular docking for predicting
the interactions between a target protein and the molecule and elucidate
its possible mechanism of action. Six protein targets were selected based
on their participation on the antioxidant activity and simulated with
11m (Fig. 6). The parameters resulting from the protein target and
ligand interaction are shown in Table 2.
PDB
Binding
energy
(kcal/
mol)
Full
Binding
interaction
(ΔG)
Bond
length
(Å)
code
Fitness
(kcal/
mol)
Superoxide
Dismutase
(SOD)
1CB4
(chain
A + B)
1DGB
(chain
A)
ꢀ 9.14
ꢀ 1369.18
LigOH-
3.62
3.70
Asp11
LigN2-Lys9
LigNH-
Although the lowest binding energy has been reached in the inter-
action with superoxide dismutase (ꢀ 9.14 kcal/mol), the ligand was
positioned out of the catalytic site, indicating that the ligand does not
interfere with the enzyme antioxidant activity. The glutathione peroxi-
dase 3, NADPH oxidase, catalase and lipoxygenase active site were not
covered by 11m as well, indicating that these proteins are not the target,
even having low binding energy and/or full fitness.
Catalase
ꢀ 8.56
ꢀ 2144.19
3.15
2.79
3.36
Arg431
LigOH-
Asp157
LigOH-
Lys349
Glutathione
peroxidase 3
2R37
(chain
A)
ꢀ 8.90
ꢀ 8.60
ꢀ 1049.03
ꢀ 2349.03
LigNH-
2.42
3.93
3.50
Thr165
LigO2-
On the other hand, 11m interacted with the peroxiredoxin-5 cata-
lytic pocket, with a strong hydrogen bond (ꢀ 8.59 kcal/mol) identified in
the side chain of the key residues Arg127, and Gly46 with low bond
length (4.29 and 2.20 Å, respectively) [4]. Thus, we have an indication
that peroxiredoxin-5 could be the molecular target for 11m.
Peroxiredoxins are thiol-dependent peroxidases, able to reduce
hydrogen peroxide, alkyl hydroperoxides, and peroxynitrite [19]. They
are well conserved enzymes, mammalian, found in both cytosol and
mitochondria, already identified in the liver and lung. These enzymes
were related to several pathologies involving tendons, cartilages, nerves,
carcinomas (ovarian, breast, colon, adrenocortical), thyroid, astrocytes,
retina, and also blood cells [25].
Arg168(NH)
LigO2-
Leu162(NH)
LigO2-
NADPH
oxidase
2CDU
(chain
B)
3.67
5.44
3.13
Arg308(NH)
LigO-Arg308
(NH2)
LigN1-
Arg305
(NH2)
Lipoxygenase
1N8Q
(chain
A)
ꢀ 6.56
ꢀ 8.59
ꢀ 3881.73
ꢀ 783.09
LigO-Arg621
(NH)
3.19
3.38
4.46
LigO-
Arg621-N
LigN2-
The Peroxiredoxin-5 (PRDX5) expression was related to acute in-
flammatory effects, as it was upregulated in rats after the injection of
lipopolysaccharides and during stress of endothelial cells and macro-
phages [24]. The increase in activity was also found in neutrophils and
monocytes [27].
Glu625(O)
LigN-Gly46
(HN)
Peroxiredoxin-
5
3MNG
(chain
A)
2.20
4.29
Arg127
(NH2)
The molecular docking was performed using the PRDX5 because it is
the most available crystalized protein isoform, the opposite of PRDX2,
related to SCD. However, the catalytic site is conserved, which indicates
that 11m could interact with the PRDX2 as well.
the intra or extracellular environment [17]. According to Lipinski’s
rules, both values are in accordance with the desirable one for a mole-
cule being a drug. These binding predictions also help the verification of
the interactive capacity of a molecule with the enzymes for the drug
metabolism (CYP P450 superfamily). The 11m does not have potential
to inhibit any of the enzymes evaluated (CYP1A2, CYP2C19, CYP2C9,
CYP2D6, CYP3A4), which indicates low hepatotoxicity and drug inter-
action. On the other hand, 11g has potential to inhibit 2 CYPs, indicating
possible drug interaction and effects on xenobiotics metabolism.
The acute oral toxicity of both molecules (11g and 11m) was esti-
mated as category III (Table 1). This parameter is important because
indicates the possible adverse effects occurring within a short time after
an oral administration of a substance and the category III is referred as
“slightly toxic and slightly irritating” according to US EPA for chemicals
[53].
In the SCD, there is a huge quantity of ROS being produced by the
RBC, such as H2O2, and consequently the activities of antioxidant en-
zymes are increased, such as superoxide dismutase, catalase, glutathione
peroxidase, and three isoforms of peroxiredoxin (I, II, and VI). So, it is
expected that the PRDX5 activity would be increased, compatible to this
́
radicals production, trying to eliminate intracellular concentration of
H2O2 [44].
Indeed, we have seen that the 11m molecule reduced the hemolysis
caused by H2O2, and besides the scavenging mechanism, it can be acting
in the enzymatic activity of PRDX5 trying to remove the ROS.
Peroxiredoxins can be inactivated during catalysis because the
cysteine residue, important for the active site, is oxidated to sulfinic acid
(Cys-SO2H), being their reactivation achieved by reduction of the sul-
finic moiety [55]. Thus, by acting on the PRDX, the 11m molecule could
act coordinating the active site of the enzyme, increasing its catalytic
efficiency, consequently increasing the H2O2 removal, as demonstrated
here.
In the classification stablished by OECD (Organization for Economic
Co-operation and Development) [40], four categories are considered
regarding lethal dose 50% (LD50) values: Category I contains com-
pounds with LD50 values less than or equal to 50 mg/kg; Category II
compounds with LD50 values >50 through 500 mg/kg; Category III
compounds with LD50 values >500 through 5000 mg/kg; Category IV
compounds with LD50 values greater than 5000 mg/kg.
4. Conclusion
The LD50 values for 11g and 11m, 500–5000 mg/kg predicted, are
comparable to hydroxyurea, which was experimentally calculated as
5760 mg/kg after an oral administration, in rats. Considering that a
chosen drug for sickle cell disease, our molecules have potential for a
drug development program [37].
After evaluating 30 molecules, we could select one (11m) according
to its antioxidant activity, evaluated by two different methods, which
was able to reduce the hemolysis caused by reactive oxygen species in
red blood cell, a model that mimics the sickle cell disease. The radical
scavenging activity of 11m can be added to the possible activation of
peroxiredoxin-5, an antioxidant enzyme, which could remove H2O2 in
both preventive and treatment action. Moreover, this molecule has
adequate pharmacokinetics proprieties, hence making it a good candi-
date as a drug to treat sickle cell disease.
This feature, added to the one hydrogen-bond acceptors surplus and
the molecular weight above 500 g/mol, made us eliminate 11g for the
study, being the 11m the chosen one for the molecular docking
evaluation.
9