Identification of novel LRRK2 inhibitors by structure-based virtual screening and alchemical free energy calculation.

The Leucine-rich repeat kinase 2 (LRRK2) target has been identified as a promising drug target for Parkinson's disease (PD) treatment. This study focuses on optimizing the activity of LRRK2 inhibitors using alchemical relative binding free energy (RBFE) calculations. Initially, we assessed various free energy calculation methods across different LRRK2 kinase inhibitor scaffolds. The results indicate that alchemical free energy calculations are promising for prospective predictions on LRRK2 inhibitors, especially for the aminopyrimidine scaffold with an RMSE of 1.15 kcal mol-1 and Rp of 0.83. Following this, we optimized a potent LRRK2 kinase inhibitor identified from previous virtual screenings, featuring a novel scaffold. Guided by RBFE predictions using alchemical methods, this optimization led to the discovery of compound LY2023-001. This compound, with a [1,2,4]triazolo[5,6-b]indole scaffold, exhibited enhanced inhibitory activity against G2019S LRRK2 (IC50 = 12.9 nM). Molecular dynamics (MD) simulations revealed that LY2023-001 formed stable hydrogen bonds with Glu1948, and Ala1950 in the G2019S LRRK2 protein. Additionally, its phenyl substituents engage in strong electrostatic interactions with Lys1906 and van der Waals interactions with Leu1885, Phe1890, Val1893, Ile1933, Met1947, Leu1949, Leu2001, Ala2016, and Asp2017. Our findings underscore the potential of computational methods in the successful optimization of small molecules, offering important insights for the development of novel LRRK2 inhibitors.

Decoding Inhibitor Egression from Wild-Type and G2019S Mutant LRRK2 Kinase: Insights into Unbinding Mechanisms for Precision Drug Design in Parkinson's Disease.

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.

A designed ankyrin-repeat protein that targets Parkinson's disease-associated LRRK2.

Leucine rich repeat kinase 2 (LRRK2) is a large multidomain protein containing two catalytic domains, a kinase and a GTPase, as well as protein interactions domains, including a WD40 domain. The association of increased LRRK2 kinase activity with both the familial and sporadic forms of Parkinson's disease has led to an intense interest in determining its cellular function. However, small molecule probes that can bind to LRRK2 and report on or affect its cellular activity are needed. Here, we report the identification and characterization of the first high-affinity LRRK2-binding designed ankyrin-repeat protein (DARPin), named E11. Using cryo-EM, we show that DARPin E11 binds to the LRRK2 WD40 domain. LRRK2 bound to DARPin E11 showed improved behavior on cryo-EM grids, resulting in higher resolution LRRK2 structures. DARPin E11 did not affect the catalytic activity of a truncated form of LRRK2 in vitro but decreased the phosphorylation of Rab8A, a LRRK2 substrate, in cells. We also found that DARPin E11 disrupts the formation of microtubule-associated LRRK2 filaments in cells, which are known to require WD40-based dimerization. Thus, DARPin E11 is a new tool to explore the function and dysfunction of LRRK2 and guide the development of LRRK2 kinase inhibitors that target the WD40 domain instead of the kinase.

Effects of bound nucleotides on the secondary structure, thermal stability, and phosphorylation of Rab3A.

Rab3A is a member of the Rab GTPase family involved in synaptic vesicle trafficking. Recent evidence has demonstrated that Rab3A is phosphorylated by leucine-rich repeat kinase 2 (LRRK2) that is implicated in both familial and sporadic forms of Parkinson's disease (PD), and an abnormal increase in Rab3A phosphorylation has been proposed as a cause of PD. Despite the potential importance of Rab3A in PD pathogenesis, its structural information is limited and the effects of bound nucleotides on its biophysical and biochemical properties remain unclear. Here, we show that GDP-bound Rab3A is preferentially phosphorylated by LRRK2 compared with GTP-bound Rab3A. The secondary structure of Rab3A, measured by circular dichroism (CD) spectroscopy, revealed that Rab3A is resistant to heat-induced denaturation at pH 7.4 or 9.0 regardless of the nucleotides bound. In contrast, Rab3A underwent heat-induced denaturation at pH 5.0 at a lower temperature in its GDP-bound form than in its GTP-bound form. The unfolding temperature of Rab3A was studied by differential scanning fluorimetry, which showed a significantly higher unfolding temperature in GTP-bound Rab3A than in GDP-bound Rab3A, with the highest at pH 7.4. These results suggest that Rab3A has unusual thermal stability under physiologically relevant conditions and that bound nucleotides influence both thermal stability and phosphorylation by LRRK2.

Subpocket Similarity-Based Hit Identification for Challenging Targets: Application to the WDR Domain of LRRK2.

We herewith applied a priori a generic hit identification method (POEM) for difficult targets of known three-dimensional structure, relying on the simple knowledge of physicochemical and topological properties of a user-selected cavity. Searching for local similarity to a set of fragment-bound protein microenvironments of known structure, a point cloud registration algorithm is first applied to align known subpockets to the target cavity. The resulting alignment then permits us to directly pose the corresponding seed fragments in a target cavity space not typically amenable to classical docking approaches. Last, linking potentially connectable atoms by a deep generative linker enables full ligand enumeration. When applied to the WD40 repeat (WDR) central cavity of leucine-rich repeat kinase 2 (LRRK2), an unprecedented binding site, POEM was able to quickly propose 94 potential hits, five of which were subsequently confirmed to bind in vitro to LRRK2-WDR.

Leucine-Rich Repeat Kinases.

Activating mutations in leucine-rich repeat kinase 2 (LRRK2) represent the most common cause of monogenic Parkinson's disease. LRRK2 is a large multidomain protein kinase that phosphorylates a specific subset of the ∼65 human Rab GTPases, which are master regulators of the secretory and endocytic pathways. After phosphorylation by LRRK2, Rabs lose the capacity to bind cognate effector proteins and guanine nucleotide exchange factors. Moreover, the phosphorylated Rabs cannot interact with their cognate prenyl-binding retrieval proteins (also known as guanine nucleotide dissociation inhibitors) and, thus, they become trapped on membrane surfaces. Instead, they gain the capacity to bind phospho-Rab-specific effector proteins, such as RILPL1, with resulting pathological consequences. Rab proteins also act upstream of LRRK2 by controlling its activation and recruitment onto membranes. LRRK2 signaling is counteracted by the phosphoprotein phosphatase PPM1H, which selectively dephosphorylates phospho-Rab proteins. We present here our current understanding of the structure, biochemical properties, and cell biology of LRRK2 and its related paralog LRRK1 and discuss how this information guides the generation of LRRK2 inhibitors for the potential benefit of patients.

Structural insights into the GTP-driven monomerization and activation of a bacterial LRRK2 homolog using allosteric nanobodies.

Roco proteins entered the limelight after mutations in human LRRK2 were identified as a major cause of familial Parkinson's disease. LRRK2 is a large and complex protein combining a GTPase and protein kinase activity, and disease mutations increase the kinase activity, while presumably decreasing the GTPase activity. Although a cross-communication between both catalytic activities has been suggested, the underlying mechanisms and the regulatory role of the GTPase domain remain unknown. Several structures of LRRK2 have been reported, but structures of Roco proteins in their activated GTP-bound state are lacking. Here, we use single-particle cryo-electron microscopy to solve the structure of a bacterial Roco protein (CtRoco) in its GTP-bound state, aided by two conformation-specific nanobodies: NbRoco1 and NbRoco2. This structure presents CtRoco in an active monomeric state, featuring a very large GTP-induced conformational change using the LRR-Roc linker as a hinge. Furthermore, this structure shows how NbRoco1 and NbRoco2 collaborate to activate CtRoco in an allosteric way. Altogether, our data provide important new insights into the activation mechanism of Roco proteins, with relevance to LRRK2 regulation, and suggest new routes for the allosteric modulation of their GTPase activity.

Rab29-dependent asymmetrical activation of leucine-rich repeat kinase 2.

Gain-of-function mutations in LRRK2, which encodes the leucine-rich repeat kinase 2 (LRRK2), are the most common genetic cause of late-onset Parkinson's disease. LRRK2 is recruited to membrane organelles and activated by Rab29, a Rab guanosine triphosphatase encoded in the PARK16 locus. We present cryo-electron microscopy structures of Rab29-LRRK2 complexes in three oligomeric states, providing key snapshots during LRRK2 recruitment and activation. Rab29 induces an unexpected tetrameric assembly of LRRK2, formed by two kinase-active central protomers and two kinase-inactive peripheral protomers. The central protomers resemble the active-like state trapped by the type I kinase inhibitor DNL201, a compound that underwent a phase 1 clinical trial. Our work reveals the structural mechanism of LRRK2 spatial regulation and provides insights into LRRK2 inhibitor design for Parkinson's disease treatment.

Role of the Residue Q1919 in Increasing Kinase Activity of G2019S LRRK2 Kinase: A Computational Study.

Mutations in multi-domain leucine-rich repeat kinase 2 (LRRK2) have been an interest to researchers as these mutations are associated with Parkinson's disease. G2019S mutation in LRRK2 kinase domain leads to the formation of additional hydrogen bonds by S2019 which results in stabilization of the active state of the kinase, thereby increasing kinase activity. Two additional hydrogen bonds of S2019 are reported separately. Here, a mechanistic picture of the formation of additional hydrogen bonds of S2019 with Q1919 (also with E1920) is presented using 'active' Roco4 kinase as a homology model and its relationship with the stabilization of the 'active' G2019S LRRK2 kinase. A conformational flipping of residue Q1919 was found which helped to form stable hydrogen bond with S2019 and made 'active' state more stable in G2019S LRRK2. Two different states were found within the 'active' kinase with respect to the conformational change (flipping) in Q1919. Two doubly-mutated systems, G2019S/Q1919A and G2019S/E1920 K, were studied separately to check the effect of Q1919 and E1920. For both cases, the stable S2 state was not formed, leading to a decrease in kinase activity. These results indicate that both the additional hydrogen bonds of S2019 (with Q1919 and E1920) are necessary to stabilize the active G2019S LRRK2.

Doubly Constrained C-terminal of Roc (COR) Domain-Derived Peptides Inhibit Leucine-Rich Repeat Kinase 2 (LRRK2) Dimerization.

Missense mutations along the leucine-rich repeat kinase 2 (LRRK2) protein are a major contributor to Parkinson's Disease (PD), the second most commonly occurring neurodegenerative disorder worldwide. We recently reported the development of allosteric constrained peptide inhibitors that target and downregulate LRRK2 activity through disruption of LRRK2 dimerization. In this study, we designed doubly constrained peptides with the objective of inhibiting C-terminal of Roc (COR)-COR mediated dimerization at the LRRK2 dimer interface. We show that the doubly constrained peptides are cell-permeant, bind wild-type and pathogenic LRRK2, inhibit LRRK2 dimerization and kinase activity, and inhibit LRRK2-mediated neuronal apoptosis, and in contrast to ATP-competitive LRRK2 kinase inhibitors, they do not induce the mislocalization of LRRK2 to skein-like structures in cells. This work highlights the significance of COR-mediated dimerization in LRRK2 activity while also highlighting the use of doubly constrained peptides to stabilize discrete secondary structural folds within a peptide sequence.

Capturing Differences in the Regulation of LRRK2 Dynamics and Conformational States by Small Molecule Kinase Inhibitors.

Mutations in the human leucine rich repeat protein kinase-2 (LRRK2) create risk factors for Parkinson's disease, and pathological functions of LRRK2 are often correlated with aberrant kinase activity. Past research has focused on developing selective LRRK2 kinase inhibitors. In this study, we combined enhanced sampling simulations with HDX-MS to characterize the inhibitor-induced dynamic changes and the allosteric communications within the C-terminal domains of LRRK2, LRRK2RCKW. We find that the binding of MLi-2 (a type I kinase inhibitor) stabilizes a closed kinase conformation and reduces the global dynamics of LRRK2RCKW, leading to a more compact LRRK2RCKW structure. In contrast, the binding of Rebastinib (a type II kinase inhibitor) stabilizes an open kinase conformation, which promotes a more extended LRRK2RCKW structure. By probing the distinct effects of the type I and type II inhibitors, key interdomain interactions are found to regulate the communication between the kinase domain and the GTPase domain. The intermediate states revealed in our simulations facilitate the efforts toward in silico design of allosteric modulators that control LRRK2 conformations and potentially mediate the oligomeric states of LRRK2 and its interactions with other proteins.

LRRK2 Quantification in Cerebrospinal Fluid of Patients with Parkinson's Disease and Atypical Parkinsonian Syndromes.

BACKGROUND: The alteration of leucine-rich repeat kinase 2 (LRRK2) kinase activity is thought to be involved in Parkinson's disease (PD) pathogenesis beyond familiar cases, and LRRK2 inhibitors are currently under investigation. Preliminary data suggest a relationship between LRRK2 alteration and cognitive impairment in PD. OBJECTIVE: To investigate cerebrospinal fluid (CSF) LRRK2 levels in PD and other parkinsonian disorders, also correlating them with cognitive impairment. METHODS: In this study, we retrospectively investigated by means of a novel highly sensitive immunoassay the levels of total and phosphorylated (pS1292) LRRK2 in CSF of cognitively unimpaired PD (n = 55), PD with mild cognitive impairment (n = 49), PD with dementia (n = 18), dementia with Lewy bodies (n = 12), atypical parkinsonian syndromes (n = 35), and neurological controls (n = 30). RESULTS: Total and pS1292 LRRK2 levels were significantly higher in PD with dementia with respect to PD with mild cognitive impairment and PD, and also showed a correlation with cognitive performances. CONCLUSIONS: The tested immunoassay may represent a reliable method for assessing CSF LRRK2 levels. The results appear to confirm an association of LRRK2 alteration with cognitive impairment in PD. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Identification of LRRK2 Inhibitors through Computational Drug Repurposing.

Parkinson's disease (PD) is the second most common neurodegenerative disorder that affects more than ten million people worldwide. However, the current PD treatments are still limited and alternative treatment strategies are urgently required. Leucine-rich repeat kinase 2 (LRRK2) has been recognized as a promising target for PD treatment. However, there are no approved LRRK2 inhibitors on the market. To rapidly identify potential drug repurposing candidates that inhibit LRRK2 kinase, we report a structure-based drug repurposing workflow that combines molecular docking, recursive partitioning model, molecular dynamics (MD) simulation, and molecular mechanics-generalized Born surface area (MM-GBSA) calculation. Thirteen compounds screened from our drug repurposing workflow were further evaluated through the experiment. The experimental results showed six drugs (Abivertinib, Aumolertinib, Encorafenib, Bosutinib, Rilzabrutinib, and Mobocertinib) with IC50 less than 5 µM that were identified as potential LRRK2 kinase inhibitors. The most potent compound Abivertinib showed potent inhibitions with IC50 toward G2019S mutation and wild-type LRRK2 of 410.3 nM and 177.0 nM, respectively. Our combination screening strategy had a 53% hit rate in this repurposing task. MD simulations and MM-GBSA free energy analysis further revealed the atomic binding mechanism between the identified drugs and G2019S LRRK2. In summary, the results showed that our drug repurposing workflow could be used to identify potent compounds for LRRK2. The potent inhibitors discovered in our work can be a starting point to develop more effective LRRK2 inhibitors.

Structure-Based Virtual Screening and De Novo Design to Identify Submicromolar Inhibitors of G2019S Mutant of Leucine-Rich Repeat Kinase 2.

Missense mutations of leucine-rich repeat kinase 2 (LRRK2), including the G2019S mutant, are responsible for the pathogenesis of Parkinson's disease. In this work, structure-based virtual screening of a large chemical library was carried out to identify a number of novel inhibitors of the G2019S mutant of LRRK2, the biochemical potencies of which ranged from the low micromolar to the submicromolar level. The discovery of these potent inhibitors was made possible due to the modification of the original protein-ligand binding energy function in order to include an accurate ligand dehydration energy term. The results of extensive molecular docking simulations indicated that the newly identified inhibitors were bound to the ATP-binding site of the G2019S mutant of LRRK2 through the multiple hydrogen bonds with backbone amide groups in the hinge region as well as the hydrophobic interactions with the nonpolar residues in the P-loop, hinge region, and interdomain region. Among 18 inhibitors derived from virtual screening, 4-(2-amino-5-phenylpyrimidin-4-yl)benzene-1,3-diol (Inhibitor 2) is most likely to serve as a new molecular scaffold to optimize the biochemical potency, because it revealed submicromolar inhibitory activity in spite of its low molecular weight (279.3 amu). Indeed, a highly potent inhibitor (Inhibitor 2n) of the G2019S mutant was derived via the structure-based de novo design using the structure of Inhibitor 2 as the molecular core. The biochemical potency of Inhibitor 2n surged to the nanomolar level due to the strengthening of hydrophobic interactions in the ATP-binding site, which were presumably caused by the substitutions of small nonpolar moieties. Due to the high biochemical potency against the G2019S mutant of LRRK2 and the putatively good physicochemical properties, Inhibitor 2n is anticipated to serve as a new lead compound for the discovery of antiparkinsonian medicines.

Roc, the G-domain of the Parkinson's disease-associated protein LRRK2.

Mutation in leucine-rich repeat (LRR) kinase 2 (LRRK2) is a common cause of Parkinson's disease (PD). Aberrant LRRK2 kinase activity is associated with disease pathogenesis and thus it is an attractive drug target for combating PD. Intense efforts in the past nearly two decades have focused on the development of small-molecule inhibitors of the kinase domain of LRRK2 and have identified potent kinase inhibitors. However, most LRRK2 kinase inhibitors have shown adverse effects; therefore, alternative-mechanism-based strategies are desperately needed. In this review, we discuss the new insights gleaned from recent cryoelectron microscope (cryo-EM) structures of LRRK2 towards understanding the mechanisms of actions of LRRK2 and explore the potential new therapeutic avenues.

LRRK2 dynamics analysis identifies allosteric control of the crosstalk between its catalytic domains.

The 2 major molecular switches in biology, kinases and GTPases, are both contained in the Parkinson disease-related leucine-rich repeat kinase 2 (LRRK2). Using hydrogen-deuterium exchange mass spectrometry (HDX-MS) and molecular dynamics (MD) simulations, we generated a comprehensive dynamic allosteric portrait of the C-terminal domains of LRRK2 (LRRK2RCKW). We identified 2 helices that shield the kinase domain and regulate LRRK2 conformation and function. One helix in COR-B (COR-B Helix) tethers the COR-B domain to the αC helix of the kinase domain and faces its activation loop, while the C-terminal helix (Ct-Helix) extends from the WD40 domain and interacts with both kinase lobes. The Ct-Helix and the N-terminus of the COR-B Helix create a "cap" that regulates the N-lobe of the kinase domain. Our analyses reveal allosteric sites for pharmacological intervention and confirm the kinase domain as the central hub for conformational control.

LRRK2 at Striatal Synapses: Cell-Type Specificity and Mechanistic Insights.

Mutations in leucine-rich repeat kinase 2 (LRRK2) cause Parkinson's disease with a similar clinical presentation and progression to idiopathic Parkinson's disease, and common variation is linked to disease risk. Recapitulation of the genotype in rodent models causes abnormal dopamine release and increases the susceptibility of dopaminergic neurons to insults, making LRRK2 a valuable model for understanding the pathobiology of Parkinson's disease. It is also a promising druggable target with targeted therapies currently in development. LRRK2 mRNA and protein expression in the brain is highly variable across regions and cellular identities. A growing body of work has demonstrated that pathogenic LRRK2 mutations disrupt striatal synapses before the onset of overt neurodegeneration. Several substrates and interactors of LRRK2 have been identified to potentially mediate these pre-neurodegenerative changes in a cell-type-specific manner. This review discusses the effects of pathogenic LRRK2 mutations in striatal neurons, including cell-type-specific and pathway-specific alterations. It also highlights several LRRK2 effectors that could mediate the alterations to striatal function, including Rabs and protein kinase A. The lessons learned from improving our understanding of the pathogenic effects of LRRK2 mutations in striatal neurons will be applicable to both dissecting the cell-type specificity of LRRK2 function in the transcriptionally diverse subtypes of dopaminergic neurons and also increasing our understanding of basal ganglia development and biology. Finally, it will inform the development of therapeutics for Parkinson's disease.

Closing the structure-to-function gap for LRRK2.

Variations in the LRRK2 gene represent one of the strongest genetic factors for Parkinson's disease (PD). It has become clear that structural knowledge of the encoded large multidomain LRRK2 protein will cast light on its biological function. The new study from Myasnikov, Zhu, et al. provides a high-resolution structure of the full-length LRRK2.

From structure to ætiology: a new window on the biology of leucine-rich repeat kinase 2 and Parkinson's disease.

Since the discovery of mutations in leucine-rich repeat kinase 2 (LRRK2) as an underlying genetic cause for the development of Parkinson's disease (PD) in 2004 (Neuron 44, 601-607; Neuron 44, 595-600), and subsequent efforts to develop LRRK2 kinase inhibitors as a therapy for Parkinson's (Expert Opin. Ther. Targets 21, 751-753), elucidating the atomic resolution structure of LRRK2 has been a major goal of research into this protein. At over 250 kDa, the large size and complicated domain organisation of LRRK2 has made this a highly challenging target for structural biologists, however, a number of recent studies using both in vitro and in situ approaches (Nature 588, 344-349; Cell 182, 1508-1518.e1516; Cell 184, 3519-3527.e3510) have provided important new insights into LRRK2 structure and the complexes formed by this protein.

Structural analysis of the full-length human LRRK2.

Mutations in leucine-rich repeat kinase 2 (LRRK2) are commonly implicated in the pathogenesis of both familial and sporadic Parkinson's disease (PD). LRRK2 regulates critical cellular processes at membranous organelles and forms microtubule-based pathogenic filaments, yet the molecular basis underlying these biological roles of LRRK2 remains largely enigmatic. Here, we determined high-resolution structures of full-length human LRRK2, revealing its architecture and key interdomain scaffolding elements for rationalizing disease-causing mutations. The kinase domain of LRRK2 is captured in an inactive state, a conformation also adopted by the most common PD-associated mutation, LRRK2G2019S. This conformation serves as a framework for structure-guided design of conformational specific inhibitors. We further determined the structure of COR-mediated LRRK2 dimers and found that single-point mutations at the dimer interface abolished pathogenic filamentation in cells. Overall, our study provides mechanistic insights into physiological and pathological roles of LRRK2 and establishes a structural template for future therapeutic intervention in PD.