Engineering Parvalbumin for the Heart: Optimizing the Mg Binding Properties of Rat β-Parvalbumin

Magnesium Ions Moderate Calcium-Induced Calcium Release in Cardiac Calcium Release Sites by Binding to Ryanodine Receptor Activation and Inhibition Sites

Ryanodine receptor channels at calcium release sites of cardiac myocytes operate on the principle of calcium-induced calcium release. In vitro experiments revealed competition of Ca2+ and Mg2+ in the activation of ryanodine receptors (RyRs) as well as inhibition of RyRs by Mg2+. The impact of RyR modulation by Mg2+ on calcium release is not well understood due to the technical limitations of in situ experiments. We turned instead to an in silico model of a calcium release site (CRS), based on a homotetrameric model of RyR gating with kinetic parameters determined from in vitro measurements. We inspected changes in the activity of the CRS model in response to a random opening of one of 20 realistically distributed RyRs, arising from Ca2+/Mg2+ interactions at RyR channels. Calcium release events (CREs) were simulated at a range of Mg2+-binding parameters at near-physiological Mg2+ and ATP concentrations. Facilitation of Mg2+ binding to the RyR activation site inhibited the formation of sparks and slowed down their activation. Impeding Mg-binding to the RyR activation site enhanced spark formation and speeded up their activation. Varying Mg2+ binding to the RyR inhibition site also dramatically affected calcium release events. Facilitation of Mg2+ binding to the RyR inhibition site reduced the amplitude, relative occurrence, and the time-to-end of sparks, and vice versa. The characteristics of CREs correlated dose-dependently with the effective coupling strength between RyRs, defined as a function of RyR vicinity, single-channel calcium current, and Mg-binding parameters of the RyR channels. These findings postulate the role of Mg2+ in calcium release as a negative modulator of the coupling strength among RyRs in a CRS, translating to damping of the positive feedback of the calcium-induced calcium-release mechanism.

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.

Calcium Buffering in the Heart in Health and Disease

Changes of intracellular Ca²⁺ concentration regulate many aspects of cardiac myocyte function. About 99% of the cytoplasmic calcium in cardiac myocytes is bound to buffers, and their properties will therefore have a major influence on Ca²⁺ signaling. This article considers the fundamental properties and identities of the buffers and how to measure them. It reviews the effects of buffering on the systolic Ca²⁺ transient and how this may change physiologically, and in heart failure and both atrial and ventricular arrhythmias, as well. It is concluded that the consequences of this strong buffering may be more significant than currently appreciated, and a fuller understanding is needed for proper understanding of cardiac calcium cycling and contractility.

Comprehensive analysis of the roles of 'black' and 'gray' clusters in structure and function of rat β-parvalbumin

Recently we found two highly conserved structural motifs in the proteins of the EF-hand calcium binding protein family. These motifs provide a supporting scaffold for the Ca²⁺ binding loops and contribute to the hydrophobic core of the EF-hand domain. Each structural motif forms a cluster of three amino acids called cluster I ('black' cluster) and cluster II ('grey' cluster). Cluster I is much more conserved and mostly incorporates aromatic amino acids. In contrast, cluster II includes a mix of aromatic, hydrophobic, and polar amino acids. The 'black' and 'gray' clusters in rat β-parvalbumin consist of F48, A100, F103 and G61, L64, M87, respectively. In the present work, we sequentially substituted these amino acids residues by Ala, except Ala100, which was substituted by Val. Physical properties of the mutants were studied by circular dichroism, scanning calorimetry, dynamic light scattering, chemical crosslinking, and fluorescent probe methods. The Ca²⁺ and Mg²⁺ binding affinities of these mutants were evaluated by intrinsic fluorescence and equilibrium dialysis methods. In spite of a rather complicated pattern of contributions of separate amino acid residues of the 'black' and 'gray' clusters into maintenance of rat β-parvalbumin structural and functional status, the alanine substitutions in the cluster I cause noticeably more pronounced changes in various structural parameters of proteins, such as hydrodynamic radius of apo-form, thermal stability of Ca²⁺/Mg²⁺-loaded forms, and total energy of Ca²⁺ binding in comparison with the changes caused by amino acid substitutions in the cluster II. These findings were further supported by the outputs of computational analysis of the effects of these mutations on the intrinsic disorder predisposition of rat β-parvalbumin, which also indicated that local intrinsic disorder propensities and the overall levels of predicted disorder were strongly affected by mutations in the cluster I, whereas mutations in cluster II had less pronounced effects. These results demonstrate that amino acids of the cluster I provide more essential contribution to the maintenance of structuraland functional properties of the protein in comparison with the residues of the cluster II.

Effects of Modified Parvalbumin EF-Hand Motifs on Cardiac Myocyte Contractile Function

Cardiac gene delivery of parvalbumin (Parv), an EF-hand Ca(2+) buffer, has been studied as a therapeutic strategy for diastolic heart failure, in which slow Ca(2+) reuptake is an important contributor. A limitation of wild-type (WT) Parv is the significant trade-off between faster relaxation and blunted contraction amplitude, occurring because WT-Parv sequesters Ca(2+) too early in the cardiac cycle and prematurely truncates sarcomere shortening in the facilitation of rapid relaxation. We recently demonstrated that an E → Q substitution (ParvE101Q) at amino acid 12 of the EF-hand Ca(2+)/Mg(2+) binding loop disrupts bidentate Ca(2+) binding, reducing Ca(2+) affinity by 99-fold and increasing Mg(2+) affinity twofold. ParvE101Q caused faster relaxation and not only preserved contractility, but unexpectedly increased it above untreated myocytes. To gain mechanistic insight into the increased contractility, we focused here on amino acid 12 of the EF-hand motif. We introduced an E → D substitution (ParvE101D) at this site, which converts bidentate Ca(2+) coordination to monodentate coordination. ParvE101D decreased Ca(2+) affinity by 114-fold and increased Mg(2+) affinity 28-fold compared to WT-Parv. ParvE101D increased contraction amplitude compared to both untreated myocytes and myocytes with ParvE101Q, with limited improvement in relaxation. Additionally, ParvE101D increased spontaneous contractions after pacing stress. ParvE101D also increased Ca(2+) transient peak height and was diffusely localized around the Z-line of the sarcomere, suggesting a Ca(2+)-dependent mechanism of enhanced contractility. Sarcoplasmic reticulum Ca(2+) load was not changed with ParvE101D, but postpacing Ca(2+) waves were increased. Together, these data show that inverted Ca(2+)/Mg(2+) binding affinities of ParvE101D increase myocyte contractility through a Ca(2+)-dependent mechanism without altering sarcoplasmic reticulum Ca(2+) load and by increasing unstimulated contractions and Ca(2+) waves. ParvE101D provides mechanistic insight into how changes in the Ca(2+)/Mg(2+) binding affinities of parvalbumin's EF-hand motif alter function of cardiac myocytes. These data are informative in developing new Ca(2+) buffering strategies for the failing heart.

Divergent Soybean Calmodulins Respond Similarly to Calcium Transients: Insight into Differential Target Regulation

Plants commonly respond to stressors by modulating the expression of a large family of calcium binding proteins including isoforms of the ubiquitous signaling protein calmodulin (CaM). The various plant CaM isoforms are thought to differentially regulate the activity of specific target proteins to modulate cellular stress responses. The mechanism(s) behind differential target activation by the plant CaMs is unknown. In this study, we used steady-state and stopped-flow fluorescence spectroscopy to investigate the strategy by which two soybean CaMs (sCaM1 and sCaM4) have evolved to differentially regulate NAD kinase (NADK), which is activated by sCaM1 but inhibited by sCaM4. Although the isolated proteins have different cation binding properties, in the presence of Mg2+ and the CaM binding domains from proteins that are differentially regulated, the two plant CaMs respond nearly identically to rapid and slow Ca2+ transients. Our data suggest that the plant CaMs have evolved to bind certain targets with comparable affinities, respond similarly to a particular Ca2+ signature, but achieve different structural states, only one of which can activate the enzyme. Understanding the basis for differential enzyme regulation by the plant CaMs is the first step to engineering a vertebrate CaM that will selectively alter the CaM signaling network.

Interaction with adenylate cyclase toxin from Bordetella pertussis affects the metal binding properties of calmodulin

Adenylate cyclase toxin domain (CyaA‐ACD) is a calmodulin (CaM)‐dependent adenylate cyclase involved in Bordetella pertussis pathogenesis. Calcium (Ca²⁺) and magnesium (Mg²⁺) concentrations impact CaM‐dependent CyaA‐ACD activation, but the structural mechanisms remain unclear. In this study, NMR, dynamic light scattering, and native PAGE were used to probe Mg²⁺‐induced transitions in CaM's conformation in the presence of CyaA‐ACD. Mg²⁺ binding was localized to sites I and II, while sites III and IV remained Ca²⁺ loaded when CaM was bound to CyaA‐ACD. 2Mg²⁺/2Ca²⁺‐loaded CaM/CyaA‐ACD was elongated, whereas mutation of site I altered global complex conformation. These data suggest that CyaA‐ACD interaction moderates CaM's Ca²⁺‐ and Mg²⁺‐binding capabilities, which may contribute to pathobiology.

Understanding Ion Binding Affinity and Selectivity in β-Parvalbumin Using Molecular Dynamics and Mean Spherical Approximation Theory

Parvalbumin (PV) is a globular calcium (Ca(2+))-selective protein expressed in a variety of biological tissues. Our computational studies of the rat β-parvalbumin (β-PV) isoform seek to elucidate the molecular thermodynamics of Ca(2+) versus magnesium (Mg(2+)) binding at the protein's two EF-hand motifs. Specifically, we have utilized molecular dynamics (MD) simulations and a mean-field electrolyte model (mean spherical approximation (MSA) theory) to delineate how the EF-hand scaffold controls the "local" thermodynamics of Ca(2+) binding selectivity over Mg(2+). Our MD simulations provide the probability density of metal-chelating oxygens within the EF-hand scaffolds for both Ca(2+) and Mg(2+), as well the conformational strain induced by Mg(2+) relative to Ca(2+) binding. MSA theory utilizes the binding domain oxygen and charge distributions to predict the chemical potential of ion binding, as well as their corresponding concentrations within the binding domain. We find that the electrostatic and steric contributions toward ion binding were similar for Mg(2+) and Ca(2+), yet the latter was 5.5 kcal/mol lower in enthalpy when internal strain within the EF hand was considered. We therefore speculate that beyond differences in dehydration energies for the Ca(2+) versus Mg(2+), strain induced in the β-PV EF hand by cation binding significantly contributes to the nearly 10,000-fold difference in binding affinity reported in the literature. We further complemented our analyses of local factors governing cation binding selectivity with whole-protein (global) contributions, such as interhelical residue-residue contacts and solvent exposure of hydrophobic surface. These contributions were found to be comparable for both Ca(2+)- and Mg(2+)-bound β-PV, which may implicate local factors, EF-hand strain, and dehydration, in providing the primary means of selectivity. We anticipate these methods could be used to estimate metal binding thermodynamics across a broad range of PV sequence homologues and EF-hand-containing, Ca(2+) binding proteins.

Designing proteins to combat disease: Cardiac troponin C as an example

Throughout history, muscle research has led to numerous scientific breakthroughs that have brought insight to a more general understanding of all biological processes. Potentially one of the most influential discoveries was the role of the second messenger calcium and its myriad of handling and sensing systems that mechanistically control muscle contraction. In this review we will briefly discuss the significance of calcium as a universal second messenger along with some of the most common calcium binding motifs in proteins, focusing on the EF-hand. We will also describe some of our approaches to rationally design calcium binding proteins to palliate, or potentially even cure cardiovascular disease. Considering not all failing hearts have the same etiology, genetic background and co-morbidities, personalized therapies will need to be developed. We predict designer proteins will open doors for unprecedented personalized, and potentially, even generalized medicines as gene therapy or protein delivery techniques come to fruition.

Parvalbumin as a metal-dependent antioxidant

Parvalbumin (PA) is a Ca(2+)-binding protein of vertebrates massively expressed in tissues with high oxygen uptake and respectively elevated level of reactive oxygen species (ROS). To characterize antioxidant properties of PA, antioxidant capacity (AOC) of intact rat α-PA has been explored. ORAC, TEAC and hydrogen peroxide AOC assays evidence conformation-dependent oxidation of the PA. AOC value for the apo-PA 4-11-fold exceeds that for the Ca(2+)-loaded protein. Despite folded conformation of apo-PA, it has AOC equivalent to that of the proteolized protein. The most populated under resting conditions PA form, Mg(2+)-bound PA, has AOC similar to that of apo-PA. ROS-induced changes in absorption spectrum of PA evidence an oxidation of PA's phenylalanines in the ORAC assay. Sensitivity of PA oxidation to its conformation enabled characterization of its metal affinity and pH-dependent behavior: a transition with pKa of 7.6 has been revealed for the Ca(2+)-loaded PA. Since total AOC of PA under in vivo conditions may reach the level of reduced glutathione, we propose that PA might modulate intracellular redox equilibria and/or signaling in a calcium-dependent manner. We speculate that the oxidation-mediated damage of some of PA-GABAergic interneurons observed in schizophrenia is due to a decline in total AOC of the reduced glutathione-PA pair.

Engineered troponin C constructs correct disease-related cardiac myofilament calcium sensitivity

Aberrant myofilament Ca(2+) sensitivity is commonly observed with multiple cardiac diseases, especially familial cardiomyopathies. Although the etiology of the cardiomyopathies remains unclear, improving cardiac muscle Ca(2+) sensitivity through either pharmacological or genetic approaches shows promise of alleviating the disease-related symptoms. Due to its central role as the Ca(2+) sensor for cardiac muscle contraction, troponin C (TnC) stands out as an obvious and versatile target to reset disease-associated myofilament Ca(2+) sensitivity back to normal. To test the hypothesis that aberrant myofilament Ca(2+) sensitivity and its related function can be corrected through rationally engineered TnC constructs, three thin filament protein modifications representing different proteins (troponin I or troponin T), modifications (missense mutation, deletion, or truncation), and disease subtypes (familial or acquired) were studied. A fluorescent TnC was utilized to measure Ca(2+) binding to TnC in the physiologically relevant biochemical model system of reconstituted thin filaments. Consistent with the pathophysiology, the restrictive cardiomyopathy mutation, troponin I R192H, and ischemia-induced truncation of troponin I (residues 1-192) increased the Ca(2+) sensitivity of TnC on the thin filament, whereas the dilated cardiomyopathy mutation, troponin T ΔK210, decreased the Ca(2+) sensitivity of TnC on the thin filament. Rationally engineered TnC constructs corrected the abnormal Ca(2+) sensitivities of the thin filament, reconstituted actomyosin ATPase activity, and force generation in skinned trabeculae. Thus, the present study provides a novel and versatile therapeutic strategy to restore diseased cardiac muscle Ca(2+) sensitivity.