390
J. Xing et al. / European Journal of Medicinal Chemistry 95 (2015) 388e399
fragment library was comprised with 305 fragments based on the
39 known ligands under the established setting rules. The infor-
mation about the spatial distribution in the target was included in
the resultant fragment library, as well as about their chemical
structures. The binding conformations of inhibitor rivaroxaban
(PDB code 2w26) and its fragments were shown in Fig. 4, and it was
found that the inhibitor was deconstructed into three component
fragments which maintained the initial orientations and spatial
positions of the ligand.
docked poses in each program. During screening, the recom-
mended criterion for similar docking results was slightly modified
from RMSD < 2.0 Å to RMSD < 3.0 Å because of limited number of
hit compounds when the proposed criteria applied. The 380 top-
ranking hit compounds from last step were independently
docked to the fXa structure (PDB code 2w26) by LigandFit and
Glide. And the top 100 compounds resulting from each docking
program were combined, and their poses were compared visually
followed by the calculation of corresponding RMSD value (RMSD_1
in Fig. 7). As a result, 43 compounds containing similar docked
poses generated from two different docking programs were
considered to be valid docking results.
2.2. Evaluation of fragment-based compound library and hits
identification
2.2.1. Construction of new molecular library
2.2.5. RMSD-based scoring function
According to the setting rules, All the new molecules should be
constructed by the fragments from different initial sets (Fig. 5). As a
result, 11858 compounds were developed from 305 fragments by
chemical space analysis.
To further evaluate the quality of the 43 compounds, the RMSD
values between the combining conformations and the predicted
binding conformations which were generated from molecular
docking were then calculated to re-rank hit compounds. The RMSD
values of the 43 compounds were showed in Fig. 7 (RMSD_2), and it
was found that most of them were distributed around 2.00 Å, with
five compounds <1.00 Å and three >3.00 Å. Table 1 exhibits all the
detailed RMSD values of the 43 potential fXa inhibitors in
ascending order, and Table 2 lists the structures of top-ten new
compounds generated by this stage. Among them, three known
2.2.2. Lipinski's rule and ADMET profile
To effectively abandon the redundant compounds and generate
a library with high quality, some virtual screening, including Lip-
inski's rule and ADMET profile, was firstly performed to filter the
new molecules. Finally, 814 compounds were selected from the
library.
compounds (1876, RMSD
¼
0.1487, Ki
¼
0.14 nM; 2398,
RMSD ¼ 0.3021, IC50 ¼ 27 nM; 643, RMSD ¼ 1.1324, Ki ¼ 90 nM)
with good fXa inhibitory activity were also constructed with the
fragments from the different initial inhibitors. From the compari-
son of RMSD values and activity data of the three compounds, we
can find that the smaller deviations between their binding and
combining conformations that they obtain, the higher inhibitory
activity they will display. This correlation is consistent with our
assumption. The RMSD values of top-ten new compounds were all
small and in the range from 0.1487 to 1.3938, indicating that all
fragments maintained their original conformation and binding
mode in the newly constructed molecules which would show good
activity against fXa with a high probability. As for the selection of
hits compounds, synthetic tractability and cost were firstly
considered. As a result, compounds 3780 and 319 were chosen as
hit compounds.
2.2.3. Structural based pharmacophore generation
Twelve fXa-inhibitor cocrystal structures which were included in
the above mentioned thirty-nine representative fXa inhibitors were
superimposed with the same coordinate. The Ludi interaction map
was used to generate pharmacophore features which consisted of
hydrogen bond acceptor (HBA), hydrogen bond donor (HBD), aro-
matic ring and hydrophobic features (Fig. 6A and B). Then some
common pharmacophoric constraints in the active site of fXa were
retained and some unrepresentative features were deleted. As a
result, six sets were considered as the key pharmacophoric features
(Fig. 6C). Subsequently, the central pharmacophore feature of each
set was defined as the common feature and constructed the Hypo1
(Fig. 6D). Hypo1 was then aligned to rivaroxaban. As were shown in
Fig. 6E and F, one HBA and HBD, one aromatic ring in the S4 pocket,
two hydrophobic cores in S1 pocket, one HBA and one HBD were all
included in the SBP model. Eventually, Hypo 1 was employed as the
query to quickly identify compounds in the library generated from
last stage, which have an excellent mapping. As a result, 380 com-
pounds were well matched this model.
The interaction mode of two hit compounds with fXa predicted
by glide in Schrodinger 2009 were illustrated in Fig. 8. For the
interaction mechanism of 3780 with fXa (Fig. 8a), 5-
chlorothiophene occupies hydrophobic S1pocket formed by
Ala190, Val213 and Tyr228 and forms a Cl-
p interaction between
the chloro atom and the phenyl ring of Tyr228, which played a part
in the interaction between inhibitors and fXa. The NH of scaffold
carboxamide forms a hydrogen bond with carbonyl oxygen of
Gly219. The P4 phenyllactam of 3780 is flanked by the phenyl
groups of Phe174 and Tyr99 in the S4 aryl binding pocket and forms
an edge to face interaction with Trp215. The carbonyl oxygen
(scaffold carboxamide) also forms hydrogen bond with Gly216. For
Compound 319 (Fig. 8B), chloropyridine occupies the S1 pocket
more deeply than 3780 because the carboxamide group extends P1
2.2.4. Consensus docking
The compounds resulting pharmacophore screening were then
analyzed by docking study. Presently, there are many potential
pitfalls and inherent limitations remaining in docking program.
There is no single docking program that is competent in all cases,
which poses a pressing question regarding the optimal protocol to
select for docking-based virtual screening. Recent studies have re-
ported a novel consensus docking strategy, which instead of
combining various docking scores, comprehensively integrated the
results of docking poses from separate docking programs [35e37].
The use of consensus docking poses could remarkably improve the
reliability of docking and remedy the drawback existing in single
docking program, which always provide similar docking scores
with no way of distinguishing which one is correct [36]. Previously,
we have successfully adopted this method in identification of
dipeptidyl peptidase IV inhibitors [37]. Now, consensus docking
strategy is also employed by using the following two different
docking programs: LigandFit in DS2.5 and glide in Schrodinger
2009. DockScore and GScore were selected to separately rank
by 2-bonds length, but Cl-p interaction with Tyr228 was not
formed because the chlorine atom deviated from the phenyl ring of
Tyr228. The other interaction mode is similar to 3780.
To confirm the reliability of the strategy, compounds 3780 and
319 were firstly synthesized and evaluated for in vitro inhibition of
fXa. To our delight, the two compounds exhibited inhibitory activity
against fXa with the IC50 value of 67 nM and 298 nM, respectively,
which encouraged us to further explore the efficacy of analogs of
these two hit compounds. As a result, a series of analogs of com-
pounds 3780 and 319 were designed, synthesized and evaluated for
their fXa inhibitory activity.