Abstract
Kinases are important targets for drug development. However, accounting for the impact of possible structural rearrangements on the binding of kinase inhibitors is complicated by the extensive flexibility of their catalytic domain. The dynamic N-lobe contains four particular mobile structural elements: the Asp-Phe-Gly (DFG) motif, the phosphate (P) positioning loop, the activation (A) loop, and the αC helix. In our previous study [Meng et al. J. Chem. Theory Comput. 2018 14, 2721-2732], we combined various simulation techniques with Markov state modeling (MSM) to explore the free energy landscape of Abl kinase beyond conformations that are known from X-ray crystallography. Here we examine the resulting Markov model in greater detail by analyzing its metastable states. A characterization of the states in terms of their DFG state, P-loop, and αC conformations is presented and compared to existing classification schemes. Several metastable states are found to be structurally close to known crystal structures of different kinases in complex with a variety of inhibitors. These results suggest that the set of conformations accessible to tyrosine kinases may be shared within the entire family and that the conformational dynamics of one kinase in the absence of any ligand can provide meaningful information about possible target conformations for inhibitors of any member of the kinase family.