Theoretical Characterization of Binding Mode of Organosilicon Inhibitor with p38: Docking, MD Simulation and MM/GBSA Free Energy Approach

P38 mitogen activated protein (MAP) kinase is an important anti-inflammatory drug target, which can be activated by responding to various stimuli such as stress and immune response. Based on the conformation of the conserved DFG loop (in or out), binding inhibitors are termed as type-I and II. Type-I inhibitors are ATP competitive, whereas type-II inhibitors bind in DFG-out conformation of allosteric pocket. It remains unclear that how these allosteric inhibitors stabilize the DFG-out conformation and interact. Organosilicon compounds provide unusual opportunity to enhance potency and diversity of drug molecules due to their low toxicity. However, very few examples have been reported to utilize this property. In this regard, we performed docking of an inhibitor (BIRB) and its silicon analog (Si-BIRB) in an allosteric binding pocket of p38. Further, molecular dynamics (MD) simulations were performed to study the dynamic behavior of the simulated complexes. The difference in the biological activity and mechanism of action of the simulated inhibitors could be explained based on the molecular mechanics/generalized Born surface area (MM/GBSA) binding free energy per residue decomposition. MM/GBSA showed that biological activities were related with calculated binding free energy of inhibitors. Analyses of the per-residue decomposed energy indicated that van der Waals and non-polar interactions were predominant in the ligand-protein interactions. Further, crucial residues identified for hydrogen bond, salt bridge and hydrophobic interactions were Tyr35, Lys53, Glu71, Leu74, Leu75, Ile84, Met109, Leu167, Asp168 and Phe169. Our results indicate that stronger hydrophobic interaction of Si-BIRB with the binding site residues could be responsible for its greater binding affinity compared with BIRB.