Abstract
Mutations in Leucine-Rich Repeat Kinase 2 (LRRK2) are a leading contributor to developing familial and idiopathic Parkinson’s disease (PD). Most PD-causing LRRK2 mutations increase the kinase activity, leading to increased phosphorylation of Rab GTPases, disrupting vesicular trafficking, cytoskeletal dynamics, and autophagy. Under homeostatic conditions, the bulk of WT and PD-mutant LRRK2 is found in the cytosol. However, exogenously expressed LRRK2 can form microtubule-associated filaments that have been shown to affect molecular transport along microtubules in vitro. While the physiological relevance of microtubule binding has not been established yet, inhibitors being designed and tested as therapeutics have been shown to either promote or prevent filament formation of LRRK2. In this study, we examine the localization and resulting molecular organization of hyperactive LRRK2-I2020T, a common PD mutant, in cells treated with type I (MLi-2) or type II (GZD-824) kinase inhibitors. Treatment with a type I kinase inhibitor results in extensive LRRK2-I2020T decoration around microtubules and microtubule bundling. Stabilization of LRRK2-I2020T filaments by type I inhibitor treatment allowed us to build a full-length closed-kinase model of LRRK2-I2020T in its cellular environment. Conversely, treatment with a type II inhibitor resulted in minimal microtubule decoration by LRRK2-I2020T compared to Type I inhibitor treated cells. This study provides a structural framework for understanding how type I and type II kinase inhibitors differentially modulate LRRK2 filament formation, demonstrating that type I inhibitor treatment promotes a distinct filament architecture, whereas such assemblies are not observed with type II inhibitors.
Introduction
Parkinson’s disease (PD) is a chronic and progressive neurodegenerative disorder affecting approximately ten million people worldwide1–4. Genome-wide association studies have identified leucine-rich repeat kinase 2 (LRRK2) as a major risk factor for developing familial and idiopathic PD5–7. LRRK2 contributes to disrupted cytoskeletal dynamics, lysosomal dysfunction, neuroinflammation and impaired dopaminergic signaling8–11.
LRRK2 is a large multidomain protein, one in four in the human genome to contain a GTPase and a kinase12–14. Its C-terminal catalytic region is composed of WD40, ROC GTPase, Kinase and COR (RCKW) domains (Fig. 1A). Pathogenic mutations in this region (e.g. I2020T, G2019S, R1441H/G/C/S, Y1699C) increase the kinase activity leading to hyperphosphorylation of Rab GTPases disrupting fusion, docking and trafficking of vesicular membranes15–18. The N-terminal region mediates protein-protein interactions, facilitating recognition of LRRK2 substrates through the Leucine rich repeats (LRR), Ankyrin (ANK) and Armadillo (ARM) domains (Fig. 1A). In the autoinhibited inactive kinase state, the N-terminus occludes catalytic domains and limits substrate access to the kinase domain19–21. Upon activation, these N-terminal repeats are proposed to undock, enabling substrate binding near the kinase moiety and subsequent phosphorylation of the target19. Recent structural studies have elucidated the architecture of LRRK2’s catalytic RCKW domains in both active and inactive conformations, as well as full-length autoinhibited and intermediate activation state structures of the protein19– 27. Resolving the full-length closed-kinase structure of LRRK2 remains a significant challenge due to the large conformational heterogeneity of the N-terminal repeats in the active state.

LRRK2IT expression in 293T cells responds differently to treatment with type I and type II kinase inhibitors.
(A) LRRK2 contains multiple functional domains including Armadillo (ARM), Ankyrin repeats (ANK), Leucine-rich repeats (LRR), Ras of complex (ROC), COR, kinase and WD40 domain. (B) Schematic of the LRRK2 protein construct used in this paper that includes a hyperactive I2020T mutation and an N-terminal GFP-tag. (C) Type I inhibitors (red) capture the LRRK2 kinase in active, closed conformation while type II inhibitors (blue) trap kinase of LRRK2 in an inactive, open kinase conformation. (D-E) Cryo-fluorescence microscopy (cryo-FM) image of a vitreous cell on an electron microscopy grid expressing LRRK2IT treated with MLi-2. The GFP-LRRK2IT signal (green) is present as filaments as well as puncta. (F-G) Cryo-FM image of a cell expressing LRRK2IT treated with GZD-824. The GFP-LRRK2IT signal (green) shows a punctate distribution of the protein throughout the cytosol of the cell. (H-I) Cryo-TEM overview of the lamella of a cell expressing LRRK2IT treated with MLi-2 (H) and GZD-824 (I). Microtubule bundles are highlighted in green with some displaying dense LRRK2 decoration. Green: microtubule bundles, red: mitochondria, pink: vesicles, blue: plasma membrane, purple: endoplasmic reticulum; Scale bars: D, F 20μm; E, G 5μm; H, I 500nm.
Given its role in PD pathogenesis, LRRK2 has been the focus of small molecule kinase inhibitor development with two drug candidates currently in clinical testing28–31. There are two major classes of LRRK2 inhibitors: type I LRRK2 specific inhibitors (e.g. MLi2, LRRK2-IN-1, HG-10-102-01) stabilize the kinase domain in its closed active-like conformation, while type II broad spectrum inhibitors (e.g. GZD-824, ponatinib, rebastinib) stabilize the kinase in its open, inactive-like conformation31–37 (Fig. 1C). Both types of inhibitors bind to the kinase domain, but type I inhibitors compete for ATP binding while type II inhibitors bind to the allosteric site adjacent to the kinase domain19,21. The conformational changes of the kinase domain when treated with different types of inhibitors have been characterized, but the effect of inhibitor treatment on LRRK2 localization and molecular architecture in cells remains unclear19,21.
In cells, LRRK2 is typically cytosolic, but under certain conditions it binds to microtubules and forms extensive filaments22,38,39. Association with microtubules disrupts molecular transport along microtubules in vitro23. Hyperactivating pathogenic mutations such as I2020T and R1441H/G/C/S increase the protein’s propensity to form microtubule associated filaments22–24. Additionally, type I, but not type II, kinase inhibitor treatment increases LRRK2 filamentation in cells24,40. While physiological significance of LRRK2 filamentation remains unclear, microtubule association provides an assay to test effects of kinase inhibitors on molecular architecture of LRRK2 in cellular models.
In this study, we investigate the effects of different types of kinase inhibitors on the cellular architecture and localization of PD-associated LRRK2-I2020T (LRRK2IT). To this end, we expressed GFP-tagged LRRK2IT in 293T cells and treated samples with either type I (MLi-2) or type II (GZD-824) inhibitors. To visualize effects of inhibitor treatment on LRRK2 architecture at nanometer resolution we used cryo-correlative light and electron microscopy (cryo-CLEM), which combines cryo-electron tomography with cryo-fluorescence microscopy (cryo-FM)41–43.
In cells treated with MLi-2, we observed highly ordered LRRK2-decorated microtubules that formed extensive microtubule bundles of 15-30 microtubules. In contrast, LRRK2 microtubule association was rarely observed in GZD-824 treated cells. These observations point to distinct mechanisms by which type I and type II kinase inhibitors influence LRRK2 filamentation in a cellular context. Besides increased filament formation, MLi-2 treatment reduced conformational heterogeneity of LRRK2IT, allowing us to characterize the architecture of the full-length microtubule-associated conformation of LRRK2IT. We present the structure of full-length LRRK2IT in a closed-kinase conformation and N-terminal repeats undocked from the RCKW domains, providing insight about hyperactive protein architecture in the cellular environment.
Results
LRRK2IT localizes to microtubules in the presence of a type I inhibitor and to the cytosol with a type II inhibitor
We used cryo-CLEM to study how cellular localization of GFP-tagged hyperactive LRRK2IT is affected by treatment with either type I (MLi-2) or type II (GZD-824) inhibitors (Fig. 1; Ext. Fig. 1). Our previous work showed that expression of LRRK2IT in 293T cells forms an extensive network of filaments around microtubules22. Here we added either MLi-2 or GZD-824 inhibitors that resulted in significant phenotypic differences in LRRK2IT localization (Fig. 1D-G; Ext. Fig. 1A-J). In cells treated with MLi-2, we observed LRRK2IT in extended filaments, puncta, and diffuse in the cytosol (Fig. 1D-E; Ext. Fig 1A-D). In contrast, when cells were treated with GZD-824, LRRK2IT was mostly localized to puncta and distributed throughout the cytosol (Fig. 1F-G; Ext. Fig 1E-H), in agreement with our previous work23,24,40.
To probe this difference at the structural level, we used cryo-FIB milling to produce thin samples of vitrified cells for cryo-ET. Cryo-FM images were aligned to cryo-TEM overviews of the lamella to localize LRRK2IT (Fig. 1H-I; Ext. Fig. 1I-J). In cells treated with MLi-2, we co-localized LRRK2IT with an extended network of microtubules (Fig. 1H-I; Ext. Fig. I-J). These microtubules were coated with a dense LRRK2IT lattice and bundled. The LRRK2IT decorated microtubule network was extensive, spanning 2.17±1.40µm in length (Fig. 1H; Ext. Fig. 1K). By contrast, microtubule bundling was not present in GZD-824 treated cells (Fig. 1I).
Type I kinase inhibitor treatment promotes bundling of LRRK2IT decorated microtubules
We acquired cryo-electron tomograms of MLi-2 and GZD-824 treated cells targeting LRRK2IT localization sites and identified all microtubules within the data. In MLi-2 treated cells, we focused on microtubules that were decorated with a LRRK2IT lattice (Fig. 2A-D; Ext. Fig. 2A-B; Ext. Fig 3). We found a substantial fraction of them to be surrounded by other decorated microtubules often closely spaced and arranged in parallel, suggesting that LRRKIT binding promotes microtubule bundling (Fig. 2A-D; Ext. Fig. 2A-B). Bundles of LRRK2IT decorated microtubules contained an average of 16±8.04 microtubules (range 5-30; Fig. 3M). The center-to-center distance between LRRK2IT-decorated microtubules was 68.7nm±5.73nm, while undecorated microtubule bundles had a center-to-center distance of 35nm±3.00nm in our data and in previous findings (Fig. 3L)22. The expanded spacing observed in the LRRK2IT decorated microtubules is due to the 15nm thick LRRK2IT layer that coats each microtubule, adding 30nm to the center-to-center distance (Fig. 3A-C, L).

Treatment with the type I kinase inhibitor MLi-2 promotes formation of higher-order LRRK2IT-decorated microtubule bundles, whereas the type II kinase inhibitor GZD-824 reduces LRRK2IT filamentation and microtubule bundling in cells.
(A, C) Representative tomogram slices through the cytosol of a cell treated with MLi-2, showing a bundle of microtubules decorated with LRRK2IT. (B, D) Segmentation of the tomogram in A, C highlighting bundles of microtubules decorated with a well-ordered LRRK2IT lattice. (E) Tomogram slice of a cell treated with GZD-824, showing sparse LRRK2IT decoration on microtubules. (F) Segmentation of the tomogram in E. LRRK2IT decoration is present on three microtubules while there are nine undecorated microtubules present in the proximity. (G) Representative image of a cell treated with GZD-824. Single microtubule is decorated with LRRK2IT. This distribution pattern contrasts sharply with the extensive microtubule-associated filaments observed following Type I inhibitor treatment. (H) Segmentation of a tomogram in G. of a cell treated with GZD-824. Six undecorated microtubules are present near a single LRRK2 decorated microtubule. Lime Green: LRRK2IT, green: microtubule, red: mitochondria, light pink: vesicles, blue: multivesicular bodies, pink: ribosomes. Scale bars: 100nm.

LRRK2IT lattice that decorates microtubules is highly ordered in cells treated with type I as compared to type II inhibitor.
(A) Segmented representation of a LRRK2IT associated microtubule bundle in MLi-2 treated cells. Individual LRRK2IT subunits are highly ordered around the microtubules. (B) Top-down view of the segmented representation of a LRRK2IT associated microtubule bundle in MLi-2 treated cells. The individual microtubules are separated from each other by 65-70nm. (C) Segmented representation of a LRRK2IT associated microtubule bundle in MLi-2 treated cells. (D-E) Close-up view of the LRRK2IT decoration on microtubules in MLi-2 treated cells. (F-G) Segmented representation of LRRK2IT associated microtubules in GZD-824 treated cells. LRRK2IT lattice is less ordered compared to MLi-2 treated lattice in (D-E). (continued) (H-I) Subtomogram average of the LRRK2IT associated with microtubules from MLi-2 (H) or GZD-824 (I) treated cells. MLi-2 treated cells exhibit a more ordered lattice of LRRK2IT than the lattice found in GZD-824 treated cells. (J-K) Nearest neighbor analysis of a central LRRK2IT subunit (highlighted in black) in MLi-2 and GZD-824 treated cells. In MLi-2 treated cells, the central LRRK2IT subunit has 16 nearest neighbors highlighted in orange. While in GZD-824 treated cells, the central LRRK2IT subunit has only 8 nearest neighbors shown in pink. (L) Evaluation of the center-to-center distance of undecorated and LRRK2IT-decorated microtubule bundles in MLi-2 and GZD-824 treated cells. (MLi-2 n=80, GZD-824 n=7, undecorated n=13; vertical line at median). (M) Quantification of LRRK2IT decorated microtubule (MT) bundles and isolated microtubules in MLi-2 and GZD-824 treated cells. (bundles quantified in MLi-2 treated cells n=12, isolated MTs quantified in MLi-2 treated cells n=16, isolated MTs quantified in GZD-824 n=6, bundle of MTs quantified in GZD-824 bundle n=1, vertical line at median). Scale bars: A-C 40nm; D-G 20nm; H-K 10nm.
In contrast, tomograms acquired from type II GZD-824 inhibitor treatment resulted in reduced LRRK2 filamentation around microtubules, with only a few rare, isolated microtubules showing LRRK2IT decoration (Fig. 2E-H; Ext. Fig. 2C-D). Interestingly, we captured a single bundle of LRRK2IT decorated microtubules in the tomograms of GZD-824 treated cells (Ext. Fig. 4). This bundle contained 15 microtubules with intermicrotubule spacing similar to LRRK2IT decorated microtubule bundles in MLi-2 treated cells. Otherwise, fluorescence microscopy of GZD-824 treated cells showed LRRK2IT was present as discrete cytosolic puncta, suggesting association with intracellular membranes (Fig. 1F-G, I; Ext. Fig. 1E-H, J). In these data, organelles appeared intact and morphologically unaltered, with no discernible densities, or organized assemblies that could be confidently attributed to LRRK2IT (Fig. 1I). Aside from the absence of microtubule bundling, the overall ultrastructure of GZD-824–treated cells closely resembled that of MLi-2–treated and of untreated cells22.
LRRK2IT microtubule decoration is highly ordered in cells treated with type I compared to type II inhibitors
To characterize the differences between LRRK2IT lattice organization in cells treated with different inhibitors, we performed subtomogram analysis of decorated microtubules and mapped the location of each individual LRRK2 dimer and quantified lattice parameters (Fig. 3). In the MLi-2 treated cells, LRRK2IT strands were organized around microtubules with a regularly spaced lattice, similar to the LRRK2IT strands in cells not treated without the inhibitor (Fig. 3A-E)22. Conversely, in GZD-824-treated cells, the LRRK2IT lattice was rare, and when existent, it was disordered (Fig. 3F-G). While average pitch. rise, and handedness of the filaments of the rare GZD-824 treated LRRK2 filaments were similar between the MLi-2 parameters (Fig. 3H-K), the relative orientation of individual dimers was variable, and the lattice was highly irregular of in GZD-824 treated cell (Fig. 3F-G). Using nearest neighbor analysis, we found that in the LRRK2IT lattice from MLi-2-treated cells, up to 16 close neighbors were placed with high precision, compared to only eight in the GZD-824-treated LRRK2IT lattice (Fig. 3J-K). Altogether this suggests that MLi-2 treatment stabilizes LRRK2 filamentation on microtubules, while GZD-824 disrupts LRRK2IT lattice weakening its association with microtubules.
Subtomogram analysis in type I inhibitor treated cells reveals the architecture of the N-terminal repeats of the full-length LRRK2 in closed-kinase conformation
Full-length LRRK2 structures with closed-kinase conformation have been determined, but the flexible N-terminal repeats remain disordered and are not resolved in these structures19–27. Given the stabilization of LRRK2IT on microtubules with MLi-2 treatment we attempted to resolve full-length closed-kinase conformation of LRRK2IT focusing on resolving the N-terminal repeat architecture.
First, we performed subtomogram analysis on regions of the microtubules that focused on a LRRK2IT dimer and its surroundings (Fig. 4). We focused refinements on the COR-COR dimerization interface, centering the reconstruction around four LRRK2 monomers, resulting in a 12Å resolution map (map A; Fig. 4A). We also refined the WD40-WD40 dimerization interface and focused refinements on four central subunits reaching 11Å resolution (map B; Ext. Fig. 5). These reconstructions allowed us to dock all the individual domains of LRRK2 into the map (Fig. 4). To determine the conformation of the kinase domain of LRRK2IT we fit the structures of active and inactive kinase of LRRK2 into map A (Fig. 4A-C). Fitting the inactive kinase structure resulted in significant clashes between COR-COR domains of neighboring LRRK2IT subunits. The closed-kinase structure of the catalytic domain fit in map A (Fig. 4A) with no clashes, confirming that the reconstruction represents LRRK2IT in its closed-kinase conformation (Fig. 4B-C). In maps A and B, we resolved a connection between a ROC domain and tubulin, present in the in-vitro microtubule-associated LRRK2 RCKW structure, but previously elusive in the in-situ maps (Fig. 4G, Ext. Fig. 5G-H)23,24.

Subtomogram analysis of LRRK2 MLi-2 reveals structural details about a N-terminal domains of closed-kinase LRRK2IT
(A) Subtomogram average map of the microtubule-bound LRRK2 from MLi-2 treated cells (threshold value 0.24). (B-C) Molecular model of two central protomers of active kinase LRRK2 RCKW fit into map A. The protomers of LRRK2 fit well within the density corresponding to the active conformation of the LRRK2 kinase. The domains are colored according to Fig. 1A. (D) Subtomogram average of LRRK2 displaying densities for the N-terminal domains emanating from the catalytic half of the protein (threshold value 0.14). (E-F) Molecular model of the two central protomers of full-length LRRK2 fit into the map shown in D. The model includes RCKW domains and N-terminal LRR-ANK domains. (G) Molecular model of a single subunit of full length LRRK2 in active conformation interacting with a microtubule (gray) via ROC domain (green). (continued) (H) Comparison between the RCKW domains of LRRK2 in intermediate active kinase conformation (gray) and fully active full-length in-cell model of LRRK2 (colored). The kinase domain in the cellular map is fully closed, demonstrated by 15Å shift between the COR domains. (I-J) Comparison of the LRR and ARM domains of LRRK2 in the in-cell model (colored) versus the intermediate kinase-active LRRK2 structure (gray). The COR– ROC domains are aligned in both structures to highlight the positional shifts in the LRR and ARM domains (K-L) Subtomogram average of microtubule-bound LRRK2 from GZD-824–treated cells. (M-N) Molecular model of the active kinase RCKW domains fitted into the map shown in panel K. The molecular model aligns well with the active kinase conformation of LRRK2. The N-terminal domains are not visible in this map. Rotated views are shown. Scale bars: A-F 5nm; G-J 2nm; K-L 10nm; M-N 5nm.
Notably, our maps have ordered density for the N-terminal LRR-ANK domains (Fig. 4D-F), and thus we set out to build a full-length model of closed-kinase LRRK2IT. We based this model on our published structure of in-vitro MLi-2 treated full-length LRRK2, where the kinase was captured in an intermediate activation state with the N-terminal repeats docked to the kinase, preventing it from adopting a fully active conformation19. We split this model in three parts: the WD40 and C-lobe of the kinase, the N-lobe of the kinase with ROC and COR domains, and the LRR and ANK domains, aligned and fitted each of these three to our map A (Fig. 4D-F). When comparing this structure with the reference model of LRRK2IT in its intermediate active conformation, we discerned key differences between the structures. First, the kinase was in the fully active conformation in the in-situ map (Fig. 4G-H). Second, the LRR-ANK domains were undocked from the kinase and WD40 domains and shifted by 19.8Å as compared to the intermediate kinase-active structure (Fig. 4I-J). In the closed-kinase form, the N-terminal repeats are known to undock from the catalytic core (Fig. 4I-J), allowing easier access of LRRK2 substrates to the kinase domain. Here, we captured and modeled for the first time the N-terminal repeats undocked from the catalytic domains with kinase in an active-like conformation.
In map C (Fig. 5A), we discerned an additional layer of protein decoration that accommodated not only LRR and ANK but also the ARM domain of LRRK2IT (Fig. 5A-C). In this map, we first fit in our in-situ model of LRRK2IT and fit an AlphaFold2 model of the ARM domain into the additional densities (Fig. 5D-F). While resolution of this map is moderate, densities that emanate from the catalytic part of the protein clearly correspond to the LRR-ANK-ARM domains (Fig. 5D-F). LRR-ANK-ARM domains all point outward, extending away from the catalytic domains of LRRK2IT (Fig. 5D-F). We hypothesize that the bundling visible in our maps could be reinforced by the connections between the LRRK2IT N-terminal ARM domains between adjacent microtubules in a bundle (Fig. 3E-F; Fig. 5D-F).

Model of full-length LRRK2IT associated with microtubules in its active-like conformation.
(A) Subtomogram average of LRRK2IT lattice decorating a microtubule. (B) Cross-section view of the LRRK2IT lattice. Catalytic domains are closest to the microtubule while outer layer represents the N-terminal protein-protein interaction domains of LRRK2. (C) Side view of the microtubule demonstrating the extension of the N-terminal layer of LRRK2. (D-F) Atomic model fit of full-length LRRK2IT in the subtomogram average map. LRRK2 domain colors are as in Fig. 1A. Scale bars: A-C 10nm; D-F 10nm.
Our attempts to characterize the molecular architecture of LRRK2IT following treatment with GZD-824 showed minimal filamentation of LRRK2IT on microtubules (Fig. 2E-H). In tomograms of GZD-824-treated cells, only six individual LRRK2-decorated microtubules and a single LRRK2 decorated microtubule bundle was detected (Fig. 2E-H, Ext. Fig. 4). Subtomogram analysis revealed a disrupted LRRK2 strand geometry, compared to that observed in MLi-2 treated cells (Fig. 3F-G). We obtained a subtomogram map of LRRK2IT decoration at significantly lower resolution (∼20 Å), in which only the RCKW domains could be resolved (Fig. 4K–N). Fitting of the RCKW domains into this map was consistent with a closed-kinase conformation of LRRK2IT (Fig. 4M–N). A likely explanation for this observation is that GZD-824 treatment was insufficient to convert all pre-existing, microtubule-bound LRRK2 molecules from an active to an inactive kinase conformation. Instead, the residual filaments detected are composed of LRRK2 molecules that formed prior to inhibitor addition and are refractory to subsequent GZD-824 binding. This interpretation is consistent with the markedly reduced frequency of filamentation observed under GZD-824 treatment.
Discussion
Utilizing cellular structural biology tools, we dissected the effects of different types of inhibitors on the localization and architecture of hyperactive LRRK2IT expressed in 293T cells. In cells treated with type I MLi-2 inhibitor, we captured a significant number of microtubules with LRRK2IT decoration, the majority of which formed extensive bundled networks of LRRK2IT filaments (Fig. 3). On the other hand, when cells were treated with GZD-824, only a few microtubules were decorated with LRRK2IT (Fig. 2, E-H). Differences in LRRK2IT propensity to form filaments can be explained by the distinct molecular mechanisms of action of inhibitor treatments. MLi-2 treatment is known to stabilize LRRK2 in its closed-kinase conformation which in turn promotes filament formation around microtubules. These observations are in stark contrast to LRRK2 localization in cells treated with GZD-824 which stabilizes LRRK2’s kinase in its inactive conformation weakening association with microtubules19,21. These findings confirm that inhibitor-induced transitions of LRRK2IT from its active to inactive kinase conformation modulates its interaction with microtubules, leading to a shift from extensive to minimal microtubule decoration in cells.
The physiological relevance of LRRK2 filamentation around microtubules remains controversial. Wild-type LRRK2 does not exhibit enriched microtubule localization, whereas several Parkinson’s disease–associated mutants (e.g., I2020T, R1441H/G/S/T) show robust microtubule association. Microtubule decoration by LRRK2IT has not been studied in cell types that endogenously express high levels of LRRK2, such as lung epithelial cells and brain-resident immune cells including microglia and macrophages44. Thus, it remains possible that aberrant LRRK2–microtubule interactions occur under physiological expression conditions, potentially disrupting homeostatic intracellular transport and being further exacerbated by type I LRRK2 inhibitors, as suggested by in vitro studies23,45.
In our previous work, we used in-cell cryo-ET to build and integrative model of the microtubule-associated catalytic ROC, COR, kinase and WD40 (RCKW) domains of LRRK2 in its closed-kinase conformation22. However, our model could not resolve the conformation of N-terminal LRR, ANK and ARM domains. Full-length LRRK2 has been captured in vitro in its intermediate active conformation19– 21,23–26. However, how the N-terminal repeats of LRRK2 are organized when the protein is in its closed-kinase conformation remained elusive. Stabilization of LRRK2 in a closed-kinase conformation by MLi-2 treatment and microtubule association reduces conformational heterogeneity to permit structure determination of full-length LRRK2IT with the N-terminal repeats undocked from the catalytic core. This structure represents one stabilized closed-kinase state within a broader conformational ensemble of undocked N-terminal repeats.
Our results clarify the relationship between kinase conformation, repeat undocking, and microtubule association. Increased microtubule association observed for I2020T mutant favors repeat undocking, a prerequisite for kinase closure and filament assembly. In contrast, the G2019S mutation enhances catalytic activity without increasing microtubule association. Microtubule binding therefore reflects conformational state rather than kinase activity per se.
Together, these findings provide a structural view of full-length LRRK2 in a closedkinase conformation and capture a resolved snapshot along its conformational continuum. Our results underscore how pathogenic mutations and distinct inhibitor classes reshape its structural landscape. More broadly, this work demonstrates the power of in-cell structural biology to illuminate dynamic regulatory transitions that are otherwise obscured by conformational heterogeneity.
Methods
Cell Culture
293T cells were cultured in Dulbecco’s Modified Eagle’s Medium (DMEM) formulated with high glucose. DMEM was supplemented by 10% (v/v) heat inactivated fetal bovine serum (FBS), 0.1mM Non-Essential Amino Acids (NEAA), 1mM Sodium Pyruvate, 1% penicillin-streptomycin (Pen-Strep), 500µg/ml Geneticin Selective antibiotic. Cells were maintained at 37ºC and 5% CO2 subcultured every 2-3 days at ∼80% confluency. Cells were detached for subculture using 0.05% trypsin-EDTA.
Transfection Experiments
293T cells were transfecting using the Lipofectamine 3000 Transfection Kit (Thermo Fisher Scientific) according to the manufacturer’s protocol. For each transfection, 1µg of GFP-tagged LRRK2 recombinant DNA was used (Addgene). Transfections were carried out in 6-well plates. Cells were incubated for 48 hours post-transfection at 37ºC and 5% CO2. Prior to sample preparation for fixation or cryo-EM grid preparation cells were treated with 500nM of MLi-2 for 2 hours, or 300nM GZD-824 for 10-30 minutes. Following treatment, transfected cells were either fixed for fluorescence microscopy analysis or processed for cryo-ET grid preparation.
Cryo-TEM Grid Preparation
Au Quantifoil R1/4 200 mesh grids were used to prepare cellular samples. Prior to cell deposition the grids were patterned to enhance cell adhesion and positioning for cryo-ET46. In brief, grids were plasma cleaned and were treated with PLL-PEG. Using UV laser custom defined micropatterns were created in the layer of PLL-PEG creating adhesive regions using the PRIMO micropatterning system (Alveole). The patterned grids were incubated with fibronectin for 1 hour at room temperature before cells were deposited on them. Detailed protocol for patterning grids can be accessed here47.
Cells were detached from culture dishes after transfection via trypsin treatment, counted and approximately 200,000 cells were deposited on four grids in a 3cm cell culture dish. The cells were then incubated on the patterned grids for 2 hours or until the cells were completely attached. Following attachment to the grids, samples were plunge frozen using a manual custom plunger assembly. The freezing medium consisted of 50/50 ethane/propane mixture. Grids were blotted for 6-8 seconds in a humidity-controlled room. Plunge frozen grids were secured in autogrids and stored in liquid nitrogen until imaging.
Cryo-Focused Ion Beam Milling
Cryo-FIB milling was performed on an Aquilos 2+ system from Thermo Fisher Scientific. Scanning electron microscope images were taken using either 2KV or 5KV at a 30pA beam current. The lamella were prepared at 11º stage angle. Milling was performed in sequential manner starting from rough milling with 0.5nA-0.1nA currents, followed by intermediate milling at 50-30pA and final polishing using 10pA current. Target lamella thickness was 150-200nm. A detailed protocol for cryo-FIB milling is described in previous publications41,48,49.
Integrated Cryo-Fluorescence Microscopy
Cryo-Fluorescence microscopy was performed using an integrated fluorescence microscopy (iFLM) system on Aquilos 2+. The 20x air objective NA 0.7 was used for imaging. For GFP tagged cellular samples excitation 470nm was used for visualization. Illumination was set to 10% intensity with 200ms exposure time to minimize photodamage. Z-stacks spanning 10-20µm in depth were acquired with step sizes ranging between 500-800nm. This allowed imaging of the whole cellular volumes with adequate resolution and limited photodamage.
Cryo-Electron Tomography
TEM image acquisition was performed on Titan Krios G3 system equipped with a K3 direct electron detector and Gatan energy filter. The milling slot of the autogrid was oriented perpendicular to the loading cartridge. Tilt-Series were acquired using Serial-EM software and Pacetomo automated tilt-series acquisition methods50,51. The images were collected at a target defocus range of 4-5 µm with a pixel size of 1.341 Å. Dose symmetric data acquisition schemes were implemented with a tilt range of ±54° degrees with 3º increments52. The pre-tilt was calculated using Pacetomo and acquisition was initiated from this position. Cumulative electron dose of approximately 140 e−/Å2 was applied across the tilt series.
Tomogram Reconstruction and Annotation
Tilt-series were processed using a Warp-Imod pipeline53–55. Initial preprocessing including frame alignment, gain correction, and motion correction was performed in Warp. Tilt-series were aligned in IMOD. The aligned tilt-series stacks were imported back into Warp 1.0.9/beta for 3D-CTF correction. Final tomograms were reconstructed in Warp and visualized in IMOD.
Tomogram Segmentation and Visualization
Cellular membrane segmentation was performed using MemBrainSeg v0.0.9 with v10_alpha pretrained segmentation model with manual clean-up performed in UCSF ChimeraX56,57. Microtubule coordinates were traced manually in IMOD and visualized in ChimeraX. Ribosomes were identified through template matching in Warp using an ab initio ribosome model, with subsequent coordinate refinement in Relion. LRRK2 subtomograms were extracted from Dynamo coordinate tables58. The segmented membranes, microtubules, LRRK2 and ribosome coordinates were inputted in ChimeraX and ArtiaX for comprehensive visualization59.
Microtubule tracing and protofilament number analysis
Tomograms were visualized in IMOD to identify LRRK2IT decorated microtubules. These microtubules were manually traced in IMOD and coordinates were saved as contours. Filaments were cropped equidistantly along the traced microtubule coordinates using filament path function in Dynamo. The cropped volumes were averaged and smoothed to generate an ab initio reference for a microtubule. Based on this ab initio reference a cylindrical featureless reference was generated. Each microtubule was cropped, and alignment was run to symmetrize each individual microtubule. The results of the alignment were evaluated to assess the protofilament number of the microtubule.
LRRK2 subtomogram analysis
LRRK2 decorated microtubule coordinates were used to extract subtomograms of LRRK2 using Dynamo filament path function. 15 subunits at 8A separation were used to crop the LRRK2 subtomograms along the microtubule path. Cropped particles were averaged to generate initial reference. Spherical mask around central 4 subunits of LRRK2 was generated and initial refinements were run in Dynamo to align cropped particles. We removed duplicate particles from these subtomograms and exported them for import into Relion3. The subtomograms were refined in Relion3 with a C1 symmetry60. Focused refinements were performed on four central subunits of LRRK2IT. 3D classification was performed following initial refinements and classes with distinct LRRK2 decorations were selected for further refinements. The refinements were either centered on COR-COR interface of the protein or WD40-WD40 interaction interface. The center of the box was shifted accordingly, and refinements were performed with newly extracted particles.
Model Building
Structures of previously published LRRK2 were used as references for model building. ChimeraX rigid body fitting was used to fit each individual domain into the density maps. The fitting was performed stepwise. The domains of WD40and kinase C-lobe were separated and fit into the maps together. Following this, kinase N-lobe, COR and ROC domains were fitted into the density for the LRRK2. LRR-ANK domains were fit separately as well accounting for the remaining electron density in the in-situ maps. ARM domain was obtained through Alphafold2 model of LRRK2 and separated from the rest of the protein. ARM domain was individually mapped into the density for the N-terminal domains in electron density maps with larger lattice of LRRK2 accounting for the extra densities.
Data availability
The datasets generated in this study, including cryo-electron microscopy maps and associated structural models, have been prepared for deposition to the Electron Microscopy Public Image Archive (EMPIAR) and the Protein Data Bank (PDB). These data will be made publicly available in the respective repositories prior to publication.
Acknowledgements
Electron microscopy data were collected at the UCSD Cryo-Electron Microscopy Facility, which was built and equipped with funds from UCSD and an initial gift from the Agouron Institute. We thank the UCSD Physics Computing for computational support. This research was funded by Aligning Science Across Parkinson’s (grant number ASAP-000519) (A.E.L. and E.V.) through the Michael J. Fox Foundation for Parkinson’s Research (MJFF) (A.E.L. and E.V.), as well as from the National Science Foundation grant DBI 1920374 (to E.V.). T.B. was supported by an American Heart Association Predoctoral Fellowship. J.H. was supported by an EMBO long-term postdoctoral Fellowship (ALTF 871-2020). S.C. was supported by Jane Coffin Childs Memorial Fund for Medical Research. E.K. was supported by Jane Coffin Childs Memorial Fund for Medical Research. E.V. is a Howard Hughes Medical Institute
Additional information
Author contributions
Conceptualization: T.B., E.V.
Sample Preparation and Data Acquisition: T.B., J.H., S.C., E.K., A.F.
Analysis: T.B., J.H.
Writing: T.B., A.E.L., E.V. with input from all the authors.
Funding
Aligning Science Across Parkinson's (ASAP) (ASAP-000519)
Elizabeth Villa
Eva P Karasmanis
Tamar Basiashvili
Siyu Chen
W Alexander Flaherty
Andres E Leschziner
Joshua Hutchings
American Heart Association (AHA)
Tamar Basiashvili
European Molecular Biology Organization (EMBO) (ALTF 871-2020)
Joshua Hutchings
Jane Coffin Childs Memorial Fund for Medical Research (JCC)
Siyu Chen
Howard Hughes Medical Institute (HHMI)
Elizabeth Villa
Michael J. Fox Foundation for Parkinson's Research (MJFF)
W Alexander Flaherty
Joshua Hutchings
Andres E Leschziner
Elizabeth Villa
Siyu Chen
Eva P Karasmanis
Tamar Basiashvili
Additional files
References
- 1.Parkinson’s diseaseThe Lancet 386:896–912https://doi.org/10.1016/S0140-6736(14)61393-3PubMedGoogle Scholar
- 2.LRRK2 in Parkinson disease: challenges of clinical trialsNat. Rev. Neurol 16:97–107https://doi.org/10.1038/s41582-019-0301-2PubMedGoogle Scholar
- 3.Mutations in LRRK2 Cause Autosomal-Dominant Parkinsonism with Pleomorphic PathologyNeuron 44:601–607https://doi.org/10.1016/j.neuron.2004.11.005PubMedGoogle Scholar
- 4.Depletion of dopamine in Parkinson’s disease and relevant therapeutic options: A review of the literatureAIMS Neurosci 10:200–231https://doi.org/10.3934/neuroscience.2023017PubMedGoogle Scholar
- 5.Genome-wide association study reveals genetic risk underlying Parkinson’s diseaseNat. Genet 41:1308–1312https://doi.org/10.1038/ng.487PubMedGoogle Scholar
- 6.The PARK8 locus in autosomal dominant parkinsonism: confirmation of linkage and further delineation of the disease-containing intervalAm. J. Hum. Genet 74:11–19https://doi.org/10.1086/380647PubMedGoogle Scholar
- 7.A new locus for Parkinson’s disease (PARK8) maps to chromosome 12p11.2-q13.1Ann. Neurol 51:296–301https://doi.org/10.1002/ana.10113PubMedGoogle Scholar
- 8.Parkinson’s disease-associated mutations in leucinerich repeat kinase 2 augment kinase activityProc. Natl. Acad. Sci. U. S. A 102:16842–16847https://doi.org/10.1073/pnas.0507360102PubMedGoogle Scholar
- 9.Pathogenic Parkinson’s disease mutations across the functional domains of LRRK2 alter the autophagic/lysosomal response to starvationBiochem. Biophys. Res. Commun 441:862–866https://doi.org/10.1016/j.bbrc.2013.10.159PubMedGoogle Scholar
- 10.LRRK2 kinase in Parkinson’s diseaseScience 360:36–37https://doi.org/10.1126/science.aar5683PubMedGoogle Scholar
- 11.Systematic proteomic analysis of LRRK2-mediated Rab GTPase phosphorylation establishes a connection to ciliogenesiseLife 6:e31012https://doi.org/10.7554/eLife.31012PubMedGoogle Scholar
- 12.Leucine-Rich Repeat KinasesAnnu. Rev. Biochem 93:261–287https://doi.org/10.1146/annurev-biochem-030122-051144PubMedGoogle Scholar
- 13.Structural biology of Parkinson’s disease-associated leucine-rich repeat kinase 2 (LRRK2)J. Biol. Chem 301:110376https://doi.org/10.1016/j.jbc.2025.110376PubMedGoogle Scholar
- 14.LRRK2 Biology from structure to dysfunction: research progresses, but the themes remain the sameMol. Neurodegener 14:49https://doi.org/10.1186/s13024-019-0344-2PubMedGoogle Scholar
- 15.LRRK2 recruitment, activity, and function in organellesFebs J 289:6871–6890https://doi.org/10.1111/febs.16099PubMedGoogle Scholar
- 16.LRRK2 activation in idiopathic Parkinson’s diseaseSci. Transl. Med 10:eaar5429https://doi.org/10.1126/scitranslmed.aar5429PubMedGoogle Scholar
- 17.The Parkinson disease causing LRRK2 mutation I2020T is associated with increased kinase activityHum. Mol. Genet 15:223–232https://doi.org/10.1093/hmg/ddi439PubMedGoogle Scholar
- 18.Phosphoproteomics reveals that Parkinson’s disease kinase LRRK2 regulates a subset of Rab GTPaseseLife 5:e12813https://doi.org/10.7554/eLife.12813PubMedGoogle Scholar
- 19.Inhibition of Parkinson’s disease–related LRRK2 by type I and type II kinase inhibitors: Activity and structuresSci. Adv 9:eadk6191https://doi.org/10.1126/sciadv.adk6191PubMedGoogle Scholar
- 20.Cryo-electron tomography reveals the microtubule-bound form of inactive LRRK2eLife 13https://doi.org/10.7554/elife.100799PubMedGoogle Scholar
- 21.Pharmacology of LRRK2 with type I and II kinase inhibitors revealed by cryo-EMCell Discov 10:1–12https://doi.org/10.1038/s41421-023-00639-8PubMedGoogle Scholar
- 22.The In Situ Structure of Parkinson’s Disease-Linked LRRK2Cell 182:1508–1518https://doi.org/10.1016/j.cell.2020.08.004PubMedGoogle Scholar
- 23.Structure of LRRK2 in Parkinson’s disease and model for microtubule interactionNature 588:344–349https://doi.org/10.1038/s41586-020-2673-2PubMedGoogle Scholar
- 24.Structural basis for Parkinson’s disease-linked LRRK2’s binding to microtubulesNat. Struct. Mol. Biol 29:1196–1207https://doi.org/10.1038/s41594-022-00863-yPubMedGoogle Scholar
- 25.Structural analysis of the full-length human LRRK2Cell 184:3519–3527https://doi.org/10.1016/j.cell.2021.05.004PubMedGoogle Scholar
- 26.Rab29-dependent asymmetrical activation of leucinerich repeat kinase 2Science 382:1404–1411https://doi.org/10.1126/science.adi9926PubMedGoogle Scholar
- 27.14-3-3 binding maintains the Parkinson’s associated kinase LRRK2 in an inactive stateNat. Commun 16:7226https://doi.org/10.1038/s41467-025-62337-1PubMedGoogle Scholar
- 28.The Development and Design Strategy of Leucine-Rich Repeat Kinase 2 Inhibitors: Promising Therapeutic Agents for Parkinson’s DiseaseJ. Med. Chem 66:2282–2307https://doi.org/10.1021/acs.jmedchem.2c01552PubMedGoogle Scholar
- 29.Type II kinase inhibitors that target Parkinson’s disease–associated LRRK2Sci. Adv 11:eadt2050https://doi.org/10.1126/sciadv.adt2050PubMedGoogle Scholar
- 30.Small-molecule LRRK2 inhibitors for PD therapy: Current achievements and future perspectivesEur. J. Med. Chem 256:115475https://doi.org/10.1016/j.ejmech.2023.115475PubMedGoogle Scholar
- 31.Identification of GZD824 as an Orally Bioavailable Inhibitor That Targets Phosphorylated and Nonphosphorylated Breakpoint Cluster Region–Abelson (Bcr-Abl) Kinase and Overcomes Clinically Acquired Mutation-Induced Resistance against ImatinibJ. Med. Chem 56:879–894https://doi.org/10.1021/jm301581yPubMedGoogle Scholar
- 32.Characterization of a selective inhibitor of the Parkinson’s disease kinase LRRK2Nat. Chem. Biol 7:203–205https://doi.org/10.1038/nchembio.538PubMedGoogle Scholar
- 33.Impact of Type II LRRK2 inhibitors on signaling and mitophagyBiochem. J 478:3555–3573https://doi.org/10.1042/bcj20210375PubMedGoogle Scholar
- 34.Discovery of a 3-(4-Pyrimidinyl) Indazole (MLi-2), an Orally Available and Selective Leucine-Rich Repeat Kinase 2 (LRRK2) Inhibitor that Reduces Brain Kinase ActivityJ. Med. Chem 60:2983–2992https://doi.org/10.1021/acs.jmedchem.7b00045PubMedGoogle Scholar
- 35.Brain Penetrant LRRK2 InhibitorACS Med Chem Lett 3:658–662https://doi.org/10.1021/ml300123aPubMedGoogle Scholar
- 36.MLi-2, a Potent, Selective, and Centrally Active Compound for Exploring the Therapeutic Potential and Safety of LRRK2 Kinase InhibitionJ. Pharmacol. Exp. Ther 355:397–409https://doi.org/10.1124/jpet.115.227587PubMedGoogle Scholar
- 37.Discovery of highly potent, selective, and brainpenetrant aminopyrazole leucine-rich repeat kinase 2 (LRRK2) small molecule inhibitorsJ. Med. Chem 57:921–936https://doi.org/10.1021/jm401654jPubMedGoogle Scholar
- 38.LRRK2 Parkinson disease mutations enhance its microtubule associationHum. Mol. Genet 21:890–899https://doi.org/10.1093/hmg/ddr526PubMedGoogle Scholar
- 39.Conformation and dynamics of the kinase domain drive subcellular location and activation of LRRK2Proc. Natl. Acad. Sci 118:e2100844118https://doi.org/10.1073/pnas.2100844118PubMedGoogle Scholar
- 40.A designed ankyrin-repeat protein that targets Parkinson’s disease-associated LRRK2J. Biol. Chem 300:107469https://doi.org/10.1016/j.jbc.2024.107469PubMedGoogle Scholar
- 41.Practical Approaches for Cryo-FIB Milling and Applications for Cellular Cryo-Electron TomographyIn:
- Gonen T.
- Nannenga B. L.
- 42.Current data processing strategies for cryoelectron tomography and subtomogram averagingBiochem. J 478:1827–1845https://doi.org/10.1042/bcj20200715PubMedGoogle Scholar
- 43.Opening windows into the cell: focused-ion-beam milling for cryo-electron tomographyCurr. Opin. Struct. Biol 23:771–777https://doi.org/10.1016/j.sbi.2013.08.006PubMedGoogle Scholar
- 44.LRRK2 levels in immune cells are increased in Parkinson’s diseaseNPJ Park Dis 3:11https://doi.org/10.1038/s41531-017-0010-8PubMedGoogle Scholar
- 45.Role of LRRK2 in axonal transport and Parkinson’s diseaseBiochem. J 482:905–919https://doi.org/10.1042/bcj20253133PubMedGoogle Scholar
- 46.Micropatterning Transmission Electron Microscopy Grids to Direct Cell Positioning within Whole-Cell Cryo-Electron Tomography WorkflowsJ. Vis. Exp. JoVE https://doi.org/10.3791/62992PubMedGoogle Scholar
- 47.Grid patterning protocol for Cryo-FIB/ET workflowprotocols.io https://doi.org/10.17504/protocols.io.e6nvwdjewlmk/v3Google Scholar
- 48.Preparing samples from whole cells using focusedion-beam milling for cryo-electron tomographyNat. Protoc 15:2041–2070https://doi.org/10.1038/s41596-020-0320-xPubMedGoogle Scholar
- 49.Cryo-focused Ion Beam Sample Preparation for Imaging Vitreous Cells by Cryo-electron TomographyBio-Protoc 5:e1575–e1575https://doi.org/10.21769/bioprotoc.1575PubMedGoogle Scholar
- 50.Automated electron microscope tomography using robust prediction of specimen movementsJ Struct Biol 152:36–51https://doi.org/10.1016/j.jsb.2005.07.007PubMedGoogle Scholar
- 51.Parallel cryo electron tomography on in situ lamellaeNat. Methods 20:131–138https://doi.org/10.1038/s41592-022-01690-1PubMedGoogle Scholar
- 52.Implementation of a cryoelectron tomography tilt-scheme optimized for high resolution subtomogram averagingJ. Struct. Biol 197:191–198https://doi.org/10.1016/j.jsb.2016.06.007PubMedGoogle Scholar
- 53.Multiparticle cryo-EM refinement with M visualizes ribosome-antibiotic complex at 3.5 Å in cellsNat. Methods 18:186–193https://doi.org/10.1038/s41592-020-01054-7PubMedGoogle Scholar
- 54.Real-time cryo-electron microscopy data preprocessing with WarpNat. Methods 16:1146–1152https://doi.org/10.1038/s41592-019-0580-yPubMedGoogle Scholar
- 55.Computer Visualization of Three-Dimensional Image Data Using IMODJ. Struct. Biol 116:71–76https://doi.org/10.1006/jsbi.1996.0013PubMedGoogle Scholar
- 56.UCSF ChimeraX: Tools for structure building and analysisProtein Sci 32:e4792https://doi.org/10.1002/pro.4792PubMedGoogle Scholar
- 57.MemBrain v2: an end-to-end tool for the analysis of membranes in cryo-electron tomographybioRxiv https://doi.org/10.1101/2024.01.05.574336Google Scholar
- 58.Dynamo: A flexible, user-friendly development tool for subtomogram averaging of cryo-EM data in high-performance computing environmentsJ. Struct. Biol 178:139–151https://doi.org/10.1016/j.jsb.2011.12.017PubMedGoogle Scholar
- 59.ArtiaX: An electron tomography toolbox for the interactive handling of sub-tomograms in UCSF ChimeraXProtein Sci 31:e4472https://doi.org/10.1002/pro.4472PubMedGoogle Scholar
- 60.Resolving macromolecular structures from electron cryo-tomography data using subtomogram averaging in RELIONNat. Protoc 11:2054–2065https://doi.org/10.1038/nprot.2016.124PubMedGoogle Scholar
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.111075. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2026, Basiashvili et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 0
- downloads
- 0
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.