Identification of Potent Small-Molecule PCSK9 Inhibitors Based on Quantitative Structure-Activity Relationship, Pharmacophore Modeling, and Molecular Docking Procedure

The leading cause of atherosclerotic cardiovascular disease (ASCVD) is elevated low-density lipoprotein cholesterol (LDL-C). Proprotein convertase subtilisin/kexin type 9 (PCSK9) attaches to the domain of LDL receptor (LDLR), diminishing LDL-C influx and LDLR cell surface presentation in hepatocytes, resulting in higher circulating LDL-C levels. PCSK9 dysfunction has been linked to lower levels of plasma LDLC and a decreased risk of coronary heart disease (CHD). Herein, using virtual screening tools, we aimed to identify a potent small-molecule PCSK9 inhibitor in compounds that are currently being studied in clinical trials. We first performed chemical absorption, distribution, metabolism, excretion, and toxicity (ADMET) filtering of 9800 clinical trial compounds obtained from the ZINC 15 database using Lipinski's rule of 5 and achieved 3853 compounds. Two-dimensional (2D) quantitative structure-activity relationship (QSAR) was initiated by computing molecular descriptors and selecting important descriptors of 23 PCSK9 inhibitors. Multivariate calibration was performed with the partial least square regression (PLS) method with 18 compounds for training to design the QSAR model and 5 compounds for the test set to assess the model. The best latent variables (LV) (LV=6) with the lowest value of Root-Mean-Square Error of Cross-Validation (RMSECV) of 0.48 and leave-one-out cross-validation correlation coefficient (R2CV) = 0.83 were obtained for the QSAR model. The low RMSEC (0.21) with high R²cal (0.966) indicates the probability of fit between the experimental data and the calibration model. Using QSAR analysis of 3853 compounds, 2635 had a pIC50<1 and were considered for pharmacophore screening. The PHASE module (a complete package for pharmacophore modeling) designed the pharmacophore hypothesis through multiple ligands. The top 14 compounds (pIC50>1) were defined as active, whereas 9 (pIC50<1) were considered as an inactive set. Three five-point pharmacophore hypotheses achieved the highest score: DHHRR1, DHHRR2, and DHRRR1. The highest and best model with survival scores (5.365) was DHHRR1, comprising 1 hydrogen donor (D), 2 hydrophobic groups (H), and 2 rings of aromatic (R) features. We selected the molecules with a higher 1.5 fitness score (257 compounds) in pharmacophore screening (DHHRR1) for molecular docking screening. Molecular docking indicates that ZINC000051951669, with a binding affinity: of -13.2 kcal/mol and 2 H-bonds, has the highest binding to the PCSK9 protein. ZINC000011726230 with energy binding: -11.4 kcal/mol and 3 H-bonds, ZINC000068248147 with binding affinity: -10.7 kcal/mol and 1 H-bond, ZINC000029134440 with a binding affinity: -10.6 kcal/mol and 4 H-bonds were ranked next, respectively. To conclude, the archived molecules identified as inhibitory PCSK9 candidates, and especially ZINC000051951669 may therefore significantly inhibit PCSK9 and should be considered in the newly designed trials.