If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Molecular Medicine Program, Taiwan International Graduate Program, Institute of Biomedical Sciences, Academia Sinica, Taipei 115, Taiwan;Institute of Biochemistry and Molecular Biology, School of Life Sciences, National Yang-Ming University, Taipei 112, Taiwan;
Molecular Medicine Program, Taiwan International Graduate Program, Institute of Biomedical Sciences, Academia Sinica, Taipei 115, Taiwan;Institute of Biochemistry and Molecular Biology, School of Life Sciences, National Yang-Ming University, Taipei 112, Taiwan;
To whom correspondence should be addressed: Institute of Biomedical Sciences, Academia Sinica, 128 Academia Rd Sec 2, Nankang, Taipei 11529, Taiwan. Tel.: 886-2-2652-3522; Fax: 886-2-2785-8847;
* This work was supported by the National Science Council, Taiwan (98-2320-B-001-017-MY3, 100-2311-B-001-002-MY3) and Academia Sinica (AS100-CDA-L05). This article contains supplemental Figs. S1 to S6 and Tables S1 to S3.
Pressure-overload stress to the heart causes pathological cardiac hypertrophy, which increases the risk of cardiac morbidity and mortality. However, the detailed signaling pathways induced by pressure overload remain unclear. Here we used phosphoproteomics to delineate signaling pathways in the myocardium responding to acute pressure overload and chronic hypertrophy in mice. Myocardial samples at 4 time points (10, 30, 60 min and 2 weeks) after transverse aortic banding (TAB) in mice underwent quantitative phosphoproteomics assay. Temporal phosphoproteomics profiles showed 360 phosphorylation sites with significant regulation after TAB. Multiple mechanical stress sensors were activated after acute pressure overload. Gene ontology analysis revealed differential phosphorylation between hearts with acute pressure overload and chronic hypertrophy. Most interestingly, analysis of the cardiac hypertrophy pathway revealed phosphorylation of the mitochondrial fission protein dynamin-related protein 1 (DRP1) by prohypertrophic kinases. Phosphorylation of DRP1 S622 was confirmed in TAB-treated mouse hearts and phenylephrine (PE)-treated rat neonatal cardiomyocytes. TAB-treated mouse hearts showed phosphorylation-mediated mitochondrial translocation of DRP1. Inhibition of DRP1 with the small-molecule inhibitor mdivi-1 reduced the TAB-induced hypertrophic responses. Mdivi-1 also prevented PE-induced hypertrophic growth and oxygen consumption in rat neonatal cardiomyocytes. We reveal the signaling responses of the heart to pressure stress in vivo and in vitro. DRP1 may be important in the development of cardiac hypertrophy.
Hypertension or acute aortic stenosis causes pressure overload in the left ventricle (LV)
, and prolonged pressure overload leads to pathological cardiac hypertrophy. Although beneficial initially, hypertrophy in the heart signals metabolic, structural and functional abnormalities and might lead to abnormal cardiac function (
). However, one of the major challenges but also the prerequisite to preventing cardiac hypertrophy is a comprehensive understanding of the mechanism to precisely prevent pathologic cardiac growth without affecting homeostasis (
The signaling pathways leading to cardiac hypertrophy with chronic pressure overload have been extensively studied by genetic and pharmacological means (
Pressure-overload hypertrophy of the developing heart reveals activation of divergent gene and protein pathways in the left and right ventricular myocardium.
Am. J. Physiol. Heart Circ. Physiol.2013; 304: H697-708
). However, the signaling pathways contributing to early changes in cardiomyocytes remain unclear.
Protein phosphorylation allows cells to quickly respond to stimuli and transmit signals by regulating enzymatic activity, protein subcellular localization, protein interaction partners, and protein stability (
). Abnormal phosphorylation of proteins has been associated with many diseases, including cardiovascular diseases. For example, hyperphosphorylation of type 2 ryanodine receptor (Ryr2) by protein kinase A leads to defective channel function in the human heart (
Despite the importance of phosphorylation-mediated regulation in heart diseases, the multiple signaling pathways of cardiac hypertrophy at the phosphoproteome scale have not been delineated. Furthermore, insights into the dynamic interplay of such pathways in vivo could greatly enhance our understanding of the molecular mechanism of acute pressure overload and cardiac hypertrophy. In this study, we used quantitative phosphoproteomics to reveal the early signaling pathways induced by acute pressure overload in the mouse LV. Low abundant phosphopeptides were enriched by immobilized metal affinity chromatography (IMAC). To facilitate accurate quantification of phosphorylation in vivo, we used a post-enrichment labeling with isobaric tag for relative and absolute quantification (iTRAQ) for quantitative phosphoproteomics and demonstrated reliable quantitation performance with ≈10% coefficient of variation (CV). This strategy provided a large-scale quantification of phosphorylation change of LV proteins at four time points (10, 30, 60 min and 2 weeks) after transverse aortic banding surgery (TAB) in mice. This study revealed potential signal pathways underlying the pressure stress response and the disease phenotypes during the progression of cardiac hypertrophy. We further demonstrated that mitochondrial fission protein dynamin-related protein 1 (DRP1) is involved in the pathological cardiac hypertrophy.
EXPERIMENTAL PROCEDURES
Animal, Transverse Aortic Banding Surgery (TAB), Mitochondrial Division Inhibitor 1 (mdivi-1) Injection, and Echocardiography Analysis
Eight-week-old C57BL/6 male mice (20–25 g) underwent pressure overload by transverse aortic banding (TAB) or sham operation as described (
). For acute TAB, experiments were repeated three times with three mice for each time tested in each replicate (Fig. 1A, Exp. Set 1). The hypertrophy experiment involved two replicates with two mice each for TAB and sham operation at 2 weeks in each replicate (Fig. 1A, Exp. Set 2).
Fig. 1Quantitative phosphoproteomics study of pressure-overloaded hearts of mice.A, Strategy for quantitative phosphoproteomic studies of mice with pressure-overloaded hearts induced by transverse aortic banding (TAB). Experimental set 1 (Exp. Set 1) and experimental set 2 (Exp. Set 2) were designed for studying the phosphorylation changes in acute-stressed and hypertrophic hearts, respectively. Phosphopeptide profiles from left-ventricle (LV) tissues of (B) acute pressure-overloaded and (C) chronic hypertrophic hearts (2 weeks after TAB) as compared with sham-operated hearts. The ratio profile of the 1:1 experiment represents the distribution of unchanged ratios. D, Immunoblotting confirmation of regulated phosphorylation sites in αB-Crystallin (Cryab), heat shock protein 27 (Hspb1), phospholamban (PLN) and Na/K ATPase α1A subunit (ATP1a1). Quantification is normalized to total protein. Data are mean ± S.E. * p < 0.05 versus sham. The position of molecular weight marker is showed to the right of immunoblot.
Mice received an intraperitoneal injection of 25 mg/Kg mdivi-1 dissolved in DMSO every 2 days. Vehicle control mice received an intraperitoneal injection of DMSO every 2 days. Before animals were killed, the pressure gradient across the banding site was checked by echocardiography to ensure the pressure overload (
LV Protein Extraction, Mitochondrial Purification, RNA Extraction and Real-Time Quantitative PCR
After the indicated time of TAB, mice were anesthetized for 3 min by isoflurane (3% in oxygen), then killed by neck dislocation. Hearts were excised and weighed, then washed with ice-cold phosphate buffered saline with a phosphatase inhibitor. Atria and right ventricles were removed. The time from neck dislocation to obtain the LV was within 5 min. The LV was frozen by use of liquid nitrogen. For LV protein extraction, 1 ml lysis buffer (20 mm Tris-HCl pH 7.5, 150 mm NaCl, 1 mm Na2EDTA, 1 mm EGTA, 1% Triton X-100, 2.5 mm sodium pyrophosphate, 1 mm beta-glycerophosphate, 1 mm Na3VO4, 1 μg/ml leupeptin) with phosphatase inhibitor mixture (Pierce) was added to each LV for homogenization by use of a Dounce homogenizer (Wheaton). The lysates were incubated on ice for 30 min. Undissolved pellets were centrifuged at 13,000 rpm for 30 min at 4 °C. For mitochondrial isolation, LV was homogenized in buffer containing 250 mm sucrose, 20 mm HEPES pH 7.5, 10 mm KCl, 1.5 mm MgCl2, 1 mm EDTA, protease inhibitor mixture (Sigma) and phosphatase inhibitor mixture. The homogenates were centrifuged first at 800 × g for 10 min at 4 °C and supernatant was further centrifuged at 10,000 × g for 20 min at 4 °C. The pellets were washed with buffer and spun at 10,000 × g for 20 min at 4 °C. The final pellets were dissolved in 1% Triton X-100 lysis buffer. Protein concentrations were determined by the Bradford method (Bio-Rad). Protein lysates were stored at −80 and analyzed within 2 months. Extraction of total RNA from LV and real-time quantitative PCR were as described (
Protein Digestion, Phosphopeptide Enrichment, and Isobaric Tag for Relative and Absolute Quantification (iTRAQ) Labeling
For protein digestion, 1 mg LV lysates from each mouse was used, and 0.75 μg α-casein and 0.25 μg β-casein were added. Gel-assisted digestion was used to improve the digestion efficiency of membrane proteins (
A multiplexed quantitative strategy for membrane proteomics: opportunities for mining therapeutic targets for autosomal dominant polycystic kidney disease.
). The IMAC elution was cleaned by use of C18 Ziptip (Millipore). Phosphopeptides were dissolved by use of 10 μl H2O, and 0.5 μl was added into 100 μl BCA assay reagent (Pierce). The standard curve was created by adding 0.0625, 0.125, 0.1875, 0.25, 0.5, and 1.0 μg BSA into 100 μl BCA assay reagent. After the weight of phosphopeptides was determined, iTRAQ reagent (
) was added to at least twofold the recommended amount. For convenience, 10 μl reconstituted iTRAQ reagent was added into solution containing 9.5 μl phosphopeptides, 2.5 μl TEABC (1 m), and 18 μl ethanol. After 2 h at room temperature, phosphopeptides labeled with different iTRAQ tags were mixed and dried by use of a vacuum pump.
Mass Spectrometry (MS), Database Searching and Phosphopeptide Quantification
Each sample was analyzed at least three times by reverse-phase ultra-performance liquid chromatography (Waters nanoACQUITY UPLC) with quadruple time-of-flight MS (Waters Q-TOF Premier) as described (
). The MS peak lists in Mascot generic format (mgf) were generated by use of Mascot Distiller v2.1.1.0 with default parameters. The combined mgf file was used for a search against the protein database IPI Mouse v3.64 (56,751 protein entries) with an in-house version of Mascot v2.2 (Matrix Science, London, UK). Parameter settings for database searching were trypsin as protease, up to 2 missed cleavages and 0.07 Da for both precursor and fragment ions tolerance. iTRAQ modifications of lysine and N-terminal were set as fixed modification. Methylthio of cysteine, oxidation of methionine, and phosphorylation of serine, threonine, and tyrosine were set as variable modifications.
Proteins were considered identified by significance threshold p < 0.05. Peptides were identified by peptide score ≥15. Peptides matched to multiple protein hits were assigned to the protein hit with highest protein score. The false positive rate for peptide identification was 0.62% for acute TAB and 1.46% for TAB at 2 weeks as estimated by Mascot searching against a randomized decoy database. The Mascot delta score was used for phosphorylation site assignment (
). In all, 1,603 MSMS spectra for the 2-week TAB experiment were manually checked. Among 1,462 MSMS spectra with a Mascot delta score ≥ 5 for the first two hits, phosphorylation sites for only one spectrum could not be determined manually. Thus, we used the phosphorylation site assigned by Mascot if the delta score for the first two hits was ≥5.
The quantification of iTRAQ labeling peptide involved use of Multi-Q (
). The mascot protein or peptide identification information was saved as an XML file. The MS and MSMS spectra information was converted from Waters (Milford, MA) MassLynx raw data to an mzXML file by use of massWolf software. The XML and mzXML files were loaded into Multi-Q for automatic calculation. The quantifiable MSMS spectra were defined as the peak area of the iTRAQ reporter ion > 25 determined by the 1:1 experiment (supplemental Fig. S2A and B).
Three externally added phosphoproteins, bovine α-S1, α-S2 and β-casein, were an internal reference for loss of LV phosphopeptides (supplemental Fig. S2C). The reproducibility of quantification for every comparison pair was estimated by the coefficient of variation (CV) of the ratio of casein phosphopeptides in each comparison (supplemental Fig. S3). The LV phosphopeptide ratio was first adjusted by the ratio of casein phosphopeptides, then adjusted by the total peptide ratio to normalize protein loading error. The total trypsin-digested peptide was also labeled with iTRAQ and then analyzed by MS, identified by a database search and quantified by Multi-Q. The signal intensity of the iTRAQ reporter ion for all peptides in each sample was summed and used to normalize the protein input error of each sample. In all, 5204 MSMS spectra were used for quantification. The ratios of MS/MS spectra matched to the same phosphopeptide were averaged. Phosphopeptide ratios of different mice were further averaged and showed as mean ± S.D. The MS/MS spectra are available at http://www.peptideatlas.org/PASS/PASS00236.
GoMiner Analysis, Motif Classification, Cardiac Phenotype Searching, and Hypertrophy Pathway Construction
Subcellular localization and biological processes of phosphoproteins and classes with significant change in phosphorylation were determined by use of the high-throughput GoMiner tool (
High-Throughput GoMiner, an ‘industrial-strength’ integrative gene ontology tool for interpretation of multiple-microarray experiments, with application to studies of Common Variable Immune Deficiency (CVID).
). p values and FDR for each GO term were calculated, and significance of changes was considered at p < 0.05. The classification of phosphorylation sequences by general kinase recognition motifs was as described (
). The cardiac phenotype of phosphoproteins was searched in the Mouse Genome Informatics (MGI) database at http://www.informatics.jax.org/. The network of phosphoproteins related to cardiac hypertrophic growth was constructed manually, and the references for events in the network are in supplemental Table S2.
Immunoblotting, Immunofluorescence, and Immunoprecipitation
Samples were incubated with the antibodies for extracellular signal-regulated kinase 1/2 (ERK1/2) T202/Y204, S6 S235/S236, heat shock protein b1 (Hspb1) S82, ATPase α1A subunit (ATP1a1) S16, DRP1 S616, DRP1, protein kinase C δ (PKC-δ), and HA-tag (Cell Signaling Technology); phospholamban (PLN) T17, tubulin, and GAPDH (Santa Cruz Biotechnology, Santa Cruz, CA); mitochondrially encoded cytochrome c oxidase (MT-CO1; Abcam, Cambridge, MA); crystallin αB (Cryab) S59 (Enzo Life Sciences, Farmingdale, NY); and α-actinin (Sigma-Aldrich). The mouse Xin antibody clone U1013 was prepared as described (
). Preparation of the rabbit polyclonal phospho-Xin S205/S208 and S295 antibodies involved corresponding synthetic peptides with phosphoserine in the indicated sites. Red fluorescence wheat germ agglutinin (WGA) was from Invitrogen. Immunoblotting, immunoprecipitation, and immunofluorescence were as described (
Primary Culture of Rat Neonatal Cardiomyocytes (rNCMs), Transfection, and Measurement of Oxygen Consumption Rate
rNCMs were isolated from ventricles of 1- to 3-day-old rats. The ventricles were digested by use of collagenase and trypsin (Worthington Biochemical Corp, Freehold, NJ) as described (
). Isolated cells were plated in dishes coated with collagen and polylysine (Sigma-Aldrich) in DMEM containing newborn bovine serum (10%), penicillin (100 units/ml) and streptomycin (100 μg/ml) and B-d-arabinofuranoside cytosine (17 μm) for 2 days. After overnight incubation in serum-free medium, rNCMs were treated with PE (10 μm) for 24 h to induce cardiac hypertrophy at 37 °C. The HA-DRP1K38A construct (a kind gift from Dr. Alexander M. Van der Bliek in UCLA) was transfected into rNCMs by use of GeneCellin (BioCellChallenge). The oxygen consumption rate (OCR) of rNCMs (5 × 104/well) was measured by use of a Seahorse XF24 Extracellular Flux Analyzer (Seahorse Bioscience, North Billerica, MA).
Statistical Analysis
Data are presented as mean ± S.E. Statistical comparisons involved Student's t test or ANOVA (simple or repeated-measures) and Tukey's post-hoc analysis. A p < 0.05 was considered statistically significant.
RESULTS
Quantification of Dynamic Phosphoproteomic Profiles Responding to Pressure Overload in Mice
We used quantitative phosphoproteomics to analyze levels of phosphorylated proteins in LVs of mice to determine protein phosphorylation events within the first hour after TAB-induced pressure overload (Fig. 1A). To assure the success of TAB in acute pressure overload, we examined the activation/phosphorylation of ERK1/2 and ribosomal protein S6, known to be acutely phosphorylated in the LV after mechanical stress (
). ERK1/2 was transiently phosphorylated at 10 min after TAB, and S6 was phosphorylated up to 60 min after TAB (supplemental Fig. S1A). Early phosphorylation of ERK1/2 and S6 indicated successful induction of pressure overload. The increased ratio of heart weight to body weight (HW/BW) also indicated successful pressure overload with TAB (supplemental Fig. S1B).
Temporal responses to pressure overload were analyzed by comparing protein phosphorylation levels at various times after TAB and sham operation (Fig. 1A, Exp. Set 1). We also compared levels of cardiac phosphoproteins with chronic pressure overload (2 weeks post-TAB (Fig. 1A, Exp. Set 2). We used a post-enrichment iTRAQ labeling strategy and externally added bovine phosphoprotein as a quantification standard. The reproducibility and accuracy of this quantitation strategy was assessed by a 1:1 experiment with four identical LV protein samples (supplemental Fig. S2). With 711 quantifiable phosphopeptides (iTRAQ reporter ion intensity ≥ 25) in the 1:1 experiment, the CV and mean absolute error (MAE) were about 20 and 16%, respectively. The CV and MAE could be reduced to about 10% with reporter ion intensity > 200 for a single phosphopeptide or with number of quantifiable spectra ≥7 (supplemental Fig. S2B and D). The reproducibility and accuracy of our method were comparable to the quantification of single peptide level by iTRAQ (
). From the 1:1 experiment, we defined a 1.5 fold-change in phosphorylation level as significant difference (p < 0.05).
We could quantify 918 unique phosphopeptides from 423 proteins (supplemental Table S1), corresponding to 872 amino acid sites (723 phosphoserines, 111 phosphothreonines, 18 phosphotyrosins and 20 ambiguous sites). In all, 679 phosphopeptides were identified in acute stress, 539 phosphopeptides were identified in chronic stress, and a total of 320 phosphopeptides were consistently identified across all four time points. We determined 164, 110, 121, and 154 differentially regulated phosphorylation sites at 10, 30, and 60 min and 2 weeks, respectively, as compared with sham operation. The overall profiling of protein samples with acute pressure overload showed a rightward shift as compared with the 1:1 experimental findings (Fig. 1B), which indicates an increase in phosphorylated peptides with acute pressure overload. Notably, the analysis revealed a dynamic and transient increase in the profile, which peaked at 10 min post-TAB and was lowest at 30 min post-TAB. The profiles of phosphopeptides with chronic pressure overload also showed a rightward shift as compared with sham operation (Fig. 1C).
To assess the biological significant change, we performed Student's t test by comparing iTRAQ signal intensities of phosphopeptides between TAB and sham mice. A total of 119 and 124 phosphopeptides for acute pressure overload and hypertrophic hearts have a p < 0.05, respectively (supplemental Table S1). We further evaluated the biological reproducibility by correlation analysis of significant changed ratio. The r2 values were 0.77 and 0.79 for acute pressure overload and hypertrophic hearts (supplemental Fig. S4).
Validation of Dynamic Phosphoproteomic Profiles Responding to Pressure Overload
Change in contractility is one of the acute responses when hearts are exposed to pressure overload (
). We observed increased phosphorylation of PLN T17 and decreased phosphorylation of ATP1a1 S16 in acute pressure-overloaded but not chronic hypertrophic hearts (supplemental Table S1). These changes were further confirmed by incubation with their site-specific phospho-antibodies (Fig. 1D). Phosphorylation of the small heat shock protein Cryab on residue S59 is sufficient to mediate its anti-apoptosis function in cardiomyocytes (
Mimicking phosphorylation of alphaB-crystallin on serine-59 is necessary and sufficient to provide maximal protection of cardiac myocytes from apoptosis.
). We detected and confirmed increased phosphorylation of Cryab S59 and Hspb1 S86 (supplemental Table S1 and Fig. 1D) in acute pressure-overloaded and chronic hypertrophic mouse hearts. These results confirmed our phosphoproteomic findings and provide evidence for the protective effect of short-term pressure overload on ischemia-reperfusion injury (
Pressure-overload hypertrophy of the developing heart reveals activation of divergent gene and protein pathways in the left and right ventricular myocardium.
Am. J. Physiol. Heart Circ. Physiol.2013; 304: H697-708
). Our study revealed only 17 of our 222 differentially phosphorylated proteins among the 152 differentially expressed proteins. In addition, our study revealed only 26 differentially phosphorylated proteins of the 1072 genes previously reported to be differentially regulated in hypertrophy hearts (
). The low rate of overlap demonstrates the unique ability of the phosphoproteomic approach to identify novel phosphorylated proteins as compared with proteomic and transcriptomic approaches. By searching the Mouse Genome Informatics (MGI) database, we compared those mutated genes involved in abnormal cardiac phenotypes with our identified phosphoproteins. We found 81 phosphoproteins associated with abnormal cardiac phenotypes; 44 and 22 showed changes in phosphorylation in acute pressure-overload and chronic hypertrophic hearts, respectively (supplemental Table S3). Thus, the phosphoproteins we identified are indeed important in cardiac functions. The differential phosphorylation of previously unknown sites suggests their novel mechanisms in cardiac functions.
Temporal Specificity of Kinase Recognition Motifs in Response to Pressure Overload
To understand the types of kinases involved in the differential phosphorylation events induced by pressure overload in mice at different times, we classified the quantified phosphorylation sites into three kinase recognition motifs: acidic (e.g. CKI, CKII, GSK3), basic (e.g. protein kinases A and C, CaMKII), and proline-directed (e.g. ERK1/2, p38 MAPK, CDK5). Compared with phosphorylation sites unchanged with TAB, at 10 min, the proportion with preference for the acidic motif was reduced from 36% to 22% and that with preference for the proline-directed motif was increased from 26% to 33%; at 30 min, the proportion with preference for the basic motif was increased from 26% to 36% (supplemental Fig. S5A). These results suggest that different kinases were activated at different times after pressure overload. For example, the pSDxD sequence (pS is phosphoserine and x is any amino acid residue) was discovered as a novel motif highly phosphorylated in hypertrophic hearts (supplemental Fig. S5B). The pSDxD motif is recognized by casein kinase II (S/TxxD/E), which suggests that the casein kinase II may play important roles in cardiac hypertrophy through phosphorylation of these sites.
Phosphorylation “Hot” Sites Responding to Pressure Overload
To maintain and regulate heart contraction, cardiomyocytes require a specified subcellular structure, including the sarcoplasmic reticulum (SR), contractile fiber, and intercalated disc (ICD) (
). Thus, we also analyzed the subcellular compartments of all quantifiable proteins for those with differential phosphorylation on acute pressure overload using GoMiner (
High-Throughput GoMiner, an ‘industrial-strength’ integrative gene ontology tool for interpretation of multiple-microarray experiments, with application to studies of Common Variable Immune Deficiency (CVID).
). Among 8 major subcellular domains, the nucleus had the most phosphoproteins, but the contractile fiber showed the most changes. The average percentage of differentially phosphorylated proteins was 39.3% in these eight subcellular domains but was 60% (21/35) in contractile fiber proteins (Fig. 2A). GoMiner analysis suggest that in addition to contractile fiber, the SR membrane, ICD, Z disc, and cell-matrix junction are hot subcellular sites responding to pressure overload stress in the heart (Table I and Fig. 2B). Many SR membrane proteins, including Ryr2, sarcoplasmic/endoplasmic reticulum calcium ATPase 2 (SERCA2), PLN and histidine rich calcium binding protein are involved in regulation of calcium transient during heart contraction (
). We also observed differential phosphorylation of Ryr2 (S2808), SERCA2 (S38), and PLN (S16 and T17) (Fig. 2B). These phosphorylation sites are important for the functions of these proteins (
Identification of Ser38 as the site in cardiac sarcoplasmic reticulum Ca(2+)-ATPase that is phosphorylated by Ca2+/calmodulin-dependent protein kinase.
), so phosphorylation of SR membrane proteins and contractile proteins may be one of the immediate responses to pressure overload in cardiac myocytes.
Fig. 2Phosphorylation “hot” sites of mouse left ventricles induced by acute pressure overload.A, Classification of quantifiable (pie chart) and differentially phosphorylated (histogram) proteins. Proteins were classified into eight major subcellular compartments according to gene ontology (GO) annotation by use of GoMiner. The number in the pie chart indicates the number of phosphoproteins in each category. The y axis of histogram is the percentage of phosphoproteins with changed expression in each compartment. B, Enriched phosphorylation changes in mechanical stress-sensing sites: Z-dsic, intercalated disc and cell-matrix junction. The fold change in phosphorylation level of indicated amino acid sites at 4 times after TAB are as in the inset. C, TAB-induced phosphorylation of Xinα validated by newly generated phospho-specific antibodies. The position of molecular weight marker is showed to the right of immunoblot.
Interestingly, three other phosphorylation hot sites—ICD, Z disc, and cell-matrix junction (also named costamere in muscle cells)—are important for sensing mechanical stretch in cardiomyocytes (
). We found many differentially phosphorylated proteins located within the costameres (dystroglycan S812; tensin 1 S792, T1015; spectrin-β S2340) as well as to the ICD, including catenin protein complex (Xinα S205, S208, S295), tight junction proteins (ZO1 S111) and the gap junction protein (connexin 43 S325, S364) with acute pressure overload (Fig. 2B). Interestingly, we have recently shown that loss of Xinα leads to cardiac hypertrophy, cardiomyopathy and conduction defects in mice (
). The phosphorylation of Xinα at S205/S208 and S295 detected in the acute TAB hearts may represent another regulatory mechanism underlying stress-induced cardiac hypertrophy. To confirm our finding of the newly identified phosphorylation sites of Xinα, we generated antibodies against phospho-S205/S208 and -S295 Xinα. We confirmed the existence of these phosphorylation events and validated the phosphorylation of the Xinα protein by pressure stress in vivo (Fig. 2C). These novel findings indicate that the phosphomodulations are launched at multiple mechanical stress sensing sites of cardiomyocytes during acute pressure overload.
Distinct Phosphorylation Patterns Between Acute Pressure-Overloaded and Hypertrophic Hearts
We compared phosphoprotein regulation in hearts with acute pressure overload and chronic hypertrophy. GoMiner analysis of biological processes showed proteins participating in heart contraction, carbohydrate metabolism, cytoskeleton organization, and heart development significantly enriched with acute pressure overload but not hypertrophy (p < 0.05, Table I). Because of their roles in heart contraction, proteins located in the contractile fiber, ICD and SR membrane were also significantly changed after acute pressure overload but not hypertrophy (Table I). In contrast, proteins involved in actin-filament organization, lymphocyte activation (including T-cell activation), and cell projection organization (including neuron projection development and axonogenesis) were enriched in hypertrophic but not acutely stressed hearts (Table I). This finding is consistent with T-cell infiltration and sympathetic nerve generation being increased in hypertrophic hearts (
). Our analyses showed that different cellular pathways were activated in cardiac myocytes in response to acute or chronic pressure overload.
Phosphoproteomics Study Suggests a Novel Mediator of Cardiac Hypertrophy
The differentially phosphorylated proteins we identified by our comparative approach allowed us to study their relationship in multiple cellular pathways related to cardiac hypertrophy. A literature search (supplemental Fig. S6; details in Supplemental Table S2) allowed for summarizing these proteins and pathways associated with protein synthesis, glycogen storage, calcium homeostasis, actin polymerization, and cell survival. Our analysis revealed the phosphorylation of many well-known prohypertrophic proteins, such as mammalian target of rapamycin (
). Among those prohypertrophic proteins, PKC-δ showed increased phosphorylation in both acute pressure-overload and chronic hypertrophic hearts (supplemental Fig. S6A). Transgene activation of PKC-δ can promote cardiac hypertrophy (
). Interestingly, acute pressure overload increased the phosphorylation of dynamin-related protein 1 (DRP1 S622/isoform1, S596/isoform 2, S579/isoform 3, S491/isoform 4), a PKC-δ substrate, and interacting protein (
). However, whether DRP1 plays a role in the development of cardiac hypertrophy is unknown.
Acute Pressure Overload Induces DRP1 S622 Phosphorylation and Mitochondrial Translocation
To investigate whether DRP1 plays a role in the development of cardiac hypertrophy, we first validated our phosphoproteomic findings by determining the level of DRP1 and phospho-DRP1 (p-DRP1) (S622) after TAB in mice. To confirm the biological reproducibility, the samples we used were independent of the samples used for phosphoproteomics. Level of p-DRP1 from LV lysates was increased at 10 min after TAB as compared with sham operation (Fig. 3A and 3B), with no difference in total DRP1 protein level. Because the interaction between DRP1 and PKC-δ in the brain is increased by hypertension (
), we examined whether the interaction between DRP1 and PKC-δ in the heart changes after acute pressure overload. We found increased level of PKC-δ on co-immunoprecipitation with anti-DRP1 antibody at 10 min after TAB (Fig. 3C). Phosphorylation of DRP1 at S622 increases its mitochondrial localization (
). With mitochondrial fractionation, the mitochondrial fraction of DRP1 increased significantly at 10 min after TAB as compared with sham operation (Fig. 3D). Thus, acute pressure overload of the mouse heart may promote the interaction between DRP1 and PKC-δ and increase the amount of S622-phosphorylated DRP1 and mitochondrial translocation of DRP1.
Fig. 3Induction of S622 phosphorylation, protein kinase C δ (PKC-δ) interaction and mitochondrial localization of dynamin-related protein 1 (DRP1) by pressure overload in mouse hearts.A, Immunoblotting of phosphorylation of DRP1 S622 in sham-operated and TAB hearts. Data are mean ± S.E. * p < 0.05 versus sham operation. The position of molecular weight marker is showed to the right of immunoblot. B, Immunofluorescence images of phospho-DRP1 S622 in cryosections of LVs from sham-operated and TAB hearts. Green: phospho-DRP1 S622; Red: wheat germ agglutinin for labeling cell-surface glycan. C, Co-immunoprecipitation of PKC-δ by anti-DRP1 antibody from LV lysates of sham-operated or TAB mice. LV lysates underwent immunoprecipitation (IP) with anti-DRP1 antibody; PKC-δ and DRP1 were detected by immunoblotting (IB) with their antibodies. Data are mean ± S.E. * p < 0.05 versus sham operation, n = 3. The position of molecular weight marker is showed to the right of immunoblot. D, Immunoblotting of DRP1 in the mitochondrial fraction of sham-operated and TAB hearts. Mitochondrial-encoded cytochrome c oxidase I (MT-CO1), a mitochondrial marker, was an internal loading control. Data are mean ± S.E. * p < 0.05 versus sham operation, n = 3. The position of molecular weight marker is showed to the right of immunoblot.
Chemical Inhibitor of DRP1 Reduces LV Hypertrophy Induced by Pressure Overload
To investigate the role of DRP1 in the development of cardiac hypertrophy, we treated mice with a chemical inhibitor of DRP1, mdivi-1, before TAB. Mdivi-1 can inhibit the GTPase activity of yeast mitochondrial division dynamin and prevent mitochondrial fission in yeast and mammalian cells (
). Our echocardiography results showed significantly increased ratio of LV mass to body weight ratio in vehicle-treated mice at 1 and 2 weeks after TAB as compared with mice under the basal condition; however, in mdivi-1-treated mice, the increased ratio was significantly suppressed at 1 and 2 weeks after TAB (Fig. 4A). Similarly, HW/BW ratio was significantly lower in mdivi-1- than vehicle-treated mice (Fig. 4B, 4C). Because cardiac hypertrophy is frequently accompanied by re-expression of cardiac fetal genes (
Rapid transcriptional activation and early mRNA turnover of brain natriuretic peptide in cardiocyte hypertrophy. Evidence for brain natriuretic peptide as an “emergency” cardiac hormone against ventricular overload.
), we compared the mRNA expression of atrial natriuretic peptide, β-myosin heavy chain and skeletal α-actin in mdivi-1- and vehicle-treated mice with TAB. Compared with vehicle-treated animals, mdivi-1-treated mice showed significantly lower re-expression of these fetal genes (Fig. 4D). TAB-induced cardiac hypertrophy results from enlarged individual myocytes, so we examined cross-sectional area of myocytes in the mouse LV and found it significantly smaller in mdivi-1- than vehicle-treated mice at 2 weeks after TAB (Fig. 4E, 4F). Thus, inhibition of DRP1 with mdivi-1 may ameliorate TAB-induced cardiac hypertrophy.
Fig. 4Inhibition of DRP-1 alleviates TAB-induced hypertrophy responses in mice. The ratio of (A) LV mass to body weight (LVM/BW) and (B) heart weight to body weight (HW/BW) for TAB mice treated with DRP1 inhibitor, mdivi-1, or vehicle. Mdivi-1 (25 mg/Kg) was given by intraperitoneal injection. The LV mass was measured by echocardiography. The heart weight was measured 2 weeks after TAB. * p < 0.05 versus basal level; # p < 0.05 versus vehicle group at indicated TAB time; n = 6 for vehicle; n = 7 for mdivi-1. C, Representative sections of vehicle- and mitochondrial division inhibitor 1 (mdivi-1)–treated hearts. Cross sections of hearts were stained with Masson's trichrome 2 weeks after TAB. D, Real-time quantitative PCR of mRNA level in vehicle- and mdivi-1-treated hearts with TAB-induced re-expression of hypertrophic marker genes 2 weeks after TAB. Marker genes are atrial natriuretic peptide (ANP), β-myosin heavy chain (βMHC), and skeletal α-actin (ACTA1). Data are mean ± S.E. # p < 0.05 versus vehicle; n = 3. E, Cross-sectional area of LV mouse cardiomyocytes 2 weeks after TAB. Sections were stained with WGA to show cell-surface glycan. The surface area of cardiomyocytes was calculated and comparative results are shown in (F). Data are mean ± S.E. # p < 0.05 versus vehicle; n = 6 for vehicle; n = 7 for mdivi-1.
Inhibition of DRP1 Prevents PE-Induced Hypertrophy in Rat Cardiomyocytes
To determine whether the reduced cardiac hypertrophy observed in mdivi-1-treated mice with TAB was attributable to the direct effect on cardiomyocytes, we isolated and cultured rNCMs and tested their hypertrophic responses to PE in vitro. Similar to observations in TAB mouse hearts, level of p-DRP1 (S622) was significantly increased in rNCMs after 1-hr PE treatment (Fig. 5A). PE significantly increased the surface area of rNCM (by 32%, p < 0.001). However, pretreatment with mdivi-1 blunted the PE-induced increase in cell-surface area (Fig. 5B). In addition to pharmacological inhibition, we also downregulated DRP1 activity by transfecting rNCMs with a dominant-negative construct, DRP1K38A (
). Similar to results with mdivi-1, the PE-induced increase in cell-surface area was blunted in DRP1K38A-transfected rNCMs (Fig. 5C). Our in vivo and in vitro studies suggest that DRP1 activity is important in the development of cardiac hypertrophy, and the DRP1 activity inhibitor mdivi-1 may help ameliorate the development of cardiac hypertrophy.
Fig. 5DRP1 inhibition prevents phenylephrine (PE)-induced hypertrophic growth and oxygen consumption of rat neonatal cardiomyocytes (rNCMs).A, Phosphorylation of DRP1 S622 in rNCMs after 1 h of stimulation with 10 μm PE. Data are mean ± S.E. * p < 0.05 versus vehicle; n = 3 for each group. The position of molecular weight marker is showed to the right of immunoblot. B, DRP1 inhibitor, mdivi-1, prevents PE-induced increase in cell surface area of rNCM. rNCMs were pretreated with 50 μm mdivi-1 or vehicle for 1 h, then treated with 10 μm PE for 24 h. rNCMs were stained with anti-α-actinin antibody (red) and DAPI (blue, for nuclei). At least 100 cells were counted for each group. Data are mean ± S.E. * p < 0.05 versus control in vehicle. C, Transfection of dominant-negative DRP1 prevents PE-induced increase in cell-surface area of rNCMs. HA-tagged dominant-negative mutant of DRP1 (HA-DRP1-K38A) was transfected into rNCMs for 2 days, then rNCMs were treated with PE for 24 h. rNCMs were incubated with anti-HA and anti-α-actinin antibodies. Untransfected groups were cells negative for HA. At least 30 cells were counted for each group. Data are mean ± S.E. * p < 0.05 versus control in vehicle. D, DRP1 inhibitor, mdivi-1, prevents PE-induced increase in oxygen consumption of rNCM. 50 μm Mdivi-1 or vehicle was added to rNCMs 20 min before measurement of real-time oxygen consumption rate (OCR). 10 μm PE and 2 μm oligomycin (oli) were automatically injected into medium at the indicated time. Oligomycin is a mitochondrial ATP synthase inhibitor and is used to estimate the amount of consumed oxygen for ATP synthesis. OCR at each time was normalized to the time just before PE addition. Data are mean ± S.E. # p < 0.05 versus control in vehicle.
). Because DRP1 mediates mitochondrial fission, we examined the impact of DRP1 inhibition on the OCR of rNCMs. PE increases the cellular OCR in hypertrophic rNCMs (
). Mdivi-1 inhibition of DRP1 caused a transient but recoverable reduction in OCR immediately after mdivi-1 treatment (Fig. 5D). PE significantly increased the OCR of rNCM but mdivi-1 pretreatment prevented this increase. DRP1-mediated mitochondrial fission may be important for metabolism changes with cardiac hypertrophy.
DISCUSSION
Quantitative Phosphoproteomic Study Reveals the Mechanism of Pressure-Overload-Induced Cardiac Hypertrophy
In this study, we used a MS-based quantitative platform to analyze temporal changes in site-specific protein phosphorylation from murine hearts subjected to TAB-induced pressure overload. This systematic and unbiased approach allowed for revealing multiple mechanisms of responses to pressure-overload stress. For example, we found increased phosphorylation of SERCA2 at S38 and PLN at S16 and T17 at 10 min after TAB. Phosphorylation of these three sites increases the activity of SERCA2 and leads to increased SR calcium reuptake and heart contraction (
Identification of Ser38 as the site in cardiac sarcoplasmic reticulum Ca(2+)-ATPase that is phosphorylated by Ca2+/calmodulin-dependent protein kinase.
). In contrast, phosphorylation of Ryr2 at S2808 was decreased 60 min after TAB. The phosphorylation of Ryr2 at S2808 increases calcium leakage from the SR in failing hearts (
). Our findings suggest that acute TAB could promote SERCA2 calcium reuptake and prevent SR calcium leakage by phosphorylation of regulatory sites of SERCA2, PLN and Ryr2. Interestingly, these changes were not observed in 2-week pressure-overloaded hypertrophic hearts. In addition, we found phosphorylation of α-tropomyosin at T282 and S283 increased in acute pressure-overloaded hearts but decreased in 2-week pressure-overloaded hearts. Dephosphorylation of α-tropomyosin is related to suppression of myocardial sarcomeric tension and ATPase activity in the hypertrophic heart (
Here we found phosphorylation changes in heart-contraction, SR-membrane and heart-development proteins with acute pressure-overloaded but not chronic hypertrophic hearts. The absence of these phosphorylation changes in the chronic hypertrophic heart may be because of reduced myocardial wall stress with increased wall thickness of the hypertrophic heart. Another possibility is the additional contraction force provided by the newly synthesized myocardium in the hypertrophic heart. However, we found increased phosphorylation of Cryab S59 and Hspb1 S82 in both acute pressure-overloaded and chronic hypertrophic hearts. The differential phosphorylation events in chronic hypertrophic hearts may be attributed to protein level change, because we and others found the protein level of Cryab up-regulated during hypertrophy (
). Increased DRP1 function may be a general response to pressure stress for the brain and heart. In addition, we found that PE, an α-adrenergic receptor (αAR) agonist, could induce DRP1 S622 phosphorylation in cardiomyocytes. Stimulation of αAR by PE causes PKC-δ activation in rNCMs (
Differential activation of protein kinase C isoforms by endothelin-1 and phenylephrine and subsequent stimulation of p42 and p44 mitogen-activated protein kinases in ventricular myocytes cultured from neonatal rat hearts.
). Therefore, TAB-induced pressure overload might trigger the αAR/PKC-δ signal pathway to activate DRP1. Although the activity of Cdk/cyclin B in adult cardiomyocytes is low (
Cell cycle regulatory molecules (cyclins, cyclin-dependent kinases and cyclin-dependent kinase inhibitors) and the cardiovascular system; potential targets for therapy?.
Long-term inhibition of Rho-kinase suppresses angiotensin II-induced cardiovascular hypertrophy in rats in vivo: effect on endothelial NAD(P)H oxidase system.
) all promote cardiac hypertrophy. DRP1 may be an important convergent point of signaling pathways for cardiac hypertrophy. Determining whether regulation of DRP1 is important for the prohypertrophic activity of these proteins would be of interest.
Our results show that pharmacological and genetic inhibition of DRP1 could prevent cardiac hypertrophy in mice and in rNCMs. The mechanism of this anti-hypertrophic effect may be related to cytosolic calcium and reactive oxygen species (ROS). Recently, DRP1 was shown to mediate the oxygen-induced ductus arteriosus smooth muscle cell (DASMC) contraction and proliferation by increasing cytosolic calcium and ROS (
), TAB may lead to DRP1 activation, which increases intracellular ROS and calcium levels and thus induces cardiac hypertrophy. Another possible mechanism for the involvement of DRP1 in cardiac hypertrophy is the regulation of mitochondrial metabolism. Mdivi-1 inhibits oxygen-induced OCR increase in DASMCs (
). We also showed that mdivi-1 prevented oxygen consumption induced by PE in rNCMs. PE may affect mitochondrial metabolic activity by regulating DRP1 and mitochondrial fission.
Mitochondrial Fission and Fusion Dynamics and Clinical Perspectives for Cardiac Hypertrophy
Our study suggests that regulation of mitochondrial fission proteins may be critical for TAB-induced cardiac hypertrophy in mice. In addition to fission, mitochondrial morphologic features are also controlled by a mitochondrial fusion protein, optic atrophy 1 (Opa1) (
). The greater hypertrophy of mice with defective mitochondrial fusion protein is consistent with our finding of less hypertrophy in mice with inhibited mitochondrial fission protein. Mitochondrial fusion proteins may suppress TAB-induced cardiac hypertrophy, but mitochondrial fission proteins may promote it.
Recently, mice carrying a single mutation in DRP1 showed dilated cardiomyopathy and congestive heart failure (CHF) (
). Myocyte hypertrophy was observed in the CHF stage of DRP1 heterozygous mutant mice. This result is in contrast with our finding that inhibition of DRP1 prevented the development of cardiac hypertrophy. One possibility for this discrepancy is that genetic manipulation in mutant DRP1+/− mice may produce developmental defects that differ from that with TAB induction. Inhibition of DRP1 has beneficial effects in several disease models (
). Our finding that inhibition of DRP1 prevented the development of cardiac hypertrophy is in agreement with these previous studies.
Traditional drug targets of cardiac hypertrophy are cell-surface receptors or ion channels, but the newer drugs target the intracellular signaling pathways, including kinases, phosphatases, and transcription modulators (
). Our phosphoproteomic study revealed novel signaling molecules that could be potential therapeutic targets for intervention in cardiomyopathy. Here, we demonstrated that mdivi-1, a small-molecule inhibitor of DRP1, could suppress pressure-overload-induced cardiac hypertrophy. Thus, DRP1 may be a new therapeutic target in cardiac hypertrophy.
To our knowledge, this is the first phosphoproteomics study to reveal the early response to acute pressure overload in the heart in vivo. Most of these dynamically activated or repressed phosphorylation sites were previously unrecognized to be regulated by pressure overload. This phosphoproteomic approach can reveal the early signal pathways in mechanosensing, metabolism, heart contraction and cardiac hypertrophy during pressure overload of the heart. We also demonstrated that inhibition of the mitochondrial fission protein DRP1 could prevent cardiac hypertrophy in mice and in rat cardiomyocytes. This study broadens our understanding of signaling responses of the heart to pressure stress. As well, interfering in DRP1 function and mitochondrial fission may help control the development of cardiac hypertrophy.
Acknowledgments
We thank the mass spectrometry center of the Institute of Chemistry, Academia Sinica.
Pressure-overload hypertrophy of the developing heart reveals activation of divergent gene and protein pathways in the left and right ventricular myocardium.
Am. J. Physiol. Heart Circ. Physiol.2013; 304: H697-708
A multiplexed quantitative strategy for membrane proteomics: opportunities for mining therapeutic targets for autosomal dominant polycystic kidney disease.
High-Throughput GoMiner, an ‘industrial-strength’ integrative gene ontology tool for interpretation of multiple-microarray experiments, with application to studies of Common Variable Immune Deficiency (CVID).
Mimicking phosphorylation of alphaB-crystallin on serine-59 is necessary and sufficient to provide maximal protection of cardiac myocytes from apoptosis.
Identification of Ser38 as the site in cardiac sarcoplasmic reticulum Ca(2+)-ATPase that is phosphorylated by Ca2+/calmodulin-dependent protein kinase.
Rapid transcriptional activation and early mRNA turnover of brain natriuretic peptide in cardiocyte hypertrophy. Evidence for brain natriuretic peptide as an “emergency” cardiac hormone against ventricular overload.
Differential activation of protein kinase C isoforms by endothelin-1 and phenylephrine and subsequent stimulation of p42 and p44 mitogen-activated protein kinases in ventricular myocytes cultured from neonatal rat hearts.
Cell cycle regulatory molecules (cyclins, cyclin-dependent kinases and cyclin-dependent kinase inhibitors) and the cardiovascular system; potential targets for therapy?.
Long-term inhibition of Rho-kinase suppresses angiotensin II-induced cardiovascular hypertrophy in rats in vivo: effect on endothelial NAD(P)H oxidase system.
34591-6&id=gr1.jpg){kind=link}
34591-6&id=gr2.jpg){kind=link}
34591-6&id=gr3.jpg){kind=link}
34591-6&id=gr4.jpg){kind=link}
34591-6&id=gr5.jpg){kind=link}