Role of ELA region in auto-activation of mutant KIT receptor: a molecular dynamics simulation insight

KIT receptor is the prime target in gastrointestinal stromal tumor (GISTs) therapy. Second generation inhibitor, Sunitinib, binds to an inactivated conformation of KIT receptor and stabilizes it in order to prevent tumor formation. Here, we investigated the dynamic behavior of wild type and mutant D816H KIT receptor, and emphasized the extended A-loop (EAL) region (805-850) by conducting molecular dynamics simulation (similar to 100 ns). We analyzed different properties such as root mean square cutoff or deviation, root mean square fluctuation, radius of gyration, solvent-accessible surface area, hydrogen bonding network analysis, and essential dynamics. Apart from this, clustering and cross-correlation matrix approach was used to explore the conformational space of the wild type and mutant EAL region of KIT receptor. Molecular dynamics analysis indicated that mutation (D816H) was able to alter intramolecular hydrogen bonding pattern and affected the structural flexibility of EAL region. Moreover, flexible secondary elements, specially, coil and turns were dominated in EAL region of mutant KIT receptor during simulation. This phenomenon increased the movement of EAL region which in turn helped in shifting the equilibrium towards the active kinase conformation. Our atomic investigation of mutant KIT receptor which emphasized on EAL region provided a better insight into the understanding of Sunitinib resistance mechanism of KIT receptor and would help to discover new therapeutics for KIT-based resistant tumor cells in GIST therapy.