Structure-based 3D-Pharmacophore modeling to discover novel interleukin 6 inhibitors: An in silico screening, molecular dynamics simulations and binding free energy calculations

Interleukin 6 (IL-6) is a cytokine with various biological functions in immune regulation, hematopoiesis, and inflammation. Elevated IL-6 levels have been identified in several severe disorders such as sepsis, acute respiratory distress syndrome (ARDS), and most recently, COVID-19. The biological activity of IL-6 relies on interactions with its specific receptor, IL-6Rα, including the membrane-bound IL-6 receptor (mIL-6R) and the soluble IL-6 receptor (sIL-6R). Thus, inhibition of the interaction between these two proteins would be a potential treatment for IL-6 related diseases. To date, no orally available small-molecule drug has been approved. This study focuses on finding potential small molecules that can inhibit protein-protein interactions between IL-6 and its receptor IL-6Rα using its crystal structure (PDB ID: 5FUC). First, two pharmacophore models were constructed based on the interactions between key residues of IL-6 (Phe74, Phe78, Leu178, Arg179, Arg182) and IL-6Rα (Phe229, Tyr230, Glu277, Glu278, Phe279). A database of approximately 22 million compounds was screened using 3D-pharmacophore models, molecular docking models, and ADMET properties. By analyzing the interactive capability of successfully docked compounds with important amino acids, 12 potential ligands were selected for further analysis via molecular dynamics simulations. Based on the stability of the complexes, the high interactions rate of each ligand with the key residues of IL-6/IL-6Rα, and the low binding free energy calculation, two compounds ZINC83804241 and ZINC02997430, were identified as the most potential IL-6 inhibitor candidates. These results will pave the way for the design and optimization of more specific compounds to combat cytokine storm in severe coronavirus patients.