Calbindin D(9k): a protein optimized for calcium binding at neutral pH

Physiology of intracellular calcium buffering

Calcium signaling underlies much of physiology. Almost all the Ca2+ in the cytoplasm is bound to buffers, with typically only ∼1% being freely ionized at resting levels in most cells. Physiological Ca2+ buffers include small molecules and proteins, and experimentally Ca2+ indicators will also buffer calcium. The chemistry of interactions between Ca2+ and buffers determines the extent and speed of Ca2+ binding. The physiological effects of Ca2+ buffers are determined by the kinetics with which they bind Ca2+ and their mobility within the cell. The degree of buffering depends on factors such as the affinity for Ca2+, the Ca2+ concentration, and whether Ca2+ ions bind cooperatively. Buffering affects both the amplitude and time course of cytoplasmic Ca2+ signals as well as changes of Ca2+ concentration in organelles. It can also facilitate Ca2+ diffusion inside the cell. Ca2+ buffering affects synaptic transmission, muscle contraction, Ca2+ transport across epithelia, and the killing of bacteria. Saturation of buffers leads to synaptic facilitation and tetanic contraction in skeletal muscle and may play a role in inotropy in the heart. This review focuses on the link between buffer chemistry and function and how Ca2+ buffering affects normal physiology and the consequences of changes in disease. As well as summarizing what is known, we point out the many areas where further work is required.

Computational Investigation of Mechanisms for pH Modulation of Human Chloride Channels

Many transmembrane proteins are modulated by intracellular or extracellular pH. Investigation of pH dependence generally proceeds by mutagenesis of a wide set of amino acids, guided by properties such as amino-acid conservation and structure. Prediction of pKas can streamline this process, allowing rapid and effective identification of amino acids of interest with respect to pH dependence. Commencing with the calcium-activated chloride channel bestrophin 1, the carboxylate ligand structure around calcium sites relaxes in the absence of calcium, consistent with a measured lack of pH dependence. By contrast, less relaxation in the absence of calcium in TMEM16A, and maintenance of elevated carboxylate sidechain pKas, is suggested to give rise to pH-dependent chloride channel activity. This hypothesis, modulation of calcium/proton coupling and pH-dependent activity through the extent of structural relaxation, is shown to apply to the well-characterised cytosolic proteins calmodulin (pH-independent) and calbindin D9k (pH-dependent). Further application of destabilised, ionisable charge sites, or electrostatic frustration, is made to other human chloride channels (that are not calcium-activated), ClC-2, GABAA, and GlyR. Experimentally determined sites of pH modulation are readily identified. Structure-based tools for pKa prediction are freely available, allowing users to focus on mutagenesis studies, construct hypothetical proton pathways, and derive hypotheses such as the model for control of pH-dependent calcium activation through structural flexibility. Predicting altered pH dependence for mutations in ion channel disorders can support experimentation and, ultimately, clinical intervention.

Benefits and constrains of covalency: the role of loop length in protein stability and ligand binding

Protein folding is governed by non-covalent interactions under the benefits and constraints of the covalent linkage of the backbone chain. In the current work we investigate the influence of loop length variation on the free energies of folding and ligand binding in a small globular single-domain protein containing two EF-hand subdomains—calbindin D_9k. We introduce a linker extension between the subdomains and vary its length between 1 to 16 glycine residues. We find a close to linear relationship between the linker length and the free energy of folding of the Ca²⁺-free protein. In contrast, the linker length has only a marginal effect on the Ca²⁺ affinity and cooperativity. The variant with a single-glycine extension displays slightly increased Ca²⁺ affinity, suggesting that the slightly extended linker allows optimized packing of the Ca²⁺-bound state. For the extreme case of disconnected subdomains, Ca²⁺ binding becomes coupled to folding and assembly. Still, a high affinity between the EF-hands causes the non-covalent pair to retain a relatively high apparent Ca²⁺ affinity. Our results imply that loop length variation could be an evolutionary option for modulating properties such as protein stability and turnover without compromising the energetics of the specific function of the protein.

Accurate Electron-Nucleus Distances from Paramagnetic Relaxation Enhancements

Measurements of paramagnetic relaxation enhancements (PREs) in ¹H NMR spectra are an important tool to obtain long-range distance information in proteins, but quantitative interpretation is easily compromised by nonspecific intermolecular PREs. Here we show that PREs generated by lanthanides with anisotropic magnetic susceptibilities offer a route to accurate calibration-free distance measurements. As these lanthanides change ¹H chemical shifts due to pseudocontact shifts, the relaxation rates in the paramagnetic and diamagnetic state can be measured with a single sample that simultaneously contains the protein labeled with a paramagnetic and a diamagnetic lanthanide ion. Nonspecific intermolecular PREs are thus automatically subtracted when calculating the PREs as the difference in nuclear relaxation rates between paramagnetic and diamagnetic protein. Although PREs from lanthanides with anisotropic magnetic susceptibilities are complicated by additional cross-correlation effects and residual dipolar couplings (RDCs) in the paramagnetic state, these effects can be controlled by the choice of lanthanide ion and experimental conditions. Using calbindin D_9k with erbium, we succeeded in measuring intramolecular PREs with unprecedented accuracy, resulting in distance predictions with a root-mean-square-deviation of <0.9 Å in the range 11-24 Å.

Na+/H+ exchanger 3 inhibitor diminishes hepcidin-enhanced duodenal calcium transport in hemizygous β-globin knockout thalassemic mice

Recent investigation has shown that the liver-derived iron-regulating hormone, hepcidin, can potentiate intestinal calcium absorption in hemizygous β-globin knockout thalassemic (BKO) mice. Since the upregulation of Fe²⁺ and H⁺ cotransporter, divalent metal transporter (DMT)-1, has been shown to correlate with thalassemia-induced intestinal calcium absorption impairment, the inhibition of the apical Na⁺/H⁺ exchanger (NHE)-3 that is essential for cytoplasmic pH regulation and transepithelial sodium absorption was hypothesized to negatively affect hepcidin action. Herein, the positive effect of hepcidin on the duodenal calcium transport was evaluated using Ussing chamber technique. The results showed that BKO mice had lower absorptive surface area and duodenal calcium transport than wild-type mice. Besides, paracellular transport of zinc in BKO mice was compromised. Hepcidin administration completely restored calcium transport. Since this hepcidin action was totally abolished by inhibitors of the basolateral calcium transporters, Na⁺/Ca²⁺ exchanger (NCX1) and plasma membrane Ca²⁺-ATPase (PMCA_1b), the enhanced calcium flux potentially occurred through the transcellular pathway rather than paracellular pathway. Interestingly, the selective NHE3 inhibitor, 100 nM tenapanor, markedly inhibited hepcidin-enhanced calcium transport. Accordingly, hepcidin is one of the promising therapeutic agents for calcium malabsorption in β-thalassemia. It mainly stimulates the transcellular calcium transport across the duodenal epithelium in an NHE3-dependent manner.

Na+/H+ exchanger 3 inhibitor diminishes the amino-acid-enhanced transepithelial calcium transport across the rat duodenum

Na⁺/H⁺ exchanger (NHE)-3 is important for intestinal absorption of nutrients and minerals, including calcium. The previous investigations have shown that the intestinal calcium absorption is also dependent on luminal nutrients, but whether aliphatic amino acids and glucose, which are abundant in the luminal fluid during a meal, similarly enhance calcium transport remains elusive. Herein, we used the in vitro Ussing chamber technique to determine epithelial electrical parameters, i.e., potential difference (PD), short-circuit current (Isc), and transepithelial resistance, as well as ⁴⁵Ca flux in the rat duodenum directly exposed on the mucosal side to glucose or various amino acids. We found that mucosal glucose exposure led to the enhanced calcium transport, PD, and Isc, all of which were insensitive to NHE3 inhibitor (100 nM tenapanor). In the absence of mucosal glucose, several amino acids (12 mM in the mucosal side), i.e., alanine, isoleucine, leucine, proline, and hydroxyproline, markedly increased the duodenal calcium transport. An inhibitor for NHE3 exposure on the mucosal side completely abolished proline- and leucine-enhanced calcium transport, but not transepithelial transport of both amino acids themselves. In conclusion, glucose and certain amino acids in the mucosal side were potent stimulators of the duodenal calcium absorption, but only amino-acid-enhanced calcium transport was NHE3-dependent.

Optogenetic approaches addressing extracellular modulation of neural excitability

The extracellular ionic environment in neural tissue has the capacity to influence, and be influenced by, natural bouts of neural activity. We employed optogenetic approaches to control and investigate these interactions within and between cells, and across spatial scales. We began by developing a temporally precise means to study microdomain-scale interactions between extracellular protons and acid-sensing ion channels (ASICs). By coupling single-component proton-transporting optogenetic tools to ASICs to create two-component optogenetic constructs (TCOs), we found that acidification of the local extracellular membrane surface by a light-activated proton pump recruited a slow inward ASIC current, which required molecular proximity of the two components on the membrane. To elicit more global effects of activity modulation on ‘bystander’ neurons not under direct control, we used densely-expressed depolarizing (ChR2) or hyperpolarizing (eArch3.0, eNpHR3.0) tools to create a slow non-synaptic membrane current in bystander neurons, which matched the current direction seen in the directly modulated neurons. Extracellular protons played contributory role but were insufficient to explain the entire bystander effect, suggesting the recruitment of other mechanisms. Together, these findings present a new approach to the engineering of multicomponent optogenetic tools to manipulate ionic microdomains, and probe the complex neuronal-extracellular space interactions that regulate neural excitability.

A microfluidic platform for quantitative measurements of effective protein charges and single ion binding in solution

Microfluidic electrophoresis enables the comparison of dry sequence and solvated protein charges, and the detection of protein–ion binding.

Endoplasmic reticulum stress or mutation of an EF-hand Ca(2+)-binding domain directs the FKBP65 rotamase to an ERAD-based proteolysis

FKBP65 is an endoplasmic reticulum (ER)-localized chaperone and rotamase, with cargo proteins that include tropoelastin and collagen. In humans, mutations in FKBP65 have recently been shown to cause a form of osteogenesis imperfecta (OI), a brittle bone disease resulting from deficient secretion of mature type I collagen. In this work, we describe the rapid proteolysis of FKBP65 in response to ER stress signals that activate the release of ER Ca(2+) stores. A large-scale screen for stress-induced cellular changes revealed FKBP65 proteins to decrease within 6-12 h of stress activation. Inhibiting IP(3)R-mediated ER Ca(2+) release blocked this response. No other ER-localized chaperone and folding mediators assessed in the study displayed this phenomenon, indicating that this rapid proteolysis of folding mediator is distinctive. Imaging and cellular fractionation confirmed the localization of FKBP65 (72 kDa glycoprotein) to the ER of untreated cells, a rapid decrease in protein levels following ER stress, and the corresponding appearance of a 30-kDa fragment in the cytosol. Inhibition of the proteasome during ER stress revealed an accumulation of FKBP65 in the cytosol, consistent with retrotranslocation and a proteasome-based proteolysis. To assess the role of Ca(2+)-binding EF-hand domains in FKBP65 stability, a recombinant FKBP65-GFP construct was engineered to ablate Ca(2+) binding at each of two EF-hand domains. Cells transfected with the wild-type construct displayed ER localization of the FKBP65-GFP protein and a proteasome-dependent proteolysis in response to ER stress. Recombinant FKBP65-GFP carrying a defect in the EF1 Ca(2+)-binding domain displayed diminished protein in the ER when compared to wild-type FKBP65-GFP. Proteasome inhibition restored mutant protein to levels similar to that of the wild-type FKBP65-GFP. A similar mutation in EF2 did not confer FKBP65 proteolysis. This work supports a model in which stress-induced changes in ER Ca(2+) stores induce the rapid proteolysis of FKBP65, a chaperone and folding mediator of collagen and tropoelastin. The destruction of this protein may identify a cellular strategy for replacement of protein folding machinery following ER stress. The implications for stress-induced changes in the handling of aggregate-prone proteins in the ER-Golgi secretory pathway are discussed. This work was supported by grants from the National Institutes of Health (R15GM065139) and the National Science Foundation (DBI-0452587).

NMR order parameters calculated in an expanding reference frame: identifying sites of short- and long-range motion

NMR order parameters are calculated from molecular dynamics computer simulations of ubiquitin and the apo (Ca(2+)-free) state of calbindin D(9k). Calculations are performed in an expanding reference frame so as to discriminate between the effects of short- and long-range motions. This approach reveals that the dominant contributions to the order parameters are short-range. Longer-range contributions are limited to specific sites, many of which have been recognized in previous studies of correlated motions. These sites are identified on the basis of an effective reorientational number, n ( eff ). Not only does this parameter identify sites of short- and long-range motion, it also provides a way of evaluating the separability condition that is key to the Lipari-Szabo model-free method. When analyzed in conjunction with the Prompers-Brüschweiler separability index, the n ( eff ) values indicate that longer-range motions play a more prominent role in apo calbindin than they do in ubiquitin.

Intracellular pH modulates inner segment calcium homeostasis in vertebrate photoreceptors

Neuronal metabolic and electrical activity is associated with shifts in intracellular pH (pH(i)) proton activity and state-dependent changes in activation of signaling pathways in the plasma membrane, cytosol, and intracellular compartments. We investigated interactions between two intracellular messenger ions, protons and calcium (Ca²(+)), in salamander photoreceptor inner segments loaded with Ca²(+) and pH indicator dyes. Resting cytosolic pH in rods and cones in HEPES-based saline was acidified by ∼0.4 pH units with respect to pH of the superfusing saline (pH = 7.6), indicating that dissociated inner segments experience continuous acid loading. Cytosolic alkalinization with ammonium chloride (NH₄Cl) depolarized photoreceptors and stimulated Ca²(+) release from internal stores, yet paradoxically also evoked dose-dependent, reversible decreases in [Ca²(+)](i). Alkalinization-evoked [Ca²(+)](i) decreases were independent of voltage-operated and store-operated Ca²(+) entry, plasma membrane Ca²(+) extrusion, and Ca²(+) sequestration into internal stores. The [Ca²(+)](i)-suppressive effects of alkalinization were antagonized by the fast Ca²(+) buffer BAPTA, suggesting that pH(i) directly regulates Ca²(+) binding to internal anionic sites. In summary, this data suggest that endogenously produced protons continually modulate the membrane potential, release from Ca²(+) stores, and intracellular Ca²(+) buffering in rod and cone inner segments.

Binding of calcium and other metal ions to the EF-hand loops of calmodulin studied by quantum chemical calculations and molecular dynamics simulations

Calcium ion binding by the four EF-hand motifs of the protein calmodulin (CaM) is a central event in Ca2+-based cellular signaling. To understand molecular details of this complex process, isolated Ca2+-binding loops can be studied, by use of both experiments and calculations. In this work, we explore the metal specificity of the four Ca2+-binding loops of CaM using density functional theory (DFT) quantum chemical calculations and molecular dynamics simulations. We study CaM complexes with the physiologically important ions of calcium (Ca2+) and magnesium (Mg2+) and also with two other ions, strontium (Sr2+) and lanthanum (La3+). The former is of interest in the area of radioactive waste bioremediation, whereas the latter is often used as a probe of Ca2+-binding sites. We obtain intrinsic metal ion-loop binding energies as well as their components: vacuum, charge-transfer, solvation, entropy, and deformation terms. A detailed analysis of the results reveals that the total binding energy depends on a delicate balance among these energy components. They, in turn, are determined by the cation's charge and size as well as the amino acid composition and flexibility of the loops and the identity of the metal-chelating residues.

Effects of pH and ionic strength on the assembly and bundling of FtsZ protofilaments: A possible role of electrostatic interactions in the bundling of protofilaments

Assembly, bundling and stability of FtsZ protofilaments are important for the formation and functioning of the cytokinetic Z-ring during bacterial division. We found that the bundling of FtsZ protofilaments decreased strongly with increasing pH from 6.0 to 7.9, while the assembly of FtsZ monomers did not decrease considerably. In addition, the disassembly of FtsZ protofilaments was strongly suppressed at pH 6.0 as compared to the elevated pHs. The far-UV circular dichroism spectra of the native FtsZ and the tryptophan emission spectra of mutated FtsZ (Y371W) did not change by increasing pH from 6 to 7.9 indicating that the structure of FtsZ was not altered significantly. Further, the inhibition of bundling of FtsZ protofilaments predominantly, and the inhibition of assembly to a lesser extent by salt indicated that electrostatic interactions are important for the assembly and bundling of FtsZ protofilaments. These observations are supported by the results of computational docking of Escherichia coli dimer structure. The results suggest that the basic intracellular pH (7.4-7.8) of E. coli may play a role in regulating the assembly dynamics of FtsZ in the Z-ring by reducing protofilament stability and bundling in bacterial cytoplasm.

Intra- versus intermolecular interactions in monellin: contribution of surface charges to protein assembly

The relative significance of weak non-covalent interactions in biological context has been much debated. Here, we have addressed the contribution of Coulombic interactions to protein stability and assembly experimentally. The sweet protein monellin, a non-covalently linked heterodimeric protein, was chosen for this study because of its ability to spontaneously reconstitute from separated fragments. The reconstitution of monellin mutants containing large surface charge perturbations was compared to the thermostability of structurally equivalent single-chain monellin containing the same sets of mutations under varying salt concentrations. The affinity between monellin fragments is found to correlate with the thermostability of single chain monellin, indicating the involvement of the same underlying Coulombic interactions. This confirms that there are no principal differences in the interactions involved in folding and binding. Based on comparison with a previous mutational study involving hydrophobic core residues, the relative contribution of Coulombic interactions to stability and affinity is modest. However, the Coulombic perturbations only affect the association rates of reconstitution in contrast to perturbations involving hydrophobic residues, which affect primarily the dissociation rates. These results indicate that Coulombic interactions are likely to be of main importance for the association of protein assembly, relevant for functions of proteins.

Solution structure and backbone dynamics of Calsensin, an invertebrate neuronal calcium-binding protein

Calsensin is an EF-hand calcium-binding protein expressed by a subset of peripheral sensory neurons that fasciculate into a single tract in the leech central nervous system. Calsensin is a 9-kD protein with two EF-hand calcium-binding motifs. Using multidimensional NMR spectroscopy we have determined the solution structure and backbone dynamics of calcium-bound Calsensin. Calsensin consists of four helices forming a unicornate-type four-helix bundle. The residues in the third helix undergo slow conformational exchange indicating that the motion of this helix is associated with calciumbinding. The backbone dynamics of the protein as measured by (15)N relaxation rates and heteronuclear NOEs correlate well with the three-dimensional structure. Furthermore, comparison of the structure of Calsensin with other members of the EF-hand calcium-binding protein family provides insight into plausible mechanisms of calcium and target protein binding.

Characterization of calretinin I-II as an EF-hand, Ca2+, H+-sensing domain

Calretinin, a neuronal protein with well-defined calcium-binding properties, has a poorly defined function. The pH dependent properties of calretinin (CR), the N-terminal (CR I-II), and C-terminal (CR III-VI) domains were investigated. A drop in pH within the intracellular range (from pH 7.5 to pH 6.5) leads to an increased hydrophobicity of calcium-bound CR and its domains as reported by fluorescence spectroscopy with the hydrophobic probe 2-(p-toluidino)-6-naphthalenesulfonic acid (TNS). The TNS data for the N- and C-terminal domains of CR are additive, providing further support for their independence within the full-length protein. Our work concentrated on CR I-II, which was found to have hydrophobic properties similar to calmodulin at lower pH. The elution of CR I-II from a phenyl-Sepharose column was consistent with the TNS data. The pH-dependent structural changes were further localized to residues 13-28 and 44-51 using nuclear magnetic resonance spectroscopy chemical shift analysis, and there appear to be no large changes in secondary structure. Protonation of His 12 and/or His 27 side chains, coupled with calcium chelation, appears to lead to the organization of a hydrophobic pocket in the N-terminal domain. CR may sense and respond to calcium, proton, and other signals, contributing to conflicting data on the proteins role as a calcium sensor or calcium buffer.

Electrostatic contributions to the kinetics and thermodynamics of protein assembly

The role of electrostatic interactions in the assembly of a native protein structure was studied using fragment complementation. Contributions of salt, pH, or surface charges to the kinetics and equilibrium of calbindin D(9k) reconstitution was measured in the presence of Ca(2+) using surface plasmon resonance and isothermal titration calorimetry. Whereas surface charge substitutions primarily affect the dissociation rate constant, the association rates are correlated with subdomain net charge in a way expected for Coulomb interactions. The affinity is reduced in all mutants, with the largest effect (260-fold) observed for the double mutant K25E+K29E. At low net charge, detailed charge distribution is important, and charges remote from the partner EF-hand have less influence than close ones. The effects of salt and pH on the reconstitution are smaller than mutational effects. The interaction between the wild-type EF-hands occurs with high affinity (K(A) = 1.3 x 10(10) M(-1); K(D) = 80 pM). The enthalpy of association is overall favorable and there appears to be a very large favorable entropic contribution from the desolvation of hydrophobic surfaces that become buried in the complex. Electrostatic interactions contribute significantly to the affinity between the subdomains, but other factors, such as hydrophobic interactions, dominate.