Characterisation of purified parvalbumin from five fish species and nucleotide sequencing of this major allergen from Pacific pilchard, Sardinops sagax

IgE-mediated allergic reaction to seafood is a common cause of food allergy including anaphylactic reactions. Parvalbumin, the major fish allergen, has been shown to display IgE cross-reactivity among fish species consumed predominantly in Europe and the Far East. However, cross-reactivity studies of parvalbumin from fish species widely consumed in the Southern hemisphere are limited as is data relating to immunological and molecular characterisation. In this study, antigenic cross-reactivity and the presence of oligomers and isomers of parvalbumin from five highly consumed fish species in Southern Africa were assessed by immunoblotting using purified parvalbumin and crude fish extracts. Pilchard (Sardinops sagax) parvalbumin was found to display the strongest IgE reactivity among 10 fish-allergic consumers. The cDNA sequence of the beta-form of pilchard parvalbumin was determined and designated Sar sa 1.0101 (accession number FM177701 EMBL/GenBank/DDBJ databases). Oligomeric forms of parvalbumin were observed in all fish species using a monoclonal anti-parvalbumin antibody and subject's sera. Isoforms varied between approximately 10-13 kDa. A highly cross-reactive allergenic isoform of parvalbumin was identified and sequenced, providing a successful primary step towards the generation of a recombinant form that could be used for diagnostic and potential therapeutic use in allergic individuals.

Physicochemical characterization of allergens: quantity, identity, purity, aggregation and conformation

Allergens and allergoids can be characterized by means of physicochemical methods, resulting in a description of the protein on different structural levels. Several techniques are available and their suitability depends on the composition of the particular sample. Current European legislation on allergen products demands characterization of final products in particular focusing on identity, degree of modification (for allergoids) and stability of the composition. Structural parameters of allergens may be used to investigate the stability of an allergen product. The challenge is to identify and optimize techniques that allow determination of protein structure in the context of a formulated pharmaceutical product. As the majority of the products currently marketed are formulated with aluminium hydroxide or aluminium phosphate as a depot, most of the methods and techniques used for protein characterization in solution are not applicable. An additional hurdle is that allergen products are based on natural extracts, comprising a mixture of proteins, both allergens and non-allergens, sometimes in the presence of other uncharacterized components from the raw material. This paper describes which methods are suitable for the different stages of allergen product manufacturing, and how these relate to the current regulatory requirements. Some of the techniques are demonstrated using a model allergen, cod parvalbumin, and a chemically modified form thereof. We conclude that a variety of methods is available for characterization of proteins in solution, and that a limited number of techniques appear to be suitable for modified allergens (allergoids). Adaptation of existing methods, e.g. mass spectroscopy and infrared spectroscopy may be helpful to obtain protein parameters of allergens in a formulated allergen product. By choosing a combination of techniques, including those additional to physicochemical approaches, relevant parameters of allergens in formulated allergen products can be assessed in order to achieve a well-characterized pharmaceutical product.

Combining molecular dynamics and machine learning to improve protein function recognition

As structural genomics efforts succeed in solving protein structures with novel folds, the number of proteins with known structures but unknown functions increases. Although experimental assays can determine the functions of some of these molecules, they can be expensive and time consuming. Computational approaches can assist in identifying potential functions of these molecules. Possible functions can be predicted based on sequence similarity, genomic context, expression patterns, structure similarity, and combinations of these. We investigated whether simulations of protein dynamics can expose functional sites that are not apparent to the structure-based function prediction methods in static crystal structures. Focusing on Ca2+ binding, we coupled a machine learning tool that recognizes functional sites, FEATURE, with Molecular Dynamics (MD) simulations. Treating molecules as dynamic entities can improve the ability of structure-based function prediction methods to annotate possible functional sites.

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.

Parvalbumin isoforms differentially accelerate cardiac myocyte relaxation kinetics in an animal model of diastolic dysfunction

The cytosolic Ca(2+)/Mg(2+)-binding protein alpha-parvalbumin (alpha-Parv) has been shown to accelerate cardiac relaxation; however, beyond an optimal concentration range, alpha-Parv can also diminish contractility. Mathematical modeling suggests that increasing Parv's Mg(2+) affinity may lower the effective concentration of Parv ([Parv]) to speed relaxation and, thus, limit Parv-mediated depressed contraction. Naturally occurring alpha/beta-Parv isoforms show divergence in amino acid primary structure (57% homology) and cation-binding affinities, with beta-Parv having an estimated 16% greater Mg(2+) affinity and approximately 200% greater Ca(2+) affinity than alpha-Parv. We tested the hypothesis that, at the same or lower estimated [Parv], mechanical relaxation rate would be more significantly accelerated by beta-Parv than by alpha-Parv. Dahl salt-sensitive (DS) rats were used as an experimental model of diastolic dysfunction. Relaxation properties were significantly slowed in adult cardiac myocytes isolated from DS rats compared with controls: time from peak contraction to 50% relaxation was 57 +/- 2 vs. 49 +/- 2 (SE) ms (P < 0.05), validating this model system. DS cardiac myocytes were subsequently transduced with alpha- or beta-Parv adenoviral vectors. Upon Parv gene transfer, beta-Parv caused significantly faster relaxation than alpha-Parv (P < 0.05), even though estimated [beta-Parv] was approximately 10% of [alpha-Parv]. This comparative analysis showing distinct functional outcomes raises the prospect of utilizing naturally occurring Parv variants to address disease-associated slowed cardiac relaxation.

Low-affinity signature of the rat beta-parvalbumin CD site. Evidence for remote determinants

Although rat beta-parvalbumin and chicken parvalbumin 3 (CPV3) are identical at 74 of 108 residues, rat beta exhibits perceptibly lower Ca2+ and Mg2+ affinities. At 25 degrees C, in Hepes-buffered saline, at pH 7.4, the overall deltadeltaG degrees ' values are 2.0 and 3.9 kcal/mol, respectively. These differences primarily reflect the disparate behavior of the CD sites in the two proteins. Their respective binding constants for Ca2+, for example, are 1.5 x 10(6) and 2.4 x 10(7) M-1. The extent to which this differential behavior is dictated by local and remote sequence differences is unknown. To explore this question, we performed mutagenesis on rat beta, substituting the corresponding CPV3 codon for residues 49, 50, 57, 58, 59, and 60. The resulting CD site is identical to CPV3 at 27 of 30 positions. The mutations were introduced in four stages, replacing residues 49 and 50 (yielding beta 49/50), then 57 and 58 (beta 49/50/57/58), then 59 (beta 49/50/57/58/59), and finally 60 (beta 49/50/57/58/59/60). Apoprotein stability was examined by scanning calorimetry and chemical denaturation and divalent ion affinity by titration calorimetry. All four variants exhibit elevated Tm values and are between 0.13 and 0.39 kcal/mol more stable at 25 degrees C. Although all four proteins display heightened divalent ion affinity, the increases are small. The maximal deltadeltaG degrees ' values, observed for 49/50/57/58/59/60, are just -0.56 and -0.96 kcal/mol for Ca2+ and Mg2+, respectively. Evidently, structural features beyond the metal ion-binding motif contribute to the unusual divalent ion-binding behavior associated with the rat beta CD site.

Potential influence of Asp in the Ca2+ coordination position 5 of parvalbumin on the calcium-binding affinity: A computational study

Parvalbumins (PV) are calcium-binding proteins, all sharing the common helix-loop-helix (EF-hand) motif. This motif contains a central twelve-residue Ca(2+)-binding loop with the flanking helices positioned roughly perpendicular to each other. The precise role of these coordination residues has been the subject of intense studies. In this work, we focus on the coordination position 5 in the CD Ca(2+)-binding site of silver hake parvalbumin isoform B (SHPV-B). The most common residue at site 5 of calcium-binding loop in canonical EF-hands is Asp [B.J. Marsden, G.S. Shaw, B.D. Sykes, Biochem. Cell Biol. 68 (1990) 587-601], but in the CD site of PV, this position is almost always serine (Ser). The substitution of Ser with Asp will add the 5th carboxylate residue in the CD coordination sphere. However, as predicted by the acid pair hypothesis, the Ca(2+)-binding affinity would be maximized in an EF-hand motif that has four carboxylate ligands paired along the +/-x, and +/-z-axes [R.E. Reid, R.S. Hodges, J. Theor. Biol. 84 (1980) 401-444]. Molecular dynamics simulations and free energy calculations were employed to investigate the influence of Ser to Asp mutation at position 5 on calcium-binding affinity. We found that the Asp variant exhibited remarkable stability during the entire molecular dynamics simulation, with not only the retention of the Ca(2+)-binding site, but also increased compactness in the coordination sphere. The S55D fragment also accommodated Ca(2+) well. We conclude that the reason why Asp which is the most common residue at site 5 of calcium-binding loop in canonical EF-hands has never been identified at this position experimentally for PVs might be related to its physiological functions.

Identification of commercial hake and grenadier species by proteomic analysis of the parvalbumin fraction

Analysis of parvalbumin fractions through proteomic methodologies allowed the differential classification of ten commercial, closely related species of the family Merlucciidae. Muscle extracts from nine hake species of the genus Merluccius including two subspecies of Merluccius australis (australis and polylepsis) and one grenadier species Macruronus novaezelandiae with two populations (novaezelandiae and magellanicus) were evaluated by 2-DE and MALDI-TOF MS. 2-DE demonstrated that the species tested displayed a low intra-specific degree of polymorphism and the isoform patterns were noticeably species-specific. MALDI-TOF mass fingerprints showed clear differences in the pattern of peptides produced by tryptic digestion between the Merluccius and the Macruronus, making the genus differentiation possible. In addition, a selective peptide mass present in the spectra from certain hakes allowed its classification in two groups: Euro-African and American hakes. Besides, some specific masses allowed a clearly individual identification for M. bilinearis, M. australis polylepsis, M. australis australis, M. productus, M. paradoxus and M. polli, while the rest of the hake species can be grouped in two clusters, comprising M. hubbsi and M. gayi in one and M. merluccius and M. capensis in the other.

Detecting fish parvalbumin with commercial mouse monoclonal anti-frog parvalbumin IgG

Parvalbumin is a calcium-binding muscle protein that is highly conserved across fish species and amphibians. It is the major cross-reactive allergen associated with both fish and frog allergy. We used two-dimensional electrophoretic and immunoblotting techniques to investigate the utility of a commercial monoclonal anti-frog parvalbumin IgG for detecting parvalbumin present in some commonly consumed fish species. The 2D electrophoresis and immunoblots revealed species-specific differences in proteins that appear to represent various numbers of isoforms of parvalbumin in carp (5), catfish (3), cod (1) and tilapia (2). No parvalbumin was detected in yellowfin tuna. Based on minor differences in relative intensities of protein staining and immunodetection, parvalbumin isoforms may have slight differences in the epitope region recognized by the anti-frog parvalbumin antibody. These results suggest that the frog anti-parvalbumin antibody can be used as a valuable tool to detect parvalbumins from the fish tested in this study, except yellowfin tuna.

Functional characterization of parvalbumin from the Arctic cod (Boreogadus saida): similarity in calcium affinity among parvalbumins from polar teleosts

Calcium dissociation constants (KD) were measured as a function of temperature for parvalbumin, a small acidic protein expressed abundantly in fast-twitch muscle, from the Arctic cod (Boreogadus saida) and compared to values previously determined for Antarctic and temperate zone teleosts. Estimates of KD were derived independently from fluorometric titrations and calorimetry. In addition, the primary structure of B. saida parvalbumin was determined. Calcium KDs for parvalbumin from B. saida were fundamentally similar to those for parvalbumins from Antarctic species (6.68+/-0.59 nM and 7.77+/-0.72 nM at 5 degrees C, respectively), but significantly different from temperate zone species (1.35+/-0.28 nM at 5 degrees C). However, estimates of KD for B. saida parvalbumin at 5 degrees C closely matched values for temperate zone fish at 25 degrees C (6.54+/-0.56 nM), recapitulating the prior observation that calcium affinity of parvalbumin is conserved at the native temperature of teleost fish. Full sequence of B. saida parvalbumin was generated using reverse-phase HPLC and RACE-PCR. The Arctic parvalbumin showed 83% homology to a carp parvalbumin. None of the 16 total substitutions between the two parvalbumins resided in the cation binding sites of the protein, indicating that the structural locus of the thermal sensitivity of function lies outside the active regions.

Parvalbumin isoforms in zebrafish

By using an analysis of existing genomic information it is concluded that in zebrafish nine genes encode parvalbumin (PV). These genes possess introns that differ in size and show nucleotide variability but they contain the same number of exons, and for each corresponding exon, the number of nucleotides therein are identical in all the paralogs. This rule also applies to the multiple PV genes of other species e.g. mammals. Each of these genes displays, however, characteristic 5' and 3' UTRs which appear highly conserved between closely related species (so that orthologs among these species can be readily identified) but which show larger numbers of mutations between species that are more distant in evolution. A tree is presented which suggests that the traditional classification of PVs as alpha or beta (based mainly on charge of the protein molecule) is not sustainable. Numbers 1-9 are assigned to the various isoforms to facilitate their identification in future studies. A bifurcation of isoforms into 1 and 4; 2 and 3; 6 and 7; 8 and 9 appears to have occurred simultaneously in more recent time, i.e. perhaps approximately 60 mys ago when primates and rodents branched.

Endogenous and exogenous Ca2+ buffers differentially modulate Ca2+-dependent inactivation of Ca(v)2.1 Ca2+ channels

Voltage-gated Ca2+ channels undergo a negative feedback regulation by Ca2+ ions, Ca2+-dependent inactivation, which is important for restricting Ca2+ signals in nerve and muscle. Although the molecular details underlying Ca2+-dependent inactivation have been characterized, little is known about how this process might be modulated in excitable cells. Based on previous findings that Ca2+-dependent inactivation of Ca(v)2.1 (P/Q-type) Ca2+ channels is suppressed by strong cytoplasmic Ca2+ buffering, we investigated how factors that regulate cellular Ca2+ levels affect inactivation of Ca(v)2.1 Ca2+ currents in transfected 293T cells. We found that inactivation of Ca(v)2.1 Ca2+ currents increased exponentially with current amplitude with low intracellular concentrations of the slow buffer EGTA (0.5 mm), but not with high concentrations of the fast Ca2+ buffer BAPTA (10 mm). However, when the concentration of BAPTA was reduced to 0.5 mm, inactivation of Ca2+ currents was significantly greater than with an equivalent concentration of EGTA, indicating the importance of buffer kinetics in modulating Ca2+-dependent inactivation of Ca(v)2.1. Cotransfection of Ca(v)2.1 with the EF-hand Ca2+-binding proteins, parvalbumin and calbindin, significantly altered the relationship between Ca2+ current amplitude and inactivation in ways that were unexpected from behavior as passive Ca2+ buffers. We conclude that Ca2+-dependent inactivation of Ca(v)2.1 depends on a subplasmalemmal Ca2+ microdomain that is affected by the amplitude of the Ca2+ current and differentially modulated by distinct Ca2+ buffers.

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).

Principal Component Analysis of the Conformational Freedom within the EF-Hand Superfamily

A database of nonredundant structures of EF-hand domains--i.e., pairs of helix-loop-helix motifs--has been assembled, and the six angles among the four helices re-determined. A principal component analysis of these angles allows us to use two such components (PC1 and PC2) to describe the system retaining 80% of the total variance. A PC2 against PC1 plot representation allows us to represent in a compact way the full range of structural diversity of EF-hand domains, their grouping into protein families, and the variation for each family upon calcium and peptide binding.

Mimotopes identify conformational epitopes on parvalbumin, the major fish allergen

Parvalbumin, the major fish allergen, is recognized by allergen-specific IgE of more than 90% of all fish-allergic patients. A detailed knowledge of allergenic structures is crucial for developing a vaccine inducing blocking antibodies specifically directed towards the IgE binding epitopes. In the present study we aimed to use the phage display technique to generate mimotopes, which mimic epitopes on parvalbumin. Parvalbumin-specific IgE was purified from sera of fish-allergic patients and used for screening of a constrained decamer phage library. After four rounds of biopanning using parvalbumin-specific IgE, five phage clones were selected which were specifically recognized by parvalbumin-specific IgE as well as IgG. DNA sequencing and peptide alignment revealed a high degree of sequence similarities between the mimotopes. Interestingly, on the surface of natural parvalbumin three regions could be defined by computational mimotope matching. In accordance, previously defined allergenic peptides of cod parvalbumin highlighted areas in close proximity or overlapping with the mimotope matching sites. From the presented data we conclude that our approach identified conformational epitopes of parvalbumin relevant for IgE and IgG binding. We suggest that these mimotopes are suitable candidates for an epitope-specific immunotherapy of fish-allergic patients.

Crystal Structure of the D94S/G98E Variant of Rat α-Parvalbumin. An Explanation for the Reduced Divalent Ion Affinity

Simultaneous replacement of Asp-94 with serine and Gly-98 with glutamate in rat alpha-parvalbumin creates a CD-site ligand array in the context of the EF-site binding loop. Previous work has shown that, relative to the wild-type CD site, this engineered site has markedly reduced Ca(2+) affinity. Seeking an explanation for this phenomenon, we have obtained the crystal structure of the alpha D94S/G98E variant. The Ca(2+) coordination within the engineered EF site of the 94/98E variant is nearly identical to that within the CD site, suggesting that the attenuated affinity of the EF site in 94/98E is not a consequence of suboptimal coordination geometry. We have also examined the divalent ion binding properties of the alpha 94/98E variant in both Na(+)- and K(+)-containing buffers. Although the Ca(2+) and Mg(2+) affinities are higher in K(+) solution, the increases are comparable to those observed for wild-type alpha. Consistent with that finding, the apparent Na(+) stoichiometry, estimated from stability studies conducted as a function of Na(+) concentration, is 1.0 +/- 0.1, identical to that of wild-type alpha. Thus, the reduced affinity for divalent ions is evidently not the result of heightened monovalent ion competition. The thermodynamic analysis indicates that the less favorable Gibbs free energy of binding reflects a substantial enthalpic penalty. Significantly, the crystal structure reveals a steric clash between Phe-57 and the C(gamma) atom of Glu-98. The consequent displacement of Phe-57 also produces a close contact with Ser-55. Thus, steric interference may be the source of the enthalpic penalty.

A united residue force-field for calcium-protein interactions

United-residue potentials are derived for interactions of the calcium cation with polypeptide chains in energy-based prediction of protein structure with a united-residue (UNRES) force-field. Specific potentials were derived for the interaction of the calcium cation with the Asp, Glu, Asn, and Gln side chains and the peptide group. The analytical expressions for the interaction energies for each of these amino acids were obtained by averaging the electrostatic interaction energy, expressed by a multipole series over the dihedral angles not considered in the united-residue model, that is, the side-chain dihedral angles chi and the dihedral angles lambda for the rotation of peptide groups about the C(alpha)...C(alpha) virtual-bond axes. For the side-chains that do not interact favorably with calcium, simple excluded-volume potentials were introduced. The parameters of the potentials were obtained from ab initio quantum mechanical calculations of model systems at the Restricted Hartree-Fock (RHF) level with the 6-31G(d,p) basis set. The energy surfaces of pairs consisting of Ca(2+)-acetate, Ca(2+)-propionate, Ca(2+)-acetamide, Ca(2+)-propionamide, and Ca(2+)-N-methylacetamide systems (modeling the Ca(2+)-Asp(-), Ca(2+)-Glu(-), Ca(2+)-Asn, Ca(2+)-Gln, and Ca(2+)-peptide group interactions) at different distances and orientations were calculated. For each pair, the restricted free energy (RFE) surfaces were calculated by numerical integration over the degrees of freedom lost when switching from the all-atom model to the united-residue model. Finally, the analytical expressions for each pair were fitted to the RFE surfaces. This force-field was able to distinguish the EF-hand motif from all potential binding sites in the crystal structures of bovine alpha-lactalbumin, whiting parvalbumin, calbindin D9K, and apo-calbindin D9K.

Solution structure of human beta-parvalbumin and structural comparison with its paralog alpha-parvalbumin and with their rat orthologs

The aim of this research was to determine the structure of human beta-parvalbumin (109 amino acids) and to compare it with its paralog and ortholog proteins. The structure was determined in solution using multinuclear and multidimensional NMR methods and refined using substitution of the EF-hand Ca(2+) ion with a paramagnetic lanthanide. The resulting family of structures had a backbone rmsd of 0.50 A. Comparison with rat oncomodulin (X-ray, 1.3 A resolution) as well as with human (NMR, backbone rmsd of 0.49 A) and rat (X-ray, 2.0 A resolution) parvalbumins reveals small but reliable local differences, often but not always related to amino acid variability. The analysis of these structures has led us to propose an explanation for the different affinity for Ca(2+) between alpha- and beta-parvalbumins and between parvalbumins and calmodulins.

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

Association of the parvalbumin AB and CD-EF domains was examined in Hepes-buffered saline, pH 7.4, employing fragments from rat alpha and beta. All of the interactions require Ca(2+). In saturating Ca(2+), the alpha AB/alpha CD-EF (alpha/alpha) complex displays an association constant of (7.6 +/- 0.4) x 10(7) M(-1). Ca(2+)-binding data for a mixture of the alpha fragments are compatible with an identical two-site model, yielding an average binding constant of (8.5 +/- 0.2) x 10(5) M(-1). The beta/beta interaction is significantly weaker, exhibiting an association constant of (3.0 +/- 0.6) x 10(6) M(-1). The Ca(2+)-binding constants for beta/beta are likewise diminished, at (1.0 +/- 0.1) x 10(5) and (2.3 +/- 0.2) x 10(4) M(-1). The magnitude of the apparent DeltaDeltaG(degree)' for Ca(2+) binding by alpha/alpha and beta/beta, at 3.4 kcal/mol, approaches that measured for the intact proteins (3.6 kcal/mol) and is substantially larger than the 1.5 kcal/mol value previously measured for the isolated CD-EF domains. This result suggests that the AB domain can modulate the Ca(2+) affinities of the CD and EF sites. Interestingly, the heterologous alpha/beta complex displays a larger association constant [(6.6 +/- 0.4) x 10(6) M(-1)] than the homologous beta/beta complex and heightened Ca(2+) affinity [binding constants of (1.3 +/- 0.1) x 10(6) and (8.8 +/- 0.2) x 10(4) M(-1)]. By contrast, beta/alpha associates more weakly than alpha/alpha and exhibits sharply reduced affinity for Ca(2+). Thus, the interaction between the beta AB domain and beta CD-EF domain may act to attenuate Ca(2+) affinity in the intact protein.

Crystal structure of a high-affinity variant of rat alpha-parvalbumin

In model peptide systems, Ca2+ affinity is maximized in EF-hand motifs containing four carboxylates positioned on the +x and -x and +z and -z axes; introduction of a fifth carboxylate ligand reduces the affinity. However, in rat beta-parvalbumin, replacement of Ser-55 with aspartate heightens divalent ion affinity [Henzl, M. T., et al. (1996) Biochemistry 35, 5856-5869]. The corresponding alpha-parvalbumin variant (S55D/E59D) likewise exhibits elevated affinity [Henzl, M. T., et al. (2003) Anal. Biochem. 319, 216-233]. To determine whether these mutations produce a variation on the archetypal EF-hand coordination scheme, we have obtained high-resolution X-ray crystallographic data for alpha S55D/E59D. As anticipated, the aspartyl carboxylate replaces the serine hydroxyl at the +z coordination position. Interestingly, the Asp-59 carboxylate abandons the role it plays as an outer sphere ligand in wild-type rat beta, rotating away from the Ca2+ and, instead, forming a hydrogen bond with the amide of Glu-62. Superficially, the coordination sphere in the CD site of alpha S55D/E59D resembles that in the EF site. However, the orientation of the Asp-59 side chain is predicted to stabilize the D-helix, which may contribute to the heightened divalent ion affinity. DSC data indicate that the alpha S55D/E59D variant retains the capacity to bind 1 equiv of Na+. Consistent with this finding, when binding measurements are conducted in K(+)-containing buffer, divalent ion affinity is markedly higher. In 0.15 M KCl and 0.025 M Hepes-KOH (pH 7.4) at 5 degrees C, the macroscopic Ca2+ binding constants are 1.8 x 10(10) and 2.0 x 10(9) M(-1). The corresponding Mg2+ binding constants are 2.7 x 10(6) and 1.2 x 10(5) M(-1).

Rat α- and β-Parvalbumins: Comparison of Their Pentacarboxylate and Site-Interconversion Variants

Introduction of a fifth carboxylate into the ligand array of the CD site (via the combined S55D and E59D mutations) or the EF site (G98D) of rat alpha-parvalbumin substantially increases divalent ion affinity. This behavior, in conflict with that seen in model peptide systems, agrees with existing data for rat beta-parvalbumin [Henzl et al. (1996) Biochemistry 35, 5856-5869]. The complete analysis of the S55D/E59D double variant necessitated characterization of alpha E59D. Whereas the D59E mutation has minimal influence on beta CD site affinity, E59D has a major impact on the alpha CD site, lowering the apparent association constant by a factor of 14. The thermodynamic consequences of exchanging the rat alpha CD and EF site ligand arrays, which differ at the +z and -x coordination positions, were also examined. When the alpha CD array is imported into the EF site, it acquires a low-affinity phenotype, in agreement with previous findings for beta [Henzl et al. (1998) Biochemistry 37, 9101-9111]. However, when the EF ligand array is introduced into the alpha CD binding loop, it retains a high-affinity signature. This result, contrary to that observed in beta, suggests that the influence of the parvalbumin CD site environment supersedes the intrinsic behavior of the ligand array, a conclusion further supported by the disparate impact of the beta D59E and alpha E59D mutations.