Thermodynamics and molecular dynamics simulations of calcium binding to the regulatory site of human cardiac troponin C: evidence for communication with the structural calcium binding sites

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.

Binding of calcium and magnesium to human cardiac troponin C

Cardiac muscle thin filaments are composed of actin, tropomyosin, and troponin that change conformation in response to Ca²⁺ binding, triggering muscle contraction. Human cardiac troponin C (cTnC) is the Ca²⁺-sensing component of the thin filament. It contains structural sites (III/IV) that bind both Ca²⁺ and Mg²⁺ and a regulatory site (II) that has been thought to bind only Ca²⁺. Binding of Ca²⁺ at this site initiates a series of conformational changes that culminate in force production. However, the mechanisms that underpin the regulation of binding at site II remain unclear. Here, we have quantified the interaction between site II and Ca²⁺/Mg²⁺ through isothermal titration calorimetry and thermodynamic integration simulations. Direct and competitive binding titrations with WT N-terminal cTnC and full-length cTnC indicate that physiologically relevant concentrations of both Ca²⁺/Mg²⁺ interacted with the same locus. Moreover, the D67A/D73A N-terminal cTnC construct in which two coordinating residues within site II were removed was found to have significantly reduced affinity for both cations. In addition, 1 mM Mg²⁺ caused a 1.4-fold lower affinity for Ca²⁺. These experiments strongly suggest that cytosolic-free Mg²⁺ occupies a significant population of the available site II. Interaction of Mg²⁺ with site II of cTnC likely has important functional consequences for the heart both at baseline as well as in diseased states that decrease or increase the availability of Mg²⁺, such as secondary hyperparathyroidism or ischemia, respectively.

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.

In depth, thermodynamic analysis of Ca2+ binding to human cardiac troponin C: Extracting buffer-independent binding parameters

BACKGROUND

Characterizing the thermodynamic parameters behind metal-biomolecule interactions is fundamental to understanding the roles metal ions play in biology. Isothermal Titration Calorimetry (ITC) is a "gold-standard" for obtaining these data. However, in addition to metal-protein binding, additional equilibria such as metal-buffer interactions must be taken into consideration prior to making meaningful comparisons between metal-binding systems.

METHODS

In this study, the thermodynamics of Ca²⁺ binding to three buffers (Bis-Tris, MES, and MOPS) were obtained from Ca²⁺-EDTA titrations using ITC. These data were used to extract buffer-independent parameters for Ca²⁺ binding to human cardiac troponin C (hcTnC), an EF-hand containing protein required for heart muscle contraction.

RESULTS

The number of protons released upon Ca²⁺ binding to the C- and N-domain of hcTnC were found to be 1.1 and 1.2, respectively. These values permitted determination of buffer-independent thermodynamic parameters of Ca²⁺-hcTnC binding, and the extracted data agreed well among the buffers tested. Both buffer and pH-adjusted parameters were determined for Ca²⁺ binding to the N-domain of hcTnC and revealed that Ca²⁺ binding under aqueous conditions and physiological ionic strength is both thermodynamically favorable and driven by entropy.

CONCLUSIONS

Taken together, the consistency of these data between buffer systems and the similarity between theoretical and experimental proton release is indicative of the reliability of the method used and the importance of extracting metal-buffer interactions in these studies.

GENERAL SIGNIFICANCE

The experimental approach described herein is clearly applicable to other metal ions and other EF-hand protein systems.

The missing links within troponin

The cardiac contraction-relaxation cycle is controlled by a sophisticated set of machinery. Of particular interest, is the revelation that allosteric networks transmit effects of binding at one site to influence troponin complex dynamics and structural-mediated signaling in often distal, functional sites in the myofilament. Our recent observations provide compelling evidence that allostery can explain the function of large-scale macromolecular events. Here we elaborate on our recent findings of interdomain communication within troponin C, using cutting-edge structural biology approaches, and highlight the importance of unveiling the unknown, distant communication networks within this system to obtain more comprehensive knowledge of how allostery impacts cardiac physiology and disease.

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.

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.

Amide hydrogens reveal a temperature-dependent structural transition that enhances site-II Ca2+-binding affinity in a C-domain mutant of cardiac troponin C

The hypertrophic cardiomyopathy-associated mutant D145E, in cardiac troponin C (cTnC) C-domain, causes generalised instability at multiple sites in the isolated protein. As a result, structure and function of the mutant are more susceptible to higher temperatures. Above 25 °C there are large, progressive increases in N-domain Ca²⁺-binding affinity for D145E but only small changes for the wild-type protein. NMR-derived backbone amide temperature coefficients for many residues show a sharp transition above 30–40 °C, indicating a temperature-dependent conformational change that is most prominent around the mutated EF-hand IV, as well as throughout the C-domain. Smaller, isolated changes occur in the N-domain. Cardiac skinned fibres reconstituted with D145E are more sensitive to Ca²⁺ than fibres reconstituted with wild-type, and this defect is amplified near body-temperature. We speculate that the D145E mutation destabilises the native conformation of EF-hand IV, leading to a transient unfolding and dissociation of helix H that becomes more prominent at higher temperatures. This creates exposed hydrophobic surfaces that may be capable of binding unnaturally to a variety of targets, possibly including the N-domain of cTnC when it is in its open Ca²⁺-saturated state. This would constitute a potential route for propagating signals from one end of TnC to the other.

Dissecting ITC data of metal ions binding to ligands and proteins

BACKGROUND

ITC is a powerful technique that can reliably assess the thermodynamic underpinnings of a wide range of binding events. When metal ions are involved, complications arise in evaluating the data due to unavoidable solution chemistry that includes metal speciation and a variety of linked equilibria.

SCOPE OF REVIEW

This paper identifies these concerns, provides recommendations to avoid common mistakes, and guides the reader through the mathematical treatment of ITC data to arrive at a set of thermodynamic state functions that describe identical chemical events and, ideally, are independent of solution conditions. Further, common metal chromophores used in biological metal sensing studies are proposed as a robust system to determine unknown solution competition.

MAJOR CONCLUSIONS

Metal ions present several complications in ITC experiments. This review presents strategies to avoid these pitfalls and proposes and experimentally validates mathematical approaches to deconvolute complex equilibria that exist in these systems.

GENERAL SIGNIFICANCE

This review discusses the wide range of complications that exists in metal-based ITC experiments. It provides a starting point for scientists new to this field and articulates concerns that will help experienced researchers troubleshoot experiments.

The exchanged EF-hands in calmodulin and troponin C chimeras impair the Ca²⁺-induced hydrophobicity and alter the interaction with Orai1: a spectroscopic, thermodynamic and kinetic study

Background

Calmodulin (CaM) plays an important role in Ca²⁺-dependent signal transduction. Ca²⁺ binding to CaM triggers a conformational change, forming a hydrophobic patch that is important for target protein recognition. CaM regulates a Ca²⁺-dependent inactivation process in store-operated Ca²⁺ entry, by interacting Orai1. To understand the relationship between Ca²⁺-induced hydrophobicity and CaM/Orai interaction, chimera proteins constructed by exchanging EF-hands of CaM with those of Troponin C (TnC) are used as an informative probe to better understand the functionality of each EF-hand.

Results

ANS was used to assess the context of the induced hydrophobic surface on CaM and chimeras upon Ca²⁺ binding. The exchanged EF-hands from TnC to CaM resulted in reduced hydrophobicity compared with wild-type CaM. ANS lifetime measurements indicated that there are two types of ANS molecules with rather distinct fluorescence lifetimes, each specifically corresponding to one lobe of CaM or chimeras. Thermodynamic studies indicated the interaction between CaM and a 24-residue peptide corresponding to the CaM-binding domain of Orail1 (Orai-CMBD) is a 1:2 CaM/Orai-CMBD binding, in which each peptide binding yields a similar enthalpy change (ΔH = −5.02 ± 0.13 kcal/mol) and binding affinity (K_a = 8.92 ± 1.03 × 10⁵ M⁻¹). With the exchanged EF1 and EF2, the resulting chimeras noted as CaM(1TnC) and CaM(2TnC), displayed a two sequential binding mode with a one-order weaker binding affinity and lower ΔH than that of CaM, while CaM(3TnC) and CaM(4TnC) had similar binding thermodynamics as CaM. The dissociation rate constant for CaM/Orai-CMBD was determined to be 1.41 ± 0.08 s⁻¹ by rapid kinetics. Stern-Volmer plots of Orai-CMBD Trp76 indicated that the residue is located in a very hydrophobic environment but becomes more solvent accessible when EF1 and EF2 were exchanged.

Conclusions

Using ANS dye to assess induced hydrophobicity showed that exchanging EFs for all Ca²⁺-bound chimeras impaired ANS fluorescence and/or binding affinity, consistent with general concepts about the inadequacy of hydrophobic exposure for chimeras. However, such ANS responses exhibited no correlation with the ability to interact with Orai-CMBD. Here, the model of 1:2 binding stoichiometry of CaM/Orai-CMBD established in solution supports the already published crystal structure.

Type 1 and Type 2 scenarios in hydrogen exchange mass spectrometry studies on protein–ligand complexes

Ligand binding to a protein can elicit a wide range of responses when studied by HDX mass spectrometry.