Stritch School of Medicine, Maywood, Illinois 60153, United States.
J Phys Chem B. 2021 Jun 24;125(24):6390-6405. doi: 10.1021/acs.jpcb.1c01269. Epub 2021 Jun 11.
Members of the parvalbumin (PV) family of calcium (Ca) binding proteins (CBPs) share a relatively high level of sequence similarity. However, their Ca affinities and selectivities against competing ions like Mg can widely vary. We conducted molecular dynamics simulations of several α-parvalbumin (αPV) constructs with micromolar to nanomolar Ca affinities to identify structural and dynamic features that contribute to their binding of ions. Specifically, we examined a D94S/G98E construct with a lower Ca 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 Mg to these isoforms, which is much weaker than Ca. 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 Mg from Ca, they did not indicate that Ca 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 Ca binding by -20 kcal/mol relative to WT versus 30 kcal/mol for D94S/G98E αPV. Meanwhile, Mg 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, Ca affinity and selectivity against Mg 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.
钙结合蛋白(CBPs)家族的成员 Parvalbumin(PV)具有相对较高的序列相似性。然而,它们对竞争离子(如 Mg)的亲和力和选择性差异很大。我们对几种具有微摩尔至纳摩尔 Ca 亲和力的α-PV(αPV)结构进行了分子动力学模拟,以确定导致其离子结合的结构和动态特征。具体来说,我们研究了一个 D94S/G98E 构建体,其 Ca 亲和力(≈-18 kcal/mol)相对于野生型(WT)(≈-22 kcal/mol)较低,以及一个 S55D/E59D 变体,其亲和力增强(≈-24 kcal/mol)。此外,我们还研究了这些同工型对 Mg 的结合,其结合强度远低于 Ca。我们使用平均球近似(MSA)理论来评估离子在蛋白质 EF 手结构域中的结合热力学,以解释离子有限大小和周围电解质组成的影响。虽然 MSA 分数区分了 Mg 和 Ca,但它们并没有表明 PV 同工型之间在结合环处的 Ca 结合亲和力不同。相反,分子力学广义 Born 表面积(MM/GBSA)近似能量,我们用它来量化蛋白质在结合离子时结构重排的热力学成本,表明 S55D/E59D αPV 相对于 WT 对 Ca 的结合亲和力有利,为-20 kcal/mol,而 D94S/G98E αPV 则为 30 kcal/mol。同时,S55D/E59D αPV 和 D94S/G98E αPV 变体对 Mg 的结合更有利,分别为-18.32 和-1.65 kcal/mol。这些能量表明,除了结合环处的氧配位外,离子结合还有其他重要贡献,这些贡献源于 α-螺旋、β-折叠特征和氢键的变化。因此,Ca 亲和力和对 Mg 的选择性是由蛋白质离子结合位点的局部效应以及其他非局部贡献共同产生的突现性质。我们的研究结果拓宽了我们对α-PV 离子结合的分子基础的理解,这可能是由广泛的 CBPs 家族成员共享的。