Myosin dynamics and myosin ATP consumption is unchanged in Ursus arctos and Ursus americanus during hibernation.

A. Representative fluorescence mant-ATP decays from single muscle fibers isolated from Ursus arctos skeletal muscle measured over 300 seconds. B-C. Percentage of myosin heads in the P1/DRX (B) or P2/SRX (C) from Ursus arctos single muscle fibers obtained during summer (active) or winter (hibernating) periods. Values were separated based on each individual fiber was MyHC type I or MyHC type II. D. T1 value in seconds denoting the ATP turnover lifetime of the DRX. E. T2 value in seconds denoting the ATP turnover lifetime in seconds of the SRX. F. Calculated myosin ATP consumption values of each single muscle fiber per minute. This was calculated using the equation shown in the materials and methods section. G. Representative fluorescence mant-ATP decays from single muscle fibers isolated from Ursus americanus skeletal muscle measured over 300 seconds. H-I. Percentage of myosin heads in the P1/DRX (G) or P2/SRX (H) from Ursus americanus single muscle fibers obtained during summer (active) or winter (hibernating) periods. Values were separated based on each individual fiber was MyHC type I or MyHC type II. J. T1 value in seconds denoting the ATP turnover lifetime of the DRX. K. T2 value in seconds denoting the ATP turnover lifetime in seconds of the SRX. L. Calculated myosin ATP consumption values of each single muscle fiber per minute. Grey circles represent the values from each individual muscle fiber which was analyzed. Colored triangles represent the mean value from an individual animal, 8-12 fibers analyzed per animal. Statistical analysis was performed upon mean values. One-way ANOVA was used for statistical testing. n = 5 individual animals per group.

Myosin ATP turnover lifetime is reduced during hibernation in small hibernators, Eliomys quercinus & Ictidomys tridecemlineatus, resulting in an increase in myosin ATP consumption at ambient temperatures.

A-B. Percentage of myosin heads in the P1/DRX (A) or P2/SRX (B) from E. quercinus single muscle fibers obtained during active, interbout arousal (IBA) or torpor periods. Values were separated based on each individual fiber was MyHC type I or MyHC type II. C. T1 value in seconds denoting the ATP turnover lifetime of the DRX in E. quercinus. D. T2 value in seconds denoting the ATP turnover lifetime in seconds of the SRX in E. quercinus. E. Calculated myosin ATP consumption values of each single muscle fiber per minute in E. quercinus. This was calculated using the equation shown in the materials and methods section. F-G. Percentage of myosin heads in the P1/DRX (F) or P2/SRX (G) from Ictidomys tridecemlineatus single muscle fibers obtained during summer active (SA), interbout arousal (IBA) or torpor periods. H. T1 value in seconds denoting the ATP turnover lifetime of the DRX in Ictidomys tridecemlineatus. I. T2 value in seconds denoting the ATP turnover lifetime in seconds of the SRX in Ictidomys tridecemlineatus. J. Calculated myosin ATP consumption values of each single muscle fiber per minute in Ictidomys tridecemlineatus. Grey circles represent the values from each individual muscle fiber which was analyzed. Colored triangles represent the mean value from an individual animal, 8-12 fibers analyzed per animal. Statistical analysis was performed upon mean values. One-way ANOVA was used to calculate statistical significance. * = p < 0.05, ** = p < 0.01. n = 5 individual animals per group.

X-ray diffraction experiments of skeletal muscle from Ictidomys tridecemlineaus demonstrate changes in M6 myosin meridional spacing during torpor.

A. Representative x-ray diffraction recordings from permeabilized skeletal muscle bundles from Ictidomys tridecemlineatus from summer active (SA), interbout arousal (IBA) and torpor. The M3 and M6 meridional reflections and the 1,0 and 1,1 equatorial reflections are indicated. B. Ratio of the 1,1 to 1,0 equatorial reflections from active, IBA and torpor skeletal muscle. C. M3 meridional spacing, measured in nm. D. Normalized intensity (A.U.) of the M3 meridional reflection. E. M6 meridional spacing, measured in nm. Colored circles represent the mean value obtained from each skeletal muscle bundle which was recorded. Data is displayed as mean ± SEM. One-way ANOVA was used to calculate statistical significance. * = p < 0.05, ** = p < 0.01, *** = p < 0.001. n = 5 individual animals per group.

Myosin dynamics of Ictidomys tridecemlineatus are protected from temperature induced change during torpor, preventing an increase in myosin ATP consumption.

A. T1 value in seconds denoting the ATP turnover lifetime of the DRX in I. tridecemlineatus at 8°C. B. T2 value in seconds denoting the ATP turnover lifetime in seconds of the SRX in I. tridecemlineatus at 8°C. C. Ratio of the T1 expressed as the mean value for each matched animal at 20°C/8°C, separated for fiber type. D. Ratio of the T2 expressed as the mean value for each matched animal at 20°C/8°C, separated for fiber type. E. Ratio of calculated myosin ATP consumption expressed as 20°C/8°C, separated for fiber type. Black triangles represent the mean ratio value for each animal. One-way ANOVA was used to calculate statistical significance. * = p < 0.05, ** = p < 0.01, *** = p < 0.001. n = 5 individual animals per group.

MYH2 protein in Ictidomys tridecemlineatus is hyper-phosphorylated during torpor, which is predicted to increase protein stability.

A. Peptide mapping of differentiated phosphorylation sites upon MYH2 protein during SA, IBA and torpor periods. Heat map demonstrates all sites observed to be differentiated following the calculation of z-scores for each site. Z-scores > 0 equal hyper-phosphorylation and z-scores < 0 equal hypo-phosphorylation for each residue. Violin plot demonstrates significantly differentiated residues using z-scores. Two-way ANOVA with Šídák’s multiple comparisons test was used to calculate statistical significance. * = p < 0.05, ** = p < 0.01, *** = p < 0.001, **** = p < 0.0001. n = 5 individual animals per group. B. Chimera of MYH2 protein created using ChimeraX software. Important regions of the protein are annotated including coiled-coil region, ATP binding domain, actin binding domain and N-terminal SH3-like domain. Also, significantly hyper-phosphorylated residues are highlighted in red. C. Schematic of MYH2 protein with regions and hyper-phosphorylated resides annotated in red. Figure made in BioRender. D. EvoEF calculations of protein stability in both wild type and phosphor-mimetic mutants. Aspartic acid was used to mimic phospho-threonine/phospho-serine due to their chemical similarity. ΔGStability indicates the stability score for the protein in its corresponding configuration. ΔΔGStability represents the change in stability in mutant proteins versus the wild type protein. ΔΔGStability of > 0 represents an increase in the stability of a mutant versus wild type.

Global proteome analysis demonstrates changes to metabolic and sarcomeric changes in skeletal muscle fibers from Ictidomys tridecemlineatus during IBA and torpor.

A. Principal component analysis for all animals analyzed during SA, IBA and torpor periods. B. Volcano plot displaying proteins which are differentially expressed during torpor vs active periods. FDR < 0.01. Red circles are upregulated proteins and blue circles are downregulated proteins. Highly differentiated proteins of interest are annotated with their respective protein name. C. Volcano plot displaying proteins which are differentially expressed during IBA vs SA periods. Red circles are upregulated proteins and blue circles are downregulated proteins. Highly differentiated proteins of interest are annotated with their respective protein name. D. Ontological associations between proteins upregulated during torpor vs SA periods. The top five association clusters are annotated on the network. A full list of clusters and the proteins lists included in clusters are available in supplementary figure 1 and supplementary table 2. E. Ontological associations between proteins downregulated during torpor vs SA periods. The top five association clusters are annotated on the network. A full list of clusters and the proteins lists included in clusters are available in supplementary figure 1 and supplementary table 2. F. Ontological associations between proteins upregulated during IBA vs SA periods. The top five association clusters are annotated on the network. A full list of clusters and the proteins lists included in clusters are available in supplementary figure 2 and supplementary table 3. G. Ontological associations between proteins downregulated during IBA vs SA periods. The top five association clusters are annotated on the network. A full list of clusters and the proteins lists included in clusters are available in supplementary figure 2 and supplementary table 3. Gene ontology networks were established using Metascape and visualized using Cytoscape. Detailed information upon the statistical testing used is available in the methods section. FDR < 0.01 significantly differentially expressed proteins were used to establish networks. Purple lines indicate a direct interaction. Circle size is determined by enrichment and color is determined by p value. n = 5 individual animals per group.