Identification of potential CRAC channel inhibitors: Pharmacophore mapping, 3D-QSAR modelling,and molecular docking approach

Upregulation of store-operated Ca²⁺ influx via ORAI1, an integral component of the CRAC channel, is responsible for abnormal cytokine release in active rheumatoid arthritis, and therefore ORAI1 has been proposed as an attractive molecular target. In this study, we attempted to predict the mechanical insights of ORAI1 inhibitors through pharmacophore modelling, 3D-QSAR, molecular docking and free energy analysis. Various hypotheses of pharmacophores were generated and from that, a pharmacophore hypothesis with two hydrogen bond acceptors, one hydrogen bond donor and two aromatic rings (AADRR) resulted in a statistically significant 3D-QSAR model (r² = 0.84 and q² = 0.74). We believe that the obtained statistical model is a reliable QSAR model for the diverse dataset of inhibitors against the IL-2 production assay. The visualization of contours in active and inactive compounds generated from the 3D-QSAR models and molecular docking studies revealed major interaction with GLN108, HIS113 and ASP114, and interestingly, these residues are located near the Ca²⁺ selectivity filter region. Free energy binding analysis revealed that Coulomb energy, van der Waals energy and non-polar solvation terms are more favourable for ligand binding. Thus, the present study provides the physical and chemical requirements for the development of novel ORAI1 inhibitors with improved biological activity.