Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a ...laudable goal is a better understanding of the structure–function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.
Purpose. Early animal studies suggest that parotid gland (PG) toxicity prediction could be improved by an accurate estimation of the radiation dose to sub-regions of the PG. Translation to clinical ...investigation requires voxel-level dose accumulation in this organ that responds volumetrically throughout treatment. To date, deformable image registration (DIR) has been evaluated for the PG using only surface alignment. We sought to develop and evaluate an advanced DIR technique capable of modeling these complex PG volume changes over the course of radiation therapy. Materials and methods. Planning and mid-treatment magnetic resonance images from 19 patients and computed tomography images from nine patients who underwent radiation therapy for head and neck cancer were retrospectively evaluated. A finite element model (FEM)-based DIR algorithm was applied between the corresponding pairs of images, based on boundary conditions on the PG surfaces only (Morfeus-spatial). To investigate an anticipated improvement in accuracy, we added a population model-based thermal expansion coefficient to simulate the dose distribution effect on the volume change inside the glands (Morfeus-spatialDose). The model accuracy was quantified using target registration error for magnetic resonance images, where corresponding anatomical landmarks could be identified. The potential clinical impact was evaluated using differences in mean dose, median dose, D98, and D50 of the PGs. Results. In the magnetic resonance images, the mean (±standard deviation) target registration error significantly reduced by 0.25 ± 0.38 mm (p = 0.01) when using Morfeus-spatialDose instead of Morfeus-spatial. In the computed tomography images, differences in the mean dose, median dose, D98, and D50 of the PGs reached 2.9 ± 0.8, 3.8, 4.1, and 3.8 Gy, respectively, between Morfeus-spatial and Morfeus-spatialDose. Conclusion. Differences between Morfeus-spatial and Morfeus-spatialDose may be impactful when considering high-dose gradients of radiation in the PGs. The proposed DIR model can allow more accurate PG alignment than the standard model and improve dose estimation and toxicity prediction modeling.
Prior papers have been inconsistent regarding how much creatinine clearance (CrCl) overestimates glomerular filtration rate (GFR). A recent cross-sectional study suggested that measurement error ...alone could entirely account for the longstanding observation that CrCl/GFR ratio is larger when GFR is lower among patients with chronic kidney disease (CKD); but there have been no validation of this in other cohorts.
To fill these gaps in knowledge regarding the relation between CrCl and GFR, we conducted cross-sectional and longitudinal analysis of the Modification of Diet in Renal Disease study (MDRD) and African American Study of Kidney Disease and Hypertension (AASK); and cross-sectional analysis of a clinical dataset from the Mayo Clinic of four different patient populations (CKD patients, kidney transplant recipients, post kidney donation subgroup and potential kidney donors). In the cross-sectional analyses (MDRD, AASK and Mayo Clinic cohort), we examined the relation between the CrCl/iothalamate GFR (iGFR) ratio at different categories of iGFR or different levels of CrCl. In the MDRD and AASK longitudinal analyses, we studied how the CrCl/iGFR ratio changed with those who had improvement in iGFR (CrCl) over time versus those who had worsening of iGFR (CrCl) over time.
Observed CrCl/iGFR ratios were generally on the lower end of the range reported in the literature for CKD (median 1.24 in MDRD, 1.13 in AASK and 1.25 in Mayo Clinic cohort). Among CKD patients in whom CrCl and iGFR were measured using different timed urine collections, CrCl/iGFR ratio were higher with lower iGFR categories but lower with lower CrCl categories. However, among CKD patients in whom CrCl and iGFR were measured using the same timed urine collections (which reduces dis-concordant measurement error), CrCl/iGFR ratio were higher with both lower iGFR categories and lower CrCl categories.
These data refute the recent suggestion that measurement error alone could entirely account for the longstanding observation that CrCl/GFR ratio increases as GFR decreases in CKD patients. They also highlight the lack of certainty in our knowledge with regard to how much CrCl actually overestimates GFR.
Electrical impulse generation and its conduction within cells or cellular networks are the cornerstone of electrophysiology. However, the advancement of the field is limited by sensing accuracy and ...the scalability of current recording technologies. Here we describe a scalable platform that enables accurate recording of transmembrane potentials in electrogenic cells. The platform employs a three-dimensional high-performance field-effect transistor array for minimally invasive cellular interfacing that produces faithful recordings, as validated by the gold standard patch clamp. Leveraging the high spatial and temporal resolutions of the field-effect transistors, we measured the intracellular signal conduction velocity of a cardiomyocyte to be 0.182 m s
, which is about five times the intercellular velocity. We also demonstrate intracellular recordings in cardiac muscle tissue constructs and reveal the signal conduction paths. This platform could provide new capabilities in probing the electrical behaviours of single cells and cellular networks, which carries broad implications for understanding cellular physiology, pathology and cell-cell interactions.
Although a truly complete understanding of whole heart activation, contraction, and deformation is well beyond our current reach, a significant amount of effort has been devoted to discovering and ...understanding the mechanisms by which myocardial structure determines cardiac function to better treat patients with cardiac disease. Several experimental studies have shown that transmural fiber strain is relatively uniform in both diastole and systole, in contrast to predictions from traditional mechanical theory. Similarly, mathematical models have largely predicted uniform fiber stress across the wall. The development of this uniform pattern of fiber stress and strain during filling and ejection is due to heterogeneous transmural distributions of several myocardial structures. This review summarizes these transmural gradients, their contributions to fiber mechanics, and the potential functional effects of their remodeling during pressure overload hypertrophy.
Mechanically-induced alterations in cardiac electrophysiology are referred to as mechano-electric feedback (MEF), and play an important role in electrical regulation of cardiac performance. The ...influence of mechanical stress and strain on electrophysiology has been investigated at all levels, however the role of MEF in arrhythmia remains poorly understood. During the normal contraction of the heart, mechano-sensitive processes are an implicit component of cardiac activity. Under abnormal mechanical events, stretch-activated mechanisms may contribute to local or global changes in electrophysiology (EP). While such mechanisms have been hypothesised to be involved in mechanically-initiated arrhythmias, the details of these mechanisms and their importance remain elusive.
We assess the theoretical role of stretch mechanisms using coupled models of cellular electrophysiology and sarcomere contraction dynamics. Using models of single ventricular myocytes, we first investigated the potential MEF contributions of stretch-activated currents (SAC), and stretch-induced myofilament calcium release, to test how strain and fibrosis may alter cellular electrophysiology. For all models investigated, SACs were alone not sufficient to create a pro-arrhythmic perturbation of the action potential with stretch. However, when combined with stretch-induced myofilament calcium release, the action potential could be shortened depending on the timing of the strain. This effect was highly model dependent, with a canine epicardial EP model being the most sensitive.
These model results suggest that known mechanisms of mechano-electric coupling in cardiac myocyte may be sufficient to be pro-arrhythmic, but only in combination and under specific strain patterns.
cTnIP82S (cTnIP83S in rodents) resides at the I-T arm of cardiac troponin I (cTnI) and was initially identified as a disease-causing mutation of hypertrophic cardiomyopathy (HCM). However, later ...studies suggested this may not be true. We recently reported that introduction of an HCM-associated mutation in either inhibitory-peptide (cTnIR146G) or cardiac-specific N-terminus (cTnIR21C) of cTnI blunts the PKA-mediated modulation on myofibril activation/relaxation kinetics by prohibiting formation of intrasubunit contacts between these regions. Here, we tested whether this also occurs for cTnIP83S. cTnIP83S increased both Ca2+ binding affinity to cTn (K Ca) and affinity of cTnC for cTnI (K C–I), and eliminated the reduction of K Ca and K C–I observed for phosphorylated-cTnIWT. In isolated myofibrils, cTnIP83S maintained maximal tension (TMAX) and Ca2+ sensitivity of tension (pCa50). For cTnIWT myofibrils, PKA-mediated phosphorylation decreased pCa50 and sped up the slow-phase relaxation (especially for those Ca2+ conditions that heart performs in vivo). Those effects were blunted for cTnIP83S myofibrils. Molecular-dynamics simulations suggested cTnIP83S moderately inhibited an intrasubunit interaction formation between inhibitory-peptide and N-terminus, but this “blunting” effect was weaker than that with cTnIR146G or cTnIR21C. In summary, cTnIP83S has similar effects as other HCM-associated cTnI mutations on troponin and myofibril function even though it is in the I-T arm of cTnI.
Familial cardiomyopathy is a precursor of heart failure and sudden cardiac death. Over the past several decades, researchers have discovered numerous gene mutations primarily in sarcomeric and ...cytoskeletal proteins causing two different disease phenotypes: hypertrophic (HCM) and dilated (DCM) cardiomyopathies. However, molecular mechanisms linking genotype to phenotype remain unclear. Here, we employ a systems approach by integrating experimental findings from preclinical studies (e.g., murine data) into a cohesive signaling network to scrutinize genotype to phenotype mechanisms. We developed an HCM/DCM signaling network model utilizing a logic-based differential equations approach and evaluated model performance in predicting experimental data from four contexts (HCM, DCM, pressure overload, and volume overload). The model has an overall prediction accuracy of 83.8%, with higher accuracy in the HCM context (90%) than DCM (75%). Global sensitivity analysis identifies key signaling reactions, with calcium-mediated myofilament force development and calcium-calmodulin kinase signaling ranking the highest. A structural revision analysis indicates potential missing interactions that primarily control calcium regulatory proteins, increasing model prediction accuracy. Combination pharmacotherapy analysis suggests that downregulation of signaling components such as calcium, titin and its associated proteins, growth factor receptors, ERK1/2, and PI3K-AKT could inhibit myocyte growth in HCM. In experiments with patient-specific iPSC-derived cardiomyocytes (MLP-W4R;MYH7-R723C iPSC-CMs), combined inhibition of ERK1/2 and PI3K-AKT rescued the HCM phenotype, as predicted by the model. In DCM, PI3K-AKT-NFAT downregulation combined with upregulation of Ras/ERK1/2 or titin or Gq protein could ameliorate cardiomyocyte morphology. The model results suggest that HCM mutations that increase active force through elevated calcium sensitivity could increase ERK activity and decrease eccentricity through parallel growth factors, Gq-mediated, and titin pathways. Moreover, the model simulated the influence of existing medications on cardiac growth in HCM and DCM contexts. This HCM/DCM signaling model demonstrates utility in investigating genotype to phenotype mechanisms in familial cardiomyopathy.
Display omitted
•HCM/DCM signaling network model links mutations to cardiomyocyte morphological changes.•The computational model predicts higher efficiency for combination pharmacotherapies in familial cardiomyopathy.•Experiments on iPSC-CMs showed rescue of HCM phenotype by combined inhibition of ERK1/2 and PI3K-AKT predicted by the model.
Pressure overload left ventricular (LV) hypertrophy is characterized by increased cardiomyocyte width and ventricle wall thickness, however the regional variation of this remodeling is unclear. ...Cardiovascular magnetic resonance (CMR) diffusion tensor imaging (DTI) may provide a non-invasive, comprehensive, and geometrically accurate method to detect regional differences in structural remodeling in hypertrophy. We hypothesized that DTI parameters, such as fractional and planar anisotropy, would reflect myocyte remodeling due to pressure overload in a regionally-dependent manner.
We investigated the regional distributions of myocyte remodeling in rats with or without transverse aortic constriction (TAC) via direct measurement of myocyte dimensions with confocal imaging of thick tissue sections, and correlated myocyte cross-sectional area and other geometric features with parameters of diffusivity from ex-vivo DTI in the same regions of the same hearts.
We observed regional differences in several parameters from DTI between TAC hearts and SHAM controls. Consistent with previous studies, helix angles from DTI correlated strongly with those measured directly from histological sections (p < 0.001, R
= 0.71). There was a transmural gradient in myocyte cross-sectional area in SHAM hearts that was diminished in the TAC group. We also found several regions of significantly altered DTI parameters in TAC LV compared to SHAM, especially in myocyte sheet angle dispersion and planar anisotropy. Among others, these parameters correlated significantly with directly measured myocyte aspect ratios.
These results show that structural remodeling in pressure overload LV hypertrophy is regionally heterogeneous, especially transmurally, with a greater degree of remodeling in the sub-endocardium compared to the sub-epicardium. Additionally, several parameters derived from DTI correlated significantly with measurements of myocyte geometry from direct measurement in histological sections. We suggest that DTI may provide a non-invasive, comprehensive method to detect regional structural myocyte LV remodeling during disease.