Molecular dynamics simulations provide insights into ULK-101 potency and selectivity toward autophagic kinases ULK1/2

Kinase domains are highly conserved within protein kinases in both sequence and structure. Many factors, including phosphorylation, amino acid substitutions or mutations, and small molecule inhibitor binding, influence conformations of the kinase domain and enzymatic activity. The serine/threonine kinases ULK1 and ULK2 are highly conserved with N- and C-terminal domains, phosphate-binding P-loops, -αC-helix, regulatory and catalytic spines, and activation loop DFG and APE motifs. Here, we performed molecular dynamics (MD) simulations to understand better the potency and selectivity of the ULK1/2 small molecule inhibitor, ULK-101. We observed stable bound states for ULK-101 to the adenosine triphosphate (ATP)-binding site of ULK2, coordinated by hydrogen bonding with the hinge backbone and the catalytic lysine sidechain. Notably, ULK-101 occupies a hydrophobic pocket associated with the N-terminus of the -αC-helix. Large movements in the P-loop are also associated with ULK-101 inhibitor binding and exit from ULK2. Our data further suggests that ULK-101 could induce a folded P-loop conformation and hydrophobic pocket reflected in its nanomolar potency and kinome selectivity.