P2Y12 receptor is an attractive target for the anti-platelet therapies, treating various thrombotic diseases. In this work, a total of 107 6-aminonicotinate-based compounds as potent P2Y12 antagonists were studies by a molecular modeling study combining three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations to explore the decisive binding conformations of these antagonists with P2Y12 and the structural features for the activity. The optimum CoMFA and CoMSIA models identified satisfactory robustness and good predictive ability, with R2 = .983, q2 = .805, [Formula: see text] = .881 for CoMFA model, and R2 = .935, q2 = .762, [Formula: see text] = .690 for CoMSIA model, respectively. The probable binding modes of compounds and key amino acid residues were revealed by molecular docking. MD simulations and MM/GBSA free energy calculations were further performed to validate the rationality of docking results and to compare the binding modes of several compound pairs with different activities, and the key residues (Val102, Tyr105, Tyr109, His187, Val190, Asn191, Phe252, His253, Arg256, Tyr259, Thr260, Val279, and Lys280) for the higher activity were pointed out. The binding energy decomposition indicated that the hydrophobic and hydrogen bond interactions play important roles for the binding of compounds to P2Y12. We hope these results could be helpful in design of potent and selective P2Y12 antagonists.
Keywords: 3D-QSAR; MD simulations; P2Y12 antagonists; anti-platelet; docking.