Factor IXa (FIXa), a blood coagulation factor, is specifically inhibited at the initiation stage of the coagulation cascade, promising an excellent approach for developing selective and safe anticoagulants. Eighty-four amidinobenzothiophene antithrombotic derivatives targeting FIXa were selected to establish three-dimensional quantitative structure-activity relationship (3D-QSAR) and three-dimensional quantitative structure-selectivity relationship (3D-QSSR) models using comparative molecular field analysis and comparative similarity indices analysis methods. Internal and external cross-validation techniques were investigated as well as region focusing and bootstrapping. The satisfactory q (2) values of 0.753 and 0.770, and r (2) values of 0.940 and 0.965 for 3D-QSAR and 3D-QSSR, respectively, indicated that the models are available to predict both the inhibitory activity and selectivity on FIXa against Factor Xa, the activated status of Factor X. This work revealed that the steric, hydrophobic, and H-bond factors should appropriately be taken into account in future rational design, especially the modifications at the 2'-position of the benzene and the 6-position of the benzothiophene in the R group, providing helpful clues to design more active and selective FIXa inhibitors for the treatment of thrombosis. On the basis of the three-dimensional quantitative structure-property relationships, 16 new potent molecules have been designed and are predicted to be more active and selective than Compound 33, which has the best activity as reported in the literature.
Keywords: 3D-QSAR; 3D-QSSR; CoMFA; CoMSIA; benzothiophene antithrombosis.