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

Structural Changes beyond the EF-Hand Contribute to Apparent Calcium Binding Affinities: Insights from Parvalbumins

Members of the parvalbumin (PV) family of calcium (Ca2+) binding proteins (CBPs) share a relatively high level of sequence similarity. However, their Ca2+ affinities and selectivities against competing ions like Mg2+ can widely vary. We conducted molecular dynamics simulations of several α-parvalbumin (αPV) constructs with micromolar to nanomolar Ca2+ affinities to identify structural and dynamic features that contribute to their binding of ions. Specifically, we examined a D94S/G98E construct with a lower Ca2+ affinity (≈-18 kcal/mol) relative to the wild type (WT) (≈-22 kcal/mol) and an S55D/E59D variant with enhanced affinity (≈-24 kcal/mol). Additionally, we also examined the binding of Mg2+ to these isoforms, which is much weaker than Ca2+. We used mean spherical approximation (MSA) theory to evaluate ion binding thermodynamics within the proteins' EF-hand domains to account for the impact of ions' finite sizes and the surrounding electrolyte composition. While the MSA scores differentiated Mg2+ from Ca2+, they did not indicate that Ca2+ binding affinities at the binding loop differed between the PV isoforms. Instead, molecular mechanics generalized Born surface area (MM/GBSA) approximation energies, which we used to quantify the thermodynamic cost of structural rearrangement of the proteins upon binding ions, indicated that S55D/E59D αPV favored Ca2+ binding by -20 kcal/mol relative to WT versus 30 kcal/mol for D94S/G98E αPV. Meanwhile, Mg2+ binding was favored for the S55D/E59D αPV and D94S/G98E αPV variants by -18.32 and -1.65 kcal/mol, respectively. These energies implicate significant contributions to ion binding beyond oxygen coordination at the binding loop, which stemmed from changes in α-helicity, β-sheet character, and hydrogen bonding. Hence, Ca2+ affinity and selectivity against Mg2+ are emergent properties stemming from both local effects within the proteins' ion binding sites as well as non-local contributions elsewhere. Our findings broaden our understanding of the molecular bases governing αPV ion binding that are likely shared by members of the broad family of CBPs.

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

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

Characterization of avian thymic hormone and chicken parvalbumin 3 target cells

Avian thymic hormone (ATH) is a β-parvalbumin produced by epithelial cells in the thymic cortex and in the eyes of chickens. Chicken parvalbumin 3 (CPV3) is a homologous protein produced in the thymus and in hair cells of the chicken ear. ATH circulates in the blood on a five-day cycle and stimulates cell-mediated immunity when administered to young chickens. We report the identification of target cells for ATH and CPV3 and the immunophenotype of target cells for ATH. Newly hatched chicks were injected intracoelomically with ATH and killed 5, 10, 15 or 20 min later. Naïve chickens also were killed at 1, 7 and 14 days of age. Various tissues were examined by EM for the presence of either ATH or CPV3 using colloidal gold labeling. Gold particles were initially present on plasma membranes of lymphocytes in T cell areas of spleen and cecal tonsils from the chicks injected with ATH, internalized within 10 min, and accumulated in nuclei by 20 min. Immunofluorescence staining also identified the presence of ATH in T cell areas of spleen and cecal tonsils. Target cells labeled for ATH were immunophenotyped by double labeling. They were positive for CD3, CD8 and the lymphocyte receptor TCR-1, a phenotype characteristic of cytotoxic γδ T cells. Some of the target cells in the spleen were TCR-3 positive. Targeting of lymphocytes by CPV3 indicated that it may also be an immunomodulating hormone.

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.

Leucine 85 is an important determinant of divalent ion affinity in rat beta-parvalbumin (Oncomodulin)

Despite 69% sequence identity with chicken parvalbumin 3 (CPV3), rat beta-parvalbumin (beta-PV) exhibits a substantially lower Ca(2+) affinity (DeltaDeltaG degrees ' = 2.0 kcal/mol). This difference largely reflects the disparate behavior of the respective CD sites. Replacement of the rat beta-PV codon with the CPV3 codon at positions 49, 50, and 57-60 produces virtual sequence identity with the CPV3 CD site. However, the resulting protein exhibits a modest (0.5 kcal/mol) improvement in Ca(2+) affinity, implying that sequence differences beyond the binding site modulate divalent ion binding behavior. The solution structure of Ca(2+)-free rat beta-PV suggested that Leu-85, phenylalanine in CPV3, might be an important determinant. Therefore, the impact of the L85F mutation on divalent ion affinity was examined in rat beta-PV, in the variant harboring all six of the aforementioned CD site mutations, and in the intermediate CD site variants. We find that the identity of residue 85, located within the E helix, strongly influences divalent ion affinity in the mammalian beta-PV isoform and that its impact is mediated by interactions with residues in the CD site. In the wild-type protein, L85F primarily affects the EF site. By contrast, in the presence of the six CD site mutations, L85F also improves the CD site performance, yielding a protein with Ca(2+) affinity rivaling that of CPV3 and markedly enhanced Mg(2+) affinity as well. The impact of L85F on CD site Ca(2+) affinity is particularly sensitive to the identities of residues 59 and 60. Interestingly, however, significant improvement in CD site Mg(2+) affinity also requires mutation of additional CD site residues.