The Dynamics of Locomotor Neuromuscular Fatigue during Ramp-Incremental Cycling to Intolerance.

INTRODUCTION: Traditional neuromuscular fatigue assessments are not task-specific and are unable to characterize neuromuscular performance decline during dynamic whole-body exercise. This study used interleaved maximal isokinetic cycling efforts to characterize the dynamics of the decline in neuromuscular performance during ramp-incremental (RI) cycle ergometry exercise to intolerance. METHODS: Eleven young healthy participants (10 male/1 female) performed two RI cycle ergometry exercise tests to intolerance: [1] RI-exercise with peak isokinetic power (Piso) at 80 rev·min-1 measured at baseline and immediately at intolerance from a maximal ~6 s effort; [2] RI-exercise where additional Piso measurements were interleaved every 90 s to characterize the decline in neuromuscular performance during the RI-test. Muscle excitation was measured using EMG during all Piso assessments, and pulmonary gas exchange was measured throughout. RESULTS: Baseline Piso was 832 ± 140 W and RI-exercise reduced Piso to 349 ± 96 W at intolerance (p = 0.001), which was not different from flywheel power at intolerance (303 ± 96 W; p = 0.292). There was no reduction in Piso between baseline cycling and gas exchange threshold (GET; baseline Piso vs. mean Piso below GET: 828 ± 146 vs. 815 ± 149 W; p = 1.00). Piso fell progressively above GET until intolerance (Piso every 90 s above GET: 759 ± 139; 684 ± 141; 535 ± 144; 374 ± 117 W; each p < 0.05 vs. baseline and mean Piso below GET). Peak muscle excitation (EMG) was also reduced only above GET (73 ± 14 % of baseline, at intolerance; p < 0.05). However, the reduction in peak Piso preceded the reduction in peak muscle excitation. CONCLUSIONS: The dynamics of the decline in neuromuscular performance (reduction in Piso and EMG) during RI-exercise are consistent with known intensity-dependent metabolic and traditional pre-post neuromuscular fatigue responses to discrete bouts of constant-power exercise.

A simple approach for image-based modelling of the heart that enables robust simulation of highly heterogeneous electrical excitation.

Remodelling of cardiac tissue structure, including intercellular electrical coupling, is a major determinant of the complex and heterogeneous excitation patterns associated with cardiac arrhythmias. Evaluation of the precise mechanisms by which local tissue structure determines global arrhythmic excitation patterns is a major challenge that may be critically important for the development of effective treatment strategies. Computational modelling is a key tool in the study of cardiac arrhythmias, yet the established approaches for organ-scale modelling are unsuitable to capture the impact of local conduction heterogeneities; a novel approach is required to provide this multi-scale mechanistic insight. We present a fundamentally simple yet powerful approach to simulate electrical excitation in highly heterogeneous whole-heart models that exploits the underlying discreteness of the myocardium. Preliminary simulations demonstrate that this approach can capture lower conduction velocities and reproduce wave breakdown and the development of re-entry in a range of conditions.

Ionic current changes underlying action potential repolarization responses to physiological pacing and adrenergic stimulation in adult rat ventricular myocytes.

This study aimed to simulate ventricular responses to elevations in myocyte pacing and adrenergic stimulation using a novel electrophysiological rat model and investigate ion channel responses underlying action potential (AP) modulations. Peak ion currents and AP repolarization to 50% and 90% of full repolarization (APD50-90 ) were recorded during simulations at 1-10 Hz pacing under control and adrenergic stimulation conditions. Further simulations were performed with incremental ion current block (L-type calcium current, ICa ; transient outward current, Ito ; slow delayed rectifier potassium current, IKs ; rapid delayed rectifier potassium current, IKr ; inward rectifier potassium current, IK1 ) to identify current influence on AP response to exercise. Simulated APD50-90 closely resembled experimental findings. Rate-dependent increases in IKs (6%-101%), IKr (141%-1339%), and ICa (0%-15%) and reductions in Ito (11%-57%) and IK1 (1%-9%) were observed. Meanwhile, adrenergic stimulation triggered moderate increases in all currents (23%-67%) except IK1 . Further analyses suggest AP plateau is most sensitive to modulations in Ito and ICa while late repolarization is most sensitive to IK1 , ICa , and IKs , with alterations in IKs predominantly stimulating the greatest magnitude of influence on late repolarization (35%-846% APD90 prolongation). The modified Leeds rat model (mLR) is capable of accurately modeling APs during physiological stress. This study highlights the importance of ICa , Ito , IK1, and IKs in controlling electrophysiological responses to exercise. This work will benefit the study of cardiac dysfunction, arrythmia, and disease, though future physiologically relevant experimental studies and model development are required.

Power Reserve at Intolerance in Ramp-Incremental Exercise Is Dependent on Incrementation Rate.

INTRODUCTION: The mechanism(s) of exercise intolerance at V˙O2max remain poorly understood. In health, standard ramp-incremental (RI) exercise is limited by fatigue-induced reductions in maximum voluntary cycling power. Whether neuromuscular fatigue also limits exercise when the RI rate is slow and RI peak power at intolerance is lower than standard RI exercise, is unknown. METHODS: In twelve healthy participants, maximal voluntary cycling power was measured during a short (~6 s) isokinetic effort at 80 rpm (Piso) at baseline and, using an instantaneous switch from cadence-independent to isokinetic cycling, immediately at the limit of RI exercise with RI rates of 50, 25, and 10 W·min-1 (RI-50, RI-25, and RI-10). Breath-by-breath pulmonary gas exchange was measured throughout. RESULTS: Baseline Piso was not different among RI rates (analysis of variance; P > 0.05). Tolerable duration increased with decreasing RI rate (RI-50, 411 ± 58 s vs RI-25, 732 ± 93 s vs RI-10, 1531 ± 288 s; P < 0.05). At intolerance, V˙O2peak was not different among RI rates (analysis of variance; P > 0.05), but RI peak power decreased with RI rate (RI-50, 361 ± 48 W vs RI-25, 323 ± 39 W vs RI-10, 275 ± 38 W; P < 0.05). Piso at intolerance was 346 ± 43 W, 353 ± 45 W, and 392 ± 69 W for RI-50, RI-25, and RI-10, respectively (P < 0.05 for RI-10 vs RI-50 and RI-25). At intolerance, in RI-50 and RI-25, Piso was not different from RI peak power (P > 0.05), thus there was no "power reserve." In RI-10, Piso was greater than RI peak power at intolerance (P < 0.001), that is, there was a "power reserve." CONCLUSIONS: In RI-50 and RI-25, the absence of a power reserve suggests the neuromuscular fatigue-induced reduction in Piso coincided with V˙O2max and limited the exercise. In RI-10, the power reserve suggests neuromuscular fatigue was insufficient to limit the exercise, and additional mechanisms contributed to intolerance at V˙O2max.

Deterministic and Stochastic Cellular Mechanisms Contributing to Carbon Monoxide Induced Ventricular Arrhythmias.

Chronic exposure to low levels of Carbon Monoxide is associated with an increased risk of cardiac arrhythmia. Microelectrode recordings from rat and guinea pig single isolated ventricular myocytes exposed to CO releasing molecule CORM-2 and excited at 0.2/s show repolarisation changes that develop over hundreds of seconds: action potential prolongation by delayed repolarisation, EADs, multiple EADs and oscillations around the plateau, leading to irreversible repolarisation failure. The measured direct effects of CO on currents in these cells, and ion channels expressed in mammalian systems showed an increase in prolonged late Na+, and a decrease in the maximal T- and L-type Ca++. peak and late Na+, ultra-rapid delayed, delayed rectifier, and the inward rectifier K+ currents. Incorporation of these CO induced changes in maximal currents in ventricular cell models; (Gattoni et al., J. Physiol., 2016, 594, 4193-4224) (rat) and (Luo and Rudy, Circ. Res., 1994, 74, 1071-1096) (guinea-pig) and human endo-, mid-myo- and epi-cardial (O'Hara et al., PLoS Comput. Biol., 2011, 7, e1002061) models, by changes in maximal ionic conductance reproduces these repolarisation abnormalities. Simulations of cell populations with Gaussian distributions of maximal conductance parameters predict a CO induced increase in APD and its variability. Incorporation of these predicted CO induced conductance changes in human ventricular cell electrophysiology into ventricular tissue and wall models give changes in indices for the probability of the initiation of re-entrant arrhythmia.

Multi-scale approaches for the simulation of cardiac electrophysiology: II - Tissue-level structure and function.

Computational models of the heart, from cell-level models, through one-, two- and three-dimensional tissue-level simplifications, to biophysically-detailed three-dimensional models of the ventricles, atria or whole heart, allow the simulation of excitation and propagation of this excitation, and have provided remarkable insight into the normal and pathological functioning of the heart. In this article we present equations for modelling cellular excitation (i.e. the cell action potential) from both a phenomenological and a biophysical perspective. Hodgkin-Huxley formalism is discussed, along with the current generation of biophysically-detailed cardiac cell models. Alternative Markovian formulations for modelling ionic currents are also presented. Equations describing propagation of this cellular excitation, through one-, two- and three-dimensional idealised or realistic tissues, are then presented. For all types of model, from cell to tissue, methods for discretisation and integration of the underlying equations are discussed. The article finishes with a discussion of two tissue-level experimental imaging techniques - diffusion tensor magnetic resonance imaging and optical imaging - that can be used to provide data for parameterisation and validation of cell- and tissue-level cardiac models.

Multi-scale approaches for the simulation of cardiac electrophysiology: I - Sub-cellular and stochastic calcium dynamics from cell to organ.

Computational models of the heart at multiple spatial scales, from sub-cellular nanodomains to the whole-organ, are a powerful tool for the simulation of cardiac electrophysiology. Application of these models has provided remarkable insight into the normal and pathological functioning of the heart. In these two articles, we present methods for modelling cardiac electrophysiology at all of these spatial scales. In part one, presented here, we discuss methods and approaches for modelling sub-cellular calcium dynamics at the whole-cell and organ scales, valuable for modelling excitation-contraction coupling and mechanisms of arrhythmia triggers.

Abnormal skeletal muscle blood flow, contractile mechanics and fibre morphology in a rat model of obese-HFpEF.

KEY POINTS: Heart failure is characterised by limb and respiratory muscle impairments that limit functional capacity and quality of life. However, compared with heart failure with reduced ejection fraction (HFrEF), skeletal muscle alterations induced by heart failure with preserved ejection fraction (HFpEF) remain poorly explored. Here we report that obese-HFpEF induces multiple skeletal muscle alterations in the rat hindlimb, including impaired muscle mechanics related to shortening velocity, fibre atrophy, capillary loss, and an impaired blood flow response to contractions that implies a perfusive oxygen delivery limitation. We also demonstrate that obese-HFpEF is characterised by diaphragmatic alterations similar to those caused by denervation - atrophy in Type IIb/IIx (fast/glycolytic) fibres and hypertrophy in Type I (slow/oxidative) fibres. These findings extend current knowledge in HFpEF skeletal muscle physiology, potentially underlying exercise intolerance, which may facilitate future therapeutic approaches. ABSTRACT: Peripheral skeletal muscle and vascular alterations induced by heart failure with preserved ejection fraction (HFpEF) remain poorly identified, with limited therapeutic targets. This study used a cardiometabolic obese-HFpEF rat model to comprehensively phenotype skeletal muscle mechanics, blood flow, microvasculature and fibre atrophy. Lean (n = 8) and obese-HFpEF (n = 8) ZSF1 rats were compared. Skeletal muscles (soleus and diaphragm) were assessed for in vitro contractility (isometric and isotonic properties) alongside indices of fibre-type cross-sectional area, myosin isoform, and capillarity, and estimated muscle PO2 . In situ extensor digitorum longus (EDL) contractility and femoral blood flow were assessed. HFpEF soleus demonstrated lower absolute maximal force by 22%, fibre atrophy by 24%, a fibre-type shift from I to IIa, and a 17% lower capillary-to-fibre ratio despite increased capillary density (all P < 0.05) with preserved muscle PO2 (P = 0.115) and isometric specific force (P > 0.05). Soleus isotonic properties (shortening velocity and power) were impaired by up to 17 and 22%, respectively (P < 0.05), while the magnitude of the exercise hyperaemia was attenuated by 73% (P = 0.012) in line with higher muscle fatigue by 26% (P = 0.079). Diaphragm alterations (P < 0.05) included Type IIx fibre atrophy despite Type I/IIa fibre hypertrophy, with increased indices of capillarity alongside preserved contractile properties during isometric, isotonic, and cyclical contractions. In conclusion, obese-HFpEF rats demonstrated blunted skeletal muscle blood flow during contractions in parallel to microvascular structural remodelling, fibre atrophy, and isotonic contractile dysfunction in the locomotor muscles. In contrast, diaphragm phenotype remained well preserved. This study identifies numerous muscle-specific impairments that could exacerbate exercise intolerance in obese-HFpEF.

Investigation of the Role of Myocyte Orientations in Cardiac Arrhythmia Using Image-Based Models.

Cardiac electrical excitation-propagation is influenced by myocyte orientations (cellular organization). Quantitatively understanding this relationship presents a significant research challenge, especially during arrhythmias in which excitation patterns become complex. Tissue-scale simulations of cardiac electrophysiology, incorporating both dynamic action potential behavior and image-based myocardial architecture, provide an approach to investigate three-dimensional (3D) propagation of excitation waves in the heart. In this study, we aimed to assess the importance of natural variation in myocyte orientations on cardiac arrhythmogenesis using 3D tissue electrophysiology simulations. Three anatomical models (i.e., describing myocyte orientations) of healthy rat ventricles-obtained using diffusion tensor imaging at 100 µm resolution-were registered to a single biventricular geometry (i.e., a single cardiac shape), in which the myocyte orientations could be represented by each of the diffusion tensor imaging data sets or by an idealized rule-based description. The Fenton-Karma cellular excitation model was modified to reproduce rat ventricular action potential duration restitution to create reaction-diffusion cardiac electrophysiology models. Over 250 3D simulations were performed to investigate the effects of myocyte orientations on the following: 1) ventricular activation, 2) location-dependent arrhythmia induction via rapid pacing, and 3) dynamics of re-entry averaged over multiple episodes. It was shown that 1) myocyte orientation differences manifested themselves in local activation times, but the influence on total activation time was small; 2) differences in myocyte orientations could critically affect the inducibility and persistence of arrhythmias for specific stimulus-location/cycle-length combinations; and 3) myocyte orientations alone could be an important determinant of scroll wave break, although no significant differences were observed in averaged arrhythmia dynamics between the four myocyte orientation scenarios considered. Our results show that myocyte orientations are an important determinant of arrhythmia inducibility, persistence, and scroll wave break. These findings suggest that where specificity is desired (for example, when predicting location-dependent, patient-specific arrhythmia inducibility), subject-specific myocyte orientations may be important.

Effects of Heart Rate and Ventricular Wall Thickness on Non-invasive Mapping: An in silico Study.

Background: Non-invasive cardiac mapping-also known as Electrocardiographic imaging (ECGi)-is a novel, painless and relatively economic method to map the electrical activation and repolarization patterns of the heart, providing a valuable tool for early identification and diagnosis of conduction abnormalities and arrhythmias. Moreover, the ability to obtain information on cardiac electrical activity non-invasively using ECGi provides the potential for a priori information to guide invasive surgical procedures, improving success rates, and reducing procedure time. Previous studies have shown the influence of clinical variables, such as heart rate, heart size, endocardial wall, and body composition on surface electrocardiogram (ECG) measurements. The influence of clinical variables on the ECG variability has provided information on cardiovascular control and its abnormalities in various pathologies. However, the effects of such clinical variables on the Body Surface Potential (BSP) and ECGi maps have yet to be systematically investigated. Methods: In this study we investigated the effects of heart size, intracardiac thickness, and heart rate on BSP and ECGi maps using a previously-developed 3D electrophysiologically-detailed ventricles-torso model. The inverse solution was solved using the three different Tikhonov regularization methods. Results: Through comparison of multiple measures of error/accuracy on the ECGi reconstructions, our results showed that using different heart geometries to solve the forward and inverse problems produced a larger estimated focal excitation location. An increase of ~2 mm in the Euclidean distance error was observed for an increase in the heart size. However, the estimation of the location of focal activity was still able to be obtained. Similarly, a Euclidean distance increase was observed when the order of regularization was reduced. For the case of activation maps reconstructed at the same ectopic focus location but different heart rates, an increase in the errors and Euclidean distance was observed when the heart rate was increased. Conclusions: Non-invasive cardiac mapping can still provide useful information about cardiac activation patterns for the cases when a different geometry is used for the inverse problem compared to the one used for the forward solution; rapid pacing rates can induce order-dependent errors in the accuracy of reconstruction.

Computational Analysis of the Mode of Action of Disopyramide and Quinidine on hERG-Linked Short QT Syndrome in Human Ventricles.

The short QT syndrome (SQTS) is a rare cardiac disorder associated with arrhythmias and sudden death. Gain-of-function mutations to potassium channels mediating the rapid delayed rectifier current, IKr, underlie SQTS variant 1 (SQT1), in which treatment with Na+ and K+ channel blocking class Ia anti-arrhythmic agents has demonstrated some efficacy. This study used computational modeling to gain mechanistic insights into the actions of two such drugs, disopyramide and quinidine, in the setting of SQT1. The O'Hara-Rudy (ORd) human ventricle model was modified to incorporate a Markov chain formulation of IKr describing wild type (WT) and SQT1 mutant conditions. Effects of multi-channel block by disopyramide and quinidine, including binding kinetics and altered potency of IKr/hERG channel block in SQT1 and state-dependent block of sodium channels, were simulated on action potential and multicellular tissue models. A one-dimensional (1D) transmural ventricular strand model was used to assess prolongation of the QT interval, effective refractory period (ERP), and re-entry wavelength (WL) by both drugs. Dynamics of re-entrant excitation waves were investigated using a 3D human left ventricular wedge model. In the setting of SQT1, disopyramide, and quinidine both produced a dose-dependent prolongation in (i) the QT interval, which was primarily due to IKr block, and (ii) the ERP, which was mediated by a synergistic combination of IKr and INa block. Over the same range of concentrations quinidine was more effective in restoring the QT interval, due to more potent block of IKr. Both drugs demonstrated an anti-arrhythmic increase in the WL of re-entrant circuits. In the 3D wedge, disopyramide and quinidine at clinically-relevant concentrations decreased the dominant frequency of re-entrant excitations and exhibited anti-fibrillatory effects; preventing formation of multiple, chaotic wavelets which developed in SQT1, and could terminate arrhythmias. This computational modeling study provides novel insights into the clinical efficacy of disopyramide and quinidine in the setting of SQT1; it also dissects ionic mechanisms underlying QT and ERP prolongation. Our findings show that both drugs demonstrate efficacy in reversing the SQT1 phenotype, and indicate that disopyramide warrants further investigation as an alternative to quinidine in the treatment of SQT1.

Trigger vs. Substrate: Multi-Dimensional Modulation of QT-Prolongation Associated Arrhythmic Dynamics by a hERG Channel Activator.

Background: Prolongation of the QT interval of the electrocardiogram (ECG), underlain by prolongation of the action potential duration (APD) at the cellular level, is linked to increased vulnerability to cardiac arrhythmia. Pharmacological management of arrhythmia associated with QT prolongation is typically achieved through attempting to restore APD to control ranges, reversing the enhanced vulnerability to Ca2+-dependent afterdepolarisations (arrhythmia triggers) and increased transmural dispersion of repolarisation (arrhythmia substrate) associated with APD prolongation. However, such pharmacological modulation has been demonstrated to have limited effectiveness. Understanding the integrative functional impact of pharmacological modulation requires simultaneous investigation of both the trigger and substrate. Methods: We implemented a multi-scale (cell and tissue) in silico approach using a model of the human ventricular action potential, integrated with a model of stochastic 3D spatiotemporal Ca2+ dynamics, and parameter modification to mimic prolonged QT conditions. We used these models to examine the efficacy of the hERG activator MC-II-157c in restoring APD to control ranges, examined its effects on arrhythmia triggers and substrates, and the interaction of these arrhythmia triggers and substrates. Results: QT prolongation conditions promoted the development of spontaneous release events underlying afterdepolarisations during rapid pacing. MC-II-157c applied to prolonged QT conditions shortened the APD, inhibited the development of afterdepolarisations and reduced the probability of afterdepolarisations manifesting as triggered activity in single cells. In tissue, QT prolongation resulted in an increased transmural dispersion of repolarisation, which manifested as an increased vulnerable window for uni-directional conduction block. In some cases, MC-II-157c further increased the vulnerable window through its effects on INa. The combination of stochastic release event modulation and transmural dispersion of repolarisation modulation by MC-II-157c resulted in an integrative behavior wherein the arrhythmia trigger is reduced but the arrhythmia substrate is increased, leading to variable and non-linear overall vulnerability to arrhythmia. Conclusion: The relative balance of reduced trigger and increased substrate underlies a multi-dimensional role of MC-II-157c in modulation of cardiac arrhythmia vulnerability associated with prolonged QT interval.

Dissociating external power from intramuscular exercise intensity during intermittent bilateral knee-extension in humans.

KEY POINTS: Continuous high-intensity constant-power exercise is unsustainable, with maximal oxygen uptake (V̇O2 max ) and the limit of tolerance attained after only a few minutes. Performing the same power intermittently reduces the O2 cost of exercise and increases tolerance. The extent to which this dissociation is reflected in the intramuscular bioenergetics is unknown. We used pulmonary gas exchange and 31 P magnetic resonance spectroscopy to measure whole-body V̇O2, quadriceps phosphate metabolism and pH during continuous and intermittent exercise of different work:recovery durations. Shortening the work:recovery durations (16:32 s vs. 32:64 s vs. 64:128 s vs. continuous) at a work rate estimated to require 110% peak aerobic power reduced V̇O2, muscle phosphocreatine breakdown and muscle acidification, eliminated the glycolytic-associated contribution to ATP synthesis, and increased exercise tolerance. Exercise intensity (i.e. magnitude of intramuscular metabolic perturbations) can be dissociated from the external power using intermittent exercise with short work:recovery durations. ABSTRACT: Compared with work-matched high-intensity continuous exercise, intermittent exercise dissociates pulmonary oxygen uptake (V̇O2) from the accumulated work. The extent to which this reflects differences in O2 storage fluctuations and/or contributions from oxidative and substrate-level bioenergetics is unknown. Using pulmonary gas-exchange and intramuscular 31 P magnetic resonance spectroscopy, we tested the hypotheses that, at the same power: ATP synthesis rates are similar, whereas peak V̇O2 amplitude is lower in intermittent vs. continuous exercise. Thus, we expected that: intermittent exercise relies less upon anaerobic glycolysis for ATP provision than continuous exercise; shorter intervals would require relatively greater fluctuations in intramuscular bioenergetics than in V̇O2 compared to longer intervals. Six men performed bilateral knee-extensor exercise (estimated to require 110% peak aerobic power) continuously and with three different intermittent work:recovery durations (16:32, 32:64 and 64:128 s). Target work duration (576 s) was achieved in all intermittent protocols; greater than continuous (252 ± 174 s; P < 0.05). Mean ATP turnover rate was not different between protocols (∼43 mm min-1 on average). However, the intramuscular phosphocreatine (PCr) component of ATP generation was greatest (∼30 mm min-1 ), and oxidative (∼10 mm min-1 ) and anaerobic glycolytic (∼1 mm min-1 ) components were lowest for 16:32 and 32:64 s intermittent protocols, compared to 64:128 s (18 ± 6, 21 ± 10 and 10 ± 4 mm min-1 , respectively) and continuous protocols (8 ± 6, 20 ± 9 and 16 ± 14 mm min-1 , respectively). As intermittent work duration increased towards continuous exercise, ATP production relied proportionally more upon anaerobic glycolysis and oxidative phosphorylation, and less upon PCr breakdown. However, performing the same high-intensity power intermittently vs. continuously reduced the amplitude of fluctuations in V̇O2 and intramuscular metabolism, dissociating exercise intensity from the power output and work done.

Data collection, handling, and fitting strategies to optimize accuracy and precision of oxygen uptake kinetics estimation from breath-by-breath measurements.

Phase 2 pulmonary oxygen uptake kinetics (ϕ2 τV̇o2P) reflect muscle oxygen consumption dynamics and are sensitive to changes in state of training or health. This study identified an unbiased method for data collection, handling, and fitting to optimize V̇o2P kinetics estimation. A validated computational model of V̇o2P kinetics and a Monte Carlo approach simulated 2 × 105 moderate-intensity transitions using a distribution of metabolic and circulatory parameters spanning normal health. Effects of averaging (interpolation, binning, stacking, or separate fitting of up to 10 transitions) and fitting procedures (biexponential fitting, or ϕ2 isolation by time removal, statistical, or derivative methods followed by monoexponential fitting) on accuracy and precision of V̇o2P kinetics estimation were assessed. The optimal strategy to maximize accuracy and precision of τV̇o2P estimation was 1-s interpolation of 4 bouts, ensemble averaged, with the first 20 s of exercise data removed. Contradictory to previous advice, we found optimal fitting procedures removed no more than 20 s of ϕ1 data. Averaging method was less critical: interpolation, binning, and stacking gave similar results, each with greater accuracy compared with analyzing repeated bouts separately. The optimal procedure resulted in ϕ2 τV̇o2P estimates for transitions from an unloaded or loaded baseline that averaged 1.97 ± 2.08 and 1.04 ± 2.30 s from true, but were within 2 s of true in only 47-62% of simulations. Optimized 95% confidence intervals for τV̇o2P ranged from 4.08 to 4.51 s, suggesting a minimally important difference of ~5 s to determine significant changes in τV̇o2P during interventional and comparative studies.NEW & NOTEWORTHY We identified an unbiased method to maximize accuracy and precision of oxygen uptake kinetics (τV̇o2P) estimation. The optimum number of bouts to average was four; interpolation, bin, and stacking averaging methods gave similar results. Contradictory to previous advice, we found that optimal fitting procedures removed no more than 20 s of phase 1 data. Our data suggest a minimally important difference of ~5 s to determine significant changes in τV̇o2P during interventional and comparative studies.

No reserve in isokinetic cycling power at intolerance during ramp incremental exercise in endurance-trained men.

During whole body exercise in health, maximal oxygen uptake (V̇o2max) is typically attained at or immediately before the limit of tolerance (LoT). At the V̇o2max and LoT of incremental exercise, a fundamental, but unresolved, question is whether maximal evocable power can be increased above the task requirement, i.e., whether there is a "power reserve" at the LoT. Using an instantaneous switch from cadence-independent (hyperbolic) to isokinetic cycle ergometry, we determined maximal evocable power at the limit of ramp-incremental exercise. We hypothesized that in endurance-trained men at LoT, maximal (4 s) isokinetic power would not differ from the power required by the task. Baseline isokinetic power at 80 rpm (Piso; measured at the pedals) and summed integrated EMG from five leg muscles (ΣiEMG) were measured in 12 endurance-trained men (V̇o2max = 4.2 ± 1.0 l/min). Participants then completed a ramp incremental exercise test (20-25 W/min), with instantaneous measurement of Piso and ΣiEMG at the LoT. Piso decreased from 788 ± 103 W at baseline to 391 ± 72 W at LoT, which was not different from the required ramp-incremental flywheel power (352 ± 58 W; P > 0.05). At LoT, the relative reduction in Piso was greater than the relative reduction in the isokinetic ΣiEMG (50 ± 9 vs. 63 ± 10% of baseline; P < 0.05). During maximal ramp incremental exercise in endurance-trained men, maximum voluntary power is not different from the power required by the task and is consequent to both central and peripheral limitations in evocable power. The absence of a power reserve suggests both the perceptual and physiological limits of maximum voluntary power production are not widely dissociated at LoT in this population.

Influence of muscle metabolic heterogeneity in determining the V̇o2p kinetic response to ramp-incremental exercise.

The pulmonary O2 uptake (V̇o2p) response to ramp-incremental (RI) exercise increases linearly with work rate (WR) after an early exponential phase, implying that a single time constant (τ) and gain (G) describe the response. However, variability in τ and G of V̇o2p kinetics to different step increments in WR is documented. We hypothesized that the "linear" V̇o2p-WR relationship during RI exercise results from the conflation between WR-dependent changes in τ and G. Nine men performed three or four repeats of RI exercise (30 W/min) and two step-incremental protocols consisting of four 60-W increments beginning from 20 W or 50 W. During testing, breath-by-breath V̇o2p was measured by mass spectrometry and volume turbine. For each individual, the V̇o2p RI response was characterized with exponential functions containing either constant or variable τ and G values. A relationship between τ and G vs. WR was determined from the step-incremental protocols to derive the variable model parameters. τ and G increased from 21 ± 5 to 98 ± 20 s and from 8.7 ± 0.6 to 12.0 ± 1.9 ml·min(-1)·W(-1) for WRs of 20-230 W, respectively, and were best described by a second-order (τ) and a first-order (G) polynomial function of WR (lowest Akaike information criterion score). The sum of squared residuals was not different (P > 0.05) when the V̇o2p RI response was characterized with either the constant or variable models, indicating that they described the response equally well. Results suggest that τ and G increase progressively with WR during RI exercise. Importantly, these relationships may conflate to produce a linear V̇o2p-WR response, emphasizing the influence of metabolic heterogeneity in determining the apparent V̇o2p-WR relationship during RI exercise.

The influence of metabolic and circulatory heterogeneity on the expression of pulmonary oxygen uptake kinetics in humans.

We examined the relationship amongst baseline work rate (WR), phase II pulmonary oxygen uptake (V̇(O2p)) time constant (τV̇(O2p)) and functional gain (G(P)=ΔV̇(O2p)/ΔWR) during moderate-intensity exercise. Transitions were initiated from a constant or variable baseline WR. A validated circulatory model was used to examine the role of heterogeneity in muscle metabolism (V̇(O2m)) and blood flow (Q̇(m)) in determining V̇(O2p) kinetics. We hypothesized that τV̇(O2p) and G(P) would be invariant in the constant baseline condition but would increase linearly with increased baseline WR. Fourteen men completed three to five repetitions of ∆40 W step transitions initiated from 20, 40, 60, 80, 100 and 120 W on a cycle ergometer. The ∆40 W step transitions from 60, 80, 100 and 120 W were preceded by 6 min of 20 W cycling, from which the progressive ΔWR transitions (constant baseline condition) were examined. The V̇(O2p) was measured breath by breath using mass spectrometry and a volume turbine. For a given ΔWR, both τV̇(O2p) (22-35 s) and G(P) (8.7-10.5 ml min(-1) W(-1)) increased (P < 0.05) linearly as a function of baseline WR (20-120 W). The τV̇(O2p) was invariant (P < 0.05) in transitions initiated from 20 W, but G(P) increased with ΔWR (P < 0.05). Modelling the summed influence of multiple muscle compartments revealed that τV̇(O2p) could appear fast (24 s), and similar to in vivo measurements (22 ± 6 s), despite being derived from τV̇(O2p) values with a range of 15-40 s and τQ̇(m) with a range of 20-45 s, suggesting that within the moderate-intensity domain phase II V̇(O2p) kinetics are slowed dependent on the pretransition WR and are strongly influenced by muscle metabolic and circulatory heterogeneity.

Systems approach to the study of stretch and arrhythmias in right ventricular failure induced in rats by monocrotaline.

We demonstrate the synergistic benefits of using multiple technologies to investigate complex multi-scale biological responses. The combination of reductionist and integrative methodologies can reveal novel insights into mechanisms of action by tracking changes of in vivo phenomena to alterations in protein activity (or vice versa). We have applied this approach to electrical and mechanical remodelling in right ventricular failure caused by monocrotaline-induced pulmonary artery hypertension in rats. We show arrhythmogenic T-wave alternans in the ECG of conscious heart failure animals. Optical mapping of isolated hearts revealed discordant action potential duration (APD) alternans. Potential causes of the arrhythmic substrate; structural remodelling and/or steep APD restitution and dispersion were observed, with specific remodelling of the Right Ventricular Outflow Tract. At the myocyte level, [Ca(2+)]i transient alternans were observed together with decreased activity, gene and protein expression of the sarcoplasmic reticulum Ca(2+)-ATPase (SERCA). Computer simulations of the electrical and structural remodelling suggest both contribute to a less stable substrate. Echocardiography was used to estimate increased wall stress in failure, in vivo. Stretch of intact and skinned single myocytes revealed no effect on the Frank-Starling mechanism in failing myocytes. In isolated hearts acute stretch-induced arrhythmias occurred in all preparations. Significant shortening of the early APD was seen in control but not failing hearts. These observations may be linked to changes in the gene expression of candidate mechanosensitive ion channels (MSCs) TREK-1 and TRPC1/6. Computer simulations incorporating MSCs and changes in ion channels with failure, based on altered gene expression, largely reproduced experimental observations.