Next Article in Journal
Hypothermia Induced by Oxcarbazepine after Transient Forebrain Ischemia Exerts Therapeutic Neuroprotection through Transient Receptor Potential Vanilloid Type 1 and 4 in Gerbils
Previous Article in Journal
Green Efficient One-Pot Synthesis and Separation of Nitrones in Water Assisted by a Self-Assembled Nanoreactor
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 1
10.3390/ijms23010235
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Omics Approach Profiling Metabolic Remodeling in Early Systolic Dysfunction and in Overt Systolic Heart Failure

by
Antoine H. Chaanine
1,2,*,
LeeAnn Higgins
3,
Todd Markowski
3,
Jarrod Harman
4,
Maureen Kachman
5,
Charles Burant
6,
L. Gabriel Navar
2,
David Busija
7,8 and
Patrice Delafontaine
1,2,7
1
Department of Medicine, Tulane University, New Orleans, LA 70112, USA
2
Department of Physiology, Tulane University, New Orleans, LA 70112, USA
3
Department of Biochemistry, Molecular Biology, and Biophysics, University of Minnesota Medical School, Minneapolis, MN 55455, USA
4
Department of Ophthalmology, Louisiana State University Health Sciences Center, New Orleans, LA 70112, USA
5
Michigan Regional Comprehensive Metabolomics Resource Core, University of Michigan, Ann Arbor, MI 48109, USA
6
Department of Internal Medicine, University of Michigan Medical School, Ann Arbor, MI 48109, USA
7
Department of Pharmacology, Tulane University, New Orleans, LA 70112, USA
8
Tulane Brain Institute, Tulane University, New Orleans, LA 70112, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(1), 235; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23010235
Submission received: 7 November 2021 / Revised: 6 December 2021 / Accepted: 22 December 2021 / Published: 26 December 2021
(This article belongs to the Section Molecular Biology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Metabolic remodeling plays an important role in the pathophysiology of heart failure (HF). We sought to characterize metabolic remodeling and implicated signaling pathways in two rat models of early systolic dysfunction (MOD), and overt systolic HF (SHF). Tandem mass tag-labeled shotgun proteomics, phospho-(p)-proteomics, and non-targeted metabolomics analyses were performed in left ventricular myocardium tissue from Sham, MOD, and SHF using liquid chromatography–mass spectrometry, n = 3 biological samples per group. Mitochondrial proteins were predominantly down-regulated in MOD (125) and SHF (328) vs. Sham. Of these, 82% (103/125) and 66% (218/328) were involved in metabolism and respiration. Oxidative phosphorylation, mitochondrial fatty acid β-oxidation, Krebs cycle, branched-chain amino acids, and amino acid (glutamine and tryptophan) degradation were highly enriched metabolic pathways that decreased in SHF > MOD. Glycogen and glucose degradation increased predominantly in MOD, whereas glycolysis and pyruvate metabolism decreased predominantly in SHF. PKA signaling at the endoplasmic reticulum–mt interface was attenuated in MOD, whereas overall PKA and AMPK cellular signaling were attenuated in SHF vs. Sham. In conclusion, metabolic remodeling plays an important role in myocardial remodeling. PKA and AMPK signaling crosstalk governs metabolic remodeling in progression to SHF.

1. Introduction

The heart is an organ with high-energy demands, which has the flexibility to utilize almost any type of substrate for adenosine triphosphate (ATP) production [1]. It cycles nearly 6 kg of ATP per day, produced in the mitochondria by oxidative phosphorylation (OXPHOS). Ninety-five percent of the generated ATP is consumed by myosin ATPase, the sodium/potassium pump, and the sarco/endoplasmic reticulum calcium ATPase (Serca2a or Atp2a2), which is necessary for the contraction and maintenance of ion homeostasis during the contraction and relaxation phases of the cardiac cycle, respectively [2]. Mitochondria are abundant in the heart and constitute about 30% of the cardiomyocyte volume [3,4]. In addition to their role in OXPHOS and ATP production, mitochondria play a role in myocardial calcium and redox signaling, ion homeostasis, and metabolism. This multiplicity of mitochondrial (mt) function is tightly coordinated and balanced to maintain cell survival [5].
Perturbations in mt function and metabolic remodeling are evident in heart failure (HF) [6,7,8]. Changes in mt morphology and dynamics take place early in pathological hypertrophy and worsen during transition to advanced systolic HF (SHF) [9,10]. Moreover, fatty acid (FA) utilization, mt FA β-oxidation, and OXPHOS are impaired early in pathological hypertrophy and gradually decline with the progression of myocardial remodeling and systolic dysfunction [6,11]. The heart initially compensates through increases in glucose uptake and glycolysis, but this is eventually impaired upon progression to SHF [6,11]. Less is known about when changes in amino acid (AA) and nucleotide metabolism and branched-chain amino acid (BCAA) catabolism occur and their contribution to myocardial remodeling. Derangements in cardiac metabolism are debated as to whether they are the cause or a marker of progression to SHF. Evidence suggests that the former is true. Changes in myocardial signaling as a result of altered workload and biomechanical and neurohumoral stress are evident in the earliest phases of cardiac stress and drive the initiation of cardiac remodeling and its progression toward SHF [12]. Kinases and phosphatases, such as protein kinase A (PKA), calcium/calmodulin-dependent protein kinase II (subunits delta (CaMKIIδ), and gamma (CaMKIIγ)), and calcineurin that drive cardiac remodeling and altered calcium signaling [13,14,15] are also implicated in metabolic reprogramming [16,17,18,19]. Moreover, classic metabolites and metabolic intermediates, summarized elsewhere [1], have been recognized as modulators of signal transduction and ion channels/calcium homeostasis, suggesting that metabolic remodeling is an important player in myocardial remodeling and the progression to overt SHF.
What remains unknown is: (1) a global understanding of metabolic disturbances, primarily changes in mitochondrial proteome and phosphoproteome in early systolic dysfunction and in transition to overt SHF, and (2) the time point at which derangements in signal transduction pathways, primarily PKA and adenosine monophosphate-activated protein kinase (AMPK) signaling modulating metabolic remodeling occur, and which of those predominantly drive metabolic remodeling in the transition to overt SHF. To gain a global understanding of metabolic remodeling at the proteome and metabolome level, and implicated signal transduction pathways, we performed tandem mass tag (TMT)-labeled shotgun proteomics, untargeted metabolomics, and shotgun phospho (p)-proteomics analysis of left ventricular (LV) myocardium tissue in two rat models of induced pressure overload (PO): (1) moderate remodeling and early systolic dysfunction (MOD), and (2) overt SHF.

2. Results

The echocardiographic parameters of the studied animals are presented in Supplementary Table S1. At week 3 post ascending aortic banding (AAB), there was evidence of concentric hypertrophy in the MOD and SHF groups with increases in LV septal and posterior wall thicknesses, and LV ejection fraction (LVEF) compared with Sham. At week 8 post-AAB, the MOD and SHF phenotypes continued to show the same degree of concentric hypertrophy, but developed increases in LV end-diastolic (LVEDV) and end-systolic volumes (LVESV) and a decrease in LVEF compared with week 3 post-AAB. Moreover, the SHF phenotype showed significantly higher LVEDV and LVESV, but lower LVEF than the MOD phenotype. The heart to body weight (BW) ratio increased in SHF > MOD vs. Sham. The LV weight to BW ratio increased equally in both MOD and SHF vs. Sham; however, right ventricular weight to BW ratio increased in SHF compared with MOD and Sham.

2.1. Visualization of the MOD and SHF Proteomic Datasets

Heat maps and PCA plots (Supplementary Figure S1A,B) show the differential log2 fold change and the variance in biological samples, respectively, for the total identified proteins in Sham and MOD and Sham and SHF proteomic runs. The PCA plots show significant separation of the Sham and MOD, and Sham and SHF biological samples, suggesting a phenotypic change in MOD and SHF groups compared with Sham. Heat maps, volcano plots, and PCA plots (Figure 1A,B and Supplementary Figure S1C), respectively, show the differential log2 fold change, the log2 fold change in group means, and the variance in biological samples, respectively, for the identified proteins that changed in MOD and SHF vs. Sham. A total of 227 and 722 proteins were down-regulated in MOD and SHF, respectively, vs. Sham. Of these, 55% (125/227) and 45% (328/722) (Figure 1C), respectively, were mt proteins (intersection of green with blue diagrams and yellow with red diagrams, respectively). A total of 318 and 783 proteins were up-regulated in MOD and SHF vs. Sham, respectively (Figure 1C). Of these, 3% (9/318) and 1.3% (10/783), respectively, were mt proteins. The corresponding analysis in Scaffold is presented in Supplementary Figure S1E–G. Heat maps, volcano plots, and PCA plots of the identified mt proteins that changed in MOD and SHF vs. Sham, analyzed in Qlucore, are presented in Supplementary Figure S2A–C. The corresponding analysis in Scaffold, including protein–protein interaction network by confidence, is presented in Supplementary Figure S2D–F. Transitioning from MOD to SHF was associated with an increase in the number of down-regulated mt proteins vs. Sham (32% (125/393) vs. 67% (328/49)). One hundred and fourteen mt proteins were downregulated in MOD and SHF vs. Sham. An additional 11 mt proteins were down-regulated in MOD vs. Sham; whereas an additional 214 (160 + 54) mt proteins were down-regulated in SHF vs. Sham (Figure 1C). A total of 160 of these were identified in the MOD proteomic run (but did not reach statistical significance), and 54 were not. A similar trend was seen in the dataset analyzed by Scaffold. Volcano plots (Supplementary Figure S2G,H), show the log2 fold change in group means for the 114 and 160 mt proteins in MOD and SHF vs. Sham that were discussed above.
Heat maps, volcano plots, and PCA plots (Figure 1D,E as well as Supplementary Figure S1D) show the differential log2 fold change and the log2 fold change in group means, as well as the variance in biological samples, respectively, for the identified proteins related to metabolism and respiration that changed in MOD and SHF vs. Sham. Transitioning from MOD to SHF was associated with an increase in the number of metabolism- and respiration-related proteins that were predominantly down-regulated, 125 and 306, in MOD and SHF, respectively, vs. Sham (Figure 1F); 82% (103/125) and 71% (218/306) of these were mt proteins. Corresponding analysis in Scaffold is presented in Supplementary Figure S3A–D. Volcano plots (Supplementary Figure S3E,F) show the log2 fold change in group means for the 113 (95 + 18) commonly identified proteins related to metabolism and respiration that were down-regulated in MOD and SHF vs. Sham, and for the 100 commonly identified mt proteins involved in metabolism and respiration that were down-regulated in SHF only. Data visualization of the identified mt proteins, not involved in metabolism or respiration, that changed in MOD and SHF vs. Sham are presented in Supplementary Figure S4A–C (analyzed in Qlucore) and D,E (analyzed in Scaffold).
Heat maps (Figure 2A and Supplementary Figure S5) show the canonical metabolic pathways enriched by z-score and p-value in the MOD and SHF vs. Sham. OXPHOS followed by mt FA β-oxidation were the highest enriched and inhibited metabolic pathways (by p-value and z-score, respectively) in SHF > MOD vs. Sham. BCAA catabolism, TCA cycle, glutamine, and tryptophan degradation were also highly enriched and inhibited, by p-value and z-score, respectively, in SHF > MOD vs. Sham. Glycogen degradation was mainly activated in MOD, whereas glucose biosynthesis and UDP-N-acetyl-D-galactosamine biosynthesis were activated in MOD > SHF vs. Sham. Glycolysis was enriched in SHF > MOD and inhibited in MOD and SHF; however, the z-score was below the cutoff value of 2. Acetate conversion to acetyl-coenzyme A (CoA) and acetyl-CoA biosynthesis were enriched and inhibited in SHF vs. MOD and Sham. Purine, guanosine, and adenosine nucleotide degradation were enriched and activated in SHF > MOD; however, the z-score was below the cutoff value of 2. Ketolysis was enriched in SHF > MOD, and was inhibited in both MOD and SHF, but the z-score was below the cutoff value of 2.
Of the upstream regulators (Figure 2A) that were highly enriched in MOD and SHF vs. Sham, the peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PPARGC1A) and the peroxisome proliferator-activated receptor alpha (PPARA) were inhibited in SHF > MOD vs. Sham, whereas the mitogen-activated protein kinase kinase 4 (MAP4K4), which is upstream of Jun N-terminal kinase (JNK), and the transforming G growth factor beta 1 (TGFB1) were activated in SHF > MOD vs. Sham. The PPARGC1A and carnitine palmitoyltransferase 1B (CPT1B) gene networks in MOD and SHF vs. Sham are presented in Supplementary Figure S6A–D, respectively. PGC-1α expression (Figure 2B) was decreased in MOD vs. Sham and further decreased in SHF vs. MOD and Sham. The expression of the mt trifunctional enzymes (HADHA and HADHB) and the mt medium-chain acyl-CoA dehydrogenase (MCAD) proteins, involved in mt FA β-oxidation, were decreased in MOD and SHF vs. Sham (Figure 2B). HADHB expression was further decreased in SHF vs. MOD and Sham. The OXPHOS network of the electron transport chain (ETC) complexes I-IV and ATP synthase (complex V) (Figure 3A) show that the ETC complexes II, IV and ATP synthase were decreased/inhibited in MOD vs. Sham, whereas all the ETC complexes I-IV and ATP synthase were down-regulated/inhibited in SHF vs. Sham. ETC complexes I and III were somewhat affected in MOD vs. Sham. ETC complexes III, IV, and ATP synthase expression were decreased in MOD vs. Sham and further decreased in SHF vs. MOD and Sham (Figure 3B). ETC complex II expression was equally decreased in MOD and SHF, whereas ETC complex I expression was decreased in SHF only (Figure 3B).

2.2. Visualization of the Untargeted Metabolomics Dataset in Sham, MOD and SHF Groups

Heat maps, volcano plots, and PCA plots (Figure 4A–C) show the differential log2 fold change, the log2 fold change in group means, and the variance in biological samples, respectively, for the identified metabolites that changed in MOD and SHF vs. Sham. The PCA plots indicate that the MOD and SHF biological samples share a similar unidirectional change for the metabolites that changed in MOD vs. Sham, whereas the biological samples in the MOD group are centered midway between the Sham and SHF groups for the metabolites that changed in SHF vs. Sham. These data suggest a phenotypic disease progression upon transitioning from MOD to an advanced SHF phenotype, similar to that shown by the proteomic datasets. A total of 45 metabolites changed in MOD (blue) and 173 in SHF (red) vs. Sham and 56 (green) changed in SHF vs. MOD (Figure 4D). Eighteen common metabolites changed in MOD and SHF vs. Sham (intersection of red and blue diagrams). Seven of these also changed in SHF vs. MOD (intersection of all three diagrams). Of the 11 commonly identified metabolites that changed in MOD and SHF vs. Sham, there was a trend toward the further down-regulation of 2′,3′ cyclic AMP (or 3′,5′ cyclic AMP) and pantetheine 4′-phosphate in SHF compared to MOD (Figure 4E). Heat maps and PCA plots of the 18 commonly identified metabolites that changed in MOD and SHF vs. Sham and of the 56 metabolites that changed in SHF vs. MOD and Sham are presented in Figure 4F,G, respectively. Figure 5A,B are a heat map and PCA plot of the fatty acyl-carnitines that changed in abundance in MOD and SHF vs. Sham. The very long- and long-chain acyl-carnitines were equally down-regulated in MOD and SHF vs. Sham, whereas the medium- and short-chain acyl-carnitines decreased in abundance in SHF vs. Sham and down-trended in MOD vs. Sham.
Heat maps, generated by Ingenuity Pathway (IPA) comparison analysis (Figure 5C), show the metabolic pathways and upstream regulators that were enriched by p-value and z-score, respectively, in MOD and SHF vs. Sham. Metabolic pathways related to CMP-N-acetylneuraminate biosynthesis, glycolysis, taurine biosynthesis, valine degradation, glutamine, and sucrose degradation were enriched in MOD vs. Sham. However, CMP-N-acetylneuraminate biosynthesis, glycolysis, isoleucine, leucine and valine degradation, UDP-N- acetyl-D-glucosamine biosynthesis, glutamine degradation, and citrulline metabolism metabolic pathways were enriched in SHF vs. Sham. Bar graphs (Supplementary Figure S7) show the metabolic pathways that were expected to be activated (orange) or inhibited (blue) in MOD and SHF vs. Sham and in SHF vs. MOD. A z-score could not be estimated due to the limited number of identified metabolites per pathway or lack of information on the expected change in certain metabolites per metabolic pathway in IPA. Transitioning from MOD to SHF was associated with the down-regulation of glycolysis, BCAA catabolism, and gluconeogenesis, as well as the up-regulation of taurine biosynthesis, CMP-N-acetylneuraminate biosynthesis, and methionine degradation. Enrichment pathway analysis (Supplementary Figure S8) shows the enrichment ratio for the identified metabolites that changed in MOD and SHF vs. Sham. Of the upstream regulators that were highly enriched, calcium and calcium-independent phospholipase A2 (PNPLA8) were decreased/inhibited in SHF vs. Sham only (Figure 5C). The calcium mechanistic network (Figure 5D) shows the associated metabolites involved in calcium cycling (2′,3′ cyclic AMP) and calcium homeostasis regulation (eicosanoids) that were down-regulated in SHF only vs. Sham.

2.3. Visualization of the MOD and SHF p-Proteomic Datasets

A total of 817 and 961 proteins, as well as 2034 and 2524 p-sites were identified in MOD and SHF p-proteomic runs, respectively (Supplementary Figure S9A). Venn diagrams (Supplementary Figure S9B) show the number of p-sites that changed in MOD and SHF vs. Sham, analyzed in Scaffold PTM and Qlucore. Heat maps, volcano plots, and PCA plots (Figure 6A–C) show the differential log2 fold change, the log2 fold change in group means, and the variance in biological samples, respectively, for the identified p-sites that changed in MOD and SHF vs. Sham. One of the Sham biological samples in the MOD p-proteomic run was an outlier and was eliminated because it skewed the data (heat map and PCA plot: Supplementary Figure S10A). Corresponding analysis in Scaffold PTM is presented in Supplementary Figure S10B,C. Heat maps, generated by IPA comparison analysis (Figure 7A and Supplementary Figure S11), show signaling pathways that were enriched by p-value and z-score in MOD and SHF vs. Sham. PKA and AMPK signaling, two important signaling pathways regulating metabolism, were enriched in SHF > MOD vs. Sham, but were inactivated in SHF only vs. Sham. These results were validated by Western blotting (Figure 7B). Phosphorylation of AMPK at threonine 172 (p-T172-AMPK) increased in MOD vs. Sham, whereas p-S485/491-AMPK increased in MOD and SHF vs. Sham. Phospholamban phosphorylation (p-S16-PLN) decreased and protein phosphatase 1a (PP1a) expression increased in MOD and SHF vs. Sham. PP1a expression up-trended further in SHF vs. MOD. Transitioning from MOD to SHF was associated with a decrease in the phosphorylation of downstream PKA targets at different cellular subcompartments. Glycogen synthase kinase-3 alpha (p-S21-GSK3α), glycogen synthase kinase-3 beta (p-S9-GSK3β), troponin I, and cardiac muscle (p-S23/24-TNNI3) were all decreased in SHF vs. MOD and Sham. p-T286-CaMKII down-trended in SHF vs. Sham.
In the next section, we will address changes in specific signaling pathways separately. Analysis of the shotgun proteomic datasets, including metabolic pathways and signaling pathways that changed in MOD and SHF vs. Sham, are presented in Supplementary Excel file S1. Analysis of the untargeted metabolomics datasets is presented in Supplementary Excel file S2. Changes in p-sites and the corresponding signaling pathways are presented in Supplementary Excel file S3.

2.4. FA Metabolism and mt FA β-Oxidation and Triacylglycerol (TAG) Metabolism

The FA receptor, Slc27a1, and the FA-binding protein 3 (Fabp3), which are involved in the FA uptake and intracellular transfer of long-chain FA and their acyl-CoA esters in the heart, respectively, were down-regulated in MOD and SHF vs. Sham (Excel file S1). Moreover, the enzyme long-chain-FA-CoA ligase 1 (Acsl1), which is involved in the activation of long-chain-FA for both the synthesis of cellular lipids and degradation via β-oxidation, was down-regulated in MOD and SHF. The Fabp3 and Acsl1 were down-regulated in SHF > MOD vs. Sham (Excel file S1). The rate-limiting enzyme, carnitine O-palmitoyltransferase 1, muscle isoform (Cpt1b), which is involved in the conversion of long-chain fatty-acyl-CoA to long-chain fatty-acyl-carnitine and their transport across the outer mitochondrial membrane, was decreased in SHF vs. Sham (Excel file S1). CPT1B was predicted to be activated in IPA in MOD and SHF, SHF > MOD, due to the number of downstream molecules detected in the dataset. However, proteomics determined that Cpt1b was down-regulated in SHF vs. Sham (log2 fold change: −0.37, p = 0.002). Moreover, the enzyme malonyl-CoA decardoxylase (Mlycd), which is involved in the conversion of malonyl-CoA to acetyl-CoA and hence the activation of Cpt1b, was down-regulated in MOD and SHF, SHF > MOD, vs. Sham (Excel file S1). Malonyl-CoA is a potent allosteric inhibitor of Cpt1b and hence FA metabolism, while acetyl-CoA carboxylase, beta isoform (Acacb), the enzyme converting acetyl-CoA to malonyl-CoA, was down-regulated in SHF only (Excel file 1). The specific acyl-CoA dehydrogenases, short-chain (Acadsb), medium-chain (Acadm), long-chain (Acadl), and family members 8 (Acad8) and 10 (Acad10), were downregulated in SHF vs. Sham (Excel file S1). This is consistent with the metabolomic data showing all forms of acyl-carnitines to have decreased in abundance in SHF vs. Sham. PGC-1α expression and signaling decreased in SHF > MOD vs. Sham. Acetyl-CoA abundance down-trended in MOD and SHF vs. Sham (Supplementary Figure S12A; acetyl-CoA relative expression decreased in one Sham biological sample and statistical significance could not be attained unless that sample was excluded). Glycerol-3-phosphate dehydrogenase (Gpd1), an enzyme needed for glycerol-3-phospate generation and de novo TAG biosynthesis, was down-regulated in SHF (Excel file S1). Metabolomics data showed that TAG abundance decreased in MOD and SHF vs. Sham. The metabolites TG(38:4) and TG(36:1) were down-regulated in SHF vs. Sham, whereas TG(44:2) was up-regulated, and DG(22:3) was down-regulated in MOD vs. Sham (Excel file S2).

2.5. Glucose, Pyruvate Metabolism, and BCAA Catabolism

The glucose uptake receptor, Glut4 (Slc2a4), the predominant isoform in the adult heart than Glut1, was equally down-regulated in MOD and SHF (Excel file S1). Glut4 activity is also enhanced through its phosphorylation by AMPK, which also triggers its translocation to the plasma membrane. AMPK kinase signaling was attenuated in SHF, suggesting a further decrease in glucose uptake in SHF. The metabolites 2-deoxy-D-glucose-6-phosphate and glyceraldehyde-3-phosphate were down-regulated in SHF vs. Sham and down-trended in MOD vs. Sham (Excel file S2). Phosphoglucomutase (Pgm1) was up-regulated in MOD vs. Sham and was mildly increased in SHF vs. Sham. Glycogen synthase 1 (Gys1) and glycogen phosphorylase (Pygm) are two nodal enzymes that regulate glycogen homeostasis in the heart by promoting glycogen synthesis and glycogen degradation, respectively. Pygm was down-regulated in SHF vs. Sham only (Excel file S1). It has been shown that Pygm is activated through its phosphorylation, either by PKA or by the calcium-activated phosphorylase kinase [20], while the serine/threonine-protein phosphatase 1 (PP1) dephosphorylates Gys1 and Pygm, promoting Gys1 activation and Pygm inhibition, respectively [21]. PP1-beta catalytic subunit (Ppp1cb) was up-regulated in SHF > MOD and PP1 expression increased in MOD and SHF. Moreover, p-S641-Gys1 was attenuated in MOD vs. Sham and p-S698-Gys1 was attenuated in SHF vs. Sham (Excel file S3). These data suggest an increase in glycogen synthesis in MOD and SHF, along with enhanced glycogen degradation in MOD, and a paradoxical net increase in glycogen content in SHF > MOD. The untargeted metabolomics data showed a similar pattern of decrease in glycogen content in MOD vs. Sham and an increase in glycogen content in SHF vs. MOD (Supplementary Figure S12B).
The phosphofructokinase (Pfkm) enzyme, which is critical in glycolysis catalyzing the initial step in the conversion of fructose-6-phosphate to fructose-1,6-biphosphate, was down-regulated in MOD and SHF, SHF > MOD (Excel file S1). Moreover, the enzymes phosphoglycerate mutase (Pgam2) and phosphopyruvate hydratase (Eno3) were down-regulated in MOD and SHF, SHF > MOD, vs. Sham (Excel file S1). Additionally, the enzymes fructose biphosphate aldolase (Aldoa) and pyruvate kinase (Pkm) were downregulated in SHF vs. Sham. Pkm, which converts phosphoenoylpyruvate to pyruvate, catalyzes the second ATP-generating step in the glycolytic pathway. These data suggest that glycolysis was impaired in MOD and SHF, SHF > MOD. Pyruvate is then translocated into the mitochondria by the mt pyruvate carriers 1 (Mpc1) and 2 (Mpc2). Mpc1, pyruvate carboxylase (Pc), and the dihydrolipoamide acetyltransferase component of the pyruvate dehydrogenase (PDH) complex were down-regulated in MOD vs. Sham. In addition to the aforementioned proteins, Mpc2 and the PDH E1 component subunits alpha (Pdha2) and beta (Pdhb) were downregulated in SHF vs. Sham (Excel file S1). These data suggest that transitioning from MOD to SHF was associated with a significant impairment in pyruvate metabolism.
The enzyme BCAA aminotransferase (Bcat), which catalyzes the initial step converting BCAA to branched-chain alpha-ketoacids (BCKAs), was down-regulated in MOD and SHF, SHF > MOD, vs. Sham (Excel file S1). BCKAs are then irreversibly decarboxylated by the BCKA dehydrogenase (Bckdh) complex, which was mainly down-regulated in SHF. The Bckdh complex is also phosphorylated and thus inhibited by the kinase Bckdk, which was about equally down-regulated in MOD and SHF. BCKAs, through a series of reactions, are eventually catabolized into acetyl-CoA for oxidation in the TCA cycle. The enzymes responsible for these reactions were down-regulated in SHF > MOD (Excel file S1). These findings were consistent with the untargeted metabolomics data suggesting that BCAAs catabolism was decreased in SHF > MOD vs. Sham.

2.6. AA, Ketone Body and Nucleotide Metabolism

The untargeted metabolomics data showed that the abundance of taurine and the N-acyl-amines, N-stearoyl-taurine and N-Docosanoyl-taurine, decreased in MOD vs. Sham. Tauropine and N-linoleoyl-taurine abundance equally decreased in MOD and SHF vs. Sham (Excel file S2). Creatine abundance decreased in MOD and further decreased in SHF vs. MOD and Sham, whereas 3-Sulfino-alanine increased in MOD and further increased in SHF vs. MOD and Sham (Excel file S2). Glutamine abundance decreased in SHF only (Excel file S2). The AA, pantetheine-4-phosphate, abundance decreased in MOD and further down-trended in SHF vs. MOD (Supplementary Figure S12D). In addition to AA, many metabolites belonging to the glycerolipid and glycerophospholipid classes were identified and their abundance decreased in SHF vs. Sham (Excel file S2), suggesting a significant alteration in cell membrane fluidity and permeability.
Metabolic pathways related to ketogenesis and ketolysis were downregulated in SHF > MOD. The enzymes acetyl-CoA acyltransferase 2 (Acaa2) and 3-hydroxybutyrate dehydrogenase 1 (Bdh1) were downregulated in MOD and SHF, SHF > MOD (Excel file S1). Hdha was about equally down-regulated in MOD and SHF, whereas Hdhb was downregulated in MOD and further decreased in SHF vs. MOD and Sham. Acetoacetic acid abundance increased in SHF vs. Sham and up-trended in MOD vs. Sham (Supplementary Figure S12F), confirming a decrease in ketone body metabolism, as the proteomics data indicated.
The enzyme, purine nucleoside phosphorylase (Pnp), was equally up-regulated in MOD and SHF. The enzyme, 5’-nucleotidase, cytosolic II (Nt5c2), was upregulated in MOD and further increased in SHF (Excel file S1). Transitioning to SHF was associated with the up-regulation of the enzymes adenosine deaminase (Ada), guanine deaminase (Gda), and xanthine dehydrogenase (Xdh) (Excel file S1). The metabolite 2’,3’ cyclic AMP decreased in MOD and SHF vs. Sham and further down-trended in SHF vs. MOD (Figure 4E and Supplementary Figure S12E).

2.7. Mitochondrial Proteome and Oxidative Capacity

Succinate dehydrogenase complexes (Sdha and Sdhb) and the succinate CoA ligases (Sucla2 and Suclg1) were mildly downregulated in MOD vs. Sham and further downregulated, including Sdhc, in SHF. Additionally, all enzymes involved in the TCA cycle were down-regulated in SHF vs. Sham (Excel file S1). Similarly, ETC complexes II, IV, and ATP synthase expression decreased in MOD vs. Sham. All the ETC complexes and ATP synthase were down-regulated in transitioning to SHF (Excel file S1). The post-translational modification (PTM) of ATP synthase was downregulated in MOD and SHF at specific p-sites, as shown in (Excel file S3). The NAD-dependent protein deacetylases (Sirt 3 and Sirt5) were down-regulated in SHF only vs. Sham (Excel file S1). Moreover, the metabolite nicotinamide ribotide, also known as nicotinamide mononucleotide (NMN), which is an intermediate in NAD+ biosynthesis, down-trended in SHF vs. Sham (p = 0.05, Supplementary Figure S12G); and the enzyme nicotinamide-nucleotide adenylyltransferase (Nmnat3), which catalyzes the reaction forming NAD+ from NMN and ATP, down-regulated in SHF vs. Sham (Excel file S1). All in all, these data suggest that transitioning from MOD to SHF is associated not only with down-regulation but also with PTM, including the phosphorylation and hyperacetylation of the OXPHOS machinery.
In addition to the changes related to OXPHOS, the mt contact site and cristae organizing system (MICOS), mt dynamics proteins mitofusin 1 (Mfn1) and optic atrophy 1 (Opa1), were down-regulated in SHF relative to Sham. The mt ribosomal proteins and the Slc25 mt carrier family proteins were down-regulated in SHF > MOD vs. Sham. Particularly, the ADP/ATP translocase 1 (Slc25a4) and the mt CoA transporter (Slc25a42) were equally down-regulated in MOD and SHF vs. Sham. The mt tricarboxylate transporter (Slc25a1), the mt phosphate carrier (Slc25a3), and the mt 2-oxoglutarate/malate carrier (Slc25a11) were down-regulated in SHF vs. Sham. Moreover, there was PTM of Slc25a4 with a predominant decrease in mt protein phosphorylation in SHF vs. Sham (Excel file S3). Proteins involved in mt calcium uptake and efflux were predominantly affected in SHF and there was a PTM (phosphorylation) of the voltage-dependent anion channel 1 (VDAC1) rather than a change in its expression. The phosphosites, p-S57-VDAC1 and p-S241-VDAC1, were down-regulated in SHF vs. Sham (Excel file S3). The mt calcium uniporter (MCU) complex and the mt proton/calcium exchanger (Letm1) were down-regulated in SHF vs. Sham (Excel file S1). Phosphosites of mt proteins that changed in MOD and SHF vs. Sham are presented in Supplementary Figure S14B and Excel file S3.

2.8. PKA and AMPK Signaling and their Implication in Metabolic Remodeling

PKA and AMPK are two important intracellular signaling pathways implicated in myocardial metabolism. PKA phosphorylation at Ser339 was mildly attenuated in both MOD and SHF vs. Sham by Scaffold analysis. The adenylate cyclase type 5 (Adcy5), the PKA type II regulatory subunit alpha (Prkar2a), and the A-kinase anchor protein 12 (AKAP12) were down-regulated in SHF vs. Sham, while AKAP13 was upregulated in SHF vs. Sham (Excel file S1). The metabolite 2′,3′ cAMP was downregulated in MOD and further down-trended in SHF. The enzyme 2’,3’-cyclic-nucleotide 3’-phosphodiesterase (Cnp) was up-regulated in SHF vs. Sham, whereas p-S152-Pde4a and p-S45-Pde7b were down-regulated in SHF vs. Sham. The p-S16-PLN and p-S101-AKAP1 were decreased in MOD vs. Sham, while p-S103-AKAP1 down-trended in MOD and SHF vs. Sham (Supplementary Figure S13). These data suggest a decrease in PKA activity at the ER–mt interface subcompartment in MOD and a more pronounced decrease in PKA signaling at multiple cellular subcompartments in SHF.
The AMPK non-catalytic subunit beta 2 (p-S107-Prkab2) was up-regulated in MOD vs. Sham. The AMPK catalytic subunit alpha 2 (p-S377-Prkaa2) was downregulated, i.e., suggested to be inhibited, in SHF relative to Sham, whereas protein phosphatase 2 regulatory subunit B’alpha (p-S41-Ppp2r5a) was up-regulated in SHF vs. Sham (Excel file S3). Previous work showed that the protein phosphatase 2 regulatory subunit B delta isoform (Ppp2r2d) dephosphorylates Thr172 AMPK residue and inactivates it [22]. The p-S25-Acaca, a known downstream target of AMPK, was down-regulated in SHF vs. Sham (Excel file S3), and p-S485/491-AMPK was enhanced in MOD and SHF vs. Sham, while p-T172-AMPK was enhanced in MOD vs. Sham and then decreased in SHF vs. MOD.

2.9. Signaling Implicated in Mitochondrial Function other than PKA and AMPK

The protein kinase c delta type, Prkcd, was upregulated in MOD and further increased in SHF vs. Sham. Additionally, protein kinase c alpha type, Prkca, was upregulated in SHF only vs. Sham (Excel file S1. The p-S661-Prkcd was up-regulated and p-T521-Prkce was down-regulated in SHF vs. Sham (Excel file S3). MAPK signaling was activated in MOD and SHF, SHF > MOD. The mTOR signaling was activated in SHF vs. Sham. Serine/threonine protein kinases and phosphatases that changed in MOD and SHF vs. Sham are presented in Excel file S1 and Supplementary Figure S14A.
Calcium cycling was mainly impaired in SHF. Calcium cycling at the ER–mt interface was attenuated in SHF > MOD. SERCA2a expression was mildly decreased in MOD (log2 fold change: −0.2, p = 0.037) and further decreased in SHF vs. Sham (Excel file S1). p-S16-PLN was attenuated in both MOD and SHF vs. Sham. MCU complex proteins (Mcu and Micu1) were downregulated in SHF only (Excel file S1). The Letm1 was mildly decreased in SHF vs. Sham. Calreticulin (Calr) was upregulated in SHF, while calsequestrin-2 (Casq2), Prkar1a, Prkar2a, and Camk2a were downregulated in SHF vs. Sham (Excel file S1). Additionally, the plasma membrane calcium entry and efflux proteins, the voltage-dependent calcium channel (Cacna2d2) that regulates L-type calcium channel activity, and the sodium/calcium exchanger (Slc8a1) were upregulated in SHF vs. Sham, respectively (Excel file S1). Moreover, p-S339-Prkaca, p-S528-Camk2d, p-S51-Vdac1, and p-S241-Vdac1 were downregulated and p-S282-Slc8a1 was upregulated in SHF vs. Sham (Excel file S3). Collectively, these data suggest predominant impairment in calcium signaling and homeostasis as well as mt calcium uptake and efflux in SHF vs. Sham.

3. Discussion

This multi-omics-based study provides a global assessment focused on metabolic remodeling and its contribution to myocardial remodeling in a rat model of moderate remodeling and early systolic dysfunction (MOD) and in progression to SHF. It highlights the complexity of signaling pathways, their cross talk, and how metabolic intermediates regulate calcium signaling, and homeostasis, cardiac metabolism, and mt function (Figure 8). We show that the mitochondria are central in the metabolic remodeling process where metabolic pathways converge to eventually undergo OXPHOS and ATP production. We have validated that the mt proteome is predominantly down-regulated or degraded in cardiac-stress-induced myocardial dysfunction, as previously shown [23]. However, our study is also unique in that we have highlighted mt proteomic changes in a less severe MOD phenotype, as well as in an overt SHF phenotype. In both MOD and SHF, mt proteins constituted about 50% of the proteins that were down-regulated relative to Sham. Transitioning from MOD to SHF was associated with at least a 1.5-fold increase in the number and the degree of down-regulation of the commonly identified mt proteins. Moreover, other mt physiological processes, such as mt calcium uptake and efflux regulation, MICOS and the mt import system, and mt translation–elongation machinery were predominantly impaired/down-regulated in SHF. These findings support the concept that the degree of mt dysfunction develops in parallel and mirrors the degree of myocardial remodeling and progression to overt systolic dysfunction, suggesting that altered mt function is the cause rather than the consequence of HF progression, as was previously described [11]. We highlight molecular mechanisms leading to progression in mt dysfunction in HF. Our findings suggest that mt pathophysiological changes in MOD are predominantly related to attenuation in mt FA β-oxidation, BCAA catabolism and OXPHOS. Upon progression to SHF, there is a worsening of the aforementioned pathophysiological processes. Additionally, there is impaired pyruvate metabolism, enhanced mt degradation and apoptosis (impaired mt calcium regulation and reactive oxygen species (ROS) production) [5,8,9,24,25], and derangements in the mt import systems. These mt pathophysiological changes constitute the drive for the worsening in mt function and myocardial remodeling, and the transition to an advanced stage of HF.
About two-thirds of the proteins involved in metabolism and respiration that were down-regulated in MOD and SHF vs. Sham were mt proteins. Impairments in OXPHOS and mt FA β-oxidation were conspicuous in SHF > MOD, as previously shown [6]. Furthermore, our data suggest that OXPHOS and mt FA β-oxidation were the most enriched metabolic pathways that were inactivated in SHF > MOD. We showed that PTM (phosphorylation and hyperacetylation) of the OXPHOS machinery and mt proteins was more predominantly evident in SHF than MOD. Our findings are consistent with previous acetylproteomics work showing extensive mt protein lysine hyperacetylation in mouse models of early-stage HF and in end-stage human HF [26,27]. Moreover, BCAA catabolism, AA metabolism (tryptophan and glutamine degradation), and the Krebs cycle were highly enriched metabolic pathways that were attenuated in SHF > MOD. Previous work showed that the accumulation of BCAAs exerts detrimental effects on the myocardium by enhancing mTOR signaling (therefore suppressing autophagy) and ROS production in the mitochondria [28,29]. Moreover, BCAA accumulation exerts an inhibitory effect on the PDH complex by inhibiting the hexosamine biosynthetic pathway, which suppresses pyruvate metabolism and sensitizes the heart to cardiac stress [30].
Glucose metabolism was differentially regulated in MOD and SHF. Glycogen and glucose-1-phosphate degradation were activated in MOD > SHF, while glucose uptake was attenuated in SHF > MOD. Glut 4 expression was equally down-regulated in MOD and SHF; however, its AMPK-mediated PTM (activation) was attenuated in SHF only. Moreover, glycolysis and pyruvate metabolism were predominantly inactivated in SHF rather than MOD. Overall, our findings indicate that glucose metabolism is mainly impaired in SHF, as previously shown [6], but with some caveats that are noteworthy to highlight. While glucose uptake was attenuated in MOD, it seems that the heart compensates at the earliest stages of systolic dysfunction by up-regulating pathways involved in endogenous glycogen degradation and glucose biosynthesis, which are then attenuated in SHF along with significant decreases in glycolysis and pyruvate metabolism. This leads to a paradoxical net increase in glycogen content. The metabolomics data suggest a pattern of decrease in glycogen content in MOD vs. Sham and an increase in glycogen content in SHF vs. MOD. Glycogen is an important energy source in the heart, which is quickly mobilized when needed, especially under acute changes in working load. Glycogen occupies about 2% of the adult and 30% of the fetal cardiomyocyte volume [31]. Unlike liver and skeletal muscle, the heart increases its glycogen content during fasting because free FA (FFA) become the predominant substrates for oxidation and inhibit glycolysis, thus routing glucose for glycogen synthesis [32]. An increase in glycogen content in advanced HF may play a detrimental role, as has been shown in high-glycogen-content ischemic rat hearts [33], and in the extreme cases of glycogen storage diseases [34]. Based on our findings, we believe that this is likely related to the dysregulation of pathways governing glycogen homeostasis and glucose metabolism, rather than just an increase in glycogen content per se.
Our data show unique findings in terms of TAG homeostasis and the accumulation of cardiolipotoxic material. Despite that mt FA β-oxidation and TAG degradation were attenuated in SHF > MOD, the abundance of ceramides and TAG-related metabolites was attenuated in both MOD and SHF. Previous work has shown that FFA uptake into the cardiac myocyte is increased in the advanced stage of HF due to enhanced sympathetic activation and an increase in plasma FFA levels [35]. Mismatch between enhanced FFA uptake and their oxidation leads to the accumulation of toxic lipid species, such as ceramides and diacylglycerol (DAG), and lipotoxicity-induced mt dysfunction and apoptosis in advanced HF [36,37]. Our findings in MOD and SHF do not support this hypothesis, but rather point toward a decrease in FA uptake and FA processing in MOD and SHF, SHF > MOD. The FA uptake receptor (Slc27a1) was equally down-regulated in MOD and SHF, while the protein, Fabp3, and the enzyme, Acsl1, involved in FA transport and activation, respectively, were down-regulated in SHF > MOD. Our findings are consistent with previous work showing a decrease in endogenous TAG content, turnover, and oxidation in rat models of PO-induced early cardiac failure [38], and volume overload-induced HF and eccentric remodeling [39].
Our data highlight important metabolites and metabolite intermediates that play a role in myocardial function and metabolism rather than in ATP production. Taurine and creatine constitute two important AAs that play a role in cardiomyocyte function, signaling, and metabolism. Taurine and its derivatives’ abundance decreased in MOD > SHF, while creatine abundance decreased in SHF > MOD. Taurine deficiency in mice, through the genetic knockout of the taurine transporter (TauT) and TauT inhibition in rats was associated with cardiomyopathy [40]. Taurine deficiency was associated with reduced glucose and FA oxidation in isolated perfused rat heart [41], and reduced ETC complexes I and III activity with increased ROS production [42], as well as altered calcium homeostasis and signaling [43]. Taurine supplementation was shown to be beneficial in HF, and was associated with improvement in LV function [44] and exercise capacity [45]. Pantetheine-4-phosphate abundance decreased in MOD and down-trended in SHF vs. MOD. Pantetheine, also known as pantothenate, is phosphorylated by the pantothenate kinase to pantetheine-4-phosphate, which is the first step in the biosynthesis of CoA and is also used to shuttle intermediates between the active sites of enzymes involved in FA peptide and polyketide synthesis, thus playing role in cell growth and metabolism [46]. Eicosanoids, which are metabolites of arachidonic acid, were down-regulated in SHF vs. Sham. These metabolites have been shown to play an important role in calcium homeostasis and cycling and cardiac contraction through the regulation/inhibition of the L-type calcium channels (LTCC) [47] and the calcium-independent phosphorylation of troponin I and myosin light chain 2 [48].
Our data suggest an enhanced nucleotide degradation, including 2′,3′ cAMP (or 3′,5′ cAMP), in SHF > MOD. Currently, we do not know whether an increase in nucleotide metabolism is an adaptive versus maladaptive process in HF. However, previous work and our data suggest that enhanced nucleotide metabolism play a detrimental rather than beneficial role and contribute to myocardial dysfunction and metabolic remodeling. Enhanced Xdh activity (Xdh expression increased in SHF) was shown to be associated with enhanced ROS and that its inhibition may be beneficial in the context of HF [49,50]. Moreover, our data show that ketone body oxidation was attenuated in SHF, contrary to what was previously shown in animal models [51] and in human end-stage HF [52].
PKA phosphorylation at Ser339 residue was mildly attenuated in both MOD and SHF vs. Sham by Scaffold analysis. Previous work has shown that PKA phosphorylation at Ser338 (reported as Ser339 in PosphoSitePlus.org) precedes PTM of the activation loop Thr197, and is required for PKA processing and maturation as well as for its phosphorylation at Thr197 [53]. In MOD, PKA signaling was more ambiguous, showing enhanced phosphorylation of some PKA targets: ryanodine receptor isoform 2 (p-T1850-Ryr2, not previously reported) and myosin light chain 3 (p-S111-Myl3, not previously reported), as well as decreased phosphorylation of other targets: p-S641-Gys1, and p-S14-Myl2 (Excel file S3). These data suggest a differential change in PKA activity at specific cellular subcompartments in MOD (as highlighted in the subsequent section), rather than a global intracellular decline in PKA activity, and that the ER–mt interface may constitute one of the important compartments where PKA activity is down-regulated at the initial stage of cardiac remodeling and systolic dysfunction.
P-AMPK-T172 was increased in MOD, while p-S485/491-AMPK was increased in MOD and SHF vs. Sham. Previous work has shown that AMPK activation is mediated through its phosphorylation at Thr172 [54,55,56], while its phosphorylation at Ser485/491 was identified as being an auto-phosphorylation site, and a target site for protein kinase B/AKT, which is responsible for the insulin-mediated inhibition of AMPK upon insulin stimulation [57,58]. PKA has also been shown to phosphorylate AMPK at Ser485/491 and to limit its activation in response to energy depletion or other regulators [56]. The collective data suggest that AMPK activity may be enhanced or unaltered in MOD and is attenuated in SHF.
Cross talk between PKA and AMPK signaling pathways exists. Both regulate common pathways involved in carbohydrate, lipid, and protein metabolism [17,59,60,61]. Our data show derangements in FA, glucose metabolism, and OXPHOS, and a decline in PGC-1α and PPARα signaling that was evident in MOD, despite the fact that AMPK signaling was not affected in MOD. This is likely explained by derangements in PKA signaling in specific cellular subcompartments in MOD due to heterogeneity in the localization or activity of the adenylate cyclases, phosphodiesterases, AKAPs, and phosphatases, such as PP1a-induced dephosphorylation of common PKA targets in specific cellular subcompartments as previously shown [16,62,63,64]. We show that PP1a expression was increased in MOD and SHF vs. Sham and up-trended further in SHF vs. MOD. The p-S101-AKAP1, a scaffold protein localizing PKA to the outer mt membrane (OMM), was attenuated in MOD vs. Sham, while p-S103-AKAP1 down-trended in both MOD and SHF vs. Sham. Previous work has shown that AMPK phosphorylates AKAP1 at Ser103 [65], while PP1 destabilizes AKAP1, promoting its degradation [66]. These data suggest a decrease in mt PKA activity at the OMM in MOD and SHF, SHF > MOD. Mitochondrially targeted PKA have been shown to be beneficial for mt function by: 1) promoting p-S-637-DRP1, i.e., inhibiting mt fission [67,68]; 2) enhancing the activity of the mt import system [69,70,71]; and 3) by directly enhancing the activity of the ETC complexes I and IV, and ATP synthase, and therefore OXPHOS and ATP production [72,73,74,75].
Besides PKA, protein kinase C isoforms: alpha (Prkca), delta (Prkcd), and epsilon (Prkce) are known to localize to mitochondria and to regulate mt function, metabolism and apoptosis [76,77]. PKC isoform activity is dependent upon their expression level, localization, and phosphorylation status [78]. Our data suggest enhanced Prkcd signaling in SHF > MOD and a decrease in Prkce signaling in SHF. Prkcd is known to contribute to stress-induced mt dysfunction in neurons [76] and is a critical mediator of post-ischemic cardiomyocyte necrosis and contractile dysfunction [79]. Moreover, Prkcd has been shown to exert a pro-hypertrophic role in the non-stressed heart [79]. Prkce plays more of a protective role in the heart and neurons [76]. Increased Prkce activity or expression in the heart has been shown to promote adaptive hypertrophy and eutrophic cardiac growth [79]. Additionally, mt-targeted Prkce has been shown to phosphorylate ETC complex IV and to increase its activity in the heart and kidney in response to stress [78,80]. Similar to Prkcd, Prkca has been shown to play a role in maladaptive hypertrophy and cardiac dysfunction [79]. MAPKs can also translocate to the OMM and modulate mt function. The MAPK, c-Jun N-terminal kinase (JNK), which is activated in cardiac stress [10], is of particular interest as it has been shown to localize to the mitochondria via its interaction with the scaffold protein Sab [81] and to phosphorylate mitofusin 2 at Ser27 (promoting its degradation and mt fission) [82], and the bcl-2 family proteins, thereby promoting mt apoptosis [10,76,83].
In conclusion, our data provide strong evidence that metabolic remodeling plays a causative role in myocardial remodeling and progression to overt SHF. We highlight changes in the mt proteome and p-proteome in early systolic dysfunction and in overt SHF. Although decline in mt OXPHOS and mt FA β-oxidation and BCAA catabolism constitute the initial mt pathophysiological processes that are impaired in early systolic dysfunction, progression to SHF is accompanied by fundamental impairment in mt function at multiple regulatory levels, including calcium and redox balance regulation, derangements in MICOS, mt import, and mt translation–elongation systems. We have highlighted the role of important AAs and metabolites in the regulation of myocardial signaling and calcium cycling and homeostasis in both MOD and SHF. Finally, we have highlighted important myocardial signaling pathways and the complexities that govern metabolic remodeling and derangements in calcium cycling in progression to SHF. Schematic drawings of important findings in MOD and SHF vs. Sham are presented in Figure 8A–C and Supplementary Figure S15A–C. Our findings are important because they provide not only a global understanding of important myocardial pathophysiological processes involved in metabolic remodeling, but also a framework for future comprehensive mechanistic studies for exploring and identifying potential therapeutic targets in HF. Strategies restoring ER–mt–PKA signaling may constitute a potential therapeutic strategy for HF that is worth investigating in the future.

4. Materials and Methods

4.1. Experimental Model of Ascending Aortic Banding for the Creation of the MOD and SHF Phenotypes

All procedures involving the handling of animals were approved by the Animal Care and Use Committee of Tulane University and adhered to the National Institutes of Health Guide for the Care and Use of Laboratory Animals (Protocol 812, approval date 25 September 2019). The MOD and SHF phenotypes were created using the AAB procedure in Sprague Dawley rats, as previously described [24,84]. Animals underwent transthoracic echocardiography as previously described [84], to assess LV size and function and assure the presence of compensated concentric hypertrophy at week 3 post-AAB, and the development of MOD and SHF phenotypes at week 8 post-AAB. The creation and characterization, by echocardiography and hemodynamics, of these two models of HF have been described elsewhere [85]. Left ventricular myocardial tissue from the same biological samples were used to perform the proteomics and metabolomics studies; n = 3 biological samples per group. Additionally, the same Sham biological samples were used for both MOD and SHF proteomics and p-proteomics runs.

4.2. Discovery-Based Proteomic and p-Proteomic Analysis Using Tandem Mass Tags and Liquid Chromatography–Mass Spectrometry

Description of the methods for protein extraction, in-solution proteolytic digestion and tandem mass tags (TMT)-labeling, global TMT sample fractionation and Orbitrap Fusion liquid chromatography–mass spectrometry (LC–MS) peptide analysis and database searching are presented in the online supplemental materials file.

4.3. Criteria for Protein Odentification

Scaffold (version 4.9, Proteome Software Inc., Portland, OR, USA) was used to validate tandem MS-based peptide and protein identifications. Peptide and protein false discovery rates (FDRs) were set to 1% and the minimum number of identified peptides was set to 2. Protein probabilities were calculated by the Protein Prophet algorithm [86]. Proteins that contained similar peptides and could not be differentiated based on MS/MS analysis alone were grouped to satisfy the principle of parsimony. Proteins sharing significant peptide evidence were grouped into clusters.

4.4. Protein Quantification in Scaffold

Scaffold Q+ (version 4.9, Proteome Software Inc., Portland, OR, USA) was used for TMT-based peptide and protein quantification. Channels were corrected according to the algorithm described in i-Tracker [87] in all samples. Normalization was then performed on intensities across samples and spectra, as described in the statistical analysis of relative labeled mass spectrometry data from complex samples using ANOVA [88]. Median values were used for averaging. Spectra data were then log2-transformed, pruned of those matched to multiple proteins, and weighted by an adaptive intensity-weighting algorithm. Of the 64,264 spectra in the experiment at the given thresholds, 36,096 (56%) were included in quantitation. Differentially expressed proteins across experimental groups were determined by applying a permutation test with Benjamini–Hochberg correction for multiple comparisons [89]. An unadjusted p value of < 0.05 was considered significant.
MS raw, log2-transformed data were exported from Scaffold and uploaded into Qlucore bioinformatics software for data visualization and presentation of the dataset in Scaffold via the generation of heat maps, PCA plots, and Venn diagrams. Additionally, the uploaded MS raw, log2-transformed datasets were statistically analyzed in Qlucore bioinformatics software by two-group comparison T-test; p < 0.05 was considered statistically significant. Heat maps, volcano plots, PCA plots, and Venn diagrams were then generated for data visualization and presentation of the statistically analyzed dataset in Qlucore.

4.5. Bioinformatic Analyses

Statistically analyzed datasets in Scaffold and Qlucore bioinformatic software programs, including all identified proteins in the proteomic TMT multiplex experiments, were uploaded into IPA bioinformatics software. A core analysis was performed for each of the two-group comparisons (MOD vs. Sham and SHF vs. Sham) using the Ingenuity Knowledge database and a cutoff p-value of < 0.05. Analyzed datasets were then compared using IPA “Comparison Analyses” function to yield the most enriched Canonical Pathways, filtered to those related to metabolism and signaling as well as upstream regulators that were shared between the two-group comparisons. The generated heat maps were presented based on a cutoff -log10 p-value of 1.3 and a z-score of ± 2. IPA analysis content information for the MOD and SHF proteomic datasets analyzed in Qlucore were the following: Analysis Creation Date: 2021-05-03, Content version: 62,089,861 (Release Date: 17 February 2021), Analysis IDs: 33,670,484 and 33,670,824, respectively. IPA analysis content information for the MOD and SHF proteomic datasets analyzed in Scaffold were the following: Analysis Creation Date: 6 February 2021, Content version: 60,467,501 (Release Date: 19 November 2020), Analysis IDs: 27,199,584 and 27,199,664, respectively.
Description of the methods for phosphorylated Peptide TMT Enrichment, Orbitrap Fusion LC-MS phosphopeptide analysis, and database searching are presented in the online supplemental materials file.

4.6. Criteria for Protein Identification

Scaffold (version 4.9, Proteome Software Inc., Portland, OR, USA) was used to validate tandem MS-based peptide and protein identifications with the ‘mudpit’ option selected for the goal of combining the results from TiO2 and FeNTA enrichments into a single report. Peptide and protein FDRs were set to 1%, and the minimum number of identified peptides was set to 2. Protein probabilities were calculated by the Protein Prophet algorithm [86]. Proteins containing similar peptides that could not be differentiated based on an MS/MS analysis were grouped to satisfy the principle of parsimony. Proteins sharing significant peptide evidence were grouped into clusters.

4.7. Protein Quantification in Scaffold

Scaffold Q+ (version 4.9, Proteome Software Inc., Portland, OR, USA) was used for TMT-based peptide and protein quantification. Channels were corrected according to the algorithm described in i-Tracker [87] by correction factors in all samples. Normalization was then performed on intensities across samples and spectra, as described in the statistical analysis of relative labeled mass spectrometry data from complex samples using ANOVA [88]. Median values were used for averaging. Spectra data were log2-transformed, pruned of those matched to multiple proteins, and weighted by an adaptive intensity-weighting algorithm. Of the 43,180 spectra in the experiment at the given thresholds, 26,377 (61%) were included in quantitation. Differentially expressed proteins across experimental groups were determined by applying a permutation test with Benjamini–Hochberg correction for multiple comparisons [89]. An unadjusted p value of < 0.05 was considered significant.

4.8. Phosphopeptide Quantification in Scaffold PTM

Scaffold PTM (version 3.2, Proteome Software Inc., Portland, OR, USA) was used for TMT-based phosphopeptide quantification with the output from Scaffold Q+. Phosphopeptide results were normalized with the global protein quantification results from the related TMT experiment using proteinXML exports from Scaffold Q+. The PTM sites derived from Scaffold were annotated in Scaffold PTM using the site localization algorithm developed by Beausoleil et al. [90]. Scaffold PTM re-analyzes MS/MS spectra identified as modified peptides and calculates Ascore values and site localization probabilities to assess the level of confidence in each PTM localization. Scaffold PTM then combines localization probabilities for all peptides containing each identified PTM site to obtain the best-estimated probability that a PTM is present at that site. For motif analysis, PTMs were scanned for over-represented patterns in the amino acids surrounding the modification sites using the method described previously [91]. The background percentage was calculated using the target human protein sequence database described earlier.
Normalized raw phosphopeptide results with the global protein quantification were exported from Scaffold PTM as log2-transformed values and uploaded into Qlucore bioinformatics software for data visualization and presentation of the statistically analyzed dataset in Scaffold PTM via heat maps, PCA plots, and Venn diagrams. Additionally, the uploaded raw log2-transformed datasets were analyzed in Qlucore bioinformatics software using the two-group comparison T-test and a cutoff of p < 0.05 for statistical significance. Heat maps, volcano plots, PCA plots, and Venn diagrams were generated for data visualization and presentation.

4.9. Bioinformatic Analyses

Bioinformatic analysis was performed in IPA as described above, except that a “Phosphorylation Core Analysis” was performed by expression log ratio and a cutoff p-value < 0.05. Analyzed datasets were compared using IPA’s “Comparison Analyses” function to yield the most enriched signaling pathways that were present in the two-group comparisons. The generated heat maps were presented based on p-value and z-score. IPA analysis content information for the MOD and SHF p-proteomic datasets analyzed in Qlucore were the following: Analysis Creation Date: 2021-04-23, Content version: 62,089,861 (Release Date: 17 February 2021), Analysis IDs: 32,911,164 and 32,911,624, respectively. IPA analysis content information for the MOD and SHF p-proteomic datasets analyzed in Scaffold were the following: Analysis Creation Date: 17 April 2021 and 16 April 2021, respectively, Content version: 62,089,861 (Release Date: 17 February 2021), Analysis IDs: 32,434,384 and 32,433,404, respectively.

4.10. Validation by Immunoblotting

Description of the method for this section is presented in the supplemental materials file.

4.11. Discovery-Based Non-Targeted Central Carbon Metabolism and Acyl-Carnitine Analysis by LC–MS

Description of the methods for sample preparation, and central carbon and acyl-carnitine LC–MS analysis are presented in the supplemental materials file.

4.11.1. Central Carbon Metabolism Data Identification and Normalization

Metabolites were identified by matching the retention time and mass (+/− 10 ppm) to authentic standards. Peak areas were integrated using Profinder v8.00 (Agilent Technologies, Santa Clara, CA, USA). Data were normalized to urine creatinine levels and Loess drift correction was applied using the area of each metabolite in the control samples (pools) for correction, as previously described using MetaboDrift 1.0 [92]. Non-targeted data analysis was performed using Agilent’s MassHunter Find by Molecular Feature workflow (v7.0) with recursion using Agilent’s Mass Profiler Pro (v8.0.) The dataset was processed using Binner [93] to remove degenerate features, and the resulting features were searched against the Metabolomics Workbench Refmet database (https://www.metabolomicsworkbench.org/databases/refmet/index.php), last accessed date 9 November 2021, to provide Metabolomics Standards Initiative (MSI) [94] Level III identifications, or to an in-house library of authentic standards to provide MSI Level I identifications. Duplicate identified metabolites that had the same mass, retention time, and similar abundance among the studied biological samples were removed. For example, both 2′,3′ cyclic AMP (cAMP) and 3′,5′ cAMP were identified that had the same mass, retention time, and similar abundance among the studied biological samples. These were thought to be duplicate metabolites. 2′,3′ cAMP was kept and 3′,5′ cAMP was removed.

4.11.2. Acyl Carnitine Data Identification

Metabolites were identified by matching the retention time and mass (+/− 10 ppm) to authentic standards. Isotope peak areas were integrated using MassHunter Quantitative Analysis vB.07.00 (Agilent Technologies, Santa Clara, CA, USA) of metabolite fluxes for central carbon, acyl-carnitine and FA metabolites.

4.11.3. Statistical Analysis of the Metabolomics MS Raw Data

MS raw data were uploaded into Qlucore bioinformatics software, and then were log2-transformed. Statistical analysis was performed in Qlucore bioinformatics software using the multi-group comparison one way ANOVA and correction by Tukey. A p-value < 0.05 was considered statistically significant. Then, heat maps, volcano plots, PCA plots, and Venn diagrams were generated for data visualization and presentation of the statistically analyzed dataset in Qlucore.

4.11.4. Bioinformatic Analyses

Bioinformatic analysis was performed using IPA as described above, except that a “Metabolomics Core Analysis” was performed using expression log ratio and a cutoff p-value < 0.05. Analyzed datasets were compared using IPA’s “Comparison Analyses” function to yield the most enriched metabolic pathways and upstream regulators that were present in the two-group comparisons. The generated heat maps were presented based on p-value for metabolic pathways and z-score for upstream regulators. IPA analysis content information for the metabolomics dataset analyzed in Qlucore were the following: Analysis Creation Date: 23 April 2021, Content version: 62,089,861 (Release Date: 17 February 2021), Analysis ID: 32,906,044.

4.12. Statistical Analysis

Statistical methods used for the analysis of the proteomics, p-proteomics and metabolomics datasets were described in detail in previous sections. A p-value < 0.05 was considered significant. Western blot data are presented as mean ± standard deviation. Statistical analyses were performed in Prism software version 9.1.0 using one-way ANOVA with Benjamini correction method. A p-value < 0.05 was considered significant.

Supplementary Materials

Supplementary Materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms23010235/s1.

Author Contributions

Conceptualization, A.H.C.; Methodology, A.H.C., L.H., T.M. and M.K.; Formal Analysis, A.H.C. and L.H.; Investigation, A.H.C.; Resources, A.H.C., L.H., T.M., M.K. and C.B.; Data Curation, A.H.C.; Writing—Original Draft Preparation, A.H.C., L.H., T.M., J.H. and M.K., Writing—Review and Editing, A.H.C., C.B., L.G.N., D.B. and P.D.; Validation, A.H.C.; Visualization, A.H.C. and J.H.; Supervision, P.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by an internal grant award to A.H.C., and the Carol Lavin Bernick Faculty Grant award to A.H.C.; an NIH/NHLBI R01-HL070241 grant award to P.D.; NIH HL-148836, and NIH AG-063345 grant awards to D.W.B., and the Louisiana Board of Regents Endowed Chairs for Eminent Scholars program to D.W.B.

Institutional Review Board Statement

All procedures involving the handling of animals were approved by the Animal Care and Use Committee of Tulane University and adhered to the National Institutes of Health Guide for the Care and Use of Laboratory Animals (Protocol 812, approval date 25 September 2019).

Informed Consent Statement

Not applicable.

Data Availability Statement

The corresponding author have all data available upon request.

Acknowledgments

The authors thank Nancy Busija, M.A., C.C.C.-SLP for a critical reading and thorough editing of the manuscript.

Conflicts of Interest

The authors have no conflict of interest to disclose.

References

  1. Doenst, T.; Nguyen, T.D.; Abel, E.D. Cardiac metabolism in heart failure: Implications beyond ATP production. Circ. Res. 2013, 113, 709–724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ingwall, J.S.; Weiss, R.G. Is the failing heart energy starved? On using chemical energy to support cardiac function. Circ. Res. 2004, 95, 135–145. [Google Scholar] [CrossRef] [PubMed]
  3. D’Erchia, A.M.; Atlante, A.; Gadaleta, G.; Pavesi, G.; Chiara, M.; De Virgilio, C.; Manzari, C.; Mastropasqua, F.; Prazzoli, G.M.; Picardi, E.; et al. Tissue-specific mtDNA abundance from exome data and its correlation with mitochondrial transcription, mass and respiratory activity. Mitochondrion 2015, 20, 13–21. [Google Scholar] [CrossRef] [PubMed]
  4. Piquereau, J.; Caffin, F.; Novotova, M.; Lemaire, C.; Veksler, V.; Garnier, A.; Ventura-Clapier, R.; Joubert, F. Mitochondrial dynamics in the adult cardiomyocytes: Which roles for a highly specialized cell? Front. Physiol. 2013, 4, 102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Tian, R.; Colucci, W.S.; Arany, Z.; Bachschmid, M.M.; Ballinger, S.W.; Boudina, S.; Bruce, J.E.; Busija, D.W.; Dikalov, S.; Dorn, G.W., II; et al. Unlocking the Secrets of Mitochondria in the Cardiovascular System: Path to a Cure in Heart Failure—A Report from the 2018 National Heart, Lung, and Blood Institute Workshop. Circulation 2019, 140, 1205–1216. [Google Scholar] [CrossRef]
  6. Neubauer, S. The failing heart—An engine out of fuel. N. Engl. J. Med. 2007, 356, 1140–1151. [Google Scholar] [CrossRef] [Green Version]
  7. Rosca, M.G.; Tandler, B.; Hoppel, C.L. Mitochondria in cardiac hypertrophy and heart failure. J. Mol. Cell Cardiol. 2013, 55, 31–41. [Google Scholar] [CrossRef] [Green Version]
  8. Chaanine, A.H.; Sreekumaran Nair, K.; Bergen, R.H., 3rd; Klaus, K.; Guenzel, A.J.; Hajjar, R.J.; Redfield, M.M. Mitochondrial Integrity and Function in the Progression of Early Pressure Overload-Induced Left Ventricular Remodeling. J. Am. Heart Assoc. 2017, 6, e005869. [Google Scholar] [CrossRef]
  9. Chaanine, A.H.; Joyce, L.D.; Stulak, J.M.; Maltais, S.; Joyce, D.L.; Dearani, J.A.; Klaus, K.; Nair, K.S.; Hajjar, R.J.; Redfield, M.M. Mitochondrial Morphology, Dynamics, and Function in Human Pressure Overload or Ischemic Heart Disease With Preserved or Reduced Ejection Fraction. Circ. Heart Fail. 2019, 12, e005131. [Google Scholar] [CrossRef]
  10. Chaanine, A.H.; Jeong, D.; Liang, L.; Chemaly, E.R.; Fish, K.; Gordon, R.E.; Hajjar, R.J. JNK modulates FOXO3a for the expression of the mitochondrial death and mitophagy marker BNIP3 in pathological hypertrophy and in heart failure. Cell Death Dis. 2012, 3, 265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Doenst, T.; Pytel, G.; Schrepper, A.; Amorim, P.; Färber, G.; Shingu, Y.; Mohr, F.W.; Schwarzer, M. Decreased rates of substrate oxidation ex vivo predict the onset of heart failure and contractile dysfunction in rats with pressure overload. Cardiovasc. Res. 2010, 86, 461–470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Hill, J.A.; Olson, E.N. Cardiac plasticity. N. Engl. J. Med. 2008, 358, 1370–1380. [Google Scholar] [CrossRef] [PubMed]
  13. Kehat, I.; Molkentin, J.D. Molecular pathways underlying cardiac remodeling during pathophysiological stimulation. Circulation 2010, 122, 2727–2735. [Google Scholar] [CrossRef] [Green Version]
  14. Kreusser, M.M.; Lehmann, L.H.; Keranov, S.; Hoting, M.O.; Oehl, U.; Kohlhaas, M.; Reil, J.C.; Neumann, K.; Schneider, M.D.; Hill, J.A.; et al. Cardiac CaM Kinase II genes δ and γ contribute to adverse remodeling but redundantly inhibit calcineurin-induced myocardial hypertrophy. Circulation 2014, 130, 1262–1273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Wang, Y.; Tandan, S.; Hill, J.A. Calcineurin-dependent ion channel regulation in heart. Trends Cardiovasc. Med. 2014, 24, 14–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Burdyga, A.; Surdo, N.C.; Monterisi, S.; Di Benedetto, G.; Grisan, F.; Penna, E.; Pellegrini, L.; Zaccolo, M.; Bortolozzi, M.; Swietach, P.; et al. Phosphatases control PKA-dependent functional microdomains at the outer mitochondrial membrane. Proc. Natl. Acad. Sci. USA 2018, 115, E6497–E6506. [Google Scholar] [CrossRef] [Green Version]
  17. Gerhart-Hines, Z.; Dominy, J.E., Jr.; Blättler, S.M.; Jedrychowski, M.P.; Banks, A.S.; Lim, J.H.; Chim, H.; Gygi, S.P.; Puigserver, P. The cAMP/PKA pathway rapidly activates SIRT1 to promote fatty acid oxidation independently of changes in NAD(+). Mol. Cell 2011, 44, 851–863. [Google Scholar] [CrossRef] [Green Version]
  18. Luczak, E.D.; Wu, Y.; Granger, J.M.; Joiner, M.A.; Wilson, N.R.; Gupta, A.; Umapathi, P.; Murphy, K.R.; Reyes Gaido, O.E.; Sabet, A.; et al. Mitochondrial CaMKII causes adverse metabolic reprogramming and dilated cardiomyopathy. Nat. Commun. 2020, 11, 4416. [Google Scholar] [CrossRef]
  19. Pfluger, P.T.; Kabra, D.G.; Aichler, M.; Schriever, S.C.; Pfuhlmann, K.; García, V.C.; Lehti, M.; Weber, J.; Kutschke, M.; Rozman, J.; et al. Calcineurin Links Mitochondrial Elongation with Energy Metabolism. Cell Metab. 2015, 22, 838–850. [Google Scholar] [CrossRef] [Green Version]
  20. Morgan, H.E.; Parmeggiani, A. Regulation of glycogenolysis in muscle. II. control of glycogen phosphorylase reaction in isolated perfused heart. J. Biol. Chem. 1964, 239, 2435–2439. [Google Scholar] [CrossRef]
  21. Kumar, G.S.; Choy, M.S.; Koveal, D.M.; Lorinsky, M.K.; Lyons, S.P.; Kettenbach, A.N.; Page, R.; Peti, W. Identification of the substrate recruitment mechanism of the muscle glycogen protein phosphatase 1 holoenzyme. Sci. Adv. 2018, 4, eaau6044. [Google Scholar] [CrossRef] [Green Version]
  22. Joseph, B.K.; Liu, H.Y.; Francisco, J.; Pandya, D.; Donigan, M.; Gallo-Ebert, C.; Giordano, C.; Bata, A.; Nickels, J.T., Jr. Inhibition of AMP Kinase by the Protein Phosphatase 2A Heterotrimer, PP2APpp2r2d. J. Biol. Chem. 2015, 290, 10588–10598. [Google Scholar] [CrossRef] [Green Version]
  23. Dai, D.F.; Hsieh, E.J.; Chen, T.; Menendez, L.G.; Basisty, N.B.; Tsai, L.; Beyer, R.P.; Crispin, D.A.; Shulman, N.J.; Szeto, H.H.; et al. Global proteomics and pathway analysis of pressure-overload-induced heart failure and its attenuation by mitochondrial-targeted peptides. Circ. Heart Fail. 2013, 6, 1067–1076. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Chaanine, A.H.; Gordon, R.E.; Kohlbrenner, E.; Benard, L.; Jeong, D.; Hajjar, R.J. Potential role of BNIP3 in cardiac remodeling, myocardial stiffness, and endoplasmic reticulum: Mitochondrial calcium homeostasis in diastolic and systolic heart failure. Circ. Heart Fail. 2013, 6, 572–583. [Google Scholar] [CrossRef] [Green Version]
  25. Narula, J.; Pandey, P.; Arbustini, E.; Haider, N.; Narula, N.; Kolodgie, F.D.; Dal Bello, B.; Semigran, M.J.; Bielsa-Masdeu, A.; Dec, G.W.; et al. Apoptosis in heart failure: Release of cytochrome c from mitochondria and activation of caspase-3 in human cardiomyopathy. Proc. Natl. Acad. Sci. USA 1999, 96, 8144–8149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Castillo, E.C.; Morales, J.A.; Chapoy-Villanueva, H.; Silva-Platas, C.; Treviño-Saldaña, N.; Guerrero-Beltrán, C.E.; Bernal-Ramírez, J.; Torres-Quintanilla, A.; García, N.; Youker, K.; et al. Mitochondrial Hyperacetylation in the Failing Hearts of Obese Patients Mediated Partly by a Reduction in SIRT3: The Involvement of the Mitochondrial Permeability Transition Pore. Cell. Physiol. Biochem. Int. J. Exp. Cell. Physiol. Biochem. Pharmacol. 2019, 53, 465–479. [Google Scholar]
  27. Horton, J.L.; Martin, O.J.; Lai, L.; Riley, N.M.; Richards, A.L.; Vega, R.B.; Leone, T.C.; Pagliarini, D.J.; Muoio, D.M.; Bedi, K.C., Jr.; et al. Mitochondrial protein hyperacetylation in the failing heart. JCI Insight 2016, 2, e84897. [Google Scholar] [CrossRef] [PubMed]
  28. Sun, H.; Olson, K.C.; Gao, C.; Prosdocimo, D.A.; Zhou, M.; Wang, Z.; Jeyaraj, D.; Youn, J.Y.; Ren, S.; Liu, Y.; et al. Catabolic Defect of Branched-Chain Amino Acids Promotes Heart Failure. Circulation 2016, 133, 2038–2049. [Google Scholar] [CrossRef] [PubMed]
  29. Sun, H.; Wang, Y. Branched chain amino acid metabolic reprogramming in heart failure. Biochim. Biophys. Acta 2016, 1862, 2270–2275. [Google Scholar] [CrossRef]
  30. Li, T.; Zhang, Z.; Kolwicz, S.C., Jr.; Abell, L.; Roe, N.D.; Kim, M.; Zhou, B.; Cao, Y.; Ritterhoff, J.; Gu, H.; et al. Defective Branched-Chain Amino Acid Catabolism Disrupts Glucose Metabolism and Sensitizes the Heart to Ischemia-Reperfusion Injury. Cell Metab. 2017, 25, 374–385. [Google Scholar] [CrossRef] [Green Version]
  31. Depre, C.; Vanoverschelde, J.L.; Taegtmeyer, H. Glucose for the heart. Circulation 1999, 99, 578–588. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Schneider, C.A.; Nguyêñ, V.T.; Taegtmeyer, H. Feeding and fasting determine postischemic glucose utilization in isolated working rat hearts. Am. J. Physiol. 1991, 260 Pt 2, H542–H548. [Google Scholar] [CrossRef]
  33. Cross, H.R.; Opie, L.H.; Radda, G.K.; Clarke, K. Is a high glycogen content beneficial or detrimental to the ischemic rat heart? A controversy resolved. Circ. Res. 1996, 78, 482–491. [Google Scholar] [CrossRef] [PubMed]
  34. Austin, S.L.; Proia, A.D.; Spencer-Manzon, M.J.; Butany, J.; Wechsler, S.B.; Kishnani, P.S. Cardiac Pathology in Glycogen Storage Disease Type III. JIMD Rep. 2012, 6, 65–72. [Google Scholar] [PubMed] [Green Version]
  35. Opie, L.H.; Knuuti, J. The adrenergic-fatty acid load in heart failure. J. Am. Coll. Cardiol. 2009, 54, 1637–1646. [Google Scholar] [CrossRef] [Green Version]
  36. Goldberg, I.J.; Trent, C.M.; Schulze, P.C. Lipid metabolism and toxicity in the heart. Cell Metab. 2012, 15, 805–812. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Wende, A.R.; Abel, E.D. Lipotoxicity in the heart. Biochim. Biophys. Acta 2010, 1801, 311–319. [Google Scholar] [CrossRef] [PubMed]
  38. O’Donnell, J.M.; Fields, A.D.; Sorokina, N.; Lewandowski, E.D. The absence of endogenous lipid oxidation in early stage heart failure exposes limits in lipid storage and turnover. J. Mol. Cell Cardiol. 2008, 44, 315–322. [Google Scholar] [CrossRef] [Green Version]
  39. Melenovsky, V.; Benes, J.; Skaroupkova, P.; Sedmera, D.; Strnad, H.; Kolar, M.; Vlcek, C.; Petrak, J.; Benes, J., Jr.; Papousek, F.; et al. Metabolic characterization of volume overload heart failure due to aorto-caval fistula in rats. Mol. Cell. Biochem. 2011, 354, 83–96. [Google Scholar] [CrossRef]
  40. Ito, T.; Kimura, Y.; Uozumi, Y.; Takai, M.; Muraoka, S.; Matsuda, T.; Ueki, K.; Yoshiyama, M.; Ikawa, M.; Okabe, M.; et al. Taurine depletion caused by knocking out the taurine transporter gene leads to cardiomyopathy with cardiac atrophy. J. Mol. Cell Cardiol. 2008, 44, 927–937. [Google Scholar] [CrossRef] [PubMed]
  41. Schaffer, S.W.; Shimada-Takaura, K.; Jong, C.J.; Ito, T.; Takahashi, K. Impaired energy metabolism of the taurine-deficient heart. Amino Acids 2016, 48, 549–558. [Google Scholar] [CrossRef]
  42. Jong, C.J.; Azuma, J.; Schaffer, S. Mechanism underlying the antioxidant activity of taurine: Prevention of mitochondrial oxidant production. Amino Acids 2012, 42, 2223–2232. [Google Scholar] [CrossRef] [PubMed]
  43. Ramila, K.C.; Jong, C.J.; Pastukh, V.; Ito, T.; Azuma, J.; Schaffer, S.W. Role of protein phosphorylation in excitation-contraction coupling in taurine deficient hearts. American journal of physiology. Heart Circ. Physiol. 2015, 308, H232–H239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Azuma, J.; Sawamura, A.; Awata, N. Usefulness of taurine in chronic congestive heart failure and its prospective application. Jpn. Circ. J. 1992, 56, 95–99. [Google Scholar] [CrossRef] [Green Version]
  45. Beyranvand, M.R.; Khalafi, M.K.; Roshan, V.D.; Choobineh, S.; Parsa, S.A.; Piranfar, M.A. Effect of taurine supplementation on exercise capacity of patients with heart failure. J. Cardiol. 2011, 57, 333–337. [Google Scholar] [CrossRef] [Green Version]
  46. Leonardi, R.; Jackowski, S. Biosynthesis of Pantothenic Acid and Coenzyme A. EcoSal Plus 2007, 2, ecosalplus.3.6.3.4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Petit-Jacques, J.; Hartzell, H.C. Effect of arachidonic acid on the L-type calcium current in frog cardiac myocytes. J. Physiol. 1996, 493 Pt 1, 67–81. [Google Scholar] [CrossRef]
  48. Damron, D.S.; Darvish, A.; Murphy, L.; Sweet, W.; Moravec, C.S.; Bond, M. Arachidonic acid-dependent phosphorylation of troponin I and myosin light chain 2 in cardiac myocytes. Circ. Res. 1995, 76, 1011–1019. [Google Scholar] [CrossRef]
  49. Hajjar, R.J.; Leopold, J.A. Xanthine oxidase inhibition and heart failure: Novel therapeutic strategy for ventricular dysfunction? Circ. Res. 2006, 98, 169–171. [Google Scholar] [CrossRef] [Green Version]
  50. Minhas, K.M.; Saraiva, R.M.; Schuleri, K.H.; Lehrke, S.; Zheng, M.; Saliaris, A.P.; Berry, C.E.; Barouch, L.A.; Vandegaer, K.M.; Li, D.; et al. Xanthine oxidoreductase inhibition causes reverse remodeling in rats with dilated cardiomyopathy. Circ. Res. 2006, 98, 271–279. [Google Scholar] [CrossRef] [Green Version]
  51. Aubert, G.; Martin, O.J.; Horton, J.L.; Lai, L.; Vega, R.B.; Leone, T.C.; Koves, T.; Gardell, S.J.; Krüger, M.; Hoppel, C.L.; et al. The Failing Heart Relies on Ketone Bodies as a Fuel. Circulation 2016, 133, 698–705. [Google Scholar] [CrossRef] [PubMed]
  52. Bedi, K.C., Jr.; Snyder, N.W.; Brandimarto, J.; Aziz, M.; Mesaros, C.; Worth, A.J.; Wang, L.L.; Javaheri, A.; Blair, I.A.; Margulies, K.B.; et al. Evidence for Intramyocardial Disruption of Lipid Metabolism and Increased Myocardial Ketone Utilization in Advanced Human Heart Failure. Circulation 2016, 133, 706–716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Keshwani, M.M.; Klammt, C.; von Daake, S.; Ma, Y.; Kornev, A.P.; Choe, S.; Insel, P.A.; Taylor, S.S. Cotranslational cis-phosphorylation of the COOH-terminal tail is a key priming step in the maturation of cAMP-dependent protein kinase. Proc. Natl. Acad. Sci. USA 2012, 109, E1221–E1229. [Google Scholar] [CrossRef] [Green Version]
  54. Hawley, S.A.; Boudeau, J.; Reid, J.L.; Mustard, K.J.; Udd, L.; Mäkelä, T.P.; Alessi, D.R.; Hardie, D.G. Complexes between the LKB1 tumor suppressor, STRAD alpha/beta and MO25 alpha/beta are upstream kinases in the AMP-activated protein kinase cascade. J. Biol. 2003, 2, 28. [Google Scholar] [CrossRef] [Green Version]
  55. Hawley, S.A.; Pan, D.A.; Mustard, K.J.; Ross, L.; Bain, J.; Edelman, A.M.; Frenguelli, B.G.; Hardie, D.G. Calmodulin-dependent protein kinase kinase-beta is an alternative upstream kinase for AMP-activated protein kinase. Cell Metab. 2005, 2, 9–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Hurley, R.L.; Barré, L.K.; Wood, S.D.; Anderson, K.A.; Kemp, B.E.; Means, A.R.; Witters, L.A. Regulation of AMP-activated protein kinase by multisite phosphorylation in response to agents that elevate cellular cAMP. J. Biol. Chem. 2006, 281, 36662–36672. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Horman, S.; Vertommen, D.; Heath, R.; Neumann, D.; Mouton, V.; Woods, A.; Schlattner, U.; Wallimann, T.; Carling, D.; Hue, L.; et al. Insulin antagonizes ischemia-induced Thr172 phosphorylation of AMP-activated protein kinase alpha-subunits in heart via hierarchical phosphorylation of Ser485/491. J. Biol. Chem. 2006, 281, 5335–5340. [Google Scholar] [CrossRef] [Green Version]
  58. Kovacic, S.; Soltys, C.L.; Barr, A.J.; Shiojima, I.; Walsh, K.; Dyck, J.R. Akt activity negatively regulates phosphorylation of AMP-activated protein kinase in the heart. J. Biol. Chem. 2003, 278, 39422–39427. [Google Scholar] [CrossRef] [Green Version]
  59. Cantó, C.; Auwerx, J. PGC-1alpha, SIRT1 and AMPK, an energy sensing network that controls energy expenditure. Curr. Opin. Lipidol. 2009, 20, 98–105. [Google Scholar] [CrossRef] [Green Version]
  60. Fernandez-Marcos, P.J.; Auwerx, J. Regulation of PGC-1α, a nodal regulator of mitochondrial biogenesis. Am. J. Clin. Nutr. 2011, 93, 884s–890s. [Google Scholar] [CrossRef] [Green Version]
  61. O’Neill, H.M.; Holloway, G.P.; Steinberg, G.R. AMPK regulation of fatty acid metabolism and mitochondrial biogenesis: Implications for obesity. Mol. Cell. Endocrinol. 2013, 366, 135–151. [Google Scholar] [CrossRef]
  62. Di Benedetto, G.; Zoccarato, A.; Lissandron, V.; Terrin, A.; Li, X.; Houslay, M.D.; Baillie, G.S.; Zaccolo, M. Protein kinase A type I and type II define distinct intracellular signaling compartments. Circ. Res. 2008, 103, 836–844. [Google Scholar] [CrossRef] [Green Version]
  63. Haj Slimane, Z.; Bedioune, I.; Lechêne, P.; Varin, A.; Lefebvre, F.; Mateo, P.; Domergue-Dupont, V.; Dewenter, M.; Richter, W.; Conti, M.; et al. Control of cytoplasmic and nuclear protein kinase A by phosphodiesterases and phosphatases in cardiac myocytes. Cardiovasc. Res. 2014, 102, 97–106. [Google Scholar] [CrossRef] [Green Version]
  64. Hinke, S.A.; Navedo, M.F.; Ulman, A.; Whiting, J.L.; Nygren, P.J.; Tian, G.; Jimenez-Caliani, A.J.; Langeberg, L.K.; Cirulli, V.; Tengholm, A.; et al. Anchored phosphatases modulate glucose homeostasis. EMBO J. 2012, 31, 3991–4004. [Google Scholar] [CrossRef] [PubMed]
  65. Hoffman, N.J.; Parker, B.L.; Chaudhuri, R.; Fisher-Wellman, K.H.; Kleinert, M.; Humphrey, S.J.; Yang, P.; Holliday, M.; Trefely, S.; Fazakerley, D.J.; et al. Global Phosphoproteomic Analysis of Human Skeletal Muscle Reveals a Network of Exercise-Regulated Kinases and AMPK Substrates. Cell Metab. 2015, 22, 922–935. [Google Scholar] [CrossRef] [Green Version]
  66. Le, A.V.; Tavalin, S.J.; Dodge-Kafka, K.L. Identification of AKAP79 as a protein phosphatase 1 catalytic binding protein. Biochemistry 2011, 50, 5279–5291. [Google Scholar] [CrossRef] [Green Version]
  67. Chang, C.R.; Blackstone, C. Cyclic AMP-dependent protein kinase phosphorylation of Drp1 regulates its GTPase activity and mitochondrial morphology. J. Biol. Chem. 2007, 282, 21583–21587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Cribbs, J.T.; Strack, S. Reversible phosphorylation of Drp1 by cyclic AMP-dependent protein kinase and calcineurin regulates mitochondrial fission and cell death. EMBO Rep. 2007, 8, 939–944. [Google Scholar] [CrossRef] [Green Version]
  69. Ginsberg, M.D.; Feliciello, A.; Jones, J.K.; Avvedimento, E.V.; Gottesman, M.E. PKA-dependent binding of mRNA to the mitochondrial AKAP121 protein. J. Mol. Biol. 2003, 327, 885–897. [Google Scholar] [CrossRef]
  70. Opalińska, M.; Meisinger, C. Mitochondrial protein import under kinase surveillance. Microb. Cell 2014, 1, 51–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Schmidt, O.; Harbauer, A.B.; Rao, S.; Eyrich, B.; Zahedi, R.P.; Stojanovski, D.; Schönfisch, B.; Guiard, B.; Sickmann, A.; Pfanner, N.; et al. Regulation of mitochondrial protein import by cytosolic kinases. Cell 2011, 144, 227–239. [Google Scholar] [CrossRef]
  72. Carlucci, A.; Lignitto, L.; Feliciello, A. Control of mitochondria dynamics and oxidative metabolism by cAMP, AKAPs and the proteasome. Trends Cell Biol. 2008, 18, 604–613. [Google Scholar] [CrossRef] [PubMed]
  73. De Rasmo, D.; Palmisano, G.; Scacco, S.; Technikova-Dobrova, Z.; Panelli, D.; Cocco, T.; Sardanelli, A.M.; Gnoni, A.; Micelli, L.; Trani, A.; et al. Phosphorylation pattern of the NDUFS4 subunit of complex I of the mammalian respiratory chain. Mitochondrion 2010, 10, 464–471. [Google Scholar] [CrossRef] [PubMed]
  74. De Rasmo, D.; Panelli, D.; Sardanelli, A.M.; Papa, S. cAMP-dependent protein kinase regulates the mitochondrial import of the nuclear encoded NDUFS4 subunit of complex I. Cell. Signal. 2008, 20, 989–997. [Google Scholar] [CrossRef]
  75. Piccoli, C.; Scacco, S.; Bellomo, F.; Signorile, A.; Iuso, A.; Boffoli, D.; Scrima, R.; Capitanio, N.; Papa, S. cAMP controls oxygen metabolism in mammalian cells. FEBS Lett. 2006, 580, 4539–4543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Lucero, M.; Suarez, A.E.; Chambers, J.W. Phosphoregulation on mitochondria: Integration of cell and organelle responses. CNS Neurosci. Ther. 2019, 25, 837–858. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Rathore, R.; Zheng, Y.M.; Li, X.Q.; Wang, Q.S.; Liu, Q.H.; Ginnan, R.; Singer, H.A.; Ho, Y.S.; Wang, Y.X. Mitochondrial ROS-PKCepsilon signaling axis is uniquely involved in hypoxic increase in [Ca2+]i in pulmonary artery smooth muscle cells. Biochem. Biophys. Res. Commun. 2006, 351, 784–790. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Malhotra, A.; Kang, B.P.; Opawumi, D.; Belizaire, W.; Meggs, L.G. Molecular biology of protein kinase C signaling in cardiac myocytes. Mol. Cell. Biochem. 2001, 225, 97–107. [Google Scholar] [CrossRef] [PubMed]
  79. Dorn, G.W., 2nd; Force, T. Protein kinase cascades in the regulation of cardiac hypertrophy. J. Clin. Investig. 2005, 115, 527–537. [Google Scholar] [CrossRef] [Green Version]
  80. Nowak, G.; Bakajsova, D.; Clifton, G.L. Protein kinase C-epsilon modulates mitochondrial function and active Na+ transport after oxidant injury in renal cells. American journal of physiology. Ren. Physiol. 2004, 286, F307–F316. [Google Scholar] [CrossRef]
  81. Win, S.; Than, T.A.; Fernandez-Checa, J.C.; Kaplowitz, N. JNK interaction with Sab mediates ER stress induced inhibition of mitochondrial respiration and cell death. Cell Death Dis. 2014, 5, e989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Leboucher, G.P.; Tsai, Y.C.; Yang, M.; Shaw, K.C.; Zhou, M.; Veenstra, T.D.; Glickman, M.H.; Weissman, A.M. Stress-induced phosphorylation and proteasomal degradation of mitofusin 2 facilitates mitochondrial fragmentation and apoptosis. Mol. Cell 2012, 47, 547–557. [Google Scholar] [CrossRef] [Green Version]
  83. Yamamoto, K.; Ichijo, H.; Korsmeyer, S.J. BCL-2 is phosphorylated and inactivated by an ASK1/Jun N-terminal protein kinase pathway normally activated at G(2)/M. Mol. Cell. Biol. 1999, 19, 8469–8478. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Chaanine, A.H.; Kohlbrenner, E.; Gamb, S.I.; Guenzel, A.J.; Klaus, K.; Fayyaz, A.U.; Nair, K.S.; Hajjar, R.J.; Redfield, M.M. FOXO3a regulates BNIP3 and modulates mitochondrial calcium, dynamics, and function in cardiac stress. Am. J. Physiol. Heart Circ. Physiol. 2016, 311, H1540–H1559. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Chaanine, A.H.; Hajjar, R.J. Characterization of the Differential Progression of Left Ventricular Remodeling in a Rat Model of Pressure Overload Induced Heart Failure. Does Clip Size Matter? Methods Mol. Biol. 2018, 1816, 195–206. [Google Scholar] [PubMed]
  86. Nesvizhskii, A.I.; Keller, A.; Kolker, E.; Aebersold, R. A statistical model for identifying proteins by tandem mass spectrometry. Anal. Chem. 2003, 75, 4646–4658. [Google Scholar] [CrossRef]
  87. Shadforth, I.P.; Dunkley, T.P.; Lilley, K.S.; Bessant, C. i-Tracker: For quantitative proteomics using iTRAQ. BMC Genom. 2005, 6, 145. [Google Scholar] [CrossRef] [Green Version]
  88. Oberg, A.L.; Mahoney, D.W.; Eckel-Passow, J.E.; Malone, C.J.; Wolfinger, R.D.; Hill, E.G.; Cooper, L.T.; Onuma, O.K.; Spiro, C.; Therneau, T.M.; et al. Statistical analysis of relative labeled mass spectrometry data from complex samples using ANOVA. J. Proteome Res. 2008, 7, 225–233. [Google Scholar] [CrossRef] [PubMed]
  89. Glickman, M.E.; Rao, S.R.; Schultz, M.R. False discovery rate control is a recommended alternative to Bonferroni-type adjustments in health studies. J. Clin. Epidemiol. 2014, 67, 850–857. [Google Scholar] [CrossRef]
  90. Beausoleil, S.A.; Villén, J.; Gerber, S.A.; Rush, J.; Gygi, S.P. A probability-based approach for high-throughput protein phosphorylation analysis and site localization. Nat. Biotechnol. 2006, 24, 1285–1292. [Google Scholar] [CrossRef]
  91. Schwartz, D.; Gygi, S.P. An iterative statistical approach to the identification of protein phosphorylation motifs from large-scale data sets. Nat. Biotechnol. 2005, 23, 1391–1398. [Google Scholar] [CrossRef] [PubMed]
  92. Thonusin, C.; IglayReger, H.B.; Soni, T.; Rothberg, A.E.; Burant, C.F.; Evans, C.R. Evaluation of intensity drift correction strategies using MetaboDrift, a normalization tool for multi-batch metabolomics data. J. Chromatogr. A 2017, 1523, 265–274. [Google Scholar] [CrossRef] [PubMed]
  93. Kachman, M.; Habra, H.; Duren, W.; Wigginton, J.; Sajjakulnukit, P.; Michailidis, G.; Burant, C.; Karnovsky, A. Deep annotation of untargeted LC-MS metabolomics data with Binner. Bioinformatics 2020, 36, 1801–1806. [Google Scholar] [CrossRef] [PubMed]
  94. Sumner, L.W.; Amberg, A.; Barrett, D.; Beale, M.H.; Beger, R.; Daykin, C.A.; Fan, T.W.; Fiehn, O.; Goodacre, R.; Griffin, J.L.; et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metab. Off. J. Metab. Soc. 2007, 3, 211–221. [Google Scholar]
Ijms 23 00235 g001a 550Ijms 23 00235 g001b 550
Figure 1. Visualization of the MOD and SHF proteomic datasets. (A) Heat maps showing differential log2 fold change in expression of identified proteins that changed in MOD vs. Sham (left), and in SHF vs. Sham (right). The same Sham biological samples were used as a reference for both TMT-label proteomic runs; n = 3 biological samples per group. (B) Volcano plots of the identified proteins that changed in MOD (left) and SHF (right) groups vs. Sham. The x-axis displays the log2 fold change in group means for down-regulated proteins (green) and up-regulated proteins (red) in MOD or SHF vs. Sham. The y-axis displays the -Log10 p-value of identified proteins that were statistically significant in MOD (left) and SHF (right) vs. Sham. A -Log10 p-value cutoff of 1.3 (p = 0.05) and higher was considered significant. (C) Venn diagrams showing the number of identified proteins that changed in SHF (red) and MOD (green) groups, divided into those whose relative expression decreased (left) or increased (right) vs. Sham. The blue and yellow diagrams show the total number of identified mitochondrial (mt) proteins in the MOD and SHF proteomic runs, respectively. (D) Heat maps showing differential log2 fold change in expression of proteins related to metabolism and respiration that changed in MOD vs. Sham (left), and in SHF vs. Sham (right). (E) Volcano plots showing the log2 fold change in group means (x-axis) and the corresponding -Log10 p-value on the y-axis of the identified proteins presented in (D). (F) Venn diagrams showing the total number of identified proteins related to metabolism and respiration that decreased (left) or increased (right) in MOD and SHF groups, green and red diagrams, respectively. They also show how many of those were mt proteins (intersection of green and blue diagrams for MOD and red and yellow diagrams for SHF) vs. non-mitochondrial proteins (numbers that fall outside the intersection of green and blue diagrams for MOD and red and yellow diagrams for SHF).
Figure 1. Visualization of the MOD and SHF proteomic datasets. (A) Heat maps showing differential log2 fold change in expression of identified proteins that changed in MOD vs. Sham (left), and in SHF vs. Sham (right). The same Sham biological samples were used as a reference for both TMT-label proteomic runs; n = 3 biological samples per group. (B) Volcano plots of the identified proteins that changed in MOD (left) and SHF (right) groups vs. Sham. The x-axis displays the log2 fold change in group means for down-regulated proteins (green) and up-regulated proteins (red) in MOD or SHF vs. Sham. The y-axis displays the -Log10 p-value of identified proteins that were statistically significant in MOD (left) and SHF (right) vs. Sham. A -Log10 p-value cutoff of 1.3 (p = 0.05) and higher was considered significant. (C) Venn diagrams showing the number of identified proteins that changed in SHF (red) and MOD (green) groups, divided into those whose relative expression decreased (left) or increased (right) vs. Sham. The blue and yellow diagrams show the total number of identified mitochondrial (mt) proteins in the MOD and SHF proteomic runs, respectively. (D) Heat maps showing differential log2 fold change in expression of proteins related to metabolism and respiration that changed in MOD vs. Sham (left), and in SHF vs. Sham (right). (E) Volcano plots showing the log2 fold change in group means (x-axis) and the corresponding -Log10 p-value on the y-axis of the identified proteins presented in (D). (F) Venn diagrams showing the total number of identified proteins related to metabolism and respiration that decreased (left) or increased (right) in MOD and SHF groups, green and red diagrams, respectively. They also show how many of those were mt proteins (intersection of green and blue diagrams for MOD and red and yellow diagrams for SHF) vs. non-mitochondrial proteins (numbers that fall outside the intersection of green and blue diagrams for MOD and red and yellow diagrams for SHF).
Ijms 23 00235 g001aIjms 23 00235 g001b
Ijms 23 00235 g002 550
Figure 2. Top-ranked Canonical Pathways related to metabolism and Upstream Regulators, resulting from IPA’s ‘Comparison Analyses’ of the experimental groups. (A) ‘Core Analyses’ of the statistically analyzed proteomic datasets, by both Qlucore (Q) and Scaffold (S) bioinformatics software, were performed in IPA for each of the two-group comparisons, i.e., MOD vs. Sham and SHF vs. Sham, with a cutoff p-value < 0.05. Analyzed datasets were then subsequently compared with each other in IPA’s ‘Comparison Analyses’ function to yield the most enriched Canonical Pathways (left) related to metabolism and Upstream Regulators (right) that were shared among the two-group comparisons. The respective z-score-based heat maps indicate shared Canonical Pathways and Upstream Regulators that were up-regulated/activated or down-regulated/inhibited in the two-group comparisons, with orange and blue color intensities representing the z-score-based extent of up-regulation/activation or down-regulation/inhibition, respectively. Canonical Pathways and Upstream Regulators with a z-score outside the cutoff value are represented by a dot on the heat map. (B) Western blot in left ventricular myocardium showing the expression of the peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC-1a), and proteins involved in mt fatty acid β-oxidation (HADHA, MCAD and HADHB), * p < 0.05 vs. Sham and † p < 0.05 vs. MOD.
Figure 2. Top-ranked Canonical Pathways related to metabolism and Upstream Regulators, resulting from IPA’s ‘Comparison Analyses’ of the experimental groups. (A) ‘Core Analyses’ of the statistically analyzed proteomic datasets, by both Qlucore (Q) and Scaffold (S) bioinformatics software, were performed in IPA for each of the two-group comparisons, i.e., MOD vs. Sham and SHF vs. Sham, with a cutoff p-value < 0.05. Analyzed datasets were then subsequently compared with each other in IPA’s ‘Comparison Analyses’ function to yield the most enriched Canonical Pathways (left) related to metabolism and Upstream Regulators (right) that were shared among the two-group comparisons. The respective z-score-based heat maps indicate shared Canonical Pathways and Upstream Regulators that were up-regulated/activated or down-regulated/inhibited in the two-group comparisons, with orange and blue color intensities representing the z-score-based extent of up-regulation/activation or down-regulation/inhibition, respectively. Canonical Pathways and Upstream Regulators with a z-score outside the cutoff value are represented by a dot on the heat map. (B) Western blot in left ventricular myocardium showing the expression of the peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC-1a), and proteins involved in mt fatty acid β-oxidation (HADHA, MCAD and HADHB), * p < 0.05 vs. Sham and † p < 0.05 vs. MOD.
Ijms 23 00235 g002
Ijms 23 00235 g003 550
Figure 3. Oxidative phosphorylation in MOD and SHF. (A) Oxidative phosphorylation (OXPHOS) network of the electron transport chain (ETC) complexes I-IV transport system and ATP synthase are presented in MOD vs. Sham (upper panel) and SHF vs. Sham (lower panel). The network highlights mt OXPHOS proteins in each ETC complex that were downregulated (green) vs. Sham, were unchanged (grey) vs. Sham, or were not identified (white) in the MOD and SHF proteomic runs. The darker the green color, the higher the degree of downregulation. (B) Expression of the ETC complexes I-IV and ATP synthase (complex V) in MOD and SHF by immunoblotting, * p < 0.05 vs. Sham and † p < 0.05 vs. MOD.
Figure 3. Oxidative phosphorylation in MOD and SHF. (A) Oxidative phosphorylation (OXPHOS) network of the electron transport chain (ETC) complexes I-IV transport system and ATP synthase are presented in MOD vs. Sham (upper panel) and SHF vs. Sham (lower panel). The network highlights mt OXPHOS proteins in each ETC complex that were downregulated (green) vs. Sham, were unchanged (grey) vs. Sham, or were not identified (white) in the MOD and SHF proteomic runs. The darker the green color, the higher the degree of downregulation. (B) Expression of the ETC complexes I-IV and ATP synthase (complex V) in MOD and SHF by immunoblotting, * p < 0.05 vs. Sham and † p < 0.05 vs. MOD.
Ijms 23 00235 g003
Ijms 23 00235 g004a 550Ijms 23 00235 g004b 550
Figure 4. Visualization of the non-targeted metabolomics dataset in Sham, MOD, and SHF groups. (A) Heat maps showing differential change in log2 fold expression of metabolites that changed in MOD vs. Sham (left), and in SHF vs. Sham (right); n = 3 samples per group. (B) Volcano plots showing the log2 fold change in group means, on the x-axis, and the corresponding -Log10 (p-value), on the y-axis, for the metabolites that were downregulated (green) or upregulated (red) in MOD (left) and SHF (right) vs. Sham. A -Log10 (p-value) cutoff of 1.3 and higher was considered significant. (C) PCA plots showing the variance in biological samples within and between the Sham, MOD, and SHF groups for the metabolites that changed in MOD (left) and SHF (right) vs. Sham. (D) Venn diagram showing the number of metabolites that changed in SHF (red) and MOD (blue) vs. Sham and those that changed in SHF vs. MOD (green). (E) Volcano plot of the 11 common identified metabolites that changed in MOD and SHF vs. Sham (intersection of red and blue diagrams excluding the green diagram). Of those, there was a trend for further downregulation of pantetheine 4′-phosphate and 2′,3′ cyclic AMP in SHF compared with MOD. (F) Heat map and PCA plot for the 17 common metabolites that changed in MOD and SHF vs. Sham. (G) Heat map and PCA plot for the metabolites that changed in SHF vs. Sham and MOD groups. Heat maps and PCA plots in both F and G show the differential log2 fold change in expression in Sham, MOD, and SHF groups and the variance in biological samples, respectively, for the metabolites presented in (F) and (G).
Figure 4. Visualization of the non-targeted metabolomics dataset in Sham, MOD, and SHF groups. (A) Heat maps showing differential change in log2 fold expression of metabolites that changed in MOD vs. Sham (left), and in SHF vs. Sham (right); n = 3 samples per group. (B) Volcano plots showing the log2 fold change in group means, on the x-axis, and the corresponding -Log10 (p-value), on the y-axis, for the metabolites that were downregulated (green) or upregulated (red) in MOD (left) and SHF (right) vs. Sham. A -Log10 (p-value) cutoff of 1.3 and higher was considered significant. (C) PCA plots showing the variance in biological samples within and between the Sham, MOD, and SHF groups for the metabolites that changed in MOD (left) and SHF (right) vs. Sham. (D) Venn diagram showing the number of metabolites that changed in SHF (red) and MOD (blue) vs. Sham and those that changed in SHF vs. MOD (green). (E) Volcano plot of the 11 common identified metabolites that changed in MOD and SHF vs. Sham (intersection of red and blue diagrams excluding the green diagram). Of those, there was a trend for further downregulation of pantetheine 4′-phosphate and 2′,3′ cyclic AMP in SHF compared with MOD. (F) Heat map and PCA plot for the 17 common metabolites that changed in MOD and SHF vs. Sham. (G) Heat map and PCA plot for the metabolites that changed in SHF vs. Sham and MOD groups. Heat maps and PCA plots in both F and G show the differential log2 fold change in expression in Sham, MOD, and SHF groups and the variance in biological samples, respectively, for the metabolites presented in (F) and (G).
Ijms 23 00235 g004aIjms 23 00235 g004b
Ijms 23 00235 g005 550
Figure 5. Fatty-acyl-carnitine abundance, top-ranked Metabolic Pathways, and Upstream Regulators in the MOD and SHF. (A,B) Heat map and PCA plot showing the differential log2 fold change in expression and variance of biological samples for fatty-acyl-carnitines that changed in MOD and SHF groups vs. Sham. (C) ‘Core Analyses’ of the statistically analyzed dataset, by Qlucore bioinformatics software, was performed in IPA for each of the two-group comparisons and then were compared in IPA’s ‘Comparison Analyses’ function. The generated heat maps show the most enriched Metabolic Pathways (upper panel) by p-value and Upstream Regulators (lower panel) by z-score. (D) Mechanistic network showing connection between calcium, as a second messenger, and metabolites implicated in calcium cycling and homeostasis. Please refer to prediction legend for details.
Figure 5. Fatty-acyl-carnitine abundance, top-ranked Metabolic Pathways, and Upstream Regulators in the MOD and SHF. (A,B) Heat map and PCA plot showing the differential log2 fold change in expression and variance of biological samples for fatty-acyl-carnitines that changed in MOD and SHF groups vs. Sham. (C) ‘Core Analyses’ of the statistically analyzed dataset, by Qlucore bioinformatics software, was performed in IPA for each of the two-group comparisons and then were compared in IPA’s ‘Comparison Analyses’ function. The generated heat maps show the most enriched Metabolic Pathways (upper panel) by p-value and Upstream Regulators (lower panel) by z-score. (D) Mechanistic network showing connection between calcium, as a second messenger, and metabolites implicated in calcium cycling and homeostasis. Please refer to prediction legend for details.
Ijms 23 00235 g005
Ijms 23 00235 g006 550
Figure 6. Visualization of the MOD and SHF phospho (p)-proteomic datasets. (A) Heat maps showing differential log2 fold change in the expression of p-sites that changed in MOD vs. Sham (left), and in SHF vs. Sham (right). The same biological samples in the Sham group were used as a reference for both TMT-label p-proteomic runs; n = 3 biological samples per group. Note: there was a technical problem in one of the Sham biological samples in the MOD p-proteomic run, which was eliminated as it was skewing the data. (B) Volcano plots showing the degree of log2 fold change, on the x-axis, with the corresponding -Log10 p-value, on the y-axis, for the p-sites that were downregulated (green) or upregulated (red) in MOD (left) and SHF (right) groups compared with Sham. (C) PCA plots showing the variance of biological samples in the Sham, MOD and SHF groups for the p-sites that changed in MOD (left) and SHF (right) vs. Sham.
Figure 6. Visualization of the MOD and SHF phospho (p)-proteomic datasets. (A) Heat maps showing differential log2 fold change in the expression of p-sites that changed in MOD vs. Sham (left), and in SHF vs. Sham (right). The same biological samples in the Sham group were used as a reference for both TMT-label p-proteomic runs; n = 3 biological samples per group. Note: there was a technical problem in one of the Sham biological samples in the MOD p-proteomic run, which was eliminated as it was skewing the data. (B) Volcano plots showing the degree of log2 fold change, on the x-axis, with the corresponding -Log10 p-value, on the y-axis, for the p-sites that were downregulated (green) or upregulated (red) in MOD (left) and SHF (right) groups compared with Sham. (C) PCA plots showing the variance of biological samples in the Sham, MOD and SHF groups for the p-sites that changed in MOD (left) and SHF (right) vs. Sham.
Ijms 23 00235 g006
Ijms 23 00235 g007 550
Figure 7. Top-ranked Signaling Pathways in the MOD and SHF groups. (A) ‘Core Analyses’ of the statistically analyzed p-proteomic datasets by Qlucore and Scaffold bioinformatics software was performed in IPA for each of the two-group comparisons (MOD vs. Sham and SHF vs. Sham), which were then compared in IPA’s ‘Comparison Analyses’ function to yield the most enriched signaling Pathways that were shared among the two-group comparisons. The generated heat maps are presented by p-value (left) and by z-score (right). (B) Western blot showing post-translational modification of AMPK and PKA downstream protein targets as well as protein phosphatase 1 (PP1) expression in the MOD and SHF groups, * p < 0.05 vs. Sham and † p < 0.05 vs. MOD.
Figure 7. Top-ranked Signaling Pathways in the MOD and SHF groups. (A) ‘Core Analyses’ of the statistically analyzed p-proteomic datasets by Qlucore and Scaffold bioinformatics software was performed in IPA for each of the two-group comparisons (MOD vs. Sham and SHF vs. Sham), which were then compared in IPA’s ‘Comparison Analyses’ function to yield the most enriched signaling Pathways that were shared among the two-group comparisons. The generated heat maps are presented by p-value (left) and by z-score (right). (B) Western blot showing post-translational modification of AMPK and PKA downstream protein targets as well as protein phosphatase 1 (PP1) expression in the MOD and SHF groups, * p < 0.05 vs. Sham and † p < 0.05 vs. MOD.
Ijms 23 00235 g007
Ijms 23 00235 g008a 550Ijms 23 00235 g008b 550
Figure 8. Schematic drawing highlighting metabolic pathways, calcium cycling, and implicated signaling pathways in Sham (A), and their directional change in MOD vs. Sham (B), and SHF vs. Sham (C). Green arrows promote pathway or activity. Rectangles show signaling pathways or proteins that were inhibited or down-regulated (green), respectively, and those that were activated or up-regulated (red), respectively. The deeper the color, the higher the intensity of activation/up-regulation or inhibition/down-regulation, as shown by the scale bar at the bottom of the figure. Oxidative phosphorylation (OXPHOS), fatty acid (FA) metabolism, and mitochondrial (mt)-FA β-oxidation were the most enriched/inhibited metabolic pathways in MOD vs. Sham, followed by BCAA catabolism and enzymes of the tricarboxylic acid (TCA) cycle. Glycolysis and pyruvate metabolism were the least enriched/inhibited in MOD vs. Sham. The metabolites taurine and 2′,3′cAMP were downregulated in MOD vs. Sham. Endoplasmic reticulum (ER)–mt interface protein kinase A (PKA) signaling and calcium cycling as well as PGC-1α signaling were attenuated in MOD vs. Sham. Glycogen degradation and protein kinase C isoform delta (PKCδ) signaling were activated in MOD vs. Sham. In SHF, there was further inhibition or activation of the aforementioned metabolic and signaling pathways, except for glycogen degradation, which decreased in SHF vs. MOD. Unique findings in SHF included decreased PKA activity throughout the cellular compartment, decreased AMPK and mt PKC isoform epsilon (PKCε), and mt calcium/calmodulin kinase 2 delta (Camk2δ) signaling. Derangements in calcium cycling and calcium homeostasis were also evident in SHF. Mitochondrial pathophysiological processes, other than OXPHOS, that are indirectly involved in or regulating OXPHOS (mt calcium uptake and efflux, MICOS complex, mt import and carrier proteins, and mt protein deacetylation) were attenuated in SHF vs. Sham. Additionally, mt translation/elongation was attenuated in SHF. Eicosanoids decreased in abundance in SHF vs. Sham, while 2′,3′cAMP down-trended in SHF vs. MOD. Abbreviations: TGFβ: transforming growth factor beta, Glut4: glucose transporter family 4, Slc27a1: long-chain FA transport protein 1, Fabp3: FA binding protein 3, CoA: Coenzyme A, BCAA: branched-chain amino acids, BCKAs: branched-chain ketoacids, Vdac1: voltage-dependent anion channel isoform 1, Mcu: mt calcium uniporter, Cpt1b: carnitine-O-palmitoyltransferase isoform b, Sirt3: Sirtuin 3, Mrpl: mt 39S ribosomal proteins, Slc25: mt solute carrier family proteins, ANT1: ADP/ATP translocase 1, Letm1: proton/calcium exchanger, NCLX: sodium/calcium/lithium exchanger, PP1: protein phosphatase 1, DRP1: dynamin-related protein 1, AKAP: A-kinase anchor protein, PLN: phospholamban, SERCA2a: Sarco/endoplasmic reticulum calcium ATPase, Calr: calreticulin, Casq2: Calsequestrin-2, Ryr2: ryanodine receptor isoform 2, NCX-1: sodium/calcium exchanger, LTCC: L-type calcium channel, VGCC: voltage-gated calcium channels, GPCR: G protein-coupled receptor, Plcd3: 1-phosphatidylinositol 4,5-bisphosphate phosphodiesterase delta-3, Adcy5: adenylate cyclase type 5, Gai: G protein subunit alpha i2, Gas: GNAS complex locus, and PDE: phosphodiesterase. Note: Both 2′, 3′ cAMP and 3′,5′ cAMP were identified and were considered duplicate metabolites (please refer to methods section for details). They are represented as 2′,3′ cAMP here.
Figure 8. Schematic drawing highlighting metabolic pathways, calcium cycling, and implicated signaling pathways in Sham (A), and their directional change in MOD vs. Sham (B), and SHF vs. Sham (C). Green arrows promote pathway or activity. Rectangles show signaling pathways or proteins that were inhibited or down-regulated (green), respectively, and those that were activated or up-regulated (red), respectively. The deeper the color, the higher the intensity of activation/up-regulation or inhibition/down-regulation, as shown by the scale bar at the bottom of the figure. Oxidative phosphorylation (OXPHOS), fatty acid (FA) metabolism, and mitochondrial (mt)-FA β-oxidation were the most enriched/inhibited metabolic pathways in MOD vs. Sham, followed by BCAA catabolism and enzymes of the tricarboxylic acid (TCA) cycle. Glycolysis and pyruvate metabolism were the least enriched/inhibited in MOD vs. Sham. The metabolites taurine and 2′,3′cAMP were downregulated in MOD vs. Sham. Endoplasmic reticulum (ER)–mt interface protein kinase A (PKA) signaling and calcium cycling as well as PGC-1α signaling were attenuated in MOD vs. Sham. Glycogen degradation and protein kinase C isoform delta (PKCδ) signaling were activated in MOD vs. Sham. In SHF, there was further inhibition or activation of the aforementioned metabolic and signaling pathways, except for glycogen degradation, which decreased in SHF vs. MOD. Unique findings in SHF included decreased PKA activity throughout the cellular compartment, decreased AMPK and mt PKC isoform epsilon (PKCε), and mt calcium/calmodulin kinase 2 delta (Camk2δ) signaling. Derangements in calcium cycling and calcium homeostasis were also evident in SHF. Mitochondrial pathophysiological processes, other than OXPHOS, that are indirectly involved in or regulating OXPHOS (mt calcium uptake and efflux, MICOS complex, mt import and carrier proteins, and mt protein deacetylation) were attenuated in SHF vs. Sham. Additionally, mt translation/elongation was attenuated in SHF. Eicosanoids decreased in abundance in SHF vs. Sham, while 2′,3′cAMP down-trended in SHF vs. MOD. Abbreviations: TGFβ: transforming growth factor beta, Glut4: glucose transporter family 4, Slc27a1: long-chain FA transport protein 1, Fabp3: FA binding protein 3, CoA: Coenzyme A, BCAA: branched-chain amino acids, BCKAs: branched-chain ketoacids, Vdac1: voltage-dependent anion channel isoform 1, Mcu: mt calcium uniporter, Cpt1b: carnitine-O-palmitoyltransferase isoform b, Sirt3: Sirtuin 3, Mrpl: mt 39S ribosomal proteins, Slc25: mt solute carrier family proteins, ANT1: ADP/ATP translocase 1, Letm1: proton/calcium exchanger, NCLX: sodium/calcium/lithium exchanger, PP1: protein phosphatase 1, DRP1: dynamin-related protein 1, AKAP: A-kinase anchor protein, PLN: phospholamban, SERCA2a: Sarco/endoplasmic reticulum calcium ATPase, Calr: calreticulin, Casq2: Calsequestrin-2, Ryr2: ryanodine receptor isoform 2, NCX-1: sodium/calcium exchanger, LTCC: L-type calcium channel, VGCC: voltage-gated calcium channels, GPCR: G protein-coupled receptor, Plcd3: 1-phosphatidylinositol 4,5-bisphosphate phosphodiesterase delta-3, Adcy5: adenylate cyclase type 5, Gai: G protein subunit alpha i2, Gas: GNAS complex locus, and PDE: phosphodiesterase. Note: Both 2′, 3′ cAMP and 3′,5′ cAMP were identified and were considered duplicate metabolites (please refer to methods section for details). They are represented as 2′,3′ cAMP here.
Ijms 23 00235 g008aIjms 23 00235 g008b
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chaanine, A.H.; Higgins, L.; Markowski, T.; Harman, J.; Kachman, M.; Burant, C.; Navar, L.G.; Busija, D.; Delafontaine, P. Multi-Omics Approach Profiling Metabolic Remodeling in Early Systolic Dysfunction and in Overt Systolic Heart Failure. Int. J. Mol. Sci. 2022, 23, 235. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23010235

AMA Style

Chaanine AH, Higgins L, Markowski T, Harman J, Kachman M, Burant C, Navar LG, Busija D, Delafontaine P. Multi-Omics Approach Profiling Metabolic Remodeling in Early Systolic Dysfunction and in Overt Systolic Heart Failure. International Journal of Molecular Sciences. 2022; 23(1):235. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23010235

Chicago/Turabian Style

Chaanine, Antoine H., LeeAnn Higgins, Todd Markowski, Jarrod Harman, Maureen Kachman, Charles Burant, L. Gabriel Navar, David Busija, and Patrice Delafontaine. 2022. "Multi-Omics Approach Profiling Metabolic Remodeling in Early Systolic Dysfunction and in Overt Systolic Heart Failure" International Journal of Molecular Sciences 23, no. 1: 235. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms23010235

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop