2
G. Maccari et al. / Bioorganic & Medicinal Chemistry Letters xxx (2017) xxx–xxx
sensus score using ComboScore as scoring function. The whole
NH
NH
complex was used to create a receptor object for the docking of
compound 1 with the Hybrid method.30 The best ten binding
modes obtained for the studied compound were further optimized
by energy minimization and calculation of the binding energy with
the Szybki and Zap toolkits respectively. The pose with the best
binding energy was saved for the consensus score. The final rank-
ing list was produced by applying a consensus scoring on the
results obtained with the shape comparison and the docking. Com-
plete results are reported in Table S1 in Supplementary data.
Among the top fifteen ranked structures the two most repre-
sented enzymes belonged to the families of the xylanase and the
chitinase. Xylanase is an enzyme involved in the degradation of
b-1,4-xylan. Chitinases are involved in the remodelling of the fun-
gal cell wall by degrading the chitin, its main component. Both
enzymes belong to the glycosidase superfamily, but the three crys-
tal structures of chitinase B from Aspergillus fumigatus were
resolved in complex with the natural compound Argifin (IC50 and
Ki respectively of 27 and 17 nM)31,32 (Fig. 2) and two oligopeptides
synthesized during a molecular simplification study,32,33 which
look similar to our compounds. We have thus hypothesized that
this enzyme could represent a possible target for our compound.
Argifin is the fermentation product of Gliocladium fungal cul-
ture.34–36 Chitin is the main component of the fungal cell wall; it
is a polymer of N-acetylglucosamine (GlcNAc) with a b-1,4 linkage
between monomers. Chitinases are the enzymes responsible for
chitin degradation and have been validated as a promising target
for the design of new therapeutic agents active against fungal
infections.37,38 Moreover, the inhibition of chitinase in insects is a
mechanism exploited in the production of new pesticides.39
To support the chitinase inhibition hypothesis, we performed
enzymatic assays against the homologue protein extracted from
Trichoderma viride. This protein was chosen for its commercial
availability and for the high degree of identity, being our attempts
to express the Candida chitinase unproductive. The enzyme was
incubated with different concentrations of inhibitor directly into
the reading plate at room temperature for three hours. At the
HN
NH
O
N
N
N
H
H
1
Fig. 1. Compound 1.
for our compounds. To confirm our hypothesis, we decided to opti-
mize the series of compounds in order to improve activity toward
the chitinase enzyme.
We setup the TF protocol by using a computational approach
consisting of the following six steps:
1. Selection of a set of ligand/protein complexes whose crystal
structures have been resolved;
2. Preparation of each complex by separating the protein from the
ligand and relaxing the system by means of energy
minimization;
3. Shape based comparison between the crystallized ligands and
our compound;
4. Docking of our compound in the binding site of the crystallo-
graphic structures;
5. Scoring of the candidate targets by means of a consensus score
rank-by-rank using both the results of the shape matching and
the docking;
6. Selection of the possible targets among the highest ranked
structures.
The set of complexes was taken from the PDB. In a first attempt,
we focused our attention on the crystal structures of protein from
the genus Candida, co-crystallized with at least one compound
(inhibitors or physiological substrates). This led to the selection
of only 95 crystal structures. To increase the chances of finding
an ideal target for our compounds, we enlarged the set of reference
structures by including all the protein of the Ascomycota phylum,
obtaining an initial database of 3636 structures. The second step,
i.e. preparation of the complexes, was performed using an in-house
python script that makes use of the OpenEye toolkits.17 The script
recognizes as ligand each molecule with at least four carbon atoms
and no more than 75 heavy atoms; compounds too far from the
protein atoms were ignored. Each ligand was processed by perceiv-
ing the bond orders and assigning the tautomerization and proto-
nation state at pH 7.4. The last two operations were performed
using the Quacpac libraries of the OpenEye toolkit.17,18 The protein
was prepared by adding the explicit hydrogens and calculating the
most probable protonation calling the external program
Reduce.19,20 Following the pre-treatment of both ligand and pro-
tein, the script recreates the complex and relaxes the system by
means of energy minimization, using the Szybki libraries included
in the OpenEye toolkit.17,21 In case of complexes with more than
one ligand two different files were saved. To improve the efficiency
of the computation the whole script was parallelized using the
mpi4py Python libraries22–24 interfaced with the openMPI
libraries.25 As a result, we obtained 2177 optimized complexes.
Steps from 3 to 5 were implemented in another python script
that makes use of the OEDocking, Shape, Quacpac, Szybki and
Zap libraries of the OpenEye toolkit.17,18,21,26–28 The script accepts
as input a list of pre-processed complexes and the studied com-
pound, assuming that the reference is a multi-conformer file con-
sisting of conformations calculated using Omega2.29 For each
complex the bound ligand was used as a reference for a shape
and colour atoms comparison against every conformer of com-
pound 1. The best values were stored to be used in the final con-
end of the incubation, the substrate (4-nitrophenyl N-acetyl-b-D-
glucosaminide) was added and the reading was conducted by
recording the values up to a maximum time of 10 min. The sub-
strate concentration was measured by reading UV absorption at
300 nm (
sis was measured at 400 nm (
result, both IC50 and Ki could be calculated.
e
M300: 11,100 MÀ1ÁcmÀ1), while the product of hydroly-
D
e
M400: 2,120 MÀ1ÁcmÀ1). As a
The assay has been performed on compound 1 and on other
compounds of the same family, either previously described14 or
prepared with the same procedure (Scheme S1 Supplementary
data). The results are reported in Table 1. The compounds were
active against chitinase with IC50 values of the same order of mag-
nitude of the MIC measured in cells.
O
O
O
N
HN
O
OH
HN
NH
O
NH
O
HN
O
NH
NH
NH
OH
O
Argifin
Fig. 2. Structures of one of the three inhibitors co-crystallized with the complex
1W9V.