Connecting Transcriptomics with Computational Modeling to Reveal Developmental Adaptations in the Human Pediatric Myocardium.

Background: Nearly 1% or 1.3 million babies are born with congenital heart disease (CHD) globally each year - many of whom will require palliative or corrective heart surgery within the first few years of life. A detailed understanding of cardiac maturation can help to expand our knowledge on cardiac diseases that develop during gestation, identify age-appropriate cardiovascular drug therapies, and inform clinical care decisions related to surgical repair, myocardial preservation, or postoperative management. Yet, to date, our knowledge of the temporal changes that cardiomyocytes undergo during postnatal development is largely limited to animal models. Methods: Right atrial tissue samples were collected from n=117 neonatal, infant, and pediatric patients undergoing correct surgery due to (acyanotic) CHD. Patients were stratified into five age groups: neonate (0-30 days), infant (31-364 days), toddler to preschool (1-5 years), school age (6-11 years), and adolescent to young adults (12-32 years). We measured age-dependent adaptations in cardiac gene expression, and used computational modeling to simulate action potential and calcium transients. Results: Enrichment of differentially expressed genes (DEG) was explored, revealing age-dependent changes in several key biological processes (cell cycle, cell division, mitosis), cardiac ion channels, and calcium handling genes. Gene-associated changes in ionic currents exhibited both linear trends and sudden shifts across developmental stages, with changes in calcium handling ( I NCX ) and repolarization ( I K1 ) most strongly associated with an age-dependent decrease in the action potential plateau potential and increase in triangulation, respectively. We also note a shift in repolarization reserve, with lower I Kr expression in younger patients, a finding likely tied to the increased amplitude of I Ks triggered by elevated sympathetic activation in pediatric patients. Conclusion: This study provides valuable insights into age-dependent changes in human cardiac gene expression and electrophysiology among patients with CHD, shedding light on molecular mechanisms underlying cardiac development and function across different developmental stages.

Cardiac-Specific Deletion of Scn8a Mitigates Dravet Syndrome-Associated Sudden Death in Adults.

BACKGROUND: Sudden unexpected death in epilepsy (SUDEP) is a fatal complication experienced by otherwise healthy epilepsy patients. Dravet syndrome (DS) is an inherited epileptic disorder resulting from loss of function of the voltage-gated sodium channel, NaV 1.1, and is associated with particularly high SUDEP risk. Evidence is mounting that NaVs abundant in the brain also occur in the heart, suggesting that the very molecular mechanisms underlying epilepsy could also precipitate cardiac arrhythmias and sudden death. Despite marked reduction of NaV 1.1 functional expression in DS, pathogenic late sodium current (INa,L) is paradoxically increased in DS hearts. However, the mechanisms by which DS directly impacts the heart to promote sudden death remain unclear. OBJECTIVES: In this study the authors sought to provide evidence implicating remodeling of Na+ - and Ca2+ -handling machinery, including NaV 1.6 and Na+/Ca2+exchanger (NCX) within transverse (T)-tubules in DS-associated arrhythmias. METHODS: The authors undertook scanning ion conductance microscopy (SICM)-guided patch clamp, super-resolution microscopy, confocal Ca2+ imaging, and in vivo electrocardiography studies in Scn1a haploinsufficient murine model of DS. RESULTS: DS promotes INa,L in T-tubular nanodomains, but not in other subcellular regions. Consistent with increased NaV activity in these regions, super-resolution microscopy revealed increased NaV 1.6 density near Ca2+release channels, the ryanodine receptors (RyR2) and NCX in DS relative to WT hearts. The resulting INa,L in these regions promoted aberrant Ca2+ release, leading to ventricular arrhythmias in vivo. Cardiac-specific deletion of NaV 1.6 protects adult DS mice from increased T-tubular late NaV activity and the resulting arrhythmias, as well as sudden death. CONCLUSIONS: These data demonstrate that NaV 1.6 undergoes remodeling within T-tubules of adult DS hearts serving as a substrate for Ca2+ -mediated cardiac arrhythmias and may be a druggable target for the prevention of SUDEP in adult DS subjects.

Indirect Correlative Light and Electron Microscopy (iCLEM): A Novel Pipeline for Multiscale Quantification of Structure From Molecules to Organs.

Correlative light and electron microscopy (CLEM) methods are powerful methods that combine molecular organization (from light microscopy) with ultrastructure (from electron microscopy). However, CLEM methods pose high cost/difficulty barriers to entry and have very low experimental throughput. Therefore, we have developed an indirect correlative light and electron microscopy (iCLEM) pipeline to sidestep the rate-limiting steps of CLEM (i.e., preparing and imaging the same samples on multiple microscopes) and correlate multiscale structural data gleaned from separate samples imaged using different modalities by exploiting biological structures identifiable by both light and electron microscopy as intrinsic fiducials. We demonstrate here an application of iCLEM, where we utilized gap junctions and mechanical junctions between muscle cells in the heart as intrinsic fiducials to correlate ultrastructural measurements from transmission electron microscopy (TEM), and focused ion beam scanning electron microscopy (FIB-SEM) with molecular organization from confocal microscopy and single molecule localization microscopy (SMLM). We further demonstrate how iCLEM can be integrated with computational modeling to discover structure-function relationships. Thus, we present iCLEM as a novel approach that complements existing CLEM methods and provides a generalizable framework that can be applied to any set of imaging modalities, provided suitable intrinsic fiducials can be identified.

Unraveling Impacts of Chamber-Specific Differences in Intercalated Disc Ultrastructure and Molecular Organization on Cardiac Conduction.

BACKGROUND: Propagation of action potentials through the heart coordinates the heartbeat. Thus, intercalated discs, specialized cell-cell contact sites that provide electrical and mechanical coupling between cardiomyocytes, are an important target for study. Impaired propagation leads to arrhythmias in many pathologies, where intercalated disc remodeling is a common finding, hence the importance and urgency of understanding propagation dependence on intercalated disc structure. Conventional modeling approaches cannot predict changes in propagation elicited by perturbations that alter intercalated disc ultrastructure or molecular organization, because of lack of quantitative structural data at subcellular through nano scales. OBJECTIVES: This study sought to quantify intercalated disc structure at these spatial scales in the healthy adult mouse heart and relate them to chamber-specific properties of propagation as a precursor to understanding the effects of pathological intercalated disc remodeling. METHODS: Using super-resolution light microscopy, electron microscopy, and computational image analysis, we provide here the first ever systematic, multiscale quantification of intercalated disc ultrastructure and molecular organization. RESULTS: By incorporating these data into a rule-based model of cardiac tissue with realistic intercalated disc structure, and comparing model predictions of electrical propagation with experimental measures of conduction velocity, we reveal that atrial intercalated discs can support faster conduction than their ventricular counterparts, which is normally masked by interchamber differences in myocyte geometry. Further, we identify key ultrastructural and molecular organization features underpinning the ability of atrial intercalated discs to support faster conduction. CONCLUSIONS: These data provide the first stepping stone to elucidating chamber-specific effects of pathological intercalated disc remodeling, as occurs in many arrhythmic diseases.

Emergent activity, heterogeneity, and robustness in a calcium feedback model of the sinoatrial node.

The sinoatrial node (SAN) is the primary pacemaker of the heart. SAN activity emerges at an early point in life and maintains a steady rhythm for the lifetime of the organism. The ion channel composition and currents of SAN cells can be influenced by a variety of factors. Therefore, the emergent activity and long-term stability imply some form of dynamical feedback control of SAN activity. We adapt a recent feedback model-previously utilized to describe control of ion conductances in neurons-to a model of SAN cells and tissue. The model describes a minimal regulatory mechanism of ion channel conductances via feedback between intracellular calcium and an intrinsic target calcium level. By coupling a SAN cell to the calcium feedback model, we show that spontaneous electrical activity emerges from quiescence and is maintained at steady state. In a 2D SAN tissue model, spatial variability in intracellular calcium targets lead to significant, self-organized heterogeneous ion channel expression and calcium transients throughout the tissue. Furthermore, multiple pacemaking regions appear, which interact and lead to time-varying cycle length, demonstrating that variability in heart rate is an emergent property of the feedback model. Finally, we demonstrate that the SAN tissue is robust to the silencing of leading cells or ion channel knockouts. Thus, the calcium feedback model can reproduce and explain many fundamental emergent properties of activity in the SAN that have been observed experimentally based on a minimal description of intracellular calcium and ion channel regulatory networks.

Unraveling Chamber-specific Differences in Intercalated Disc Ultrastructure and Molecular Organization and Their Impact on Cardiac Conduction.

During each heartbeat, the propagation of action potentials through the heart coordinates the contraction of billions of individual cardiomyocytes and is thus, a critical life process. Unsurprisingly, intercalated discs, which are cell-cell contact sites specialized to provide electrical and mechanical coupling between adjacent cardiomyocytes, have been the focus of much investigation. Slowed or disrupted propagation leads to potentially life-threatening arrhythmias in a wide range of pathologies, where intercalated disc remodeling is a common finding. Hence, the importance and urgency of understanding intercalated disc structure and its influence on action potential propagation. Surprisingly, however, conventional modeling approaches cannot predict changes in propagation elicited by perturbations that alter intercalated disc ultrastructure or molecular organization, owing to lack of quantitative structural data at subcellular through nano scales. In order to address this critical gap in knowledge, we sought to quantify intercalated disc structure at these finer spatial scales in the healthy adult mouse heart and relate them to function in a chamber-specific manner as a precursor to understanding the impacts of pathological intercalated disc remodeling. Using super-resolution light microscopy, electron microscopy, and computational image analysis, we provide here the first ever systematic, multiscale quantification of intercalated disc ultrastructure and molecular organization. By incorporating these data into a rule-based model of cardiac tissue with realistic intercalated disc structure, and comparing model predictions of electrical propagation with experimental measures of conduction velocity, we reveal that atrial intercalated discs can support faster conduction than their ventricular counterparts, which is normally masked by inter-chamber differences in myocyte geometry. Further, we identify key ultrastructural and molecular organization features underpinning the ability of atrial intercalated discs to support faster conduction. These data provide the first stepping stone to elucidating chamber-specific impacts of pathological intercalated disc remodeling, as occurs in many arrhythmic diseases.

Modeling incomplete penetrance in long QT syndrome type 3 through ion channel heterogeneity: an in silico population study.

Many cardiac diseases are characterized by an increased late sodium current, including heart failure, hypertrophic cardiomyopathy, and inherited long QT syndrome type 3 (LQT3). The late sodium current in LQT3 is caused by a gain-of-function mutation in the voltage-gated sodium channel Nav1.5. Despite a well-defined genetic cause of LQT3, treatment remains inconsistent because of incomplete penetrance of the mutation and variability of antiarrhythmic efficacy. Here, we investigate the relationship between LQT3-associated mutation incomplete penetrance and variability in ion channel expression, simulating a population of 1,000 individuals with the O'Hara-Rudy model of the human ventricular myocyte. We first simulate healthy electrical activity (i.e., in the absence of a mutation) and then incorporate heterozygous expression for three LQT3-associated mutations (Y1795C, I1768V, and ΔKPQ), to directly compare the effects of each mutation on individuals across a diverse population. For all mutations, we find that susceptibility, defined by either the presence of an early afterdepolarization (EAD) or prolonged action potential duration (APD), primarily depends on the balance between the conductance of IKr and INa, for which individuals with a higher IKr-to-INa ratio are less susceptible. Furthermore, we find distinct differences across the population, observing individuals susceptible to zero, one, two, or all three mutations. Individuals tend to be less susceptible with an appropriate balance of repolarizing currents, typically via increased IKs or IK1. Interestingly, the more critical repolarizing current is mutation specific. We conclude that the balance between key currents plays a significant role in mutant-specific presentation of the disease phenotype in LQT3.NEW & NOTEWORTHY An in silico population approach investigates the relationship between variability in ion channel expression and gain-of-function mutations in the voltage-gated sodium channel associated with the congenital disorder long QT syndrome type 3 (LQT3). We find that ion channel variability can contribute to incomplete penetrance of the mutation, with mutant-specific differences in ion channel conductances leading to susceptibility to proarrhythmic action potential duration prolongation or early afterdepolarizations.

A mathematical model of immunomodulatory treatment in myocardial infarction.

A heart attack, or acute myocardial infarction (MI) is caused by the acute occlusion of a coronary artery. MI is associated with 30% mortality; approximately half of the deaths occur prior to arrival at the hospital. Reperfusion therapy in the hospital is a medical treatment to restore blood flow through the blocked artery; treatment includes drugs and surgery. However, the damage to the heart muscles through the infarct area is permanent and there is additional damage around the infarct area due to inflammation or insufficient oxygen supply. Approximately half of the patients who survive MI are hospitalized again within one year after reperfusion treatment. In this paper we develop a mathematical model of MI and use it to assess the efficacy of drugs used, post reperfusion, to reduce the damage caused by inflammation in a region of the left ventricular wall surrounding the infarct area. The mathematical model, represented by a system of partial differential equations. The model variables include myocytes, endothelial cells, neutrophils, macrophages, fibroblasts and cytokines that play a role in the interactions among these cells. The drugs used to in the model include IL-1, TNF-α and TGF-ß inhibitors, and the delivery of VEGF. The model is based on mice data. In particular, we find that immunomodulatory treatment with TNF-α and IL-1 inhibitors can significantly increase the low density of myocytes bordering the infarct area by 50-60% and decrease the abnormally high density of ECM in a region surrounding the infarct area.

Distributed synthesis of sarcolemmal and sarcoplasmic reticulum membrane proteins in cardiac myocytes.

It is widely assumed that synthesis of membrane proteins, particularly in the heart, follows the classical secretory pathway with mRNA translation occurring in perinuclear regions followed by protein trafficking to sites of deployment. However, this view is based on studies conducted in less-specialized cells, and has not been experimentally addressed in cardiac myocytes. Therefore, we undertook direct experimental investigation of protein synthesis in cardiac tissue and isolated myocytes using single-molecule visualization techniques and a novel proximity-ligated in situ hybridization approach for visualizing ribosome-associated mRNA molecules for a specific protein species, indicative of translation sites. We identify here, for the first time, that the molecular machinery for membrane protein synthesis occurs throughout the cardiac myocyte, and enables distributed synthesis of membrane proteins within sub-cellular niches where the synthesized protein functions using local mRNA pools trafficked, in part, by microtubules. We also observed cell-wide distribution of membrane protein mRNA in myocardial tissue from both non-failing and hypertrophied (failing) human hearts, demonstrating an evolutionarily conserved distributed mechanism from mouse to human. Our results identify previously unanticipated aspects of local control of cardiac myocyte biology and highlight local protein synthesis in cardiac myocytes as an important potential determinant of the heart's biology in health and disease.

Intercalated disk nanoscale structure regulates cardiac conduction.

The intercalated disk (ID) is a specialized subcellular region that provides electrical and mechanical connections between myocytes in the heart. The ID has a clearly defined passive role in cardiac tissue, transmitting mechanical forces and electrical currents between cells. Recent studies have shown that Na+ channels, the primary current responsible for cardiac excitation, are preferentially localized at the ID, particularly within nanodomains such as the gap junction-adjacent perinexus and mechanical junction-associated adhesion-excitability nodes, and that perturbations of ID structure alter cardiac conduction. This suggests that the ID may play an important, active role in regulating conduction. However, the structures of the ID and intercellular cleft are not well characterized and, to date, no models have incorporated the influence of ID structure on conduction in cardiac tissue. In this study, we developed an approach to generate realistic finite element model (FEM) meshes replicating nanoscale of the ID structure, based on experimental measurements from transmission electron microscopy images. We then integrated measurements of the intercellular cleft electrical conductivity, derived from the FEM meshes, into a novel cardiac tissue model formulation. FEM-based calculations predict that the distribution of cleft conductances is sensitive to regional changes in ID structure, specifically the intermembrane separation and gap junction distribution. Tissue-scale simulations predict that ID structural heterogeneity leads to significant spatial variation in electrical polarization within the intercellular cleft. Importantly, we found that this heterogeneous cleft polarization regulates conduction by desynchronizing the activation of postjunctional Na+ currents. Additionally, these heterogeneities lead to a weaker dependence of conduction velocity on gap junctional coupling, compared with prior modeling formulations that neglect or simplify ID structure. Further, we found that disruption of local ID nanodomains can either slow or enhance conduction, depending on gap junctional coupling strength. Our study therefore suggests that ID nanoscale structure can play a significant role in regulating cardiac conduction.

A mathematical model of the multiple sclerosis plaque.

Multiple sclerosis is an autoimmune disease that affects white matter in the central nervous system. It is one of the primary causes of neurological disability among young people. Its characteristic pathological lesion is called a plaque, a zone of inflammatory activity and tissue destruction that expands radially outward by destroying the myelin and oligodendrocytes of white matter. The present paper develops a mathematical model of the multiple sclerosis plaques. Although these plaques do not provide reliable information of the clinical disability in MS, they are nevertheless useful as a primary outcome measure of Phase II trials. The model consists of a system of partial differential equations in a simplified geometry of the lesion, consisting of three domains: perivascular space, demyelinated plaque, and white matter. The model describes the activity of various pro- and anti-inflammatory cells and cytokines in the plaque, and quantifies their effect on plaque growth. We show that volume growth of plaques are in qualitative agreement with reported clinical studies of several currently used drugs. We then use the model to explore treatments with combinations of such drugs, and with experimental drugs. We finally consider the benefits of early vs. delayed treatment.

Rheumatoid arthritis - a mathematical model.

Rheumatoid arthritis (RA) is a common autoimmune disease that mainly affects the joints. It is characterized by synovial inflammation, which may result in cartilage and bone destruction. The present paper develops a mathematical model of chronic RA. The model is represented by a system of partial differential equations (PDEs) in the synovial fluid, the synovial membrane, and the cartilage. The model characterizes the progression of the disease in terms of the degradation of the cartilage. More precisely, we assume a simplified geometry in which the synovial membrane and the cartilage are planar layers adjacent to each other. We then quantify the state of the disease by how much the cartilage layer has decreased, or, equivalently, how much the synovial layer has increased. The model is used to evaluate treatments of RA by currently used drugs, as well as by experimental drugs.

Correlation between lifetime heterogeneity and kinetics heterogeneity during chlorophyll fluorescence induction in leaves: 2. Multi-frequency phase and modulation analysis evidences a loosely connected PSII pigment-protein complex.

We report the first direct decomposition of the fluorescence lifetime heterogeneity during multiphasic fluorescence induction in dark-adapted leaves by multi-frequency phase and modulation fluorometry (PMF). A very fast component, assigned to photosystem I (PSI), was found to be constant in lifetime and yield, whereas the two slow components, which are strongly affected by the closure of the reaction centers by light, were assigned to PSII. Based on a modified "reversible radical pair" kinetic model with three compartments, we showed that a loosely connected pigment complex, which is assumed to be the CP47 complex, plays a specific role with respect to the structure and function of the PSII: (i) it explains the heterogeneity of PSII fluorescence lifetime as a compartmentation of excitation energy in the antenna, (ii) it is the site of a conformational change in the first second of illumination, and (iii) it is involved in the mechanisms of nonphotochemical quenching (NPQ). On the basis of the multi-frequency PMF analysis, we reconciled two apparently antagonistic aspects of chlorophyll a fluorescence in vivo: it is heterogeneous with respect to the kinetic structure (several lifetime components) and homogeneous with respect to average quantities (quasi-linear mean tau-Phi relationship).

Correlation between lifetime heterogeneity and kinetics heterogeneity during chlorophyll fluorescence induction in leaves: 1. Mono-frequency phase and modulation analysis reveals a conformational change of a PSII pigment complex during the IP thermal phase.

The relationship between the fluorescence lifetime (tau) and yield (Phi) obtained in phase and modulation fluorometry at 54 MHz during the chlorophyll fluorescence induction in dark-adapted leaves under low actinic light has been investigated. Three typical phases have been identified: (i) linear during the OI photochemical rise, (ii) convex curvature during the subsequent IP thermal rise, and (iii) linear during the PS slow decay. A similar relationship has been obtained in the fluorescence induction for the fluorescence yield measured at 685 nm plotted versus the fluorescence yield measured at 735 nm. A spectrally resolved analysis shows that the curvature of the tau-Phi relationship is not due to chlorophyll fluorescence reabsorption effects. Several other hypotheses are discussed and we conclude that the curvature of the tau-Phi relationship is due to a variable and transitory nonphotochemical quenching. We tentatively propose that this quenching results from a conformational change of a pigment-protein complex of Photosystem II core antenna during the IP phase and could explain both spectral and temporal transitory changes of the fluorescence. A variable blue shift of the 685 nm peak of the fluorescence spectrum during the IP phase has been observed, supporting this hypothesis.

Early light-induced proteins protect Arabidopsis from photooxidative stress.

The early light-induced proteins (ELIPs) belong to the multigenic family of light-harvesting complexes, which bind chlorophyll and absorb solar energy in green plants. ELIPs accumulate transiently in plants exposed to high light intensities. By using an Arabidopsis thaliana mutant (chaos) affected in the posttranslational targeting of light-harvesting complex-type proteins to the thylakoids, we succeeded in suppressing the rapid accumulation of ELIPs during high-light stress, resulting in leaf bleaching and extensive photooxidative damage. Constitutive expression of ELIP genes in chaos before light stress resulted in ELIP accumulation and restored the phototolerance of the plants to the wild-type level. Free chlorophyll, a generator of singlet oxygen in the light, was detected by chlorophyll fluorescence lifetime measurements in chaos leaves before the symptoms of oxidative stress appeared. Our findings indicate that ELIPs fulfill a photoprotective function that could involve either the binding of chlorophylls released during turnover of pigment-binding proteins or the stabilization of the proper assembly of those proteins during high-light stress.