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

Modulation of Structure and Dynamics of Cardiac Troponin by Phosphorylation and Mutations Revealed by Molecular Dynamics Simulations

Adrenaline acts on β1 receptors in the heart muscle to enhance contractility, increase the heart rate, and increase the rate of relaxation (lusitropy) via activation of the cyclic AMP-dependent protein kinase, PKA. Phosphorylation of serines 22 and 23 in the N-terminal peptide of cardiac troponin I is responsible for lusitropy. Mutations associated with cardiomyopathy suppress the phosphorylation-dependent change. Key parts of troponin responsible for this modulatory system are disordered and cannot be resolved by conventional structural approaches. We performed all-atom molecular dynamics simulations (5 × 1.5 μs runs) of the troponin core (419 amino acids) in the presence of Ca2+ in the bisphosphorylated and unphosphorylated states for both wild-type troponin and the troponin C (cTnC) G159D mutant. PKA phosphorylation affects troponin dynamics. There is significant rigidification of the structure involving rearrangement of the cTnI(1-33)-cTnC interaction and changes in the distribution of the cTnC helix A/B angle, troponin I (cTnI) switch peptide (149-164) docking, and the angle between the regulatory head and ITC arm domains. The familial dilated cardiomyopathy cTnC G159D mutation whose Ca2+ sensitivity is not modulated by cTnI phosphorylation exhibits a structure inherently more rigid than the wild type, with phosphorylation reversing the direction of all metrics relative to the wild type.

Targeting Troponin C with Small Molecules Containing Diphenyl Moieties: Calcium Sensitivity Effects on Striated Muscles and Structure-Activity Relationship

Despite large investments from academia and industry, heart failure, which results from a disruption of the contractile apparatus, remains a leading cause of death. Cardiac muscle contraction is a calcium-dependent mechanism, which is regulated by the troponin protein complex (cTn) and specifically by the N-terminal domain of its calcium-binding subunit (cNTnC). There is an increasing need for the development of small molecules that increase calcium sensitivity without altering the systolic calcium concentration, thereby strengthening the cardiac function. Here, we examined the effect of our previously identified calcium-sensitizing small molecule, ChemBridge compound 7930079, in the context of several homologous muscle systems. The effect of this molecule on force generation in isolated cardiac trabeculae and slow skeletal muscle fibers was measured. Furthermore, we explored the use of Gaussian accelerated molecular dynamics in sampling highly predictive receptor conformations based on NMR-derived starting structures. Additionally, we took a rational computational approach for lead optimization based on lipophilic diphenyl moieties. This integrated structural-biochemical-physiological approach led to the identification of three novel low-affinity binders, which had similar binding affinities to the known positive inotrope trifluoperazine. The most potent identified calcium sensitizer was compound 16 with an apparent affinity of 117 ± 17 μM.

Targeting Troponin C with Small Molecules Containing Diphenyl Moieties: Calcium Sensitivity Effects on Striated Muscle and Structure Activity Relationship

Despite large investments from academia and industry, heart failure, which results from a disruption of the contractile apparatus, remains a leading cause of death. Cardiac muscle contraction is a calcium-dependent mechanism, which is regulated by the troponin protein complex (cTn) and specifically by the N-terminal domain of its calcium binding subunit (cNTnC). There is an increasing need for the development of small molecules that increase calcium sensitivity without altering systolic calcium concentration, thereby strengthening cardiac function. Here, we examined the effect of our previously identified calcium sensitizing small molecule, ChemBridge compound 7930079, in the context of several homologous muscle systems. The effect of this molecule on force generation in isolated cardiac trabeculae and slow skeletal muscle fibers was measured. Furthermore, we explored the use of Gaussian accelerated molecular dynamics in sampling highly predictive receptor conformations based on NMR derived starting structures. Additionally, we took a rational computational approach for lead optimization based on lipophilic diphenyl moieties. This led to the identification of three novel low affinity binders, which had similar binding affinities to known positive inotrope trifluoperazine. The most potent identified calcium sensitizer was compound 16 with an apparent affinity of 117 ± 17 μM .

Myofilament-associated proteins with intrinsic disorder (MAPIDs) and their resolution by computational modeling

The cardiac sarcomere is a cellular structure in the heart that enables muscle cells to contract. Dozens of proteins belong to the cardiac sarcomere, which work in tandem to generate force and adapt to demands on cardiac output. Intriguingly, the majority of these proteins have significant intrinsic disorder that contributes to their functions, yet the biophysics of these intrinsically disordered regions (IDRs) have been characterized in limited detail. In this review, we first enumerate these myofilament-associated proteins with intrinsic disorder (MAPIDs) and recent biophysical studies to characterize their IDRs. We secondly summarize the biophysics governing IDR properties and the state-of-the-art in computational tools toward MAPID identification and characterization of their conformation ensembles. We conclude with an overview of future computational approaches toward broadening the understanding of intrinsic disorder in the cardiac sarcomere.

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.

Umbrella Sampling Simulations Measure Switch Peptide Binding and Hydrophobic Patch Opening Free Energies in Cardiac Troponin

The cardiac troponin (cTn) complex is an important regulatory protein in heart contraction. Upon binding of Ca2+, cTn undergoes a conformational shift that allows the troponin I switch peptide (cTnISP) to be released from the actin filament and bind to the troponin C hydrophobic patch (cTnCHP). Mutations and modifications to this complex can change its sensitivity to Ca2+ and alter the energetics of the transition from the Ca2+-unbound, cTnISP-unbound form to the Ca2+-bound, cTnISP-bound form. We utilized targeted molecular dynamics (TMD) to obtain a trajectory of this transition pathway, followed by umbrella sampling to estimate the free energy associated with the cTnISP-cTnCHP binding and the cTnCHP opening events for wild-type (WT) cTn. We were able to reproduce experimental values for the cTnISP-cTnCHP binding event and obtain cTnCHP opening free energies in agreement with previous computational measurements of smaller cTnC systems. This excellent agreement for WT cTn demonstrated the strength of computational methods in studying the dynamics and energetics of the cTn complex. We then introduced mutations to the cTn complex that cause cardiomyopathy or alter its Ca2+ sensitivity and observed a general decrease in the free energy of opening the cTnCHP. For these same mutations, we observed no general trend in the effect on the cTnISP-cTnCHP binding event. Our method sets the stage for future computational studies on this system that predict the consequences of yet uncharacterized mutations on cTn dynamics and energetics.

Free-Energy Surfaces of Two Cardiac Thin Filament Conformational Changes during Muscle Contraction

The troponin core is an important regulatory complex in cardiac sarcomeres. Contraction is initiated by a calcium ion binding to cardiac troponin C (cTnC), initiating a conformational shift within the protein, altering its interactions with cardiac troponin I (cTnI). The change in cTnC-cTnI interactions prompts the C-terminal domain of cTnI to dissociate from actin, allowing tropomyosin to reveal myosin-binding sites on actin. Each of the concerted movements in the cardiac thin filament (CTF) is crucial for allowing the contraction of cardiomyocytes, yet little is known about the free energy associated with each transition, which is vital for understanding contraction on a molecular level. Using metadynamics, we calculated the free-energy surface of two transitions in the CTF: cTnC opening in the presence and absence of Ca²⁺ and cTnI dissociating from actin with both open and closed cTnC. These results not only provide the free-energy surface of the transitions but will also be shown to determine if the order of transitions in the contraction cycle is important. From our calculations, we found that the calcium ion helps stabilize the open conformation of cTnC and that the C-terminus of cTnI is stabilized by cTnC in the open conformation when dissociating from the actin surface.

Computational Methods Elucidate Consequences of Mutations and Post-translational Modifications on Troponin I Effective Concentration to Troponin C

Ca²⁺ binding to cardiac troponin C (cTnC) causes a conformational shift that exposes a hydrophobic patch (cTnC_HP) for binding of the cTnI switch peptide (cTnI_SP), ultimately resulting in contraction of the heart. The inhibitory peptide (cTnI_IP), attached at the N-terminal end of the cTnI_SP, serves as a tether for the cTnI_SP to the rest of the troponin complex. Due to this tethered nature, the cTnI_SP remains within proximity of the hydrophobic patch region, resulting in the cTnC_HP experiencing an elevated "effective concentration" of the cTnI_SP. Mutations to the cTnI_IP region have been hypothesized to cause disease by affecting the ability of the cTnI_SP to "find" the hydrophobic patch, resulting in alterations to the heart's ability to contract normally. We tested this hypothesis using molecular dynamics (MD) simulations of the troponin complex using a model that contained all three subunits of troponin: C, I, and T. We developed methods that allowed us to quantitatively measure the effective concentration of the cTnI_SP from the simulations. A significant reduction in the cTnI_SP effective concentration was observed when the cTnI_IP was removed from the system, showcasing the importance of a tethered cTnI_SP. Through accelerated MD methods, we proposed the minimum effective concentration of a tethered cTnI_SP to be approximately 21 mM. Modification of the cTnI_IP via PKC-mediated phosphorylation of T143 was shown to significantly increase the estimated effective concentration of cTnI_SP, help the cTnI_SP find the cTnC_HP more effectively, and maintain the relative shape of the cTnI_IP when compared to the native model. All of these data indicate that pT143 may be able to help promote binding of cTnI_SP to the cTnC_HP. We then tested six mutations within the cTnI_IP region that are known cTnC Ca²⁺-sensitizing mutations and have been linked with cardiomyopathy. We did not observe a significant reduction in the effective concentration upon the introduction of these mutations; however, we did observe increased variability in the flexibility and dynamics of the cTnI_IP region when compared to native. Our observations led us to hypothesize that the mechanism by which these cardiomyopathic mutations affect Ca²⁺ sensitivity is by altering the off rate of cTnI_SP from the hydrophobic patch.

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

Members of the parvalbumin (PV) family of calcium (Ca2+) binding proteins (CBPs) share a relatively high level of sequence similarity. However, their Ca2+ affinities and selectivities against competing ions like Mg2+ can widely vary. We conducted molecular dynamics simulations of several α-parvalbumin (αPV) constructs with micromolar to nanomolar Ca2+ affinities to identify structural and dynamic features that contribute to their binding of ions. Specifically, we examined a D94S/G98E construct with a lower Ca2+ 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 Mg2+ to these isoforms, which is much weaker than Ca2+. 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 Mg2+ from Ca2+, they did not indicate that Ca2+ 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 Ca2+ binding by -20 kcal/mol relative to WT versus 30 kcal/mol for D94S/G98E αPV. Meanwhile, Mg2+ 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, Ca2+ affinity and selectivity against Mg2+ 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.

Molecular dynamics provides new insights into the mechanism of calcium signal transduction and interdomain interactions in cardiac troponin

Understanding the regulation of cardiac muscle contraction at a molecular level is crucial for the development of therapeutics for heart conditions. Despite the availability of atomic structures of the protein components of cardiac muscle thin filaments, detailed insights into their dynamics and response to calcium are yet to be fully depicted. In this study, we used molecular dynamics simulations of the core domains of the cardiac muscle protein troponin to characterize the equilibrium dynamics of its calcium-bound and calcium-free forms, with a focus on elements of cardiac muscle contraction activation and deactivation, that is, calcium binding to the cardiac troponin Ca²⁺ -binding subunit (TnC) and the release of the switch region of the troponin inhibitory subunit (TnI) from TnC. The process of calcium binding to the TnC binding site is described as a three-step process commencing with calcium capture by the binding site residues, followed by cooperative residue interplay bringing the calcium ion to the binding site, and finally, calcium-water exchange. Furthermore, we uncovered a set of TnC-TnI interdomain interactions that are critical for TnC N-lobe hydrophobic pocket dynamics. Absence of these interactions allows the closure of the TnC N-lobe hydrophobic pocket while the TnI switch region remains expelled, whereas if the interactions are maintained, the hydrophobic pocket remains open. Modification of these interactions may fine-tune the ability of the TnC N-lobe hydrophobic pocket to close or remain open, modulate cardiac contractility and present potential therapy-relevant targets.

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.

Anomalous structural dynamics of minimally frustrated residues in cardiac troponin C triggers hypertrophic cardiomyopathy

Cardiac TnC (cTnC) is highly conserved among mammals, and genetic variants can result in disease by perturbing Ca2+-regulation of myocardial contraction. Here, we report the molecular basis of a human mutation in cTnC's αD-helix (TNNC1-p.C84Y) that impacts conformational dynamics of the D/E central-linker and sampling of discrete states in the N-domain, favoring the "primed" state associated with Ca2+ binding. We demonstrate cTnC's αD-helix normally functions as a central hub that controls minimally frustrated interactions, maintaining evolutionarily conserved rigidity of the N-domain. αD-helix perturbation remotely alters conformational dynamics of the N-domain, compromising its structural rigidity. Transgenic mice carrying this cTnC mutation exhibit altered dynamics of sarcomere function and hypertrophic cardiomyopathy. Together, our data suggest that disruption of evolutionary conserved molecular frustration networks by a myofilament protein mutation may ultimately compromise contractile performance and trigger hypertrophic cardiomyopathy.

Discovery of Novel Small-Molecule Calcium Sensitizers for Cardiac Troponin C: A Combined Virtual and Experimental Screening Approach

Heart failure is a leading cause of death throughout the world and is triggered by a disruption of the cardiac contractile machinery. This machinery is regulated in a calcium-dependent manner by the protein complex troponin. Calcium binds to the N-terminal domain of cardiac troponin C (cNTnC) setting into motion the cascade of events leading to muscle contraction. Because of the severity and prevalence of heart failure, there is a strong need to develop small-molecule therapeutics designed to increase the calcium sensitivity of cardiac troponin in order to treat this devastating condition. Molecules that are able to stabilize an open configuration of cNTnC and additionally facilitate the binding of the cardiac troponin I (cTnI) switch peptide have the potential to enable increased calcium sensitization and strengthened cardiac function. Here, we employed a high throughput virtual screening methodology built upon the ability of computational docking to reproduce known experimental results and to accurately recognize cNTnC conformations conducive to small molecule binding using a receiver operator characteristic curve analysis. This approach combined with concurrent stopped-flow kinetic experimental verification led to the identification of a number of sensitizers, which slowed the calcium off-rate. An initial hit, compound 4, was identified with medium affinity (84 ± 30 μM). Through refinement, a calcium sensitizing agent, compound 5, with an apparent affinity of 1.45 ± 0.09 μM was discovered. This molecule is one of the highest affinity calcium sensitizers known to date.

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.

Mechanism of Cardiac Troponin C Calcium Sensitivity Modulation by Small Molecules Illuminated by Umbrella Sampling Simulations

Cardiac troponin C (cTnC) binds intracellular calcium and subsequently cardiac troponin I (cTnI), initiating cardiac muscle contraction. Due to its role in contraction, cTnC has been a therapeutic target in the search for small molecules to treat conditions that interfere with normal muscle contraction like the heritable cardiomyopathies. Structural studies have shown the binding location of small molecules such as bepridil, dfbp-o, 3-methyldiphenylamine (DPA), and W7 to be a hydrophobic pocket in the regulatory domain of cTnC (cNTnC) but have not shown the influence of these small molecules on the energetics of opening this domain. Here we describe an application of an umbrella sampling method used to elucidate the impact these calcium sensitivity modulators have on the free energy of cNTnC hydrophobic patch opening. We found that all these molecules lowered the free energy of opening in the absence of the cTnI, with bepridil facilitating the least endergonic transformation. In the presence of cTnI, however, we saw a stabilization of the open configuration due to DPA and dfbp-o binding, and a destabilization of the open configuration imparted by bepridil and W7. Predicted poor binding molecule NSC34337 left the hydrophobic patch in under 3 ns in conventional MD simulations suggesting that only hydrophobic patch binders stabilized the open conformation. In conclusion, this study presents a novel approach to study the impact of small molecules on hydrophobic patch opening through umbrella sampling, and it proposes mechanisms for calcium sensitivity modulation.

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.

Successful Identification of Cardiac Troponin Calcium Sensitizers Using a Combination of Virtual Screening and ROC Analysis of Known Troponin C Binders

Calcium-dependent cardiac muscle contraction is regulated by the protein complex troponin. Calcium binds to the N-terminal domain of troponin C (cNTnC) which initiates the process of contraction. Heart failure is a consequence of a disruption of this process. With the prevalence of this condition, a strong need exists to find novel compounds to increase the calcium sensitivity of cNTnC. Desirable are small chemical molecules that bind to the interface between cTnC and the cTnI switch peptide and exhibit calcium sensitizing properties by possibly stabilizing cTnC in an open conformation. To identify novel drug candidates, we employed a structure-based drug discovery protocol that incorporated the use of a relaxed complex scheme (RCS). In preparation for the virtual screening, cNTnC conformations were identified based on their ability to correctly predict known cNTnC binders using a receiver operating characteristics analysis. Following a virtual screen of the National Cancer Institute's Developmental Therapeutic Program database, a small number of molecules were experimentally tested using stopped-flow kinetics and steady-state fluorescence titrations. We identified two novel compounds, 3-(4-methoxyphenyl)-6,7-chromanediol (NSC600285) and 3-(4-methylphenyl)-7,8-chromanediol (NSC611817), that show increased calcium sensitivity of cTnC in the presence of the regulatory domain of cTnI. The effects of NSC600285 and NSC611817 on the calcium dissociation rate was stronger than that of the known calcium sensitizer bepridil. Thus, we identified a 3-phenylchromane group as a possible key pharmacophore in the sensitization of cardiac muscle contraction. Building on this finding is of interest to researchers working on development of drugs for calcium sensitization.

Characterization of Zebrafish Cardiac and Slow Skeletal Troponin C Paralogs by MD Simulation and ITC

Zebrafish, as a model for teleost fish, have two paralogous troponin C (TnC) genes that are expressed in the heart differentially in response to temperature acclimation. Upon Ca(2+) binding, TnC changes conformation and exposes a hydrophobic patch that interacts with troponin I and initiates cardiac muscle contraction. Teleost-specific TnC paralogs have not yet been functionally characterized. In this study we have modeled the structures of the paralogs using molecular dynamics simulations at 18°C and 28°C and calculated the different Ca(2+)-binding properties between the teleost cardiac (cTnC or TnC1a) and slow-skeletal (ssTnC or TnC1b) paralogs through potential-of-mean-force calculations. These values are compared with thermodynamic binding properties obtained through isothermal titration calorimetry (ITC). The modeled structures of each of the paralogs are similar at each temperature, with the exception of helix C, which flanks the Ca(2+) binding site; this region is also home to paralog-specific sequence substitutions that we predict have an influence on protein function. The short timescale of the potential-of-mean-force calculation precludes the inclusion of the conformational change on the ΔG of Ca(2+) interaction, whereas the ITC analysis includes the Ca(2+) binding and conformational change of the TnC molecule. ITC analysis has revealed that ssTnC has higher Ca(2+) affinity than cTnC for Ca(2+) overall, whereas each of the paralogs has increased affinity at 28°C compared to 18°C. Microsecond-timescale simulations have calculated that the cTnC paralog transitions from the closed to the open state more readily than the ssTnC paralog, an unfavorable transition that would decrease the ITC-derived Ca(2+) affinity while simultaneously increasing the Ca(2+) sensitivity of the myofilament. We propose that the preferential expression of cTnC at lower temperatures increases myofilament Ca(2+) sensitivity by this mechanism, despite the lower Ca(2+) affinity that we have measured by ITC.

Changes in the dynamics of the cardiac troponin C molecule explain the effects of Ca2+-sensitizing mutations

Cardiac troponin C (cTnC) is the regulatory protein that initiates cardiac contraction in response to Ca²⁺ TnC binding Ca²⁺ initiates a cascade of protein-protein interactions that begins with the opening of the N-terminal domain of cTnC, followed by cTnC binding the troponin I switch peptide (TnI_SW). We have evaluated, through isothermal titration calorimetry and molecular-dynamics simulation, the effect of several clinically relevant mutations (A8V, L29Q, A31S, L48Q, Q50R, and C84Y) on the Ca²⁺ affinity, structural dynamics, and calculated interaction strengths between cTnC and each of Ca²⁺ and TnI_SW Surprisingly the Ca²⁺ affinity measured by isothermal titration calorimetry was only significantly affected by half of these mutations including L48Q, which had a 10-fold higher affinity than WT, and the Q50R and C84Y mutants, each of which had affinities 3-fold higher than wild type. This suggests that Ca²⁺ affinity of the N-terminal domain of cTnC in isolation is insufficient to explain the pathogenicity of these mutations. Molecular-dynamics simulation was used to evaluate the effects of these mutations on Ca²⁺ binding, structural dynamics, and TnI interaction independently. Many of the mutations had a pronounced effect on the balance between the open and closed conformations of the TnC molecule, which provides an indirect mechanism for their pathogenic properties. Our data demonstrate that the structural dynamics of the cTnC molecule are key in determining myofilament Ca²⁺ sensitivity. Our data further suggest that modulation of the structural dynamics is the underlying molecular mechanism for many disease mutations that are far from the regulatory Ca²⁺-binding site of cTnC.

SEEKR: Simulation Enabled Estimation of Kinetic Rates, A Computational Tool to Estimate Molecular Kinetics and Its Application to Trypsin-Benzamidine Binding

We present the Simulation Enabled Estimation of Kinetic Rates (SEEKR) package, a suite of open-source scripts and tools designed to enable researchers to perform multiscale computation of the kinetics of molecular binding, unbinding, and transport using a combination of molecular dynamics, Brownian dynamics, and milestoning theory. To demonstrate its utility, we compute the kon, koff, and ΔGbind for the protein trypsin with its noncovalent binder, benzamidine, and examine the kinetics and other results generated in the context of the new software, and compare our findings to previous studies performed on the same system. We compute a kon estimate of (2.1 ± 0.3) × 107 M-1 s-1, a koff estimate of 83 ± 14 s-1, and a ΔGbind of -7.4 ± 0.1 kcal·mol-1, all of which compare closely to the experimentally measured values of 2.9 × 107 M-1 s-1, 600 ± 300 s-1, and -6.71 ± 0.05 kcal·mol-1, respectively.

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.

Myofilament Calcium Sensitivity: Consequences of the Effective Concentration of Troponin I

Control of calcium binding to and dissociation from cardiac troponin C (TnC) is essential to healthy cardiac muscle contraction/relaxation. There are numerous aberrant post-translational modifications and mutations within a plethora of contractile, and even non-contractile, proteins that appear to imbalance this delicate relationship. The direction and extent of the resulting change in calcium sensitivity is thought to drive the heart toward one type of disease or another. There are a number of molecular mechanisms that may be responsible for the altered calcium binding properties of TnC, potentially the most significant being the ability of the regulatory domain of TnC to bind the switch peptide region of TnI. Considering TnI is essentially tethered to TnC and cannot diffuse away in the absence of calcium, we suggest that the apparent calcium binding properties of TnC are highly dependent upon an “effective concentration” of TnI available to bind TnC. Based on our previous work, TnI peptide binding studies and the calcium binding properties of chimeric TnC-TnI fusion constructs, and building upon the concept of effective concentration, we have developed a mathematical model that can simulate the steady-state and kinetic calcium binding properties of a wide assortment of disease-related and post-translational protein modifications in the isolated troponin complex and reconstituted thin filament. We predict that several TnI and TnT modifications do not alter any of the intrinsic calcium or TnI binding constants of TnC, but rather alter the ability of TnC to “find” TnI in the presence of calcium. These studies demonstrate the apparent consequences of the effective TnI concentration in modulating the calcium binding properties of TnC.

Molecular Basis of S100A1 Activation at Saturating and Subsaturating Calcium Concentrations

The S100A1 protein mediates a wide variety of physiological processes through its binding of calcium (Ca(2+)) and endogenous target proteins. S100A1 presents two Ca(2+)-binding domains: a high-affinity "canonical" EF (cEF) hand and a low-affinity "pseudo" EF (pEF) hand. Accumulating evidence suggests that both Ca(2+)-binding sites must be saturated to stabilize an open state conducive to peptide recognition, yet the pEF hand's low affinity limits Ca(2+) binding at normal physiological concentrations. To understand the molecular basis of Ca(2+) binding and open-state stabilization, we performed 100 ns molecular dynamics simulations of S100A1 in the apo/holo (Ca(2+)-free/bound) states and a half-saturated state, for which only the cEF sites are Ca(2+)-bound. Our simulations indicate that the pattern of oxygen coordination about Ca(2+) in the cEF relative to the pEF site contributes to the former's higher affinity, whereas Ca(2+) binding strongly reshapes the protein's conformational dynamics by disrupting β-sheet coupling between EF hands. Moreover, modeling of the half-saturated configuration suggests that the open state is unstable and reverts toward a closed state in the absence of the pEF Ca(2+) ion. These findings indicate that Ca(2+) binding at the cEF site alone is insufficient to stabilize opening; thus, posttranslational modification of the protein may be required for target peptide binding at subsaturating intracellular Ca(2+) levels.

The Importance of Intrinsically Disordered Segments of Cardiac Troponin in Modulating Function by Phosphorylation and Disease-Causing Mutations

Troponin plays a central role in regulation of muscle contraction. It is the Ca²⁺ switch of striated muscles including the heart and in the cardiac muscle it is physiologically modulated by PKA-dependent phosphorylation at Ser22 and 23. Many cardiomyopathy-related mutations affect Ca²⁺ regulation and/or disrupt the relationship between Ca²⁺ binding and phosphorylation. Unlike the mechanism of heart activation, the modulation of Ca²⁺-sensitivity by phosphorylation of the cardiac specific N-terminal segment of TnI (1–30) is structurally subtle and has proven hard to investigate. The crystal structure of cardiac troponin describes only the relatively stable core of the molecule and the crucial mobile parts of the molecule are missing including TnI C-terminal region, TnI (1–30), TnI (134–149) (“inhibitory” peptide) and the C-terminal 28 amino acids of TnT that are intrinsically disordered. Recent studies have been performed to answer this matter by building structural models of cardiac troponin in phosphorylated and dephosphorylated states based on peptide NMR studies. Now these have been updated by more recent concepts derived from molecular dynamic simulations treating troponin as a dynamic structure. The emerging model confirms the stable core structure of troponin and the mobile structure of the intrinsically disordered segments. We will discuss how we can describe these segments in terms of dynamic transitions between a small number of states, with the probability distributions being altered by phosphorylation and by HCM or DCM-related mutations that can explain how Ca²⁺-sensitivity is modulated by phosphorylation and the effects of mutations.

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.

Cardiac troponin structure-function and the influence of hypertrophic cardiomyopathy associated mutations on modulation of contractility

Cardiac troponin (cTn) acts as a pivotal regulator of muscle contraction and relaxation and is composed of three distinct subunits (cTnC: a highly conserved Ca(2+) binding subunit, cTnI: an actomyosin ATPase inhibitory subunit, and cTnT: a tropomyosin binding subunit). In this mini-review, we briefly summarize the structure-function relationship of cTn and its subunits, its modulation by PKA-mediated phosphorylation of cTnI, and what is known about how these properties are altered by hypertrophic cardiomyopathy (HCM) associated mutations of cTnI. This includes recent work using computational modeling approaches to understand the atomic-based structural level basis of disease-associated mutations. We propose a viewpoint that it is alteration of cTnC-cTnI interaction (rather than the Ca(2+) binding properties of cTn) per se that disrupt the ability of PKA-mediated phosphorylation at cTnI Ser-23/24 to alter contraction and relaxation in at least some HCM-associated mutations. The combination of state of the art biophysical approaches can provide new insight on the structure-function mechanisms of contractile dysfunction resulting cTnI mutations and exciting new avenues for the diagnosis, prevention, and even treatment of heart diseases.

Troponin structure: its modulation by Ca2+and phosphorylation studied by molecular dynamics simulations

The only available crystal structure of the human cardiac troponin molecule (cTn) in the Ca²⁺activated state does not include crucial segments, including the N-terminus of the cTn inhibitory subunit (cTnI).

Effects of HCM cTnI mutation R145G on troponin structure and modulation by PKA phosphorylation elucidated by molecular dynamics simulations

Cardiac troponin (cTn) is a key molecule in the regulation of human cardiac muscle contraction. The N-terminal cardiac-specific peptide of the inhibitory subunit of troponin, cTnI (cTnI(1-39)), is a target for phosphorylation by protein kinase A (PKA) during β-adrenergic stimulation. We recently presented evidence indicating that this peptide interacts with the inhibitory peptide (cTnl(137-147)) when S23 and S24 are phosphorylated. The inhibitory peptide is also the target of the point mutation cTnI-R145G, which is associated with hypertrophic cardiomyopathy (HCM), a disease associated with sudden death in apparently healthy young adults. It has been shown that both phosphorylation and this mutation alter the cTnC-cTnI (C-I) interaction, which plays a crucial role in modulating contractile activation. However, little is known about the molecular-level events underlying this modulation. Here, we computationally investigated the effects of the cTnI-R145G mutation on the dynamics of cTn, cTnC Ca(2+) handling, and the C-I interaction. Comparisons were made with the cTnI-R145G/S23D/S24D phosphomimic mutation, which has been used both experimentally and computationally to study the cTnI N-terminal specific effects of PKA phosphorylation. Additional comparisons between the phosphomimic mutations and the real phosphorylations were made. For this purpose, we ran triplicate 150 ns molecular dynamics simulations of cTnI-R145G Ca(2+)-bound cTnC(1-161)-cTnI(1-172)-cTnT(236-285), cTnI-R145G/S23D/S24D Ca(2+)-bound cTnC(1-161)-cTnI(1-172)-cTnT(236-285), and cTnI-R145G/PS23/PS24 Ca(2+)-bound cTnC(1-161)-cTnI(1-172)-cTnT(236-285), respectively. We found that the cTnI-R145G mutation did not impact the overall dynamics of cTn, but stabilized crucial Ca(2+)-coordinating interactions. However, the phosphomimic mutations increased overall cTn fluctuations and destabilized Ca(2+) coordination. Interestingly, cTnI-R145G blunted the intrasubunit interactions between the cTnI N-terminal extension and the cTnI inhibitory peptide, which have been suggested to play a crucial role in modulating troponin function during β-adrenergic stimulation. These findings offer a molecular-level explanation for how the HCM mutation cTnI-R145G reduces the modulation of cTn by phosphorylation of S23/S24 during β-adrenergic stimulation.

Multiscale Estimation of Binding Kinetics Using Brownian Dynamics, Molecular Dynamics and Milestoning

The kinetic rate constants of binding were estimated for four biochemically relevant molecular systems by a method that uses milestoning theory to combine Brownian dynamics simulations with more detailed molecular dynamics simulations. The rate constants found using this method agreed well with experimentally and theoretically obtained values. We predicted the association rate of a small charged molecule toward both a charged and an uncharged spherical receptor and verified the estimated value with Smoluchowski theory. We also calculated the kon rate constant for superoxide dismutase with its natural substrate, O2-, in a validation of a previous experiment using similar methods but with a number of important improvements. We also calculated the kon for a new system: the N-terminal domain of Troponin C with its natural substrate Ca2+. The kon calculated for the latter two systems closely resemble experimentally obtained values. This novel multiscale approach is computationally cheaper and more parallelizable when compared to other methods of similar accuracy. We anticipate that this methodology will be useful for predicting kinetic rate constants and for understanding the process of binding between a small molecule and a protein receptor.

Ca2+-induced PRE-NMR changes in the troponin complex reveal the possessive nature of the cardiac isoform for its regulatory switch

The interaction between myosin and actin in cardiac muscle, modulated by the calcium (Ca²⁺) sensor Troponin complex (Tn), is a complex process which is yet to be fully resolved at the molecular level. Our understanding of how the binding of Ca²⁺ triggers conformational changes within Tn that are subsequently propagated through the contractile apparatus to initiate muscle activation is hampered by a lack of an atomic structure for the Ca²⁺-free state of the cardiac isoform. We have used paramagnetic relaxation enhancement (PRE)-NMR to obtain a description of the Ca²⁺-free state of cardiac Tn by describing the movement of key regions of the troponin I (cTnI) subunit upon the release of Ca²⁺ from Troponin C (cTnC). Site-directed spin-labeling was used to position paramagnetic spin labels in cTnI and the changes in the interaction between cTnI and cTnC subunits were then mapped by PRE-NMR. The functionally important regions of cTnI targeted in this study included the cTnC-binding N-region (cTnI57), the inhibitory region (cTnI143), and two sites on the regulatory switch region (cTnI151 and cTnI159). Comparison of ¹H-¹⁵N-TROSY spectra of Ca²⁺-bound and free states for the spin labeled cTnC-cTnI binary constructs demonstrated the release and modest movement of the cTnI switch region (∼10 Å) away from the hydrophobic N-lobe of troponin C (cTnC) upon the removal of Ca²⁺. Our data supports a model where the non-bound regulatory switch region of cTnI is highly flexible in the absence of Ca²⁺ but remains in close vicinity to cTnC. We speculate that the close proximity of TnI to TnC in the cardiac complex is favourable for increasing the frequency of collisions between the N-lobe of cTnC and the regulatory switch region, counterbalancing the reduction in collision probability that results from the incomplete opening of the N-lobe of TnC that is unique to the cardiac isoform.

Computational studies of the effect of the S23D/S24D troponin I mutation on cardiac troponin structural dynamics

During β-adrenergic stimulation, cardiac troponin I (cTnI) is phosphorylated by protein kinase A (PKA) at sites S23/S24, located at the N-terminus of cTnI. This phosphorylation has been shown to decrease KCa and pCa50, and weaken the cTnC-cTnI (C-I) interaction. We recently reported that phosphorylation results in an increase in the rate of early, slow phase of relaxation (kREL,slow) and a decrease in its duration (tREL,slow), which speeds up the overall relaxation. However, as the N-terminus of cTnI (residues 1-40) has not been resolved in the whole cardiac troponin (cTn) structure, little is known about the molecular-level behavior within the whole cTn complex upon phosphorylation of the S23/S24 residues of cTnI that results in these changes in function. In this study, we built up the cTn complex structure (including residues cTnC 1-161, cTnI 1-172, and cTnT 236-285) with the N-terminus of cTnI. We performed molecular-dynamics (MD) simulations to elucidate the structural basis of PKA phosphorylation-induced changes in cTn structure and Ca(2+) binding. We found that introducing two phosphomimic mutations into sites S23/S24 had no significant effect on the coordinating residues of Ca(2+) binding site II. However, the overall fluctuation of cTn was increased and the C-I interaction was altered relative to the wild-type model. The most significant changes involved interactions with the N-terminus of cTnI. Interestingly, the phosphomimic mutations led to the formation of intrasubunit interactions between the N-terminus and the inhibitory peptide of cTnI. This may result in altered interactions with cTnC and could explain the increased rate and decreased duration of slow-phase relaxation seen in myofibrils.

Regulatory domain of troponin moves dynamically during activation of cardiac muscle

Heart muscle is activated by Ca²⁺ to generate force and shortening, and the signaling pathway involves allosteric mechanisms in the thin filament. Knowledge about the structure-function relationship among proteins in the thin filament is critical in understanding the physiology and pathology of the cardiac function, but remains obscure. We investigate the conformation of the cardiac troponin (Tn) on the thin filament and its response to Ca²⁺ activation and propose a molecular mechanism for the regulation of cardiac muscle contraction by Tn based uniquely on information from in situ protein domain orientation. Polarized fluorescence from bifunctional rhodamine is used to determine the orientation of the major component of Tn core domain on the thin filaments of cardiac muscle. We show that the C-terminal lobe of TnC (CTnC) does not move during activation, suggesting that CTnC, together with the coiled coil formed by the TnI and TnT chains (IT arm), acts as a scaffold that holds N-terminal lobe of TnC (NTnC) and the actin binding regions of troponin I. The NTnC, on the other hand, exhibits multiple orientations during both diastole and systole. By combining the in situ orientation data with published in vitro measurements of intermolecular distances, we construct a model for the in situ structure of the thin filament. The conformational dynamics of NTnC plays an important role in the regulation of cardiac muscle contraction by moving the C-terminal region of TnI from its actin-binding inhibitory location and enhancing the movement of tropomyosin away from its inhibitory position.

•

In situ conformational changes of troponin in myocardium were investigated.

•

A model for the cardiac thin filament was constructed based on the in situ data.

•

The IT arm of cardiac troponin acts as a scaffold that holds the regulatory domain.

•

The regulatory domain of cardiac troponin moves dynamically during activation.

•

The dynamics of regulatory domain is important in cardiac muscle regulation.

Effects of calcium binding and the hypertrophic cardiomyopathy A8V mutation on the dynamic equilibrium between closed and open conformations of the regulatory N-domain of isolated cardiac troponin C

Troponin C (TnC) is the calcium-binding subunit of the troponin complex responsible for initiating striated muscle contraction in response to calcium influx. In the skeletal TnC isoform, calcium binding induces a structural change in the regulatory N-domain of TnC that involves a transition from a closed to open structural state and accompanying exposure of a large hydrophobic patch for troponin I (TnI) to subsequently bind. However, little is understood about how calcium primes the N-domain of the cardiac isoform (cTnC) for interaction with the TnI subunit as the open conformation of the regulatory domain of cTnC has been observed only in the presence of bound TnI. Here we use paramagnetic relaxation enhancement (PRE) to characterize the closed to open transition of isolated cTnC in solution, a process that cannot be observed by traditional nuclear magnetic resonance methods. Our PRE data from four spin-labeled monocysteine constructs of isolated cTnC reveal that calcium binding triggers movement of the N-domain helices toward an open state. Fitting of the PRE data to a closed to open transition model reveals the presence of a small population of cTnC molecules in the absence of calcium that possess an open conformation, the level of which increases substantially upon Ca(2+) binding. These data support a model in which calcium binding creates a dynamic equilibrium between the closed and open structural states to prime cTnC for interaction with its target peptide. We also used PRE data to assess the structural effects of a familial hypertrophic cardiomyopathy point mutation located within the N-domain of cTnC (A8V). The PRE data show that the Ca(2+) switch mechanism is perturbed by the A8V mutation, resulting in a more open N-domain conformation in both the apo and holo states.