A novel mutant cardiac troponin C disrupts molecular motions critical for calcium binding affinity and cardiomyocyte contractility

Troponin C gene mutations on cardiac muscle cell and skeletal Regulation: A comprehensive review

BACKGROUND

The troponin complex plays a crucial role in regulating skeletal and cardiac contraction. Congenital myopathies can occur due to several mutations in genes that encode skeletal troponin. Moreover, there is limited information regarding the composition of skeletal troponin. This review specifically examines a comprehensive review of the TNNC gene mutations on cardiac and skeletal regulations.

MAIN BODY

Troponin C (TNNC) has been linked to a newly discovered inherited muscle disorder. Genetic variations in genes that encode skeletal troponin can impair the function of sarcomeres. Various treatment approaches have been employed to mitigate the impact of variations, including the use of troponin activators, the injection of wild-type protein via AAV gene therapy, and myosin modification to enhance muscle contraction. The processes responsible for the pathophysiological implications of the variations in genes that encode skeletal troponin are not fully understood.

CONCLUSION

This comprehensive review will contribute to the understanding of the relationship between human cardiomyopathy and TNNC mutations and will guide the development of therapy approaches.

Are fetal microchimerism and circulating fetal extracellular vesicles important links between spontaneous preterm delivery and maternal cardiovascular disease risk?

Trafficking and persistence of fetal microchimeric cells (fMCs) and circulating extracellular vesicles (EVs) have been observed in animals and humans, but their consequences in the maternal body and their mechanistic contributions to maternal physiology and pathophysiology are not yet fully defined. Fetal cells and EVs may help remodel maternal organs after pregnancy-associated changes, but the cell types and EV cargos reaching the mother in preterm pregnancies after exposure to various risk factors can be distinct from term pregnancies. As preterm delivery-associated maternal complications are rising, revisiting this topic and formulating scientific questions for future research to reduce the risk of maternal morbidities are timely. Epidemiological studies report maternal cardiovascular risk as one of the major complications after preterm delivery. This paper suggests a potential link between fMCs and circulating EVs and adverse maternal cardiovascular outcomes post-pregnancies, the underlying mechanisms, consequences, and methods for and how this link might be assessed.

MicroBundleCompute: Automated segmentation, tracking, and analysis of subdomain deformation in cardiac microbundles

Advancing human induced pluripotent stem cell derived cardiomyocyte (hiPSC-CM) technology will lead to significant progress ranging from disease modeling, to drug discovery, to regenerative tissue engineering. Yet, alongside these potential opportunities comes a critical challenge: attaining mature hiPSC-CM tissues. At present, there are multiple techniques to promote maturity of hiPSC-CMs including physical platforms and cell culture protocols. However, when it comes to making quantitative comparisons of functional behavior, there are limited options for reliably and reproducibly computing functional metrics that are suitable for direct cross-system comparison. In addition, the current standard functional metrics obtained from time-lapse images of cardiac microbundle contraction reported in the field (i.e., post forces, average tissue stress) do not take full advantage of the available information present in these data (i.e., full-field tissue displacements and strains). Thus, we present "MicroBundleCompute," a computational framework for automatic quantification of morphology-based mechanical metrics from movies of cardiac microbundles. Briefly, this computational framework offers tools for automatic tissue segmentation, tracking, and analysis of brightfield and phase contrast movies of beating cardiac microbundles. It is straightforward to implement, runs without user intervention, requires minimal input parameter setting selection, and is computationally inexpensive. In this paper, we describe the methods underlying this computational framework, show the results of our extensive validation studies, and demonstrate the utility of exploring heterogeneous tissue deformations and strains as functional metrics. With this manuscript, we disseminate "MicroBundleCompute" as an open-source computational tool with the aim of making automated quantitative analysis of beating cardiac microbundles more accessible to the community.

Altering Calcium Sensitivity in Heart Failure: A Crossroads of Disease Etiology and Therapeutic Innovation

Heart failure (HF) presents a significant clinical challenge, with current treatments mainly easing symptoms without stopping disease progression. The targeting of calcium (Ca2+) regulation is emerging as a key area for innovative HF treatments that could significantly alter disease outcomes and enhance cardiac function. In this review, we aim to explore the implications of altered Ca2+ sensitivity, a key determinant of cardiac muscle force, in HF, including its roles during systole and diastole and its association with different HF types-HF with preserved and reduced ejection fraction (HFpEF and HFrEF, respectively). We further highlight the role of the two rate constants kon (Ca2+ binding to Troponin C) and koff (its dissociation) to fully comprehend how changes in Ca2+ sensitivity impact heart function. Additionally, we examine how increased Ca2+ sensitivity, while boosting systolic function, also presents diastolic risks, potentially leading to arrhythmias and sudden cardiac death. This suggests that strategies aimed at moderating myofilament Ca2+ sensitivity could revolutionize anti-arrhythmic approaches, reshaping the HF treatment landscape. In conclusion, we emphasize the need for precision in therapeutic approaches targeting Ca2+ sensitivity and call for comprehensive research into the complex interactions between Ca2+ regulation, myofilament sensitivity, and their clinical manifestations in HF.

Mouse Models of Cardiomyopathies Caused by Mutations in Troponin C

Cardiac muscle contraction is regulated via Ca2+ exchange with the hetero-trimeric troponin complex located on the thin filament. Binding of Ca2+ to cardiac troponin C, a Ca2+ sensing subunit within the troponin complex, results in a series of conformational re-arrangements among the thin filament components, leading to an increase in the formation of actomyosin cross-bridges and muscle contraction. Ultimately, a decline in intracellular Ca2+ leads to the dissociation of Ca2+ from troponin C, inhibiting cross-bridge cycling and initiating muscle relaxation. Therefore, troponin C plays a crucial role in the regulation of cardiac muscle contraction and relaxation. Naturally occurring and engineered mutations in troponin C can lead to altered interactions among components of the thin filament and to aberrant Ca2+ binding and exchange with the thin filament. Mutations in troponin C have been associated with various forms of cardiac disease, including hypertrophic, restrictive, dilated, and left ventricular noncompaction cardiomyopathies. Despite progress made to date, more information from human studies, biophysical characterizations, and animal models is required for a clearer understanding of disease drivers that lead to cardiomyopathies. The unique use of engineered cardiac troponin C with the L48Q mutation that had been thoroughly characterized and genetically introduced into mouse myocardium clearly demonstrates that Ca2+ sensitization in and of itself should not necessarily be considered a disease driver. This opens the door for small molecule and protein engineering strategies to help boost impaired systolic function. On the other hand, the engineered troponin C mutants (I61Q and D73N), genetically introduced into mouse myocardium, demonstrate that Ca2+ desensitization under basal conditions may be a driving factor for dilated cardiomyopathy. In addition to enhancing our knowledge of molecular mechanisms that trigger hypertrophy, dilation, morbidity, and mortality, these cardiomyopathy mouse models could be used to test novel treatment strategies for cardiovascular diseases. In this review, we will discuss (1) the various ways mutations in cardiac troponin C might lead to disease; (2) relevant data on mutations in cardiac troponin C linked to human disease, and (3) all currently existing mouse models containing cardiac troponin C mutations (disease-associated and engineered).

Umbrella Sampling Simulations of Cardiac Thin Filament Reveal Thermodynamic Consequences of Troponin I Inhibitory Peptide Mutations

The cardiac thin filament comprises F-actin, tropomyosin, and troponin (cTn). cTn is composed of three subunits: troponin C (cTnC), troponin I (cTnI), and troponin T (cTnT). To computationally study the effect of the thin filament on cTn activation events, we employed targeted molecular dynamics followed by umbrella sampling using a model of the thin filament to measure the thermodynamics of cTn transition events. Our simulations revealed that the thin filament causes an increase in the free energy required to open the cTnC hydrophobic patch and causes a more favorable interaction between this region and the cTnI switch peptide. Mutations to the cTn complex can lead to cardiomyopathy, a collection of diseases that present clinically with symptoms of hypertrophy or dilation of the cardiac muscle, leading to impairment of the heart's ability to function normally and ultimately myocardial infarction or heart failure. Upon introduction of cardiomyopathic mutations to R145 of cTnI, we observed a general decrease in the free energy of opening the cTnC hydrophobic patch, which is on par with previous experimental results. These mutations also exhibited a decrease in electrostatic interactions between cTnI-R145 and actin-E334. After introduction of a small molecule to the wild-type cTnI-actin interface to intentionally disrupt intersubunit contacts, we successfully observed similar thermodynamic consequences and disruptions to the same protein-protein contacts as observed with the cardiomyopathic mutations. Computational studies utilizing the cTn complex in isolation would have been unable to observe these effects, highlighting the importance of using a more physiologically relevant thin-filament model to investigate the global consequences of cardiomyopathic mutations to the cTn complex.

Computational Exploration and Characterization of Potential Calcium Sensitizing Mutations in Cardiac Troponin C

Calcium-dependent heart muscle contraction is regulated by the cardiac troponin protein complex (cTn) and specifically by the N-terminal domain of its calcium binding subunit (cNTnC). cNTnC contains one calcium binding site (site II), and altered calcium binding in this site has been studied for decades. It has been previously shown that cNTnC mutants, which increase calcium sensitization may have therapeutic benefits, such as restoring cardiac muscle contractility and functionality post-myocardial infarction events. Here, we computationally characterized eight mutations for their potential effects on calcium binding affinity in site II of cNTnC. We utilized two distinct methods to estimate calcium binding: adaptive steered molecular dynamics (ASMD) and thermodynamic integration (TI). We observed a sensitizing trend for all mutations based on the employed ASMD methodology. The TI results showed excellent agreement with experimentally known calcium binding affinities in wild-type cNTnC. Based on the TI results, five mutants were predicted to increase calcium sensitivity in site II. This study presents an interesting comparison of the two computational methods, which have both been shown to be valuable tools in characterizing the impacts of calcium sensitivity in mutant cNTnC systems.

Lower troponin expression in the right ventricle of rats explains interventricular differences in E-C coupling

Despite distinctive functional and anatomic differences, a precise understanding of the cardiac interventricular differences in excitation–contraction (E–C) coupling mechanisms is still lacking. Here, we directly compared rat right and left cardiomyocytes (RVCM and LVCM). Whole-cell patch clamp, the IonOptix system, and fura-2 fluorimetry were used to measure electrical properties (action potential and ionic currents), single-cell contractility, and cytosolic Ca²⁺ ([Ca²⁺]_i), respectively. Myofilament proteins were analyzed by immunoblotting. RVCM showed significantly shorter action potential duration (APD) and higher density of transient outward K⁺ current (I_to). However, the triggered [Ca²⁺]_i change (Ca²⁺ transient) was not different, while the decay rate of the Ca²⁺ transient was slower in RVCM. Although the relaxation speed was also slower, the sarcomere shortening amplitude (ΔSL) was smaller in RVCM. SERCA activity was ∼60% lower in RVCM, which is partly responsible for the slower decay of the Ca²⁺ transient. Immunoblot analysis revealed lower expression of the cardiac troponin complex (cTn) in RVCM, implying a smaller Ca²⁺ buffering capacity (κ_S), which was proved by in situ analysis. The introduction of these new levels of cTn, I_to, and SERCA into a mathematical model of rat LVCM reproduced the similar Ca²⁺ transient, slower Ca²⁺ decay, shorter APD, and smaller ΔSL of RVCM. Taken together, these data show reduced expression of cTn proteins in the RVCM, which provides an explanation for the interventricular difference in the E–C coupling kinetics.

Microengineered platforms for characterizing the contractile function of in vitro cardiac models

Emerging heart-on-a-chip platforms are promising approaches to establish cardiac cell/tissue models in vitro for research on cardiac physiology, disease modeling and drug cardiotoxicity as well as for therapeutic discovery. Challenges still exist in obtaining the complete capability of in situ sensing to fully evaluate the complex functional properties of cardiac cell/tissue models. Changes to contractile strength (contractility) and beating regularity (rhythm) are particularly important to generate accurate, predictive models. Developing new platforms and technologies to assess the contractile functions of in vitro cardiac models is essential to provide information on cell/tissue physiologies, drug-induced inotropic responses, and the mechanisms of cardiac diseases. In this review, we discuss recent advances in biosensing platforms for the measurement of contractile functions of in vitro cardiac models, including single cardiomyocytes, 2D monolayers of cardiomyocytes, and 3D cardiac tissues. The characteristics and performance of current platforms are reviewed in terms of sensing principles, measured parameters, performance, cell sources, cell/tissue model configurations, advantages, and limitations. In addition, we highlight applications of these platforms and relevant discoveries in fundamental investigations, drug testing, and disease modeling. Furthermore, challenges and future outlooks of heart-on-a-chip platforms for in vitro measurement of cardiac functional properties are discussed.

A micromachined force sensing apparatus and method for human engineered cardiac tissue and induced pluripotent stem cell characterization

Induced pluripotent stem cell derived-cardiomyocytes (iPSC-CMs) have great potential for cell therapy, drug assessment, and for understanding the pathophysiology and genetic underpinnings of cardiac diseases. Contraction forces are one of the most important characteristics of cardiac function and are predictors of healthy and diseased states. Cantilever techniques, such as atomic force microscopy, measure the vertical force of a single cell, while systems designed to more closely resemble the physical heart function, such as engineered cardiac tissue held by end-posts, measure the axial force. One important question is how do these two force measurements correlate? By establishing a correlation of the axial and vertical force, we will be one step closer in being able to use single cell iPSC-CMs as models. A novel micromachined sensor for measuring force contractions of engineered tissue has been developed. Using this novel sensor, a correlation between axial force and vertical force is experimentally established. This finding supports the use of vertical measurements as an alternative to tissue measurements.

Cardiomyocyte Dysfunction in Inherited Cardiomyopathies

Inherited cardiomyopathies form a heterogenous group of disorders that affect the structure and function of the heart. Defects in the genes encoding sarcomeric proteins are associated with various perturbations that induce contractile dysfunction and promote disease development. In this review we aimed to outline the functional consequences of the major inherited cardiomyopathies in terms of myocardial contraction and kinetics, and to highlight the structural and functional alterations in some sarcomeric variants that have been demonstrated to be involved in the pathogenesis of the inherited cardiomyopathies. A particular focus was made on mutation-induced alterations in cardiomyocyte mechanics. Since no disease-specific treatments for familial cardiomyopathies exist, several novel agents have been developed to modulate sarcomere contractility. Understanding the molecular basis of the disease opens new avenues for the development of new therapies. Furthermore, the earlier the awareness of the genetic defect, the better the clinical prognostication would be for patients and the better the prevention of development of the disease.

Targeted Drug Delivery for Cardiovascular Disease: Modeling of Modulated Extracellular Vesicle Release Rates

Invasive and medical therapy has led to major improvements in cardiovascular disease management, but important challenges remain open. The discovery of a nano-sized system of extracellular vesicles (EVs) is opening new possibilities for reprogramming malfunctioning cells and indicates that EVs can be employed in therapeutic biomedical applications as engineered drug vehicles. Molecular communication (MC) has applications for treating cells with directed drug delivery, employing special targeting transmembrane proteins. In this paper, we propose a novel drug delivery system for cardiovascular diseases using an EV-mediated MC platform and exemplify the potential use in hypertrophic cardiomyopathy. We utilize intracellular calcium signaling as a natural mediator of EVs released from synthetic cells and model the release rate. We propose to use the cells as a therapeutic release system with a control signal input which modulates the EVs release rate as the output signal. We also study the frequency domain of the proposed model and estimate the transfer function of the therapeutic release system model numerically where the root-mean-square error for two separate estimated output signals are 0.0353 and 0.0124. The proposed EV-mediated targeted drug delivery system can make breakthroughs in future healthcare, in cardiovascular and other diseases where targeting is required.

Calcium and Heart Failure: How Did We Get Here and Where Are We Going?

The occurrence and prevalence of heart failure remain high in the United States as well as globally. One person dies every 30 s from heart disease. Recognizing the importance of heart failure, clinicians and scientists have sought better therapeutic strategies and even cures for end-stage heart failure. This exploration has resulted in many failed clinical trials testing novel classes of pharmaceutical drugs and even gene therapy. As a result, along the way, there have been paradigm shifts toward and away from differing therapeutic approaches. The continued prevalence of death from heart failure, however, clearly demonstrates that the heart is not simply a pump and instead forces us to consider the complexity of simplicity in the pathophysiology of heart failure and reinforces the need to discover new therapeutic approaches.

Structural Changes beyond the EF-Hand Contribute to Apparent Calcium Binding Affinities: Insights from Parvalbumins

Members of the parvalbumin (PV) family of calcium (Ca²⁺) binding proteins (CBPs) share a relatively high level of sequence similarity. However, their Ca²⁺ affinities and selectivities against competing ions like Mg²⁺ can widely vary. We conducted molecular dynamics simulations of several α-parvalbumin (αPV) constructs with micromolar to nanomolar Ca²⁺ affinities to identify structural and dynamic features that contribute to their binding of ions. Specifically, we examined a D94S/G98E construct with a lower Ca²⁺ affinity (≈-18 kcal/mol) relative to the wild type (WT) (≈-22 kcal/mol) and an S55D/E59D variant with enhanced affinity (≈-24 kcal/mol). Additionally, we also examined the binding of Mg²⁺ to these isoforms, which is much weaker than Ca²⁺. We used mean spherical approximation (MSA) theory to evaluate ion binding thermodynamics within the proteins' EF-hand domains to account for the impact of ions' finite sizes and the surrounding electrolyte composition. While the MSA scores differentiated Mg²⁺ from Ca²⁺, they did not indicate that Ca²⁺ binding affinities at the binding loop differed between the PV isoforms. Instead, molecular mechanics generalized Born surface area (MM/GBSA) approximation energies, which we used to quantify the thermodynamic cost of structural rearrangement of the proteins upon binding ions, indicated that S55D/E59D αPV favored Ca²⁺ binding by -20 kcal/mol relative to WT versus 30 kcal/mol for D94S/G98E αPV. Meanwhile, Mg²⁺ binding was favored for the S55D/E59D αPV and D94S/G98E αPV variants by -18.32 and -1.65 kcal/mol, respectively. These energies implicate significant contributions to ion binding beyond oxygen coordination at the binding loop, which stemmed from changes in α-helicity, β-sheet character, and hydrogen bonding. Hence, Ca²⁺ affinity and selectivity against Mg²⁺ are emergent properties stemming from both local effects within the proteins' ion binding sites as well as non-local contributions elsewhere. Our findings broaden our understanding of the molecular bases governing αPV ion binding that are likely shared by members of the broad family of CBPs.

Adaptative Steered Molecular Dynamics Study of Mutagenesis Effects on Calcium Affinity in the Regulatory Domain of Cardiac Troponin C

Calcium-dependent cardiac muscle contraction is regulated by the protein complex troponin (cTn) and specifically by the regulatory N-terminal domain (N-cTnC) which contains one active Ca²⁺ binding site (site II). It has been previously shown that cardiac muscle contractility and functionality is affected by mutations in N-cTnC which alter calcium binding affinity. Here, we describe the application of adaptive steered molecular dynamics to characterize the influence of N-cTnC mutations on site II calcium binding affinity. We observed the correct trends for all of the studied calcium sensitizing and desensitizing mutants, in conjunction with loop II perturbations. Additionally, the potential of mean force accuracy was shown to increase substantially with increasingly slower speeds and using fewer trajectories. This study presents a novel approach to computationally estimate the Ca²⁺ binding affinity of N-cTnC structures and is a valuable potential tool to support the design and characterization of novel mutations with potential therapeutic benefits.

Dilated Cardiomyopathy Mutations and Phosphorylation disrupt the Active Orientation of Cardiac Troponin C

Cardiac troponin (cTn) is made up of three subunits, cTnC, cTnI, and cTnT. The regulatory N-terminal domain of cTnC (cNTnC) controls cardiac muscle contraction in a calcium-dependent manner. We show that calcium-saturated cNTnC can adopt two different orientations, with the "active" orientation consistent with the 2020 cryo-EM structure of the activated cardiac thin filament by Yamada et al. Using solution NMR ¹⁵N R₂ relaxation analysis, we demonstrate that the two domains of cTnC tumble independently (average R₂ 10 s^-1), being connected by a flexible linker. However, upon addition of cTnI_1-77, the complex tumbles as a rigid unit (R₂ 30 s^-1). cTnI phosphomimetic mutants S22D/S23D, S41D/S43D and dilated cardiomyopathy- (DCM-)associated mutations cTnI K35Q, cTnC D75Y, and cTnC G159D destabilize the active orientation of cNTnC, with intermediate ¹⁵N R₂ rates (R₂ 17-23 s^-1). The active orientation of cNTnC is stabilized by the flexible tails of cTnI, cTnI_1-37 and cTnI_135-209. Surprisingly, when cTnC is incorporated into complexes lacking these tails (cTnC-cTnI_38-134, cTnC-cTnT_223-288, or cTnC-cTnI_38-134-cTnT_223-288), the cNTnC domain is still immobilized, revealing a new interaction between cNTnC and the IT-arm that stabilizes a "dormant" orientation. We propose that the calcium sensitivity of the cardiac troponin complex is regulated by an equilibrium between active and dormant orientations, which can be shifted through post-translational modifications or DCM-associated mutations.

A comprehensive guide to genetic variants and post-translational modifications of cardiac troponin C

Familial cardiomyopathy is an inherited disease that affects the structure and function of heart muscle and has an extreme range of phenotypes. Among the millions of affected individuals, patients with hypertrophic (HCM), dilated (DCM), or left ventricular non-compaction (LVNC) cardiomyopathy can experience morphologic changes of the heart which lead to sudden death in the most detrimental cases. TNNC1, the gene that codes for cardiac troponin C (cTnC), is a sarcomere gene associated with cardiomyopathies in which probands exhibit young age of presentation and high death, transplant or ventricular fibrillation events relative to TNNT2 and TNNI3 probands. Using GnomAD, ClinVar, UniProt and PhosphoSitePlus databases and published literature, an extensive list to date of identified genetic variants in TNNC1 and post-translational modifications (PTMs) in cTnC was compiled. Additionally, a recent cryo-EM structure of the cardiac thin filament regulatory unit was used to localize each functionally studied amino acid variant and each PTM (acetylation, glycation, s-nitrosylation, phosphorylation) in the structure of cTnC. TNNC1 has a large number of variants (> 100) relative to other genes of the same transcript size. Surprisingly, the mapped variant amino acids and PTMs are distributed throughout the cTnC structure. While many cardiomyopathy-associated variants are localized in α-helical regions of cTnC, this was not statistically significant χ² (p = 0.72). Exploring the variants in TNNC1 and PTMs of cTnC in the contexts of cardiomyopathy association, physiological modulation and potential non-canonical roles provides insights into the normal function of cTnC along with the many facets of TNNC1 as a cardiomyopathic gene.

Increased myocyte calcium sensitivity in end-stage pediatric dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is the most common cause of heart failure (HF) in children, resulting in high mortality and need for heart transplantation. The pathophysiology underlying pediatric DCM is largely unclear; however, there is emerging evidence that molecular adaptations and response to conventional HF medications differ between children and adults. To gain insight into alterations leading to systolic dysfunction in pediatric DCM, we measured cardiomyocyte contractile properties and sarcomeric protein phosphorylation in explanted pediatric DCM myocardium (N = 8 subjects) compared with nonfailing (NF) pediatric hearts (N = 8 subjects). Force-pCa curves were generated from skinned cardiomyocytes in the presence and absence of protein kinase A. Sarcomeric protein phosphorylation was quantified with Pro-Q Diamond staining after gel electrophoresis. Pediatric DCM cardiomyocytes demonstrate increased calcium sensitivity (pCa₅₀ =5.70 ± 0.0291), with an associated decrease in troponin (Tn)I phosphorylation compared with NF pediatric cardiomyocytes (pCa₅₀ =5.59 ± 0.0271, P = 0.0073). Myosin binding protein C and TnT phosphorylation are also lower in pediatric DCM, whereas desmin phosphorylation is increased. Pediatric DCM cardiomyocytes generate peak tension comparable to that of NF pediatric cardiomyocytes [DCM 29.7 mN/mm², interquartile range (IQR) 21.5-49.2 vs. NF 32.8 mN/mm², IQR 21.5-49.2 mN/mm²; P = 0.6125]. In addition, cooperativity is decreased in pediatric DCM compared with pediatric NF (Hill coefficient: DCM 1.56, IQR 1.31-1.94 vs. NF 1.94, IQR 1.36-2.86; P = 0.0425). Alterations in sarcomeric phosphorylation and cardiomyocyte contractile properties may represent an impaired compensatory response, contributing to the detrimental DCM phenotype in children.NEW & NOTEWORTHY Our study is the first to demonstrate that cardiomyocytes from infants and young children with dilated cardiomyopathy (DCM) exhibit increased calcium sensitivity (likely mediated by decreased troponin I phosphorylation) compared with nonfailing pediatric cardiomyocytes. Compared with published values in adult cardiomyocytes, pediatric cardiomyocytes have notably decreased cooperativity, with a further reduction in the setting of DCM. Distinct adaptations in cardiomyocyte contractile properties may contribute to a differential response to pharmacological therapies in the pediatric DCM population.

Computational Studies of Cardiac and Skeletal Troponin

Troponin is a key regulatory protein in muscle contraction, consisting of three subunits troponin C (TnC), troponin I (TnI), and troponin T (TnT). Calcium association to TnC initiates contraction by causing a series of dynamic and conformational changes that allow the switch peptide of TnI to bind and subsequently cross bridges to form between the thin and thick filament of the sarcomere. Owing to its pivotal role in contraction regulation, troponin has been the focus of numerous computational studies over the last decade. These studies elegantly supplemented a large volume of experimental work and focused on the structure, dynamics and function of the whole troponin complex, individual subunits, and even on segments of the thin filament. Molecular dynamics, Brownian dynamics, and free energy simulations have been used to elucidate the conformational dynamics and underlying free energy landscape of troponin, calcium, and switch peptide binding, as well as the effect of disease mutations, small molecules and post-translational modifications such as phosphorylation. Frequently, simulations have been used to confirm or explain experimental observations. Computer-aided drug discovery tools have been employed to identify novel potential calcium sensitizing agents binding to the TnC-TnI interface. Finally, Markov modeling has contributed to simulating contraction within the sarcomere on the mesoscale. Here we are reviewing and classifying the existing computational work on troponin and its subunits, outline current gaps in simulations elucidating troponin's role in contraction and suggest potential future developments in the field.

Desensitizing mouse cardiac troponin C to calcium converts slow muscle towards a fast muscle phenotype

KEY POINTS

The Ca2+ -desensitizing D73N mutation in slow skeletal/cardiac troponin C caused dilatated cardiomyopathy in mice, but the consequences of this mutation in skeletal muscle were not known. The D73N mutation led to a rightward shift in the force versus pCa (-log [Ca]) relationship in slow-twitch mouse fibres. The D73N mutation led to a rightward shift in the force-stimulation frequency relationship and reduced fatigue resistance of mouse soleus muscle. The D73N mutation led to reduced cross-sectional area of slow-twitch fibres in mouse soleus muscle without affecting fibre type composition of the muscle. The D73N mutation resulted in significantly shorter times to peak force and to relaxation during isometric twitches and tetani in mouse soleus muscle. The D73N mutation led to major changes in physiological properties of mouse soleus muscle, converting slow muscle toward a fast muscle phenotype.

ABSTRACT

The missense mutation, D73N, in mouse cardiac troponin C has a profound impact on cardiac function, mediated by a decreased myofilament Ca2+ sensitivity. Mammalian cardiac muscle and slow skeletal muscle normally share expression of the same troponin C isoform. Therefore, the objective of this study was to determine the consequences of the D73N mutation in skeletal muscle, as a potential mechanism that contributes to the morbidity associated with heart failure or other conditions in which Ca2+ sensitivity might be altered. Effects of the D73N mutation on physiological properties of mouse soleus muscle, in which slow-twitch fibres are prevalent, were examined. The mutation resulted in a rightward shift of the force-stimulation frequency relationship, and significantly faster kinetics of isometric twitches and tetani in isolated soleus muscle. Furthermore, soleus muscles from D73N mice underwent a significantly greater reduction in force during a fatigue test. The mutation significantly reduced slow fibre mean cross-sectional area without affecting soleus fibre type composition. The effects of the mutation on Ca2+ sensitivity of force development in soleus skinned slow and fast fibres were also examined. As expected, the D73N mutation did not affect the Ca2+ sensitivity of force development in fast fibres but resulted in substantially decreased Ca2+ sensitivity in slow fibres. The results demonstrate that a point mutation in a single constituent of myofilaments (slow/cardiac troponin C) led to major changes in physiological properties of skeletal muscle and converted slow muscle toward a fast muscle phenotype with reduced fatigue resistance and Ca2+ sensitivity of force generation.

Molecular Dynamics and Umbrella Sampling Simulations Elucidate Differences in Troponin C Isoform and Mutant Hydrophobic Patch Exposure

Troponin C (TnC) facilitates muscle contraction through calcium-binding within its N-terminal region (NTnC). As previously observed using molecular dynamics (MD) simulations, this calcium-binding event leads to an increase in the dynamics of helices lining a hydrophobic patch on TnC. Simulation times of multiple microseconds were required to even see a partial opening of the hydrophobic patch, limiting the ability to thoroughly and quantitatively investigate these rare events. Here we describe the application of umbrella sampling to probe the TnC hydrophobic patch opening in a more targeted and quantitative fashion. Umbrella sampling was utilized to investigate the differences in the free energy of opening between cardiac (cTnC) and fast skeletal TnC (sTnC). We found that, in agreement with previous reports, holo (calcium-bound) sTnC had a lower free energy of opening compared with holo cTnC. Additionally, differences in the free energy of opening of hypertrophic (HCM) and dilated cardiomyopathy (DCM) cTnC systems were investigated. MD simulations and umbrella sampling revealed a lower free energy of opening for the HCM mutations A8V and A31S, as well as the calcium-sensitizing mutation L48Q. The DCM mutations, Y5H, Q50R, and E59D/D75Y, all exhibited a higher free energy of opening. An umbrella sampling simulation of cTnI-bound holo cTnC exhibited the lowest free energy in the open configuration, in agreement with experimental data. In conclusion, this study presents a novel and successful protocol for applying umbrella sampling simulations to quantitatively study the molecular basis of muscle contraction and proposes a mechanism by which HCM and DCM-associated mutations influence contraction.

Insights and Challenges of Multi-Scale Modeling of Sarcomere Mechanics in cTn and Tm DCM Mutants-Genotype to Cellular Phenotype

Dilated Cardiomyopathy (DCM) is a leading cause of sudden cardiac death characterized by impaired pump function and dilatation of cardiac ventricles. In this review we discuss various in silico approaches to elucidating the mechanisms of genetic mutations leading to DCM. The approaches covered in this review focus on bridging the spatial and temporal gaps that exist between molecular and cellular processes. Mutations in sarcomeric regulatory thin filament proteins such as the troponin complex (cTn) and Tropomyosin (Tm) have been associated with DCM. Despite the experimentally-observed myofilament measures of contractility in the case of these mutations, the mechanisms by which the underlying molecular changes and protein interactions scale up to organ failure by these mutations remains elusive. The review highlights multi-scale modeling approaches and their applicability to study the effects of sarcomeric gene mutations in-silico. We discuss some of the insights that can be gained from computational models of cardiac biomechanics when scaling from molecular states to cellular level.

Fluorescent Protein-Based Ca2+ Sensor Reveals Global, Divalent Cation-Dependent Conformational Changes in Cardiac Troponin C

Cardiac troponin C (cTnC) is a key effector in cardiac muscle excitation-contraction coupling as the Ca²⁺ sensing subunit responsible for controlling contraction. In this study, we generated several FRET sensors for divalent cations based on cTnC flanked by a donor fluorescent protein (CFP) and an acceptor fluorescent protein (YFP). The sensors report Ca²⁺ and Mg²⁺ binding, and relay global structural information about the structural relationship between cTnC’s N- and C-domains. The sensors were first characterized using end point titrations to decipher the response to Ca²⁺ binding in the presence or absence of Mg²⁺. The sensor that exhibited the largest responses in end point titrations, CTV-TnC, (Cerulean, TnC, and Venus) was characterized more extensively. Most of the divalent cation-dependent FRET signal originates from the high affinity C-terminal EF hands. CTV-TnC reconstitutes into skinned fiber preparations indicating proper assembly of troponin complex, with only ~0.2 pCa unit rightward shift of Ca²⁺-sensitive force development compared to WT-cTnC. Affinity of CTV-TnC for divalent cations is in agreement with known values for WT-cTnC. Analytical ultracentrifugation indicates that CTV-TnC undergoes compaction as divalent cations bind. C-terminal sites induce ion-specific (Ca²⁺ versus Mg²⁺) conformational changes in cTnC. Our data also provide support for the presence of additional, non-EF-hand sites on cTnC for Mg²⁺ binding. In conclusion, we successfully generated a novel FRET-Ca²⁺ sensor based on full length cTnC with a variety of cellular applications. Our sensor reveals global structural information about cTnC upon divalent cation binding.

Molecular Effects of cTnC DCM Mutations on Calcium Sensitivity and Myofilament Activation-An Integrated Multiscale Modeling Study

Mutations in cardiac troponin C (D75Y, E59D, and G159D), a key regulatory protein of myofilament contraction, have been associated with dilated cardiomyopathy (DCM). Despite reports of altered myofilament function in these mutants, the underlying molecular alterations caused by these mutations remain elusive. Here we investigate in silico the intramolecular mechanisms by which these mutations affect myofilament contraction. On the basis of the location of cardiac troponin C (cTnC) mutations, we tested the hypothesis that intramolecular effects can explain the altered myofilament calcium sensitivity of force development for D75Y and E59D cTnC, whereas altered cardiac troponin C-troponin I (cTnC-cTnI) interaction contributes to the reported contractile effects of the G159D mutation. We employed a multiscale approach combining molecular dynamics (MD) and Brownian dynamics (BD) simulations to estimate cTnC calcium association and hydrophobic patch opening. We then integrated these parameters into a Markov model of myofilament activation to compute the steady-state force-pCa relationship. The analysis showed that myofilament calcium sensitivity with D75Y and E59D can be explained by changes in calcium binding affinity of cTnC and the rate of hydrophobic patch opening, if a partial cTnC interhelical opening angle (110°) is sufficient for cTnI switch peptide association to cTnC. In contrast, interactions between cTnC and cTnI within the cardiac troponin complex must also be accounted for to explain contractile alterations due to G159D. In conclusion, this is the first multiscale in silico study to elucidate how direct molecular effects of genetic mutations in cTnC translate to altered myofilament contractile function.

Expression signature of lncRNAs and their potential roles in cardiac fibrosis of post-infarct mice

The present study aimed to investigate whether long non-coding RNAs (lncRNAs) are involved in cardiac fibrogenesis induced by myocardial infarction (MI). The differentially expressed lncRNAs and mRNAs in peri-infarct region of mice 4 weeks after MI were selected for bioinformatic analysis including gene ontology (GO) enrichment, pathway and network analysis. Left ventricular tissue levels of lncRNAs and mRNAs were compared between MI and sham control mice, using a false discovery rate (FDR) of <5%. Out of 55000 lncRNAs detected, 263 were significantly up-regulated and 282 down-regulated. Out of 23000 mRNAs detected, 142 were significantly up-regulated and 67 down-regulated. Among the differentially expressed lncRNAs, 53 were up-regulated by ≥2.0-fold change and 37 down-regulated by ≤0.5-fold change. Nine up-regulated and five down-regulated lncRNAs were randomly selected for quantitative real-time PCR (qRT-PCR) verification. GO and pathway analyses revealed 173 correlated lncRNA-mRNA pairs for 57 differentially expressed lncRNAs and 20 differentially expressed genes which are related to the development of cardiac fibrosis. We identified TGF-β3 as the top-ranked gene, a critical component of the transforming growth factor-β (TGF-β) and mitogen activated protein kinase (MAPK) signalling pathways in cardiac fibrosis. NONMMUT022554 was identified as the top-ranked lncRNA, positively correlated with six up-regulated genes, which are involved in the extracellular matrix (ECM)-receptor interactions and the phosphoinositid-3 kinase/protein kinase B (PI3K-Akt) signalling pathway. Our study has identified the expression signature of lncRNAs in cardiac fibrosis induced by MI and unravelled the possible involvement of the deregulated lncRNAs in cardiac fibrosis and the associated pathological processes.

Structure and function of cardiac troponin C (TNNC1): Implications for heart failure, cardiomyopathies, and troponin modulating drugs

In striated muscle, the protein troponin complex turns contraction on and off in a calcium-dependent manner. The calcium-sensing component of the complex is troponin C, which is expressed from the TNNC1 gene in both cardiac muscle and slow-twitch skeletal muscle (identical transcript in both tissues) and the TNNC2 gene in fast-twitch skeletal muscle. Cardiac troponin C (cTnC) is made up of two globular EF-hand domains connected by a flexible linker. The structural C-domain (cCTnC) contains two high affinity calcium-binding sites that are always occupied by Ca(2+) or Mg(2+) under physiologic conditions, stabilizing an open conformation that remains anchored to the rest of the troponin complex. In contrast, the regulatory N-domain (cNTnC) contains a single low affinity site that is largely unoccupied at resting calcium concentrations. During muscle activation, calcium binding to cNTnC favors an open conformation that binds to the switch region of troponin I, removing adjacent inhibitory regions of troponin I from actin and allowing muscle contraction to proceed. Regulation of the calcium binding affinity of cNTnC is physiologically important, because it directly impacts the calcium sensitivity of muscle contraction. Calcium sensitivity can be modified by drugs that stabilize the open form of cNTnC, post-translational modifications like phosphorylation of troponin I, or downstream thin filament protein interactions that impact the availability of the troponin I switch region. Recently, mutations in cTnC have been associated with hypertrophic or dilated cardiomyopathy. A detailed understanding of how calcium sensitivity is regulated through the troponin complex is necessary for explaining how mutations perturb its function to promote cardiomyopathy and how post-translational modifications in the thin filament affect heart function and heart failure. Troponin modulating drugs are being developed for the treatment of cardiomyopathies and heart failure.

Cell biology of sarcomeric protein engineering: disease modeling and therapeutic potential

The cardiac sarcomere is the functional unit for myocyte contraction. Ordered arrays of sarcomeric proteins, held in stoichiometric balance with each other, respond to calcium to coordinate contraction and relaxation of the heart. Altered sarcomeric structure-function underlies the primary basis of disease in multiple acquired and inherited heart disease states. Hypertrophic and restrictive cardiomyopathies are caused by inherited mutations in sarcomeric genes and result in altered contractility. Ischemia-mediated acidosis directly alters sarcomere function resulting in decreased contractility. In this review, we highlight the use of acute genetic engineering of adult cardiac myocytes through stoichiometric replacement of sarcomeric proteins in these disease states with particular focus on cardiac troponin I. Stoichiometric replacement of disease causing mutations has been instrumental in defining the molecular mechanisms of hypertrophic and restrictive cardiomyopathy in a cellular context. In addition, taking advantage of stoichiometric replacement through gene therapy is discussed, highlighting the ischemia-resistant histidine-button, A164H cTnI. Stoichiometric replacement of sarcomeric proteins offers a potential gene therapy avenue to replace mutant proteins, alter sarcomeric responses to pathophysiologic insults, or neutralize altered sarcomeric function in disease.

A device for rapid and quantitative measurement of cardiac myocyte contractility

Cardiac contractility is the hallmark of cardiac function and is a predictor of healthy or diseased cardiac muscle. Despite advancements over the last two decades, the techniques and tools available to cardiovascular scientists are limited in their utility to accurately and reliably measure the amplitude and frequency of cardiomyocyte contractions. Isometric force measurements in the past have entailed cumbersome attachment of isolated and permeabilized cardiomyocytes to a force transducer followed by measurements of sarcomere lengths under conditions of submaximal and maximal Ca(2+) activation. These techniques have the inherent disadvantages of being labor intensive and costly. We have engineered a micro-machined cantilever sensor with an embedded deflection-sensing element that, in preliminary experiments, has demonstrated to reliably measure cardiac cell contractions in real-time. Here, we describe this new bioengineering tool with applicability in the cardiovascular research field to effectively and reliably measure cardiac cell contractility in a quantitative manner. We measured contractility in both primary neonatal rat heart cardiomyocyte monolayers that demonstrated a beat frequency of 3 Hz as well as human embryonic stem cell-derived cardiomyocytes with a contractile frequency of about 1 Hz. We also employed the β-adrenergic agonist isoproterenol (100 nmol l(-1)) and observed that our cantilever demonstrated high sensitivity in detecting subtle changes in both chronotropic and inotropic responses of monolayers. This report describes the utility of our micro-device in both basic cardiovascular research as well as in small molecule drug discovery to monitor cardiac cell contractions.

Gene regulation, alternative splicing, and posttranslational modification of troponin subunits in cardiac development and adaptation: a focused review

Troponin plays a central role in regulating the contraction and relaxation of vertebrate striated muscles. This review focuses on the isoform gene regulation, alternative RNA splicing, and posttranslational modifications of troponin subunits in cardiac development and adaptation. Transcriptional and posttranscriptional regulations such as phosphorylation and proteolysis modifications, and structure-function relationships of troponin subunit proteins are summarized. The physiological and pathophysiological significances are discussed for impacts on cardiac muscle contractility, heart function, and adaptations in health and diseases.

Biochemical characterisation of Troponin C mutations causing hypertrophic and dilated cardiomyopathies

Cardiac muscle contraction occurs through an interaction of the myosin head with the actin filaments, a process which is regulated by the troponin complex together with tropomyosin and is Ca(2+) dependent. Mutations in genes encoding sarcomeric proteins are a common cause of familial hypertrophic and dilated cardiomyopathies. The scope of this review is to gather information from studies regarding the in vitro characterisation of six HCM and six DCM mutations on the cardiac TnC gene and to suggest, if possible, how they may lead to dysfunction. Since TnC is the subunit responsible for Ca(2+) binding, mutations in the TnC could possibly have a strong impact on Ca(2+) binding affinities. Furthermore, the interactions of mutant TnCs with their binding partners could be altered. From the characterisation studies available to date, we can conclude that the HCM mutations on TnC increase significantly the Ca(2+) sensitivity of force development or of ATPase activity, producing large pCa shifts in comparison to WT TnC. In contrast, the DCM mutations on TnC have a tendency to decrease the Ca(2+) sensitivity of force development or of ATPase activity in comparison to WT TnC. Furthermore, the DCM mutants of TnC are not responsive to the TnI phosphorylation signal resulting in filaments that preserve their Ca(2+) sensitivity in contrast to WT filaments that experience a decrease in Ca(2+) sensitivity upon TnI phosphorylation.

Thin filament incorporation of an engineered cardiac troponin C variant (L48Q) enhances contractility in intact cardiomyocytes from healthy and infarcted hearts

Many current pharmaceutical therapies for systolic heart failure target intracellular [Ca(2+)] ([Ca(2+)]i) metabolism, or cardiac troponin C (cTnC) on thin filaments, and can have significant side-effects, including arrhythmias or adverse effects on diastolic function. In this study, we tested the feasibility of directly increasing the Ca(2+) binding properties of cTnC to enhance contraction independent of [Ca(2+)]i in intact cardiomyocytes from healthy and myocardial infarcted (MI) hearts. Specifically, cardiac thin filament activation was enhanced through adenovirus-mediated over-expression of a cardiac troponin C (cTnC) variant designed to have increased Ca(2+) binding affinity conferred by single amino acid substitution (L48Q). In skinned cardiac trabeculae and myofibrils we and others have shown that substitution of L48Q cTnC for native cTnC increases Ca(2+) sensitivity of force and the maximal rate of force development. Here we introduced L48Q cTnC into myofilaments of intact cardiomyocytes via adeno-viral transduction to deliver cDNA for the mutant or wild type (WT) cTnC protein. Using video-microscopy to monitor cell contraction, relaxation, and intracellular Ca(2+) transients (Fura-2), we report that incorporation of L48Q cTnC significantly increased contractility of cardiomyocytes from healthy and MI hearts without adversely affecting Ca(2+) transient properties or relaxation. The improvements in contractility from L48Q cTnC expression are likely the result of enhanced contractile efficiency, as intracellular Ca(2+) transient amplitudes were not affected. Expression and incorporation of L48Q cTnC into myofilaments was confirmed by Western blot analysis of myofibrils from transduced cardiomyocytes, which indicated replacement of 18±2% of native cTnC with L48Q cTnC. These experiments demonstrate the feasibility of directly targeting cardiac thin filament proteins to enhance cardiomyocyte contractility that is impaired following MI.

Myofilament dysfunction as an emerging mechanism of volume overload heart failure

Two main hemodynamic overload mechanisms [i.e., volume and pressure overload (VO and PO, respectively] result in heart failure (HF), and these two mechanisms have divergent pathologic alterations and different pathophysiological mechanisms. Extensive evidence from animal models and human studies of PO demonstrate a clear association with alterations in Ca(2+) homeostasis. By contrast, emerging evidence from animal models and patients with regurgitant valve disease and dilated cardiomyopathy point toward a more prominent role of myofilament dysfunction. With respect to VO HF, key features of excitation-contraction coupling defects, myofilament dysfunction, and extracellular matrix composition will be discussed.

Inherited cardiomyopathies caused by troponin mutations

Genetic investigations of cardiomyopathy in the recent two decades have revealed a large number of mutations in the genes encoding sarcomeric proteins as a cause of inherited hypertrophic cardiomyopathy (HCM), dilated cardiomyopathy (DCM), or restrictive cardiomyopathy (RCM). Most functional analyses of the effects of mutations on cardiac muscle contraction have revealed significant changes in the Ca(2+)-regulatory mechanism, in which cardiac troponin (cTn) plays important structural and functional roles as a key regulatory protein. Over a hundred mutations have been identified in all three subunits of cTn, i.e., cardiac troponins T, I, and C. Recent studies on cTn mutations have provided plenty of evidence that HCM- and RCM-linked mutations increase cardiac myofilament Ca(2+) sensitivity, while DCM-linked mutations decrease it. This review focuses on the functional consequences of mutations found in cTn in terms of cardiac myofilament Ca(2+) sensitivity, ATPase activity, force generation, and cardiac troponin I phosphorylation, to understand potential molecular and cellular pathogenic mechanisms of the three types of inherited cardiomyopathy.

Structural and functional consequences of cardiac troponin C L57Q and I61Q Ca(2+)-desensitizing variants

Two cTnC variants, L57Q and I61Q, both of which are located on helix C within the N domain of cTnC, were originally reported in the skeletal muscle system [Tikunova, Davis, J. Biol. Chem. 279 (2004) 35341-35352], as the analogous L58Q and I62Q sTnC, and demonstrated a decreased Ca(2+) binding affinity. Here, we provide detailed characterization of structure-function relationships for these two cTnC variants, to determine if they behave differently in the cardiac system and as a framework for determining similarities and differences with other cTnC mutations that have been associated with DCM. We have used an integrative approach to study the structure and function of these cTnC variants both in solution and in silico, to understand how the L57Q and I61Q mutations influence Ca(2+) binding at site II, the subsequent effects on the interaction with cTnI, and the structural changes which are associated with these changes. Steady-state and stopped flow fluorescence spectroscopy confirmed that a decrease in Ca(2+) affinity for recombinant cTnC and cTn complexes containing the L57Q or I61Q variants. The L57Q variant was intermediate between WT and I61Q cTnC and also did not significantly alter cTnC-cTnI interaction in the absence of Ca(2+), but did decrease the interaction in the presence of Ca(2+). In contrast, I61Q decreased the cTnC-cTnI interaction in both the absence and presence of Ca(2+). This difference in the absence of Ca(2+) suggests a greater structural change in cNTnC may occur with the I61Q mutation than the L57Q mutation. MD simulations revealed that the decreased Ca(2+) binding induced by I61Q may result from destabilization of the Ca(2+) binding site through interruption of intra-molecular interactions when residue 61 forms new hydrogen bonds with G70 on the Ca(2+) binding loop. The experimentally observed interruption of the cTnC-cTnI interaction caused by L57Q or I61Q is due to the disruption of key hydrophobic interactions between helices B and C in cNTnC. This study provides a molecular basis of how single mutations in the C helix of cTnC can reduce Ca(2+) binding affinity and cTnC-cTnI interaction, which may provide useful insights for a better understanding of cardiomyopathies and future gene-based therapies.

Molecular basis of calcium-sensitizing and desensitizing mutations of the human cardiac troponin C regulatory domain: a multi-scale simulation study

Troponin C (TnC) is implicated in the initiation of myocyte contraction via binding of cytosolic Ca²⁺ and subsequent recognition of the Troponin I switch peptide. Mutations of the cardiac TnC N-terminal regulatory domain have been shown to alter both calcium binding and myofilament force generation. We have performed molecular dynamics simulations of engineered TnC variants that increase or decrease Ca²⁺ sensitivity, in order to understand the structural basis of their impact on TnC function. We will use the distinction for mutants that are associated with increased Ca²⁺ affinity and for those mutants with reduced affinity. Our studies demonstrate that for GOF mutants V44Q and L48Q, the structure of the physiologically-active site II Ca²⁺ binding site in the Ca²⁺-free (apo) state closely resembled the Ca²⁺-bound (holo) state. In contrast, site II is very labile for LOF mutants E40A and V79Q in the apo form and bears little resemblance with the holo conformation. We hypothesize that these phenomena contribute to the increased association rate, k(on), for the GOF mutants relative to LOF. Furthermore, we observe significant positive and negative positional correlations between helices in the GOF holo mutants that are not found in the LOF mutants. We anticipate these correlations may contribute either directly to Ca²⁺ affinity or indirectly through TnI association. Our observations based on the structure and dynamics of mutant TnC provide rationale for binding trends observed in GOF and LOF mutants and will guide the development of inotropic drugs that target TnC.

Expression and assembly of active human cardiac troponin in Escherichia coli

Cardiomyopathy-related mutations in human cardiac troponin subunits, including troponin C (hcTnC), troponin I (hcTnI), and troponin T (hcTnT), are well-documented. Recently, it has been recognised that human cardiac troponin (hcTn) is a sophisticated allosteric system. Therefore, the effect of drugs on this protein complex should be studied with assembled hcTn rather than a short fragment of a subunit or the subunit itself. Here, we describe the expression and assembly of active hcTn in Escherichia coli, a novel method that is rapid and simple, and produces large amounts of functional hcTn.

Beyond the cardiac myofilament: hypertrophic cardiomyopathy- associated mutations in genes that encode calcium-handling proteins

Traditionally regarded as a genetic disease of the cardiac sarcomere, hypertrophic cardiomyopathy (HCM) is the most common inherited cardiovascular disease and a significant cause of sudden cardiac death. While the most common etiologies of this phenotypically diverse disease lie in a handful of genes encoding critical contractile myofilament proteins, approximately 50% of patients diagnosed with HCM worldwide do not host sarcomeric gene mutations. Recently, mutations in genes encoding calcium-sensitive and calcium-handling proteins have been implicated in the pathogenesis of HCM. Among these are mutations in TNNC1- encoded cardiac troponin C, PLN-encoded phospholamban, and JPH2-encoded junctophilin 2 which have each been associated with HCM in multiple studies. In addition, mutations in RYR2-encoded ryanodine receptor 2, CASQ2-encoded calsequestrin 2, CALR3-encoded calreticulin 3, and SRI-encoded sorcin have been associated with HCM, although more studies are required to validate initial findings. While a relatively uncommon cause of HCM, mutations in genes that encode calcium-handling proteins represent an emerging genetic subset of HCM. Furthermore, these naturally occurring disease-associated mutations have provided useful molecular tools for uncovering novel mechanisms of disease pathogenesis, increasing our understanding of basic cardiac physiology, and dissecting important structure-function relationships within these proteins.

Long-timescale molecular dynamics simulations elucidate the dynamics and kinetics of exposure of the hydrophobic patch in troponin C

Troponin (Tn) is an important regulatory protein in the thin-filament complex of cardiomyocytes. Calcium binding to the troponin C (TnC) subunit causes a change in its dynamics that leads to the transient opening of a hydrophobic patch on TnC's surface, to which a helix of another subunit, troponin I (TnI), binds. This process initiates contraction, making it an important target for studies investigating the detailed molecular processes that underlie contraction. Here we use microsecond-timescale Anton molecular dynamics simulations to investigate the dynamics and kinetics of the opening transition of the TnC hydrophobic patch. Free-energy differences for opening are calculated for wild-type Ca(2+)-bound TnC (∼8 kcal/mol), V44Q Ca(2+)-bound TnC (3.2 kcal/mol), E40A Ca(2+)-bound TnC (∼12 kcal/mol), and wild-type apo TnC (∼20 kcal/mol). These results suggest that the mutations have a profound impact on the frequency with which the hydrophobic patch presents to TnI. In addition, these simulations corroborate that cardiac wild-type TnC does not open on timescales relevant to contraction without calcium being bound.

Cardiomyopathy-Related Mutations in Cardiac Troponin C, L29Q and G159D, Have Divergent Effects on Rat Cardiac Myofiber Contractile Dynamics

Previous studies of cardiomyopathy-related mutations in cardiac troponin C (cTnC)-L29Q and G159D-have shown diverse findings. The link between such mutant effects and their divergent impact on cardiac phenotypes has remained elusive due to lack of studies on contractile dynamics. We hypothesized that a cTnC mutant-induced change in the thin filament will affect global myofilament mechanodynamics because of the interactions of thin filament kinetics with both Ca(2+) binding and crossbridge (XB) cycling kinetics. We measured pCa-tension relationship and contractile dynamics in detergent-skinned rat cardiac papillary muscle fibers reconstituted with the recombinant wild-type rat cTnC (cTnC(WT)), cTnC(L29Q), and cTnC(G159D) mutants. cTnC(L29Q) fibers demonstrated a significant decrease in Ca(2+) sensitivity, but cTnC(G159D) fibers did not. Both mutants had no effect on Ca(2+)-activated maximal tension. The rate of XB recruitment dynamics increased in cTnC(L29Q) (26%) and cTnC(G159D) (25%) fibers. The rate of XB distortion dynamics increased in cTnC(G159D) fibers (15%). Thus, the cTnC(L29Q) mutant modulates the equilibrium between the non-cycling and cycling pool of XB by affecting the on/off kinetics of the regulatory units (Tropomyosin-Troponin); whereas, the cTnC(G159D) mutant increases XB cycling rate. Different effects on contractile dynamics may offer clue regarding how cTnC(L29Q) and cTnC(G159D) cause divergent effects on cardiac phenotypes.

Dynamics and calcium association to the N-terminal regulatory domain of human cardiac troponin C: a multiscale computational study

Troponin C (TnC) is an important regulatory molecule in cardiomyocytes. Calcium binding to site II in TnC initiates a series of molecular events that result in muscle contraction. The most direct change upon Ca(2+) binding is an opening motion of the molecule that exposes a hydrophobic patch on the surface allowing for Troponin I to bind. Molecular dynamics simulations were used to elucidate the dynamics of this crucial protein in three different states: apo, Ca(2+)-bound, and Ca(2+)-TnI-bound. Dynamics between the states are compared, and the Ca(2+)-bound system is investigated for opening motions. On the basis of the simulations, NMR chemical shifts and order parameters are calculated and compared with experimental observables. Agreement indicates that the simulations sample the relevant dynamics of the system. Brownian dynamics simulations are used to investigate the calcium association of TnC. We find that calcium binding gives rise to correlative motions involving the EF hand and collective motions conducive of formation of the TnI-binding interface. We furthermore indicate the essential role of electrostatic steering in facilitating diffusion-limited binding of Ca(2+).

Multi-Scale Continuum Modeling of Biological Processes: From Molecular Electro-Diffusion to Sub-Cellular Signaling Transduction

This article provides a brief review of multi-scale modeling at the molecular to cellular scale, with new results for heart muscle cells. A finite element-based simulation package (SMOL) was used to investigate the signaling transduction at molecular and sub-cellular scales (http://mccammon.ucsd.edu/smol/, http://FETK.org) by numerical solution of time-dependent Smoluchowski equations and a reaction-diffusion system. At the molecular scale, SMOL has yielded experimentally-validated estimates of the diffusion-limited association rates for the binding of acetylcholine to mouse acetylcholinesterase using crystallographic structural data. The predicted rate constants exhibit increasingly delayed steady-state times with increasing ionic strength and demonstrate the role of an enzyme's electrostatic potential in influencing ligand binding. At the sub-cellular scale, an extension of SMOL solves a non-linear, reaction-diffusion system describing Ca²⁺ ligand buffering and diffusion in experimentally-derived rodent ventricular myocyte geometries. Results reveal the important role for mobile and stationary Ca²⁺ buffers, including Ca²⁺ indicator dye. We found that the alterations in Ca²⁺-binding and dissociation rates of troponin C (TnC) and total TnC concentration modulate subcellular Ca²⁺ signals. Model predicts that reduced off-rate in whole troponin complex (TnC, TnI, TnT) versus reconstructed thin filaments (Tn, Tm, actin) alters cytosolic Ca²⁺ dynamics under control conditions or in disease-linked TnC mutations. The ultimate goal of these studies is to develop scalable methods and theories for integration of molecular-scale information into simulations of cellular-scale systems.

Correlation of molecular and functional effects of mutations in cardiac troponin T linked to familial hypertrophic cardiomyopathy: an integrative in silico/in vitro approach

Nearly 70% of all of the known cTnT mutations that cause familial hypertrophic cardiomyopathy fall within the TNT1 region that is critical to cTn-Tm binding. The high resolution structure of this domain has not been determined, and this lack of information has hindered structure-function analysis. In the current study, a coupled computational experimental approach was employed to correlate changes in cTnT dynamics to basic function using the regulated in vitro motility assay (R-IVM). An in silico approach to calculate forces in terms of a bending coordinate was used to precisely identify decreases in bending forces at residues 105 and 106 within the proposed cTnT "hinge" region. Significant functional changes were observed in multiple functional properties, including a decrease in the cooperativity of calcium activation, the calcium sensitivity of sliding speed, and maximum sliding speed. Correlation of the computational and experimental findings revealed an association between TNT1 flexibility and the cooperativity of thin filament calcium activation where an increase in flexibility led to a decrease in cooperativity. Further analysis of the primary sequence of the TNT1 region revealed a unique pattern of conserved charged TNT1 residues altered by the R92W and R92L mutations and may represent the underlying "structure" modulating this central functional domain. These data provide a framework for further integrated in silico/in vitro approaches that may be extended into a high-throughput predictive screen to overcome the current structural limitations in linking molecular phenotype to genotype in thin filament cardiomyopathies.

Functional characterization of TNNC1 rare variants identified in dilated cardiomyopathy

TNNC1, which encodes cardiac troponin C (cTnC), remains elusive as a dilated cardiomyopathy (DCM) gene. Here, we report the clinical, genetic, and functional characterization of four TNNC1 rare variants (Y5H, M103I, D145E, and I148V), all previously reported by us in association with DCM (Hershberger, R. E., Norton, N., Morales, A., Li, D., Siegfried, J. D., and Gonzalez-Quintana, J. (2010) Circ. Cardiovasc. Genet. 3, 155-161); in the previous study, two variants (Y5H and D145E) were identified in subjects who also carried MYH7 and MYBPC3 rare variants, respectively. Functional studies using the recombinant human mutant cTnC proteins reconstituted into porcine papillary skinned fibers showed decreased Ca(2+) sensitivity of force development (Y5H and M103I). Furthermore, the cTnC mutants diminished (Y5H and I148V) or abolished (M103I) the effects of PKA phosphorylation on Ca(2+) sensitivity. Only M103I decreased the troponin activation properties of the actomyosin ATPase when Ca(2+) was present. CD spectroscopic studies of apo (absence of divalent cations)-, Mg(2+)-, and Ca(2+)/Mg(2+)-bound states indicated that all of the cTnC mutants (except I148V in the Ca(2+)/Mg(2+) condition) decreased the α-helical content. These results suggest that each mutation alters the function/ability of the myofilament to bind Ca(2+) as a result of modifications in cTnC structure. One variant (D145E) that was previously reported in association with hypertrophic cardiomyopathy and that produced results in vivo in this study consistent with prior hypertrophic cardiomyopathy functional studies was found associated with the MYBPC3 P910T rare variant, likely contributing to the observed DCM phenotype. We conclude that these rare variants alter the regulation of contraction in some way, and the combined clinical, molecular, genetic, and functional data reinforce the importance of TNNC1 rare variants in the pathogenesis of DCM.

Thin filament mutations: developing an integrative approach to a complex disorder

Sixteen years ago, mutations in cardiac troponin (Tn)T and α-tropomyosin were linked to familial hypertrophic cardiomyopathy, thus transforming the disorder from a disease of the β-myosin heavy chain to a disease of the cardiac sarcomere. From the outset, studies suggested that mutations in the regulatory thin filament caused a complex, heterogeneous pattern of ventricular remodeling with wide variations in clinical expression. To date, the clinical heterogeneity is well matched by an extensive array of nearly 100 independent mutations in all components of the cardiac thin filament. Significant advances in our understanding of the biophysics of myofilament activation, coupled to the emerging evidence that thin filament linked cardiomyopathies are progressive, suggests that a renewed focus on the most proximal events in both the molecular and clinical pathogenesis of the disease will be necessary to achieve the central goal of using genotype information to manage affected patients. In this review, we examine the existing biophysical and clinical evidence in support of a more proximal definition of thin filament cardiomyopathies. In addition, new high-resolution, integrated approaches are presented to help define the way forward as the field works toward developing a more robust link between genotype and phenotype in this complex disorder.

Iterative Optimization of Molecular Mechanics Force Fields from NMR Data of Full-Length Proteins

High quality molecular mechanics force fields of proteins are key for the quantitative interpretation of experimental data and the predictive understanding of protein function based on computer simulations. A strategy is presented for the optimization of protein force fields based on full-length proteins in their native environment that is guided by experimental NMR chemical shifts and residual dipolar couplings (RDCs). An energy-based reweighting approach is applied to a long molecular dynamics trajectory, performed with a parent force field, to efficiently screen a large number of trial force fields. The force field that yields the best agreement with the experimental data is then used as the new parent force field, and the procedure is repeated until no further improvement is obtained. This method is demonstrated for the optimization of the backbone φ,ψ dihedral angle potential of the Amber ff99SB force field using six trial proteins and another 17 proteins for cross-validation using (13)C chemical shifts with and without backbone RDCs. The φ,ψ dihedral angle potential is systematically improved by the inclusion of correlation effects through the addition of up to 24 bivariate Gaussian functions of variable height, width, and tilt angle. The resulting force fields, termed ff99SB_φψ(g24;CS) and ff99SB_φψ(g8;CS,RDC), perform significantly better than their parent force field in terms of both NMR data reproduction and Cartesian coordinate root-mean-square deviations between the MD trajectories and the X-ray crystal structures. The strategy introduced here represents a powerful addition to force field optimization approaches by overcoming shortcomings of methods that are solely based on quantum-chemical calculations of small molecules and protein fragments in the gas phase.

Molecular dynamics studies on troponin (TnI-TnT-TnC) complexes: insight into the regulation of muscle contraction

Mutations of any subunit of the troponin complex may lead to serious disorders. Rational approaches to managing these disorders require knowledge of the complex interactions among the three subunits that are required for proper function. Molecular dynamics (MD) simulations were performed for both skeletal (sTn) and cardiac (cTn) troponin. The interactions and correlated motions among the three components of the troponin complex were analyzed using both Molecular Mechanics-Generalized Born Surface Area (MMGBSA) and cross-correlation techniques. The TnTH2 helix was strongly positively correlated with the two long helices of TnI. The C domain of TnC was positively correlated with TnI and TnT. The N domain of TnC was negatively correlated with TnI and TnT in cTn, but not in sTn. The two C-domain calcium-binding sites of TnC were dynamically correlated. The two regulatory N-domain calcium-binding sites of TnC were dynamically correlated, even though the calcium-binding site I is dysfunctional. The strong interaction residue pairs and the strong dynamically correlated residues pairs among the three components of troponin complexes were identified. These correlated motions are consistent with the idea that there is a high degree of cooperativity among the components of the regulatory complex in response to Ca(2+) and other effectors. This approach may give insight into the mechanism by which mutations of troponin cause disease. It is interesting that some observed disease causing mutations fall within regions of troponin that are strongly correlated or interacted.

Calcium binding kinetics of troponin C strongly modulate cooperative activation and tension kinetics in cardiac muscle

Tension development and relaxation in cardiac muscle are regulated at the thin filament via Ca(2+) binding to cardiac troponin C (cTnC) and strong cross-bridge binding. However, the influence of cTnC Ca(2+)-binding properties on these processes in the organized structure of cardiac sarcomeres is not well-understood and likely differs from skeletal muscle. To study this we generated single amino acid variants of cTnC with altered Ca(2+) dissociation rates (k(off)), as measured in whole troponin (cTn) complex by stopped-flow spectroscopy (I61Q cTn>WT cTn>L48Q cTn), and exchanged them into cardiac myofibrils and demembranated trabeculae. In myofibrils at saturating Ca(2+), L48Q cTnC did not affect maximum tension (T(max)), thin filament activation (k(ACT)) and tension development (k(TR)) rates, or the rates of relaxation, but increased duration of slow phase relaxation. In contrast, I61Q cTnC reduced T(max), k(ACT) and k(TR) by 40-65% with little change in relaxation. Interestingly, k(ACT) was less than k(TR) with I61Q cTnC, and this difference increased with addition of inorganic phosphate, suggesting that reduced cTnC Ca(2+)-affinity can limit thin filament activation kinetics. Trabeculae exchanged with I61Q cTn had reduced T(max), Ca(2+) sensitivity of tension (pCa(50)), and slope (n(H)) of tension-pCa, while L48Q cTn increased pCa(50) and reduced n(H). Increased cross-bridge cycling with 2-deoxy-ATP increased pCa(50) with WT or L48Q cTn, but not I61Q cTn. We discuss the implications of these results for understanding the role of cTn Ca(2+)-binding properties on the magnitude and rate of tension development and relaxation in cardiac muscle.