People with obesity exhibit losses in muscle proteostasis that are partly improved by exercise training.

This pilot experiment examines if a loss in muscle proteostasis occurs in people with obesity and whether endurance exercise positively influences either the abundance profile or turnover rate of proteins in this population. Men with (n = 3) or without (n = 4) obesity were recruited and underwent a 14-d measurement protocol of daily deuterium oxide (D2 O) consumption and serial biopsies of vastus lateralis muscle. Men with obesity then completed 10-weeks of high-intensity interval training (HIIT), encompassing 3 sessions per week of cycle ergometer exercise with 1 min intervals at 100% maximum aerobic power interspersed by 1 min recovery periods. The number of intervals per session progressed from 4 to 8, and during weeks 8-10 the 14-d measurement protocol was repeated. Proteomic analysis detected 352 differences (p < 0.05, false discovery rate < 5%) in protein abundance and 19 (p < 0.05) differences in protein turnover, including components of the ubiquitin-proteasome system. HIIT altered the abundance of 53 proteins and increased the turnover rate of 22 proteins (p < 0.05) and tended to benefit proteostasis by increasing muscle protein turnover rates. Obesity and insulin resistance are associated with compromised muscle proteostasis, which may be partially restored by endurance exercise.

A type IV functional response with different shapes in a predator-prey model.

Group defense is a phenomenon that occurs in many predator-prey systems. Different functional responses with substantially different properties representing such a mechanism exist. Here, we develop a functional response using timescale separation. A prey-dependent catch rate represents the group defense. The resulting functional response contains a single parameter that controls whether the group defense functional response is saturating or dome-shaped. Based on that, we show that the catch rate must not increase monotonically with increasing prey density to lead to a dome-shaped functional response. We apply bifurcation analysis to show that non-monotonic group defense is usually more successful. However, we also find parameter regions in which a paradox occurs. In this case, higher group defense can give rise to a stable limit cycle, while for lower values, the predator would go extinct. The study does not only provide valuable insight on how to include functional responses representing group defense in mathematical models, but it also clarifies under which circumstances the usage of different functional responses is appropriate.

Mathematical modelling indicates that lower activity of the haemostatic system in neonates is primarily due to lower prothrombin concentration.

Haemostasis is governed by a highly complex system of interacting proteins. Due to the central role of thrombin, thrombin generation and specifically the thrombin generation curve (TGC) is commonly used as an indicator of haemostatic activity. Functional characteristics of the haemostatic system in neonates and children are significantly different compared with adults; at the same time plasma levels of haemostatic proteins vary considerably with age. However, relating one to the other has been difficult, both due to significant inter-individual differences for individuals of similar age and the complexity of the biochemical reactions underlying haemostasis. Mathematical modelling has been very successful at representing the biochemistry of blood clotting. In this study we address the challenge of large inter-individual variability by parameterising the Hockin-Mann model with data from individual patients, across different age groups from neonates to adults. Calculating TGCs for each patient of a specific age group provides us with insight into the variability of haemostatic activity across that age group. From our model we observe that two commonly used metrics for haemostatic activity are significantly lower in neonates than in older patients. Because both metrics are strongly determined by prothrombin and prothrombin levels are considerably lower in neonates we conclude that decreased haemostatic activity in neonates is due to lower prothrombin availability.

Modelling modal gating of ion channels with hierarchical Markov models.

Many ion channels spontaneously switch between different levels of activity. Although this behaviour known as modal gating has been observed for a long time it is currently not well understood. Despite the fact that appropriately representing activity changes is essential for accurately capturing time course data from ion channels, systematic approaches for modelling modal gating are currently not available. In this paper, we develop a modular approach for building such a model in an iterative process. First, stochastic switching between modes and stochastic opening and closing within modes are represented in separate aggregated Markov models. Second, the continuous-time hierarchical Markov model, a new modelling framework proposed here, then enables us to combine these components so that in the integrated model both mode switching as well as the kinetics within modes are appropriately represented. A mathematical analysis reveals that the behaviour of the hierarchical Markov model naturally depends on the properties of its components. We also demonstrate how a hierarchical Markov model can be parametrized using experimental data and show that it provides a better representation than a previous model of the same dataset. Because evidence is increasing that modal gating reflects underlying molecular properties of the channel protein, it is likely that biophysical processes are better captured by our new approach than in earlier models.

Examination of the Effects of Heterogeneous Organization of RyR Clusters, Myofibrils and Mitochondria on Ca2+ Release Patterns in Cardiomyocytes.

Spatio-temporal dynamics of intracellular calcium, [Ca2+]i, regulate the contractile function of cardiac muscle cells. Measuring [Ca2+]i flux is central to the study of mechanisms that underlie both normal cardiac function and calcium-dependent etiologies in heart disease. However, current imaging techniques are limited in the spatial resolution to which changes in [Ca2+]i can be detected. Using spatial point process statistics techniques we developed a novel method to simulate the spatial distribution of RyR clusters, which act as the major mediators of contractile Ca2+ release, upon a physiologically-realistic cellular landscape composed of tightly-packed mitochondria and myofibrils. We applied this method to computationally combine confocal-scale (~ 200 nm) data of RyR clusters with 3D electron microscopy data (~ 30 nm) of myofibrils and mitochondria, both collected from adult rat left ventricular myocytes. Using this hybrid-scale spatial model, we simulated reaction-diffusion of [Ca2+]i during the rising phase of the transient (first 30 ms after initiation). At 30 ms, the average peak of the simulated [Ca2+]i transient and of the simulated fluorescence intensity signal, F/F0, reached values similar to that found in the literature ([Ca2+]i ≈1 µM; F/F0≈5.5). However, our model predicted the variation in [Ca2+]i to be between 0.3 and 12.7 µM (~3 to 100 fold from resting value of 0.1 µM) and the corresponding F/F0 signal ranging from 3 to 9.5. We demonstrate in this study that: (i) heterogeneities in the [Ca2+]i transient are due not only to heterogeneous distribution and clustering of mitochondria; (ii) but also to heterogeneous local densities of RyR clusters. Further, we show that: (iii) these structure-induced heterogeneities in [Ca2+]i can appear in line scan data. Finally, using our unique method for generating RyR cluster distributions, we demonstrate the robustness in the [Ca2+]i transient to differences in RyR cluster distributions measured between rat and human cardiomyocytes.

A kinetic model for type I and II IP3R accounting for mode changes.

Based upon an extensive single-channel data set, a Markov model for types I and II inositol trisphosphate receptors (IP(3)R) is developed. The model aims to represent accurately the kinetics of both receptor types of IP(3)R depending on the concentrations of inositol trisphosphate (IP(3)), adenosine trisphosphate (ATP), and intracellular calcium (Ca(2+)). In particular, the model takes into account that for some combinations of ligands the IP(3)R switches between extended periods of inactivity alternating with intervals of bursting activity (mode changes). In a first step, the inactive and active modes are modeled separately. It is found that, within modes, both receptor types are ligand-independent. In a second step, the submodels are connected by transition rates. Ligand-dependent regulation of the channel activity is achieved by modulating these transitions between active and inactive modes. As a result, a compact representation of the IP(3)R is obtained that accurately captures stochastic single-channel dynamics including mode changes in a model with six states and 10 rate constants, only two of which are ligand-dependent.

MCMC can detect nonidentifiable models.

Continuous-time Markov models have been considered the best representation for the stochastic dynamics of ion channels for more than thirty years. For most single-channel data sets, several open and closed states are required for accurately representing the dynamics. However, each data point only shows if the channel is open or closed but not in which state it is. Consequently, some model structures are inherently overparameterized and therefore, in principle, unsuitable for representing any data--those models are called "nonidentifiable". As of this writing, it seems to be poorly understood which continuous-time Markov models are identifiable and which are not, therefore the unconscious use of a nonidentifiable model is a considerable concern. To address this problem, an improved variant of a recently published Markov-chain Monte Carlo method is presented. The algorithm is tested using test data as well as experimental data. We demonstrate that, opposed to a widely used maximum-likelihood estimator, it gives clear warning signs when a nonidentifiable model is used for fitting. Furthermore, for test data that was generated from a nonidentifiable model, the Markov-chain Monte Carlo results recover much more information from the data than maximum-likelihood estimation.

MCMC estimation of Markov models for ion channels.

Ion channels are characterized by inherently stochastic behavior which can be represented by continuous-time Markov models (CTMM). Although methods for collecting data from single ion channels are available, translating a time series of open and closed channels to a CTMM remains a challenge. Bayesian statistics combined with Markov chain Monte Carlo (MCMC) sampling provide means for estimating the rate constants of a CTMM directly from single channel data. In this article, different approaches for the MCMC sampling of Markov models are combined. This method, new to our knowledge, detects overparameterizations and gives more accurate results than existing MCMC methods. It shows similar performance as QuB-MIL, which indicates that it also compares well with maximum likelihood estimators. Data collected from an inositol trisphosphate receptor is used to demonstrate how the best model for a given data set can be found in practice.

An extension of the Beretta-Kuang model of viral diseases.

A model for the complete life cycle of marine viruses is presented. The Beretta-Kuang model introduces an explicit equation for viral particles but the replication process of viral particles in their hosts is not considered. The extended model keeps the structure of the original model. This makes it possible to estimate the growth parameters of the viruses for a given parametrisation of the Beretta-Kuang model.

Predation may defeat spatial spread of infection.

A model of a phytoplankton-zooplankton prey-predator system with viral infection of phytoplankton is investigated. Virus particles (V) are taken into account by an explicit equation. Phytoplankton is split into a susceptible (S) and an infected (I) class. A lytic infection is considered, thus, infected phytoplankton cells stop reproducing as soon as the infection starts and die at an increased mortality rate. Zooplankton (Z) is grazing on both susceptible and infected phytoplankton following a Holling-type II functional response. After the local dynamics of the V-S-I-Z system is analysed, numerical solutions of a stochastic reaction-diffusion model of the four species are presented. These show a spatial competition between zooplankton and viruses, although these two species are not explicitly coupled by the model equations.