Nemaline myopathies are the most common form of congenital myopathies. Variants in ACTA1 (NEM3) comprise 15–25% of all nemaline myopathy cases. Patients harboring variants in ACTA1 present with a heterogeneous disease course characterized by stable or progressive muscle weakness and, in severe cases, respiratory failure and death. To date, no specific treatments are available. Since NEM3 is an actin-based thin filament disease, we tested the ability of tirasemtiv, a fast skeletal muscle troponin activator, to improve skeletal muscle function in a mouse model of NEM3, harboring the patient-based p.Asp286Gly variant in Acta1. Acute and long-term tirasemtiv treatment significantly increased muscle contractile capacity at submaximal stimulation frequencies in both fast-twitch extensor digitorum longus and gastrocnemius muscle, and intermediate-twitch diaphragm muscle in vitro and in vivo. Additionally, long-term tirasemtiv treatment in NEM3 mice resulted in a decreased respiratory rate with preserved minute volume, suggesting more efficient respiration. Altogether, our data support the therapeutic potential of fast skeletal muscle troponin activators in alleviating skeletal muscle weakness in a mouse model of NEM3 caused by the Acta1:p.Asp286Gly variant.
Introduction
Nemaline myopathies (NEM) are among the most common forms of congenital myopathies (Colombo et al., 2015) and compose a group of genetically heterogeneous disorders characterized by skeletal muscle weakness (Wallgren-Pettersson et al., 2011). NEM disease severity ranges from neonatally severe to a mildly progressive phenotype (Sanoudou and Beggs, 2001; Ryan et al., 2001). Histopathologically, NEM patient’s muscle biopsies are characterized by the accumulation of nemaline bodies, a fiber type shift toward type 1 fibers, and in severe cases, disarrangement of the myofibrillar architecture (Yin et al., 2014). To date, 13 genes have been implicated in NEM. α-Skeletal-actin 1 (ACTA1) (Ilkovski et al., 2001); α- and β-tropomyosin (TPM3 and TPM2) (Donner et al., 2002; Laing et al., 1995); nebulin (NEB) (Pelin et al., 1999); leiomodin-3 (LMOD3) (Yuen et al., 2015); troponin T, slow skeletal type (TNNT1) (Johnston et al., 2000), and Troponin T3, fast skeletal type (TNNT3) (Sandaradura et al., 2018); cofilin 2 (CFL2) (Agrawal et al., 2007); unconventional myosin 18B (MYO18B) (Malfatti et al., 2015); myopalladin (MYPN) (Miyatake et al., 2017); Kelch family members 40 (KLHL40) and 41 (KLHL41) (Gupta et al., 2013; Ravenscroft et al., 2013); and Kelch repeat and BTB (POZ) domain containing 13 (KBTBD13) (Sambuughin et al., 2010). 12 of these genes code for proteins that are associated with the thin filament, a key component of the skeletal muscle sarcomere. Variants in these genes result in sarcomeric dysfunction (Ochala, 2008; de Winter and Ottenheijm, 2017; Ottenheijm et al., 2010; Ochala et al., 2010, Winter et al., 2016, de Winter et al., 2020). No specific treatments are currently available.
Nemaline myopathy type 3 (NEM3) is caused by variants in the skeletal muscle α-actin gene (ACTA1) and comprises the second most common form of NEM (15–25% of all cases [Wallgren-Pettersson et al., 2011]). In vitro studies showed that variants in ACTA1 affect thin filament structure and function (Vandamme et al., 2009; Marston et al., 2004), and the variants contribute to contractile weakness as assessed in vitro and in vivo (Gineste et al., 2013b; Ochala et al., 2015; Ochala et al., 2012; Lindqvist et al., 2013; Ravenscroft et al., 2011a; Joureau et al., 2018; Nguyen et al., 2011). Thus, restoring thin filament function may prove as a useful therapy to ameliorate muscle weakness in NEM3 patients.
Tirasemtiv (CK-2017357) is a fast skeletal muscle troponin activator, which increases the binding of calcium to the troponin complex (Li et al., 2021) and augments the force-generating capacity at submaximal nerve stimulations in rodents and humans (Russell et al., 2012; Hansen et al., 2014; Andrews et al., 2018). In vitro treatment with tirasemtiv has been shown to increase muscle strength in NEM2 (caused by variants in NEB) and in other neuromuscular disorders (Lee et al., 2019; Hwee et al., 2014; van de Locht et al., 2021; de Winter et al., 2013). Recently, we investigated the acute and long-term effect of tirasemtiv in a mouse model for severe NEM3 due to the heterozygous NM_009606.3 (Acta1):p.His40Tyr variant in Acta1 (referred to as Acta1(H40Y)). Tirasemtiv had a positive acute and long-term effect on muscle force with a reduction in the energetic cost of contraction (de Winter et al., 2021).
It is unknown whether tirasemtiv has a similar positive inotropic effect on muscles in milder forms of NEM3, which have a smaller therapeutic window. Therefore, in this study, we tested the in vitro and in vivo effect of tirasemtiv on skeletal muscle function in a mouse model that recapitulates a mild form of NEM3 caused by the patient-based (Acta1):p.Asp286Gly variant in Acta1 (Joureau et al., 2018). We hypothesized that both acute and long-term tirasemtiv treatment would increase skeletal muscle force at submaximal stimulation frequencies. Furthermore, we studied protein oxidation and pathways controlling proteostasis, e.g., ubiquitin-proteosome pathway, and mitochondrial dynamics to assess whether they were altered in this variant of NEM3 and could be involved in tirasemtiv impact on skeletal muscle contractility.
Throughout the manuscript, variant (Acta1):p.Asp286Gly is referred to as Acta1(D286G) to be consistent with existing literature (Ravenscroft et al., 2011a). Our results show that Acta1(D286G) mice had significantly lower hind limb muscle mass as compared to wild type mice. Acute and long-term tirasemtiv treatment significantly increased in vitro force production at submaximal stimulation frequencies in extensor digitorum longus (EDL) and diaphragm and reduced fatigability in the EDL muscle. Long-term tirasemtiv treatment showed improved in vivo contractility of gastrocnemius (GAS) muscle at submaximal stimulation frequencies. Myofibrillar proteins and myosin oxidation changes could not explain the impact of tirasemtiv on contractility. Together, these findings highlight the therapeutic potential of fast skeletal troponin activators to alleviate skeletal muscle weakness in mild forms of NEM3.
Materials and methods
Acta1(D286G) mouse model
Transgenic 8–12 mo-old male mice expressing the ACTA1 protein harboring the D286G variant p.Asp286Gly mutation in Acta1 (referred to as Acta1(D286G)) and wild type littermates were used for the experiments (Ravenscroft et al., 2011b). Experiments were conducted in agreement with the French and Dutch guidelines for animal care. All animal experiments were approved by the Institutional Animal Care Committee of Aix-Marseille University (#15–14052012) and by the local animal ethics committee at VU University (AVD114002016501). Mice were housed in an environment-controlled facility (12–12 h light-dark cycle, 22°C), and received water and standard food ad libitum. Mice were identified through PCR genotyping from mouse tail DNA. All mice survived the study and had a normal life span.
Tirasemtiv treatment
For the experiments in which the acute effect of tirasemtiv was studied, mice were I.P. injected with vehicle or 3 mg/kg tirasemtiv (Hwee et al., 2014). Experiments were performed ∼30 min after injection. For the experiments in which muscle contractility was studied in vitro, 3 µM tirasemtiv was added to the experimental solutions. For the studies in which the long-term effects of tirasemtiv were studied, 8-mo-old mice were first fed for 1 wk with custom-made mouse pellets (BioServ). After 1 wk, mice were switched to the same pellets containing tirasemtiv (600 ppm) or the same pellets without tirasemtiv. Mice were kept on the chow for 4 wk. Whole-muscle contractility and respiratory function were tested before treatment and after 4 wk of tirasemtiv-enriched diet or regular diet. After 4 wk, mice were euthanized and tissues were collected, i.e., tibialis anterior (TA), EDL, soleus (SOL), GAS, and diaphragm. Blood was collected to determine plasma concentrations of tirasemtiv, but we did not measure its concentration in skeletal muscles. Therefore, we do not know tirasemtiv exposure levels in muscles during the long-term treatment.
In vitro EDL and diaphragm muscle characterization
In vitro characterization of intact muscle was performed as described previously (de Winter et al., 2020). The experimental protocols consisted of a full tetanus at 150 Hz (diaphragm) and 200 Hz (EDL) and a force-frequency protocol. For the force-frequency protocol, the muscle was stimulated with incremental stimulation frequencies (diaphragm: 1, 5, 10, 20, 40, 60, 80, 100, and 150 Hz; and EDL: 1, 5, 10, 20, 40, 60, 80, 100, 150, and 200 Hz). Data were discarded when stimulation at 150 Hz (diaphragm) or 200 Hz (EDL) rendered a force that was <90% of the force generated during the first maximal tetanic stimulation. Stimuli were applied with a train duration of 600 ms. The resting interval was 30 s between the stimulations at 1 and 10 Hz; 60 s after stimulation at 20 Hz; 90 s after stimulation at 30 Hz; and 120 s between stimulations at 60, 80, 100, 150, and 200 Hz. Then a fatigue protocol (80 contractions; 40 Hz; 2 s on and 1 s off) was performed. The peak force of each contraction was measured and normalized to muscle mass. A fatigue index corresponding to the ratio between the last five and the first five contractions was determined. After completion of the contractility measurements, the length and mass of the muscles were determined. Cross-sectional area (CSA; in mm2) was calculated by dividing muscle mass (g) by muscle length (cm) multiplied by specific density (1.056 g/ml) × 100.
In vivo plantar flexor muscle measurements
Animal preparation: Mice were anesthetized and individually placed supine in a home-built cradle specially designed for the strictly non-invasive functional investigation of the left hind limb muscles as described previously (Gineste et al., 2013a). Force output measurements: Non-invasive transcutaneous electrical stimulation was first elicited with square-wave pulses (0.5-ms duration) on the plantar flexor muscles. The individual maximal stimulation intensity was determined by progressively increasing the stimulus intensity until there was no further peak twitch force increase. Plantar flexion force was assessed in response to stimulation trains delivered at 20 and 100 Hz (train duration = 750 ms). We chose those specific stimulation frequencies in vivo to match the relative force production we observed in vitro at 40 and 150 Hz. The resulting force was divided by the sum of GAS/plantaris + SOL muscle mass (in mN/mg). Then, a fatigue protocol (80 stimulations; 40 Hz; 1.5 s on and 6 s off) was performed. The peak force of each contraction was measured, normalized to muscle mass, and the corresponding values were averaged every five contractions. A fatigue index corresponding to the ratio between the last five and the first five contractions was determined.
Muscle fiber CSA analysis
Muscle fiber CSA was determined in the mid-belly region of GAS and EDL muscles and in a portion of the diaphragm muscle as previously described (Gineste et al., 2013b). Briefly, muscle serial transverse sections (10-μm thick) were obtained from each muscle and were immunostained with monoclonal antibodies against myosin heavy chain (MHC) isoforms (BA-F8 against MHC-1 and SC-71 against MHC-2A). The cryosections were incubated with a primary antibody for 1 h at 37°C, rinsed with PBS buffer, and incubated in a secondary rabbit anti-mouse antibody conjugated with peroxidase (DAKO) for 1 h at 37°C. After washing in PBS buffer, the stain was visualized by using a DAB (3,3′-diaminobenzidine) solution. Fibers negative for the two antibodies were considered type 2X or 2B. Images of the stained sections were captured from a light microscope (LeicaDMLS; ocular magnification: 10×, objective magnification: 10×, the numerical aperture of the objective lenses: 0.22) equipped with a camera (DFC450C; Leica), using the software LASV3.8. Fiber CSA was measured using the ImageJ 1.5 analysis software (Schneider et al., 2012) (NIH) and expressed in micrometers squared. Given the low number of fibers positive for antibodies against MHC-1 and MHC-2A, both in EDL and GAS muscles, CSA was assessed on the negative fibers (type 2X and 2B), which were the most prominent ones in both muscles. In the diaphragm muscle, CSA was measured without differentiation between fiber types.
Plethysmography analysis
Mice were placed in an unrestrained whole-body plethysmography chamber for 30 min for acclimatization. After acclimatization, tidal volume (VT), respiratory rate (RR), and minute ventilation (MV) were monitored for 15 min at room air. After 15 min, mice were exposed to a 5% CO2 gas mixture for 30 min and monitored. After the 5% CO2 exposure, mice were reexposed to room air for 15 min and monitored (de Winter et al., 2021). The data were analyzed by Buxco Finepointe Respirometry and Inhalation software (Data Science International). FinePointe software was used to acquire and analyze the respirometry data that were collected using the Buxco Small Animal Whole Body Plethysmography system. RR and VT were provided by the Finepointe software. MV was calculated as the product of RR and VT.
MHC composition analysis
Frozen EDL and GAS muscles were pulverized in a steel mortar with liquid nitrogen to obtain a powder that was immediately resuspended in Laemmli solution (Laemmli, 1970). The samples were incubated in ice for 20 min and finally spun at 18,000 g for 30 min. Protein concentration in the dissolved samples was determined with a protein assay kit (RC DC Biorad). About 10 µg of proteins for each sample was loaded on 8% SDS-PAGE polyacrylamide gels and the electrophoresis was run overnight at 250 V. Following a Coomassie stain, four bands corresponding to MHC isoforms were separated and their densitometric analysis was performed to assess the relative percentage of MHC isoforms (MHC-1, MHC-2A, MHC-2B, and MHC-2X) in the samples (Pellegrino et al., 2003).
Western blot analysis
Frozen GAS muscle samples were pulverized and immediately resuspended in a lysis buffer (20 mM Tris-HCl, 1% Triton X100, 10% glycerol, 150 mM NaCl, 5 mM EDTA, 100 mM NaF, and 2 mM NaPPi supplemented with 1× protease, phosphatase inhibitors [Sigma-Aldrich] and 1 mM PMSF). The homogenate obtained was kept on ice for 20 min and then centrifuged at 18,000 g for 20 min at 4°C. The supernatant was stored at −80°C until ready to use. Protein concentration was evaluated for each sample and equal amounts of muscle samples were loaded on gradient precast gels purchased from Bio-Rad (AnyKd). After the gel run, proteins were electro-transferred to PVDF membranes at 35 mA overnight. The membranes were incubated in 5% milk for 2 h, rinsed with TBST buffer (0.02 M Tris and 0.05 M NaCl, pH 7.4–7.6), and subsequently probed with specific primary antibodies (see below). Thereafter, the membranes were incubated in HRP-conjugated secondary antibodies. The protein bands were visualized by an enhanced chemiluminescence method in which luminol was excited by peroxidase in the presence of H2O2 (ECL Select, GE Healthcare). The content of each protein investigated was assessed by determining the brightness–area product of the protein band (Cannavino et al., 2015). Antibodies used were as follows: anti-rabbit GAPDH (catalog ab9484, 1:2,000; Abcam); anti-rabbit Catalase (catalog ab52477, 1:1,000; Abcam); anti-mouse SOD1 (catalog ab16831, 1:1,000; Abcam); anti-rabbit FIS1 (catalog ab71498, 1:1,000; Abcam); anti-rabbit p-DRP1(ser616) (catalog #3455, 1:1,000; Cell Signalling); anti-rabbit p-DRP1(ser637) (catalog #4867, 1:1,000; Cell Signalling); anti-rabbit DRP1 (catalog #8570, 1:3,000; Cell Signalling); anti-mouse MFN1 (catalog ab57602, 1:1,000; Abcam); anti-rabbit MFN2 (catalog ab50843, 1:1,000; Abcam); anti-rabbit OPA1 (catalog ab157457, 1:3,000; Abcam); anti-mouse IgG (catalog P026, 1:5,000; Dako North America Inc.); and anti-rabbit IgG (catalog #7074, 1:10,000; Cell Signalling).
Carbonylated protein analysis
Frozen GAS samples from each subject group were suspended in a lysis buffer (50 mM Tris-HCl pH 7.6, 250 mM NaCl, 5 mM EDTA protease inhibitor cocktail, and phosphatase inhibitor cocktail), left on ice for 20 min, and finally centrifuged at 18,000 g for 20 min at 4°C. Protein concentration was determined using the RC DC TM protein assay kit (Biorad product). The protein carbonylation level was detected using the OxyBlot TM Kit (Millipore), which provides reagents for sensitive immunodetection of carbonyl groups. Carbonyl groups in the protein side chains were derivatized to 2,4-dinitrophenylhydrazone (DNP) by reaction with 2,4-dinitrophenylhydrazine (DNPH). In detail, 10 μg of proteins for each muscle sample was denatured with SDS solution at a final concentration of 6%. DNPH solution was added to obtain the derivation; the reaction was stopped using a neutralization solution (90449; Millipore) after 10 min of incubation at room temperature. The DNP-derivatized protein samples were separated by polyacrylamide gel electrophoresis (Anykd Biorad gels) followed by Western blotting. Proteins were transferred to nitrocellulose membranes at 100 V for 2 h, stained with Ponceau Red (Sigma-Aldrich), and then scanned. The membranes were blocked by incubation with 3% bovine serum albumin for 1 h and then incubated with rabbit anti-DNP antibody overnight at 4°C and subsequently with a horseradish peroxidase-antibody conjugate (goat anti-rabbit IgG). The positive bands were visualized by using a chemiluminescent reagent (ECL advance as described previously [Cannavino et al., 2015]). The total protein carbonylation level and the MHC carbonylation level were analyzed quantitatively by comparison of the signal intensity of immune-positive proteins normalized on the total protein amount loaded on gels (ponceau staining signal) (Cannavino et al., 2015).
Gene expression analysis
Total RNA was extracted from GAS muscles using an SV Total RNA isolation kit (Promega). The RNA concentration was measured using a Nano Drop instrument (Thermo Fisher Scientific) and 400 ng was used to generate cDNA with SuperScript III reverse transcriptase (Invitrogen). The cDNA was analyzed by quantitative RT-PCR (AB7500; Applied Biosystems) using a SYBR Green PCR kit (Applied Biosystems), and the data were normalized to Gapdh content. Oligonucleotide primers were provided by Sigma-Aldrich and were MuRF-1 (FP: 5′-ACCTGCTGGTGGAAAACATC-3′, RP: 5′-ACCTGCTGGTGGAAAACATC-3′) and Atrogin-1 (FP: 5′-GCAAACACTGCCACATTCTCTC-3′, RP: 5′-CTTGAGGGGAAAGTGAGACG-3′). Differentially expressed genes were determined using a default threshold of 0.6. The difference between Ct (cycle threshold) values was calculated for each mRNA by taking the mean Ct of duplicate reactions and subtracting the mean Ct of duplicate reactions for the reference RNA measured on an aliquot from the same RT reaction (ΔCt = Ct target gene − Ct reference gene). All samples were then normalized to the ΔCt value of a calibrator sample to obtain a ΔΔCt value (ΔCt target − ΔCt calibrator) (comparative method) (Cannavino et al., 2015).
Statistics
Data are presented as mean ± standard error of the mean. For mouse model characterization and long-term tirasemtiv effect data analysis, an unpaired t test or two-way ANOVA followed by Sidak’s multiple comparison test (for force-frequency data and baseline MHC composition) was used. For acute tirasemtiv effect data analysis, paired t test or repeated measures two-way ANOVA followed by Sidak’s multiple comparison test (for force-frequency data) was used. All data analysis was performed with GraphPad Prism 9 software (version 9.5.1 [733]). A probability value <0.05 was considered statistically significant.
Online supplemental material
Fig. S1 describes the effect of long-term tirasemtiv administration on protein expression of muscle mitochondrial dynamic markers. All mitochondrial markers were measured in GAS muscles of wild type and Acta1(D286G) tirasemtiv- and vehicle-treated mice.
Results
Characterization of Acta1(D286G) mouse skeletal muscles
In line with previous work (Ravenscroft et al., 2011a), homozygous Acta1(D286G) mice had lower whole-body weight as well as lower hind limb muscle mass when compared with wild type mice (Fig. 1 A). Since tirasemtiv targets exclusively fast-twitch muscle fibers, we determined the morphological and functional characteristics of the fast glycolytic EDL and GAS muscles. Both EDL and GAS muscles of Acta1(D286G) mice had a significantly higher proportion of MHC type IIx and a significantly lower proportion of type MHC IIb compared with wild type mice. The proportion of MHC type I and IIa was comparable between groups (Fig. 1 B, top and bottom panels). Representative images of EDL and GAS MHC fiber type using SDS-PAGE polyacrylamide gels are shown in Fig. 1 C, bottom panel.
Muscle atrophy is often observed in NEM. Since we observed lower muscle mass in TA, SOL, EDL, and GAS in Acta1(D286G) mice, we measured EDL physiological CSA and CSA of the most prominent fiber types, 2X and 2B, of both EDL and GAS muscle to assess atrophy. Whole EDL CSA and single fiber CSA of EDL and GAS were significantly lower in Acta1(D286G) mice (Fig. 1 D). Since muscle atrophy is usually associated with a decrease in force production, we measured the contractile parameters of EDL at incremental stimulation frequencies, which reflect muscle function during different daily life activities (Tikkanen et al., 2013; Jasmin and Gardiner, 1987). Relative forces of EDL at submaximal stimulation frequencies (60 Hz) were significantly lower in Acta1(D286G) mice (Fig. 1 D). As a consequence, the FF50, which represents the stimulation frequency required to generate 50% of maximal force, was significantly higher in Acta1(D286G) mice (Fig. 1 E; inset). The maximal force (200 Hz) of EDL was significantly lower in Acta1(D286G) mice; however, the maximal normalized force (i.e., tension) was comparable with wild type (Fig. 1 F). Thus, muscles of Acta1(D286G) mice have atrophy and a reduced maximal absolute force in addition to lower relative force in response to submaximal stimulation frequencies.
Acute effect of tirasemtiv on skeletal muscle function
We first investigated the acute effect of tirasemtiv treatment on EDL muscles in vitro. EDL muscles were isolated and mounted in an in vitro whole muscle intact mechanics system and electrically stimulated (schematic representation, Fig. 2 A). Wild type and Acta1(D286G) EDL muscles were exposed to 3 µM of tirasemtiv. This dose was chosen based on previous work that showed maximal contractility gain without impairing muscle relaxation kinetics (de Winter et al., 2021). Tirasemtiv exposure resulted in a leftward shift of the force-frequency curve with significantly higher relative tension at stimulation frequencies of 20, 40, 60, 80, 100, and 150 Hz in wild type and in Acta1(D286G) mice (Fig. 2 B; top and middle panels). The FF50 was significantly lower in wild type and Acta1(D286G) tirasemtiv-treated muscles (Fig. 2 B, top and middle panel insets). Fig. 2 B, bottom panel, shows tensions generated by EDL during the force-frequency protocol for all groups and experimental conditions. As anticipated, tirasemtiv did not affect the maximal tension in wild type and Acta1(D286G) EDL muscles (Fig. 2 C). Acute tirasemtiv forces are shown in Table 1. Furthermore, we assessed the effect of tirasemtiv on the fatigability of EDL. We repeatedly stimulated wild type and Acta1(D286G) EDL muscles (Fig. 2 D; top row) and determined the fatigue index. Acta1(D286G) EDL muscles exposed to tirasemtiv showed significantly higher resistance to fatigue development (Fig. 2 D; bottom right panel). No effect was observed on wild type muscles (Fig. 2 D, bottom left panel). Thus, acute tirasemtiv treatment increased force production at submaximal stimulation frequencies and increased the resistance to fatigue in Acta1(D286G) EDL muscles.
Effect of long-term tirasemtiv treatment on skeletal muscle morphology and contractility
To investigate the long-term effects of tirasemtiv on contractility of EDL (in vitro) and GAS (in vivo), wild type and Acta1(D286G) mice were fed tirasemtiv-enriched chow for 4 wk. After 4 wk of treatment, plasma levels of tirasemtiv ranged from 2,020 to 8,100 ng/ml. No tirasemtiv was detected in the plasma of mice fed standard chow (Fig. 3 A). Tirasemtiv treatment did not result in any significant body or muscle mass changes in wild type and Acta1(D286G) mice compared with vehicle-treated mice (Fig. 3 B). Additionally, long-term tirasemtiv treatment did not affect whole EDL CSA or GAS single fiber CSA (Fig. 3, C and D). We then assessed GAS fiber type composition. GAS of tirasemtiv-treated wild type mice had a significantly higher percentage of MHC type 2X fibers as compared to vehicle-treated mice (Fig. 3 E). We did not observe significant changes in fiber type composition in tirasemtiv-treated Acta1(D286G) mice. Since we observed a positive ionotropic effect of acute tirasemtiv treatment on EDL contractility, we next evaluated the long-term effect of tirasemtiv on in vitro EDL. Note that after EDL isolation and prior to the contractility assay, to wash out tirasemtiv and rule out acute effects on contractility, EDL was bathed for ∼20 min in tirasemtiv-free Ringer solution according to previous work (de Winter et al., 2021; Hwee et al., 2014). Similar to the observed acute effect of tirasemtiv in EDL muscles, long-term tirasemtiv treatment resulted in a leftward shift of the force-frequency curve with significantly higher relative tensions at stimulation frequencies of 40, 60, and 80 Hz in wild type and 20 and 40 Hz in Acta1(D286G) mice (Fig. 4 A; top and middle panel). The FF50 was significantly lower in the wild type but not in Acta1(D286G) tirasemtiv-treated muscles (Fig. 4 A; top and middle panel, insets). Fig. 4 A, bottom panel, shows tensions generated by EDL during the force-frequency protocol for all groups and experimental conditions. No effect of tirasemtiv was observed at maximal stimulation frequencies (200 Hz) in wild type mice and Acta1(D286G) EDL muscles (Fig. 4 B; left and right panel, accordingly). Furthermore, long-term tirasemtiv treatment did not have an effect on EDL fatigability in both wild type and Acta1(D286G) mice (Fig. 4 C). Long-term in vitro tirasemtiv force values are shown in Table 2.
Next, we evaluated the in vivo effect of long-term tirasemtiv administration on GAS contractile function. Whole GAS function was evaluated in vivo using an experimental setup previously described (Giannesini et al., 2010) (schematic representation, Fig. 4 D). Note that during these measurements, GAS muscles were not washed in tirasemtiv-free ringer solution as was done with EDL muscles. Similar to what we observed in EDL, at submaximal stimulation frequency, Acta1(D286G) tirasemtiv-treated mice produced significantly higher tension as compared to vehicle-treated mice (Fig. 4 E; top left panel). No significant differences were observed in tension at maximal stimulation frequency (Fig. 4 E; top right panel) or GAS response to fatigue after tirasemtiv treatment (Fig. 4 E, bottom left and right panel). The long-term in vivo tirasemtiv force values are shown in Table 3.
Thus, these data indicate that long-term tirasemtiv administration enhances submaximal muscle tension in both EDL and GAS muscles of Acta1(D286G) mice. The sustained in vitro effect of tirasemtiv on EDL muscles after removal by washing suggests that in addition to its direct effect, long-term tirasemtiv treatment may induce structural changes within the muscle that improve force generation.
Characterization of Acta1(D286G) mouse diaphragm muscle function
Respiratory dysfunction is often observed in NEM patients. Thus, we investigated the effect of tirasemtiv administration on respiratory function and diaphragm contractility in vitro and in vivo. First, we characterized diaphragm morphology and function in untreated wild type and Acta1(D286G) mice. In contrast to EDL and GAS, we did not observe significant differences in diaphragm fiber CSA between wild type and Acta1(D286G) mice (Fig. 5 A). We then characterized the MHC composition of diaphragm muscle strips. Acta1(D286G) and wild type mouse diaphragm were predominantly composed of MHC type 2A and 2X (∼90%). The diaphragm of Acta1(D286G) mice had a significantly lower percentage of MHC type I compared with wild type mice (1.2 versus 6.9%, respectively, Fig. 5 B). Next, we characterized diaphragm contractility in vitro. No significant changes were observed in the force-frequency relation and FF50 between both groups (Fig. 5 C). Maximal tension was significantly reduced in diaphragm strips of Acta1(D286G) mice compared with wild type (Fig. 5 D, left panel); however, unlike EDL and GAS, no changes were observed at submaximal stimulation frequencies (Fig. 5 D, right panel). We then characterized in vivo respiratory function in wild type and Acta1(D286G) mice using whole-body plethysmography. RR was comparable between groups (Fig. 5 E; left panel); however, tidal volume (VT; the volume of air moved in and out the lungs during one respiratory cycle) and minute volume (VE; the product of RR and VT), were significantly reduced in Acta1(D286G) mice (Fig. 5 E; middle and right panel). Thus, the diaphragm of Acta1(D286G) mice shows contractile weakness and impaired respiratory function.
Acute effect of tirasemtiv on respiratory function
Next, we tested the acute effect of tirasemtiv on diaphragm contractility in vitro and on respiratory function in vivo. Acute tirasemtiv treatment resulted in a leftward shift in the force-frequency relation with significantly higher force production at submaximal stimulation frequencies of 1, 10, 20, 40, and 60 Hz for wild type and 20 and 40 Hz for Acta1(D286G) mice (Fig. 6 A, top and bottom panel). FF50 was significantly reduced in both wild type and Acta1(D286G) diaphragm muscle strips (Fig. 6 A, top and bottom insets). The positive inotropic effect of tirasemtiv was not present at maximal stimulation frequencies (200 Hz) in wild type and Acta1(D286G) diaphragm muscle strips (Fig. 6 B, left and right panels). Acute tirasemtiv force values are shown in Table 4.
To assess the in vivo effect of tirasemtiv treatment on respiratory function, we performed whole-body plethysmography during acute and long-term tirasemtiv administration. Note that during these assays, tirasemtiv was present in the diaphragm muscles. Acute tirasemtiv treatment did not have any observable effects on RR, VT, or VE (Fig.6, C and D). Acute tirasemtiv plethysmography values are shown in Table 5. Thus, acute tirasemtiv exposure did not affect respiratory function in Acta1(D286G) mice.
Effect of long-term tirasemtiv treatment on respiratory function
Next, we assessed the effect of long-term tirasemtiv treatment on diaphragm muscle strips. We did not observe any changes in force production at any stimulation frequency nor differences in FF50 values (Fig. 7 A). Furthermore, no changes in tension were observed at maximal stimulation frequency in both groups (Fig. 7 B). Long-term tirasemtiv force values are shown in Table 6. Furthermore, long-term tirasemtiv treatment did not affect tidal volume and minute volume (Fig. 7, C and D, top and middle panels) but did significantly reduce the RR in Acta1(D286G) mice (Fig. 7 D, top panel). Long-term Plethysmography values are shown in Table 7. Thus, these data suggest that tirasemtiv improves ventilatory efficiency by maintaining VT and VE with fewer respiratory cycles.
Long-term effect of tirasemtiv on gene and protein expression and posttranslational modifications
We then assessed the effect of tirasemtiv on pathways involved in the regulation of muscle protein turnover. Muscle RING finger protein 1 (MuRF1) and atrogin-1 are muscle-specific E3 ligases that are involved in mediating muscle atrophy. We measured mRNA levels of both E3 ligases in GAS muscle after long-term tirasemtiv exposure. Tirasemtiv exposure significantly downregulated MurRF-1 and Atrogin-1 in wild type mice. In Acta1(D286G), tirasemtiv resulted in the upregulation of atrogin-1 (Fig. 8 A), with no changes in MuRF1 levels.
As dysfunctional mitochondria is a major source of reactive oxygen species, we aimed to assess whether tirasemtiv reduces the stress on such organelles. To this end, we assessed antioxidant levels in GAS muscle. No significant changes were observed in the protein levels of antioxidant enzymes superoxide dismutase 1 (SOD1) and Catalase (Fig. 8 B). Furthermore, we measured GAS total protein and MHC carbonylation. Such posttranslational modifications are implicated in skeletal muscle dysfunction in various chronic disorders (Barreiro et al., 2005a, 2005b; Li et al., 2015; Marin-Corral et al., 2009; Coirault et al., 2007). We observed that long-term tirasemtiv treatment had no effect on total protein carbonylation; however, it significantly reduced MHC carbonylation in both wild type and Acta1(D286G) mice (Fig. 8 C, bottom panel). Finally, we measured mitochondrial dynamic markers. Mitochondrial biogenesis, fission and fusion markers, such as peroxisome proliferator-activated receptor gamma coactivator 1-α (PGC-1α), Mitofusin-1 (MFN1), Mitofusin2 (MFN2), dynamin-like 120 kD protein (OPA1), dynamin-related protein 1 (DRP1), and mitochondrial fission 1 (FIS1) protein levels remain unchanged in both wild type and Acta1(D286G) mice after long-term tirasemtiv exposure (Fig. S1, A–C) for the exception of phosphorylated DRP1 (DRP1p616), which was significantly higher in tirasemtiv-treated wild type mice (Fig. 1 C, top panel).
Discussion
Here, we tested the ability of the fast skeletal muscle troponin activator tirasemtiv to improve muscle function in Acta1(D286G) mice. We show that acute and long-term tirasemtiv treatment significantly increases muscle contractile performance in vitro and in vivo at submaximal stimulation frequencies in both fast-twitch EDL and GAS and intermediate-twitch diaphragm muscles. The positive inotropic effect following acute and long-term tirasemtiv treatment is clinically relevant as most daily life activities, as well as tidal breathing, require submaximal muscle activation (Tikkanen et al., 2013; Jasmin and Gardiner, 1987). Altogether, our data support the potential positive therapeutic effect of fast skeletal muscle troponin activators in alleviating muscle weakness in a mouse model of NEM3.
Tirasemtiv improves skeletal muscle force-generating capacity and respiratory function in Acta1(D286G) mice
NEM is one of the most common forms of congenital myopathies, and, to date, no specific treatments are available. NEM causative genes encode proteins that are associated with the thin filament. Variants in thin filament genes such as ACTA1 (NEM3) are usually associated with a severe congenital form often characterized by profound skeletal muscle weakness (antenatal or neonatal), myofibrillar disarrangement, and, in severe cases, death due to respiratory insufficiency (Labasse et al., 2022). However, milder forms of NEM3 have also been reported, with a mild disease progression (Agrawal et al., 2004; Witting et al., 2016; Feng and Marston, 2009; Laing et al., 2009). NEM3 caused by the p.Asp286Gly (D286G) variant is associated with a milder form of the disease. Previous work has shown that this variant reduces the calcium sensitivity of force in mouse myofibers (Ravenscroft et al., 2011a) and directly affects the binding of myosin cross-bridges to actin, thereby reducing the fraction of acto-myosin interactions in the strong binding state and contractile force (Ochala et al., 2012). Thus, increasing thin filament function may be an attractive approach to augment muscle function.
Our results are in line with our previous work performed on a more severe NEM3 mouse model with an Acta1 p.His40Tyr variant (Acta1(H40Y)) (de Winter et al., 2021). Here, we show that acute tirasemtiv treatment results in a >50% increase in EDL muscle contractility in both wild type and Acta1(D286G) mice and a 45% and 64% increase in contractility in diaphragm muscle strips of Acta1(D286G) and wild type mice, respectively, at submaximal stimulation frequencies. Furthermore, we show that long-term tirasemtiv treatment via chow results in a significant increase in contractile force in EDL and GAS muscles at submaximal stimulations. This suggests that long-term administration of tirasemtiv does not desensitize fast twitch muscles. It should be noted, however, that we cannot rule out that the washout of tirasemtiv after long-term treatment was complete and did not confound our results.
In contrast to our previous work on an Acta1(H40Y) mouse model, here, we did not observe any significant improvements in EDL, GAS, or diaphragm muscle contractile function at maximal stimulation frequencies (de Winter et al., 2021). This may be explained by GAS and EDL muscles presenting with a more severe baseline phenotype in the Acta1(H40Y) model, with a significant reduction in maximal tension compared with the Acta1(D286G) model. Overall, acute and long-term tirasemtiv treatment had a significant positive inotropic effect on skeletal muscle contractile function at submaximal stimulation frequencies.
In addition to peripheral skeletal muscle weakness, NEM patients also present with weakness of the respiratory muscles, i.e., intercostal muscles and the diaphragm (van Kleef et al., 2022), which can cause insufficient respiration and in severe cases respiratory failure and death. Our results show that acute tirasemtiv treatment has a strong positive inotropic effect on diaphragm contractility at submaximal stimulation frequencies. This effect is in line with previous work performed on other myopathic mouse models (Lee et al., 2019; Hwee et al., 2014). Additionally, at a clinical level, tirasemtiv has been shown to improve spirometric parameters in patients diagnosed with amyotrophic lateral sclerosis and myasthenia gravis (Shefner et al., 2013, 2019; Sanders et al., 2015). Unlike the results shown in previous work (Hwee et al., 2014; de Winter et al., 2021), here we did not see any improvements in VT, VE, or RR after acute tirasemtiv treatment. However, long-term tirasemtiv treatment in Acta1(D286G) mice showed a significant decrease in RR without affecting their VT or VE, suggesting perhaps a more efficient respiratory cycle.
Next to respiratory function, other important outcome parameters to assess treatment effects are the energy cost of contraction and resistance to fatigue. We have previously shown that tirasemtiv treatment significantly lowers the energy cost and improves resistance to fatigability in skeletal muscles in a more severe NEM3 mouse model (de Winter et al., 2021). This positive effect has also been observed in healthy rats treated with CK-2066260, a structural analog of tirasemitiv. CK-2066260 treatment significantly increased muscle endurance and resistance to fatigue without an increase in ATP demand or glycogen usage (Cheng et al., 2019). This effect might have potential benefits for NEM patients who present with insufficient respiration and in severe cases, respiratory failure. Future studies should address the effect of tirasemtiv on other respiratory muscles. In the present study we focused on the diaphragm; however, accessory respiratory muscles may benefit from fast skeletal muscle troponin activators as well. Previous work characterizing fiber type composition in human postmortem muscle tissue and patient’s muscle biopsies showed that not only the diaphragm but also the accessory respiratory muscles (internal and external intercostal, scalene, and sternocleidomastoids muscles) and extra-thoracic muscles (rectus abdominis and transverse abdominis muscles) have a substantial percentage of fast-twitch muscle fibers on which tirasemtiv could potentially exert a positive ionotropic effect (Levine et al., 1997; Mizuno, 1991). Thus, tirasemtiv may result in a more efficient respiratory cycle and thereby alleviate respiratory muscle weakness in NEM patients.
Long-term tirasemtiv treatment reduces MHC carbonylation
A remarkable finding was the improved submaximal tension generation in EDL and GAS muscles of both wild type and Acta1(D286G) mice following long-term tirasemtiv treatment. Previous work on Acta1(D286G) mice has shown that the acto-myosin interaction is altered (Ochala et al., 2012), which can in part explain the reduction in contractility. However, other changes in muscle, such as protein carbonylation, have been associated with loss of muscle mass and muscle dysfunction in human chronic diseases and in animal models for different chronic conditions (Barreiro, 2016; Marin-Corral et al., 2010). Furthermore, altered redox balance may impair the acto–myosin interaction (Prochniewicz et al., 2008), and dysfunctional mitochondria can alter redox balance in mouse models of both moderate and mild forms of NEM3 (Tinklenberg et al., 2023), resulting in an increased oxidative stress and consequently, impaired muscle function, i.e., altered acto–myosin kinetics (Prochniewicz et al., 2008) and calcium sensitivity of force (Heunks et al., 2001). Interestingly, here, we show that tirasemtiv did not affect markers of mitochondrial biogenesis and antioxidant status. Furthermore, no changes in myofibrillar protein carbonylation following tirasemtiv treatment were found. We did however observe a significant reduction in MHC carbonylation in wild type and Acta1(D286G) mice after long-term tirasemtiv treatment. It appears that the effect of tirasemtiv on muscle contractility is independent of altered protein turnover and protein oxidation. Although speculative, the reduced MHC carbonylation may be a result of tirasemtiv’s effect on skeletal muscle energetics. Tirasemtiv reduces muscle energetic cost during activation (Cheng et al, 2019), which results in less free radicals, i.e., reactive oxygen and nitrogen species production. These radicals are known to impair muscle homeostasis and contractile function (Westerblad and Allen, 2011; Debold, 2015; Lian et al, 2022). Thus, we propose that tirasemtiv protects myofilament integrity by reducing free-radical production and consequently reducing MHC carbonylation. Future work should address the effect of tirasemtiv on other posttranslational modifications such as ubiquitination, acetylation, citrullination, S-glutathiolation, and S-nitrosylation, which are known to also regulate muscle function and metabolism (Fert-Bober et al, 2015; Collins and Goldberg, 2017; Jeong et al, 2018; Kramer et al, 2018; Liang et al, 2022).
Clinical relevance
NEMs are among the most common forms of congenital myopathies. In the past years, our knowledge regarding the molecular mechanisms driving various subtypes of NEM has significantly advanced (Ravenscroft et al., 2018; Laitila and Wallgren-Pettersson, 2021, Sewry et al., 2019). However, there is still an unmet need for the development of specific therapies. Our data provide a proof-of-concept supporting the positive role of fast skeletal muscle troponin activators in improving muscle function during daily life activities that require submaximal activation. To date, fast skeletal muscle troponin activators have not been tested in clinical trials for NEM patients. We propose that fast skeletal muscle troponin activators may aid NEM patients by increasing muscle tension at submaximal stimulations and support coughing and respiratory muscle function during episodes of insufficient respiration. A phase III clinical trial showed poor tolerability of tirasemtiv in patients with amyotrophic lateral sclerosis (Shefner et al., 2019). However, in this phase III trial, authors did note that participants who tolerated the designated dose showed a trend toward increased respiratory function. The development of tirasemtiv analogs may further contribute to the development of specific treatments for NEM.
Data availability
The data are available from the corresponding author upon reasonable request.
Acknowledgments
Henk L. Granzier served as editor.
This study was supported by the ERA-NET E-Rare2 grant TREAT-NEMMYOP to C.A.C. Ottenheijm, J. Gondin, and R. Bottinelli; The Netherlands Organisation for Health Research and Development (ZonMW) VICI Grant 91819613 to C.A.C. Ottenheijm; and ZonMW-VENI Grant 09150161910168 to J.M. de Winter.
Author contributions: Conceptualization: C.A.C. Ottenheijm, J. Gondin, R. Bottinelli, and J.M. de Winter. Data curation: C.A.C. Ottenheijm, J.M. de Winter, R.A. Galli, T.C. Borsboom, J. Gondin, C. Gineste, R. Bottinelli, M.-A. Pellegrino, and M. Rossi. Formal analysis: R.A. Galli, T.C. Borsboom, C. Gineste, M. Rossi, and L. Brocca. Funding acquisition: C.A.C. Ottenheijm, R. Bottinelli, J. Gondin, and J.M. de Winter. Investigation: C.A.C. Ottenheijm, J.M. de Winter, T.C. Borsboom, C. Gineste, L. Brocca, M. Rossi, D.T. Hwee, F.I. Malik, R. Bottinelli, J. Gondin, and M.-A. Pellegrino. Methodology: C.A.C. Ottenheijm, J.M. de Winter, R. Bottinelli, M.-A. Pellegrino, and J. Gondin. Project administration: C.A.C. Ottenheijm and J.M. de Winter. Resources: C.A.C. Ottenheijm, J. Gondin, R. Bottinelli, J.M. de Winter, D.T. Hwee, and F.I. Malik. Supervision: C.A.C. Ottenheijm, J.M. de Winter, R. Bottinelli, J. Gondin, M.-A. Pellegrino, and L. Brocca. Visualization: R.A. Galli, C.A.C. Ottenheijm, and J.M. de Winter. Writing—original draft: R.A. Galli, C.A.C. Ottenheijm, and J.M. de Winter. Writing—review & editing: R.A. Galli, C.A.C. Ottenheijm, J.M. de Winter, J. Gondin, R. Bottinelli, C. Gineste, M.-A. Pellegrino, D.T. Hwee, and F.I. Malik.
References
Author notes
Disclosures: D.T. Hwee reported other from Cytokinetics during the conduct of the study and other from Cytokinetics outside the submitted work. In addition, D.T. Hwee had a patent to tirasemtiv use issued for Cytokinetics. F.I. Malik reported personal fees from Cytokinetics, Inc. during the conduct of the study and personal fees from Cytokinetics, Inc. outside the submitted work. F.I. Malik also had a patent to tirasemtiv patent licensed to Cytokinetics. No other disclosures were reported.