Statistical analysis of structural characteristics of protein Ca2+-binding sites

High-Throughput Identification of Calcium Regulated Proteins Across Diverse Proteomes

Calcium ions play important roles in nearly every biological process, yet whole-proteome analysis of calcium effectors has been hindered by lack of high-throughput, unbiased, and quantitative methods to identify proteins-calcium engagement. To address this, we adapted protein thermostability assays in the budding yeast, human cells, and mouse mitochondria. Based on calcium-dependent thermostability, we identified 2884 putative calcium-regulated proteins across human, mouse, and yeast proteomes. These data revealed calcium engagement of novel signaling hubs and cellular processes, including metabolic enzymes and the spliceosome. Cross-species comparison of calcium-protein engagement and mutagenesis experiments identified residue-specific cation engagement, even within well-known EF-hand domains. Additionally, we found that the dienoyl-CoA reductase DECR1 binds calcium at physiologically-relevant concentrations with substrate-specific affinity, suggesting direct calcium regulation of mitochondrial fatty acid oxidation. These unbiased, proteomic analyses of calcium effectors establish a key resource to dissect cation engagement and its mechanistic effects across multiple species and diverse biological processes.

Drastic alterations in the loop structure around colchicine upon complex formation with an engineered lipocalin indicate a conformational selection mechanism

Using Anticalin technology, a lipocalin protein dubbed Colchicalin, with the ability to bind the toxic plant alkaloid colchicine with picomolar affinity, has previously been engineered, thus offering a potential antidote in vivo and also allowing its sensitive detection in biological samples. To further analyze the mode of ligand recognition, the crystal structure of Colchicalin is now reported in its unliganded form and is compared with the colchicine complex. A superposition of the protein structures revealed major rearrangements in the four structurally variable loops of the engineered lipocalin. Notably, the binding pocket in the unbound protein is largely occupied by the inward-bent loop #3, in particular Ile97, as well as by the phenylalanine side chain at position 71 in loop #2. Upon binding of colchicine, a dramatic shift of loop #3 by up to 11.1 Å occurs, in combination with a side-chain flip of Phe71, thus liberating the necessary space within the ligand pocket. Interestingly, the proline residue at the neighboring position 72, which arose during the combinatorial engineering of Colchicalin, remained in a cis configuration in both structures. These findings provide a striking example of a conformational adaptation mechanism, which is a long-known phenomenon for antibodies in immunochemistry, during the recognition of a small ligand by an engineered lipocalin, thus illustrating the general similarity between the mode of antigen/ligand binding by immunoglobulins and lipocalins.

Machine learning-based modulation of Ca2+-binding affinity in EF-hand proteins and comparative structural insights into site-specific cooperative binding

Ca2+-binding proteins are present in almost all living organisms and different types display different levels of binding affinities for the cation. Here, we report two new scoring schemes enabling the user to estimate and manipulate the calcium binding affinities in EF hand containing proteins. To validate this, we designed a unique EF-hand loop capable of binding calcium with high affinity by altering five residues. The N-terminal domain of Entamoeba histolytica calcium-binding protein1 (NtEhCaBP1) is used for site-directed mutagenesis to incorporate the designed loop sequence into the second EF-hand motif of this protein, referred as Nt-EhCaBP1-EF2 mutant. The binding isotherms calculated using ITC calorimetry showed that Nt-EhCaBP1-EF2 mutant site binds Ca2+ with higher affinity than Wt-Nt-EhCaBP1, by ∼600 times. The crystal structure of the mutant displayed more compact Ca2+-coordination spheres in both of its EF loops than the structure of the wildtype protein. The compact coordination sphere of EF-2 causes the bend in the helix-3, which leads to the formation of unexpected hexamer of NtEhCaBP1-EF2 mutant structure. Further dynamic correlation analysis revealed that the mutation in the second EF loop changed the entire residue network of the monomer, resulting in stronger coordination of Ca2+ even in another EF-hand loop.

Self-assembled peptide P11-4 interacts with the type I collagen C-terminal telopeptide domain and calcium ions

OBJECTIVES

Evaluate molecularly the role of P11-4 self-assembly peptide in dentin remineralization and its interaction with collagen I.

METHODS

The calcium-responsive P11-4 peptide was analyzed by intrinsic fluorescence emission spectrum, circular dichroism spectrum (CD), and atomic force microscope (AFM). Differential light scattering was used to monitor the nucleation growth rate of calcium phosphate nanocrystals in the absence or in the presence of P11-4. AFM was used to analyze the radial size (nm) of calcium phosphate nanocrystals formed in the absence or in the presence of P11-4, as well as to verify the spatial structure of P11-4 in the absence or in the presence of Ca2+.

RESULTS

The interaction of Ca2+ with the P11-4 (KD = 0.58 ± 0.06 mM) promotes the formation of β-sheet antiparallel structure, leads to its precipitation in saturated solutions of Ca/P = 1.67 and induces the formation of parallel large fibrils (0.6 - 1.5 µm). P11-4 organized the HAP nucleation by reducing both the growth rate and size variability of nanocrystals, analyzed by the F test (p < 0.0001, N = 30). P11-4 interacts (KD = 0.75 ± 0.06 μM) with the KGHRGFSGL motif present at the C-terminal collagen telopeptide domain. P11-4 also increased the amount of HAP and collagen in the MDPC-23 cells.

SIGNIFICANCE

The presented data propose a mechanism that will help future clinical and/or basic research to better understand a molecule able to inhibit structural collagen loss and help the impaired tissue to remineralize.

Mechanistic insights into the role of calcium in the allosteric regulation of the calmodulin-regulated death-associated protein kinase

Calcium (Ca²⁺) signaling plays an important role in the regulation of many cellular functions. Ca²⁺-binding protein calmodulin (CaM) serves as a primary effector of calcium function. Ca²⁺/CaM binds to the death-associated protein kinase 1 (DAPK1) to regulate intracellular signaling pathways. However, the mechanism underlying the influence of Ca²⁺ on the conformational dynamics of the DAPK1-CaM interactions is still unclear. Here, we performed large-scale molecular dynamics (MD) simulations of the DAPK1-CaM complex in the Ca²⁺-bound and-unbound states to reveal the importance of Ca²⁺. MD simulations revealed that removal of Ca²⁺ increased the anti-correlated inter-domain motions between DAPK1 and CaM, which weakened the DAPK1-CaM interactions. Binding free energy calculations validated the decreased DAPK1-CaM interactions in the Ca²⁺-unbound state. Structural analysis further revealed that Ca²⁺ removal caused the significant conformational changes at the DAPK1-CaM interface, especially the helices α1, α2, α4, α6, and α7 from the CaM and the basic loop and the phosphate-binding loop from the DAPK1. These results may be useful to understand the biological role of Ca²⁺ in physiological processes.

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

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.

Bioinorganic insights of the PQQ-dependent alcohol dehydrogenases

Among the several alcohol dehydrogenases, PQQ-dependent enzymes are mainly found in the α, β, and γ-proteobacteria. These proteins are classified into three main groups. Type I ADHs are localized in the periplasm and contain one Ca²⁺-PQQ moiety, being the methanol dehydrogenase (MDH) the most representative. In recent years, several lanthanide-dependent MDHs have been discovered exploding the understanding of the natural role of lanthanide ions. Type II ADHs are localized in the periplasm and possess one Ca²⁺-PQQ moiety and one heme c group. Finally, type III ADHs are complexes of two or three subunits localized in the cytoplasmic membrane and possess one Ca²⁺-PQQ moiety and four heme c groups, and in one of these proteins, an additional [2Fe-2S] cluster has been discovered recently. From the bioinorganic point of view, PQQ-dependent alcohol dehydrogenases have been revived recently mainly due to the discovery of the lanthanide-dependent enzymes. Here, we review the three types of PQQ-dependent ADHs with special focus on their structural features and electron transfer processes. The PQQ-Alcohol dehydrogenases are classified into three main groups. Type I and type II ADHs are located in the periplasm, while type III ADHs are in the cytoplasmic membrane. ADH-I have a Ca-PQQ or a Ln-PQQ, ADH-II a Ca-PQQ and one heme-c and ADH-III a Ca-PQQ and four hemes-c. This review focuses on their structural features and electron transfer processes.

Structural Aspects and Prediction of Calmodulin-Binding Proteins

Calmodulin (CaM) is an important intracellular protein that binds Ca²⁺ and functions as a critical second messenger involved in numerous biological activities through extensive interactions with proteins and peptides. CaM’s ability to adapt to binding targets with different structures is related to the flexible central helix separating the N- and C-terminal lobes, which allows for conformational changes between extended and collapsed forms of the protein. CaM-binding targets are most often identified using prediction algorithms that utilize sequence and structural data to predict regions of peptides and proteins that can interact with CaM. In this review, we provide an overview of different CaM-binding proteins, the motifs through which they interact with CaM, and shared properties that make them good binding partners for CaM. Additionally, we discuss the historical and current methods for predicting CaM binding, and the similarities and differences between these methods and their relative success at prediction. As new CaM-binding proteins are identified and classified, we will gain a broader understanding of the biological processes regulated through changes in Ca²⁺ concentration through interactions with CaM.

International Union of Basic and Clinical Pharmacology. CVIII. Calcium-Sensing Receptor Nomenclature, Pharmacology, and Function

The calcium-sensing receptor (CaSR) is a class C G protein-coupled receptor that responds to multiple endogenous agonists and allosteric modulators, including divalent and trivalent cations, L-amino acids, γ-glutamyl peptides, polyamines, polycationic peptides, and protons. The CaSR plays a critical role in extracellular calcium (Ca2+ o) homeostasis, as demonstrated by the many naturally occurring mutations in the CaSR or its signaling partners that cause Ca2+ o homeostasis disorders. However, CaSR tissue expression in mammals is broad and includes tissues unrelated to Ca2+ o homeostasis, in which it, for example, regulates the secretion of digestive hormones, airway constriction, cardiovascular effects, cellular differentiation, and proliferation. Thus, although the CaSR is targeted clinically by the positive allosteric modulators (PAMs) cinacalcet, evocalcet, and etelcalcetide in hyperparathyroidism, it is also a putative therapeutic target in diabetes, asthma, cardiovascular disease, and cancer. The CaSR is somewhat unique in possessing multiple ligand binding sites, including at least five putative sites for the "orthosteric" agonist Ca2+ o, an allosteric site for endogenous L-amino acids, two further allosteric sites for small molecules and the peptide PAM, etelcalcetide, and additional sites for other cations and anions. The CaSR is promiscuous in its G protein-coupling preferences, and signals via Gq/11, Gi/o, potentially G12/13, and even Gs in some cell types. Not surprisingly, the CaSR is subject to biased agonism, in which distinct ligands preferentially stimulate a subset of the CaSR's possible signaling responses, to the exclusion of others. The CaSR thus serves as a model receptor to study natural bias and allostery. SIGNIFICANCE STATEMENT: The calcium-sensing receptor (CaSR) is a complex G protein-coupled receptor that possesses multiple orthosteric and allosteric binding sites, is subject to biased signaling via several different G proteins, and has numerous (patho)physiological roles. Understanding the complexities of CaSR structure, function, and biology will aid future drug discovery efforts seeking to target this receptor for a diversity of diseases. This review summarizes what is known to date regarding key structural, pharmacological, and physiological features of the CaSR.

Design of Calcium-Binding Proteins to Sense Calcium

Calcium controls numerous biological processes by interacting with different classes of calcium binding proteins (CaBP’s), with different affinities, metal selectivities, kinetics, and calcium dependent conformational changes. Due to the diverse coordination chemistry of calcium, and complexity associated with protein folding and binding cooperativity, the rational design of CaBP’s was anticipated to present multiple challenges. In this paper we will first discuss applications of statistical analysis of calcium binding sites in proteins and subsequent development of algorithms to predict and identify calcium binding proteins. Next, we report efforts to identify key determinants for calcium binding affinity, cooperativity and calcium dependent conformational changes using grafting and protein design. Finally, we report recent advances in designing protein calcium sensors to capture calcium dynamics in various cellular environments.

Crystal structure of proteolyzed VapBC and DNA-bound VapBC from Salmonella enterica Typhimurium LT2 and VapC as a putative Ca2+ -dependent ribonuclease

Bacterial toxin-antitoxin (TA) system has gained attention for its essential roles in cellular maintenance and survival under harsh environmental conditions such as nutrient deficiency and antibiotic treatment. There are at least 14 TA systems in Salmonella enterica serovar Typhimurium LT2, a pathogenic bacterium, and none of the structures of these TA systems have been determined. We determined the crystal structure of the VapBC TA complex from S. Typhimurium LT2 in proteolyzed and DNA-bound forms at 2.0 Å and 2.8 Å resolution, respectively. The VapC toxin possesses a pilT N-terminal domain (PIN-domain) that shows ribonuclease activity, and the VapB antitoxin has an AbrB-type DNA binding domain. In addition, the structure revealed details of interaction mode between VapBC and the cognate promoter DNA, including the inhibition of VapC by VapB and linear conformation of bound DNA in the VapBC complex. The complexation of VapBC with the linear DNA is not consistent with known structures of VapBC homologs in complex with bent DNA. We also identified VapC from S. Typhimurium LT2 as a putative Ca²⁺ -dependent ribonuclease, which differs from previous data showing that VapC homologs have Mg²⁺ or Mn²⁺ -dependent ribonuclease activities. The present studies could provide structural understanding of the physiology of VapBC systems and foundation for the development of new antibiotic drugs against Salmonella infection.

Structure of the Human ACP-ISD11 Heterodimer

In recent years, the mammalian mitochondrial protein complex for iron-sulfur cluster assembly has been the focus of important studies. This is partly because of its high degree of relevance in cell metabolism and because mutations of the involved proteins are the cause of several human diseases. Cysteine desulfurase NFS1 is the key enzyme of the complex. At present, it is well-known that the active form of NFS1 is stabilized by the small protein ISD11. In this work, the structure of the human mitochondrial ACP-ISD11 heterodimer was determined at 2.0 Å resolution. ACP-ISD11 forms a cooperative unit stabilized by several ionic interactions, hydrogen bonds, and apolar interactions. The 4'-phosphopantetheine-acyl chain, which is covalently bound to ACP, interacts with several residues of ISD11, modulating together with ACP the foldability of ISD11. Recombinant human ACP-ISD11 was able to interact with the NFS1 desulfurase, thus yielding an active enzyme, and the NFS1/ACP-ISD11 core complex was activated by frataxin and ISCU proteins. Internal motions of ACP-ISD11 were studied by molecular dynamics simulations, showing the persistence of the interactions between both protein chains. The conformation of the dimer is similar to that found in the context of the (NFS1/ACP-ISD11)₂ supercomplex core, which contains the Escherichia coli ACP instead of the human variant. This fact suggests a sequential mechanism for supercomplex consolidation, in which the ACP-ISD11 complex may fold independently and, after that, the NFS1 dimer would be stabilized.

A theoretical and experimental approach to evaluate zein-calcium interaction in nixtamalization process

The possible interactions between α-zein and Ca²⁺ in nixtamalization process were analyzed from a multidisciplinary approach, considering the effect of these interactions on the thermal properties of the nixtamalized flour. SDS-PAGE under reducing and non-reducing conditions did not reveal differences between patterns of zeins from nixtamalized and control samples. However, analysis from affinity capillary electrophoresis indicated an increment in protein volume when calcium is added to zein extracted from nixtamalized flour. In addition, the binding constant for the zein-calcium interaction was calculated indicating a higher affinity for calcium by zein from nixtamalized samples. Molecular dynamics simulations indicated that the interaction α-zein-Ca²⁺ through C-ter was more favorable than Glu48. However, in excess of Ca²⁺ ions, each site could bind one calcium atom at the same time, confirming that aggregation of α-zein through calcium bridges is possible, expanding the technological applications of this protein.

Designing Calcium-Binding Proteins for Molecular MR Imaging

Early diagnosis, noninvasive detection, and staging of various diseases, remain one of the major clinical barriers to effective medical treatment and prevention of disease progression toward major clinical consequences. Molecular imaging technologies play an indispensable role in the clinical field in overcoming these major barriers. The increasing application of imaging techniques and agents in early detection of different diseases such as cancer has resulted in improved treatment response and clinical patient management. In this chapter we will first introduce criteria for the design and engineering of calcium-binding protein (CaBP) parvalbumin as a protein Gd-MRI contrast agent (ProCA) with unprecedented metal selectivity for Gd³⁺ over physiological metal ions. We will then discuss the further development of targeted MRI contrast agent for molecular imaging of PSMA biomarker for early detection of prostate cancer.

Calcium-Binding Proteins with Disordered Structure and Their Role in Secretion, Storage, and Cellular Signaling

Calcium is one of the most important second messengers and its intracellular signaling regulates many aspects of cell physiology. Calcium ions, like phosphate ions, are highly charged and thus are able to alter protein conformation upon binding; thereby they constitute key factors in signal transduction. One of the most common calcium-binding structural motifs is the EF-hand, a well-defined helix-loop-helix structural domain, present in many calcium-binding proteins (CBPs). Nonetheless, some CBPs contain non-canonical, disordered motifs, which usually bind calcium with high capacity and low affinity, and which represent a subset of proteins with specific functions, but these functions rarely involve signaling. When compared with phosphorylation-mediated signal transduction, the role of intrinsic disorder in calcium signaling is significantly less prominent and not direct. The list of known examples of intrinsically disordered CBPs is relatively short and the disorder in these examples seems to be linked to secretion and storage. Calcium-sensitive phosphatase calcineurin is an exception, but it represents an example of transient disorder, which is, nevertheless, vital to the functioning of this protein. The underlying reason for the different role of disordered proteins in the two main cellular signaling systems appears to be linked to the gradient of calcium concentration, present in all living cells.

Defining potential roles of Pb(2+) in neurotoxicity from a calciomics approach

Metal ions play crucial roles in numerous biological processes, facilitating biochemical reactions by binding to various proteins. An increasing body of evidence suggests that neurotoxicity associated with exposure to nonessential metals (e.g., Pb(2+)) involves disruption of synaptic activity, and these observed effects are associated with the ability of Pb(2+) to interfere with Zn(2+) and Ca(2+)-dependent functions. However, the molecular mechanism behind Pb(2+) toxicity remains a topic of debate. In this review, we first discuss potential neuronal Ca(2+) binding protein (CaBP) targets for Pb(2+) such as calmodulin (CaM), synaptotagmin, neuronal calcium sensor-1 (NCS-1), N-methyl-d-aspartate receptor (NMDAR) and family C of G-protein coupled receptors (cGPCRs), and their involvement in Ca(2+)-signalling pathways. We then compare metal binding properties between Ca(2+) and Pb(2+) to understand the structural implications of Pb(2+) binding to CaBPs. Statistical and biophysical studies (e.g., NMR and fluorescence spectroscopy) of Pb(2+) binding are discussed to investigate the molecular mechanism behind Pb(2+) toxicity. These studies identify an opportunistic, allosteric binding of Pb(2+) to CaM, which is distinct from ionic displacement. Together, these data suggest three potential modes of Pb(2+) activity related to molecular and/or neural toxicity: (i) Pb(2+) can occupy Ca(2+)-binding sites, inhibiting the activity of the protein by structural modulation, (ii) Pb(2+) can mimic Ca(2+) in the binding sites, falsely activating the protein and perturbing downstream activities, or (iii) Pb(2+) can bind outside of the Ca(2+)-binding sites, resulting in the allosteric modulation of the protein activity. Moreover, the data further suggest that even low concentrations of Pb(2+) can interfere at multiple points within the neuronal Ca(2+) signalling pathways to cause neurotoxicity.

Calcium Dynamics Mediated by the Endoplasmic/Sarcoplasmic Reticulum and Related Diseases

The flow of intracellular calcium (Ca2+) is critical for the activation and regulation of important biological events that are required in living organisms. As the major Ca2+ repositories inside the cell, the endoplasmic reticulum (ER) and the sarcoplasmic reticulum (SR) of muscle cells are central in maintaining and amplifying the intracellular Ca2+ signal. The morphology of these organelles, along with the distribution of key calcium-binding proteins (CaBPs), regulatory proteins, pumps, and receptors fundamentally impact the local and global differences in Ca2+ release kinetics. In this review, we will discuss the structural and morphological differences between the ER and SR and how they influence localized Ca2+ release, related diseases, and the need for targeted genetically encoded calcium indicators (GECIs) to study these events.

Molecular Basis for Modulation of Metabotropic Glutamate Receptors and Their Drug Actions by Extracellular Ca2

Metabotropic glutamate receptors (mGluRs) associated with the slow phase of the glutamatergic signaling pathway in neurons of the central nervous system have gained importance as drug targets for chronic neurodegenerative diseases. While extracellular Ca²⁺ was reported to exhibit direct activation and modulation via an allosteric site, the identification of those binding sites was challenged by weak binding. Herein, we review the discovery of extracellular Ca²⁺ in regulation of mGluRs, summarize the recent developments in probing Ca²⁺ binding and its co-regulation of the receptor based on structural and biochemical analysis, and discuss the molecular basis for Ca²⁺ to regulate various classes of drug action as well as its importance as an allosteric modulator in mGluRs.

Molecular Basis of S100A1 Activation at Saturating and Subsaturating Calcium Concentrations

The S100A1 protein mediates a wide variety of physiological processes through its binding of calcium (Ca(2+)) and endogenous target proteins. S100A1 presents two Ca(2+)-binding domains: a high-affinity "canonical" EF (cEF) hand and a low-affinity "pseudo" EF (pEF) hand. Accumulating evidence suggests that both Ca(2+)-binding sites must be saturated to stabilize an open state conducive to peptide recognition, yet the pEF hand's low affinity limits Ca(2+) binding at normal physiological concentrations. To understand the molecular basis of Ca(2+) binding and open-state stabilization, we performed 100 ns molecular dynamics simulations of S100A1 in the apo/holo (Ca(2+)-free/bound) states and a half-saturated state, for which only the cEF sites are Ca(2+)-bound. Our simulations indicate that the pattern of oxygen coordination about Ca(2+) in the cEF relative to the pEF site contributes to the former's higher affinity, whereas Ca(2+) binding strongly reshapes the protein's conformational dynamics by disrupting β-sheet coupling between EF hands. Moreover, modeling of the half-saturated configuration suggests that the open state is unstable and reverts toward a closed state in the absence of the pEF Ca(2+) ion. These findings indicate that Ca(2+) binding at the cEF site alone is insufficient to stabilize opening; thus, posttranslational modification of the protein may be required for target peptide binding at subsaturating intracellular Ca(2+) levels.

Molecular Basis of the Extracellular Ligands Mediated Signaling by the Calcium Sensing Receptor

Ca²⁺-sensing receptors (CaSRs) play a central role in regulating extracellular calcium concentration ([Ca²⁺]_o) homeostasis and many (patho)physiological processes in multiple organs. This regulation is orchestrated by a cooperative response to extracellular stimuli such as small changes in Ca²⁺, Mg²⁺, amino acids, and other ligands. In addition, CaSR is a pleiotropic receptor regulating several intracellular signaling pathways, including calcium mobilization and intracellular calcium oscillation. Nearly 200 mutations and polymorphisms have been found in CaSR in relation to a variety of human disorders associated with abnormal Ca²⁺ homeostasis. In this review, we summarize efforts directed at identifying binding sites for calcium and amino acids. Both homotropic cooperativity among multiple calcium binding sites and heterotropic cooperativity between calcium and amino acid were revealed using computational modeling, predictions, and site-directed mutagenesis coupled with functional assays. The hinge region of the bilobed Venus flytrap (VFT) domain of CaSR plays a pivotal role in coordinating multiple extracellular stimuli, leading to cooperative responses from the receptor. We further highlight the extensive number of disease-associated mutations that have also been shown to affect CaSR's cooperative action via several types of mechanisms. These results provide insights into the molecular bases of the structure and functional cooperativity of this receptor and other members of family C of the G protein-coupled receptors (cGPCRs) in health and disease states, and may assist in the prospective development of novel receptor-based therapeutics.

UCP3 is associated with Hax-1 in mitochondria in the presence of calcium ion

Uncoupling protein 3 (UCP3) is known to regulate energy dissipation, proton leakage, fatty acid oxidation, and oxidative stress. To identify the putative protein regulators of UCP3, we performed yeast two-hybrid screens. Here we report that UCP3 interacted with HS-1 associated protein X-1 (Hax-1), an anti-apoptotic protein that was localized in the mitochondria, and is involved in cellular responses to Ca(2+). The hydrophilic sequences within loop 2, and the matrix-localized hydrophilic domain of mouse UCP3, were necessary for binding to Hax-1 at the C-terminal domain, adjacent to the mitochondrial inner membrane. Interestingly, interaction of these proteins occurred in a calcium-dependent manner. Moreover, the NMR spectrum of the C-terminal domain of Hax-1 was dramatically changed by removal of Ca(2+), suggesting that the C-terminal domain of Hax-1 underwent a Ca(2+)-induced conformational change. In the Ca(2+)-free state, the C-terminal Hax-1 tended to unfold, suggesting that Ca(2+) binding may induce protein folding of the Hax-1 C-terminus. These results suggested that the UCP3-Hax-1 complex may regulate mitochondrial functional changes caused by mitochondrial Ca(2+).

An electrostatic mechanism for Ca(2+)-mediated regulation of gap junction channels

Gap junction channels mediate intercellular signalling that is crucial in tissue development, homeostasis and pathologic states such as cardiac arrhythmias, cancer and trauma. To explore the mechanism by which Ca(2+) blocks intercellular communication during tissue injury, we determined the X-ray crystal structures of the human Cx26 gap junction channel with and without bound Ca(2+). The two structures were nearly identical, ruling out both a large-scale structural change and a local steric constriction of the pore. Ca(2+) coordination sites reside at the interfaces between adjacent subunits, near the entrance to the extracellular gap, where local, side chain conformational rearrangements enable Ca(2+)chelation. Computational analysis revealed that Ca(2+)-binding generates a positive electrostatic barrier that substantially inhibits permeation of cations such as K(+) into the pore. Our results provide structural evidence for a unique mechanism of channel regulation: ionic conduction block via an electrostatic barrier rather than steric occlusion of the channel pore.

Fast kinetics of calcium signaling and sensor design

Fast calcium signaling is regulated by numerous calcium channels exhibiting high spatiotemporal profiles which are currently measured by fluorescent calcium sensors. There is still a strong need to improve the kinetics of genetically encoded calcium indicators (sensors) to capture calcium dynamics in the millisecond time frame. In this review, we summarize several major fast calcium signaling pathways and discuss the recent developments and application of genetically encoded calcium indicators to detect these pathways. A new class of genetically encoded calcium indicators designed with site-directed mutagenesis on the surface of beta-barrel fluorescent proteins to form a pentagonal bipyramidal-like calcium binding domain dramatically accelerates calcium binding kinetics. Furthermore, novel genetically encoded calcium indicators with significantly increased fluorescent lifetime change are advantageous in deep-field imaging with high light-scattering and notable morphology change.

Myoplasmic resting Ca2+ regulation by ryanodine receptors is under the control of a novel Ca2+-binding region of the receptor

Passive SR (sarcoplasmic reticulum) Ca²⁺ leak through the RyR (ryanodine receptor) plays a critical role in the mechanisms that regulate [Ca²⁺]_rest (intracellular resting myoplasmic free Ca²⁺ concentration) in muscle. This process appears to be isoform-specific as expression of either RyR1 or RyR3 confers on myotubes different [Ca²⁺]_rest. Using chimaeric RyR3–RyR1 receptors expressed in dyspedic myotubes, we show that isoform-dependent regulation of [Ca²⁺]_rest is primarily defined by a small region of the receptor encompassing amino acids 3770–4007 of RyR1 (amino acids 3620–3859 of RyR3) named as the CLR (Ca²⁺ leak regulatory) region. [Ca²⁺]_rest regulation by the CLR region was associated with alteration of RyRs’ Ca²⁺-activation profile and changes in SR Ca²⁺-leak rates. Biochemical analysis using Tb³⁺-binding assays and intrinsic tryptophan fluorescence spectroscopy of purified CLR domains revealed that this determinant of RyRs holds a novel Ca²⁺-binding domain with conformational properties that are distinctive to each isoform. Our data suggest that the CLR region provides channels with unique functional properties that modulate the rate of passive SR Ca²⁺ leak and confer on RyR1 and RyR3 distinctive [Ca²⁺]_rest regulatory properties. The identification of a new Ca²⁺-binding domain of RyRs with a key modulatory role in [Ca²⁺]_rest regulation provides new insights into Ca²⁺-mediated regulation of RyRs.

This paper reports the finding of a new class of Ca²⁺-binding domain of intracellular Ca²⁺ channels from muscle cells. This domain provides channels with distinctive properties that result in channel-specific modulation of the intracellular resting Ca²⁺ concentration.

Design of ProCAs (protein-based Gd(3+) MRI contrast agents) with high dose efficiency and capability for molecular imaging of cancer biomarkers

Magnetic resonance imaging (MRI) is the leading imaging technique for disease diagnostics, providing high resolution, three-dimensional images noninvasively. MRI contrast agents are designed to improve the contrast and sensitivity of MRI. However, current clinically used MRI contrast agents have relaxivities far below the theoretical upper limit, which largely prevent advancing molecular imaging of biomarkers with desired sensitivity and specificity. This review describes current progress in the development of a new class of protein-based MRI contrast agents (ProCAs) with high relaxivity using protein design to optimize the parameters that govern relaxivity. Further, engineering with targeting moiety allows these contrast agents to be applicable for molecular imaging of prostate cancer biomarkers by MRI. The developed protein-based contrast agents also exhibit additional in vitro and in vivo advantages for molecular imaging of disease biomarkers, such as high metal-binding stability and selectivity, reduced toxicity, proper blood circulation time, and higher permeability in tumor tissue in addition to improved relaxivities.

Ca2+-binding motif of βγ-crystallins

βγ-Crystallin-type double clamp (N/D)(N/D)XX(S/T)S motif is an established but sparsely investigated motif for Ca(2+) binding. A βγ-crystallin domain is formed of two Greek key motifs, accommodating two Ca(2+)-binding sites. βγ-Crystallins make a separate class of Ca(2+)-binding proteins (CaBP), apparently a major group of CaBP in bacteria. Paralleling the diversity in βγ-crystallin domains, these motifs also show great diversity, both in structure and in function. Although the expression of some of them has been associated with stress, virulence, and adhesion, the functional implications of Ca(2+) binding to βγ-crystallins in mediating biological processes are yet to be elucidated.

Complex structure of type VI peptidoglycan muramidase effector and a cognate immunity protein

The type VI secretion system (T6SS) is a bacterial protein-export machine that is capable of delivering virulence effectors between Gram-negative bacteria. The T6SS of Pseudomonas aeruginosa transports two lytic enzymes, Tse1 and Tse3, to degrade cell-wall peptidoglycan in the periplasm of rival bacteria that are competing for niches via amidase and muramidase activities, respectively. Two cognate immunity proteins, Tsi1 and Tsi3, are produced by the bacterium to inactivate the two antibacterial effectors, thereby protecting its siblings from self-intoxication. Recently, Tse1-Tsi1 has been structurally characterized. Here, the structure of the Tse3-Tsi3 complex is reported at 1.9 Å resolution. The results reveal that Tse3 contains a C-terminal catalytic domain that adopts a soluble lytic transglycosylase (SLT) fold in which three calcium-binding sites were surprisingly observed close to the catalytic Glu residue. The electrostatic properties of the substrate-binding groove are also distinctive from those of known structures with a similar fold. All of these features imply that a unique catalytic mechanism is utilized by Tse3 in cleaving glycosidic bonds. Tsi3 comprises a single domain showing a β-sandwich architecture that is reminiscent of the immunoglobulin fold. Three loops of Tsi3 insert deeply into the groove of Tse3 and completely occlude its active site, which forms the structural basis of Tse3 inactivation. This work is the first crystallographic report describing the three-dimensional structure of the Tse3-Tsi3 effector-immunity pair.

Cd2+ as a Ca2+ surrogate in protein-membrane interactions: isostructural but not isofunctional

Due to its favorable spectroscopic properties, Cd(2+) is frequently used as a probe of Ca(2+) sites in proteins. We investigate the ability of Cd(2+) to act as a structural and functional surrogate of Ca(2+) in protein-membrane interactions. C2 domain from protein kinase Cα (C2α) was chosen as a paradigm for the Ca(2+)-dependent phosphatidylserine-binding peripheral membrane domains. We identified the Cd(2+)-binding sites of C2α using NMR spectroscopy, determined the 1.6 Å crystal structure of Cd(2+)-bound C2α, and characterized metal-ion-dependent interactions between C2α and phospholipid membranes using fluorescence spectroscopy and ultracentrifugation experiments. We show that Cd(2+) forms a tight complex with the membrane-binding loops of C2α but is unable to support its membrane-binding function. This is in sharp contrast with Pb(2+), which is almost as effective as Ca(2+) in driving the C2α-membrane association process. Our results provide the first direct evidence for the specific role of divalent metal ions in mediating protein-membrane interactions, have important implications for metal substitution studies in proteins, and illustrate the potential diversity of functional responses caused by toxic metal ions.

Calciomics: integrative studies of Ca2+-binding proteins and their interactomes in biological systems

Calcium ion (Ca(2+)), the fifth most common chemical element in the earth's crust, represents the most abundant mineral in the human body. By binding to a myriad of proteins distributed in different cellular organelles, Ca(2+) impacts nearly every aspect of cellular life. In prokaryotes, Ca(2+) plays an important role in bacterial movement, chemotaxis, survival reactions and sporulation. In eukaryotes, Ca(2+) has been chosen through evolution to function as a universal and versatile intracellular signal. Viruses, as obligate intracellular parasites, also develop smart strategies to manipulate the host Ca(2+) signaling machinery to benefit their own life cycles. This review focuses on recent advances in applying both bioinformatic and experimental approaches to predict and validate Ca(2+)-binding proteins and their interactomes in biological systems on a genome-wide scale (termed "calciomics"). Calmodulin is used as an example of Ca(2+)-binding protein (CaBP) to demonstrate the role of CaBPs on the regulation of biological functions. This review is anticipated to rekindle interest in investigating Ca(2+)-binding proteins and Ca(2+)-modulated functions at the systems level in the post-genomic era.

How calcium inhibits the magnesium-dependent kinase gsk3β: a molecular simulation study

Glycogen synthase kinase 3β (GSK3β) is a ubiquitous serine/threonine kinase that plays a pivotal role in many biological processes. GSK3β catalyzes the transfer of γ-phosphate of ATP to the unique substrate Ser/Thr residues with the assistance of two natural activating cofactors Mg(2+). Interestingly, the biological observation reveals that a non-native Ca(2+) ion can inhibit the GSK3β catalytic activity. Here, the inhibitory mechanism of GSK3β by the displacement of native Mg(2+) at site 1 by Ca(2+) was investigated by means of 80 ns comparative molecular dynamics (MD) simulations of the GSK3β···Mg(2+)-2/ATP/Mg(2+) -1 and GSK3β···Mg(2+)-2/ATP/Ca(2+)-1 systems. MD simulation results revealed that using the AMBER point charge model force field for Mg(2+) was more appropriate in the reproduction of the active site architectural characteristics of GSK3β than using the magnesium-cationic dummy atom model force field. Compared with the native Mg(2+) bound system, the misalignment of the critical triphosphate moiety of ATP, the erroneous coordination environments around the Mg(2+) ion at site 2, and the rupture of the key hydrogen bond between the invariant Lys85 and the ATP O(β2) atom in the Ca(2+) substituted system were observed in the MD simulation due to the Ca(2+) ion in active site in order to achieve its preferred sevenfold coordination geometry, which adequately abolish the enzymatic activity. The obtained results are valuable in understanding the possible mechanism by why Ca(2+) inhibits the GSK3β activity and also provide insights into the mechanism of Ca(2+) inhibition in other structurally related protein kinases.

Structural comparison of ColH and ColG collagen-binding domains from Clostridium histolyticum

Clostridium histolyticum secretes collagenases, ColG and ColH, that cause extensive tissue destruction in myonecrosis. The C-terminal collagen-binding domain (CBD) of collagenase is required for insoluble collagen fibril binding and subsequent collagenolysis. The high-resolution crystal structures of ColG-CBD (s3b) and ColH-CBD (s3) are reported in this paper. The new X-ray structure of s3 was solved at 2.0-Å resolution (R = 17.4%; R(free) = 23.3%), while the resolution of the previously determined s3b was extended to 1.4 Å (R = 17.9%; R(free) = 21.0%). Despite sharing only 30% sequence identity, the molecules resemble one another closely (root mean square deviation [RMSD] C(α) = 1.5 Å). All but one residue, whose side chain chelates with Ca(2+), are conserved. The dual Ca(2+) binding site in s3 is completed by an unconserved aspartate. Differential scanning calorimetric measurements showed that s3 gains thermal stability, comparable to s3b, by binding to Ca(2+) (holo T(m) = 94.1°C; apo T(m) = 70.2°C). holo s3 is also stabilized against chemical denaturants urea and guanidine HCl. The three most critical residues for collagen interaction in s3b are conserved in s3. The general shape of the binding pocket is retained by altered loop structures and side chain positions. Small-angle X-ray scattering data revealed that s3 also binds asymmetrically to minicollagen. Besides the calcium-binding sites and the collagen-binding pocket, architecturally important hydrophobic residues and the hydrogen-bonding network around the cis-peptide bond are well conserved within the metallopeptidase subfamily M9B. CBDs were previously shown to bind to the extracellular matrix of various tissues. Compactness and extreme stability in physiological Ca(2+) concentration possibly make both CBDs suitable for targeted growth factor delivery.

Detection of highly curved membrane surfaces using a cyclic peptide derived from synaptotagmin-I

The generation of highly curved membranes is essential to cell growth, division, and movement. Recent research in the field is focused to answer questions related to the consequences of changes in the topology of the membrane once it is created, broadly termed as membrane curvature sensing. Most probes that are used to study curvature sensing are intact membrane active proteins such as DP1/Yop1p, ArfGAP1, BAR domains, and Synaptotagmin-I (Syt1). Taking a cue from nature, we created the cyclic peptide C2BL3C based on the membrane penetration C2B loop 3 of Syt1 via "Click" chemistry. Using a combination of spectroscopic techniques, we investigated the peptide-lipid interactions of this peptide with synthetic phospholipid vesicles and exosomes from rat blood plasma. We found that the macrocycle peptide probe was selective for lipid vesicles with highly curved surfaces (d < 100 nm). These results suggested that C2BL3C functions as a selective detector of highly curved phospholipid bilayers.

Predicting Ca2+ -binding sites using refined carbon clusters

Identifying Ca(2+) -binding sites in proteins is the first step toward understanding the molecular basis of diseases related to Ca(2+) -binding proteins. Currently, these sites are identified in structures either through X-ray crystallography or NMR analysis. However, Ca(2+) -binding sites are not always visible in X-ray structures due to flexibility in the binding region or low occupancy in a Ca(2+) -binding site. Similarly, both Ca(2+) and its ligand oxygens are not directly observed in NMR structures. To improve our ability to predict Ca(2+) -binding sites in both X-ray and NMR structures, we report a new graph theory algorithm (MUG(C) ) to predict Ca(2+) -binding sites. Using carbon atoms covalently bonded to the chelating oxygen atoms, and without explicit reference to side-chain oxygen ligand co-ordinates, MUG(C) is able to achieve 94% sensitivity with 76% selectivity on a dataset of X-ray structures composed of 43 Ca(2+) -binding proteins. Additionally, prediction of Ca(2+) -binding sites in NMR structures was obtained by MUG(C) using a different set of parameters, which were determined by the analysis of both Ca(2+) -constrained and unconstrained Ca(2+) -loaded structures derived from NMR data. MUG(C) identified 20 of 21 Ca(2+) -binding sites in NMR structures inferred without the use of Ca(2+) constraints. MUG(C) predictions are also highly selective for Ca(2+) -binding sites as analyses of binding sites for Mg(2+) , Zn(2+) , and Pb(2+) were not identified as Ca(2+) -binding sites. These results indicate that the geometric arrangement of the second-shell carbon cluster is sufficient not only for accurate identification of Ca(2+) -binding sites in NMR and X-ray structures but also for selective differentiation between Ca(2+) and other relevant divalent cations.

Structure and properties of the Ca(2+)-binding CUB domain, a widespread ligand-recognition unit involved in major biological functions

CUB domains are 110-residue protein motifs exhibiting a β-sandwich fold and mediating protein-protein interactions in various extracellular proteins. Recent X-ray structural and mutagenesis studies have led to the identification of a particular CUB domain subset, cbCUB (Ca(2+)-binding CUB domain). Unlike other CUB domains, these harbour a homologous Ca(2+)-binding site that underlies a conserved binding site mediating ionic interaction between two of the three conserved acidic Ca(2+) ligands and a basic (lysine or arginine) residue of a protein ligand, similar to the interactions mediated by the low-density lipoprotein receptor family. cbCUB-mediated protein-ligand interactions usually involve multipoint attachment through several cbCUBs, resulting in high-affinity binding through avidity, despite the low affinity of individual interactions. The aim of the present review is to summarize our current knowledge about the structure and functions of cbCUBs, which represent the majority of the known CUB repertoire and are involved in a variety of major biological functions, including immunity and development, as well as in various cancer types. Examples discussed in the present review include a wide range of soluble and membrane-associated human proteins, as well as some archaeal and invertebrate proteins. The fact that these otherwise unrelated proteins share a common Ca(2+)-dependent ligand-binding ability suggests a mechanism inherited from very primitive ancestors. The information provided in the present review should stimulate further investigations on the crucial interactions mediated by cbCUB-containing proteins.

Analysis and prediction of calcium-binding pockets from apo-protein structures exhibiting calcium-induced localized conformational changes

Calcium binding in proteins exhibits a wide range of polygonal geometries that relate directly to an equally diverse set of biological functions. The binding process stabilizes protein structures and typically results in local conformational change and/or global restructuring of the backbone. Previously, we established the MUG program, which utilized multiple geometries in the Ca(2+)-binding pockets of holoproteins to identify such pockets, ignoring possible Ca(2+)-induced conformational change. In this article, we first report our progress in the analysis of Ca(2+)-induced conformational changes followed by improved prediction of Ca(2+)-binding sites in the large group of Ca(2+)-binding proteins that exhibit only localized conformational changes. The MUG(SR) algorithm was devised to incorporate side chain torsional rotation as a predictor. The output from MUG(SR) presents groups of residues where each group, typically containing two to five residues, is a potential binding pocket. MUG(SR) was applied to both X-ray apo structures and NMR holo structures, which did not use calcium distance constraints in structure calculations. Predicted pockets were validated by comparison with homologous holo structures. Defining a "correct hit" as a group of residues containing at least two true ligand residues, the sensitivity was at least 90%; whereas for a "correct hit" defined as a group of residues containing at least three true ligand residues, the sensitivity was at least 78%. These data suggest that Ca(2+)-binding pockets are at least partially prepositioned to chelate the ion in the apo form of the protein.

Elucidation of a novel extracellular calcium-binding site on metabotropic glutamate receptor 1{alpha} (mGluR1{alpha}) that controls receptor activation

Metabotropic glutamate receptor 1α (mGluR1α) exerts important effects on numerous neurological processes. Although mGluR1α is known to respond to extracellular Ca(2+) ([Ca(2+)](o)) and the crystal structures of the extracellular domains (ECDs) of several mGluRs have been determined, the calcium-binding site(s) and structural determinants of Ca(2+)-modulated signaling in the Glu receptor family remain elusive. Here, we identify a novel Ca(2+)-binding site in the mGluR1α ECD using a recently developed computational algorithm. This predicted site (comprising Asp-318, Glu-325, and Asp-322 and the carboxylate side chain of the receptor agonist, Glu) is situated in the hinge region in the ECD of mGluR1α adjacent to the reported Glu-binding site, with Asp-318 involved in both Glu and calcium binding. Mutagenesis studies indicated that binding of Glu and Ca(2+) to their distinct but partially overlapping binding sites synergistically modulated mGluR1α activation of intracellular Ca(2+) ([Ca(2+)](i)) signaling. Mutating the Glu-binding site completely abolished Glu signaling while leaving its Ca(2+)-sensing capability largely intact. Mutating the predicted Ca(2+)-binding residues abolished or significantly reduced the sensitivity of mGluR1α not only to [Ca(2+)](o) and [Gd(3+)](o) but also, in some cases, to Glu. The dual activation of mGluR1α by [Ca(2+)](o) and Glu has important implications for the activation of other mGluR subtypes and related receptors. It also opens up new avenues for developing allosteric modulators of mGluR function that target specific human diseases.

The structure of an Arf-ArfGAP complex reveals a Ca2+ regulatory mechanism

Arfs are small G proteins that have a key role in vesicle trafficking and cytoskeletal remodeling. ArfGAP proteins stimulate Arf intrinsic GTP hydrolysis by a mechanism that is still unresolved. Using a fusion construct we solved the structure of the ArfGAP ASAP3 in complex with Arf6 in the transition state. This structure clarifies the ArfGAP catalytic mechanism and shows a glutamine((Arf6)) and an arginine finger((ASAP3)) as the important catalytic residues. Unexpectedly the structure shows a calcium ion, liganded by both proteins in the complex interface, stabilizing the interaction and orienting the catalytic machinery. Calcium stimulates the GAP activity of ASAPs, but not other members of the ArfGAP family. This type of regulation is unique for GAPs and any other calcium-regulated processes and hints at a crosstalk between Ca(2+) and Arf signaling.

Integration of Diverse Research Methods to Analyze and Engineer Ca-Binding Proteins: From Prediction to Production

In recent years, increasingly sophisticated computational and bioinformatics tools have evolved for the analyses of protein structure, function, ligand interactions, modeling and energetics. This includes the development of algorithms to recursively evaluate side-chain rotamer permutations, identify regions in a 3D structure that meet some set of search parameters, calculate and minimize energy values, and provide high-resolution visual tools for theoretical modeling. Here we discuss the interdependency between different areas of bioinformatics, the evolution of different algorithm design approaches, and finally the transition from theoretical models to real-world design and application as they relate to Ca(2+)-binding proteins. Within this context, it has become evident that significant pre-experimental design and calculations can be modeled through computational methods, thus eliminating potentially unproductive research and increasing our confidence in the correlation between real and theoretical models. Moving from prediction to production, it is anticipated that bioinformatics tools will play an increasingly significant role in research and development, improving our ability to both understand the physiological roles of Ca(2+) and other metals and to extend that knowledge to the design of function-specific synthetic proteins capable of fulfilling different roles in medical diagnostics and therapeutics.

Inverse tuning of metal binding affinity and protein stability by altering charged coordination residues in designed calcium binding proteins

Ca(2+ )binding proteins are essential for regulating the role of Ca(2+ )in cell signaling and maintaining Ca(2+ )homeostasis. Negatively charged residues such as Asp and Glu are often found in Ca(2+ )binding proteins and are known to influence Ca(2+ )binding affinity and protein stability. In this paper, we report a systematic investigation of the role of local charge number and type of coordination residues in Ca(2+ )binding and protein stability using de novo designed Ca(2+ )binding proteins. The approach of de novo design was chosen to avoid the complications of cooperative binding and Ca(2+)-induced conformational change associated with natural proteins. We show that when the number of negatively charged coordination residues increased from 2 to 5 in a relatively restricted Ca(2+)-binding site, Ca(2+ )binding affinities increased by more than 3 orders of magnitude and metal selectivity for trivalent Ln(3+ )over divalent Ca(2+ )increased by more than 100-fold. Additionally, the thermal transition temperatures of the apo forms of the designed proteins decreased due to charge repulsion at the Ca(2+ )binding pocket. The thermal stability of the proteins was regained upon Ca(2+ )and Ln(3+ )binding to the designed Ca(2+ )binding pocket. We therefore observe a striking tradeoff between Ca(2+)/Ln(3+ )affinity and protein stability when the net charge of the coordination residues is varied. Our study has strong implications for understanding and predicting Ca(2+)-conferred thermal stabilization of natural Ca(2+ )binding proteins as well as for designing novel metalloproteins with tunable Ca(2+ )and Ln(3+ )binding affinity and selectivity.PACS codes: 05.10.-a.

Towards predicting Ca2+-binding sites with different coordination numbers in proteins with atomic resolution

Ca(2+)-binding sites in proteins exhibit a wide range of polygonal geometries that directly relate to an equally-diverse set of biological functions. Although the highly-conserved EF-Hand motif has been studied extensively, non-EF-Hand sites exhibit much more structural diversity which has inhibited efforts to determine the precise location of Ca(2+)-binding sites, especially for sites with few coordinating ligands. Previously, we established an algorithm capable of predicting Ca(2+)-binding sites using graph theory to identify oxygen clusters comprised of four atoms lying on a sphere of specified radius, the center of which was the predicted calcium position. Here we describe a new algorithm, MUG (MUltiple Geometries), which predicts Ca(2+)-binding sites in proteins with atomic resolution. After first identifying all the possible oxygen clusters by finding maximal cliques, a calcium center (CC) for each cluster, corresponding to the potential Ca(2+) position, is located to maximally regularize the structure of the (cluster, CC) pair. The structure is then inspected by geometric filters. An unqualified (cluster, CC) pair is further handled by recursively removing oxygen atoms and relocating the CC until its structure is either qualified or contains fewer than four ligand atoms. Ligand coordination is then determined for qualified structures. This algorithm, which predicts both Ca(2+) positions and ligand groups, has been shown to successfully predict over 90% of the documented Ca(2+)-binding sites in three datasets of highly-diversified protein structures with 0.22 to 0.49 A accuracy. All multiple-binding sites (i.e. sites with a single ligand atom associated with multiple calcium ions) were predicted, as were half of the low-coordination sites (i.e. sites with less than four protein ligand atoms) and 14/16 cofactor-coordinating sites. Additionally, this algorithm has the flexibility to incorporate surface water molecules and protein cofactors to further improve the prediction for low-coordination and cofactor-coordinating Ca(2+)-binding sites.