Dataset Information

An Investigation of Molecular Docking and Molecular Dynamic Simulation on Imidazopyridines as B-Raf Kinase Inhibitors.

ABSTRACT: In the recent cancer treatment, B-Raf kinase is one of key targets. Nowadays, a group of imidazopyridines as B-Raf kinase inhibitors have been reported. In order to investigate the interaction between this group of inhibitors and B-Raf kinase, molecular docking, molecular dynamic (MD) simulation and binding free energy (?Gbind) calculation were performed in this work. Molecular docking was carried out to identify the key residues in the binding site, and MD simulations were performed to determine the detail binding mode. The results obtained from MD simulation reveal that the binding site is stable during the MD simulations, and some hydrogen bonds (H-bonds) in MD simulations are different from H-bonds in the docking mode. Based on the obtained MD trajectories, ?Gbind was computed by using Molecular Mechanics Generalized Born Surface Area (MM-GBSA), and the obtained energies are consistent with the activities. An energetic analysis reveals that both electrostatic and van der Waals contributions are important to ?Gbind, and the unfavorable polar solvation contribution results in the instability of the inhibitor with the lowest activity. These results are expected to understand the binding between B-Raf and imidazopyridines and provide some useful information to design potential B-Raf inhibitors.

SUBMITTER: Xie H

PROVIDER: S-EPMC4661884 | biostudies-literature | 2015

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

An Investigation of Molecular Docking and Molecular Dynamic Simulation on Imidazopyridines as B-Raf Kinase Inhibitors.

Xie Huiding H Li Yupeng Y Yu Fang F Xie Xiaoguang X Qiu Kaixiong K Fu Jijun J

International journal of molecular sciences 20151116 11

In the recent cancer treatment, B-Raf kinase is one of key targets. Nowadays, a group of imidazopyridines as B-Raf kinase inhibitors have been reported. In order to investigate the interaction between this group of inhibitors and B-Raf kinase, molecular docking, molecular dynamic (MD) simulation and binding free energy (ΔGbind) calculation were performed in this work. Molecular docking was carried out to identify the key residues in the binding site, and MD simulations were performed to determin ...[more]

PMID: 26580609

Similar Datasets

Project description:Protein-protein interactions are indispensable physiological processes regulating several biological functions. Despite the availability of structural information on protein-protein complexes, deciphering their complex topology remains an outstanding challenge. Raf kinase inhibitory protein (RKIP) has gained substantial attention as a favorable molecular target for numerous pathologies including cancer and Alzheimer's disease. RKIP interferes with the RAF/MEK/ERK signaling cascade by endogenously binding with C-Raf (Raf-1 kinase) and preventing its activation. In the current investigation, the binding of RKIP with C-Raf was explored by knowledge-based protein-protein docking web-servers including HADDOCK and ZDOCK and a consensus binding mode of C-Raf/RKIP structural complex was obtained. Molecular dynamics (MD) simulations were further performed in an explicit solvent to sample the conformations for when RKIP binds to C-Raf. Some of the conserved interface residues were mutated to alanine, phenylalanine and leucine and the impact of mutations was estimated by additional MD simulations and MM/PBSA analysis for the wild-type (WT) and constructed mutant complexes. Substantial decrease in binding free energy was observed for the mutant complexes as compared to the binding free energy of WT C-Raf/RKIP structural complex. Furthermore, a considerable increase in average backbone root mean square deviation and fluctuation was perceived for the mutant complexes. Moreover, per-residue energy contribution analysis of the equilibrated simulation trajectory by HawkDock and ANCHOR web-servers was conducted to characterize the key residues for the complex formation. One residue each from C-Raf (Arg398) and RKIP (Lys80) were identified as the druggable "hot spots" constituting the core of the binding interface and corroborated by additional long-time scale (300 ns) MD simulation of Arg398Ala mutant complex. A notable conformational change in Arg398Ala mutant occurred near the mutation site as compared to the equilibrated C-Raf/RKIP native state conformation and an essential hydrogen bonding interaction was lost. The thirteen binding sites assimilated from the overall analysis were mapped onto the complex as surface and divided into active and allosteric binding sites, depending on their location at the interface. The acquired information on the predicted 3D structural complex and the detected sites aid as promising targets in designing novel inhibitors to block the C-Raf/RKIP interaction.

Project description:The aurora kinase is a key enzyme that is implicated in tumor growth. Research revealed that small molecules that target aurora kinase have beneficial effects as anticancer agents. In the present study, in order to identify potential antibreast cancer agents with aurora kinase inhibitory activity, we employed QSARINS software to perform the quantitative structure-activity relationship (QSAR). The statistical values resulted from the study include R2 = 0.8902, CCCtr = 0.7580, Q2 LOO = 0.7875, Q2LMO = 0.7624, CCCcv = 0.7535, R2ext = 0.8735, and CCCext = 0.8783. Among the four generated models, the two best models encompass five important variables, including PSA, EstateVSA5, MoRSEP3, MATSp5, and RDFC24. The parameters including the atomic volume, atomic charges, and Sanderson's electronegativity played an important role in designing newer lead compounds. Based on the above data, we have designed six series of compounds including 1a-e, 2a-e, 3a-e, 4a-e, 5a-e, and 6a-e. All these compounds were subjected to molecular docking studies by using AutoDock v4.2.6 against the aurora kinase protein (1MQ4). Among the above 30 compounds, the 2-amino thiazole derivatives 1a, 2a, 3e, 4d, 5d, and 6d have excellent binding interactions with the active site of 1MQ4. Compound 1a had the highest docking score (-9.67) and hence was additionally subjected to molecular dynamic simulation investigations for 100 ns. The stable binding of compound 1a with 1MQ4 was verified by RMSD, RMSF, RoG, H-bond, molecular mechanics-generalized Born surface area (MM-GBSA), free binding energy calculations, and solvent-accessible surface area (SASA) analyses. Furthermore, newly designed compound 1a exhibited excellent ADMET properties. Based on the above findings, we propose that the designed compound 1a may be utilized as the best theoretical lead for future experimental research of selective inhibition of aurora kinase, therefore assisting in the creation of new antibreast cancer drugs.

Dataset Information

An Investigation of Molecular Docking and Molecular Dynamic Simulation on Imidazopyridines as B-Raf Kinase Inhibitors.

Publications

An Investigation of Molecular Docking and Molecular Dynamic Simulation on Imidazopyridines as B-Raf Kinase Inhibitors.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets