Diffusion restrictions surrounding mitochondria: a mathematical model of heart muscle fibers.
Ontology highlight
ABSTRACT:
SUBMITTER: Ramay HR
PROVIDER: S-EPMC2711342 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Many long skeletal muscles are comprised of fibers that terminate intrafascicularly. Force from terminating fibers can be transmitted through shear within the endomysium that surrounds fibers or through tension within the endomysium that extends from fibers to the tendon; however, it is unclear which pathway dominates in force transmission from terminating fibers. The purpose of this work was to develop mathematical models to (i) compare the efficacy of lateral (through shear) and longitudinal (through tension) force transmission in intrafascicularly terminating fibers, and (ii) determine how force transmission is affected by variations in the structure and properties of fibers and the endomysium. The models demonstrated that even though the amount of force that can be transmitted from an intrafascicularly terminating fiber is dependent on fiber resting length (the unstretched length at which passive stress is zero), endomysium shear modulus, and fiber volume fraction (the fraction of the muscle cross-sectional area that is occupied by fibers), fibers that have values of resting length, shear modulus, and volume fraction within physiologic ranges can transmit nearly all of their peak isometric force laterally through shearing of the endomysium. By contrast, the models predicted only limited force transmission ability through tension within the endomysium that extends from the fiber to the tendon. Moreover, when fiber volume fraction decreases to unhealthy ranges (less than 50%), the force-transmitting potential of terminating fibers through shearing of the endomysium decreases significantly. The models presented here support the hypothesis that lateral force transmission through shearing of the endomysium is an effective mode of force transmission in terminating fibers.
| S-EPMC3134549 | biostudies-literature
Project description:The mammalian heart must function as an efficient pump while simultaneously conducting electrical signals to drive the contraction process. In the ventricles, electrical activation begins at the insertion points of the Purkinje network in the endocardium. How does the diffusion component of the subsequent excitation wave propagate from the endocardium in a healthy heart wall without creating directional biases? We show that this is a consequence of the particular geometric organization of myocytes in the heart wall. Using a generalized helicoid to model fiber orientation, we treat the myocardium as a curved space via Riemannian geometry, and then use stochastic calculus to model local signal diffusion. Our analysis shows that the helicoidal arrangement of myocytes minimizes the directional biases that could lead to aberrant propagation, thereby explaining how electrophysiological principles are consistent with local measurements of cardiac fiber geometry. We discuss our results in the context of the need to balance electrical and mechanical requirements for heart function.
| S-EPMC5940931 | biostudies-literature
Project description:Efflux time courses of endogenous cytosolic proteins were obtained from rabbit psoas muscle fibers skinned in oil and transferred to physiological salt solution. Proteins were separated by gel electrophoresis and compared to load-matched standards for quantitative analysis. A radial diffusion model incorporating the dissociation and dissipation of supramolecular complexes accounts for an initial lag and subsequent efflux of glycolytic and glycogenolytic enzymes. The model includes terms representing protein crowding, myofilament lattice hindrance, and binding to the cytomatrix. Optimization algorithms returned estimates of the apparent diffusion coefficients, D(r,t), that were very low at the onset of diffusion (∼10(-10) cm(2) s(-1)) but increased with time as cytosolic protein density, which was initially high, decreased. D(r,t) at later times ranged from 2.11 × 10(-7) cm(2) s(-1) (parvalbumin) to 0.20 × 10(-7) cm(2) s(-1) (phosphofructose kinase), values that are 3.6- to 12.3-fold lower than those predicted in bulk water. The low initial values are consistent with the presence of complexes in situ; the higher later values are consistent with molecular sieving and transient binding of dissociated proteins. Channeling of metabolic intermediates via enzyme complexes may enhance production of adenosine triphosphate at rates beyond that possible with randomly and/or sparsely distributed enzymes, thereby matching supply with demand.
| S-EPMC3944625 | biostudies-literature
Project description:Variations of free calcium concentration ([Ca2+]) are powerful intracellular signals, controlling contraction as well as metabolism in muscle cells. To fully understand the role of calcium redistribution upon excitation and contraction in skeletal muscle cells, the local [Ca2+] in different compartments needs to be taken into consideration. Fluorescent probes allow the determination of [Ca2+] in the cytosol where myofibrils are embedded, the lumen of the sarcoplasmic reticulum (SR) and the mitochondrial matrix. Previously, models have been developed describing intracellular calcium handling in skeletal and cardiac muscle cells. However, a comprehensive model describing the kinetics of the changes in free calcium concentration in these three compartments is lacking. We designed a new 3D compartmental model of the half sarcomere with radial symmetry, which accounts for diffusion of Ca2+ into the three compartments and simulates its dynamics at rest and at various rates of stimulation in mice skeletal muscle fibers. This model satisfactorily reproduces both the amplitude and time course of the variations of [Ca2+] in the three compartments in mouse fast fibers. As an illustration of the applicability of the model, we investigated the effects of Calsequestrin (CSQ) ablation. CSQ is the main Ca2+ buffer in the SR, localized in close proximity of its calcium release sites and near to the mitochondria. CSQ knock-out mice muscles still preserve a near-normal contractile behavior, but it is unclear whether this is caused by additional SR calcium buffering or a significant contribution of calcium entry from extracellular space, via stored-operated calcium entry (SOCE). The model enabled quantitative assessment of these two scenarios by comparison to measurements of local calcium in the cytosol, the SR and the mitochondria. In conclusion, the model represents a useful tool to investigate the impact of protein ablation and of pharmacological interventions on intracellular calcium dynamics in mice skeletal muscle.
| S-EPMC6062086 | biostudies-literature
Project description:Multiple cellular organelles tightly orchestrate intracellular calcium (Ca2+) dynamics to regulate cellular activities and maintain homeostasis. The interplay between the endoplasmic reticulum (ER), a major store of intracellular Ca2+, and mitochondria, an important source of adenosine triphosphate (ATP), has been the subject of much research, as their dysfunction has been linked with metabolic diseases. Interestingly, throughout the cell's cytosolic domain, these two organelles share common microdomains called mitochondria-associated ER membranes (MAMs), where their membranes are in close apposition. The role of MAMs is critical for intracellular Ca2+ dynamics as they provide hubs for direct Ca2+ exchange between the organelles. A recent experimental study reported correlation between obesity and MAM formation in mouse liver cells, and obesity-related cellular changes that are closely associated with the regulation of Ca2+ dynamics. We constructed a mathematical model to study the effects of MAM Ca2+ dynamics on global Ca2+ activities. Through a series of model simulations, we investigated cellular mechanisms underlying the altered Ca2+ dynamics in the cells under obesity. We predict that, as the dosage of stimulus gradually increases, liver cells from obese mice will reach the state of saturated cytosolic Ca2+ concentration at a lower stimulus concentration, compared to cells from healthy mice.
| S-EPMC6726250 | biostudies-literature
Project description:Thrombin is an enzyme that plays many important roles in the blood clotting process; it activates platelets, cleaves coagulation proteins within feedback loops, and cleaves fibrinogen into fibrin, which polymerizes into fibers to form a stabilizing gel matrix in and around growing clots. Thrombin also binds to the formed fibrin matrix, but this interaction is not well understood. Thrombin-fibrin binding is often described as two independent, single-step binding events, one high-affinity and one low-affinity. However, kinetic schemes describing these single-step binding events do not explain experimentally-observed residency times of fibrin-bound thrombin. In this work, we study a bivalent, sequential-step binding scheme as an alternative to the high-affinity event and, in addition to the low-affinity one. We developed mathematical models for the single- and sequential-step schemes consisting of reaction-diffusion equations to compare to each other and to experimental data. We then used Bayesian inference, in the form of Markov chain Monte Carlo, to learn model parameter distributions from previously published experimental data. For the model to best fit the data, we made an additional assumption that thrombin was irreversibly sequestered; we hypothesized that this could be due to thrombin becoming physically trapped within fibrin fibers as they formed. We further estimated that ?30% of thrombin in the experiments to which we compare our model output became physically trapped. The notion of physically trapped thrombin may provide new insights into conflicting observations regarding the speed of fibrinolysis. Finally, we show that our new model can be used to further probe scenarios dealing with thrombin allostery.
| S-EPMC6817721 | biostudies-literature
Project description:Williams-Beuren syndrome (WBS) is a rare neurodevelopmental disorder that, together with a rather characteristic neurocognitive profile, presents a strong cardiovascular phenotype. The cardiovascular features of WBS are mainly related to a gene dosage effect due to hemizygosity of the elastin (ELN) gene; however, the phenotypic variability between WBS patients indicates the presence of important modulators of the clinical impact of elastin deficiency. Recently, two genes within the WBS region have been linked to mitochondrial dysfunction. Numerous cardiovascular diseases are related to mitochondrial dysfunction; therefore, it could be a modulator of the phenotype present in WBS. Here, we analyze mitochondrial function and dynamics in cardiac tissue from a WBS complete deletion (CD) model. Our research reveals that cardiac fiber mitochondria from CD animals have altered mitochondrial dynamics, accompanied by respiratory chain dysfunction with decreased ATP production, reproducing alterations observed in fibroblasts from WBS patients. Our results highlight two major factors: on the one hand, that mitochondrial dysfunction is probably a relevant mechanism underlying several risk factors associated with WBS disease; on the other, the CD murine model mimics the mitochondrial phenotype of WBS and could be a great model for carrying out preclinical tests on drugs targeting the mitochondria.
| S-EPMC10298715 | biostudies-literature
Project description:In a previous study, the three-dimensional structures of mitochondria in type I and type IIb muscle fibers of chicken were analyzed. The study reported differences in the shape of the mitochondria and the distribution of lipid droplets. In this study, we three-dimensionally analyzed mitochondria and lipid droplets of type II muscle fiber subtypes IIa, IIb, and IIc of chicken lateral iliotibial muscle in the same field of view using correlative light electron microscopy (CLEM) and array tomography methods. The reconstructed images showed that the mitochondria of type IIa muscle fiber were thick and aligned along the myofibrils, and many lipid droplets were embedded in the mitochondria. The mitochondria of type IIb muscle fibers were intermittent, aligned along the myofibrils, and showed contact between adjacent horizontal mitochondria. No lipid droplets were observed in type IIb muscle fiber. In type IIc muscle fiber, we observed irregularly shaped mitochondria with small diameters aligned along the myofibrils. Lipid droplets not only were embedded in the mitochondria but also existed independently in some cases. The combination of array tomography and CLEM methods enabled three-dimensional electron microscopic observation of mitochondria in different subtypes of type II muscle fibers. The subtypes of type II muscle fibers differed in mitochondrial occupancy and morphology and in lipid droplet distribution, and characteristics that had been demonstrated biochemically were also demonstrated ultrastructurally.
| S-EPMC9246695 | biostudies-literature
Project description:This article provides models and code for numerically simulating muscle-fluid-structure interactions (FSIs). This work was presented as part of the symposium on Leading Students and Faculty to Quantitative Biology through Active Learning at the society-wide meeting of the Society for Integrative and Comparative Biology in 2015. Muscle mechanics and simple mathematical models to describe the forces generated by muscular contractions are introduced in most biomechanics and physiology courses. Often, however, the models are derived for simplifying cases such as isometric or isotonic contractions. In this article, we present a simple model of the force generated through active contraction of muscles. The muscles' forces are then used to drive the motion of flexible structures immersed in a viscous fluid. An example of an elastic band immersed in a fluid is first presented to illustrate a fully-coupled FSI in the absence of any external driving forces. In the second example, we present a valveless tube with model muscles that drive the contraction of the tube. We provide a brief overview of the numerical method used to generate these results. We also include as Supplementary Material a MATLAB code to generate these results. The code was written for flexibility so as to be easily modified to many other biological applications for educational purposes.
| S-EPMC6296398 | biostudies-literature
Project description:The unexpectable variations of the diagnosed disease symptoms are quite often observed during medical diagnosis. In stochastics, such behavior is called "grey swan" or "black swan" as synonyms of sudden, unpredictable change. Evolution of the disease's symptoms is usually described by means of Markov processes, where dependency on process history is neglected. The common expectation is that such processes are Gaussian. It is demonstrated here that medical observation can be described as a Markov process and is non-Gaussian. Presented non-Gaussian processes have "fat tail" probability density distribution (pdf). "Fat tail" permits a slight change of probability density distribution and triggers an unexpectable big variation of the diagnosed parameter. Such "fat tail" solution is delivered by the anomalous diffusion model applied here to describe disease evolution and to explain the possible presence of "swans" mentioned above. The proposed model has been obtained as solution of the Fractal Fokker-Planck equation (FFPE). The paper shows a comparison of the results of the theoretical model of anomalous diffusion with experimental results of clinical studies using bioimpedance measurements in cardiology. This allows us to consider the practical usefulness of the proposed solutions.
| S-EPMC9777956 | biostudies-literature