Discovery of LRRK2 inhibitors using sequential in silico joint pharmacophore space (JPS) and ensemble docking

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.

Computationally Designed Molecules Modulate ALS-Related Amyloidogenic TDP-43307-319 Aggregation

Abnormal cytosolic aggregation of TAR DNA-binding protein of 43 kDa (TDP-43) is observed in multiple diseases, including amyotrophic lateral sclerosis (ALS), frontotemporal lobar degeneration, and Alzheimer's disease. Previous studies have shown that TDP-43307-319 located at the C-terminal of TDP-43 can form higher-order oligomers and fibrils. Of particular interest are the hexamers that adopt a cylindrin structure that has been strongly correlated to neurotoxicity. In this study, we use the joint pharmacophore space (JPS) model to identify and generate potential TDP-43 inhibitors. Five JPS-designed molecules are evaluated using both experimental and computational methods: ion mobility mass spectrometry, thioflavin T fluorescence assay, circular dichroism spectroscopy, atomic force microscopy, and molecular dynamics simulations. We found that all five molecules can prevent the amyloid fibril formation of TDP-43307-319, but their efficacy varies significantly. Furthermore, among the five molecules, [AC0101] is the most efficient in preventing the formation of higher-order oligomers and dissociating preformed higher-order oligomers. Molecular dynamics simulations show that [AC0101] both is the most flexible and forms the most hydrogen bonds with the TDP-43307-319 monomer. The JPS-designed molecules can insert themselves between the β-strands in the hexameric cylindrin structure of TDP-43307-319 and can open its structure. Possible mechanisms for JPS-designed molecules to inhibit and dissociate TDP-43307-319 oligomers on an atomistic scale are proposed.

Computationally Designed Small Molecules Disassemble Both Soluble Oligomers and Protofibrils of Amyloid β-Protein Responsible for Alzheimer's Disease

Alzheimer's disease (AD) is one of the world's most pressing health crises. AD is an incurable disease affecting more than 6.5 million Americans, predominantly the elderly, and in its later stages, leads to memory loss, dementia, and death. Amyloid β (Aβ) protein aggregates have been one of the pathological hallmarks of AD since its initial characterization. The early stages of Aβ accumulation and aggregation involve the formation of oligomers, which are considered neurotoxic and play a key role in further aggregation into fibrils that eventually appear in the brain as amyloid plaques. We have recently shown by combining ion mobility mass spectrometry (IM-MS) and atomic force microscopy (AFM) that Aβ42 rapidly forms dodecamers (12-mers) as the terminal oligomeric state, and these dodecamers seed the early formation of Aβ42 protofibrils. The link between soluble oligomers and fibril formation is one of the essential aspects for understanding the root cause of the disease state and is critical to developing therapeutic interventions. Utilizing a joint pharmacophore space (JPS) method, potential drugs have been designed specifically for amyloid-related diseases. These small molecules were generated based on crucial chemical features necessary for target selectivity. In this paper, we utilize our combined IM-MS and AFM methods to investigate the impact of three second-generation JPS small-molecule inhibitors, AC0201, AC0202, and AC0203, on dodecamer as well as fibril formation in Aβ42. Our results indicate that AC0201 works well as an inhibitor and remodeler of both dodecamers and fibril formation, AC0203 behaves less efficiently, and AC0202 is ineffective.

Small-molecule LRRK2 inhibitors for PD therapy: Current achievements and future perspectives

Leucine-rich repeat kinase 2 (LRRK2) is a multifunctional protein that orchestrates a diverse array of cellular processes, including vesicle transport, autophagy, lysosome degradation, neurotransmission, and mitochondrial activity. Hyperactivation of LRRK2 triggers vesicle transport dysfunction, neuroinflammation, accumulation of α-synuclein, mitochondrial dysfunction, and the loss of cilia, ultimately leading to Parkinson's disease (PD). Therefore, targeting LRRK2 protein is a promising therapeutic strategy for PD. The clinical translation of LRRK2 inhibitors was historically impeded by issues surrounding tissue specificity. Recent studies have identified LRRK2 inhibitors that have no effect on peripheral tissues. Currently, there are four small-molecule LRRK2 inhibitors undergoing clinical trials. This review provides a summary of the structure and biological functions of LRRK2, along with an overview of the binding modes and structure-activity relationships (SARs) of small-molecule inhibitors targeting LRRK2. It offers valuable references for developing novel drugs targeting LRRK2.

Imaging Leucine-Rich Repeat Kinase 2 In Vivo with 18F-Labeled Positron Emission Tomography Ligand

Leucine-rich repeat kinase 2 (LRRK2) has been demonstrated to be closely involved in the pathogenesis of Parkinson's disease (PD), and pharmacological blockade of LRRK2 represents a new opportunity for therapeutical treatment of PD and other related neurodegenerative conditions. The development of an LRRK2-specific positron emission tomography (PET) ligand would enable a target occupancy study in vivo and greatly facilitate LRRK2 drug discovery and clinical translation as well as provide a molecular imaging tool for studying physiopathological changes in neurodegenerative diseases. In this work, we present the design and development of compound 8 (PF-06455943) as a promising PET radioligand through a PET-specific structure-activity relationship optimization, followed by comprehensive pharmacology and ADME/neuroPK characterization. Following an efficient 18F-labeling method, we have confirmed high brain penetration of [18F]8 in nonhuman primates (NHPs) and validated its specific binding in vitro by autoradiography in postmortem NHP brain tissues and in vivo by PET imaging studies.

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.

Discovery of G2019S-Selective Leucine Rich Repeat Protein Kinase 2 inhibitors with in vivo efficacy

Mutations in the Leucine Rich Repeat Protein Kinase 2 gene (LRRK2) are the most common genetic causes of Parkinson's Disease (PD). The G2019S mutation is the most common inherited LRRK2 mutation, occurs in the kinase domain, and results in increased kinase activity. We report the discovery and development of compound 38, an indazole-based, G2019S-selective (>2000-fold vs. WT) LRRK2 inhibitor capable of entering rodent brain (K_p = 0.5) and selectively inhibiting G2019S-LRRK2. The compounds disclosed herein present a starting point for further development of brain penetrant G2019S selective inhibitors that hopefully reduce lung phenotype side-effects and pave the way to providing a precision medicine for people with PD who carry the G2019S mutation.

LRRK2 Targeting Strategies as Potential Treatment of Parkinson's Disease

Parkinson's Disease (PD) affects millions of people worldwide with no cure to halt the progress of the disease. Leucine-rich repeat kinase 2 (LRRK2) is the most common genetic cause of PD and, as such, LRRK2 inhibitors are promising therapeutic agents. In the last decade, great progress in the LRRK2 field has been made. This review provides a comprehensive overview of the current state of the art, presenting recent developments and challenges in developing LRRK2 inhibitors, and discussing extensively the potential targeting strategies from the protein perspective. As currently there are three LRRK2-targeting agents in clinical trials, more developments are predicted in the upcoming years.

Glial Nrf2 signaling mediates the neuroprotection exerted by Gastrodia elata Blume in Lrrk2-G2019S Parkinson's disease

The most frequent missense mutations in familial Parkinson's disease (PD) occur in the highly conserved LRRK2/PARK8 gene with G2019S mutation. We previously established a fly model of PD carrying the LRRK2-G2019S mutation that exhibited the parkinsonism-like phenotypes. An herbal medicine, Gastrodia elata Blume (GE), has been reported to have neuroprotective effects in toxin-induced PD models. However, the underpinning molecular mechanisms of GE beneficiary to G2019S-induced PD remain unclear. Here, we show that these G2019S flies treated with water extracts of GE (WGE) and its bioactive compounds, gastrodin and 4-HBA, displayed locomotion improvement and dopaminergic neuron protection. WGE suppressed the accumulation and hyperactivation of G2019S proteins in dopaminergic neurons and activated the antioxidation and detoxification factor Nrf2 mostly in the astrocyte-like and ensheathing glia. Glial activation of Nrf2 antagonizes G2019S-induced Mad/Smad signaling. Moreover, we treated LRRK2-G2019S transgenic mice with WGE and found that the locomotion declines, the loss of dopaminergic neurons, and the number of hyperactive microglia were restored. WGE also suppressed the hyperactivation of G2019S proteins and regulated the Smad2/3 pathways in the mice brains. We conclude that WGE prevents locomotion defects and the neuronal loss induced by G2019S mutation via glial Nrf2/Mad signaling, unveiling a potential therapeutic avenue for PD.

Modulating ALS-Related Amyloidogenic TDP-43307-319 Oligomeric Aggregates with Computationally Derived Therapeutic Molecules

TDP-43 aggregates are a salient feature of amyotrophic lateral sclerosis (ALS), frontotemporal dementia (FTD), and a variety of other neurodegenerative diseases, including Alzheimer's disease (AD). With an anticipated growth in the most susceptible demographic, projections predict neurodegenerative diseases will potentially affect 15 million people in the United States by 2050. Currently, there are no cures for ALS, FTD, or AD. Previous studies of the amyloidogenic core of TDP-43 have demonstrated that oligomers greater than a trimer are associated with toxicity. Utilizing a joint pharmacophore space (JPS) method, potential drugs have been designed specifically for amyloid-related diseases. These molecules were generated on the basis of key chemical features necessary for blood-brain barrier permeability, low adverse side effects, and target selectivity. Combining ion-mobility mass spectrometry and atomic force microscopy with the JPS computational method allows us to more efficiently evaluate a potential drug's efficacy in disrupting the development of putative toxic species. Our results demonstrate the dissociation of higher-order oligomers in the presence of these novel JPS-generated inhibitors into smaller oligomer species. Additionally, drugs approved by the Food and Drug Administration for the treatment of ALS were also evaluated and demonstrated to maintain higher-order oligomeric assemblies. Possible mechanisms for the observed action of the JPS molecules are discussed.

Inhibiting and Remodeling Toxic Amyloid-Beta Oligomer Formation Using a Computationally Designed Drug Molecule That Targets Alzheimer's Disease

Alzheimer's disease (AD) is rapidly reaching epidemic status among a burgeoning aging population. Much evidence suggests the toxicity of this amyloid disease is most influenced by the formation of soluble oligomeric forms of amyloid β-protein, particularly the 42-residue alloform (Aβ42). Developing potential therapeutics in a directed, streamlined approach to treating this disease is necessary. Here we utilize the joint pharmacophore space (JPS) model to design a new molecule [AC0107] incorporating structural characteristics of known Aβ inhibitors, blood-brain barrier permeability, and limited toxicity. To test the molecule's efficacy experimentally, we employed ion mobility mass spectrometry (IM-MS) to discover [AC0107] inhibits the formation of the toxic Aβ42 dodecamer at both high (1:10) and equimolar concentrations of inhibitor. Atomic force microscopy (AFM) experiments reveal that [AC0107] prevents further aggregation of Aβ42, destabilizes preformed fibrils, and reverses Aβ42 aggregation. This trend continues for long-term interaction times of 2 days until only small aggregates remain with virtually no fibrils or higher order oligomers surviving. Pairing JPS with IM-MS and AFM presents a powerful and effective first step for AD drug development. Graphical Abstract.

Molecular docking, 3D-QSAR and structural optimization on imidazo-pyridine derivatives dually targeting AT1 and PPARg

Telmisartan, a bifunctional agent of blood pressure lowering and glycemia reduction, was previously reported to antagonize angiotensin II type 1 (AT1) receptor and partially activate peroxisome proliferator-activated receptor γ (PPARγ) simultaneously. Through the modification to telmisartan, researchers designed and obtained imidazo-\pyridine derivatives with the IC₅₀s of 0.49∼94.1 nM against AT1 and EC₅₀s of 20∼3640 nM towards PPARγ partial activation. For minutely inquiring the interaction modes with the relevant receptor and analyzing the structure-activity relationships, molecular docking and 3D-QSAR (Quantitative structure-activity relationships) analysis of these imidazo-\pyridines on dual targets were conducted in this work. Docking approaches of these derivatives with both receptors provided explicit interaction behaviors and excellent matching degree with the binding pockets. The best CoMFA (Comparative Molecular Field Analysis) models exhibited predictive results of q²=0.553, r²=0.954, SEE=0.127, r²_pred=0.779 for AT1 and q²=0.503, r²=1.00, SEE=0.019, r²_pred=0.604 for PPARγ, respectively. The contour maps from the optimal model showed detailed information of structural features (steric and electrostatic fields) towards the biological activity. Combining the bioisosterism with the valuable information from above studies, we designed six molecules with better predicted activities towards AT1 and PPARγ partial activation. Overall, these results could be useful for designing potential dual AT1 antagonists and partial PPARγ agonists.

Design of Leucine-Rich Repeat Kinase 2 (LRRK2) Inhibitors Using a Crystallographic Surrogate Derived from Checkpoint Kinase 1 (CHK1)

Mutations in leucine-rich repeat kinase 2 (LRRK2), such as G2019S, are associated with an increased risk of developing Parkinson's disease. Surrogates for the LRRK2 kinase domain based on checkpoint kinase 1 (CHK1) mutants were designed, expressed in insect cells infected with baculovirus, purified, and crystallized. X-ray structures of the surrogates complexed with known LRRK2 inhibitors rationalized compound potency and selectivity. The CHK1 10-point mutant was preferred, following assessment of surrogate binding affinity with LRRK2 inhibitors. Fragment hit-derived arylpyrrolo[2,3-b]pyridine LRRK2 inhibitors underwent structure-guided optimization using this crystallographic surrogate. LRRK2-pSer935 HEK293 IC₅₀ data for 22 were consistent with binding to Ala2016 in LRRK2 (equivalent to Ala147 in CHK1 10-point mutant structure). Compound 22 was shown to be potent, moderately selective, orally available, and brain-penetrant in wild-type mice, and confirmation of target engagement was demonstrated, with LRRK2-pSer935 IC₅₀ values for 22 in mouse brain and kidney being 1.3 and 5 nM, respectively.

Recent Advancements in Docking Methodologies

Nowadays molecular docking has become an important methodology in CADD (Computer-Aided Drug Design)-assisted drug discovery process. It is an important computational tool widely used to predict binding mode, binding affinity and binding free energy of a protein-ligand complex. The important factors responsible for accurate results in docking studies are correct binding site prediction, use of suitable small-molecule databases, consistent docking pose, high dock score with good MD (Molecular Dynamics), clarity whether the compound is an inhibitor or agonist, etc. However, still there are several limitations which make it difficult to obtain accurate results from docking studies. In this chapter, the main focus is on recent advancements in various aspects of molecular docking such as ligand sampling, protein flexibility, scoring functions, fragment docking, post-processing, docking into homology models and protein-protein docking.

Screening for chemical modulators for LRRK2

After the discovery of leucine-rich repeat kinase 2 (LRRK2) as a risk factor for sporadic Parkinson's disease (PD) and mutations in LRRK2 as a cause of some forms of familial PD, there has been substantial interest in finding chemical modulators of LRRK2 function. Most of the pathogenic mutations in LRRK2 are within the enzymatic cores of the protein; therefore, many screens have focused on finding chemical modulators of this enzymatic activity. There are alternative screening approaches that could be taken to investigate compounds that modulate LRRK2 cellular functions. These screens are more often phenotypic screens. The preparation for a screen has to be rigorous and enable high-throughput accurate assessment of a compound's activity. The pipeline to beginning a drug screen and some LRRK2 inhibitor and phenotypic screens will be discussed.

LRRK2 inhibitors and their potential in the treatment of Parkinson's disease: current perspectives

Major advances in understanding how genetics underlies Parkinson's disease (PD) have provided new opportunities for understanding disease pathogenesis and potential new targets for therapeutic intervention. One such target is leucine-rich repeat kinase 2 (LRRK2), an enigmatic enzyme implicated in both familial and idiopathic PD risk. Both academia and industry have promoted the development of potent and selective inhibitors of LRRK2, and these are currently being employed to assess the safety and efficacy of such compounds in preclinical models of PD. This review examines the evidence that LRRK2 kinase activity contributes to the pathogenesis of PD and outlines recent progress on inhibitor development and early results from preclinical safety and efficacy testing. This review also looks at some of the challenges remaining for translation of LRRK2 inhibitors to the clinic, if indeed this is ultimately warranted. As a disease with no current cure that is increasing in prevalence in line with an aging population, there is much need for developing new treatments for PD, and targeting LRRK2 is currently a promising option.

A multicomponent pharmacophore fragment-decoration approach to identify selective LRRK2-targeting probes

Herein we report the development of a new versatile chemical tool for the rapid identification of LRRK2-targeting probes as potential anti-Parkinson's agents.