Leucine-rich repeat kinase 2 (LRRK2)
remains a viable target for
drug development since the discovery of the association of its mutations
with Parkinson’s disease (PD). G2019S (in the kinase domain)
is the most common mutation for LRRK2-based PD. Though various types
of inhibitors have been developed for the kinase domain to reduce
the effect of the mutation, understanding the working of these inhibitors
at the molecular level is still ongoing. This study focused on the
exploration of the dissociation mechanism (pathways) of inhibitors
from (WT and G2019S) LRRK2 kinase (using homology model CHK1 kinase),
which is one of the crucial aspects in drug discovery. Here, two ATP-competitive
type I inhibitors, PF-06447475 and MLi-2 (Comp1 and Comp2 ), and
one non-ATP-competitive type II inhibitor, rebastinib (Comp3), were
considered for this investigation. To study the unbinding process,
random accelerated molecular dynamics simulations were performed.
The binding free energies of the three inhibitors for different egression
paths were determined using umbrella sampling. This work found four
major egression pathways that were adopted by the inhibitors Comp1
(path1, path2, and path3), Comp2 (path1, path2 and path3), and Comp3
(path3 and path4). Also, the mechanism of unbinding for each path
and key residues involved in unbinding were explored. Mutation was
not observed to impact the preference of the particular egression
pathways for both LRRK2-Comp1 and -Comp2 systems. However, the findings
suggested that the size of the inhibitor molecules might have an effect
on the preference of the egression pathways. The binding energy and
residence time of the inhibitors followed a similar trend to experimental
observations. The findings of this work might provide insight into
designing more potent inhibitors for the G2019S LRRK2 kinase.