If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
∗ Wen-jing Ding and Xue-hui Li contributed equally to this work.
Wen-jing Ding
Footnotes
∗ Wen-jing Ding and Xue-hui Li contributed equally to this work.
Affiliations
Department of Geriatric Medicine, Qilu Hospital of Shandong University, Jinan 250012, ChinaKey Laboratory of Cardiovascular Proteomics of Shandong Province, Qilu Hospital of Shandong University, Jinan 250012, China
∗ Wen-jing Ding and Xue-hui Li contributed equally to this work.
Xue-hui Li
Footnotes
∗ Wen-jing Ding and Xue-hui Li contributed equally to this work.
Affiliations
Department of Geriatric Medicine, Qilu Hospital of Shandong University, Jinan 250012, ChinaKey Laboratory of Cardiovascular Proteomics of Shandong Province, Qilu Hospital of Shandong University, Jinan 250012, China
Department of Geriatric Medicine, Qilu Hospital of Shandong University, Jinan 250012, ChinaKey Laboratory of Cardiovascular Proteomics of Shandong Province, Qilu Hospital of Shandong University, Jinan 250012, China
Department of Geriatric Medicine, Qilu Hospital of Shandong University, Jinan 250012, ChinaInstitute of Basic Medical Sciences, Qilu Hospital of Shandong University, Jinan 250012, China
Department of Geriatric Medicine, Qilu Hospital of Shandong University, Jinan 250012, ChinaKey Laboratory of Cardiovascular Proteomics of Shandong Province, Qilu Hospital of Shandong University, Jinan 250012, China
Department of Geriatric Medicine, Qilu Hospital of Shandong University, Jinan 250012, ChinaKey Laboratory of Cardiovascular Proteomics of Shandong Province, Qilu Hospital of Shandong University, Jinan 250012, ChinaInstitute of Basic Medical Sciences, Qilu Hospital of Shandong University, Jinan 250012, China
Department of Geriatric Medicine, Qilu Hospital of Shandong University, Jinan 250012, ChinaKey Laboratory of Cardiovascular Proteomics of Shandong Province, Qilu Hospital of Shandong University, Jinan 250012, China
# Zhen Zhang and Yan-qiu Xing are co correspondents.
Zhen Zhang
Correspondence
Correspondence to: Zhen Zhang, Department of Geriatric Medicine, Qilu Hospital of Shandong University, 107 Wen Hua Xi Rd, Jinan, Shandong 250012, China,
# Zhen Zhang and Yan-qiu Xing are co correspondents.
Affiliations
Department of Geriatric Medicine, Qilu Hospital of Shandong University, Jinan 250012, ChinaKey Laboratory of Cardiovascular Proteomics of Shandong Province, Qilu Hospital of Shandong University, Jinan 250012, China
# Zhen Zhang and Yan-qiu Xing are co correspondents.
Yan-qiu Xing
Correspondence
Correspondence to: Yan-qiu Xing, Department of Geriatric Medicine, Qilu Hospital of Shandong University, 107 Wen Hua Xi Rd, Jinan, Shandong 250012, China,
# Zhen Zhang and Yan-qiu Xing are co correspondents.
Affiliations
Department of Geriatric Medicine, Qilu Hospital of Shandong University, Jinan 250012, ChinaKey Laboratory of Cardiovascular Proteomics of Shandong Province, Qilu Hospital of Shandong University, Jinan 250012, China
Alteration in myocardial ultrastructure and diastolic dysfunction of AMPKα2 knockout mice.
•
Proteomic and β-hydroxybutyrylation modification omics analysis of cardiac tissues in AMPKα2 knockout mice and wild-type mice.
•
Insights from β-hydroxybutyrylation modification omics on AMPK regulation of fatty acid degradation and TCA cycle.
ABSTRACT
AMP-activated protein kinase alpha 2 (AMPKα2) regulates energy metabolism, protein synthesis, and glucolipid metabolism myocardial cells. Ketone bodies (KB) produced by fatty acid β-oxidation, especially β-hydroxybutyrate (β-OHB), are fatty energy-supplying substances for the heart, brain, and other organs during fasting and long-term exercise. They also regulate metabolic signaling for multiple cellular functions. Lysine β-hydroxybutyrylation (Kbhb) is a β-OHB mediated protein post-translational modification (PTMs). Histone Kbhb has been identified in yeast, mouse, and human cells. However, whether AMPK regulates protein Kbhb is yet unclear. Hence, the present study explored the changes in proteomics and Kbhb modification omics in the hearts of AMPKα2 knockout mice using a comprehensive quantitative proteomic analysis. Based on mass spectrometry (LC-MS/MS) analysis, the number of 1181 Kbhb modified sites in 455 proteins were quantified between AMPKα2 knockout (AK) mice and wild-type (WT) mice; 244 Kbhb sites in 142 proteins decreased or increased after AMPKα2 knockout (fold change >1.5 or <1/1.5, P<0.05). The regulation of Kbhb sites in 26 key enzymes of fatty acid degradation and tricarboxylic acid cycle (TCA cycle) was noted in AMPKα2 knockout mouse cardiomyocytes. These findings, for the first time, identified proteomic features and Kbhb modification of cardiomyocytes after AMPKα2 knockout, suggesting that AMPKα2 regulates energy metabolism by modifying protein Kbhb.
]. AMPK consists of three subunits in various isoforms: the catalytic α (α1 and α2) subunits and the regulatory β (β1, β2) and γ (γ1, γ2, and γ3) subunit. In mouse hearts, α subunit 2 (AMPKα2) is the dominant catalytic subtype that accounts for 70–80% of AMPK activity [
], regulating energy metabolism, cell growth, oxidative stress, cell growth, glucose and fatty acid oxidation. AMPK is the main signaling molecule regulating myocardial energy. It acts on many metabolism-related transporters and enzymes; for example, hydroxymethylglutaryl-CoA reductase (HMGCR), glucose transporter 1 (GLUT1) and GLUT4, glycogen synthase (GS), and acetyl-CoA carboxylase (ACC) [
]. Previous studies have shown that through phosphorylation of key substrates, AMPK can promote the synthesis of acetyl-CoA and ketone bodies by suppressing the expression of carnitine palmitoyltransferase 1 (CPT-1), which limits the rate-limiting enzymes in mitochondrial fatty acid β-oxidation[
Intermediates were formed from fatty acid β-oxidation: ketone bodies consisting of β-OHB, acetoacetate (Acac), and acetone. β-OHB and others were produced in the mitochondria of hepatocytes that are subsequently secreted into the blood and utilized by organs with ketolytic mechanisms, such as the heart [
]. The classical role of β-OHB is an energy source; however, recent studies have shown that β-OHB induces Kbhb in cells, and this modification may contribute to tumor progression and energy regulation [
]. In 2021, Koronowski et al. found that ketogenic diet leads to fatty acid oxidation to derive acetyl-CoA, produces excessive β-OHB, and increases S-adenosyl-L-homocysteine hydrolase (AHCY) Kbhb, which in turn decreases AHCY activity and alters the methionine cycle [
]. To date, Kbhb is mostly described in histone proteins and barely described in non-histone proteins, especially metabolism proteins. Thus, we conducted proteomic and post-translational modification omics studies in AMPKα2 knockout mouse hearts to comprehensively demonstrate the regulatory mechanism of AMPKα2 in fatty acid oxidation, TCA cycle, and other processes at the protein level. The current findings showed that AMPKα2-specific knockout did not alter weight,grip, heart weight, tibia length in mice but altered the myocardium mitochondrial ultrastructure and cardiac diastolic function. Finally, we identified 1582 Kbhb modified sites in 585 proteins, and 244 Kbhb modified sites in 142 proteins were significantly altered (fold change >1.5 or <1/1.5, P<0.05) after AMPKα2 knockout. Functional enrichment analysis showed that dramatically downregulated Kbhb sites were related to fatty acid metabolism, amino acid metabolism, citric acid cycle, cGMP-PKG signaling pathway, and congenital heart disease pathway. Upregulated Kbhb sites were involved in arginine and proline metabolism, antigen processing and presentation, and dilated cardiomyopathy.
Furthermore, we observed that a subset of differential Kbhb sites were located on key proteins for fatty acid degradation and TCA cycle, suggesting that AMPKα2 probably affects the functions of the crucial proteins through Kbhb at these sites.
MATERIALS AND METHODS
Experimental Design and Statistical Rationale
Our study was based on AK mice and WT mice. Heart samples were analyzed from three mice per genotype (biological replicates). Based on three-dimensional assessments of reproducibility (PCC, PCA, RSD), statistical power was deemed to be sufficient (Fig. S1). Student’s t tests were used to calculate the P value, and P value < 0.05 was considered as the significance index. For physiological data, independent Student’s t test or one-way analysis of variance, were used to assess differences between WT mice and AK mice.
Generation of AMPKα2 knockout mice
CRISPR/Cas9 gene-editing technology was used to cut the protein-coding region of the C57BL/6J mice target gene AMPKα2, and the mouse fertilized egg cells were repaired by non-homologous end-joining (NHEJ), resulting in fragment deletion in the protein coding region, making the AMPKα2 protein ineffective, thereby achieving the gene knockout. The genotype of the AMPKα2 knockout was verified by PCR testing. Briefly stated, mouse tail tissue was used to extract DNA. Using the primers listed below: AMPKα2 wild type (Fig. 1A, a) (forward: 5′-TGACATCCTGTGGTGCTGAA-3′, reverse: 5′-CTGCCTAGTGCTGACTCTGA-3′) and AMPK α2 knockout (Fig. 1A, b) (forward: 5’-GCAGAGGCAGGCGAATTTC-3′, reverse:5′-GATTGTTCACTGGCTAATCTTAAGC-3′). For electrophoresis, the DNA products were placed on a 2% agarose gel and pictures were acquired. Only DNA band near 582 bp indicates knock out homozygote, and just one DNA band around 468 bp denotes heterozygote homozygote. A wild-type genotype is indicated by both the 1500-bp and 582-bp DNA bands obtained using PCR. (Fig 1A).
Fig. 1Phenotype of AK and WT mice. A, Genotyping PCR (a) AMPKα2 wild-type band. (b) AMPK α2 knockout band. Just one DNA band near 582 bp denotes knockout homozygote. B, Western blot was used to ascertain the expression of proteins such AMPKα 2 and AMPK α1 in the heart. C, Mice appearance and cardiac structure did not alter significantly after AMPKα2 knockout. D, Histogram indicated that weight, grip, heart weight (HW), Tibia length (TL), HW/BW, HW/TL were not significant in WT mice and AK mice. (n=6).
Experimental mouse relative programs were approved to the Laboratory Animal Ethical and Welfare Committee of Shandong University Cheeloo College of Medicine (Approval No. 20157). Wild-type C57BL/6J mice (WT mice) were purchased from Jinan Tengli Trade Co. Ltd, and AK mice (C57BL/6J) were generated and identified by Beijing Viewsolid Biotechnology Co. Ltd. A total of 10 WT mice and 10 AK mice were raised in the specific pathogen-free animal chamber of Shandong University Laboratory Animal Center. At the age of 3 months, the body weight (BW), grip and echocardiographic of the animals were measured and the changes in the activity, reaction, and fur were recorded. The room temperature was maintained at 20–25 °C, the humidity was 70%, and light was regulated (12 h of light/12 h of dark). The mice were allowed free access to the feed and potable water.
Grip Strength
A grip strength meter (Jiangsu Science Biological Technology Co. Ltd, Jiangsu, China) was used to measure the front paws grip strength of WT and AK mice. To test the mice's front paw grip, we took hold of the mouse at the base of the tail and lowered it vertically toward the bar. The mouse's tail pulled it slightly backwards as its two forelimbs (paws) grabbed the bar. The peak force in newton (N) was recorded at the time the mouse released its paws from the bar in each measurement. With 30 seconds of rest in between each trial, mice received 5 trials. A final grip per mouse was determined by averaging the top three out of five trials.
Echocardiographic imaging
After three months, 2% isoflurane anesthesia was used to perform transthoracic echocardiography on a heated platform utilizing the Vevo2100 imaging equipment (VisualSonics, Toronto, Canada). M-mode echocardiography was used to evaluate the left ventricular ejection fraction (LVEF), left ventricular end-diastolic diameter (LVEDD), and fractional shortening (FS) in the parasternal long-axis view. Pulsed Doppler was used to measure the early (E) and late (A) diastolic mitral flow velocities in the apical four-chamber view, and the ratio of E/A was calculated. The ratio of E'/A' was calculated by measuring the early (E') and late (A') diastolic mitral annulus velocities using tissue Doppler imaging in the apical four-chamber view.
Tissue preparation
Mice were weighed and anesthetized by intraperitoneal injection with 10% chloral hydrate. After weighing the heart, a part of the myocardial tissue was stored at −80 °C, while the remaining was used for preparing light and electron microscopy sections.
Blood parameter determinations
Blood was drawn for testing. Serum concentrations of triglycerides (TG) and total cholesterol were determined using an automatized chemistry analyzer (Rayto, Shenzhen, China). Free fatty acids test kit was used to measure the fatty acids in the serum of mice (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). Serum β-hydroxybutyrate were assessed with Beta-Hydroxybutyric acid ELISA kits (USCN Life Science, Wuhan, China). All data were expressed as mean ±SD. P-value <0.05 was considered statistically significant.
Evaluation of myocardial morphology
The mice cardiac tissue was fixed in 4% formaldehyde solution (3-5d), dehydrated in ethanol (Kelong Chemical Reagent Factory, Chengdu), paraffin embedded, and paraffin sectionalized (LEICA, RM2235). Sections were cut into 4μm thick and stained with hematoxylin-eosin (hematoxylin, Sigma, 041M0014V; eosin, Maikun Chemical Co. Ltd, 20120831). Digital microscope (×400) (OLYMPUS, DX45) was used to capture images of the heart slice.
Electron microscopy
Mice cardiac tissue was sliced into 1-mm3 organization blocks on ice and instantly immersed in 3% glutaraldehyde in cacodilate buffer at 4 °C for 2 h before post-fixation in 1% osmium-tetroxide phosphate buffer about 2 h. Subsequently, the sections were embedded in epoxide resin after dehydration in a graded ethanol series with acetone. Using a Marada-G2 type CCD photo ultrastructure of pictures, the semithin slices were cut into ultrathin sections (about 70 nm). Finally, the sections were examined under a JEM-1200EX (JEOL, Japan) transmission electron microscope after dying with uranyl acetate and lead citrate.
Western Blot
WT and AK mice heart lysates were separated on sodium dodecyl-polyacrylamide gel electrophoresis (SDS-PAGE) gels and transferred to nitrocellulose filter membrane (0.22 μm, New York, USA). Then, we blocked with 5% non-fat milk in TBST. The membrane was incubated with the primary antibody at 4°C overnight and then with the secondary antibody at room temperature for 1 h. Chemiluminescence was used to detect proteins (apparatus: Chemiluminescent Imaging System,shanghai, China). The following antibodies were used:AMPK alpha 1 polyclonal antibody (proteintech,Chicago, USA), AMPK Alpha 2 Polyclonal antibody (proteintech,Chicago, USA); anti-Kbhb antibody (PTMBioLab, Hangzhou, China), anti-succinylated lysine (PTM Biolab, Hangzhou, China),anti-malonylated lysine (PTM Biolab, Hangzhou, China),anti-crotonylated lysine (PTM Biolab, Hangzhou, China),anti-lactylated lysine (PTM Biolab, Hangzhou, China), mouse IgG (H+L) secondary antibody (Pierce, Rockford, USA), rabbit IgG (H+L) secondary antibody (Pierce, Rockford, USA).
Protein extraction and trypsin degradation
Briefly, AK and WT mice heart samples were ground into powder, followed by the addition of lysis buffer (8 M urea, 1% protease inhibitors, 3 μM TSA, and 50 mM NAM) and centrifugation at 12,000 g for 10 minutes at 4 °C. The supernatant was then gathered.
To analyze the proteomics of the heart tissue, the protein solution was alkylated with 11 mM iodoacetamide at room temperature in the dark after being reduced with 5 mM dithiothreitol at 56 °C. The protein sample was then diluted by adding 100 mM TEAB tourea concentration less than 2 M. Finally, trypsin was added for the first digestion (1:50 trypsin-to-protein mass ratio, overnight) and for a further digestion (1:100 trypsin-to-protein mass ratio, 4 h). The C18 SPE column desalted the peptides in the end. The peptides were separated using a reverse phase HPLC column (Agilent 300Extend C18, 5 m particle size, 4.6 mm id, 250 mm length). LC/MS analysis was performed after vacuum freeze-drying. For Kbhb, the supernatants of AK and WT mouse cardiac protein samples were mixed with 20% TCA (m/v) (Sigma-Aldrich, St. Louis, Missouri, USA) to deposit proteins and washed with pre-cooled acetone. Subsequently, the protein pellets were resuspended in 200 mM triethylammonium bicarbonate buffer (Sigma-Aldrich, St. Louis, Missouri, USA), dispersed by sonication, and subjected to trypsinization. The samples were reduced with 5 mM dithiothreitol (Sigma-Aldrich, St. Louis, Missouri, USA) and alkylated with iodoacetamide (Sigma-Aldrich, St. Louis, Missouri, USA). Finally, the peptides were desalted by Strata X SPE column.
Kbhb-modified enrichment
To enrich Kbhb peptides, mouse cardiac samples were mixed with NETN buffer (100 mM NaCl, 1 mM EDTA, 50 mM Tris-HCl, 0.5% NP-40, pH 8.0), which was subsequently incubated with antibody beads (PTM Bio,Catalog # PTM1204) at 4°C overnight. Then, the mixture was rinsed twice with H2O (ThermoFisher Scientific, Massachusetts, USA) and four times with NETN buffer solution. The attached peptides were washed out of the antibody beads with 0.1% trifluoroacetic acid (Sigma-Aldrich, St. Louis, Missouri, USA). Then, the eluted portions were pooled and vacuum-dried. The cardiac samples were desalted for LC-MS/MS analysis using C18 ZipTips (Millipore, Massachusetts, USA), according to the manufacturer’s instructions.
LC-MS/MS Analysis
The tryptic peptides were added in solvent A (0.1% formic acid, 2% acetonitrile/in water) and injected on a home-made reverse-phase analytical column with integratedspray tip (100 μm i.d. × 25 cm) packed with 1.9 μm/120 Å ReproSil-PurC18 resins (Dr. Maisch GmbH, Ammerbuch, Germany). The peptides were evaluated using a capillary source using Bruker Daltonics timsTOF Pro (Bruker Daltonics, USA) mass spectrometry. The TOF detectors were used to examine the precursors and fragments, with a scan range of 100–1700 m/z on the MS/MS. In PASEF mode, the timsTOF Pro was employed (parallel accumulation serial fragmentation). Following the selection of precursors with charge states 0–5 for fragmentation, 10 PASEF-MS/MS images were collected per cycle.
The final data of MS/MS were analyzed the MaxQuant search engine (v.1.6.15.0). Tandem mass spectra were retrieved against the Mus musculus SwissProt database (updated on 21/07/21, 17089 entries) concatenated with the reverse decoy database. The cleavage enzyme, trypsin/P, was selected and allowed up to four missed cleavages (proteomics allowed up two missed cleavages). In the initial and primary searches, the mass tolerance for precursor ions was set to 20 ppm, while the fragment ion mass tolerance was set to 20 ppm. Carbamidomethyl on Cys was specified as fixed modification, but oxidation of Met, Kbhb of Lys, and acetylation of the protein’s N-terminal were variable (acetylation of the protein’s N-terminal and oxidation of Met were variable in proteomics). False discovery rate <1%.
Data analysis
Proteins/peptides in the potential contaminant database and reverse decoy database were excluded. For the proteomic analysis, Label-free quantification (LFQ) of MaxQuant was used to determine intensities and normalize protein quantities. LFQ intensity values were transformed to the relative quantitative value after centralization (the LFQ intensity of each sample divided by the mean of all samples). Proteins that were detected in at least two unique peptide sampleswere included for further analysis. For Kbhb proteomics, the cutoff of site with localization probability estimated by MaxQuant was required to be 0.75 or higher. Intensity values were transformed to the relative quantitative value after centralization (each sample's intensity of modified peptides divided by the average of all samples). the relative quantitative value of the modified peptide is typically divided by the relative quantitative value of the corresponding protein to remove the influence from protein expression of modifications. Proteins/sites with a p value < 0.05 were deemed significant using the student t test. The fold-change (mean values of AK mice/WT mice) > 1.5 or 1/1.5, respectively, was used to define up- and down-regulated proteins/sites.
Bioinformatics Methods
Gene Ontology (GO) annotated proteome from UniProt-GOA database (http://www.ebi.ac.uk/GOA/). Then, the differential proteins and differential Kbhb proteins were classified according to the three categories of biological process (BPs), cellular component (CCs) and molecular function (MFs) using gene ontology annotation. Kyoto Encyclopedia of Genes and Genomes (KEGG) connects known information on molecular interaction networks. In this study, we analyzed the metabolic pathways of differential proteins and differential Kbhb proteins to elucidate the effect of AMPKα2 knockout in mice heart based on KEGG enrichment. Differential Kbhb proteins and differential proteins were searched against the STRING database (v. 11.0) for protein-protein interactions between AK mice and WT mice. Followed by visualization using Cytoscape software (v. 3.7.2). In the outer circle, the ratio value following log2 treatment was displayed using the Omics Visualizer 1.3.0 plug-in. The color is more yellow the higher the value and more purple the lower the value. The protein's location with the biggest difference served as the ratio value for the modified group (ranked by the absolute value of ratio after log2 treatment). Confidence score ≥ 0.7.
Statistical analysis
Heart samples were analyzed from three mice per genotype, and reproducibility was assessed using relative standard deviation (RSD), Principal Component Analysis (PCA), and Pearson's Correlation Coefficient (PCC)(Fig. S1). T-test was used to analyze the differences in body weight (BW), grip, heart weight (HW), tibia length (TL), HW/BW and HW/TL between AK and WT mice, and graphs were drawn in GraphPad Prism 8 software. Quantitative values are presented as mean ± standard error of mean (SEM). The statistical significance was indicated by P<0.05.
RESULTS
General condition of AK and WT mice
Both WT and AK mice were sacrificed at 3 months of age. The genotype and the phenotype of the AK mice was verified by PCR testing and western blot (WB) (Fig 1A, Fig 1B). The appearance of AK mice was indistinguishable from that of the WT mice (Fig 1C). Compared to WT mice, no conspicuous difference was detected in the BW, grip, HW, TL, HW/BW and HW/TL after AMPKa2 knockout (Fig 1D). AK mice were characterized by higher plasma FFA concentrations than were found in WT mice, while triglyceride (TG), free fatty acid (FFA) and β-OHB levels were no difference (Table 1)
Table 1Blood parameters of WT mice and AK mice
WT mice
AK mice
β- hydroxybutyrate(μg/ml)
0.03±0.01
0.04±0.02
FFA(mM)
2.15±0.29
2.72±0.17*
TG(mM)
0.65±0.13
0.70±0.12
Total cholesterol(mM)
2.01±0.07
1.95±0.04
FFA: free fatty acids; TG:triglyceride. All data are expressed as mean ± SEM (n = 6 per group). *p<0.05
The histological myocardial alterations in WT mice and AK mice under a light microscope are depicted in Fig. 2Aa and 2Ab. In the WT group, the color of myocardium was uniform, the morphology was consistent, and the boundaries were clear. A few cardiac cell nucleus in the AK mice had swollen, along with the loose cytoplasm. The ultrastructure of myocardial cells was displayed in Fig. 2Ac and 2Ad. In the WT group, the myocardium mitochondria, round or oval shape, were arranged orderly, and the cristae in the whole mitochondrial cavity were dense and regularly distributed (Fig. 2A, c). However, compared to the WT group, the mitochondria of AK mice were distributed randomly, and the mitochondrial crest flattened. (Fig. 2A, d).
Fig. 2Cardiac characteristics of AK mice and WT mice. A, AMPKα2 knockout destroyed the myocardial structure in mice. Representative images of HE staining of myocardium from WTY (A, a) and AKY (A, b) at ×400 magnification. Black scale bar represents 50 mm. Electron microscopic images of myocardium from WT mice (A, c) and AK mice (A, d). Black scale bar represents 1 mm. B, (b1) left ventricular long-axis views on a two-dimensional echocardiography, with scale bar in millimeters. (b2) M-mode echocardiograms displaying the size of the left ventricle, a scale bar in millimeters, and a time stamp in seconds on the right and bottom. (b3) Pulse-wave Doppler echocardiograms depicting mitral inflow velocities, time stamp in seconds is at the bottom, with the scale bar in mm/s on the right. (b4) Mitral annular velocities are shown on tissue Doppler echocardiograms with time stamp in seconds on the bottoms and a scale bar in mm/s on the right.
To reveal functional changes induced by AMPKα2 knockout, we performed echocardiography in AK mice and WT mice (Fig. 2B, Table 2). To determine the left ventricular systolic function, FS and LVEF were measured by M-mode measurements. All these parameters were no significantly decreased after AMPKα2 knockout. Worthwhile, after AMPKα2 knockout, an apparent decrease in the LV diastolic functional parameter, such as E/A was observed. LVEDD and E'/A' exhibited no differences after AMPKα2 knockout. Taken together, these results indicate that AMPKα2 knockout slightly impairs ventricular diastolic function.
Table 2Effects of AMPKα2 knockout on mouse left ventricular function and dimension in mice
WT mice
AK mice
LVEF (%)
81.06±6.99
69.933±17.45
FS (%)
48.85±8.26
40.08±13.92
LVEDD (mm)
2.90±0.09
3.18±0.07
E/A
1.79±0.25
0.675±0.15*
E'/A'
1.55±0.84
1.24±0.07
LVEF: left ventricular ejection fraction; FS: fractional shortening; LVEDD: Left ventricular end-diastolic diameter; E/A: the ratio of the early to late diastolic mitral inflow velocities; E'/A': the ratio of the early to late diastolic mitral annular velocities. *p<0.05,n=5 per group.
Quantification and identification of the myocardial proteome
A total of 3659 proteins were quantified in AK and WT mice. Finally, 522 differential proteins were obtained by defining the truncation value as 1.5-fold. Among, these, 57% (299/522) and 43% (223/522) were upregulated and downregulated differential proteins after AMPKα2 knockout (Fig. 3B). Supplementary data summarizes the relevant information (name, molecular weight, subcellular localization, and main functions) of the proteins.
Fig. 3Whole proteome analysis of AK mice and WT mice. A, Flowchart for the identification of omics in myocardial samples from WT and AK mice. B, Volcano plot of differential protein quantify for proteomics. Yellow represents upregulated proteins and blue represents downregulated proteins. C, Subcellular structure diagram of differential proteins. D, GO analysis results of differential proteins (yellow: biological process; green: cellular component; purple: molecular function). E, COG/KOG pathway of differential proteins (purple: cellular processes and signaling, green: information storage and processing, yellow: metabolism, blue: poorly characterized). F, Top enriched items for KEGG pathway analysis in differential proteins. The vertical axis is the functional category or pathway, and the horizontal axis is the proportion of differentially expressed modified proteins in this functional type compared to the proportion of identified proteins. Gradation from yellow to purple indicates a decreasing P value.
In order to understand the effect of AMPKα2 on the mouse heart, differential proteins were divided into various subcellular locations, biological processes, and molecular functions based on GO analysis, clusters of Orthologous Groups (COG/KOG) functional classification, and KEGG pathway analysis. The subcellular location of all the differentially expressed proteins is annotated in Fig. 3C. The largest proportion of differential proteins was localized in the cytoplasm. Most of the differential proteins were in the cytoplasm (159, 30.46%), followed by nucleus (102, 19.54%) and extracellular differential proteins (99, 18.97%). In addition, only 14.56% of the differential proteins were distributed in mitochondria (Fig. 3C). GO analysis revealed that differential proteins were aggregated into multiple biological processes, related to biological processes (BPs) (56%), cellular components (CCs) (26%), and molecular functions (MFs) (18%) (Fig. 3D).
All differential proteins in the WT and AK mice were categorized into several cellular functions, associated with energy production and conversion, glycolipid metabolism, genetic material processing and modification, and regulation of cell cycle (Fig. 3E). Notably, after AMPK knockout, about 49% of the differential proteins were related to cellular processes and signal transduction (Fig. 3E). About 38% of proteins are involved in metabolism and information storage processing. In addition to classifying the functional sorts using GO analysis, proteins with prominent differential enrichment were noted based on KEGG analysis. The complement and coagulation cascades and Staphylococcus aureus infection were significant in the KEGG pathway analysis for differential proteins. Circadian rhythm and systemic lupus erythematosus were also related to the differential proteins (Fig. 3F).
Kbhb modification omics of WT and AK
PTMs are crucial for protein function (activity and stability) and are closely related to various cardiovascular diseases [
]. Based on western blot (WB) analysis, we showed five post-translational modifications (PTMs), including succinylation, crotonylation, malonylation, β-hydroxybutyrylation, and lactylation (Fig. 4). Following AMPKα2 knockout, the levels of succinylation, crotonylation, malonylation, β-hydroxybutyrylation were all noticeably reduced. β-OHB mediated Kbhb has attracted our attention. KB, considered efficient fuel, is meaningful in heart's energy supply[
Trimetazidine Ameliorates Myocardial Metabolic Remodeling in Isoproterenol-Induced Rats Through Regulating Ketone Body Metabolism via Activating AMPK and PPAR alpha.
]. At present, β-OHB mediated Kbhb has not been reported in the heart. To further explore the effect of AMPK on the Kbhb of cardiac tissue, we conducted a proteomic analysis of Kbhb.
Fig. 4Lysine succinylation (K-succ), crotonylation (K-Cr), malonylation (K-Mal), β-hydroxybutyrylation (K-Bhb), and lactylation (K-La) levels of WT mice and AK mice.
A total of 1582 Kbhb sites were identified from 585 proteins. Among these, 313 (53.5%) had only one Kbhb site, and 272 (46.5%) had >2 Kbhb sites (Fig. 5A), while, some proteins had >10 Kbhb sites and were heavily modified (Fig. 5B). Furthermore, Kbhb substantially modified some proteins. Titin (Ttn) protein is a key component of vertebrate rhabdoid muscle assembly and function, consisting of 180 Kbhb sites. Myosin 6 (Myh 6), involved in muscle contraction, contains 48 Kbhb sites, while isocitrate dehydrogenase consisting of 16 Kbhb site is involved in the TCA cycle and acts in the intermediary metabolism and energy production (Fig. 5B). Interestingly, several Kbhb sites in these proteins might have an influence on their activity and functions.
Fig. 5Quantitative Kbhb-modified proteomics in AK and WT mice. A, Distribution of protein Kbhb sites. 313 (53.5%) had only one Kbhb site. B, Statistical map of proteins with >10 Kbhb sites. C, Volcano plot of Kbhb-modified differential sites. D, Subcellular distribution of Kbhb protein. E, Functional enrichment of Kbhb-modified proteins by the Gene Ontology database (yellow: biological process; green: cellular component; purple: molecular function). F, COG/KOG functional classification chart of differential Kbhb proteins. G, Top enriched items for KEGG pathway analysis. (Orange: up KEGG pathway, purple: down KEGG pathway).
Western blot analysis confirmed that AMPKα2 knockout significantly reduces the overall Kbhb level in the heart. A 1.5-fold change resulted in 244 differential Kbhb sites of 142 proteins. Subsequently, 140 differential Kbhb sites of 84 proteins were downregulated, while 104 differential Kbhb sites of 70 proteins were upregulated after AMPKα2 knockout (Fig. 5C). The subcellular localization of all differentially Kbhb proteins is annotated in Fig. 5D. Most of the differential Kbhb proteins were located in the mitochondria (50/35.46%) and cytoplasm (49/34.75%), suggesting that Kbhb mainly acts through the mitochondria (Fig. 5D).
The differential Kbhb modified proteins were sorted according to BPs, CCs, and MFs by GO annotation (Fig. 5E). The top two biological processes were cellular process (134, 20%) and metabolic processes (107, 16%). Further functional analysis based on COG/KOG category revealed that several Kbhb proteins were associated with metabolism, especially lipid metabolism (19, 25%) (Fig. 5F). KEGG pathway analysis showed that prominent downregulated Kbhb modified sites were associated with amino acid and fatty acid metabolism, cGMP-PKG signaling, and TCA cycle, and congenital heart disease pathways. The upregulated Kbhb sites were involved in arginine and proline metabolism, antigen processing and presentation, and dilated cardiomyopathy (Fig. 5G). Next, the flanking region of the site was analyzed, and we found that glycine (G), tryptophan (W), leucine (L), phenylalanine (F) and tyrosine (Y) were overexpressed at −1 and +1 positions around the Lys Kbhb sites (Fig. 6A).
Fig. 6Characterizing Kbhb in mouse myocardium. A, The map shows the frequency change of amino acids close to the Kbhb site. Glycine (G), tryptophan (W), leucine (L), phenylalanine (F) and tyrosine (Y) were overexpressed, −1 and + 1 positions around the Lys Kbhb sites. B-C, All Kbhb were equally divided into four groups according to their fold difference (Q1<0.5, Q2: 0.5–1/1.5, Q3: 1.5–2.0, Q4:>2.0). B, Statistical chart of quantity of each group. C, KEGG pathway of the Kbhb proteome. Fat digestion and absorption, TCA cycle, cholesterol metabolism, and glyoxylate and dicarboxylate metabolism pathway were enriched with low AK/WT ratio. D, Interaction network of differential Kbhb and proteome between TCA cycle and fatty acid degradation pathway. The Kbhb group is blue, while the protein group is red. The number of interacting nodes determines the size of the circle. In the outer circle, the ratio value following log2 treatment was displayed using the Omics Visualizer 1.3.0 plug-in. The color is more yellow the higher the value and more purple the lower the value. The protein's location with the biggest difference served as the ratio value for the modified group (ranked by the absolute value of ratio after log2 treatment).
To further study the effect of Kbhb after AMPKα2 gene knockout, Kbhb sites were divided into four groups based on differential expression (Q1<0.5, Q2: 0.5–1/1.5, Q3: 1.5–2.0, Q4>2.0), followed by KEGG pathway analysis for each group (Fig. 6B). In the KEGG category, pathways related to RNA degradation, HIF-1 signaling pathway, and glycolysis/gluconeogenesis were enriched with high AK/WT ratio. Conversely, fat digestion and absorption, TCA cycle, cholesterol metabolism, and glyoxylate and dicarboxylate metabolism pathways were significantly downregulated with low AK/WT ratio (Fig. 6C).
In this study, KEGG category showed that the part of significantly downregulated Kbhb sites in AK mice were located in proteins interrelated to mitochondrial fatty acid degradation and TCA cycle (Table 3, Table 4). Therefore, the interaction network of differential Kbhb proteins and differential proteins was analyzed based on the interaction gene/protein retrieval tool (STRING v.11.0) database, including 26 proteins related to lipid degradation and TCA cycle (Fig. 6D, Fig. S2). Notably, 10 of the key enzymes had >2 different Kbhb sites. CS had 5 differential Kbhb sites, and CPT-1 had 2 differential Kbhb sites. Multiple differential Kbhb sites in these fatty acid oxidation- and TCA cycle-related proteins might have significant effects on energy metabolism (Supplementary data 2).
Table 3Down-regulated Kbhb sites on mitochondrial fatty acid degradation
ATP production in the mouse heart is mainly dependent on mitochondrial fatty acid β-oxidation, with ketone body (KB) and others supplying only fraction of the precursor [
Trimetazidine Ameliorates Myocardial Metabolic Remodeling in Isoproterenol-Induced Rats Through Regulating Ketone Body Metabolism via Activating AMPK and PPAR alpha.
]. AMPK is an energy regulator involved in various physiological processes, such as restoring energy balance, protein synthesis, and glucose and fatty acid metabolism [
]. β-hydroxybutyrate (β-OHB) along with acetoacetate (AcAc) and acetone are the constituents of ketone bodies. β-OHB is an intermediate in mitochondrial fatty acid β-oxidation and is defined as a substrate that maintains metabolic homeostasis and a metabolic signal-regulating lipolysis, oxidative stress, and cellular function [
]. Herein, we showed that the AMPKα2 knockout disrupts the normal ultrastructure of the cardiomyocytes with ventricular diastolic function impaired as detected by echocardiography and reported a proteomic analysis of Kbhb in WT mice and AK mice. Based on HPLC-MS/MS, we identified 522 up- or downregulated differential proteins in cardiac tissues, while 244 differential Kbhb modified sites were detected after AMPKα2 knockout. This study, for the first time, observed that the Kbhb sites of key enzymes of energy metabolism, for instance, mitochondrial fatty acid β-oxidation and TCA cycle, were significantly downregulated after AMPKα2 knockout.
The current data indicated that AMPKα2 knockout did not elicit significant changes in appearance, body weight, grip and heart weight. (Fig. 1C, 1 D), indicating that specific knock out of AMPKα2 did not affect the general condition in young mice. We also observed that AMPKα2 knockout provokes unfavorable changes in the myocardium, along with the ventricular diastolic function impaired, in line with ultrastructural disarray and mitochondria distributed randomly (Fig. 2, Table 2). Notably, AMPK knockout induced flattened and disappeared mitochondrial crest in mouse cardiac tissue, suggesting that the knockout affects energy metabolism, which is the primary event in the mitochondria.
PTMs are crucial for protein function (activity and stability) and are closely related to various cardiovascular diseases [
]. Hence, we performed proteomics and Kbhb proteomics analyses to elucidate the potential molecular mechanisms of AMPK energy regulators. Quantitative Kbhb proteomics analysis discovered 244 Kbhb significantly altered protein sites, of which 42.6% (104/244) were upregulated and 57.4% (144/244) were downregulated after AMPKα2 knockout (Fig.5C). The Kbhb modified proteins were mainly distributed in the mitochondria (35.45%). Further analysis showed that pathway were enriched for Kbhb, including macronutrient, energy metabolism, and congenital heart disease. Our data investigated the functions and roles in AMPKα2 knockout progression by GO and KEGG pathway enrichment analyses. Several nutrient pathways were enriched for Kbhb-containing glycolysis, amino acid metabolism, fatty acid β-oxidation, and the TCA cycle.
The richest KB, β-hydroxybutyrate, has diverse bioactive properties, not only as an energy source but also as an inhibitor of histone deacetylases (HDACs). Shimazu et al. showed that ketogenic diet or fasting and other conditions connected with increasing β-OHB abundance, and corresponding increase in overall histone acetylation in mouse tissues, is similar to the Kbhb-induced conditions [
]. KBs, especially β-OHB, represent a transported form of acetyl-CoA. In mouse myocardial mitochondria, acetyl-CoA could be formed from pyruvate by PDH and fatty acids through β-oxidation [
]. Acetyl-CoA levels, together with acetyltransferase and deacetylase, adjust protein in acetylation through mitochondrial fatty acid β-oxidation and TCA cycle pathways. As a regulatory center of energy metabolism, AMPK regulates a variety of metabolic enzymes. Previous studies have shown that the activation of the AMPK/ACC/CPT-1 pathway promotes fatty acid β-oxidation and increases the yield of β-OHB (KB), which in turn increases the overall histone acetylation in mouse tissues[
]. According to Koronowski et al., as the concentration of β-OHB increased, its activated CoA formed β-OHB-CoA, which acted as the substrate of Kbhb, alter the overall Kbhb level. AMPK facilitates the transition from pyruvate metabolism to the TCA cycle by tight regulation of PDH activity [
]. In our findings, the Kbhb sites of key enzymes in β-oxidation and TCA cycle (for instance CPT-1 and PDHA1) were significantly downregulated after AMPKα2 knockout, which appears to overlap significantly with the acetylation target pathways. Rardin et al. showed that the proteins and sites of kac and Kbhb overlapped 75% when mitochondrial analysis was limited [
]. Since AMPK has been linked to abnormal metabolic pathway activation and epigenetic adjustment, the characteristics of carcinoma, affecting AMPK could provide novel targets for carcinoma therapy [
Previously, it was thought that tumor tissue could bypass the TCA cycle and exploit glycolysis (Warburg effect); however, recent evidence suggested that a large number of cancer cells depend on the TCA cycle for supply energy [
]. Zhen et al. demonstrated that AMPK can convert the metabolic processes to the TCA cycle by regulating PDH activity, thereby preventing oxidative stress and metabolism-induced cell death, which ultimately leads to cancer cell metastasis [
]. The present study showed that after AMPKα2 knockout, the site of the key enzyme PDH in the TCA cycle was downregulated by 3.73-fold. Recent studies have reported that AMPK sustains the activity of the PDH complex (PDHc) by phosphorylating PDHA [
]. Whether AMPK can directly affect the activity of PDHc by mediating the Kbhb PDHA, thereby regulating the TCA cycle, requires further investigation.
Taken together, our study is the first analysis of Kbhb substrates post-AMPKα2 knockout and provides a dataset of cardiac Kbhb in mammals. Further studies on the function of protein Kbhb modifications in energy metabolism will provide an in-depth insight into protein PTMs, discuss the distinctive characteristic of AMPK in regulating energy metabolism, and provide a clue for future drug exploits for AMPK-targeted cancer therapy.
For the first time, Kbhb modification was reported in the cardiac tissue of AMPKα2 knockout mice; nevertheless, the present study had some limitations. Although a considerable number of metabolites and Kbhb peptides have been identified in this study, the exact modification site function remains to be elucidated.
Data availability
Annotated spectra for identified β-hydroxybutyrylation modification have been deposited on MS-viewer (https://msviewer.ucsf.edu/cgi-bin/msform.cgi?form=msviewer) with search key "vau9y04iwv". Mass spectrometry proteomics data have been deposited on the ProteomeXchange Consortium via the iProX partner repository. Project name: The role of β-hydroxybutyrylation in cardiac energy metabolism using AMPKα2 knockout mouse model. Dataset identifier: PXD034269.
Conflict of Interest
The authors declare no competing interests.
Acknowledgments
This work was supported by the Key Technology Research and Development Program of Shandong [2022CXGC010510]; the Natural Science Foundation of Shandong Province [ZR2020MH313]; the Key Technology Research and Development Program of Shandong [2017GSF218101]; the National Natural Science Foundation of China [81700725]; the National Major Science and Technology projects of China [2012ZX09303016-003]; the Natural Science Foundation of Shandong Province [ZR2017BH003].
Trimetazidine Ameliorates Myocardial Metabolic Remodeling in Isoproterenol-Induced Rats Through Regulating Ketone Body Metabolism via Activating AMPK and PPAR alpha.
Wen-jing Ding: Data curation, Software, Validation, Writing-Original draft preparation; Xue-hui Li: Data curation, Software, Validation, Editing; Cong-min Tang: Visualization, Software; Xue-chun Yang: Investigation, Software; Yan Sun: Data curation, Investigation; Yi-ping Song: Data curation, Visualization; Ming-ying Ling: Investigation, Visualization; Rong Yan: Software, Data curation; Hai-qing Gao: Data curation; Wen-hua Zhang: Resources; Na Yu: Software; Jun-chao Feng: Visualization; Zhen Zhang: Methodology, Software; Yan-qiu Xing: Supervision, Conceptualization, Methodology, Software.
In Brief
Quantitative β-hydroxybutyrylation modification omics was used to detect AMPKα2 knockout mice cardiac tissue. The β-hydroxybutyrylation sites of myocardial fatty acid degradation and TCA cycle associated enzymes were dramatically downregulated after AMPKα2 knockout. Protein β-hydroxybutyrylation may have a role in AMPK-mediated energy metabolism regulation, providing a clue for future drug exploits for AMPK-related therapy.