Research article | Open | Published:
The PINK1 p.I368N mutation affects protein stability and ubiquitin kinase activity
Molecular Neurodegenerationvolume 12, Article number: 32 (2017)
Mutations in PINK1 and PARKIN are the most common causes of recessive early-onset Parkinson’s disease (EOPD). Together, the mitochondrial ubiquitin (Ub) kinase PINK1 and the cytosolic E3 Ub ligase PARKIN direct a complex regulated, sequential mitochondrial quality control. Thereby, damaged mitochondria are identified and targeted to degradation in order to prevent their accumulation and eventually cell death. Homozygous or compound heterozygous loss of either gene function disrupts this protective pathway, though at different steps and by distinct mechanisms. While structure and function of PARKIN variants have been well studied, PINK1 mutations remain poorly characterized, in particular under endogenous conditions. A better understanding of the exact molecular pathogenic mechanisms underlying the pathogenicity is crucial for rational drug design in the future.
Here, we characterized the pathogenicity of the PINK1 p.I368N mutation on the clinical and genetic as well as on the structural and functional level in patients’ fibroblasts and in cell-based, biochemical assays.
Under endogenous conditions, PINK1 p.I368N is expressed, imported, and N-terminally processed in healthy mitochondria similar to PINK1 wild type (WT). Upon mitochondrial damage, however, full-length PINK1 p.I368N is not sufficiently stabilized on the outer mitochondrial membrane (OMM) resulting in loss of mitochondrial quality control. We found that binding of PINK1 p.I368N to the co-chaperone complex HSP90/CDC37 is reduced and stress-induced interaction with TOM40 of the mitochondrial protein import machinery is abolished. Analysis of a structural PINK1 p.I368N model additionally suggested impairments of Ub kinase activity as the ATP-binding pocket was found deformed and the substrate Ub was slightly misaligned within the active site of the kinase. Functional assays confirmed the lack of Ub kinase activity.
Here we demonstrated that mutant PINK1 p.I368N can not be stabilized on the OMM upon mitochondrial stress and due to conformational changes in the active site does not exert kinase activity towards Ub. In patients’ fibroblasts, biochemical assays and by structural analyses, we unraveled two pathomechanisms that lead to loss of function upon mutation of p.I368N and highlight potential strategies for future drug development.
Parkinson's disease (PD) is the second most common neurodegenerative disease and is characterized by mitochondrial dysfunctions and loss of dopaminergic neurons [1, 2]. While most cases of PD are sporadic, mutations in the PTEN-induced putative kinase 1 (PINK1) and PARKIN (PARK2) are the most common causes of early-onset forms of PD [3, 4]. It is known that the mitochondrial kinase PINK1 and the cytosolic E3 ubiquitin (Ub) ligase PARKIN functionally cooperate during mitochondrial quality control to identify, label and remove damaged organelles .
In healthy mitochondria, newly translated full-length PINK1 (~63 kDa) is constitutively imported through the TOM complex, the translocase of the outer mitochondrial membrane (OMM) and then through the TIM complex, the translocase of the inner mitochondrial membrane (IMM). Upon import, the N-terminal mitochondrial targeting sequence (MTS) of PINK1 is cleaved off by the matrix processing peptidase (MPP) and the intermediate isoform (~60 kDa) is further processed by the Presenilin-associated rhomboid-like protein (PARL) in the IMM to generate a 52 kDa PINK1 fragment. This cleaved form of PINK1 is exported back to the cytosol by an unknown mechanism and degraded by the Ub/proteasome system (UPS) [6–9].
Upon mitochondrial depolarization, full-length PINK1 is no longer imported into mitochondria, but accumulates, with the kinase domain facing the cytosol, on the OMM and forms a dimeric structure associated with the TOM complex [10–12]. This allows phosphorylation of its substrates, the small modifier protein Ub and its E3 ligase PARKIN, at a conserved Ser65 residue [13–15]. Both PINK1-mediated phosphorylation events fully activate the auto-inhibited enzymatic functions of PARKIN and further facilitate its recruitment from the cytosol . Then, PINK1 and PARKIN together cooperatively label damaged mitochondria in a feed forward mechanism with phosphorylated poly-Ub chains that serve as the “mitophagy tag” for their removal via the autophagy/lysosome system [17–21]. Interference at any of these steps abrogates protection through PINK1/PARKIN-directed mitochondrial quality control.
A deeper understanding of the particular pathomechanisms of individual PARKIN mutations has further allowed delineating the regulation and sequence of events along this process . These detailed structure-function studies have already spurred efforts for a rationalized drug design to activate PARKIN through different strategies . Similarly, missense mutations in PINK1 can interfere with mitochondrial quality control through different molecular mechanisms . For instance PINK1 p.Q456X results in instability of its transcript through non-sense mediated decay, leading to a complete loss-of-function on the protein level . In contrast, the mutant p.G411S is expressed and upon damage forms a dimer on the OMM similar to PINK1 wild type (WT) . However, in addition to partial loss of its kinase activity, in a heterodimer with PINK1 WT, p.G411S also exerts a partial, dominant-negative effect. This lowers overall protection through mitochondrial quality control and enhances risk for PD even in heterozygous carriers.
Towards a complete dissection of PINK1 regulation and a rationalized drug design in the future, we characterized here the novel p.I368N mutation in a combined clinical-genetic and structure-function approach. In patients’ fibroblasts, we found that stabilization of PINK1 p.I368N on the OMM is greatly impaired upon mitochondrial stress. Additional structural analysis revealed changes in the active site of PINK1 p.I368N that should affect ATP and substrate binding. Indeed, in cells we observed loss of kinase activity even under conditions of forced expression. Altogether, our study identified two distinct molecular mechanisms that can drive pathogenicity of PINK1 missense mutations and might help guiding future drug design.
The Polish family with the PINK1 p.I368N mutation was recruited to participate in this research study from the movement disorders clinic of the Department of Neurology of the Medical University of Warsaw in Warsaw, Poland. The diagnosis of PD was made based on the UK Parkinson’s Disease Society Brain Bank clinical diagnostic criteria . At the time of blood and skin specimens retrieval, the research subjects were independently clinically examined by two board certified neurologists (DK and ZKW) utilizing standardized narrative Mayo Clinic (MC) history, standardized MC narrative family history, and standardized MC neurological examination . In addition the Hoehn and Yahr and Mini Mental State examination scales were obtained [29, 30]. The blood samples and 3 mm single skin specimens were retrieved in routine fashion by utilizing the standard universal precautions. Both subjects were right handed and parkinsonian features were more expressed on the dominant side. The response to the standard antiparkinsonian therapy including levodopa was excellent and there was no need to proceed with surgical interventions such as for example a deep brain stimulation.
Material for analysis of genomic DNA was extracted from blood and fibroblasts using a standard protocol. To confirm the genotype (c.T1103A/c.T1103A; p.I368N/p.I368N) and exclude other mutations in PINK1, all coding exons were directly sequenced by Sanger method and analyzed against the reference sequence NM_032409.2 (NCBI36/hg18). As far as we know this is the only family carrying this specific mutation in the Polish population. Though the analyzed siblings both carry a homozygous mutation, there was no consanguineous marriage of their parents. Gene rearrangements were analyzed by multiple ligation-dependent probe amplification (MLPA) technique using available tests SALSA P051 PD probes (MRC Holland) – probes for all PINK1 exons. Mutations in PINK1 and PARKIN were excluded in control fibroblast cells by sequencing of all coding exons.
Modeling PINK1 structures and refinement
The modeling of the full-length human PINK1 protein, NP_115785.1 (581 amino acid residues), has been described recently . In brief, each individual domain was modeled as a separate unit and built into a composite full-length structure. The hybrid model is derived from consensus between the programs PRIME (Prime v3.0, Schrödinger, LLC, New York, NY) [31, 32], YASARA SSP/Homology/PSSM Method [31–38], DISTILL (Porter) and TASSER [39–44] and combines homology, threading, ab initio, and compositing techniques.
Modification and mutation modeling
Mutations of amino acids were completed using the Maestro within the Schrödinger suite and parameterized using OPLS2.0 and Amber FF. Additionally, Maestro was used for placing mutated residues (or extending) automatically within an existing peptide chain. Also, MacroModel features within Maestro allow for the quick minimization of the structure for local geometry fixes to correct stereochemistry and packing of the amino acids. Further minimization was completed on the post-Schrödinger model using OPLS2005 within YASARA2, which has an AMBER set and can be used to further parameterize modifications to import into existing molecular dynamics integration engines, as the parameters for the modification are well documented for YASARA and Schrödinger [32, 37, 45–47].
Prediction of binding energy changes in proteins using FoldX algorithm
For each proposed mutation locus, we equilibrated the corresponding flexibility zone in the wild type complex and evaluated the free energy of unfolding (ΔGwt) with the KB potential FoldX . Combined with our own codes, we have a method to scan across trajectories and examine interaction zones of interest. We repeated the calculation for the mutant (p.I368N) complex to obtain the respective free energy of unfolding (ΔGmutant). The difference between these two quantities is our estimate of the experimental change in binding free energy induced by the mutation : ΔΔG = ΔGwt- ΔGmutant.
Molecular dynamics simulation and methods
MDS was completed on each model for conformational sampling, using methods previously described [50–53]. Each PINK1 system was minimized with relaxed restraints using Steepest Descent and Conjugate Gradient PR, and equilibrated in solvent with physiological salt conditions, as described in the literature [50–53]. After equilibration was established, each system was allowed to run an additional MD production length of >250 nanoseconds. The primary purpose of MD for this study was conformational stability, refinement, and interaction calculations that may occur at the active site or dimer interface. The protocol for refinement include the following steps: (1) Minimization with explicit water molecules and ions, (2) Energy minimization of the entire system, and (3) MDS for >2 ns to relax to the force field (AMBER03). Following the refinement protocol, production simulations were completed to collect data.
Molecular dynamics protocol
AMBER03 force fields were used with the current release of NAnoscale Molecular Dynamics 2 engine [54, 55]. The PINK1 systems simulated, including hydrogens, consist of 1.74 × 104 atoms prior to solvation with TIP3P water and ions. In all cases, we neutralized with counter-ions, and then created a solvent with 150 mM Na + Cl- to recreate physiological strength. TIP3P water molecules were added around the protein at a depth of 12–15 Å from the edge of the molecule depending upon the side . Our protocol has been previously described in the literature [50, 51, 53, 57]. Solvated protein simulations consist of a box with 2.23–3.35 × 105 atoms including nucleic acids, proteins, counter-ions, solvent ions, and solvent waters. Simulations were carried out using the particle mesh Ewald technique with repeating boundary conditions with a 9 Å nonbonded cut-off, using SHAKE with a 2-fs timestep. Pre-equilibration was started with three stages of minimization with 10,000 steps of SD, PRCG, relaxing restraints, then followed by 1000 ps of heating under MD, with the atomic positions of nucleic and protein fixed. Then, two cycles of minimization (5000 steps each) and heating (1000–5000 ps) were carried out with soft restraints of 10 and 5 kcal/(mol · Å2) applied to all nucleic and protein backbone atoms. Next, 5000-steps of minimization were performed with solute restraints reduced to 1 kcal/(mol · Å2). Following that, 2000–4000 ps of MDS were completed using relaxing restraints (1 kcal/(mol · Å2)) until all atoms are unrestrained, while the system was slowly heated from 1 to 310 K using velocity rescaling upon reaching the desired 310 K during this equilibration phase. Additionally, NPT equilibration based MD was used with velocity rescaling for 1–2 ns. Finally, production runs of MD were carried out with constant pressure boundary conditions (relaxation time of 1.0 ps). A constant temperature of 310 K was maintained using the Berendsen weak-coupling algorithm with a time constant of 1.0 ps. SHAKE constraints were applied to all hydrogens to eliminate X-H vibrations, which yielded a longer simulation time step (2 fs). Our methods for equilibration and production run protocols are in the literature [50, 51, 53, 57]. Equilibration was determined from a global flattening of RMSD over time after an interval of 1–3 ns or similarly with total energy. Translational and rotational center-of-mass motions were initially removed. Periodically, simulations were interrupted to have the center-of-mass removed again by a subtraction of velocities to account for the “flying ice-cube” effect . Following the simulation, the individual frames were superposed back to the origin, to remove rotation and translation effects. Archived snapshots were recorded every 2 ps, yielding >2000 snapshots per simulation.
Docking of kinase PINK1 and substrate ubiquitin
Using mixture of Piper protein-protein docking for generating 40,000 conformers of random seed orientation over the PINK1 and ubiquitin proteins, we applied restriction parameters as follows: only kinase domains from PINK1 (N-lobe and C-lobe) and Ser65 face of Ub were weighted with increased attraction parameters +0.25, whilst we also did multiple LC-MOD/MC generated side chain re-orientations during the final combination with an simulated annealing protocol to allow for settling of interactions before applying a PRCG minimization protocol. We examined the docking interfaces with brief bursts of simulated annealing and MDS for further refinement for taking the top seeded structures and elimination of the poor scorers. Top surviving poses were scored and best from each series was retained [50, 52, 53, 59, 60]. The best pairs generated were refined and studied using MDS calculations.
Primary human dermal fibroblasts were grown in Dulbecco’s modified Eagle medium (DMEM, Invitrogen) supplemented with 10% fetal bovine serum (FBS, BioWest), 1% PenStrep and 1% non-essential amino acids (both Invitrogen). Control fibroblasts were received from Cell Applications Inc., mutations in PINK1 and PARKIN genes were excluded by sequencing and PINK1-dependent responses in those cells have been characterized before [18, 25, 26]. HeLa cells (ATCC) were cultured in DMEM supplemented with 10% FBS. Stable HeLa cell clones expressing EGFP-Parkin or 3xFLAG-Parkin C431S have been described before . All cells were maintained at 37 °C, 5% CO2 in a humidified atmosphere.
Chemical treatments, siRNA and DNA transfection of cells
Carbonyl cyanide m-chlorophenyl hydrazone (CCCP), epoxomicin and cycloheximide were purchased from Sigma-Aldrich, valinomycin from Axxora. siRNA transfections were performed with 20nM control (all stars negative control) or PINK1-specific siRNA (5′ GACGCTGTTCCTCGTTATGAA-3′) using HiPerfect (all from Qiagen) or Lipofectamine2000 (Invitrogen). siRNA resistant PINK1-V5 constructs have been described before [5, 24]. The mutations p.F104A, p.L347P and p.I368N were introduced by site-directed mutagenesis and verified by Sanger sequencing. The standard transfection protocol for DNA was as follows: for one well of a 12-well plate 1 μg of DNA and 2.5 μl of Lipofectamine2000 (Invitrogen) were each mixed with 100 μl of Opti-MEM media (Invitrogen), incubated for 5 min at room temperature, mixed together and incubated 20 min before adding to the cells. For low-level expression DNA and lipofectamine amounts were reduced by 50–75% (total DNA per 12 well: 0.5 or 0.25 μg).
RNA extraction and real-time PCR analysis
Total mRNA from fibroblasts was prepared using the RNeasy mini kit (Qiagen) according to the manufacturer’s protocol and then reverse-transcribed into cDNA using a High Fidelity cDNA kit (Roche Applied Science). cDNA was amplified by real-time PCR on a LightCycler 480 system (Roche) using Taqman probes (BioRad PINK1 Id qHsaCEP0053094 and RPL27 Id qHsaCEP0051648). Levels were determined using the second derivative method.
The following primary antibodies were used for Western blot (WB) and/or immunofluorescence (IF): rabbit anti-CDC37 (WB, 1:1,000, 4793, Cell Signaling), mouse anti-HSP90 (WB, 1:2,000, 610419, BD), mouse anti-FLAG (WB, 1:150,000, F3165, Sigma Aldrich), mouse anti-GAPDH (WB, 1:150,000, H86504M, Meridian Life science), mouse anti-Mitofusin 1 (WB, 1:5,000, ab57602, Abcam), rabbit anti-p38 MAPK (WB, 1:1,000, #9212, Cell Signaling), rabbit anti-PINK1 (WB, 1:2,000, BC100-494, Novus Biologicals), rabbit anti-PINK1 (WB, 1:2,000, IF, 1:100, #6946, Cell Signaling), mouse anti-Parkin (WB, 1:3,000, #4211, Prk8, Cell Signaling), mouse anti-TOM20 (IF, 1:250, sc-17764, Santa Cruz Biotechnology), rabbit anti-TOM20 (WB, 1:20,000, 11802-1-AP, ProteinTech Group), rabbit anti-TOM40 (WB, 1:1,000, 18409-1-AP, ProteinTech Group), mouse anti-Ub (WB, 1:1,000, MAB1510, ubi-1, Millipore), rabbit anti-p-Ser65-Ub (WB, 1:5,000-1:15,000, IF, 1:250, in house ), mouse anti-V5 (WB, 1:5,000, R960-25, Invitrogen), rabbit anti-V5 (WB, 1:5,000, ab9115, Abcam), mouse anti-vinculin (WB, 1:100,000-500,000, V9131, Sigma). Streptavidin HRP (WB, 1:150,0000, Pierce).
Protein extraction, SDS-PAGE and Western blot (WB)
Cells were collected and lysed in RIPA buffer (50 mM Tris pH 8.0, 150 mM NaCl, 1% NP-40, 0.5% deoxycholate, 0.1% SDS) plus Complete protease and PhosSTOP phosphatase inhibitors (Roche Applied Science) on ice for 30 min. Proteins were collected as supernatants after centrifugation at 20,000 g for 15 min at 4 °C. SDS–PAGE was performed using 8–16% or 16% Tris Glycine gels (Invitrogen). Proteins were transferred onto PVDF membranes and detected using standard immunoblotting procedures. To detect phosphorylated PINK1, 8% Tris Glycine gels containing 50 μM Phos-tag acrylamide (Wako chemicals) and 100 μM ZnCl2 were used . After electrophoresis, Phos-tag acrylamide gels were washed using transfer buffer with 0.01% SDS and 1 mM EDTA for 20 min and then with transfer buffer containing 0.01% SDS without EDTA for 10 min to remove excess Zn2+ before transfer.
Subcellular fractionation and solubility analysis
Cells were harvested in fractionation buffer containing 20 mM HEPES-KOH (pH 7.4), 1 mM EDTA, 250 mM sucrose and cocktail of protease and phosphatase inhibitors (Roche Applied Science). After homogenization with a 22 G syringe (20 strokes), nuclei and cell debris were removed by centrifugation at 1500 g for 5 min. The post-nuclear supernatant (PNS) was transferred into a new tube and centrifuged at 12,000 g for 10 min. The mitochondria-enriched pellet (‘mitochondrial fraction’) was dissolved in RIPA buffer. Cytoplasmic fractions were isolated as supernatant after ultracentrifugation at 100,000 g for 30 min. For solubility studies, the PNS was dissolved in RIPA buffer and pelleted at 14,000 g for 15 min. The supernatant was used as the “soluble” fraction, the pellet was further washed and solubilized in 8 M urea and used as the “insoluble” fraction.
Analysis of ubiquitin charging of PARKIN C431S
Transfected HeLa 3xFLAG-Parkin C431S cells  were harvested in preheated (95 °C) SDS lysis buffer (50 mM Tris pH 7.6, 150 mM NaCl, 1% SDS). Lysates were homogenized by 10 strokes through a 20G needle. To verify the band shift by oxyester formation, aliquots of lysates were treated with or without NaOH (final concentration: 100 mM) for 1 h at 37 °C before loading onto SDS gels.
Fibroblasts were seeded onto PDL (P0899, Sigma-Aldrich) coated glass coverslips treated with valinomycin, CCCP or left untreated. Cells were then fixed with 4% (w/v) paraformaldehyde, permeabilized with 1% Triton-X-100, blocked with 10% goat serum and incubated with primary antibodies followed by secondary antibodies (anti-rabbit Alexa Fluor-488, anti-mouse IgG Alexa Fluor-568, 1:1000, Invitrogen). Nuclei were stained with Hoechst 33342 (1:5000, Invitrogen). For staining of PINK1 tyramide signal amplification (T20922, Invitrogen) was used. Coverslips were mounted onto microscope slides using fluorescence mounting medium (Dako). High resolution confocal images were taken with an AxioObserver microscope equipped with an ApoTome Imaging System (Zeiss).
High content imaging (HCI)
To quantify p-Ser65-Ub levels in valinomycin treated fibroblasts and EGFP-PARKIN translocation in CCCP treated HeLa cells, we used the BD Pathway 855 HCI platform (BD Biosciences) as described earlier [18, 26, 61]. To collect signal only from cells transfected with siRNA and various PINK1-V5 constructs, co-expressed mCherry fluorescent intensity was used as a restriction parameter for the analysis of PARKIN translocation.
HeLa cells were transiently transfected with V5-6xHis-tagged PINK1 WT, p.I368N or p.L347P for 48 h and were then treated with 15 μM CCCP for 2 h or left untreated. Cells were lysed in co-IP buffer (50 mM HEPES pH 7.5, 150 mM NaCl, 10 mM KCl, 1.5 mM MgCl2, 1 mM EDTA, 0.5 mM EGTA, 10% glycerol and 0.2% NP-40) supplemented with protease and phosphatase inhibitors (Roche Applied Science) and briefly sonicated. For PINK1-V5 immunoprecipitation (IP), total cell lysate was incubated overnight at 4 °C with mouse anti-V5-agarose beads (A7345, Sigma-Aldrich). Formed immuno-complexes were spun down, washed 3x with co-IP buffer and eluted from beads with 50 μl 2x SDS sample buffer and boiled at 95 °C. 10 μg of input cell lysates and 10 μl of immunoprecipitates were analyzed by Western blot.
In vitro kinase assay
Immunoprecipitated WT and p.I368N PINK1-V5 from CCCP treated and untreated HeLa cells were washed twice with co-IP buffer and once with in vitro phosphorylation buffer (20 mM HEPES pH 7.4, 50 mM NaCl, 5 mM MgCl2, 0.1 mM EGTA, 0.01% Triton X-100). After complete removal of the buffer 100 μl of the in vitro reaction mixture was added. The reaction contained 1 μg of N-terminally biotinylated Ub (Boston Biochem), 1 mM TCEP (Gold Bio) and 2 mM ATP trisodium salt (Sigma-Aldrich) in phosphorylation buffer. Reactions were carried out at 37 °C under constant agitation for 24 h, stopped by addition of 6x SDS sample buffer, heated to 56 °C for 15 min and analyzed by immunoblotting using p-Ser65-Ub antibody. Detection of biotinylated Ub with streptavidin-HRP (Pierce) served as loading control.
Densitometry, PINK1 half-life calculation and statistical methods
Densitometric analysis of Western blots was performed with the Image Studio software (LiCor). The half-life of PINK1 proteins was calculated by non-linear one phase decay curve fitting using GraphPad Prism version 6.
All quantitative results are expressed as mean ± SEM from at least three independent experiments. Statistical comparisons were performed using one-way and two-way ANOVA, respectively, with Tukey’s post-hoc test (*, p < 0.05; **, p < 0.005; ***, p < 0.0005). All analyses were performed using GraphPad Prism version 6.
Clinical features of PINK1 p.I368N are similar to sporadic late-onset PD
Here we analyzed two patients that are homozygous carriers of a recently identified PINK1 c.T1103A (p.I368N) mutation . Detailed clinical data of the affected members are provided in Additional file 1: Table S1. Both patients have early-onset and slowly progressive PD that started at relatively young age (28 and 33 years, respectively), consistent with the typical clinical presentation of homozygous PINK1 loss-of-function . Some clinical features (e.g. laterality of parkinsonism, resting tremor, bradykinesia and good response to levodopa) are indistinguishable from sporadic late-onset PD. Yet, REM sleep behavior disorder (RBD), which is a clinical biomarker of synucleinopathies, was not detected in these patients . Both patients show mild autonomic failure (e.g. constipation), but no dementia, and no depression or anxiety, which are occasionally recognized as a non-motor symptom of PINK1-associated PD .
Structural changes of the p.I368N mutation distort the ATP-binding site of PINK1
The residue Ile368 is located at the +1 position of the catalytic loop of the consensus motif HRDLKxxN in Ser/Thr kinases (Fig. 1a). Using homology modeling, we have recently generated a structural model of full-length PINK1 at an all-atom resolution . To examine how the p.I368N mutation would affect protein structure we performed molecular modeling and free, unbiased molecular dynamics simulations (MDS) (Additional file 2: Movie S1 and Additional file 3: Movie S2 for PINK1 WT and p.I368N mutant, respectively). The calculated Root Mean Square Fluctuations (RMSF) provide the movement of each individual residue of PINK1 (581 amino acids) as a time-averaged sum over the course of the simulation. Though global conformational changes were not observed, overall RMSF of PINK1 p.I368N were slightly lower than WT, yet with few higher peaks around residues of the ATP-binding pocket (Fig. 1b). Overlay of the WT and mutant kinase domains together with Root Mean Square Deviations (RMSD), as a measure of structural change, revealed an immediate but minor difference before and after MDS (Fig. 1c, d). However, specifically the orientation of helices surrounding the ATP-binding pocket was altered in the mutant compared to PINK1 WT (Fig. 1e) and resulted in local changes of RMSD within the ATP-binding site (Fig. 1f). In addition, surface analyses revealed a distorted ATP pocket shape for p.I368N (Fig. 1g). Together with FoldX measurements of Gibb’s free energy that indicated decreased local stability (Additional file 4: Table S2) and an increase in the radius of gyration (Fig. 1h) this indicated a negative effect on the catalytic activity of PINK1.
PINK1 p.I368N decreases binding efficiency and mispositions the substrate Ub
To study the enzymatic activity, we had recently docked the substrate Ub into the kinase domain of PINK1 WT . We could show that auto-phosphorylation of PINK1 at Ser228 and Ser402 increased interactions with Ub that not only stabilized the kinase/substrate complex but also positively affected atomic distances for a more efficient phosphate transfer. Here, we docked Ub to a p.I368N homodimer to investigate binding and alignment of the substrate within the active site of the kinase (Fig. 2a). Indeed, closer inspection of the mutant enzymatic complex, revealed a slight mispositioning of the Ub substrate (Fig. 2b). However, distances measured between key atoms of PINK1, ATP (terminal phosphate), and Ub (oxygen of Ser65) were only slightly altered, and still within the range for effective catalysis. In summary, structural analyses suggested a kinase defect for the PINK1 p.I368N mutation would likely be caused by impaired ATP stabilization in the active site and potentially altered binding of the substrate Ub.
PINK1 p.I368N mutant protein is not properly stabilized upon dissipation of the mitochondrial membrane potential
To study the effects of endogenous mutant PINK1 in cells, we analyzed primary skin fibroblasts from controls and both patients carrying the homozygous p.I368N mutation. We quantified PINK1 mRNA levels before and after mitochondrial depolarization with the K+ionophore valinomycin by quantitative reverse transcriptase PCR (qRT-PCR), but found no significant difference between WT and mutant PINK1 fibroblasts (Fig. 3a). Next, we examined PINK1 protein levels under basal conditions and upon mitochondrial damage (Fig. 3b). In the absence of stress, PINK1 protein was hardly detectable by Western blot (WB) under endogenous conditions. While full-length PINK1 WT robustly accumulated in controls upon treatment with the mitochondrial uncoupler CCCP, stabilization of p.I368N mutant PINK1 protein was significantly impaired (Fig. 3b, c). To confirm these findings, fibroblasts were subjected to immunofluorescence (IF) with antibodies against PINK1 and the mitochondrial marker protein TOM20. Consistently, PINK1 signal was considerably weaker in p.I368N mutant cells compared to WT controls (Additional file 5: Figure S1A). Despite significantly reduced protein levels, localization of the mutant PINK1 p.I368N to damaged mitochondria seemed unaltered. This was further corroborated by subcellular fractionations, where both WT and mutant full-length PINK1 were detected in the mitochondrial fraction upon depolarization (Additional file 5: Figure S1B). Taken together, p.I368N mutant cells showed no changes in mRNAs levels, but significantly reduced levels of the full-length PINK1 protein on mitochondria upon depolarization.
The PINK1 p.I368N mutation affects kinase activity and PARKIN activation
To determine the enzymatic activity of the PINK1 p.I368N mutation under endogenous conditions we used recently developed antibodies specific to the PINK1-phosphorylated form of Ub (p-Ser65-Ub) [18, 19, 26]. While Ub was strongly and time-dependently phosphorylated upon CCCP treatment in control cells, no p-Ser65-Ub signal was detected by WB in cells from patient’s carrying the PINK1 p.I368N mutation (Fig. 3b). The same results although with an overall faster stress response were obtained upon depolarization of mitochondria with valinomycin (Additional file 5: Figure S1C). To further confirm this by an independent method, fibroblasts were treated with valinomycin and analyzed by IF (Fig. 3d). In line, PINK1 mutant cells showed no induction of specific p-Ser65-Ub signal above background, which was quantified by high content imaging (HCI) (Fig. 3e and Additional file 5: Figure S1D). Consistent with the lack of p-Ser65-Ub signal in patients’ cells, PARKIN activation/modification and Mitofusin1 (MFN1) ubiquitylation were not observed (see Fig. 3b). Overall, these results show that the PINK1 p.I368N mutation reduces protein stabilization upon stress resulting in failure to activate PARKIN and mitochondrial quality control.
Cleaved forms of PINK1 p.I368N protein are stabilized upon proteasome inhibition
In order to reveal whether differences in PINK1 protein levels are caused by aberrant processing or protein instability, we co-treated cells with a mitochondrial uncoupler and a proteasome inhibitor to stabilize the cleaved forms of PINK1. WB analysis showed that although full-length PINK1 protein (~63 kDa) was significantly reduced in p.I368N mutant cells upon mitochondrial depolarization, the PARL-cleaved form (~52 kDa) was stabilized upon epoxomicin treatment to the same extent as PINK1 WT in the control cells (Fig. 4a, b). Full-length PINK1 and in particularly the PARL-cleaved form has been reported before as an insoluble protein [7, 66, 67]. Thus, we separately analyzed soluble and insoluble fractions of control and p.I368N mutant fibroblasts that had been treated with a combination of valinomycin and epoxomicin. Full-length WT PINK1 was mostly found in the soluble fraction, while its cleaved form strongly accumulated in the insoluble fraction upon proteasome inhibition (Fig. 4c). While the PARL-cleaved form of p.I368N accumulated in the insoluble fraction similar to PINK1 WT, full-length mutant PINK1 remained barely detectable in either fraction. Yet, we could not identify any distinct, additional cleavage product for p.I368N that would indicate aberrant processing of mutant PINK1. Accordingly, several computational prediction tools failed to report de novo protease cleavage sites for PINK1 p.I368N (data not shown).
Protein instability of PINK1 p.I368N is masked upon overexpression
To examine whether full-length PINK1 p.I368N is indeed an inherently unstable protein, we compared WT and mutant levels upon overexpression in HeLa cells. In contrast to endogenous PINK1 in fibroblasts, expression levels of full-length and cleaved p.I368N mutant were surprisingly similar to PINK1 WT upon standard transfection (1 μg DNA per well of a 12-well plate), regardless of mitochondrial depolarization (Additional file 6: Figure S2A). However, transfection of only half of the DNA amounts (0.5 μg versus 1 μg), resulted in decreased levels of PINK1 p.I368N compared to WT protein comparable to the endogenous condition in patients’ fibroblasts. In line with previous studies that had raised caution, our data show that high-level PINK1 overexpression can override important regulatory mechanisms and may mask defects of PINK1 mutations .
Although the exact reasons are unclear, epoxomicin treatment negatively affected stabilization of full-length PINK1 on damaged mitochondria in line with the literature . To prevent any confounding effect of proteasome inhibition, we made use of an engineered mutation that specifically blocks degradation of cleaved PINK1. It is known that PINK1, once imported through the IMM, is cleaved by PARL in its transmembrane region between Ala103 and Phe104, resulting in its degradation by the N-end rule pathway . The F104A mutation is able to block degradation of the PARL- cleaved form of PINK1 independently of the proteasome activity. Upon knockdown of endogenous PINK1 in HeLa cells, we expressed PINK1 siRNA-resistant variants with single or double p.I386N and p.F104A mutations at low levels (Fig. 4d). Under these conditions, full-length PINK1 p.I368N was less stable than PINK1 WT, but similar levels were reached for the p.F104A mutant. Cleaved PINK1 was barely detectable for WT or p.I368N, but strongly accumulated in the single and double mutants harboring the p.F104A mutation. Strikingly, also full-length PINK1 p.I368N was significantly stabilized upon CCCP treatment by the presence of the p.F104A mutation (Fig. 4e). Altogether these results show that while generation, processing and stability of PINK1 p.I368N are not impaired under basal conditions, specifically the mutant full-length form does not properly accumulate on the OMM upon mitochondrial stress.
PINK1 p.I368N mutant has a reduced protein half-life
To better characterize the instability of PINK1 p.I368N over time, we performed experiments in the presence of a protein synthesis inhibitor. PINK1 siRNA knockdown HeLa cells were transfected with only 25% of the recommended DNA amounts (0.25 μg) to ensure low expression of PINK1-V5 close to endogenous levels. Cells were left untreated or treated with cycloheximide (CHX) and the decay of full-length PINK1 protein was monitored in the absence or presence of CCCP, respectively (Fig. 5a and b). Under basal conditions, without mitochondrial damage, the half-life of PINK1-V5 WT was approximately 29 min, while PINK1 p.I368N was stable for only about 7 min (Fig. 5a). Mitochondrial damage dramatically increased the half-life of PINK1-V5 WT (>4 h) and slightly increased also the stability of the p.I368N mutant to 23 min (Fig. 5b). The calculated half-life for PINK1 WT was in good agreement with a previous report . Together this confirms that PINK1 p.I368N is normally expressed, similar to PINK1 WT, but is destabilized on the protein level, in particular during mitochondrial stress.
Full-length PINK1 p.I368N is not stabilized on the OMM/by the HSP90/CDC37 co-chaperone complex
To elucidate the reasons for the destabilization of the p.I368N mutation, we investigated key protein interactions of PINK1 by co-immunoprecipitation (Fig. 5c and d). We focused on the HSP90/CDC37 chaperone complex that supports folding of several kinases in the cytosol including PINK1 . The PD-associated PINK1 mutation p.L347P had been reported before as unstable due to its of lower binding affinity for HSP90/CDC37 [66, 71] and thus was included as a control. While PINK1-V5 WT co-immunoprecipitated substantial amounts of the two chaperones with and without CCCP treatment, both proteins were significantly reduced in immunoprecipitates of PINK1-V5 p.I368N and p.L347P, respectively. In addition, we examined the interaction of PINK1 with the integral OMM protein TOM40 that forms the channel of the protein import complex . It is known that full-length PINK1 accumulates into a higher molecular weight protein complex with the import machinery [10, 73], which regulates its insertion into the OMM, stability and enzymatic function towards PARKIN activation. In line with this, PINK1-V5 WT was found to interact with TOM40 specifically upon CCCP treatment, while the PINK1 mutations p.I368N or p.L347P failed, despite comparable amounts of immunoprecipitated PINK1-V5. Together, these data suggest that the stability of the PINK1 p.I368N mutation is altered through reduced interactions with stabilizing chaperone and the mitochondrial protein import complex.
The PINK1 p.I368N mutant affects kinase activity
To investigate kinase activity of the PINK1 mutation independent of protein instability, HeLa cells were transfected with WT or p.I368N mutant using standard methods (1 μg DNA). At equal expression levels, only PINK1 WT was auto-phosphorylated upon mitochondrial depolarization as shown by separation and detection on phos-tag gel (Fig. 6a). Accordingly, p-Ser65-Ub levels were only increased in HeLa cells transfected with PINK1 WT, but not in cells expressing the p.I368N mutant or the empty vector, also when endogenous PINK1 WT was silenced (Additional file 6: Figure S2B). To directly confirm lack of kinase activity of PINK1 p.I368N, we performed in vitro kinase assays with recombinant Ub. PINK1-V5 WT and p.I368N were immunoprecipitated from overexpressing cells that had been treated with or without CCCP and samples were incubated with biotinylated mono-Ub as a substrate. Phosphorylated Ub was detected exclusively in presence of PINK1 WT obtained from CCCP-treated cells (Fig. 6b).
To confirm a biological downstream effect, we analyzed translocation of GFP-PARKIN to damaged mitochondria in stable HeLa cells by HCI . As expected, siRNA-mediated knockdown of endogenous PINK1 fully abrogated PARKIN translocation to damaged mitochondria (Fig. 6c), while exogenous expression of PINK1 WT restored it but the PINK1 p.I368N mutant failed. To biochemically verify failure to activate PARKIN, PINK1 WT or p.I368N mutant were transiently transfected into HeLa cells stably expressing an active-site mutant of PARKIN. PARKIN C431S traps Ub in an oxyester bond (8 kDa band shift) and facilitates the analysis of the E3 ligase’s Ub-charging as a read-out for its activation . Consistently, significant PARKIN C431S ~ Ub levels were only observed upon re-transfection of cells with PINK1 WT but not with the p.I368N mutant or an empty vector (Fig. 6d, e). Altogether, these data confirm the lack of p.I368N kinase activity even at protein levels comparable to PINK1 WT, and thus its inability to activate PARKIN and the protective mitochondrial quality control.
While EOPD-associated missense mutations in PINK1 and PARKIN typically result in biological loss-of-function, their particular pathomechanisms may differ. Several studies have highlighted specific molecular defects of individual missense variants along a sequential and complex-regulated mitochondrial quality control pathway [5, 22, 24]. In fact, detailed analyses of PARKIN variants have significantly contributed to a better understanding of the pathway itself, the molecular events and their sequence. This is of particular interest not only towards a complete understanding of the PINK1/PARKIN pathway, but also towards a rationalized drug design in the future.
By now numerous missense mutations have been identified in PINK1 from both patients and controls . However, the molecular pathomechanisms of PINK1 mutations have not been well investigated, in particular not on the endogenous level in patients’ cells. We have just recently shown that various mutations disrupt the enzymatic and biological functions of PINK1 at different steps [24, 25]. In addition to simple loss-of-function, specific PINK1 mutations can also exert a partial, dominant-negative function, simply through formation of a heterodimer with PINK1 WT . This results in a significantly increased risk for PD even in heterozygotes. In this study, we revealed two additional molecular mechanisms that can lead to loss of PINK1 function and to increased disease risk.
PINK1 p.I368N was normally expressed, imported into and processed in mitochondria under non-stress conditions, in contrast to PINK1 WT, however, stabilization of the full-length mutant protein on the OMM was strongly impaired. Mitochondrial depolarization combined with proteasome inhibition stabilized similar levels of the PARL-cleaved forms of both PINK1 WT and p.I368N mutant. Yet, proteasome inhibition blocks progression of mitochondrial quality control and N-terminally processed PINK1 itself, which tends to aggregate in the insoluble fraction [7, 66, 67], and has been suggested to block PARKIN translocation from the cytosol . To prevent potentially confounding effects, we employed the PINK1 p.F104A mutation that blocks degradation of the PARL-cleaved form through the N-end rule pathway even without proteasome inhibition . Interestingly, a p.F104A/p.I368N double mutant showed slightly reduced levels of full-length and slightly, but significantly, enhanced levels of cleaved PINK1, compared to the p.F104A mutation alone. This could indicate an overall enhanced and potentially continuous import of PINK1 p.I368N that is coupled to PARL cleavage and degradation by the UPS in the cytosol. It has previously been suggested that the N-terminus of PINK1 is still imported into the inter-membrane space even after depolarization .
In addition to this, we found significantly reduced interactions of PINK1 p.I368N with the cytosolic chaperones HSP90 and CDC37. This chaperone system promotes folding of many kinases and has been shown to regulate stability and subcellular distribution of PINK1 [69, 70, 77, 78]. Mutations such as p.L347P that showed reduced binding to HSP90/CDC37 [71, 79] were suggested before to be subject to misfolding and increased turnover through either the Ub/proteasome system or another protease . Similarly, protein half-life of full-length PINK1 p.I368N was reduced, both before and after mitochondrial stress. While we cannot formally exclude aberrant cleavage of the PINK1 p.I368N mutant in either compartment, a particular instability of the protein in the cytosolic vs. mitochondrial milieu was obvious.
It has been suggested that import of PINK1 into mitochondria is directly regulated by binding to HSP90 and CDC37. The chaperone complex not only assists in folding, stability, and thus kinase activity of PINK1 [66, 69], but also contributes to retention and accumulation of PINK1 on the OMM. As such, a reduced interaction would increase PINK1 inside mitochondria , which in turn would provide more substrate to MPP in the matrix and PARL in the IMM. Of note, loss of those proteases resulted in OMM accumulation of endogenous, full-length PINK1, which induced mitophagy even in the absence of membrane depolarization [8, 9, 76].
In this regard it is important to mention that PINK1 WT interacted with the mitochondrial import channel TOM40 specifically upon CCCP treatment, but both of the mutations that had reduced HSP90/CDC37 binding, failed. It is known that upon mitochondrial depolarization PINK1 auto-phosphorylates [73, 80] and dimerizes into a higher molecular weight protein complex together with the TOM machinery in the OMM [10, 12, 81]. However, it is unclear whether PINK1 then remains in the TOM complex or integrates laterally into the lipid phase of the OMM [9, 10, 76] once it has acquired its orientation to phosphorylate Ub and PARKIN.
In the structure of PINK1, Ile368 is positioned next to the catalytic loop residues 360-367 of the kinase domain. The catalytic loop contains a conserved catalytic base residue Asp362 that is essential for transfer of phosphate from ATP to the substrate . Based on our structural analyses, we conclude that the mutation p.I368N induces local perturbations such that Asp362 cannot adopt the appropriate orientation needed for phosphate transfer. In addition, the de-formed ATP-binding pocket may de-stabilize the interaction which could result in dissociation of ATP from PINK1. All of this together with altered Ub binding would additively decrease the kinase activity of PINK1 p.I368N. While recombinant human PINK1 is notoriously difficult to purify as an active protein in sufficient amounts, it will be important to validate our computational predictions by structural and biophysical comparisons of WT and mutant forms.
Nevertheless, similar defects might be expected to underlie the pathogenicity of other PD associated mutations that are nearby, such as p.N367S or p.L369P. Importantly, a recent study had identified the ATP-analog kinetin to rescue a kinase-defective mutant PINK1 p.G309D . Future rationalized drug design that is based on a detailed understanding of the pathomechanisms should furnish more such small molecules. In addition to active site binders, the complex and as of yet poorly understood life cycle of PINK1 that includes its import, cleavage and re-routing between both mitochondrial membranes might reveal further allosteric sites or strategies to enhance activity and/or retention/stabilization on the OMM.
Some studies have already raised a note of caution as massive overexpression of PINK1 can result in its stabilization on the OMM and subsequent PARKIN recruitment to even healthy mitochondria [11, 16, 68, 80]. In line, our data indicate that subtle defects seen with physiological levels of PINK1 can be masked or overrun by overexpression. Here, we showed that full-length PINK1 p.I368N is unstable on the OMM under endogenous conditions, but this effect was not seen in the overexpression situation. In addition, in a recent study we have found a partial loss-of-function for the PINK1 p.G411S mutation under endogenous conditions , which was missed in an earlier overexpression study .
In summary, we have shown that the PINK1 p.I368N mutation affects protein stabilization specifically on the OMM and results in complete loss of kinase activity towards Ub likely caused by distortion of the ATP-binding site. Together with earlier studies, this highlights the presence of several distinct pathomechanisms that inactivate PINK1 and emphasizes at novel therapeutic avenues that could be pursued in a structure-function based drug design.
Mutations in PINK1 and PARKIN are the most common forms of EOPD. Both proteins play critical roles in stress-induced mitochondrial quality control that protects cells from death. In this study, we examined the PINK1 p.I368N mutation on the clinical and genetic as well as the structural and functional level. We identified two distinct molecular mechanisms that can lead to inactivation of PINK1 function and thus result in loss of neuroprotection. Our research highlights future avenues for drug design that are based on an understanding of the pathogenic events underlying loss of PINK1 activity.
Carbonyl cyanide m-chlorophenyl hydrazone
Early-onset Parkinson’s disease
High content imaging
Inner mitochondrial membrane
Matrix processing peptidase
Molecular dynamics simulation
Mitochondrial targeting sequence
Outer mitochondrial membrane
Presenilin-associated rhomboid-like protein
PTEN-induced putative kinase 1
PINK1 phosphorylated ubiquitin (Serine-65)
Quantitative reverse transcriptase polymerase chain reaction
Root mean square deviation
Root mean square fluctuation
Translocase of the inner mitochondrial membrane
Translocase of the outer mitochondrial membrane
Corti O, Lesage S, Brice A. What Genetics Tells Us About the Causes and Mechanisms of Parkinson's Disease. Physiol Rev. 2011;91:1161–218.
Dawson TM, Ko HS, Dawson VL. Genetic Animal Models of Parkinson's Disease. Neuron. 2010;66:646–61.
Valente EM, Abou-Sleiman PM, Caputo V, Muqit MM, Harvey K, Gispert S, Ali Z, Del Turco D, Bentivoglio AR, Healy DG, et al. Hereditary early-onset Parkinson's disease caused by mutations in PINK1. Science. 2004;304:1158–60.
Kitada T, Asakawa S, Hattori N, Matsumine H, Yamamura Y, Minoshima S, Yokochi M, Mizuno Y, Shimizu N. Mutations in the parkin gene cause autosomal recessive juvenile parkinsonism. Nature. 1998;392:605–8.
Geisler S, Holmstrom KM, Skujat D, Fiesel FC, Rothfuss OC, Kahle PJ, Springer W. PINK1/Parkin-mediated mitophagy is dependent on VDAC1 and p62/SQSTM1. Nat Cell Biol. 2010;12:119–31.
Muqit MM, Abou-Sleiman PM, Saurin AT, Harvey K, Gandhi S, Deas E, Eaton S, Payne Smith MD, Venner K, Matilla A, et al. Altered cleavage and localization of PINK1 to aggresomes in the presence of proteasomal stress. J Neurochem. 2006;98:156–69.
Yamano K, Youle RJ. PINK1 is degraded through the N-end rule pathway. Autophagy. 2013;9:1758–69.
Deas E, Plun-Favreau H, Gandhi S, Desmond H, Kjaer S, Loh SH, Renton AE, Harvey RJ, Whitworth AJ, Martins LM, et al. PINK1 cleavage at position A103 by the mitochondrial protease PARL. Hum Mol Genet. 2011;20:867–79.
Greene AW, Grenier K, Aguileta MA, Muise S, Farazifard R, Haque ME, McBride HM, Park DS, Fon EA. Mitochondrial processing peptidase regulates PINK1 processing, import and Parkin recruitment. EMBO Rep. 2012;13:378–85.
Lazarou M, Jin SM, Kane LA, Youle RJ. Role of PINK1 binding to the TOM complex and alternate intracellular membranes in recruitment and activation of the E3 ligase Parkin. Dev Cell. 2012;22:320–33.
Okatsu K, Oka T, Iguchi M, Imamura K, Kosako H, Tani N, Kimura M, Go E, Koyano F, Funayama M, et al. PINK1 autophosphorylation upon membrane potential dissipation is essential for Parkin recruitment to damaged mitochondria. Nat Commun. 2012;3:1016.
Okatsu K, Kimura M, Oka T, Tanaka K, Matsuda N. Unconventional PINK1 localization to the outer membrane of depolarized mitochondria drives Parkin recruitment. J Cell Sci. 2015;128:964–78.
Kondapalli C, Kazlauskaite A, Zhang N, Woodroof HI, Campbell DG, Gourlay R, Burchell L, Walden H, Macartney TJ, Deak M, et al. PINK1 is activated by mitochondrial membrane potential depolarization and stimulates Parkin E3 ligase activity by phosphorylating Serine 65. Open Biology. 2012;2(5):120080.
Kane LA, Lazarou M, Fogel AI, Li Y, Yamano K, Sarraf SA, Banerjee S, Youle RJ. PINK1 phosphorylates ubiquitin to activate Parkin E3 ubiquitin ligase activity. J Cell Biol. 2014;205:143–53.
Koyano F, Okatsu K, Kosako H, Tamura Y, Go E, Kimura M, Kimura Y, Tsuchiya H, Yoshihara H, Hirokawa T, et al. Ubiquitin is phosphorylated by PINK1 to activate parkin. Nature. 2014;510:162–6.
Okatsu K, Koyano F, Kimura M, Kosako H, Saeki Y, Tanaka K, Matsuda N. Phosphorylated ubiquitin chain is the genuine Parkin receptor. J Cell Biol. 2015;209:111–28.
Ordureau A, Sarraf SA, Duda DM, Heo JM, Jedrychowski MP, Sviderskiy VO, Olszewski JL, Koerber JT, Xie T, Beausoleil SA, et al. Quantitative proteomics reveal a feedforward mechanism for mitochondrial PARKIN translocation and ubiquitin chain synthesis. Mol Cell. 2014;56:360–75.
Fiesel FC, Ando M, Hudec R, Hill AR, Castanedes-Casey M, Caulfield TR, Moussaud-Lamodiere EL, Stankowski JN, Bauer PO, Lorenzo-Betancor O, et al. (Patho-)physiological relevance of PINK1-dependent ubiquitin phosphorylation. EMBO Rep. 2015;16:1114–30.
Fiesel FC, Springer W. Disease relevance of phosphorylated ubiquitin (p-S65-Ub). Autophagy. 2015;11:2125–6.
Heo JM, Ordureau A, Paulo JA, Rinehart J, Harper JW. The PINK1-PARKIN Mitochondrial Ubiquitylation Pathway Drives a Program of OPTN/NDP52 Recruitment and TBK1 Activation to Promote Mitophagy. Mol Cell. 2015;60:7–20.
Lazarou M, Sliter DA, Kane LA, Sarraf SA, Wang C, Burman JL, Sideris DP, Fogel AI, Youle RJ. The ubiquitin kinase PINK1 recruits autophagy receptors to induce mitophagy. Nature. 2015;524:309–14.
Fiesel FC, Caulfield TR, Moussaud-Lamodiere EL, Ogaki K, Dourado DF, Flores SC, Ross OA, Springer W. Structural and Functional Impact of Parkinson Disease-Associated Mutations in the E3 Ubiquitin Ligase Parkin. Hum Mutat. 2015;36:774–86.
Caulfield TR, Fiesel FC, Springer W. Activation of the E3 ubiquitin ligase Parkin. Biochem Soc Trans. 2015;43:269–74.
Geisler S, Holmstrom KM, Treis A, Skujat D, Weber SS, Fiesel FC, Kahle PJ, Springer W. The PINK1/Parkin-mediated mitophagy is compromised by PD-associated mutations. Autophagy. 2010;6:871–8.
Siuda J, Jasinska-Myga B, Boczarska-Jedynak M, Opala G, Fiesel FC, Moussaud-Lamodiere EL, Scarffe LA, Dawson VL, Ross OA, Springer W, et al. Early-onset Parkinson's disease due to PINK1 p.Q456X mutation--clinical and functional study. Parkinsonism Relat Disord. 2014;20:1274–8.
Puschmann A, Fiesel FC, Caulfield TR, Hudec R, Ando M, Truban D, Hou X, Ogaki K, Heckman MG, James ED, et al. Heterozygous PINK1 p.G411S increases risk of Parkinson's disease via a dominant-negative mechanism. Brain. 2017;140(Pt 1):98-117.
Hughes AJ, Daniel SE, Kilford L, Lees AJ. Accuracy of clinical diagnosis of idiopathic Parkinson's disease: a clinico-pathological study of 100 cases. J Neurol Neurosurg Psychiatry. 1992;55:181–4.
Mayo Clinic. Department of Neurology. Mayo Clinic examinations in neurology. 7th ed. St. Louis: Mosby; 1998.
Hoehn MM, Yahr MD. Parkinsonism: onset, progression and mortality. Neurology. 1967;17:427–42.
Folstein MF, Folstein SE, McHugh PR. "Mini-mental state". A practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res. 1975;12:189–98.
Jacobson MP, Friesner RA, Xiang Z, Honig B. On the role of the crystal environment in determining protein side-chain conformations. J Mol Biol. 2002;320:597–608.
Krieger E, Joo K, Lee J, Lee J, Raman S, Thompson J, Tyka M, Baker D, Karplus K. Improving physical realism, stereochemistry, and side-chain accuracy in homology modeling: Four approaches that performed well in CASP8. Proteins. 2009;77 Suppl 9:114–22.
Altschul SF, Koonin EV. Iterated profile searches with PSI-BLAST--a tool for discovery in protein databases. Trends Biochem Sci. 1998;23:444–7.
Hooft RW, Sander C, Scharf M, Vriend G. The PDBFINDER database: a summary of PDB, DSSP and HSSP information with added value. Comput Appl Biosci. 1996;12:525–9.
Hooft RW, Vriend G, Sander C, Abola EE. Errors in protein structures. Nature. 1996;381:272.
King RD, Sternberg MJ. Identification and application of the concepts important for accurate and reliable protein secondary structure prediction. Protein Sci. 1996;5:2298–310.
Krieger E, Dunbrack Jr RL, Hooft RW, Krieger B. Assignment of protonation states in proteins and ligands: combining pKa prediction with hydrogen bonding network optimization. Methods Mol Biol. 2012;819:405–21.
Qiu J, Elber R. SSALN: an alignment algorithm using structure-dependent substitution matrices and gap penalties learned from structurally aligned protein pairs. Proteins. 2006;62:881–91.
Zhou H, Pandit SB, Lee SY, Borreguero J, Chen H, Wroblewska L, Skolnick J. Analysis of TASSER-based CASP7 protein structure prediction results. Proteins. 2007;69 Suppl 8:90–7.
Zhou H, Pandit SB, Skolnick J. Performance of the Pro-sp3-TASSER server in CASP8. Proteins. 2009;77 Suppl 9:123–7.
Zhou H, Skolnick J. Ab initio protein structure prediction using chunk-TASSER. Biophys J. 2007;93:1510–8.
Zhou H, Skolnick J. Protein structure prediction by pro-Sp3-TASSER. Biophys J. 2009;96:2119–27.
Zhou H, Skolnick J. Improving threading algorithms for remote homology modeling by combining fragment and template comparisons. Proteins. 2010;78:2041–8.
Zhou H, Skolnick J. Template-based protein structure modeling using TASSER(VMT). Proteins. 2011;80(2):352-61. doi:10.1002/prot.23183.
Maestro-9.1. Maestro 9.1. 91 Mth ed. New York: Schrödinger, LLC; 2010.
Schrödinger. Suite 2012. BioLuminate, version 10th ed. New York: Schrödinger, LLC; 2012.
Schrödinger. Biologics Suite. In: Schrodinger, editor. BioLuminate, version 11. BioLuminate, version 11th ed. New York: Schrödinger, LLC; 2013.
Guerois R, Nielsen JE, Serrano L. Predicting changes in the stability of proteins and protein complexes: a study of more than 1000 mutations. J Mol Biol. 2002;320:369–87.
Brandsdal BO, Aqvist J, Smalas AO. Computational analysis of binding of P1 variants to trypsin. Protein Sci. 2001;10:1584–95.
Caulfield T, Devkota B. Motion of transfer RNA from the A/T state into the A-site using docking and simulations. Proteins. 2012;80:2489–500.
Caulfield T, Medina-Franco JL. Molecular dynamics simulations of human DNA methyltransferase 3B with selective inhibitor nanaomycin A. J Struct Biol. 2011;176:185–91.
Caulfield TR. Inter-ring rotation of apolipoprotein A-I protein monomers for the double-belt model using biased molecular dynamics. J Mol Graph Model. 2011;29:1006–14.
Caulfield TR, Devkota B, Rollins GC. Examinations of tRNA Range of Motion Using Simulations of Cryo-EM Microscopy and X-Ray Data. J Biophys. 2011;2011:219515.
Case DA, Cheatham III TE, Darden T, Gohlke H, Luo R, Merz Jr KM, Onufriev A, Simmerling C, Wang B, Woods R. The Amber biomolecular simulation programs. J Chemical Chemistry. 2005;26:1668–88.
Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, Chipot C, Skeel R, Kale L, Schulten K. Scalable molecular dynamics with NAMD. J Comput Chem. 2005;26(16):1781-802.
Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. Comparison of simple potential functions for simulating liquid water. J Chemical Physics. 1983;79:926–35.
Reblova K, Fadrna E, Sarzynska J, Kulinski T, Kulhanek P, Ennifar E, Koca J, Sponer J. Conformations of flanking bases in HIV-1 RNA DIS kissing complexes studied by molecular dynamics. Biophys J. 2007;93:3932–49.
Cheatham III TE, Young MA. Molecular dynamics simulation of nucleic acids: Successes, limitations and promise. Biopolymers. 2001;56:232–56.
Caulfield TR, Fiesel FC, Moussaud-Lamodiere EL, Dourado DF, Flores SC, Springer W. Phosphorylation by PINK1 releases the UBL domain and initializes the conformational opening of the E3 ubiquitin ligase Parkin. PLoS Comput Biol. 2014;10:e1003935.
Kuck D, Caulfield T, Lyko F, Medina-Franco JL. Nanaomycin A selectively inhibits DNMT3B and reactivates silenced tumor suppressor genes in human cancer cells. Mol Cancer Ther. 2010;9:3015–23.
Fiesel FC, Moussaud-Lamodiere EL, Ando M, Springer W. A specific subset of E2 ubiquitin-conjugating enzymes regulate Parkin activation and mitophagy differently. J Cell Sci. 2014;127:3488–504.
Koziorowski D, Hoffman-Zacharska D, Slawek J, Jamrozik Z, Janik P, Potulska-Chromik A, Roszmann A, Tataj R, Bal J, Friedman A. Incidence of mutations in the PARK2, PINK1, PARK7 genes in Polish early-onset Parkinson disease patients. Neurol Neurochir Pol. 2013;47:319–24.
Valente EM, Ferraris A. Blue Books of Neurology. Chapter 5 PINK1 (PARK6) and Parkinson's Disease. In: Schapira AHV ,Lang AET, Fahn S, editors. Elsevier. 2010;34:66–82. http://doi.org/10.1016/B978-1-4160-6641-5.00005-2. ISBN: 978-1-4160-6641-5.
Postuma RB, Gagnon JF, Bertrand JA, Genier Marchand D, Montplaisir JY. Parkinson risk in idiopathic REM sleep behavior disorder: preparing for neuroprotective trials. Neurology. 2015;84:1104–13.
Chaudhuri KR, Tolosa E, Schapira AHV, Poewe W. Non-Motor Symptoms of Parkinson's Disease. USA: Oxford University Press; 2014.
Beilina A, Van Der Brug M, Ahmad R, Kesavapany S, Miller DW, Petsko GA, Cookson MR. Mutations in PTEN-induced putative kinase 1 associated with recessive parkinsonism have differential effects on protein stability. Proc Natl Acad Sci U S A. 2005;102:5703–8.
Silvestri L, Caputo V, Bellacchio E, Atorino L, Dallapiccola B, Valente EM, Casari G. Mitochondrial import and enzymatic activity of PINK1 mutants associated to recessive parkinsonism. Hum Mol Genet. 2005;14:3477–92.
Narendra DP, Jin SM, Tanaka A, Suen DF, Gautier CA, Shen J, Cookson MR, Youle RJ. PINK1 is selectively stabilized on impaired mitochondria to activate Parkin. PLoS Biol. 2010;8:e1000298.
Lin W, Kang UJ. Characterization of PINK1 processing, stability, and subcellular localization. J Neurochem. 2008;106:464–74.
Weihofen A, Ostaszewski B, Minami Y, Selkoe DJ. Pink1 Parkinson mutations, the Cdc37/Hsp90 chaperones and Parkin all influence the maturation or subcellular distribution of Pink1. Hum Mol Genet. 2008;17:602–16.
Moriwaki Y, Kim YJ, Ido Y, Misawa H, Kawashima K, Endo S, Takahashi R. L347P PINK1 mutant that fails to bind to Hsp90/Cdc37 chaperones is rapidly degraded in a proteasome-dependent manner. Neurosci Res. 2008;61:43–8.
Gabriel K, Egan B, Lithgow T. Tom40, the import channel of the mitochondrial outer membrane, plays an active role in sorting imported proteins. EMBO J. 2003;22:2380–6.
Okatsu K, Uno M, Koyano F, Go E, Kimura M, Oka T, Tanaka K, Matsuda N. A dimeric PINK1-containing complex on depolarized mitochondria stimulates Parkin recruitment. J Biol Chem. 2013;288:36372–84.
Truban D, Hou X, Caulfield TR, Fiesel FC, Springer W. PINK1, Parkin, and mitochondrial quality control – What can we learn about Parkinson's disease pathobiology? J Parkinsons Dis. 2017;7(1):13-29.
Fedorowicz MA, de Vries-Schneider RL, Rub C, Becker D, Huang Y, Zhou C, Alessi Wolken DM, Voos W, Liu Y, Przedborski S. Cytosolic cleaved PINK1 represses Parkin translocation to mitochondria and mitophagy. EMBO Rep. 2014;15:86–93.
Meissner C, Lorenz H, Hehn B, Lemberg MK. Intramembrane protease PARL defines a negative regulator of PINK1- and PARK2/Parkin-dependent mitophagy. Autophagy. 2015;11:1484–98.
Sim CH, Gabriel K, Mills RD, Culvenor JG, Cheng HC. Analysis of the regulatory and catalytic domains of PTEN-induced kinase-1 (PINK1). Hum Mutat. 2012;33:1408–22.
Cardona F, Sanchez-Mut JV, Dopazo H, Perez-Tur J. Phylogenetic and in silico structural analysis of the Parkinson disease-related kinase PINK1. Hum Mutat. 2011;32:369–78.
Lin W, Kang UJ. Structural determinants of PINK1 topology and dual subcellular distribution. BMC Cell Biol. 2010;11:90.
Aerts L, Craessaerts K, De Strooper B, Morais VA. PINK1 kinase catalytic activity is regulated by phosphorylation on serines 228 and 402. J Biol Chem. 2015;290:2798–811.
Liu W, Vives-Bauza C, Acin-Perez R, Yamamoto A, Tan Y, Li Y, Magrane J, Stavarache MA, Shaffer S, Chang S, et al. PINK1 defect causes mitochondrial dysfunction, proteasomal deficit and alpha-synuclein aggregation in cell culture models of Parkinson's disease. PLoS One. 2009;4:e4597.
Hertz NT, Berthet A, Sos ML, Thorn KS, Burlingame AL, Nakamura K, Shokat KM. A neo-substrate that amplifies catalytic activity of parkinson's-disease-related kinase PINK1. Cell. 2013;154:737–47.
We thank all patients and controls that participated in our research studies.
W.S. is partially supported by the National Institutes of Health (NIH)/National Institute of Neurological Disorders and Stroke (NINDS) [R01 #NS085070], the Michael J. Fox Foundation for Parkinson’s Research and the Foundation for Mitochondrial Medicine, Mayo Clinic Center for Individualized Medicine (CIM), Center for Regenerative Medicine (CRM) and the Center for Biomedical Discovery (CBD), the Marriott Family Foundation, and a Gerstner Family Career Development Award. F.C.F. is the recipient of a fellowship from the Younkin Scholar Program and has been supported by the American Parkinson Disease Association (APDA). M.A. acknowledges support from the Uehara Memorial Foundation. Mayo Clinic Florida is a Morris K. Udall Parkinson’s Disease Research Center of Excellence [NIH/NINDS P50 #NS072187 to Z.K.W. and O.A.R.]. O.A.R. is supported by NIH/NINDS [R01 #NS078086]. Z.K.W. and O.A.R are supported by the Mayo Clinic Center for Individualized Medicine (CIM) and Z.K.W. is supported by the Mayo Clinic Center for Regenerative Medicine (CRM). Z.K.W., O.A.R., and W.S. are recipients of a Mayo Clinic Neuroscience Focused Research Team Award. GB is supported by the Mayo Clinic Center for Regenerative Medicine (CRM) Neuroregeneration Lab. V.L.D is supported by grants from the NIH/NINDS NS38377, MDSCRF 2007-MSCRFI-0420-00, 2009-MSCRFII-0125-00, MDSCRF 2012-MSCRFII-0268-00, MDSCRF 2013-MSCRFII-0105-00. T.M.D. is the Leonard and Madlyn Abramson Professor in Neurodegenerative Diseases, and is supported by grants from the NIH/NINDS NS38377, MDSCRF 2007-MSCRFI-0420-00, 2009-MSCRFII-0125-00, MDSCRF 2012-MSCRFII-0268-00, MDSCRF 2013-MSCRFII-0105-00, and the JPB Foundation. T.M.D. and V.L.D. acknowledge the joint participation by the Adrienne Helis Malvin Medical Research Foundation through their direct engagement in the continuous active conduct of medical research in conjunction with the Johns Hopkins Hospital and the Johns Hopkins University School of Medicine and the Foundation’s Parkinson’s Disease Program M-1.
Availability of data and material
The datasets supporting the conclusions of this article are included within the article and its additional files.
MA, FCF, RH, and TRC performed experiments and analyzed the data. KO, PG-S, DK, AF, LC, VLD, TMD, GB, OAR and ZKW contributed clinical and genetic information as well as materials of the patients. WS conceived and designed the study. FCF and WS wrote the manuscript. All authors contributed to drafting or revising the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The research project was approved by the Ethical Committee of the Medical University of Warsaw and by the IRB Committee of the Mayo Clinic, and both Committees approved the collaboration between institutions and researchers. Both research subjects signed the informed consent to participate in this study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Clinical features of PD patients homozygous for the PINK1 p.I368N mutation are summarized. Abbreviations: RBD: REM sleep behavior disorder; MMSE: Mini Mental State Examination; DBS: deep brain stimulation. (DOCX 16 kb)
FoldX measurement for mutational energy effect on stability of binding pocket. Using the FoldX algorithm, which compares the WT structure with the mutant, calculations on WT and mutant at t = 0 and after 100 ns were performed. While WT is nearly zero over that time interval, the N368 mutant increases energy (∆∆G) by 1.5–2 kcal/mol*Å2. This modest increase supports that the mutation distorts the local region via an increase of Gibb’s free energy. (DOCX 16 kb)
Reduced full-length PINK1 levels in p.I368N mutant fibroblasts upon treatment with valinomycin. (A) Representative confocal IF images of control and two PINK1 p.I368N fibroblasts showing mitochondrial localization but greatly reduced levels of the mutant protein. Cells were left untreated (-) or treated with 10 μM CCCP for 6 days (+) and stained with antibodies against PINK1 (green) and TOM20 (mitochondria, red). Nuclei were counterstained with Hoechst (blue). Scale bars represent 10 μm. (B) Subcellular fractionation of WT control and PINK1 p.I368N fibroblasts treated with or without 1 μM valinomycin for 24 h. Despite different protein levels, both WT and p.I368N mutant full-length PINK1 localized to the mitochondrial fraction. A shift of PARKIN into higher molecular weights species indicative of PINK1-dependent activation was not observed in lysates from PINK1 p.I368N mutant cells. Purity of the mitochondrial and cytosolic fractions was determined using antibodies recognizing TOM20 and p38 MAPK, respectively. PNS denotes post-nuclear supernatant. (C) Control and two PINK1 p.I368N fibroblasts were treated with 1 μM valinomycin for 0, 2, 4 or 8 h and total lysates were analyzed by WB with indicated antibodies. GAPDH served as a loading control. Similar to CCCP treatment (Fig. 3b), rapid stabilization of full-length PINK1 along with an increase of p-Ser65-Ub levels was observed in control cells, but not in PINK1 p.I368N mutant fibroblasts. (D) Representative Images obtained with a 20x magnification on the BD pathway 855 system. Control fibroblasts and PINK1 I368N cells were seeded in 96-well imaging plates and treated with valinomycin as indicated. Cells were fixed and stained with p-S65-Ub antibodies (green) and Hoechst (blue). Scale bars indicate 10 μm. (TIF 13839 kb)
Figure S2. Expression of PINK1-V5 WT and p.I368N mutation at different levels in HeLa cells. (A) HeLa cells were transfected with PINK1 siRNA and with different amounts (0.5 or 1 μg per well of a 12-well plate) of V5 empty vector, PINK1-V5 WT or p.I368N, as indicated. Cells were left untreated or incubated with 10 μM CCCP for 4 h and lysates were analyzed on WB with anti-V5 and PINK1 antibodies. GAPDH was used as a loading control. V5/GAPDH ratios were determined by densitometry and normalized to values of CCCP treated PINK1 WT samples transfected with 1 μg of DNA. (B) Denitometric analysis of data from (A). Shown is the mean ± SEM from four independent experiments. Statistical significance was assessed by two-way ANOVA with Tukey’s post hoc; *, p < 0.05; ***, p < 0.0005; ns, not significant. (C) HeLa cells were transfected with control or PINK1 siRNA and with V5 empty vector, PINK1-V5 WT or p.I368N. Cells were treated with CCCP for 4 h and proteins were extracted and analyzed by WB using the indicated antibodies. p-Ser65-Ub signal was detected only in PINK1-V5 WT but not in p.I368N transfected cells upon mitochondrial stress despite similar PINK1 proteins levels. (TIF 3542 kb)