Quantification of skeletal muscle mitochondrial function by 31P magnetic resonance spectroscopy techniques: a quantitative review.

Magnetic resonance spectroscopy (MRS) can give information about cellular metabolism in vivo which is difficult to obtain in other ways. In skeletal muscle, non-invasive (31) P MRS measurements of the post-exercise recovery kinetics of pH, [PCr], [Pi] and [ADP] contain valuable information about muscle mitochondrial function and cellular pH homeostasis in vivo, but quantitative interpretation depends on understanding the underlying physiology. Here, by giving examples of the analysis of (31) P MRS recovery data, by some simple computational simulation, and by extensively comparing data from published studies using both (31) P MRS and invasive direct measurements of muscle O2 consumption in a common analytical framework, we consider what can be learnt quantitatively about mitochondrial metabolism in skeletal muscle using MRS-based methodology. We explore some technical and conceptual limitations of current methods, and point out some aspects of the physiology which are still incompletely understood.

Pioglitazone treatment restores in vivo muscle oxidative capacity in a rat model of diabetes.

AIM: To determine the effect of pioglitazone treatment on in vivo and ex vivo muscle mitochondrial function in a rat model of diabetes. METHODS: Both the lean, healthy rats and the obese, diabetic rats are Zucker Diabetic Fatty (ZDF) rats. The homozygous fa/fa ZDF rats are obese and diabetic. The heterozygous fa/+ ZDF rats are lean and healthy. Diabetic Zucker Diabetic Fatty rats were treated with either pioglitazone (30 mg/kg/day) or water as a control (n = 6 per group), for 2 weeks. In vivo ¹H and ³¹P magnetic resonance spectroscopy was performed on skeletal muscle to assess intramyocellular lipid (IMCL) content and muscle oxidative capacity, respectively. Ex vivo muscle mitochondrial respiratory capacity was evaluated using high-resolution respirometry. In addition, several markers of mitochondrial content were determined. RESULTS: IMCL content was 14-fold higher and in vivo muscle oxidative capacity was 26% lower in diabetic rats compared with lean rats, which was, however, not caused by impairments of ex vivo mitochondrial respiratory capacity or a lower mitochondrial content. Pioglitazone treatment restored in vivo muscle oxidative capacity in diabetic rats to the level of lean controls. This amelioration was not accompanied by an increase in mitochondrial content or ex vivo mitochondrial respiratory capacity, but rather was paralleled by an improvement in lipid homeostasis, that is lowering of plasma triglycerides and muscle lipid and long-chain acylcarnitine content. CONCLUSION: Diminished in vivo muscle oxidative capacity in diabetic rats results from mitochondrial lipid overload and can be alleviated by redirecting the lipids from the muscle into adipose tissue using pioglitazone treatment.

Combined in vivo and in silico investigations of activation of glycolysis in contracting skeletal muscle.

The hypothesis was tested that the variation of in vivo glycolytic flux with contraction frequency in skeletal muscle can be qualitatively and quantitatively explained by calcium-calmodulin activation of phosphofructokinase (PFK-1). Ischemic rat tibialis anterior muscle was electrically stimulated at frequencies between 0 and 80 Hz to covary the ATP turnover rate and calcium concentration in the tissue. Estimates of in vivo glycolytic rates and cellular free energetic states were derived from dynamic changes in intramuscular pH and phosphocreatine content, respectively, determined by phosphorus magnetic resonance spectroscopy ((31)P-MRS). Computational modeling was applied to relate these empirical observations to understanding of the biochemistry of muscle glycolysis. Hereto, the kinetic model of PFK activity in a previously reported mathematical model of the glycolytic pathway (Vinnakota KC, Rusk J, Palmer L, Shankland E, Kushmerick MJ. J Physiol 588: 1961-1983, 2010) was adapted to contain a calcium-calmodulin binding sensitivity. The two main results were introduction of regulation of PFK-1 activity by binding of a calcium-calmodulin complex in combination with activation by increased concentrations of AMP and ADP was essential to qualitatively and quantitatively explain the experimental observations. Secondly, the model predicted that shutdown of glycolytic ATP production flux in muscle postexercise may lag behind deactivation of PFK-1 (timescales: 5-10 s vs. 100-200 ms, respectively) as a result of accumulation of glycolytic intermediates downstream of PFK during contractions.

In vivo postprandial lipid partitioning in liver and skeletal muscle in prediabetic and diabetic rats.

AIMS/HYPOTHESIS: Insulin resistance and type 2 diabetes have been associated with ectopic lipid deposition. This study investigates the derangements in postprandial lipid handling in liver and skeletal muscle tissue at different stages during the pathogenesis of type 2 diabetes in a rat model. METHODS: Four groups (n = 6) of male Zucker diabetic fatty rats were used for this study: prediabetic fa/fa rats and healthy fa/+ littermates at the age of 6 weeks, and diabetic fa/fa rats and healthy fa/+ littermates at the age of 12 weeks. In vivo (1)H-[(13)C] magnetic resonance spectroscopy measurements were performed in liver and tibialis anterior muscle at baseline and 4, 24 and 48 h after oral administration of 1.5 g [U-(13)C]algal lipid mixture per kilogram body weight. Total and (13)C-labelled intracellular lipid contents were determined from the magnetic resonance spectra. RESULTS: In both prediabetic and diabetic rats, total lipid contents in muscle and liver were substantially higher than in healthy controls and this was accompanied by a 2.3-fold greater postprandial lipid uptake in the liver (p < 0.001). Interestingly, in prediabetic rats, skeletal muscle appeared to be protected from excess lipid uptake whereas after developing overt diabetes muscle lipid uptake was 3.4-fold higher than in controls (p < 0.05). Muscle lipid use was significantly lower in prediabetic and diabetic muscle, indicative of impairments in lipid oxidation. CONCLUSIONS/INTERPRETATION: In vivo postprandial lipid handling is disturbed in both liver and skeletal muscle tissue in prediabetic and diabetic rats, but the uptake of dietary lipids in muscle is only increased after the development of overt diabetes.

Physical activity is the key determinant of skeletal muscle mitochondrial function in type 2 diabetes.

CONTEXT: Conflicting data exist on mitochondrial function and physical activity in type 2 diabetes mellitus (T2DM) development. OBJECTIVE: The aim was to assess mitochondrial function at different stages during T2DM development in combination with physical exercise in longstanding T2DM patients. DESIGN AND METHODS: We performed cross-sectional analysis of skeletal muscle from 12 prediabetic 11 longstanding T2DM male subjects and 12 male controls matched by age and body mass index. INTERVENTION: One-year intrasubject controlled supervised exercise training intervention was done in longstanding T2DM patients. MAIN OUTCOME MEASUREMENTS: Extensive ex vivo analyses of mitochondrial quality, quantity, and function were collected and combined with global gene expression analysis and in vivo ATP production capacity after 1 yr of training. RESULTS: Mitochondrial density, complex I activity, and the expression of Krebs cycle and oxidative phosphorylation system-related genes were lower in longstanding T2DM subjects but not in prediabetic subjects compared with controls. This indicated a reduced capacity to generate ATP in longstanding T2DM patients only. Gene expression analysis in prediabetic subjects suggested a switch from carbohydrate toward lipid as an energy source. One year of exercise training raised in vivo skeletal muscle ATP production capacity by 21 ± 2% with an increased trend in mitochondrial density and complex I activity. In addition, expression levels of ß-oxidation, Krebs cycle, and oxidative phosphorylation system-related genes were higher after exercise training. CONCLUSIONS: Mitochondrial dysfunction is apparent only in inactive longstanding T2DM patients, which suggests that mitochondrial function and insulin resistance do not depend on each other. Prolonged exercise training can, at least partly, reverse the mitochondrial impairments associated with the longstanding diabetic state.

Comparison of in vivo postexercise phosphocreatine recovery and resting ATP synthesis flux for the assessment of skeletal muscle mitochondrial function.

(31)P magnetic resonance spectroscopy (MRS) has been used to assess skeletal muscle mitochondrial function in vivo by measuring 1) phosphocreatine (PCr) recovery after exercise or 2) resting ATP synthesis flux with saturation transfer (ST). In this study, we compared both parameters in a rat model of mitochondrial dysfunction with the aim of establishing the most appropriate method for the assessment of in vivo muscle mitochondrial function. Mitochondrial dysfunction was induced in adult Wistar rats by daily subcutaneous injections with the complex I inhibitor diphenyleneiodonium (DPI) for 2 wk. In vivo (31)P MRS measurements were supplemented by in vitro measurements of oxygen consumption in isolated mitochondria. Two weeks of DPI treatment induced mitochondrial dysfunction, as evidenced by a 20% lower maximal ADP-stimulated oxygen consumption rate in isolated mitochondria from DPI-treated rats oxidizing pyruvate plus malate. This was paralleled by a 46% decrease in in vivo oxidative capacity, determined from postexercise PCr recovery. Interestingly, no significant difference in resting, ST-based ATP synthesis flux was observed between DPI-treated rats and controls. These results show that PCr recovery after exercise has a more direct relationship with skeletal muscle mitochondrial function than the ATP synthesis flux measured with (31)P ST MRS in the resting state.

Increased mitochondrial content rescues in vivo muscle oxidative capacity in long-term high-fat-diet-fed rats.

Mitochondria are thought to play a crucial role in the etiology of muscle insulin resistance (IR). The aim of this study was to gain more insight into the timing and nature of mitochondrial adaptations during the development of high-fat-diet (HFD)-induced IR. Adult Wistar rats were fed HFD or normal chow for 2.5 and 25 wk. Intramyocellular lipids (IMCLs) were quantified in vivo using (1)H magnetic resonance spectroscopy (MRS). Muscle oxidative capacity was assessed in vivo using (31)P MRS and in vitro by measuring mitochondrial DNA copy number and oxygen consumption in isolated mitochondria. MRS in tibialis anterior muscle revealed 3.3-fold higher IMCL content and 1.2-fold increased oxidative capacity after 2.5 wk of HFD feeding. The latter result could be fully accounted for by increased mitochondrial content. After 25 wk of HFD, maximal ADP-stimulated oxygen consumption in isolated mitochondria oxidizing pyruvate plus malate remained unaffected, while IMCL and mitochondrial content had further increased compared to controls (5.1-fold and 1.4-fold, respectively). Interestingly, in vivo oxidative capacity at this time point was identical to controls. These results show that skeletal muscle in HFD-induced IR accompanied by IMCL accumulation requires a progressively larger mitochondrial pool size to maintain normal oxidative capacity in vivo.

Early or advanced stage type 2 diabetes is not accompanied by in vivo skeletal muscle mitochondrial dysfunction.

OBJECTIVE: Several lines of evidence support a potential role of skeletal muscle mitochondrial dysfunction in the pathogenesis of insulin resistance and/or type 2 diabetes. However, it remains to be established whether mitochondrial dysfunction represents either cause or consequence of the disease. We examined in vivo skeletal muscle mitochondrial function in early and advanced stages of type 2 diabetes, with the aim to gain insight in the proposed role of mitochondrial dysfunction in the aetiology of insulin resistance and/or type 2 diabetes. METHODS: Ten long-standing, insulin-treated type 2 diabetes patients, 11 subjects with impaired fasting glucose, impaired glucose tolerance and/or recently diagnosed type 2 diabetes, and 12 healthy, normoglycaemic controls, matched for age and body composition and with low habitual physical activity levels were studied. In vivo mitochondrial function of the vastus lateralis muscle was evaluated from post-exercise phosphocreatine (PCr) recovery kinetics using (31)P magnetic resonance spectroscopy (MRS). Intramyocellular lipid (IMCL) content was assessed in the same muscle using single-voxel (1)H MRS. RESULTS: IMCL content tended to be higher in the type 2 diabetes patients when compared with normoglycaemic controls (P=0.06). The(31)P MRS parameters for mitochondrial function, i.e. PCr and ADP recovery time constants and maximum aerobic capacity, did not differ between groups. CONCLUSIONS: The finding that in vivo skeletal muscle oxidative capacity does not differ between long-standing, insulin-treated type 2 diabetes patients, subjects with early stage type 2 diabetes and sedentary, normoglycaemic controls suggests that mitochondrial dysfunction does not necessarily represent either cause or consequence of insulin resistance and/or type 2 diabetes.

31P MR spectroscopy and in vitro markers of oxidative capacity in type 2 diabetes patients.

BACKGROUND: Skeletal muscle mitochondrial function in type 2 diabetes (T2D) is currently being studied intensively. In vivo (31)P magnetic resonance spectroscopy ((31)P MRS) is a noninvasive tool used to measure mitochondrial respiratory function (MIFU) in skeletal muscle tissue. However, microvascular co-morbidity in long-standing T2D can interfere with the (31)P MRS methodology. AIM: To compare (31)P MRS-derived parameters describing in vivo MIFU with an in vitro assessment of muscle respiratory capacity and muscle fiber-type composition in T2D patients. METHODS: (31)P MRS was applied in long-standing, insulin-treated T2D patients. (31)P MRS markers of MIFU were measured in the M. vastus lateralis. Muscle biopsy samples were collected from the same muscle and analyzed for succinate dehydrogenase activity (SDH) and fiber-type distribution. RESULTS: Several (31)P MRS parameters of MIFU showed moderate to good correlations with the percentage of type I fibers and type I fiber-specific SDH activity (Pearson's R between 0.70 and 0.75). In vivo and in vitro parameters of local mitochondrial respiration also correlated well with whole-body fitness levels (VO (2peak)) in these patients (Pearson's R between 0.62 and 0.90). CONCLUSION: Good correlations exist between in vivo and in vitro measurements of MIFU in long-standing insulin-treated T2D subjects, which are qualitatively and quantitatively consistent with previous results measured in healthy subjects. This justifies the use of (31)P MRS to measure MIFU in relation to T2D.

Reorientational eigenmode dynamics: a combined MD/NMR relaxation analysis method for flexible parts in globular proteins.

An approach is presented for the interpretation of heteronuclear NMR spin relaxation data in mobile protein parts in terms of reorientational eigenmode dynamics. The method is based on the covariance matrix of the spatial functions of the nuclear spin interactions that cause relaxation expressed as spherical harmonics of rank 2. The approach was applied to characterize the dynamics of a loop region of ubiquitin. The covariance matrix was determined from a conformational ensemble generated by a 5 ns molecular dynamics simulation. It was found that the time correlation functions of the dominant eigenmodes decay in good approximation with a single correlation time. From the reorientational eigenmodes, their eigenvalues, and correlation times, NMR relaxation data were calculated in accordance with Bloch-Wangsness-Redfield relaxation theory and directly compared with experimental (15)N relaxation parameters. Using a fitting procedure, agreement between calculated and experimental data was improved significantly by adjusting eigenvalues and correlation times of the dominant modes. The presented procedure provides detailed information on correlated reorientational dynamics of flexible parts in globular proteins. The covariance matrix was linked to the covariance matrix of backbone dihedral angle fluctuations, allowing one to study the motional behavior of these degrees of freedom on nano- and subnanosecond time scales.

Model-free approach to the dynamic interpretation of residual dipolar couplings in globular proteins.

The effects of internal motions on residual dipolar NMR couplings of proteins partially aligned in a liquid-crystalline environment are analyzed using a 10 ns molecular dynamics (MD) computer simulation of ubiquitin. For a set of alignment tensors with different orientations and rhombicities, MD-averaged dipolar couplings are determined and subsequently interpreted for different scenarios in terms of effective alignment tensors, average orientations of dipolar vectors, and intramolecular reorientational vector distributions. Analytical relationships are derived that reflect similarities and differences between motional scaling of dipolar couplings and scaling of dipolar relaxation data (NMR order parameters). Application of the self-consistent procedure presented here to dipolar coupling measurements of biomolecules aligned in different liquid-crystalline media should allow one to extract in a "model-free" way average orientations of dipolar vectors and specific aspects of their motions.

Collective reorientational motion and nuclear spin relaxation in proteins.

Significant progress in NMR methodology for measuring spin-relaxation data at many different 15N and 13C sites in proteins demands new and increasingly sophisticated ways of data interpretation. Recent work of our group concerning the use of anisotropic and reorientational collective motional models for spin-relaxation interpretation is briefly reviewed and a number of important aspects of collective reorientational motional models are discussed at the example of a 11 ns molecular dynamics computer simulation of the protein ubiquitin.

Characterization of NMR relaxation-active motions of a partially folded A-state analogue of ubiquitin.

The dominant dynamics of a partially folded A-state analogue of ubiquitin that give rise to NMR 15N spin relaxation have been investigated using molecular dynamics (MD) computer simulations and reorientational quasiharmonic analysis. Starting from the X-ray structure of native ubiquitin with a protonation state corresponding to a low pH, the A-state analogue was generated by a MD simulation of a total length of 33 ns in a 60%/40% methanol/water mixture using a variable temperature scheme to control and speed up the structural transformation. The N-terminal half of the A-state analogue consists of loosely coupled native-like secondary structural elements, while the C-terminal half is mostly irregular in structure. Analysis of dipolar N-H backbone correlation functions reveals reorientational amplitudes and time-scale distributions that are comparable to those observed experimentally. Thus, the trajectory provides a realistic picture of a partially folded protein that can be used for gaining a better understanding of the various types of reorientational motions that are manifested in spin-relaxation parameters of partially folded systems. For this purpose, a reorientational quasiharmonic reorientational analysis was performed on the final 5 ns of the trajectory of the A-state analogue, and for comparison on a 5 ns trajectory of native ubiquitin. The largest amplitude reorientational modes show a markedly distinct behavior for the two states. While for native ubiquitin, such motions have a more local character involving loops and the C-terminal end of the polypeptide chain, the A-state analogue shows highly collective motions in the nanosecond time-scale range corresponding to larger-scale movements between different segments. Changes in reorientational backbone entropy between the A-state analogue and the native state of ubiquitin, which were computed from the reorientational quasiharmonic analyses, are found to depend significantly on motional correlation effects.

Thermal unfolding of a llama antibody fragment: a two-state reversible process.

Camelids produce functional "heavy chain" antibodies which are devoid of light chains and CH1 domains [Hamers-Casterman, C., et al. (1993) Nature 363, 446-448]. It has been shown that the variable domains of these heavy chain antibodies (the V(HH) fragments) are functional at or after exposure to high temperatures, in contrast to conventional antibodies [Linden van der, R. H. J., et al. (1999) Biochim. Biophys. Acta 1431, 37-44]. For a detailed understanding of the higher thermostability of these V(HH) fragments, knowledge of their structure and conformational dynamics is required. As a first step toward this goal, we report here the essentially complete (1)H and (15)N NMR backbone resonance assignments of a llama V(HH) antibody fragment, and an extensive analysis of the structure at higher temperatures. The H-D exchange NMR data at 300 K indicate that the framework of the llama V(HH) fragment is highly protected with a DeltaG(ex) of >5.4 kcal/mol, while more flexibility is observed for surface residues, particularly in the loops and the two outer strands (residues 4-7, 10-13, and 58-60) of the beta-sheet. The CD data indicate a reversible, two-state unfolding mechanism with a melting transition at 333 K and a DeltaH(m) of 56 kcal/mol. H-D exchange studies using NMR and ESI-MS show that below 313 K exchange occurs through local unfolding events whereas above 333 K exchange mainly occurs through global unfolding. The lack of a stable core at high temperatures, observed for V(HH) fragments, has also been observed for conventional antibody fragments. The main distinction between the llama V(HH) fragment and conventional antibody fragments is the reversibility of the thermal unfolding process, explaining its retained functionality after exposure to high temperatures.

Tryptophan mediated photoreduction of disulfide bond causes unusual fluorescence behaviour of Fusarium solani pisi cutinase.

The fluorescence signal of the single tryptophan residue (Trp69) of Fusarium solani pisi cutinase is highly quenched. However, prolonged irradiation of the enzyme in the tryptophan absorption band causes an increase of the tryptophan fluorescence quantum yield by an order of magnitude. By using a combination of NMR spectroscopy and chemical detection of free thiol groups with a sulfhydryl reagent we could unambiguously show that the unusual fluorescence behaviour of Trp69 in cutinase is caused by the breaking of the disulfide bond between Cys31 and Cys109 upon irradiation, while the amide-aromatic hydrogen bond between Ala32 and Trp69 remains intact. This is the first example of tryptophan mediated photoreduction of a disulfide bond in proteins.

NMR studies of Fusarium solani pisi cutinase in complex with phosphonate inhibitors.

The backbone dynamics of Fusarium solani pisi cutinase in complex with a phosphonate inhibitor has been studied by a variety of nuclear magnetic resonance experiments to probe internal motions on different time scales. The results have been compared with dynamical studies performed on free cutinase. In solution, the enzyme adopts its active conformation only upon binding the inhibitor. While the active site Ser120 is rigidly attached to the stable alpha/beta core of the protein, the remainder of the binding site is very flexible in the free enzyme. The other two active site residues Asp175 and His188 as well as the oxyanion hole residues Ser42 and Gln121 are only restrained into their proper positions upon binding of the substrate-like inhibitor. The flap helix, which opens and closes the binding site in the free molecule, is also fixed in the cutinase-inhibitor complex. Our results are in contrast with the X-ray analysis results, namely that in the protein crystal, free cutinase has a well-defined active site and a preformed oxyanion hole and that it does not need any rearrangements to bind its substrate. Our solution studies show that cutinase does need conformational rearrangements to bind its substrate, which may form the rate-limiting step in catalysis.

Backbone dynamics of Fusarium solani pisi cutinase probed by nuclear magnetic resonance: the lack of interfacial activation revisited.

The backbone dynamics of Fusarium solani pisi cutinase has been studied by a variety of nuclear magnetic resonance experiments to probe internal motions on different time scales. The core of cutinase appears to be highly rigid. The binding site, including the oxyanion hole, is mobile on the microsecond to millisecond time scale, in contrast to the well-defined active site and preformed oxyanion hole elucidated by X-ray crystallography [Martinez, C., de Geus, P., Lauwereys, M., Matthyssens, G., and Cambillau, C. (1992) Nature 356, 615-618]. In this crystal structure, cutinase has a rather open conformation, in which the hydrophobic binding site is exposed. The observed mobility in solution most likely represents the interconversion between open and more closed conformations, like in a true lipase. The opening and closing motions are on a time scale which corresponds with the kinetics of the hydrolysis reaction, i.e., the millisecond range, which suggests that these conformational rearrangements form the rate-limiting step in catalysis. We conclude that the crystal structure probably represents one of the multiple conformations present in solution, which fortuitously is the active conformation. The implications of our findings are discussed with particular reference to the explanation of the lack of interfacial activation as found for cutinase.

1H, 13C, and 15N resonance assignments of Fusarium solani pisi cutinase and preliminary features of the structure in solution.

Essentially complete (96%) sequence-specific assignments were made for the backbone and side-chain 1H, 13C, and 15N resonances of Fusarium solani pisi cutinase, produced as a 214-residue heterologous protein in Escherichia coli, using heteronuclear NMR techniques. Three structural features were noticed during the assignment. (1) The secondary structure in solution corresponds mostly with the structure from X-ray diffraction, suggesting that both structures are globally similar. (2) The HN of Ala32 has a strongly upfield-shifted resonance at 3.97 ppm, indicative of an amide-aromatic hydrogen bond to the indole ring of Trp69 that stabilizes the N-terminal side of the parallel beta-sheet. (3) The NMR data suggest that the residues constituting the oxyanion hole are quite mobile in the free enzyme in solution, in contrast to the existence of a preformed oxyanion hole as observed in the crystal structure. Apparently, cutinase forms its oxyanion hole upon binding of the substrate like true lipases.

Refined solution structure of the Tyr41-->His mutant of the M13 gene V protein. A comparison with the crystal structure.

The three-dimensional solution structure of mutant Tyr41-->His of the single-stranded DNA binding protein encoded by gene V of the filamentous bacteriophage M13 has been refined in two stages. The first stage involved the collection of additional NOE-based distance constraints, which were then used in eight cycles of back-calculations and structure calculations. The structures of the gene V protein dimers were calculated using simulated annealing, employing restrained molecular dynamics with a geometric force field. In the second stage of the refinement procedure, solvent was explicitly included during the dynamic calculations. A total of 30 structures was calculated for the protein, representing its solution structure in water. The first calculation step significantly improved the convergence of the structures, whereas the subsequent simulations in water made the structures physically more realistic. This is, for instance, illustrated by the number of hydrogen bonds formed in the molecule, which increased considerably upon going to aqueous solution. It is shown that the solution structure of the mutant gene V protein is nearly identical to the crystal structure of the wild-type molecule, except for the DNA-binding loop (residues 16-28). This antiparallel beta-hairpin is twisted and partially folded back towards the core of the protein in the NMR structure, whereas it is more extended and points away from the rest of the molecule in the X-ray structure. Unrestrained molecular dynamics calculations suggest that this latter conformation is energetically unstable in solution.