Crystal structure of the unique parvalbumin component from muscle of the leopard shark (Triakis semifasciata). The first X-ray study of an alpha-parvalbumin

Comprehensive Sequence Analysis of Parvalbumins in Fish and Their Comparison with Parvalbumins in Tetrapod Species

Parvalbumins are small molecules with important functions in Ca²⁺ signaling, but their sequence comparisons to date, especially in fish, have been relatively poor. We here, characterize sequence motifs that distinguish parvalbumin subfamilies across vertebrate species, as well as those that distinguish individual parvalbumins (orthologues) in fish, and map them to known parvalbumin structures. As already observed by others, all classes of jawed vertebrates possess parvalbumins of both the α-parvalbumin and oncomodulin subfamilies. However, we could not find convincing phylogenetic support for the common habit of classifying all non-α-parvalbumins together as "β-parvalbumins." In teleost (modern bony) fish, we here distinguish parvalbumins 1-to-10, of which the gene copy number can differ between species. The genes for α-parvalbumins (pvalb6 and pvalb7) and oncomodulins (pvalb8 and pvalb9) are well conserved between teleost species, but considerable variation is observed in their copy numbers of the non-α/non-oncomodulin genes pvalb1-to-5 and pvalb10. Teleost parvalbumins 1-to-4 are hardly distinguishable from each other and are highly expressed in muscle, and described allergens belong to this subfamily. However, in some fish species α-parvalbumin expression is also high in muscle. Pvalb5 and pvalb10 molecules form distinct lineages, the latter even predating the origin of teleosts, but have been lost in some teleost species. The present study aspires to be a frame of reference for future studies trying to compare different parvalbumins.

Purification, Characterization, and Crystal Structure of Parvalbumins, the Major Allergens in Mustelus griseus

Fish play important roles in human nutrition and health, but also trigger allergic reactions in some population. Parvalbumin (PV) represents the major allergen of fish. While IgE cross-reactivity to PV in various bony fish species has been well characterized, little information is available about allergens in cartilaginous fish. In this study, two shark PV isoforms (named as SPV-I and SPV-II) from Mustelus griseus were purified. Their identities were further confirmed by mass spectroscopic analysis. IgE immunoblot analysis showed that sera from fish-allergic patients reacted to both SPV-I and SPV-II, but the majority of sera reacted more intensely to SPV-I than SPV-II. Thermal denaturation monitored by CD spectrum showed that both of the SPV allergens are highly thermostable. SPV-I maintained its IgE-binding capability after heat denaturation, while the IgE-binding capability of SPV-II was reduced. The results of crystal structure showed that SPV-I and SPV-II were similar in their overall tertiary structure, but their amino acid sequences shared lower similarities, indicating that the differences in the IgE-binding capabilities of SPV-I and SPV-II might be due to differential antigen epitopes in these two isoforms.

Seafood allergy: A comprehensive review of fish and shellfish allergens

Seafood refers to several distinct groups of edible aquatic animals including fish, crustacean, and mollusc. The two invertebrate groups of crustacean and mollusc are, for culinary reasons, often combined as shellfish but belong to two very different phyla. The evolutionary and taxonomic diversity of the various consumed seafood species poses a challenge in the identification and characterisation of the major and minor allergens critical for reliable diagnostics and therapeutic treatments. Many allergenic proteins are very different between these groups; however, some pan-allergens, including parvalbumin, tropomyosin and arginine kinase, seem to induce immunological and clinical cross-reactivity. This extensive review details the advances in the bio-molecular characterisation of 20 allergenic proteins within the three distinct seafood groups; fish, crustacean and molluscs. Furthermore, the structural and biochemical properties of the major allergens are described to highlight the immunological and subsequent clinical cross-reactivities. A comprehensive list of purified and recombinant allergens is provided, and the applications of component-resolved diagnostics and current therapeutic developments are discussed.

Allergenicity of bony and cartilaginous fish - molecular and immunological properties

Allergy to bony fish is common and probably increasing world-wide. The major heat-stable pan-fish allergen, parvalbumin (PV), has been identified and characterized for numerous fish species. In contrast, there are very few reports of allergic reactions to cartilaginous fish despite widespread consumption. The molecular basis for this seemingly low clinical cross-reactivity between these two fish groups has not been elucidated. PV consists of two distinct protein lineages, α and β. The α-lineage of this protein is predominant in muscle tissue of cartilaginous fish (Chondrichthyes), while β-PV is abundant in muscle tissue of bony fish (Osteichthyes). The low incidence of allergic reactions to ingested rays and sharks is likely due to the lack of molecular similarity, resulting in reduced immunological cross-reactivity between the two PV lineages. Structurally and physiologically, both protein lineages are very similar; however, the amino acid homology is very low with 47-54%. Furthermore, PV from ancient fish species such as the coelacanth demonstrates 62% sequence homology to leopard shark α-PV and 70% to carp β-PV. This indicates the extent of conservation of the PV isoforms lineages across millennia. This review highlights prevalence data on fish allergy and sensitization to fish, and details the molecular diversity of the two protein lineages of the major fish allergen PV among different fish groups, emphasizing the immunological and clinical differences in allergenicity.

Immunological cross-reactivity between four distant parvalbumins-Impact on allergen detection and diagnostics

Fish are the largest and most diverse group of vertebrates. Fish are also a part of the eight food groups that cause the majority of IgE mediated food reactions. Detection tools for fish allergens are however limited due to the great diversity of fish species, despite fish allergy and its major allergen parvalbumin being well documented. The most commonly studied fish are frequently consumed in North America and Europe. However, much less is known about fish allergens in the Australasian region although fish is widely consumed in this region. A comprehensive phylogenetic analysis was performed of known parvalbumin amino acid sequences to determine possible candidate antigens for new cross-reactive antibodies to be used to detect most fish parvalbumins. Polyclonal rabbit antibodies were raised against parvalbumins from frequently consumed barramundi (Lates calcarifer), basa (Pangasius bocourti), pilchard (Sardinops sagax) and Atlantic salmon (Salmo salar). These were evaluated for cross-reactivity against a panel of 45 fish extracts (raw, heated and canned fish). Anti-barramundi parvalbumin proved to be the most cross-reactive antibody, detecting 87.5% of the 40 species analyzed, followed by anti-pilchard and anti-basa antibody. In contrast the anti-salmon antibody was very specific and only reacted to salmonidae and a few other fish. All analyzed fish species, except mahi mahi, swordfish, yellowfin tuna and all 5 canned fish had parvalbumin detected in raw extracts. However antibody reactivity to many fish was heat liable or susceptible to denaturation, demonstrating that some parvalbumins have most likely conformational epitopes, which lose antibody reactivity after heat treatment. We have demonstrated the generation of highly cross-reactive anti-parvalbumin antibodies that could be used for the detection of allergenic fish parvalbumin in contaminated food products. This cross-reactivity study thus shows processing of fish, especially canning, can have on impact on antibody recognition by ELISA, possibly similar to IgE-binding in vivo.

Coordination to divalent cations by calcium-binding proteins studied by FTIR spectroscopy

We review the Fourier-transform infrared (FTIR) spectroscopy of side-chain COO(-) groups of Ca(2+)-binding proteins: parvalbumins, bovine calmodulin, akazara scallop troponin C and related calcium binding proteins and peptide analogues. The COO(-) stretching vibration modes can be used to identify the coordination modes of COO(-) groups of Ca(2+)-binding proteins to metal ions: bidentate, unidentate, and pseudo-bridging. FTIR spectroscopy demonstrates that the coordination structure of Mg(2+) is distinctly different from that of Ca(2+) in the Ca(2+)-binding site in solution. The interpretation of COO(-) stretches is ensured on the basis of the spectra of calcium-binding peptide analogues. The implication of COO(-) stretches is discussed for Ca(2+)-binding proteins. This article is part of a Special Issue entitled: FTIR in membrane proteins and peptide studies.

Purification and characterization of parvalbumins, the major allergens in red stingray (Dasyatis akajei)

Fish has received increasing attention because it induces IgE-mediated food allergy. Parvalbumin (PV) represents the major allergen of fish, and IgE cross-reactivity to PV in various teleost fish species has been shown, while little information is available about allergens in elasmobranch fish. In this study, two PV isoforms (named as PV-I and PV-II) from red stingray (Dasyatis akajei) were purified to homogeneity by a series of procedures including ammonium sulfate precipitation and column chromatographies of DEAE-Sepharose and Sephacryl S-200. Purified PVs revealed a single band on tricine-sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The molecular masses of PV-I and PV-II were 12.29 and 11.95 kDa, respectively, as determined by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Western blot using antifrog PV monoclonal antibody (PARV-19) showed positive reactions to the two proteins, confirming that they were PVs, although their immunological reactivities were weaker than those of PV from silver carp. The N-terminal amino acid sequence of PV-I was determined, and comparison with PVs from other fish species showed low homology between teleost and elasmobranch fish. The isoelectric points of PV-I and PV-II were 5.4 and 5.0, respectively, as determined by two-dimensional electrophoresis (2-DE), suggesting that both isoforms belong to the α-group. IgE immunoblotting analysis showed that sera from fish-allergic patients reacted to both PV-I and PV-II from red stingray. Thermal stability revealed that PV-I easily formed oligomers than PV-II, which might contribute to the maintenance of its allerginicity during heat processing.

Structure of avian thymic hormone, a high-affinity avian beta-parvalbumin, in the Ca2+-free and Ca2+-bound states

Originally isolated on the basis of its capacity to stimulate T-cell maturation and proliferation, avian thymic hormone (ATH) is nevertheless a parvalbumin, one of two beta-lineage isoforms expressed in birds. We recently learned that addition of Ca(2+)-free ATH to a solution of 8-anilinonaphthalene-1-sulfonate (ANS) markedly increases ANS emission. This behavior, not observed in the presence of Ca(2+), suggests that apolar surface area buried in the Ca(2+)-bound state becomes solvent accessible upon Ca(2+) removal. In order to elucidate the conformational alterations that accompany Ca(2+) binding, we have obtained the solution structure of the Ca(2+)-free protein using NMR spectroscopy and compared it to the Ca(2+)-loaded protein, solved by X-ray crystallography. Although the metal-ion-binding (CD-EF) domains are largely coincident in the superimposed structures, a major difference is observed in the AB domains. The tight association of helix B with the E and F helices in the Ca(2+)-bound state is lost upon removal of Ca(2+), producing a deep hydrophobic cavity. The B helix also undergoes substantial rotation, exposing the side chains of F24, Y26, F29, and F30 to solvent. Presumably, the increase in ANS emission observed in the presence of unliganded ATH reflects the interaction of these hydrophobic residues with the fluorescent probe. The increased solvent exposure of apolar surface area in the Ca(2+)-free protein is consistent with previously collected scanning calorimetry data, which indicated an unusually low change in heat capacity upon thermal denaturation. The Ca(2+)-free structure also provides added insight into the magnitude of ligation-linked conformational alteration compatible with a high-affinity metal-ion-binding signature. The exposure of substantial apolar surface area suggests the intriguing possibility that ATH could function as a reverse Ca(2+) sensor.

Expression levels of parvalbumins determine allergenicity of fish species

BACKGROUND

Parvalbumins are the most important fish allergens. Polysensitization to various fish species is frequently reported and linked to the cross-reactivity of their parvalbumins. Studies on cross-reactivity and its association to the allergenicity of purified natural parvalbumins from different fish species are still lacking. In addition, some studies indicate that dark muscled fish such as tuna are less allergenic.

METHODS

Total protein extracts and purified parvalbumins from cod, whiff, and swordfish, all eaten frequently in Spain, were tested for their IgE-binding properties with 16 fish allergic patients' sera from Madrid. The extent of cross-reactivity of these parvalbumins was investigated by IgE ELISA inhibition assays. Additionally, the cDNA sequences of whiff and swordfish parvalbumins were determined.

RESULTS

Extractable amounts of parvalbumins from cod were 20 times and from whiff 30 times higher than from swordfish. Parvalbumins were recognized by 94% of the patients in extracts of cod and whiff, but only by 60% in swordfish extracts. Nevertheless, a high cross-reactivity was determined for all purified parvalbumins by IgE inhibition. The amino acid sequence identities of the three parvalbumins were in a range of 62-74%.

CONCLUSIONS

The parvalbumins of cod, whiff and swordfish are highly cross-reactive. The high amino acid sequence identity among cod, whiff and swordfish parvalbumins results in the observed IgE cross-reactivity. The low allergenicity of swordfish is due to the low expression levels of its parvalbumin.

Identification of low-molecular-weight protein (SCP1) from shark cartilage with anti-angiogenesis activity and sequence similarity to parvalbumin

Cartilage was considered as a possible natural source of anti-angiogenesis compounds due to its known avascular nature. In this study, a low-molecular-weight protein with an anti-angiogenesis activity was isolated from shark cartilage using a mild extraction procedure. The protein was purified to homogeneity by gel filtration and electroelution techniques and its N-terminal amino acid sequence was determined. The purified protein, designated as SCP1, represented a molecular weight of 13.7 kDa, pI of 6.9-7 and its N-terminal sequence revealed sequence similarity to alpha parvalbumin family. The protein inhibited angiogenesis when subjected to microvessel sprouting of collagen-embedded rat aortic ring assay. It is suggested that SCP1 could be considered as a new angiogenesis inhibitor derived from shark cartilage.

Molecular cloning, expression and phylogenetic analyses of parvalbumin in tilapia, Oreochromis mossambicus

The gene expression of parvalbumin (Pvalb), a high-affinity calcium-binding protein and the major fish allergen, was significantly increased in the tilapia fry treated with methyltestosterone (MT) as examined using a subtractive hybridization assay. Using the real-time quantitative PCR, we further confirmed the increased Pvalb expression in the MT-treated tilapia fry. The 568 base pairs (bp) tilapia Pvalb (tPvalb) cDNA clone was fully sequenced and found to contain a coding region of 330 bp, which encodes a 108 amino acids protein with a molecular weight of 11,370.5 and an calculated isoelectric point of 4.56. The predicted secondary structure of tPvalb is comprised of seven alpha helices. It contains two characteristic EF-hand calcium-binding motifs, one PKC and five casein kinase II consensus phosphorylation sites. The tPvalb is highly homologous to the selected fish Pvalbs at a similarity ranging from 53% to 80%. The phylogenetic tree analysis showed that the tPvalb is closest to the Scomber japonicus Pvalb. The tPvalb was found to express in the heart, muscle, gill, kidney, brain and ovary of adult fish by RT-PCR analysis. In situ hybridization also revealed that the tPvalb was highly expressed in the hypothalamus and sarcoplasmic reticulum. A tPvalb glutathione S-transferase (GST) fusion protein was generated and digested by thrombin to remove the GST moiety. Further Western analysis showed that the tPvalb protein was cross-reacted to an anti-rat Pvalb antibody. Those results suggest that Pvalb is evolutionally conserved in tilapia.

Calcium-binding proteins and their role in allergic diseases

Calcium-binding proteins (CBPs) are ubiquitous pollen allergens and important food allergens in fish and amphibians. Calcium-binding allergens containing two EF-hands (polcalcins) have been detected and characterized in pollen from trees, grasses, and weeds. Timothy grass Phl p 7 is the most cross-reactive allergen among polcalcins. Although there is cross-reactivity described within the subfamilies of calcium-binding allergens, there are no strong indications for IgE cross-reactivity between CBPs from plants, fish, and humans. Therefore, Phl p 7 could be used as marker to identify multiple pollen-sensitized patients, whereas cod Gad c 1 or carp Cyp c 1 could be selected for the diagnosis of fish allergy. Hom s 4, a calcium-binding autoantigen, might be an interesting candidate to monitor chronic skin inflammation in atopic and nonatopic individuals. Diagnostic tests containing these molecules could allow the identification of most patients sensitized to calcium-binding allergens/antigens. In general, IgE recognition of calcium-binding allergens is influenced by binding or release of calcium ions. This knowledge could be used to engineer hypoallergenic CBPs for specific immunotherapy.

Exploring structurally conserved solvent sites in protein families

Protein-bound water molecules are important components of protein structure, and therefore, protein function and energetics. Although structural conservation of solvent has been studied in a few protein families, a lack of suitable computational tools has hindered more comprehensive analyses. Herein we present a semiautomated computational approach for identifying solvent sites that are conserved among proteins sharing a common three-dimensional structure. This method is tested on six protein families: (1) monodomain cytochrome c, (2) fatty-acid binding protein, (3) lactate/malate dehydrogenase, (4) parvalbumin, (5) phospholipase A2, and (6) serine protease. For each family, the method successfully identified previously known conserved solvent sites. Moreover, the method discovered 22 novel conserved solvent sites, some of which have higher degrees of conservation than the previously known sites. All six families studied had solvent sites with more than 90% conservation and these sites were invariably located in regions of the protein with very high sequence conservation. These results suggest that highly conserved solvent sites, by virtue of their proximity to conserved residues, should be considered as one of the defining three-dimensional structural characteristics of protein families and folds.

Crystal structure of rat alpha-parvalbumin at 1.05 Angstrom resolution

The crystal structure of rat alpha-parvalbumin has been determined at 1.05 Angstrom resolution, using synchrotron data collected at Advanced Photon Source beamline 19-ID. After refinement with SHELX, employing anisotropic displacement parameters and riding hydrogen atoms, R = 0.132 and R(free) = 0.162. The average coordinate estimated standard deviations are 0.021 Angstrom and 0.038 Angstrom for backbone atoms and side-chain atoms, respectively. Besides providing a more precise view of the alpha-isoform than previously available, these data permit comparison with the 0.91 Angstrom structure determined for pike beta-parvalbumin. Visualization of the anisotropic displacement parameters as thermal ellipsoids yields insight into the atomic motion within the Ca(2+)-binding sites. The asymmetric unit includes three parvalbumin (PV) molecules. Interestingly, the EF site in one displays uncharacteristic flexibility. The ellipsoids for Asp-92 are particularly large and non-spherical, and the shape of the Ca(2+) ellipsoid implies significant vibrational motion perpendicular to the plane defined by the four y and z ligands. The relative dearth of crystal-packing interactions in this site suggests that the heightened flexibility may be the result of diminished intermolecular contacts. The implication is that, by impeding conformational mobility, crystal-packing forces may cause serious overestimation of EF-hand rigidity. The high quality of the data permitted 11 residues to be modeled in alternative side-chain conformations, including the two core residues, Ile-97 and Leu-105. The discrete disorder observed for Ile-97 may have functional ramifications, providing a mechanism for communicating binding status between the CD and EF binding loops and between the PV metal ion-binding domain and the N-terminal AB region.

Structure of rat parvalbumin with deleted AB domain: implications for the evolution of EF hand calcium-binding proteins and possible physiological relevance

Among the EF-hand Ca(2+)-binding proteins, parvalbumin (PV) and calbindin D9k (CaB) have the function of Ca(2+) buffers. They evolved from an ancestor protein through two phylogenetic pathways, keeping one pair of EF-hands. They differ by the extra helix-loop-helix (AB domain) found in PV and by the linker between the binding sites. To investigate whether the deletion of AB in PV restores a CaB-like structure, we prepared and solved the structure of the truncated rat PV (PVratDelta37) by X-ray and NMR. PVratDelta37 keeps the PV fold, but is more compact, having a well-structured linker, which differs remarkably from CaB. PvratDelta37 has no stable apo-form, has lower affinity for Ca(2+) than full-length PV, and does not bind Mg(2+), in contrast to CaB. Structural differences of the hydrophobic core are partially responsible for lowering the calcium-binding affinity of the truncated protein. It can be concluded that the AB domain, like the linker of CaB, plays a role in structural stabilization. The AB domain of PV protects the hydrophobic core, and is required to maintain high affinity for divalent cation binding. Therefore, the AB domain possibly modulates PV buffer function.

X-Ray crystal structure and molecular dynamics simulations of silver hake parvalbumin (Isoform B)

Parvalbumins constitute a class of calcium-binding proteins characterized by the presence of several helix-loop-helix (EF-hand) motifs. In a previous study (Revett SP, King G, Shabanowitz J, Hunt DF, Hartman KL, Laue TM, Nelson DJ, 1997, Protein Sci 7:2397-2408), we presented the sequence of the major parvalbumin isoform from the silver hake (Merluccius bilinearis) and presented spectroscopic and structural information on the excised "EF-hand" portion of the protein. In this study, the X-ray crystal structure of the silver hake major parvalbumin has been determined to high resolution, in the frozen state, using the molecular replacement method with the carp parvalbumin structure as a starting model. The crystals are orthorhombic, space group C2221, with a = 75.7 A, b = 80.7 A, and c = 42.1 A. Data were collected from a single crystal grown in 15% glycerol, which served as a cryoprotectant for flash freezing at -188 degrees C. The structure refined to a conventional R-value of 21% (free R 25%) for observed reflections in the range 8 to 1.65 A [1 > 2sigma(I)]. The refined model includes an acetylated amino terminus, 108 residues (characteristic of a beta parvalbumin lineage), 2 calcium ions, and 114 water molecules per protein molecule. The resulting structure was used in molecular dynamics (MD) simulations focused primarily on the dynamics of the ligands coordinating the Ca2+ ions in the CD and EF sites. MD simulations were performed on both the fully Ca2+ loaded protein and on a Ca2+ deficient variant, with Ca2+ only in the CD site. There was substantial agreement between the MD and X-ray results in addressing the issue of mobility of key residues in the calcium-binding sites, especially with regard to the side chain of Ser55 in the CD site and Asp92 in the EF site.

Crystal structure of the EF-hand parvalbumin at atomic resolution (0.91 A) and at low temperature (100 K). Evidence for conformational multistates within the hydrophobic core

Several crystal structures of parvalbumin (Parv), a typical EF-hand protein, have been reported so far for different species with the best resolution achieving 1.5 A. Using a crystal grown under microgravity conditions, cryotechniques (100 K), and synchrotron radiation, it has now been possible to determine the crystal structure of the fully Ca2+-loaded form of pike (component pI 4.10) Parv.Ca2 at atomic resolution (0.91 A). The availability of such a high quality structure offers the opportunity to contribute to the definition of the validation tools useful for the refinement of protein crystal structures determined to lower resolution. Besides a better definition of most of the elements in the protein three-dimensional structure than in previous studies, the high accuracy thus achieved allows the detection of well-defined alternate conformations, which are observed for 16 residues out of 107 in total. Among them, six occupy an internal position within the hydrophobic core and converge toward two small buried cavities with a total volume of about 60 A3. There is no indication of any water molecule present in these cavities. It is probable that at temperatures of physiological conditions there is a dynamic interconversion between these alternate conformations in an energy-barrier dependent manner. Such motions for which the amplitudes are provided by the present study will be associated with a time-dependent remodeling of the void internal space as part of a slow dynamics regime (millisecond timescales) of the parvalbumin molecule. The relevance of such internal dynamics to function is discussed.

Hydration-coupled dynamics in proteins studied by neutron scattering and NMR: the case of the typical EF-hand calcium-binding parvalbumin

The influence of hydration on the internal dynamics of a typical EF-hand calciprotein, parvalbumin, was investigated by incoherent quasi-elastic neutron scattering (IQNS) and solid-state 13C-NMR spectroscopy using the powdered protein at different hydration levels. Both approaches establish an increase in protein dynamics upon progressive hydration above a threshold that only corresponds to partial coverage of the protein surface by the water molecules. Selective motions are apparent by NMR in the 10-ns time scale at the level of the polar lysyl side chains (externally located), as well as of more internally located side chains (from Ala and Ile), whereas IQNS monitors diffusive motions of hydrogen atoms in the protein at time scales up to 20 ps. Hydration-induced dynamics at the level of the abundant lysyl residues mainly involve the ammonium extremity of the side chain, as shown by NMR. The combined results suggest that peripheral water-protein interactions influence the protein dynamics in a global manner. There is a progressive induction of mobility at increasing hydration from the periphery toward the protein interior. This study gives a microscopic view of the structural and dynamic events following the hydration of a globular protein.

Finding local structural similarities among families of unrelated protein structures: a generic non-linear alignment algorithm

We have developed a generic tool for the automatic identification of regions of local structural similarity in unrelated proteins having different folds, as well as for defining more global similarities that result from homologous protein structures. The computer program GENFIT has evolved from the genetic algorithm-based three-dimensional protein structure comparison program GA_FIT. GENFIT, however, can locate and superimpose regions of local structural homology regardless of their position in a pair of structures, the fold topology, or the chain direction. Furthermore, it is possible to restrict the search to a volume centered about a region of interest (e.g., catalytic site, ligand-binding site) in two protein structures. We present a number of examples to illustrate the function of the program, which is a parallel processing implementation designed for distribution to multiple machines over a local network or to run on a single multiprocessor computer.

Conformational stabilities of the rat alpha- and beta-parvalbumins

It is widely believed that beta-parvalbumin (PV) isoforms are intrinsically less stable than alpha-parvalbumins, due to greater electrostatic repulsion and an abbreviated C-terminal helix. However, when examined by differential scanning calorimetry, the apo-form of the rat beta-PV (i.e. oncomodulin) actually displays greater thermal stability than the alpha-PV. Whereas the melting temperature of the a isoform is 45.8 degrees C at physiological pH and ionic strength, the Tm for the beta isoform is more than 7 degrees higher (53.6 degrees C). This result suggests that factors besides net charge and C-terminal helix length strongly influence parvalbumin conformational stability. Extension of the F helix in the beta-PV, by insertion of Ser-109, has a modest stabilizing effect, raising the Tm, by 1.1 degrees. Truncation of the alpha-PV F helix, by removal of Glu-108, has a more profound impact, lowering the Tm by 4.0 degrees.

Ca2+/Mg2+ exchange in parvalbumin and other EF-hand proteins. A theoretical study

A remarkable conformational rearrangement occurs upon Ca2+/Mg2+ exchange in the C-terminal EF-hand site (labelled site EF or EF-4) of parvalbumin, as initially established by X-ray crystallography. Such a conformational rearrangement is characterised as follows: (i) the co-ordination number decreases from seven oxygen atoms in the Ca-loaded form to six oxygen atoms in the Mg-loaded form, the heptaco-ordination of Ca2+ corresponding with a skewed pentagonal bipyramid configuration of the seven oxygen atoms, whereas the hexaco-ordination of Mg2+ corresponds with a regular octahedral configuration of the six oxygen atoms; and (ii) Glu101, at the relative position 12 in the EF-hand loop sequence (labelled "Glu12"), acts as a bidentate ligand in the Ca-loaded form and as a monodentate ligand in the Mg-loaded form. As part of the conformational rearrangement, the chi1 dihedral angle undergoes a gauche(+) to gauche(-) transition upon substitution of Ca2+ by Mg2+, whereas the chi2 angle remains practically unchanged and the chi3 angles in both forms adopt a nearly mirror image relationship. In order to understand the molecular mechanisms underlying such a conformational rearrangement, we undertook a theoretical study using the free energy perturbation (FEP) method, starting from high-resolution crystal structures of the same parvalbumin (pike 4. 10 isoform) differing by the substitution of their two cationic sites EF-3 (or CD) and EF-4 (or EF), i.e. the 1pal structure with EF-3(Ca2+) and EF-4(Ca2+), the 4pal structure with EF-3(Ca2+) and EF-4(Mg2+). When Mg2+ is "alchemically" transformed into Ca2+ within the EF-4 site of 4pal, the conformational rearrangement of Glu12 is correctly predicted by the FEP calculation. When Ca2+ is transformed into Mg2+ within the EF-3 site of 4pal, the FEP calculation predicts the topology of the fully Mg-loaded form for which no crystallographic data is presently available. As expected, Glu62 (at the relative position 12 in EF-3 loop) is predicted to be a monodentate residue within a regular octahedral arrangement of six oxygen atoms around Mg2+. We also investigated the behaviour during Ca2+/Mg2+ exchange of two other typical EF-hand proteins, troponin C (TnC) and calmodulin (CaM), for which no three-dimensional structure of their Mg-loaded forms is available so far. It is also predicted that the EF-3 site of TnC and the EF-1 site of CaM have their invariant Glu12 residues switching from the bidentate to the monodentate configuration when Ca2+ is substituted by Mg2+, with six oxygen atoms being observed in the co-ordination sphere of the alchemically generated Mg2+ cation.

Mode of binding of anti-P-glycoprotein antibody MRK-16 to its antigen. A crystallographic and molecular modeling study

Monoclonal antibody MRK-16 recognizes a discontinuous extracellular epitope on the multidrug resistance-associated ATP-binding cassette transporter, P-glycoprotein. The atomic basis for specificity of this antibody is of interest because of its potential as a modulator of P-glycoprotein activity. The crystal structure of Fab MRK-16 is reported to a resolution of 2.8 A. A structure for a portion of the epitope was derived by comparison to regions of solved structures with similar primary sequence. This has permitted a proposal for the mode of binding of the peptide epitope to the antibody, in which the peptide makes specific contacts with complementarity-determining regions H1, H2, and H3 from the heavy chain and L3 from the light chain. These interactions are consistent with epitope mapping studies and with the observation that MRK-16 is specific for human class I P-glycoprotein. This result identifies side chains in MRK-16 that would be amenable to alteration in antibody engineering experiments to derive improved multidrug resistance inhibitors for clinical use during chemotherapy. In particular, Arg-H97 contacts both Glu-746 and Asp-744 of the peptide, Arg-L96 contacts Asp-743, and Thr-H33 interacts with Thr-747. All of these epitope residues were implicated in mediating specificity by epitope mapping studies.

The cation-binding domain from the alpha subunit of integrin alpha5 beta1 is a minimal domain for fibronectin recognition

The cation-binding domain from the alpha subunit of human integrin alpha5beta1 was produced as a recombinant protein, alpha5-(229-448). This protein displays a well defined fold with a content of 30-35% alpha-helix and 20-25% beta-strand, based on circular dichroism. The binding of Ca2+ or Mg2+ to alpha5-(229-448) results in a biphasic conformational rearrangement consistent with the occurrence of two classes of cation-binding sites differing by their affinities. The two classes of sites are located in two conformationally independent lobes, as established by a parallel study of two recombinant half-domains (N- and C-terminal) that also adopt stable folds. Upon saturation with divalent cations, alpha5-(229-448) binds an Arg-Gly-Asp (RGD)-containing fibronectin ligand to form a 1:1 complex. Complex formation is associated with a specific conformational adaptation of the ligand, suggesting an induced fit mechanism. In contrast, neither of the half-domains is competent for ligand binding. The alpha5-(229-448)-fibronectin complex is dissociated in the presence of an RGD peptide, as well as of a simple carboxylic acid, suggesting that the RGD aspartyl carboxylate is an essential element that directly interacts with the alpha5 cation-binding domain.

15N NMR relaxation studies of calcium-loaded parvalbumin show tight dynamics compared to those of other EF-hand proteins

Dynamics of the rat alpha-parvalbumin calcium-loaded form have been determined by measurement of 15N nuclear relaxation using proton-detected heteronuclear NMR spectroscopy. The relaxation data were analyzed using spectral density functions and the Lipari-Szabo formalism. The major dynamic features for the rat alpha-parvalbumin calcium-loaded form are (1) the extreme rigidity of the helix-loop-helix EF-hand motifs and the linker segment connecting them, (2) the N and C termini of the protein being restricted in their mobility, (3) a conformational exchange occurring at the kink of helix D, and (4) the residue at relative position 2 in the Ca2+-binding sites having an enhanced mobility. Comparison of the Ca2+-binding EF-hand domains of alpha-parvalbumin-Ca2+, calbindin-Ca2+, and calmodulin-Ca2+ shows that parvalbumin is probably the most rigid of the EF-hand proteins. It also illustrates the dynamical properties which are conserved in the EF-hand domains from different members of this superfamily: (1) a tendency toward higher mobility of NH vectors at relative position 2 in the Ca2+-binding loop, (2) a restricted mobility for the other residues in the binding loop, and (3) an overall rigidity for the helices of EF-hand motifs. The differences in mobility between parvalbumin and the two EF-hand proteins occur mainly at the linker connecting the pair of EF hands and also at the C terminus of the last helix. In parvalbumin-Ca2+, these two regions are characterized by a pronounced rigidity compared to the corresponding more mobile regions in calbindin-Ca2+ and calmodulin-Ca2+.

Ca2+ coordination to backbone carbonyl oxygen atoms in calmodulin and other EF-hand proteins: 15N chemical shifts as probes for monitoring individual-site Ca2+ coordination

Examination of the NMR 15N chemical shifts of a number of EF-hand proteins shows that the shift value for the amido nitrogen of the residue in position 8 of a canonical EF-hand loop (or position 10 of a pseudo EF-hand loop) provides a good indication of metal occupation of that site. The NH of the residue in position 8 is covalently bonded to the carbonyl of residue 7, the only backbone carbonyl that coordinates to the metal ion in a canonical EF-hand loop. Upon metal coordination to this carbonyl, there is an appreciable deshielding of the 15N nucleus at position 8 (+4 to +8 ppm) due to the polarization of the O(7)=C(7)-N(8) amido group and the corresponding reduction in the electron density of the nitrogen atom. This deshielding effect is effectively independent of the binding of metal to the other site of an EF-hand pair, allowing the 15N shifts to be used as probes for site-specific occupancy of metal binding sites. In addition, a Ca2+-induced change in side-chain Halpha-Calpha-Cbeta-Hbeta torsion angle for isoleucine or valine residues in position 8 can also contribute to the deshielding of the amide 15N nucleus. This conformational effect occurs only in sites I or III and takes place upon binding a Ca2+ ion to the other site of an EF-hand pair (site II or IV) regardless of whether the first site is occupied. The magnitude of this effect is in the range +5 to +7 ppm. A Ca2+ titration of 15N-labeled apo-calmodulin was performed using 2D 1H-15N HSQC NMR spectra. The changes in the 15N chemical shifts and intensities for the peaks corresponding to the NH groups of residues in position 8 of the EF-hand loops allowed the amount of metal bound at sites II, III and IV to be monitored directly at partial degrees of saturation. The peak corresponding to site I could only be monitored at the beginning and end of the titration because of line broadening effects in the intermediate region of the titration. Sites III and IV both titrate preferentially and the results demonstrate clearly that sites in either domain fill effectively in parallel, consistent with a significant positive intradomain cooperativity of calcium binding.

Isolation of parvalbumin isotypes by preparative HPLC techniques

Parvalbumins are highly stable Ca2+ binding proteins, present in large quantities in the sarcoplasmic reticulum in the white muscle of most lower vertebrates and fish. The properties of these proteins make them promising antigens for the use as a specific biomarker for fish species identification. Parvalbumin isotypes were isolated, on a preparative scale level, by use of size exclusion chromatography (SEC) and anion exchange HPLC. The utility of this technique, with regard to maximizing purified isotypes, is discussed.

Conformational effects of calcium release from parvalbumin: comparison of computational simulations with spectroscopic investigations

The effect of Ca2+ binding to parvalbumin was monitored by probes of conformation including absorption, fluorescence, circular dichroism (CD), infrared (IR) spectroscopy and differential scanning calorimetry. These experimental studies were compared with molecular dynamics computations on the structures of the Ca-bound and Ca-free forms of cod parvalbumin. The UV CD spectra show that removal of calcium results in a decrease in the alpha-helical content of the protein. The IR amide I' and III' regions are very much affected by Ca removal and are indicative of significant perturbation of secondary structure. The fluorescence of tryptophan, the IR markers, and UV ellipticity all show changes with temperature, pointing to a lowering of protein stability upon Ca removal. These results are consistent with the structures obtained for both the Ca-bound and Ca-free proteins after 200 ps of solvated molecular dynamics simulations which show a decrease in the secondary structure upon Ca removal.

The isolation of parvalbumin isoforms from the tail muscle of the American alligator (Alligator mississipiensis)

Multiple parvalbumin isoforms have been detected in the tail (skeletal) muscle of the American alligator (Alligator mississipiensis). One of these isoforms (APV-1) has been highly purified and partially characterized. Protein purification involved mainly gel filtration and anion exchange chromatography, and characterization included gel electrophoresis, amino acid composition analysis, metal ion analysis, MALDI-TOF and ESI mass spectrometry, ultraviolet and fluorescence spectroscopy, and one- and two-dimensional 500 MHz proton NMR spectroscopy. The alligator isoforms are rich in phenylalanine and deficient in the other aromatic residues as is typical for parvalbumins. In fact, the one highly purified isoform that forms the basis of this study has only phenyl-alanine as an aromatic residue. Ion exchange chromatography further indicates that this isoform has a relatively high isoelectric point (pl approximately 5.0), indicating that it is an alpha-lineage parvalbumin. This alligator parvalbumin isoform is unusual in that it has an atypically high Ca2+ content (almost 3.0 mole of Ca2+ per mole of protein) following purification, a fact supported by terbium fluorescence titration experiments. Preliminary comparative analysis of the highly purified alligator parvalbumin isoform (in the Ca2-loaded state) by two-dimensional 1H-NMR (2D 1H TOCSY and 2D 1H NOESY) indicates that there is considerable similarity in structure between the alligator protein and a homologous protein obtained from the silver hake (a saltwater fish species).

Oncomodulin is abundant in the organ of Corti

A small, acidic Ca(2+)-binding protein (CBP-15) was recently detected in extracts of the mammalian auditory receptor organ, the organ of Corti [Senarita et al. (1995) Hear. Res. 90, 169-175]. N-terminal sequence data for CBP-15 [Thalmann et al. (1995) Biochem. Biophys. Res. Commun. 215, 142-147] implied membership in the parvalbumin family and possible identity with the mammalian beta-parvalbumin oncomodulin. As shown herein, the latter conclusion is supported by strong cross-reactivity between CBP-15 and isoform-specific antibodies to oncomodulin. Moreover, we have succeeded in amplifying the guinea pig CBP-15 coding sequence from organ of Corti cDNA using degenerate oligonucleotide primers based on the rat oncomodulin sequence. The deduced amino acid sequence of guinea pig CBP-15 displays 90%, 92%, and 98% identity with mouse, rat, and human oncomodulin isoforms. Demonstration of the presence of oncomodulin in the organ of Corti is the first documentation of this substance in a postnatal mammalian tissue.

Comparative modeling of the three-dimensional structure of the calmodulin-related TCH2 protein from Arabidopsis

Plants adapt to various stresses by developmental alterations that render them less easily damaged. Expression of the TCH2 gene of Arabidopsis is strongly induced by stimuli such as touch and wind. The gene product, TCH2, belongs to the calmodulin (CaM) family of proteins and contains four highly conserved Ca(2+)-binding EF-hands. We describe here the structure of TCH2 in the fully Ca(2+)-saturated form, constructed using comparative molecular modeling, based on the x-ray structure of paramecium CaM. Like known CaMs, the overall structure consists of two globular domains separated by a linker helix. However, the linker region has added flexibility due to the presence of 5 glycines within a span of 6 residues. In addition, TCH2 is enriched in Lys and Arg residues relative to other CaMs, suggesting a preference for targets which are more negatively charged. Finally, a pair of Cys residues in the C-terminal domain, Cys126 and Cys131, are sufficiently close in space to form a disulfide bridge. These predictions serve to direct future biochemical and structural studies with the overall aim of understanding the role of TCH2 in the cellular response of Arabidopsis to environmental stimuli.

The role of glycine (residue 89) in the central helix of EF-hand protein troponin-C exposed following amino-terminal alpha-helix deletion

Because an N-terminal alpha-helical (N-helix) arm and a KGK-triplet (residues 88KGK90) in the central helix of troponin-C (TnC) are missing in calmodulin, several recent studies have attempted to elucidate the structure-function correlations of these units. Presently, with a family of genetically manipulated derivatives especially developed for this study and tested on permeabilized isolated single skeletal muscle fiber segments, we explored the specificities of the amino acid residues within the N-helix and the KGK-triplet in TnC. Noticeably, the amino acid compositions vary between the N-helices of the cardiac and skeletal TnC isoforms. On the other hand, the KGK-triplet is located similarly in both TnC isoforms. We previously indicated that deletion of the N-helix (mutant delta Nt) diminishes the tension obtained on activation with maximal calcium, but the contractile function is revived by the superimposed deletion of the 88KGK90-triplet (mutant delta Nt delta KGK; see Gulati J, Babu A, Su H, Zhang YF, 1993, J Biol Chem 268:11685-11690). Using this functional test, we find that replacement of Gly-89 with a Leu or an Ala could also overcome the contractile defect associated with N-helix deletion. On the other hand, replacement of the skeletal TnC N-helix with cardiac type N-helix was unable to restore contractile function. The findings indicate a destabilizing influence of Gly-89 residue in skeletal TnC and suggest that the N-terminal arm in normal TnC serves to moderate this effect. Moreover, specificity of the N-helix between cardiac and skeletal TnCs raises the possibility that resultant structural disparities are also important for the functional distinctions of the TnC isoforms.

Molecular tuning of ion binding to calcium signaling proteins

Intracellular calcium plays an essential role in the transduction of most hormonal, neuronal, visual, and muscle stimuli. (Recent reviews include Putney, 1993; Berridge, 1993a,b; Tsunoda, 1993; Gnegy, 1993; Bachset al.1992; Hanson & Schulman, 1992; Villereal & Byron, 1992; Premack & Gardner, 1992; Meanset al.1991).

Oligomerization of an avian thymic parvalbumin. Chemical evidence for a Ca(2+)-specific conformation

CPV3, the third parvalbumin isoform to be identified in the chicken, is produced exclusively in the thymus gland. Although parvalbumins are typically cysteine-deficient, CPV3 contains two cysteine residues, at positions 18 and 72. The reported three-dimensional parvalbumin structures suggest that the side chain of cysteine-72 should be solvent-accessible. Accordingly, we find that CPV3 readily forms disulfide-linked oligomers in the absence of reducing agents. The reaction, employing either O2 or ferricyanide ion as the oxidant, is apparently restricted to the Ca(2+)-bound form of the protein. The differing reactivity of the Ca2+, Mg2+, and apo-forms has significant structural implications.

Human alpha and beta parvalbumins. Structure and tissue-specific expression

alpha and beta parvalbumins are Ca(2+)-binding proteins of the EF-hand type. We determined the protein sequence of human brain alpha parvalbumin by mass spectrometry and cloned human beta parvalbumin (or oncomodulin) from genomic DNA and preterm placental cDNA. beta parvalbumin differs in 54 positions from alpha parvalbumin and lacks the C-terminal amino acid 109. From MS analyses of alpha and beta parvalbumins we conclude that parvalbumins generally lack posttranslational modifications. alpha and beta parvalbumins were differently expressed in human tissues when analyzed by immunoblotting and polymerase-chain-reaction techniques. Whereas alpha parvalbumin was found in a number of adult human tissues, beta parvalbumin was restricted to preterm placenta. The pattern of alpha parvalbumin expression also differs in man compared to other vertebrates. For example, in rat, alpha parvalbumin was found in extrafusal and intrafusal skeletal-muscle fibres whereas, in man, alpha parvalbumin was restricted to the muscle spindles. Different functions for alpha and beta parvalbumins are discussed.

Kinetic control of Ca(II) signaling: tuning the ion dissociation rates of EF-hand Ca(II) binding sites

EF-hand Ca(II) binding sites share a conserved architecture and are prevalent in Ca(II) signaling pathways. The ion binding kinetics of these sites are carefully tuned to provide the physiologically appropriate activation and inactivation time scales. Here we examine kinetic tuning by the side chain at the ninth position of the EF-loop. A model is proposed in which both the size and charge of the side chain contribute to kinetic tuning. To test this model, the ninth loop position of the EF-hand-like site in the Escherichia coli D-galactose binding protein has been engineered and the Tb(III) dissociation kinetics of the resulting sites have been analyzed. Substitutions at this position are observed to generate up to 10(4)-fold changes in Tb(III) dissociation rates, with little effect on Tb(III) affinity. Furthermore, the observed pattern of rate changes confirm the model's predictions; long side chains at the ninth loop position yield slow dissociation kinetics as predicted for a steric block, whereas acidic side chains yield slow dissociation kinetics as expected for an electrostatic barrier.

Symmetrical rearrangement of the cation-binding sites of parvalbumin upon Ca2+/Mg2+ exchange. A study by 1H 2D NMR

Two forms of parvalbumin, i.e., the fully Ca-loaded form PaCa2 and the fully Mg-loaded form PaMg2, are investigated by 2D 1H NMR in solution. A detailed analysis of the resonances, which belong to residues involved in direct coordination of Ca2+ and Mg2+, establishes that the sixth ligand, a highly conserved Glu residue at the relative position 12 in both cation-binding sites CD and EF, undergoes a conformational rearrangement through a 120 degrees rotation of its side chain about the C alpha-C beta bond with PaMg2 adopting the less energetically favored g- conformation, as inferred from scalar coupling constants and dipole-dipole contacts measured on the COSY and NOESY spectra, respectively. Similarly, chemical shift effects, which selectively involve NH and C alpha H resonances (as well as side-chain resonances) in both CD and EF sites, point to a symmetrical behavior of both cation-binding sites upon Ca2+/Mg2+ exchange.

Evolution of EF-hand calcium-modulated proteins. II. Domains of several subfamilies have diverse evolutionary histories

In the first report in this series we described the relationships and evolution of 152 individual proteins of the EF-hand subfamilies. Here we add 66 additional proteins and define eight (CDC, TPNV, CLNB, LPS, DGK, 1F8, VIS, TCBP) new subfamilies and seven (CAL, SQUD, CDPK, EFH5, TPP, LAV, CRGP) new unique proteins, which we assume represent new subfamilies. The main focus of this study is the classification of individual EF-hand domains. Five subfamilies--calmodulin, troponin C, essential light chain, regulatory light chain, CDC31/caltractin--and three uniques--call, squidulin, and calcium-dependent protein kinase--are congruent in that all evolved from a common four-domain precursor. In contrast calpain and sarcoplasmic calcium-binding protein (SARC) each evolved from its own one-domain precursor. The remaining 19 subfamilies and uniques appear to have evolved by translocation and splicing of genes encoding the EF-hand domains that were precursors to the congruent eight and to calpain and to SARC. The rates of evolution of the EF-hand domains are slower following formation of the subfamilies and establishment of their functions. Subfamilies are not readily classified by patterns of calcium coordination, interdomain linker stability, and glycine and proline distribution. There are many homoplasies indicating that similar variants of the EF-hand evolved by independent pathways.