Association of the AB and CD-EF domains from rat alpha- and beta-parvalbumin

What Is Parvalbumin for?

Parvalbumin (PA) is a small, acidic, mostly cytosolic Ca2+-binding protein of the EF-hand superfamily. Structural and physical properties of PA are well studied but recently two highly conserved structural motifs consisting of three amino acids each (clusters I and II), which contribute to the hydrophobic core of the EF-hand domains, have been revealed. Despite several decades of studies, physiological functions of PA are still poorly known. Since no target proteins have been revealed for PA so far, it is believed that PA acts as a slow calcium buffer. Numerous experiments on various muscle systems have shown that PA accelerates the relaxation of fast skeletal muscles. It has been found that oxidation of PA by reactive oxygen species (ROS) is conformation-dependent and one more physiological function of PA in fast muscles could be a protection of these cells from ROS. PA is thought to regulate calcium-dependent metabolic and electric processes within the population of gamma-aminobutyric acid (GABA) neurons. Genetic elimination of PA results in changes in GABAergic synaptic transmission. Mammalian oncomodulin (OM), the β isoform of PA, is expressed mostly in cochlear outer hair cells and in vestibular hair cells. OM knockout mice lose their hearing after 3–4 months. It was suggested that, in sensory cells, OM maintains auditory function, most likely affecting outer hair cells’ motility mechanisms.

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.

The Highly Conservative Cysteine of Oncomodulin as a Feasible Redox Sensor

Oncomodulin (Ocm), or parvalbumin β, is an 11-12 kDa Ca²⁺-binding protein found inside and outside of vertebrate cells, which regulates numerous processes via poorly understood mechanisms. Ocm consists of two active Ca²⁺-specific domains of the EF-hand type ("helix-loop-helix" motif), covered by an EF-hand domain with inactive EF-hand loop, which contains a highly conservative cysteine with unknown function. In this study, we have explored peculiarities of the microenvironment of the conservative Cys18 of recombinant rat Ocm (rWT Ocm), redox properties of this residue, and structural/functional sensitivity of rWT Ocm to the homologous C18S substitution. We have found that pK_a of the Cys18 thiol lays beyond the physiological pH range. The measurement of redox dependence of rWT Ocm thiol-disulfide equilibrium (glutathione redox pair) showed that redox potential of Cys18 for the metal-free and Ca²⁺-loaded protein is of -168 mV and -176 mV, respectively. Therefore, the conservative thiol of rWT Ocm is prone to disulfide dimerization under physiological redox conditions. The C18S substitution drastically reduces α-helices content of the metal-free and Mg²⁺-bound Ocm, increases solvent accessibility of its hydrophobic residues, eliminates the cooperative thermal transition in the apo-protein, suppresses Ca²⁺/Mg²⁺ affinity of the EF site, and accelerates Ca²⁺ dissociation from Ocm. The distinct structural and functional consequences of the minor structural modification of Cys18 indicate its possible redox sensory function. Since some other EF-hand proteins also contain a conservative redox-sensitive cysteine located in an inactive EF-hand loop, it is reasonable to suggest that in the course of evolution, some of the EF-hands attained redox sensitivity at the expense of the loss of their Ca²⁺ affinity.

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.

Site-directed mutagenesis of rat α-parvalbumin: replacement of canonical CD-site residues with their non-consensus counterparts from rat β-parvalbumin

Rat β-parvalbumin (β-PV) displays low divalent-ion affinity. Its CD site is distinguished by six non-consensus residues--the "CD-loop residues"--at positions 49, 50, 57-60. Additionally, leucine occupies position 85, rather than phenylalanine, the β-lineage-consensus residue. Replacement of the CD-loop residues in rat β with the canonical residues was previously found to have little effect on divalent-ion affinity, unless L85 is replaced by phenylalanine. Herein, we replace the canonical CD-loop residues in rat α-PV with their rat β-PV counterparts. Although the mutations have a generally modest impact on affinity, E59D confers Ca(2+)-specificity on the CD site, in the presence or absence of the other mutations. Despite their minimal impact on ΔG, several CD-loop mutations markedly alter ΔH, evidently by perturbing the apo-protein conformation. The L85F mutation was also examined. In wild-type rat α, L85F increases EF-site Ca(2+) affinity. In the CD-loop variants, the mutation leaves the ΔG for Ca(2+)-binding largely unaffected. However, several variants display highly exothermic binding enthalpies, indicative of ligation-linked protein-folding. Consistent with that idea, scanning-calorimetry data confirm that L85F has significantly destabilized those proteins.

Solution and high-pressure NMR studies of the structure, dynamics, and stability of the cross-reactive allergenic cod parvalbumin Gad m 1

Beta-parvalbumins from different fish species have been identified as the main elicitors of IgE-mediated reactions in fish-allergic individuals. Here, we report for the first time the NMR determination of the structure and dynamics of the major Atlantic cod (Gadus morhua) allergen Gad m 1 and compare them with other known parvalbumins. Although the Gad m 1 structure and accessibility of putative IgE epitopes are similar to parvalbumins in mackerel and carp, the charge distribution at the putative epitopes is different. The determination of the Gad m 1 structure contributes to a better understanding of cross-reactivity among fish parvalbumins. In addition, the high-pressure NMR and temperature variation experiments revealed the important contribution of the AB motif and other regions to the protein folding. This structural information could assist the future identification of hot spots for targeted mutations to develop hypoallergenic Ca(2+) -free forms for potential use in immunotherapy.

Sequence microheterogeneity of parvalbumin, the major fish allergen

The microheterogeneity of amino acid sequence observed in various allergens may affect immune response, but incidence of sequence microheterogeneity in allergens and its relation to their allergenicity are unclear. The occurrence of sequence microheterogeneity in major fish allergen, parvalbumin (PA), has been explored using bioinformatics approaches. 44% of 111 species with known PA sequence have PA isoforms. 41% of these species exhibit from 1 to 4 cases of PA sequence microheterogeneity, i.e. unique pairs of PA isoforms with sequence identity above 90%. 29% of 210 PA sequences studied are characterized by microheterogeneity. The occurrence of allergens among them is 2.5-fold higher than among other PAs. The incidence of sequence microheterogeneity within more allergenic β isoform of PA is 2.0-fold lower than that for its less allergenic α isoform. 39 residues affected by PA microheterogeneity are concentrated in the region of helices A, B, F, while helices D and E are the most conservative region. 44% and 11% of the microheterogeneous substitutions are located in the species-specific and cross-reactive IgE-binding epitopes of PAs, respectively. 45% and 48% of the substitution cases are predicted to cause notable changes in protein disorder propensity and protein stability, respectively. Hence, the increased allergenicity rate among PAs having microheterogeneous isoforms can be related to differences in their IgE-binding caused directly or in an allosteric manner. Overall, sequence microheterogeneity is shown to be inherent to many of PAs and likely contributes to PA allergenicity.

Resurrecting prehistoric parvalbumins to explore the evolution of thermal compensation in extant Antarctic fish parvalbumins

Parvalbumins (PVs) from Antarctic notothenioid fishes display a pattern of thermal adaptation that likely reflects evolutionary changes in protein conformational flexibility. We have used ancestral sequence reconstruction and homology modeling to identify two amino acid changes that could potentially account for the present thermal sensitivity pattern of Antarctic fish PVs compared with a PV from a theoretical warm-adapted ancestral fish. To test this hypothesis, ancient PVs were resurrected in the lab using PV from the notothenioid Gobionotothen gibberifrons as a platform for introducing mutations comparable to the reconstructed ancestral PV sequences. The wild-type PV (WT) as well as three mutant expression constructs were engineered: lysine 8 to asparagine (K8N), lysine 26 to asparagine (K26N) and a double mutant (DM). Calcium equilibrium dissociation constants (K(d)) versus temperature curves for all mutants were right-shifted, as predicted, relative to that of WT PV. The K(d) values for the K8N and K26N single mutants were virtually identical at all temperatures and showed an intermediate level of thermal sensitivity. The DM construct displayed a full conversion of thermal sensitivity pattern to that of a PV from a warm/temperate-adapted fish. Additionally, the K(d) versus temperature curve for the WT construct revealed greater thermal sensitivity compared with the mutant constructs. Measurements of the rates of Ca(2+) dissociation (k(off)) showed that all mutants generally had slower k(off) values than WT at all temperatures. Calculated rates of Ca(2+) binding (k(on)) for the K8N and K26N mutants were similar to values for the WT PV at all temperatures. In contrast, the calculated k(on) values for the DM PV were faster, providing mechanistic insights into the nature of potentially adaptive changes in Ca(2+) binding in this PV. The overall results suggest that the current thermal phenotype of Antarctic PVs can be recapitulated by just two amino acid substitutions.

Metal-specific structural changes in parvalbumin

Parvalbumin is a small protein of EF-hand family whose main role is considered to be metal buffering. Recent evidences indicate that parvalbumin also fulfills more complicated functions, which may be determined by the diversity in structural changes in response to the binding of different metal cations. In the present work the conformations of α and β isoforms of pike parvalbumin in the Ca(2+)- and Mg(2+)-loaded state were studied by intrinsic fluorescence, circular dichroism and bis-ANS extrinsic fluorescence. We have determined the structural region causing different spectral response on the binding of Mg(2+)- and Ca(2+) ions in pike β-parvalbumin. Our data reveal similarity of the metal-bound forms of α-parvalbumin. In contrast, those of β isoform differ significantly in the tyrosine spectral range. We also discuss the possible physiological consequences of the structural rearrangements accompanied Mg(2+)/Ca(2+) exchange in pike β-parvalbumin.

Inverse tuning of metal binding affinity and protein stability by altering charged coordination residues in designed calcium binding proteins

Ca(2+ )binding proteins are essential for regulating the role of Ca(2+ )in cell signaling and maintaining Ca(2+ )homeostasis. Negatively charged residues such as Asp and Glu are often found in Ca(2+ )binding proteins and are known to influence Ca(2+ )binding affinity and protein stability. In this paper, we report a systematic investigation of the role of local charge number and type of coordination residues in Ca(2+ )binding and protein stability using de novo designed Ca(2+ )binding proteins. The approach of de novo design was chosen to avoid the complications of cooperative binding and Ca(2+)-induced conformational change associated with natural proteins. We show that when the number of negatively charged coordination residues increased from 2 to 5 in a relatively restricted Ca(2+)-binding site, Ca(2+ )binding affinities increased by more than 3 orders of magnitude and metal selectivity for trivalent Ln(3+ )over divalent Ca(2+ )increased by more than 100-fold. Additionally, the thermal transition temperatures of the apo forms of the designed proteins decreased due to charge repulsion at the Ca(2+ )binding pocket. The thermal stability of the proteins was regained upon Ca(2+ )and Ln(3+ )binding to the designed Ca(2+ )binding pocket. We therefore observe a striking tradeoff between Ca(2+)/Ln(3+ )affinity and protein stability when the net charge of the coordination residues is varied. Our study has strong implications for understanding and predicting Ca(2+)-conferred thermal stabilization of natural Ca(2+ )binding proteins as well as for designing novel metalloproteins with tunable Ca(2+ )and Ln(3+ )binding affinity and selectivity.PACS codes: 05.10.-a.

Metal-controlled interdomain cooperativity in parvalbumins

Conformational behavior of five homologous proteins, parvalbumins (PAs) from northern pike (alpha and beta isoforms), Baltic cod, and rat (alpha and beta isoforms), was studied by scanning calorimetry, circular dichroism, and bis-ANS fluorescence. The mechanism of the temperature-induced denaturation of these proteins depends dramatically on both the peculiarities of their amino acid sequences and on their interaction with metal ions. For example, the pike alpha-PA melting can be described by two successive two-state transitions with mid-temperatures of 90 and 120 degrees C, suggesting the presence of two thermodynamic domains. The intermediate state populated at the end of the first transition was shown to bind Ca(2+) ions, and was characterized by the largely preserved secondary structure and increased solvent exposure of hydrophobic groups. Mg(2+)- and Na(+)-loaded forms of pike alpha-PA demonstrated a single two-state transition. Therefore, the mechanism of the PA thermal denaturation is controlled by metal binding. It ranged from the absence of detectable first-order transition (apo-form of pike PA), to the two-state transition (e.g., Mg(2+)- and Na(+)-loaded forms of pike alpha-PA), to the more complex mechanisms (Ca(2+)-loaded PAs) involving at least one partially folded intermediate. Analysis of isolated cavities in the protein structures revealed that the interface between the CD and EF subdomains of Ca(2+)-loaded pike alpha-PA is much more loosely packed compared with PAs manifesting single heat-sorption peak. The impairment of interactions between CD and EF subdomains may cause a loss of structural cooperativity and appearance of two separate thermodynamic domains. One more peculiar feature of pike alpha-PA is that depending on its interactions with metal ions, it can be an intrinsically disordered protein (apo-form), an ordered protein of mesophilic (Na(+)-bound state), thermophilic (Mg(2+)-form), or even of the hyperthermophilic origin (Ca(2+)-form).

Sequence microheterogeneity of parvalbumin pI 5.0 of pike: a mass spectrometric study

Parvalbumin (PA) is a muscle and neuronal calcium-binding protein, the major fish and frog allergen. Its characteristic feature is the presence of multiple isoforms with significantly different amino acid sequences. Here we show that the major isoform of northern pike muscle PA (pI 5.0, alpha-PA) exhibits microheterogeneity of amino acid sequence. ESI Q-TOF mass-spectrometry (MS) analysis of alpha-PA sample showed the presence of two components with mass difference of 71 Da. Analysis of tryptic and endoproteinase Asp-N digests of alpha-PA by MALDI-TOF MS revealed peptides, corresponding to two different amino acid sequences. The sequence differences between variant proteins are limited to AB-domain and include substitutions K27A and L31K, and an extra Leu residue between K11 and K12. Since the affected residues comprise a cluster on the surface of PA, an involvement of the identified region into target recognition is suggested. The substitutions at positions 27 and 31 are located in the region of previously identified epitopes of parvalbumin relevant for PA-specific IgE and IgG binding, which suggests different immunoactivities of the variants. The found microheterogeneity of PA is suggested to be of importance for physiological adaptation of the propulsive musculature to developmental and/or environmental requirements and may contribute to PA allergenicity.

Solution structure of Ca2+-free rat beta-parvalbumin (oncomodulin)

Relative to other parvalbumin isoforms, the mammalian beta-parvalbumin (oncomodulin) displays attenuated divalent ion affinity. High-resolution structural data for the Ca(2+)-bound protein have provided little insight into the physical basis for this behavior, prompting an examination of the unliganded state. This article describes the solution structure and peptide backbone dynamics of Ca(2+)-free rat beta-parvalbumin (beta-PV). Ca(2+) removal evidently provokes significant structural alterations. Interaction between the D helix and the AB domain in the Ca(2+)-bound protein is greatly diminished in the apo-form, permitting the D helix to straighten. There is also a significant reorganization of the hydrophobic core and a concomitant remodeling of the interface between the AB and CD-EF domains. These modifications perturb the orientation of the C and D helices, and the energetic penalty associated with their reversal could contribute to the low-affinity signature of the CD site. By contrast, Ca(2+) removal causes a comparatively minor perturbation of the E and F helices, consistent with the more typical divalent ion affinity observed for the EF site. Ca(2+)-free rat beta-PV retains structural rigidity on the picosecond-nanosecond timescale. At 20 degrees C, the majority of amide vectors show no evidence for motion on timescales above 20 ps, and the average order parameter for the entire molecule is 0.92.

Reconstitution of calmodulin from domains and subdomains: Influence of target peptide

Reconstitution studies of a protein from domain fragments can furnish important insights into the distinctive role of particular domain interactions and how they affect biophysical properties important for function. Using isothermal titration calorimetry (ITC) and a number of spectroscopic and chromatographic tools, including CD, fluorescence and NMR spectroscopy, size-exclusion chromatography and non-denaturing agarose gel electrophoresis, we have investigated the reconstitution of the ubiquitous Ca2+-sensor protein calmodulin (CaM) and its globular domains from fragments comprising one or two EF-hands. The studies were carried out with and without the target peptide from smooth muscle myosin light chain kinase (smMLCKp). The CaM-target complex can be reconstituted from the three components consisting of the target peptide and the globular domains TR1C and TR2C. In the absence of peptide, there is no evidence for association of the globular domains. The globular domains can further be reconstituted from their corresponding native subdomains. The dissociation constant, K(D), in 2 mM Tris-HCl (pH 7.5), for the subdomain complexes, EF1:EF2 and EF3:EF4, was determined with ITC to 9.3 x 10(-7) M and 5.9 x 10(-8) M, respectively. Thus, the affinity between the two C-terminal subdomains, located within TR2C, is stronger by a factor of 16 than that between the corresponding subdomains within TR1C. These observations are corroborated by the spectroscopic and chromatographic investigations.

Divalent ion-binding properties of the two avian β-parvalbumins

Birds express three parvalbumins, one alpha isoform and two beta isoforms. The latter are known as avian thymic hormone (ATH) and avian parvalbumin 3. Although both were discovered in thymus tissue, and presumably function in T-cell maturation, they have been detected in other tissue settings. We have conducted detailed Ca2+- and Mg2+-binding studies on recombinant ATH and the C72S variant of CPV3, employing global analysis of isothermal titration calorimetry data. In Hepes-buffered saline, ATH binds Ca2+ with apparent microscopic binding constants of 2.4 +/- 0.2 x 10(8) and 1.0 +/- 0.1 x 10(8) M(-1). The corresponding values for CPV3-C72S are substantially lower, 4.5 +/- 0.5 x 10(7) and 2.4 +/- 0.2 x 10(7) M(-1), a 1.9-kcal/mol difference in binding free energy. Thus, the beta-parvalbumin lineage displays a spectrum of Ca2+-binding affinity, with ATH and the mammalian beta isoform at the high- and low-affinity extremes and CPV3 in the middle. Interestingly, despite its decreased Ca2+ affinity, CPV3-C72S exhibits increased affinity for Mg2+, relative to ATH. Whereas the latter displays Mg2+-binding constants of 2.2 +/- 0.2 x 10(4) and 1.2 +/- 0.1 x 10(4) M(-1), CPV3-C72S yields values of 5.0 +/- 0.8 x 10(4) and 2.1 +/- 0.3 x 10(4) M(-1).