Predicting the effects of dATP on cardiac contraction using multiscale modeling of the sarcomere.
Ontology highlight
ABSTRACT:
SUBMITTER: McCabe KJ
PROVIDER: S-EPMC7721865 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We developed a model of cardiac sarcomere contraction to study the calcium-tension relationship in cardiac muscle. Calcium mediates cardiac contraction through its interactions with troponin (Tn) and subsequently tropomyosin molecules. Experimental studies have shown that a slight increase in intracellular calcium concentration leads to a rapid increase in sarcomeric tension. Though it is widely accepted that the rapid increase is not possible without the concept of cooperativity, the mechanism is debated. We use the hypothesis that there exists a base level of cooperativity intrinsic to the thin filament that is boosted by mechanical tension, i.e. a high level of mechanical tension in the thin filament impedes the unbinding of calcium from Tn. To test these hypotheses, we developed a computational model in which a set of three parameters and inputs of calcium concentration and sarcomere length result in output tension. Tension as simulated appeared in good agreement with experimentally measured tension. Our results support the hypothesis that high tension in the thin filament impedes Tn deactivation by increasing the energy required to detach calcium from the Tn. Given this hypothesis, the model predicted that the areas with highest tension, i.e. closest to the Z-disk end of the single overlap region, show the largest concentration of active Tn's.
| S-EPMC5055322 | biostudies-other
Project description:Cardiac development is a complex process that is strictly controlled by various factors, including PcG protein complexes. Several studies have reported the critical role of PRC2 in cardiogenesis. However, little is known about the regulation mechanism of PRC1 in embryonic heart development. To gain more insight into the mechanistic role of PRC1 in cardiogenesis, we generated a PRC1 loss-of-function zebrafish line by using the CRISPR/Cas9 system targeting rnf2, a gene encoding the core subunit shared by all PRC1 subfamilies. Our results revealed that Rnf2 is not involved in cardiomyocyte differentiation and heart tube formation, but that it is crucial to maintaining regular cardiac contraction. Further analysis suggested that Rnf2 loss-of-function disrupted cardiac sarcomere assembly through the ectopic activation of non-cardiac sarcomere genes in the developing heart. Meanwhile, Rnf2 deficiency disrupts the construction of the atrioventricular canal and the sinoatrial node by modulating the expression of bmp4 and other atrioventricular canal marker genes, leading to an impaired cardiac conduction system. The disorganized cardiac sarcomere and defective cardiac conduction system together contribute to defective cardiac contraction. Our results emphasize the critical role of PRC1 in the cardiac development.
| S-EPMC8583368 | biostudies-literature
Project description:Understanding the dynamics of a cardiac muscle twitch contraction is complex because it requires a detailed understanding of the kinetic processes of the Ca2+ transient, thin-filament activation, and the myosin-actin cross-bridge chemomechanical cycle. Each of these steps has been well defined individually, but understanding how all three of the processes operate in combination is a far more complex problem. Computational modeling has the potential to provide detailed insight into each of these processes, how the dynamics of each process affect the complexity of contractile behavior, and how perturbations such as mutations in sarcomere proteins affect the complex interactions of all of these processes. The mechanisms involved in relaxation of tension during a cardiac twitch have been particularly difficult to discern due to nonhomogeneous sarcomere lengthening during relaxation. Here we use the multiscale MUSICO platform to model rat trabecular twitches. Validation of computational models is dependent on being able to simulate different experimental datasets, but there has been a paucity of data that can provide all of the required parameters in a single experiment, such as simultaneous measurements of force, intracellular Ca2+ transients, and sarcomere length dynamics. In this study, we used data from different studies collected under similar experimental conditions to provide information for all the required parameters. Our simulations established that twitches either in an isometric sarcomere or in fixed-length, multiple-sarcomere trabeculae replicate the experimental observations if models incorporate a length-tension relationship for the nonlinear series elasticity of muscle preparations and a scheme for thick-filament regulation. The thick-filament regulation assumes an off state in which myosin heads are parked onto the thick-filament backbone and are unable to interact with actin, a state analogous to the super-relaxed state. Including these two mechanisms provided simulations that accurately predict twitch contractions over a range of different conditions.
| S-EPMC7852458 | biostudies-literature
Project description:The instantaneous sarcomere length (SL) is regarded as an important indicator of the functional properties of striated muscle. Previously, we found greater sarcomere elongations at the distal end compared to the mid-portion in the mouse tibialis anterior (TA) when the muscle was stretched passively. Here, we wanted to see if SL dispersions increase with activation, as has been observed in single myofibrils, and if SL dispersions differ for different locations in a muscle. Sarcomere lengths were measured at a mid- and a distal location of the TA in live mice using second harmonic generation imaging. Muscle force was measured using a tendon force transducer. We found that SL dispersions increased substantially from the passive to the active state, and were the same for the mid- and distal portions of TA. Sarcomere length non-uniformities within a segment of ~30 serial sarcomeres were up to 1.0?µm. We conclude from these findings that passive, mean SLs obtained from a single location are not necessarily representative of the distribution of SL in active muscle, and thus may be misinterpreted when deriving muscle mechanical properties, such as the force-length relationship. In view of these findings, it seems crucial to determine how SL distributions within a muscle relate to the most fundamental properties of muscle, such as the maximal isometric force.
| S-EPMC6189036 | biostudies-literature
Project description:Sarcomeres, the basic contractile units of striated muscle, produce the forces driving muscular contraction through cross-bridge interactions between actin-containing thin filaments and myosin II-based thick filaments. Until now, direct visualization of the molecular architecture underlying sarcomere contractility has remained elusive. Here, we use in situ cryo-electron tomography to unveil sarcomere contraction in frozen-hydrated neonatal rat cardiomyocytes. We show that the hexagonal lattice of the thick filaments is already established at the neonatal stage, with an excess of thin filaments outside the trigonal positions. Structural assessment of actin polarity by subtomogram averaging reveals that thin filaments in the fully activated state form overlapping arrays of opposite polarity in the center of the sarcomere. Our approach provides direct evidence for thin filament sliding during muscle contraction and may serve as a basis for structural understanding of thin filament activation and actomyosin interactions inside unperturbed cellular environments.
| S-EPMC8253822 | biostudies-literature
Project description:Cardiovascular disease is often related to defects of subcellular components in cardiac myocytes, specifically in the dyadic cleft, which include changes in cleft geometry and channel placement. Modeling of these pathological changes requires both spatially resolved cleft as well as whole cell level descriptions. We use a multiscale model to create dyadic structure-function relationships to explore the impact of molecular changes on whole cell electrophysiology and calcium cycling. This multiscale model incorporates stochastic simulation of individual L-type calcium channels and ryanodine receptor channels, spatially detailed concentration dynamics in dyadic clefts, rabbit membrane potential dynamics, and a system of partial differential equations for myoplasmic and lumenal free Ca2+ and Ca2+-binding molecules in the bulk of the cell. We found action potential duration, systolic, and diastolic [Ca2+] to respond most sensitively to changes in L-type calcium channel current. The ryanodine receptor channel cluster structure inside dyadic clefts was found to affect all biomarkers investigated. The shape of clusters observed in experiments by Jayasinghe et al. and channel density within the cluster (characterized by mean occupancy) showed the strongest correlation to the effects on biomarkers.
| S-EPMC6990380 | biostudies-literature
Project description:Theoretical prediction of electronic absorption spectra without input from experiments is no easy feat, as it requires addressing all of the factors that affect line shapes. In practice, however, the methodologies are limited to treat these ingredients only to a certain extent. Here, we present a multiscale protocol that addresses the temperature, solvent, and nuclear quantum effects as well as anharmonicity and the reconstruction of the final spectra from individual transitions. First, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics is conducted to obtain trajectories of solute-solvent configurations, from which the corresponding quantum-corrected ensembles are generated through the generalized smoothed trajectory analysis (GSTA). The optical spectra of the ensembles are then produced by calculating vertical transitions using time-dependent density-functional theory (TDDFT) with implicit solvation. To obtain the final spectral shapes, the stick spectra from TDDFT are convoluted with Gaussian kernels where the half-widths are determined by a statistically motivated strategy. We have tested our method by calculating the UV-vis spectra of a recently discovered acridine photocatalyst in two redox states. Vibronic progressions and broadenings due to the finite lifetime of the excited states are not included in the methodology yet. Nuclear quantization affects the relative peak intensities and widths, which is necessary to reproduce the experimental spectrum. We have also found that using only the optimized geometry of each molecule works surprisingly well if a proper empirical broadening factor is applied. This is explained by the rigidity of the conjugated chromophore moieties of the selected molecules, which are mainly responsible for the excitations in the spectra. In contrast, we have also shown that other parts of the molecules are flexible enough to feature anharmonicities that impair the use of other techniques such as Wigner sampling.
| S-EPMC8515811 | biostudies-literature
Project description:Motivation: Monitoring, assessment and prediction of environmental risks that chemicals pose demand rapid and accurate diagnostic assays. A variety of toxicological effects have been associated with explosive compounds TNT, RDX and HMX. One important goal of microarray experiments is to discover novel biomarker genes for quantitative phenotypic prediction. We have developed an earthworm microarray containing 15,208 unique oligo probes. Our objective was to identify biomarker genes that can be used to quantitatively predict earthworm tissue residues of the explosives compounds that they were exposed to and took in from the HMX-spiked soil. Results: We collected a large microarray gene expression and earthworm tissue residue dataset. First, differentially expressed genes were identified for each exposure duration (4, 14 and 28 days). These genes were used in multivariate regression modeling for HMX residue prediction. Eighteen different regression models were tested and compared. The best performing model was able to achieve very high prediction accuracies with R2 values of 0.715, 0.728 and 0.822 for 4 days, 14 days and 28 days exposures, separately. Conclusions: This study demonstrated that multivariate regression coupled with high throughput microarray gene expression was a promising approach to quantitative phenotypic prediction. Adult earthworms (Eisenia fetida) were exposed in soil spiked with HMX (0, 8, 16, 32, 64, or 128 mg/kg) for 4, 14 or 28 days. Each treatment originally had 10 replicate worms with all 10 worms survived at the end of exposure. Total RNA was isolated from the surviving worms. A total of 120 worm RNA samples were hybridized to a custom-designed oligo array using AgilentM-bM-^@M-^Ys one-color Low RNA Input Linear Amplification Kit. The array contains 15,208 non-redundant 60-mer probes, each targeting a unique E. fetida transcript. After hybridization and scanning, gene expression data were acquired using AgilentM-bM-^@M-^Ys Feature Extraction Software (v.9.1.3). The 120-array dataset consists of three exposure groups (4 day, 14 day and 28 day) with each group having the following 8 treatments: day 0 control, blank control, solvent control, 8, 16, 32, 64, and 128 mg HMX/g soil. Each treatment had 5 biological replicates (i.e., five individual worms).
2013-08-10 | E-GEOD-42866 | biostudies-arrayexpress
Project description:Background and purposeHeart failure can reflect impaired contractile function at the myofilament level. In healthy hearts, myofilaments become more sensitive to Ca2+ as cells are stretched. This represents a fundamental property of the myocardium that contributes to the Frank-Starling response, although the molecular mechanisms underlying the effect remain unclear. Mavacamten, which binds to myosin, is under investigation as a potential therapy for heart disease. We investigated how mavacamten affects the sarcomere-length dependence of Ca2+ -sensitive isometric contraction to determine how mavacamten might modulate the Frank-Starling mechanism.Experimental approachMulticellular preparations from the left ventricular-free wall of hearts from organ donors were chemically permeabilized and Ca2+ activated in the presence or absence of 0.5-μM mavacamten at 1.9 or 2.3-μm sarcomere length (37°C). Isometric force and frequency-dependent viscoelastic myocardial stiffness measurements were made.Key resultsAt both sarcomere lengths, mavacamten reduced maximal force and Ca2+ sensitivity of contraction. In the presence and absence of mavacamten, Ca2+ sensitivity of force increased as sarcomere length increased. This suggests that the length-dependent activation response was maintained in human myocardium, even though mavacamten reduced Ca2+ sensitivity. There were subtle effects of mavacamten reducing force values under relaxed conditions (pCa 8.0), as well as slowing myosin cross-bridge recruitment and speeding cross-bridge detachment under maximally activated conditions (pCa 4.5).Conclusion and implicationsMavacamten did not eliminate sarcomere length-dependent increases in the Ca2+ sensitivity of contraction in myocardial strips from organ donors at physiological temperature. Drugs that modulate myofilament function may be useful therapies for cardiomyopathies.
| S-EPMC7707091 | biostudies-literature
Project description:Thick-filament sarcomere mutations are a common cause of hypertrophic cardiomyopathy (HCM), a disorder of heart muscle thickening associated with sudden cardiac death and heart failure, with unclear mechanisms. We engineered four isogenic induced pluripotent stem cell (iPSC) models of ?-myosin heavy chain and myosin-binding protein C3 mutations, and studied iPSC-derived cardiomyocytes in cardiac microtissue assays that resemble cardiac architecture and biomechanics. All HCM mutations resulted in hypercontractility with prolonged relaxation kinetics in proportion to mutation pathogenicity, but not changes in calcium handling. RNA sequencing and expression studies of HCM models identified p53 activation, oxidative stress, and cytotoxicity induced by metabolic stress that can be reversed by p53 genetic ablation. Our findings implicate hypercontractility as a direct consequence of thick-filament mutations, irrespective of mutation localization, and the p53 pathway as a molecular marker of contraction stress and candidate therapeutic target for HCM patients.
| S-EPMC6335568 | biostudies-literature