Molecular mechanisms of calcium and magnesium binding to parvalbumin

Insights into modulation of calcium signaling by magnesium in calmodulin, troponin C and related EF-hand proteins

The Ca(2+)-binding helix-loop-helix structural motif called "EF-hand" is a common building block of a large family of proteins that function as intracellular Ca(2+)-receptors. These proteins respond specifically to micromolar concentrations of Ca(2+) in the presence of ~1000-fold excess of the chemically similar divalent cation Mg(2+). The intracellular free Mg(2+) concentration is tightly controlled in a narrow range of 0.5-1.0mM, which at the resting Ca(2+) levels is sufficient to fully or partially saturate the Ca(2+)-binding sites of many EF-hand proteins. Thus, to convey Ca(2+) signals, EF-hand proteins must respond differently to Ca(2+) than to Mg(2+). In this review the structural aspects of Mg(2+) binding to EF-hand proteins are considered and interpreted in light of the recently proposed two-step Ca(2+)-binding mechanism (Grabarek, Z., J. Mol. Biol., 2005, 346, 1351). It is proposed that, due to stereochemical constraints imposed by the two-EF-hand domain structure, the smaller Mg(2+) ion cannot engage the ligands of an EF-hand in the same way as Ca(2+) and defaults to stabilizing the apo-like conformation of the EF-hand. It is proposed that Mg(2+) plays an active role in the Ca(2+)-dependent regulation of cellular processes by stabilizing the "off state" of some EF-hand proteins, thereby facilitating switching off their respective target enzymes at the resting Ca(2+) levels. Therefore, some pathological conditions attributed to Mg(2+) deficiency might be related to excessive activation of underlying Ca(2+)-regulated cellular processes. This article is part of a Special Issue entitled: 11th European Symposium on Calcium.

Metal cation complexation with natural organic matter in aqueous solutions: molecular dynamics simulations and potentials of mean force

Natural organic matter (NOM, or humic substance) has a known tendency to form colloidal aggregates in aqueous environments, with the composition and concentration of cationic species in solution, pH, temperature, and the composition of the NOM itself playing important roles. Strong interaction of carboxylic groups of NOM with dissolved metal cations is thought to be the leading chemical interaction in NOM supramolecular aggregation. Computational molecular dynamics (MD) study of the interactions of Na(+), Mg(2+), and Ca(2+) with the carboxylic groups of a model NOM fragment and acetate anions in aqueous solutions provides new quantitative insight into the structure, energetics, and dynamics of the interactions of carboxylic groups with metal cations, their association, and the effects of cations on the colloidal aggregation of NOM molecules. Potentials of mean force and the equilibrium constants describing overall ion association and the distribution of metal cations between contact ion pairs and solvent-separated ions pairs were computed from free MD simulations and restrained umbrella sampling calculations. The results provide insight into the local structural environments of metal-carboxylate association and the dynamics of exchange among these sites. All three cations prefer contact ion pair to solvent-separated ion pair coordination, and Na(+) and Ca(2+) show a strong preference for bidentate contact ion pair formation. The average residence time of a Ca(2+) ion in a contact ion pair with the carboxylic groups is of the order of 0.5 ns, whereas the corresponding residence time of a Na(+) ion is only between 0.02 and 0.05 ns. The average residence times of a Ca(2+) ion in a bidentate coordinated contact ion pair vs a monodentate coordinated contact ion pair are about 0.5 and 0.08 ns, respectively. On the 10 ns time scale of our simulations, aggregation of the NOM molecules occurs in the presence of Ca(2+) but not Na(+) or Mg(2+). These results agree with previous experimental observations and are explained by both Ca(2+) ion bridging between NOM molecules and decreased repulsion between the NOM molecules due to the reduced net charge of the NOM-metal complexes. Simulations on a larger scale are needed to further explore the relative importance of the different aggregation mechanisms and the stability of NOM aggregates.

Recombinant expression and affinity purification of snake venom gland parvalbumin in Escherichia coli

Parvalbumins (PV) are small, acidic, water soluble and calcium-binding proteins generally present in muscular and nervous tissues. In the present study, we identified and characterized a cDNA clone encoding PV, named AplPV, from a snake (Agkistrodon piscivorus leucostoma) venom gland cDNA library. AplPV belongs to EF-hand proteins with six alpha-helices constituting three EF-hand domains. The deduced amino acid sequence of AplPV is 91% and 68% identical to the previously characterized PVs of Boa constrictor and Cyprinus carpio, respectively. The full-length cDNA was subcloned into the expression vector pGEX and transformed into Escherichia coli (E.coli) to produce recombinant protein. The bacterially expressed GST-AplPV fusion protein was highly expressed, and effectively purified by Glutathione-Sepharose affinity chromatography. A high concentration of thrombin protease specifically cleaved and removed the GST tag from fusion protein, and further purified by Benzamidine column for removal of thrombin protease. As a result, the 12 kDa AplPV recombinant protein alone was purified. To investigate the tissue-specific biological occurrence of AplPV, a polyclonal antibody (anti-AplPV-antibody) was raised against GST-AplPV fusion protein in rabbit. Western blot analysis revealed that immunoreactive bands were exhibited in both recombinant protein and samples of venom glands, but not in any crude venom. This specific occurrence indicates a specialized function of AplPV in snake venom glands.

Comparison of natural and recombinant forms of the major fish allergen parvalbumin from cod and carp

Allergic reaction following fish consumption can trigger life-threatening reactions in predisposed individuals. Parvalbumins from different species have been identified as the major fish allergens. There are two distinct phylogenetic lineages of parvalbumins, alpha and beta. Most allergic reactions are caused by beta-parvalbumins. We cloned and expressed cDNAs encoding cod (Gadus morhua) and carp (Cyprinus carpio) beta-parvalbumins and purified natural cod beta-parvalbumin. CD spectra of the purified proteins showed that their overall secondary structure contents were very similar. No differences in thermal stability were monitored in the calcium-bound or calcium-depleted form of natural cod parvalbumin. IgE reactivity was assessed using 26 sera of fish allergic patients from Spain, The Netherlands, and Greece in immunoblot and ELISA experiments. Twenty-five of the 26 patients with IgE reactivity to native and recombinant cod parvalbumin also reacted to the recombinant carp parvalbumin. IgE inhibition assays were performed using cod and carp extracts and purified recombinant parvalbumin of cod and carp. High crossreactivity among cod and carp parvalbumins was observed in immunoblots as well as in fluid phase assays. Natural and recombinant parvalbumins gave comparable results when performing various in vitro diagnostic assays.

Structure-function analysis of plant aquaporin AtPIP2;1 gating by divalent cations and protons

Water channel proteins, AQPs (aquaporins), of the PIP (plasma membrane intrinsic protein) subfamily, provide a means for fine and quick adjustments of the plant water status. A molecular model for gating of PIPs by cytosolic protons (H(+)) and divalent cations was derived from the atomic structure of spinach SoPIP2;1 (Spinacia oleracea PIP2;1) in an open- and a closed-pore conformation. In the present study, we produced the Arabidopsis AtPIP2;1 (Arabidopsis thaliana PIP2;1) homologue in Pichia pastoris, either WT (wild-type) or mutations at residues supposedly involved in gating. Stopped-flow spectrophotometric measurements showed that, upon reconstitution in proteoliposomes, all forms function as water channels. The first functional evidence for a direct gating of PIPs by divalent (bivalent) cations was obtained. In particular, cadmium and manganese were identified, in addition to calcium (Ca(2+)) and H(+) as potent inhibitors of WT AtPIP2;1. Our results further show that His(199), the previously identified site for H(+) sensing, but also N-terminal located Glu(31), and to a lesser extent Asp(28), are involved in both divalent-cation- and H(+)-mediated gating. In contrast, mutation of Arg(124) rendered AtPIP2;1 largely insensitive to Ca(2+) while remaining fully sensitive to H(+). The role of these residues in binding divalent cations and/or stabilizing the open or closed pore conformations is discussed.

Molecular dynamics study of intrinsic stability in six RNA terminal loop motifs

Single stranded RNA molecules can assume a wide range of tertiary structures beyond the canonical A-form double helix. Certain sequences, termed motifs, are more common than a random distribution would suggest. The existence of such motifs can be rationalized in structural terms. In this study, we have investigated the intrinsic structural stability of RNA terminal loop motifs using multiple MD simulations in explicit water. Representative loops were chosen from the major tetraloop motifs, including also the U-turn motif. Not all loops retain their folded starting structure, but lowering the temperature to 277 K, or adding adjacent base pairs from the stem to which the motif is attached, helps stabilizing the folded loop structure.

Activation and inhibition of photoreceptor guanylyl cyclase by guanylyl cyclase activating protein 1 (GCAP-1): the functional role of Mg2+/Ca2+ exchange in EF-hand domains

Guanylyl cyclase activating protein 1 (GCAP-1), a Ca(2+)/Mg(2+) sensor protein that accelerates retinal guanylyl cyclase (RetGC) in the light and decelerates it in the dark, is inactive in cation-free form. Binding of Mg(2+) in EF-hands 2 and 3 was essential for RetGC activation in the conditions mimicking light adaptation. Mg(2+) binding in EF-hand 2 affected the conformation of a neighboring non-metal binding domain, EF-hand-1, and increased GCAP-1 affinity for RetGC nearly 40-fold compared with the metal-free EF-hand 2. Mg(2+) binding in EF-hand 3 increased GCAP-1 affinity for RetGC 5-fold and its maximal RetGC stimulation 2-fold. Mg(2+) binding in EF-hand 4 affected neither GCAP-1 affinity for RetGC, nor RetGC activation. Inactivation of Ca(2+) binding in EF-hand 4 was sufficient to render GCAP-1 a constitutive activator of RetGC, whereas the EF-hand 3 role in Ca(2+)-dependent deceleration of RetGC was likely to be through the neighboring EF-hand 4. Inactivation of Ca(2+) binding in EF-hand 2 affected cooperativity of RetGC inhibition by Ca(2+), but did not prevent the inhibition. We conclude that 1) Mg(2+) binding in EF-hands 2 and 3, but not EF-hand 4, is essential for the ability of GCAP-1 to activate RetGC in the light; 2) Mg(2+) or Ca(2+) binding in EF-hand 3 and especially in EF-hand 2 is required for high-affinity interaction with the cyclase and affects the conformation of the neighboring EF-hand 1, a domain required for targeting RetGC; and 3) RetGC inhibition is likely to be primarily caused by Ca(2+) binding in EF-hand 4.

Magnesium Promotes Structural Integrity and Conformational Switching Action of a Calcium Sensor Protein

Calcium binding proteins carry out various signal transduction processes upon binding to Ca2+. In general, these proteins perform their functions in a high background of Mg2+. Here, we report the role of Mg2+ on a calcium sensor protein from Entamoeba histolytica (EhCaBP), containing four Ca2+-binding sites. Mg2+-bound EhCaBP exists as a monomer with a conformation different from that of the holo- and apo-EhCaBP. NMR and biophysical data on EhCaBP demonstrate that Mg2+ stabilizes the closed conformation of the apo form. In the presence of Mg2+, the partially collapsed apo-EhCaBP gains stability and structural integrity. Mg2+ binds to only 3 out of 4 calcium binding sites in EhCaBP. The Ca2+ binding affinity and cooperativity of the conformational switching from the "closed" to the "open" state is significantly modulated by the presence of Mg2+. This fine-tuning of the Ca2+ concentration to switch its conformation is essential for CaBPs to carry out the signal transduction process efficiently.

Genome-wide identification and analyses of the rice calmodulin and related potential calcium sensor proteins

BACKGROUND

A wide range of stimuli evoke rapid and transient increases in [Ca2+]cyt in plant cells which are transmitted by protein sensors that contain EF-hand motifs. Here, a group of Oryza sativa L. genes encoding calmodulin (CaM) and CaM-like (CML) proteins that do not possess functional domains other than the Ca2+-binding EF-hand motifs was analyzed.

RESULTS

By functional analyses and BLAST searches of the TIGR rice database, a maximum number of 243 proteins that possibly have EF-hand motifs were identified in the rice genome. Using a neighbor-joining tree based on amino acid sequence similarity, five loci were defined as Cam genes and thirty two additional CML genes were identified. Extensive analyses of the gene structures, the chromosome locations, the EF-hand motif organization, expression characteristics including analysis by RT-PCR and a comparative analysis of Cam and CML genes in rice and Arabidopsis are presented.

CONCLUSION

Although many proteins have unknown functions, the complexity of this gene family indicates the importance of Ca2+-signals in regulating cellular responses to stimuli and this family of proteins likely plays a critical role as their transducers.

Effect of carboxylate-binding mode on metal binding/selectivity and function in proteins

We delineate the factors governing the carboxylate-binding mode (monodentate vs bidentate) in metalloproteins. We reveal how the carboxylate-binding mode affects the binding affinity and selectivity of a metal ion as well as the function of a metalloprotein using Ca2+-binding proteins and enzymes (ribonuclease H1, phosphoserine phosphatase, and ribonucleotide reductase) as examples. The collected data indicate that a carboxylate monodentate left arrow over right arrow bidentate switch, in addition to other structural factors, could be used to fine tune the metal-binding site affinity and/or selectivity, thus modifying the function/properties of the metalloprotein.

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.

Ca2+ and Mg2+ Binding Properties of GCAP-1

Guanylyl cyclase-activating protein 1 (GCAP-1) is an EF-hand protein that activates retinal guanylyl cyclase (RetGC) in photoreceptors at low free Ca2+ in the light and inhibits it in the dark when Ca2+ concentrations rise. We present the first direct evidence that Mg2+-bound form of GCAP-1, not its cation-free form, is the true activator of RetGC-1 under physiological conditions. Of four EF-hand structures in GCAP-1, three bound Ca2+ ions and could exchange Ca2+ for Mg2+. At concentrations of free Ca2+ and Mg2+ typical for the light-adapted photoreceptors, all three metal-binding EF-hands were predominantly occupied by Mg2, and the presence of bound Mg2+ in GCAP-1 was essential for its ability to stimulate RetGC-1. In the Mg2+-bound form of GCAP-1 all three Trp residues became more exposed to the polar environment compared with its apo form. The replacement of Mg2+ by Ca2+ in the EF-hands 2 and 3 further exposed Trp-21 to the solution in a non-metal-binding EF-hand domain 1 that interacts with RetGC. Contrary to that, replacement of Mg2+ by Ca2+ in the EF-hand 4 moved Trp-94 in the entering alpha-helix of the EF-hand 3 back to the non-polar environment. Our results demonstrate that Mg2+ regulates GCAP-1 not only by adjusting its Ca2+ sensitivity to the physiological conditions in photoreceptors but also by creating the conformation required for RetGC stimulation.

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.

Handling calcium signaling: Arabidopsis CaMs and CMLs

The Arabidopsis genome harbors seven calmodulin (CAM) and 50 CAM-like (CML) genes that encode potential calcium sensors. The CAMs encode only four protein isoforms. Selective pressure to maintain multiple CAMs indicates nonredundancy. Sequence divergence, even in the EF hand calcium-binding motif, exists among the CMLs and, therefore, divergent functions are likely to have evolved. Expression data recently available from Massively Parallel Signature Sequencing and Genevestigator compilation of microarrays are reviewed. The seven Arabidopsis CAMs are highly and relatively uniformly expressed. Differential expression is evident among the distinct CMLs over developmental stages, in various organs and in response to many different stimuli. In spite of the potential importance in mediating plant calcium signaling, the physiological functions of the Arabidopsis CaMs and CMLs remain largely unknown.

Mg2+ and Ca2+ differentially regulate DNA binding and dimerization of DREAM

DREAM (calsenilin/KChIP3) is an EF-hand calcium-binding protein that represses transcription of prodynorphin and c-fos genes. Here we present structural and binding studies on single-site mutants of DREAM designed to disable Ca(2+) binding to each of the functional EF-hands (EF-2: D150N; EF-3: E186Q; and EF-4: E234Q). Isothermal titration calorimetry (ITC) analysis of Ca(2+) binding to the various mutants revealed that, in the absence of Mg(2+), Ca(2+) binds independently and sequentially to EF-3 (DeltaH = -2.4 kcal/mol), EF-4 (DeltaH = +5.2 kcal/mol), and EF-2 (DeltaH = +1 kcal/mol). By contrast, only two Ca(2+) bind to DREAM in the presence of physiological levels of Mg(2+) for both wild-type and D150N, suggesting that EF-2 binds constitutively to Mg(2+). ITC measurements demonstrate that one Mg(2+) binds enthalpically with high affinity (K(d) = 13 mum and DeltaH = -0.79 kcal/mol) and two or more Mg(2+) bind entropically in the millimolar range. Size-exclusion chromatography studies revealed that Mg(2+) stabilizes DREAM as a monomer, whereas Ca(2+) induces protein dimerization. Electrophoretic mobility shift assays indicated that Mg(2+) is essential for sequence-specific binding of DREAM to DNA response elements (DREs) in prodynorphin and c-fos genes. The EF-hand mutants bind specifically to DRE, suggesting they are functionally intact. None of the EF-hand mutants bind DRE at saturating Ca(2+) levels, suggesting that binding of a single Ca(2+) at either EF-3 or EF-4 is sufficient to drive conformational changes that abolish DNA binding. NMR structural analysis indicates that metal-free DREAM adopts a folded yet flexible molten globule-like structure. Both Ca(2+) and Mg(2+) induce distinct conformational changes, which stabilize tertiary structure of DREAM. We propose that Mg(2+) binding at EF-2 may structurally bridge DREAM to DNA targets and that Ca(2+)-induced protein dimerization disrupts DNA binding.

Temperature sensitivity of calcium binding for parvalbumins from Antarctic and temperate zone teleost fishes

Parvalbumin (PV) is a soluble calcium-binding protein that is especially abundant in fast-twitch muscles of fish and other lower vertebrates. Despite its prevalence in ectothermic taxa, few data address the effects of temperature on PV binding function. In this study, calcium dissociation constants (KD) were measured as a function of temperature (0-25 degrees C) for PV from two Antarctic (Gobionotothen gibberifrons and Chaenocephalus aceratus) and two temperate zone fish species (Cyprinus carpio and Micropterus salmoides). Measurements by fluorometric competitive binding assay show that KD values for PVs from the Antarctic species were significantly higher at all assay temperatures and were less sensitive to temperature relative to carp and bass. However, estimates of KD are fundamentally similar for PVs from the Antarctic and temperate zone species when examined at their native physiological temperature. Variation in pH and ionic strength within a physiologically relevant range had only modest effects on KD. Thermodynamics of calcium binding to PV from G. gibberifrons and C. carpio was measured by isothermal microcalorimetry. When measured at 15 degrees C, the Gibbs free energy change (deltaG) was significantly greater for calcium binding to PV from G. gibberifrons than from carp (-43.4+/-1.5 kJ mol(-1) and -46.6+/-3.0 kJ mol(-1), respectively), and the relative contribution of entropy to deltaG for calcium binding to PV from the Antarctic species was about twice that of carp (deltaS=16.0+/-0.8 J degrees C(-1) mol(-1) for G. gibberifrons; deltaS=7.5+/-0.8 J degrees C(-1) mol(-1) for C. carpio).

The endogenous calcium ions of horseradish peroxidase C are required to maintain the functional nonplanarity of the heme

Horseradish peroxidase C (HRPC) binds 2 mol calcium per mol of enzyme with binding sites located distal and proximal to the heme group. The effect of calcium depletion on the conformation of the heme was investigated by combining polarized resonance Raman dispersion spectroscopy with normal coordinate structural decomposition analysis of the hemes extracted from models of Ca(2+)-bound and Ca(2+)-depleted HRPC generated and equilibrated using molecular dynamics simulations. Results show that calcium removal causes reorientation of heme pocket residues. We propose that these rearrangements significantly affect both the in-plane and out-of-plane deformations of the heme. Analysis of the experimental depolarization ratios are clearly consistent with increased B(1g)- and B(2g)-type distortions in the Ca(2+)-depleted species while the normal coordinate structural decomposition results are indicative of increased planarity for the heme of Ca(2+)-depleted HRPC and of significant changes in the relative contributions of three of the six lowest frequency deformations. Most noteworthy is the decrease of the strong saddling deformation that is typical of all peroxidases, and an increase in ruffling. Our results confirm previous work proposing that calcium is required to maintain the structural integrity of the heme in that we show that the preferred geometry for catalysis is lost upon calcium depletion.