3D-QSAR-Based Pharmacophore Modeling, Virtual Screening, and Molecular Dynamics Simulations for the Identification of Spleen Tyrosine Kinase Inhibitors

Abstract

Spleen tyrosine kinase (SYK) is an essential mediator of immune cell signaling and has been anticipated as a therapeutic target for autoimmune diseases, notably rheumatoid arthritis, allergic rhinitis, asthma, and cancers. Significant attempts have been undertaken in recent years to develop SYK inhibitors; however, limited success has been achieved due to poor pharmacokinetics and adverse effects of inhibitors. The primary goal of this research was to identify potential inhibitors having high affinity, selectivity based on key molecular interactions, and good drug-like properties than the available inhibitor, fostamatinib. In this study, a 3D-QSAR model was built for SYK based on known inhibitor IC₅₀ values. The best pharmacophore model was then used as a 3D query to screen a drug-like database to retrieve hits with novel chemical scaffolds. The obtained compounds were subjected to binding affinity prediction using the molecular docking approach, and the results were subsequently validated using molecular dynamics (MD) simulations. The simulated compounds were ranked according to binding free energy (ΔG), and the binding affinity was compared with fostamatinib. The binding mode analysis of selected compounds revealed that the hit compounds form hydrogen bond interactions with hinge region residue Ala451, glycine-rich loop residue Lys375, Ser379, and DFG motif Asp512. Identified hits were also observed to form a desirable interaction with Pro455 and Asn457, the rare feature observed in SYK inhibitors. Therefore, we argue that identified hit compounds ZINC98363745, ZINC98365358, ZINC98364133, and ZINC08789982 may help in drug design against SYK.

Keywords: 3D QSAR; MD simulation; SYK inhibitor; autoimmune diseases; molecular docking; pharmacophore.