Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Two Structural Motifs within Canonical EF-Hand Calcium-Binding Domains Identify Five Different Classes of Calcium Buffers and Sensors

Abstract

Proteins with EF-hand calcium-binding motifs are essential for many cellular processes, but are also associated with cancer, autism, cardiac arrhythmias, and Alzheimer's, skeletal muscle and neuronal diseases. Functionally, all EF-hand proteins are divided into two groups: (1) calcium sensors, which function to translate the signal to various responses; and (2) calcium buffers, which control the level of free Ca2+ ions in the cytoplasm. The borderline between the two groups is not clear, and many proteins cannot be described as definitive buffers or sensors. Here, we describe two highly-conserved structural motifs found in all known different families of the EF-hand proteins. The two motifs provide a supporting scaffold for the DxDxDG calcium binding loop and contribute to the hydrophobic core of the EF hand domain. The motifs allow more precise identification of calcium buffers and calcium sensors. Based on the characteristics of the two motifs, we could classify individual EF-hand domains into five groups: (1) Open static; (2) Closed static; (3) Local dynamic; (4) Dynamic; and (5) Local static EF-hand domains.

Introduction

Calcium is essential for life [1] and plays at least two major roles in living organisms: structural and regulatory [2], [3]. Calcium is primarily found outside of cells, where it is complexed with phosphates or carbonates to form exo- and endoskeletons, serving as macro-scale structural scaffolds while also functioning to buffer the approximately 10−3 M extracellular levels of Ca2+ ions. In contrast, intracellular calcium concentrations are at least 104 times lower, and require control mechanisms to be maintained at the appropriately low levels. Any failure of these control mechanisms may lead to sustained calcium overload and eventual cell and organ malfunction. The disparity in extracellular and cytoplasmic concentrations of Ca2+ ions supports the unique signaling and regulatory roles of calcium within the cell.

Calcium regulates many important aspects of cell activity, beginning with fertilization and ending with the apoptotic suicide of cells at the end of their life cycle [2][5]. Calcium ions are traditionally considered as secondary messengers liberated from intracellular and extracellular stores even though calcium itself may function to release Ca2+ ions from these stores. Calcium can also act as an extracellular primary messenger, thus taking on the role of a near-universal signaling molecule recognized by a wide variety of calcium-binding proteins in eukaryotes, prokaryotes and even viruses [2][6]. At the level of the protein structure, calcium may play important structural roles at the molecular level, required for maintaining appropriate conformations of individual proteins, such as at the β-propeller inter-domain interface in integrin α subunits and the homologous domains found in bacteria [7].

The scope of calcium ion functions extends to the regulation of contraction of all types of muscles, where binding of calcium to troponin C triggers the interaction of actin and myosin [8], [9]. As in most other cells, the concentrations of free Ca2+ ions in resting muscle is about 0.05–0.1 µM, but increases by one to two orders of magnitude when an external signal stimulates the cell. These Ca2+ ions are released either from intracellular stores (i.e. from the sarcoplasmic reticulum) or pumped up through sarcolemma from the intercellular space. In addition, increases in intracellular calcium concentrations play a central role in the function of neurons, triggering neurotransmitter release [10], [11]. Many aspects of neuronal activity, ranging from rapid modulation of channel function within a millisecond timescale to long-term switches in gene expression, are controlled by changes in the cytosolic calcium concentration. All of these various actions of Ca2+ ions are mediated by calcium-binding proteins that, in turn, interact with their target proteins. For example, regulation of transcription is coupled with numerous intracellular signaling processes often mediated by secondary messengers. Growing evidence points to the importance of calcium, as one of the most versatile second messengers, in activating or inhibiting gene transcription through actions frequently mediated by members of the EF-hand superfamily of calcium-binding proteins [12], [13].

Currently, there are eleven known calcium binding consensus motifs that proteins use to recognize Ca2+ ions, and whose sequence-variation profiles are established within the PROSITE database (http://prosite.expasy.org) [14]. The most commonly observed consensus motif for calcium binding is a characteristic DxDxDG calcium-binding loop [15], [16]. Among the 18 folds of proteins with the DxDxDG calcium-binding loop, there is a canonical EF-hand Ca2+-binding helix-loop-helix domain, where two antiparallel DxDxDG calcium-binding loops are flanked by interacting incoming and exiting α-helices [17][19]. The EF-hand domain is currently found in more than 71 non-redundant protein structures (see Homologous Superfamily EF-hand (1.10.238.10), the CATH database [20]) and is observed throughout all of the domains of life [17], [21][26]. We refer the reader to Gifford et al. (2007; Table 1 of the supplementary material in their article) for a listing of the range of functions performed by EF-hand containing proteins [17]. Given the wide range of functions of proteins with EF-hand calcium-binding domains, it is not surprising that these domains have been associated with numerous health issues, including Alzheimer's disease [27], cancer [28][30], neuronal disease [31], and disorders involving sodium ion channels in epilepsy, skeletal muscle disease, autism, and cardiac arrhythmias [32], [33].

thumbnail
Table 1. List of eleven non-redundant, representative calcium-bound X-ray and NMR protein complexes, which represent eleven different families of EF-hand domains.

https://doi.org/10.1371/journal.pone.0109287.t001

Conventionally, all EF-hand -containing proteins are divided into two groups: (1) Calcium sensors that include calmodulin, recoverin, and most of the other known EF-hand proteins, which function to translate the signal of a change in concentration of Ca2+ ions to various responses; and (2) Calcium buffers, represented by parvalbumin, calbindin D9k, calbindin D28k and calretinin that serve to modulate calcium signals as they bind free Ca2+ ions [17]. Calbindin-D28k and possibly calretinin, oncomodulin, and the mammalian β parvalbumin may have additional calcium sensor functions, leaving parvalbumin and calbindin-D9k as the only “pure” calcium buffers [34]. A typical hallmark for sensors is their relatively large calcium-dependent conformational changes, which are cooperative [35] and often accompanied by exposure of hydrophobic surfaces allowing the interaction with their target proteins [36][40].

During recent years, we have been examining the structural determinants that govern metal-controlled structural cooperativity, intrinsic disorder, dimer formation and various structural and chemical properties of parvalbumins and S100 proteins, which constitute the two largest sub-families of proteins with the Ca2+-binding EF-hand domain fold [41][44]. Here, we report two unique structural motifs found in all known families of the EF-hand proteins. These motifs provide a supporting scaffold for the calcium-binding DxDxDG sequence motif of the EF-hand. Each structural motif incorporates a cluster of three amino acids whose structure and structural rearrangement on calcium binding can serve to classify EF-hand domains into five separate classes.

Results and Discussion

Two Structural Motifs Outside of the Calcium Binding Regions Stabilize the EF-Hand Domain

Nearly 40 years ago, Kretsinger and Nockolds (1973) described what came to be known as the canonical “EF-hand” Ca2+-binding supersecondary structure [45], the well-known helix-loop-helix structural unit that many proteins use for binding of calcium ions [17][19]. Since then, many different subfamilies of EF-hand proteins have been recognized [46] and the Protein Data Bank (PDB) currently lists more than 280 crystal structures and 130 NMR structures that contain one or several domains with the “EF hand-like” fold belonging to the “EF-hand” structural superfamily [47]. Table 1 lists representative crystal or NMR structures from the eleven different structural families of the EF-hand containing proteins (SCOP; Fold: EF hand-like; Superfamily: EF-hand) [48], where each domain consists of two antiparallel helix-loop-helix substructures. By convention, the substructures are named as the “Odd” (the substructure appearing first along the sequence; N-terminal) and “Even” (the substructure appearing second along the sequence; C-terminal) motifs and the N-terminal and C-terminal flanking α-helices in each helix-loop-helix substructure are respectively called the “incoming” and “exiting” α-helix (Figure 1A) [19], [46].

thumbnail
Figure 1. The structural elements of the EF-hand fold.

Conserved elements of the EF-hand domains include (A) the flanking “incoming” and “exiting” α-helices (helix I to helix IV), the two interacting nodes, where the flanking α-helices interact and form interacting clusters I and II (shown in black and grey circles), position of the DxDxDG calcium binding loops, and the central β-sheet. (B) Structural alignment of the Odd and Even EF-hand motifs from the eleven fold representative structures. (C) Structural superposition of amino acids from Figure 1B. In 1K94_A, B, Figure 1B, “G3” designates two different three amino acid long insertions after Gly201 in 1K94_A and after Gly501 in 1K94_B. The structural frame of reference, “X”, “Y”, “Z”, “-X”, “-Y”, “-Z” and “X-4”, designate the seven key equivalent structural positions within all EF-hand domains, as numbered in Kretsinger and Nockolds [40]. The Odd EF-hand helix-loop-helix supersecondary structure is shown in grey. In panel A, B and C, the residues that belong to cluster I and cluster II are respectively highlighted in black and grey.

https://doi.org/10.1371/journal.pone.0109287.g001

On the basis of their sequences, databases such as Pfam generally lump the Odd and Even motifs together as a single entity [60]. However, structurally, the motifs are different. While they have equivalent characteristic DxDxDG sequences within the Ca2+ ion binding loops and the β-strands forming the antiparallel β-sheet linker of the EF-hand domain (Figure 1A), their differences lie within the flanking α-helices.

Here, we have treated the Odd and Even EF-hand motifs as two separate structural entities, first examining 11 representative structures individually in order to identify equivalent structural positions shared within the Odd motif and within the Even motif (Figure 1B). Then, equivalent positions were used to superpose the EF-hand domains from the eleven different structural families listed in Table 1 (Figure 1C). The comparisons showed that the flanking α-helices are not only different for the Odd and Even EF-hand motifs, but that they also interact differently with each other, forming two different interacting nodes, each with a cluster of three residues, located at opposite ends of the central β-sheet linker (Figure 1A). The distinctive structural differences within the flanking α-helices of the Odd motif and of the Even motif contribute to the interactions that play a critical role in the formation, conformational dynamics and function of the EF-hand domain.

Kretsinger and Nockolds (1973) labeled the six residues that bind calcium in the carp muscle calcium-binding protein as X, Y, Z and -X, -Y, -Z [45], and we have used these positions as a convenient frame of reference for labeling other amino acid locations within the EF-hand motif. Because there are two helix-loop-helix EF-hand motifs – the Odd and Even motifs – interacting with each other in every EF-hand domain, each domain contains two sets of X, Y, Z and -X, -Y, -Z positions (Figure 1A).

When we examined the structural superposition of the eleven representatives of the EF-hand containing proteins we found conserved characteristic clusters of interacting residues located around positions X-4. Because the Odd and Even EF-hand motifs are antiparallel, there are two such clusters located on the opposite sides of the antiparallel β-sheet linker, which include interacting residues from equivalent positions X-4, -X+1 and -Z+1 (Figure 1A). Cluster I (solid black circles, Figure 1A) locks the ends of the EF-hand domain, i.e. the incoming α-helix of the Odd EF-hand motif (residue X-4) and the exiting α-helix of the Even EF-hand motif (residues -X+1 and -Z+1), while cluster II (solid grey circles, Figure 1A) locks the exiting α-helix of the Odd EF-hand motif (residues -X+1 and -Z+1) and the incoming α-helix of the Even EF-hand motif (residue X-4).

Conservation of Amino Acids within the Two Clusters

The amino acids that constitute the two clusters are different. Cluster I is much less variable in sequence and mostly incorporates aromatic amino acids (phenylalanine, tryptophan, tyrosine): 15 out of 15 at X-4, 11/15 at -X+1 and 15/15 at -Z+1 (highlighted in black, Figure 1B). In contrast, cluster II includes a mix of aromatic, hydrophobic and polar amino acids of different sizes with 8/15 aromatic residues at X-4 (the rest are non-polar), 2/15 at -X+1 (others are variable in size and charge, e.g. glycine, lysine), and 4/15 at -Z+1 (the other residues are non-polar leucine and valine) (highlighted in grey, Figure 1B). In the one symmetric representative where one of three EF-hand domains is formed from homodimeric interactions (Ca2+-loaded human grancalcin; PDB code: 1K94), clusters I and II are obviously identical, fully aromatic, and no connecting loops are present (Figure 1B). The other two EF-hand domains in grancalcin are formed from a single-chain and are topologically asymmetric, similar to the other ten representative domains listed in Table 1.

Because of the differences in amino acid composition between the two clusters, we initially analyzed all interactions within them, with respect to whether they are considered stabilizing or destabilizing (Table S1 in File S1) according to the Automated Analysis of Interatomic Contacts software (LPC/CSU software) [61]. The LPC/CSU analysis reports that an equivalent number of non-polar stabilizing interactions could be seen in both clusters. However, cluster I, made up primarily of aromatic residues, almost completely lacked destabilizing interactions. In contrast, cluster II, composed of a wider mix of residue types, lacks some stabilizing interactions, made a larger number of destabilizing interactions, and even contained amino acids that did not directly interact with each other – i.e. an incomplete triad (Table S1 in File S1). For example, in the EF-hand of nucleoside diphosphate kinase (3BUX) tryptophan W202 (X-4) physically shields F205 from K225, preventing their interaction in cluster II. In osteonectin (1SRA), the interactions between H232 and F253 are absent as a result of a unique disulphide bridge, C256-C272, which joins the incoming and exiting α-helices of the even EF-hand helix-loop-helix motif. As a result of the disulphide bridge, the main-chain conformation and that of the F253 side chain are affected, resulting in a cavity filled with water (Wat429, Wat432 and Wat470 in 1SRA).

Clusters I and II are thus different even though they are formed in both instances by the symmetrically similar positioned incoming and exiting α-helices. Unlike cluster II, cluster I has a predominant aromatic mini-core that is clearly stabilized by a set of linked CH-π and CH-O hydrogen bonds present in ten of the eleven representative structures (Figure 2; Tables S2 and S3 in File S1). The observed differences between the two clusters may reflect both structural and functional requirements: Reduced variability in cluster I may reflect a need for additional stabilization at the domains' open end – where the chains enter and exit the EF-hand – whereas cluster II may be inherently more stabile because a loop links and restrains the two ends of the EF-hand thus allowing for a more variable core. The difference also likely reflects functional and conformational requirements related to dynamic cooperative binding and release of Ca2+ ions and protein ligands.

thumbnail
Figure 2. Example interactions within cluster I.

(A) and (B) illustrate two types of interactions between the flanking α-helices I and IV (cluster I in Figure 1). The interactions occur via amino acids at positions -X+1, X-4 and -Z+1, which are also shown as black circles in Figures 1A and 1C and highlighted in black in the alignment in Figure 1B. Because in cluster I, positions X-4 and -Z+1 are purely aromatic (Φ) in all EF-hand representative structures, cluster I is called aromatic. (C) Contains the description of interactions for the eleven EF-hand fold representatives. The pattern “(-X+1)helixIV CH-π Φ(X-4)helixI CH-π Φ(-Z+1)helixIV weak HB (-X+1)helixIV” indicates a circular interaction, where a side-chain atom of the –X+1 residue from the flanking α-helix IV forms a CH-π interaction with the ring of the X-4 aromatic amino acid from helix I, which, in turn, forms a CH-π interaction with the ring of the –Z+1 aromatic amino acid from helix IV, which interacts with the initial –X+1 residue from the flanking α-helix IV by means of a weak CH-O hydrogen bond.

https://doi.org/10.1371/journal.pone.0109287.g002

In order to investigate the role of the two conserved but different clusters of amino acids stabilizing the structure of the EF-hand domain and the effect of calcium binding on the domain conformation, we have applied two different approaches: (1) For the eleven family representative structures, we calculated the total interacting surface areas among the amino acids within each cluster, as well as the interacting surface area between the two clusters within the domain (Table 2); (2) We also calculated the (i) interacting surface areas (Table 3; Table S4 in File S1) and (ii) the root mean-squared deviations of superpositioned clusters (RMSD; Table 4; Table S5 in File S1) for comparisons of all EF-hand domains found in the PDB (at the time when analyzed) whose structures are known for both the calcium-bound (holo; some also with bound target protein) and calcium-unbound (apo) forms. Individual values for all compared proteins are shown in Tables S4 and S5 in File S1.

thumbnail
Table 2. Summarized areas of interacting surfaces among the amino acids within clusters I and II (black and grey clusters in Figure 1, respectively), and the interacting area between the two clusters (I/II), within the EF-domains of the eleven representative structures.

https://doi.org/10.1371/journal.pone.0109287.t002

thumbnail
Table 3. Effects of calcium binding on the conformation of all EF-hand domains whose structures are known in the apo-form and with bound Ca2+ ions, and target protein ligands.

https://doi.org/10.1371/journal.pone.0109287.t003

thumbnail
Table 4. Effects of calcium binding on the conformation of all EF-hand domains, whose structures are known in the apo-form and with bound Ca2+ ions, and target protein ligands (where they exist).

https://doi.org/10.1371/journal.pone.0109287.t004

EF-Hand Family-Dependent Variation of the Two Conserved Clusters

The interacting surface area data from the fold representative structures (Table 2) show that the two clusters are distinctly different. Cluster I (with the exception of one domain in Cbp40 – plasmodial specific calcium-binding protein LAV1–2; 1IJ5) appears to be more compact, consistently having a larger interacting surface area than cluster II. In the representative fold structures, the interacting surface areas of cluster I are similar, regardless of whether the structure contains two (2PVB_A), one (1F8H_A) or no (1WLM_A) bound calcium atoms. In contrast, the interacting surface areas seen for cluster II vary significantly, from values closer to those of cluster I, reflecting a compact cluster, to values as low as 34.1 Å2.

The interacting area between clusters I and II varies from no contact in calmodulin to 34.2 Å2 in grancalcin and 45.8 Å2 in the EF-hand multi-domain protein Cbl. Where interactions occur between clusters I and II, the amino acid at positions -Z+1 from both clusters are involved and the exiting α-helices of the two EF-hand motifs are consequently positioned closer to each other.

Effects of Calcium Binding on the Conformation of Individual Clusters Across Structures With and Without Bound Calcium

When considering all (19 in total) available pairs of EF-hand domains whose three-dimensional structures are known with and without bound calcium, the interacting surface area (Table 3) and RMSD (Table 4) values consistently show that interactions within cluster I are not altered dramatically on the binding of calcium or calcium plus a target protein (contact area with and without bound calcium averaging 89±12 Å2 and 84±11 Å2, respectively; RMSD for superpositioned Cα atoms of 0.6±0.2 Å; rows A in Tables 3, 4). Indeed, for both cluster I and cluster II the presence or absence of a protein ligand had little or no impact on the interactions within the cluster. In comparison to cluster I, the maximum contact surface area within cluster II averages about 20 Å2 less and in some cases the interacting surface area increases considerably on calcium binding. Cluster II can be divided into two groups. In one group, calcium binding produces little change in both the contact surface area within the cluster (68±12 Å2 → 66±15 Å2; row B in Table 3) and the local conformation (RMSD, 0.6±0.3 Å; row B in Table 4). In the other group, calcium binding leads to dramatically increased interactions among the three residues of cluster II (11±8 Å2 → 66±10 Å2; row C, Table 3) along with a notable conformational difference (RMSD, 1.6±0.1 Å; row C in Table 4).

Effects of Calcium Binding on the Overall Fold of the EF-Hand Domain

The effects of binding of Ca2+ ions on the overall fold can be seen through the differences in the interacting surface areas between clusters I and II (column “B–G” in Table S4 in File S1), and RMSD values for the superpositioning of the two clusters as a unit (column “B–G” in Table S5 in File S1). Both sets of comparisons indicate that all EF-hand domains are divided into either those that do not change the relative position between clusters I and II: 2±2 Å2 (apo) → 2±2 Å2 (Ca2+) (row E in Table 3), and 32±2 Å2 (apo) → 28±7 Å2 (Ca2+) (row F in Table 3), and with an RMSD of 0.8±0.3 Å (row E in Table 4); or those that change their conformation upon calcium binding, 33±14 Å2 (apo) → 3±3 Å2 (Ca2+) (row D in Table 3), and the RMSD of 1.7±0.2 Å (row D in Table 4).

Effects of Calcium Binding on the Conformation of Individual Clusters and the Overall Fold of the EF-Hand Domains, which bind single Ca2+ ion

All 19 pairs of EF-hand domains, whose conformational changes upon calcium binding are analyzed and shown above in Tables 3 and 4, require two Ca2+ ions in two defined calcium binding sites. However, there is also a small group of domains with the single calcium binding EF-hand motifs. Currently, there are only four such proteins, whose three-dimensional structures are known with and without bound calcium. These four proteins incorporate six different EF-hand domains shown in Table 5.

thumbnail
Table 5. Single calcium binding EF-hand domains with known apo- and holo-form structures.

https://doi.org/10.1371/journal.pone.0109287.t005

Similarly to the EF-hand domains with two Ca2+ ions, the single calcium binding EF-hand motifs show that interactions within cluster I are not altered dramatically on the binding of calcium or calcium plus a target protein (row A in Table 3; and Table 5). Moreover, with the exception of recoverin (C-domain), also the conformational changes within the cluster II and the relative conformation between clusters I and II do fall within the categories of double calcium binding EF-hand motifs (Table 3). For example, in all known EF-hand domains of grancalcin and the myosin essential light chain, the conformation of cluster II and the relative conformation between clusters I and II do not significantly change upon calcium binding, as seen, for example, in the calbindin D9K calcium buffer (rows B and E in Table 3; and Table 5). On the other hand, recoverin (N-domain) and troponin C (cardiac) do rearrange the relative conformation between clusters I and II and undergo domain opening, similar to the N- and C-domains of calmodulin (rows B and D in Table 3; and Table 5).

A notable exception from the behavior of double calcium binding EF-hand motifs is shown by the recoverin, C-domain. Calcium binding to this protein results in significant weakening of interactions within the cluster II as opposed to all the other known single or double calcium binding EF-hand motifs descried above, where cluster II becomes tighter and starts resembling cluster I. The domain conformation of recoverin (C-domain) is already open without bound Ca2+ ion, and remains open after the Ca2+ ion is bound.

Clusters I and II and a Short β-Sheet Linker Form the Hydrophobic Core of the EF-Hand Domain

Clusters I and II join the incoming and exiting helices of the EF-hand domain structure and predominantly consist of non-polar and hydrophobic amino acids. The short β-sheet linker of the EF-hand domain is also mainly hydrophobic and is formed from two β-strands bound together by two hydrogen bonds (shown as dashed lines, “H-bonds” in Figure 1A). Both clusters interact with the β-sheet linker and together contribute to the hydrophobic core of the EF-hand. The side chain of the central amino acid (X-4) of cluster I (e.g. F47 in 2PVB) in all eleven representative structures of the EF-hand domain is always positioned in close proximity and perpendicular to the plane of the β-sheet linker (I58 and I97 in 2PVB) (Figure 3). The side-chain ring of residue X-4 forms two CH-O close contacts with the β-strand and the contact distance and angle parameters of the contacts satisfy criteria for classic weak CH-O hydrogen bonds (Table S6 in File S1) [63].

thumbnail
Figure 3. Interaction of cluster I with the central β sheet.

The side chain of the aromatic residue at position X-4 in cluster I is perpendicular and directly interacts with the β-sheet linker of the EF-hand domain, contributing to its hydrophobic core. (A) Interactions between residue X-4 of cluster I (F47) and the β-sheet linker in 2PVB. (B) The interactions in the fold-representative structures include two classic, weak CH-O hydrogen bonds, whose detailed parameters are given in Table S6 in File S1. The (X-4)odd – O/(-Y-1)odd interaction shows the weak CH-O hydrogen bond between residue X-4 of cluster I and residue -Y-1 (G56 in 2PVB) of the “odd” EF-hand motif, while the (X-4)odd – O/(-Y+1)even interaction shows the weak CH-O hydrogen bond between the same residue X-4 of cluster I and residue -Y-1 (I97 in 2PVB) of the “even” EF-hand motif.

https://doi.org/10.1371/journal.pone.0109287.g003

Because the amino acid at X-4 in cluster II is not strictly aromatic, but sometimes hydrophobic or even polar (Figure 1B), several types of interactions take place with the β sheet. If the central amino acid X-4 is an aromatic residue, interactions with the short β-sheet linker are similar to those for cluster I, and both clusters contribute to the EF-hand hydrophobic core. Where residue X-4 is not aromatic, the interactions with the β-sheet linker vary from being polar (e.g. multi-domain protein Cbl, 3BUX_B, where the side-chain NZ atom of K225 interacts with the side-chain OD1 atom of D234 through the HOH445 water molecule) to being weak Van der Waals interactions (e.g. parvalbumin, 2PVB, where the L86 side chain interacts with the Cα atom of G95 at 4.5 Å). In the representative structures the side chain of residue X-4 of cluster II points directly towards and interacts with the β-sheet linker, contributing to hydrophobic core's stability (Figure 3).

Biological Implications

EF-Hand Domain Clusters I and II are Structurally and Functionally Non-Equivalent.

Regardless of the outward symmetry in the folding pattern of the EF-hand domains (Figure 1A), the two halves of the domain have been shown experimentally to be non-equivalent and non-interchangeable, and thus asymmetrical. Lakowski et al. (2007) engineered and refolded calmodulin creating a “reversed” calmodulin fragment starting with the even helix-loop-helix EF-hand motif from the N-terminal EF-hand domain, including the inter-domain region, and ending with the odd helix-loop-helix EF-hand motif from the C-terminal EF-hand domain [72]. As a result, this EF-hand construct would have the fold of the EF-hand domain, as in calmodulin, with similar flanking α-helices and the short β-strand linker, but the locations of clusters I and II would be swapped. These structural changes led to ∼100 times weaker binding of the first calcium atom and 3000 times weaker binding of the second calcium atom, rendering the calmodulin fragment inactive [72]. It may be a general requirement of EF-hands that cluster I must have a more static structure with a high level of mutual interactions because it is located where the polypeptide chains enter and leave the domain. In contrast, cluster II has a constraining 5–15 residue loop [19] joining the two helices where residues of cluster II are located and may account for the higher degree of amino acid variability in comparison with cluster I. Variability at cluster II may also play an important role and be necessary in order to allow for different degrees of conformational change on calcium binding and release that are directly associated with the molecular and biological function.

Clusters I and II, and Nearby Residues are Critical for Both Structure and Function.

The functional importance of the region adjacent to the calcium-binding motif can also be seen through the effects of naturally-occurring and engineered mutations. There are a number of EF-hand containing proteins that are linked with health issues [27], but only a few mutations have been clearly established as being directly linked to a particular condition. For example, in guanylate cyclase activating protein 1 (GCAP1) [73] polymorphisms have been directly associated with autosomal dominant cone dystrophy: i.e. Y99C at position X-1 of the odd EF-hand motif [74] and the replacement of asparagine-threonine for I143 (X-1) of the even EF-hand motif [75].

Engineered mutations have also pointed to the functional importance of clusters I and II and nearby residues. Koltzscher and Gerke (2000) have shown that for the calcium buffer S100P an F15A mutation at X-4 in cluster I completely abolished dimerization required for activation [76]. In Calbindin D9K (fold family 1 in Table 1; 1IG5_A in Figure 1B), the mutation of amino acids to alanine at X-4 (F10) and -Z+1 (F66) in cluster I and -Z+1 (L28) in cluster II significantly reduced functional calcium binding whereas the mutation of phenylalanine at -Z+1 to the larger aromatic tryptophan increased calcium binding 25-fold (Kragelund et al., 1998) [77]. Tikunova et al. (2002) examined the role of 27 hydrophobic residues within the EF-hand domain by replacing these residues with polar glutamine in the F29W activating mutant of the N-domain of troponin C (N-TnC, calmodulin-like family, fold family 6 in Table 1) [78]. The largest observed effects were the 123-fold decrease in the Kd for F26QTnCF29W at X-4 in cluster I, followed by I37QTnCF29W within the β-strand linker (↓24%), I62QTnCF29W at X-4 in cluster II (↓12%), and F78QTnCF29W at -Z+1 in cluster I (↓8%), and similar 6 to 7 -fold changes at I73, V80 and M81, but small effects on -X+1 (F75) in cluster I and -Z+1 (L42) in cluster II. Increases in the Kd were seen for V45QTnCF29W (↑19%) and L49QTnCF29W (↑19%) also located the exiting helix of the odd EF hand motif. The packing of the hydrophobic triad in the vicinity of cluster I of the C-terminal calmodulin-like domain (fold family 6, Table 1; PDB code 4DCK) was shown to be critical for the normal functioning of the sodium channel of NaV1.5 [33]. The degree of likely disruption to the local hydrophobic environment was associated with the severity of effects on channel gating by mutations of Y1795 at X-7 (incoming helix of the odd motif), for which inherited mutations cause congenital cardiac arrhythmias, and X-4 (W1798) and -Z+1 (I1853) within cluster I. In the ALG-2 protein (PDB code: 2ZN9; penta-EF-hand protein family; family fold 9, Table 1), the F148S mutation at cluster I alters the wild-type binding site recognition for the protein ligand phospholipid scramblase protein PLSCR3 [79].

Clusters I and II Can Serve as Markers for EF-Hand Domain Types.

EF-hand domains have been classified as buffers and as sensors [36], [80], [81]. Sensors such as calmodulin and troponin C [34] as well as S100A1, S100A4, S100A6 and S100B [82] would involve communication with a protein ligand and involve a calcium-sensitive conformational change in order to expose hydrophobic residues and thus mediate the interaction. In contrast, buffers (signal modulators) [81] would not, and a lack of a conformational change within the EF-hand domain would indicate that a domain was a buffer. For example, S100A16, whose structure did not show any significant conformational change upon calcium binding [83]. Similarly, no significant structural changes on calcium binding were observed for the calpain C-terminal domain [84], Cbp40 [85] and S100P [86] and they would be considered as likely buffers on this basis.

The structural changes refer to the transition between “closed” and “open” protein conformation, which may or may not occur on calcium binding, and the resulting calcium-dependent/independent or calcium-sensitive/insensitive mechanism of protein function. Originally, it was shown that in calmodulin and troponin C, distances and angles between the flanking α-helices did change upon calcium binding, while in calbindin 9K they did not [19]. This calcium dependent or independent behavior of EF-hand domains was later extended to other EF-hand proteins, where the domains were examined to see whether they would undergo structural rearrangements similar to those seen in calmodulin and troponin C (calcium-sensitive or sensor mechanism) or would show no structural changes as seen in calbindin 9K (calcium-insensitive or buffer mechanism) [19], [82]. However, as Schwaller (2009) notes, “pure” calcium buffers like calbindin D9K and parvalbumin are few, as the buffers calbindin-28K, oncomodulin and calrectinin may also have calcium sensor functions [34].

The differences observed among pairs of structures with and without bound calcium have allowed us to effectively separate and classify EF-hand domains into five separate groups (Figure 4). The quasi-symmetry of the EF-hand domain belies the fact that the component parts, with the rare exception of homodimeric interactions, are not identical and not truly symmetrical and that the domain seen across different proteins can have a largely static structure or one that rearranges considerably on binding Ca2+ ions. As a key underlying component of the EF-hand domain structure, the three-residue clusters, I and II, are critical for the structure of the EF-hand domains they are found within and consequently the overall function of the protein. Cluster I is primarily aromatic and as a unit it appears to be the more static and stabile linchpin for the overall structure since it maintains higher mutual surface area contacts that do not rearrange dramatically on calcium binding in any of the available structures, varying between 80±11 Å2 and 89±12 Å2 (row A, Table 3). In contrast, cluster II is more variable in sequence, the interaction surface area is lower than for cluster I, ranging between 11±8 Å2 (row C, Table 3) and 68±12 Å2 (row B, Table 3), and the contacts can reorganize significantly when Ca2+ ions are bound or released. Thus, it would appear that the makeup of cluster II is important for the functional requirements of the individual EF-hand.

thumbnail
Figure 4. Five distinguishable groups of EF-hand domains based on the degree of structural rearrangements within clusters I and II, and the inter-cluster interactions that take place upon calcium binding: (1) Open Static (open domain conformation and no conformational changes); (2) Closed Static (closed domain conformation and no conformational changes); (3) Local Dynamic (simultaneous conformational changes in clusters I and II and the entire domain); (4) Dynamic (only global conformational changes, but not in clusters I and II); and (5) Local Static (stable open domain conformation, conformational changes only in clusters I and II, but not the entire domain).

Domain level conformational changes, are shown by the small (“closed” domains) and large (“open” domains) distance between clusters I and II (triangles, ellipses and inverted triangles). Domain types (3) and (4) do undergo domain opening, while domain types (1), (2) and (5) do not. The conformation of the domains of type (2) remains “closed”, while the conformation the domains of types (1) and (5) remains open. Local conformational changes in clusters I and II are shown by normal triangles (compact conformation of cluster I), inverted triangles (compact conformation of cluster II), and ellipses (less compact, more open conformation of cluster II). The conformation of cluster I does not change in all of the known structures and is the same in buffers and sensors, such as in calbindin 9K, calmodulin and troponin C. The conformation of cluster II does change from being less compact to more compact in the domains of types (3) and (5). In group (4), EF-hand domains undergo domain opening, but the conformations of the conserved clusters I and II remain intact.

https://doi.org/10.1371/journal.pone.0109287.g004

The interaction surface area for cluster II would correlate with the function of the EF-hand according to the buffer-sensor classification as follows. If the interaction area within cluster II, apo form, is small (row C, Table 3), then this points towards the domain being a calcium sensor. If the interaction surface area between clusters I and II for the holo form is small (rows E and F, Table 3), then it is likely that the domain is a calcium buffer.

Individual EF-Hand Domains can be Classified into Five Groups Based on the Characteristics of Clusters I and II.

Altogether, the EF-hand domains fall into five distinguishable groups based on the degree of structural rearrangements that take place upon calcium binding (Figure 4). In all five types the three residues of cluster I form mutually compact interactions.

Type 1. Open static EF-hand domains – where the domain conformation is open, cluster II has compact interactions, and clusters I and II do not interact. These features are not affected by the presence or absence of calcium. Observed for the structures of EF-hand domains considered to function as buffers: calbindin 9K, S100P, parvalbumin, oncomodulin and Cbp40 (N-terminal domain).

Type 2. Closed static EF-hand domains – where the domain conformation is closed, cluster II is compact, and clusters I and II are close to each other and interact with each other. These features are not affected by the presence or absence of calcium. These features are shared by the calcium-free structures of the Type 4 set (see below). Includes the buffers S100A16, Cbp40 (C-terminal domain) and calpain (N-terminal domain).

Type 3. Local dynamic EF-hand domains – where a closed-to-open transition takes place on calcium binding in which the cluster II becomes more compact and interactions between clusters I and II are eliminated. Includes the calcium sensors, such as S100A1, S100A4-A6 and S100B.

Type 4. Dynamic EF-hand domains – where cluster II is compact in the absence of calcium but where the inter-cluster distance increases on calcium binding. Includes the classical calcium sensors such as calmodulin (N- and C- domains) and troponin C, and polcalcin.

Type 5. Local static EF-hand domains – where the conformation of cluster II becomes more compact on calcium binding but the inter-cluster distance remains essentially unaltered. Includes the calcium sensors S100A13 and calpain (C-domain).

With the exception of the Type 2 domains, all other domains on calcium binding have clusters I and II that are compact and are separated from each other (Figure 4). The clearly different functional role of clusters I and II in calcium buffers and sensors explains why the engineered swap of the clusters in the EF-hand domain construct of calmodulin led to the loss of function described by Lakowski et al. [72]. As a consequence of these comparisons of EF-hand domain structures, it should be possible to classify some EF-hands based on a single example of a three-dimensional structure. For structures with bound calcium, closed static EF-hand domains of Type 2 can be differentiated from the other four types; whereas, in the unbound state it would be possible to identify domains belonging to Types 1, 3 and 5. The overall conformation of EF-hand domains are highly variable [37], [81]; but, when considered in terms of the two local, three-residue clusters that help maintain the domain structure, it is possible to classify individual proteins into discrete groupings. Detailed summary on the datasets after classifying them into five groups, including characteristic features of the classes, as well as PDB and sequence IDs (http://www.uniprot.org) is shown in Table S7 in File S1.

Materials and Methods

Currently, there are eleven known calcium binding consensus motifs [14], whose sequence-variation profiles are established within the PROSITE database [87]. The EF-hand motif is currently found in more than 71 non-redundant representative domains (see Homologous Superfamily EF-hand (1.10.238.10), the CATH database [20]). We have extracted and analyzed more than 280 (all) crystal structures and 130 NMR structures currently listed in the Protein Data Bank (PDB) [47] that contain one or several domains with the “EF hand-like” fold belonging to the “EF-hand” structural superfamily, according to the SCOP database (SCOP; Fold: EF hand-like; Superfamily: EF-hand) [48].

From 410 available X-ray and NMR structures of EF-hand domains, we had to create three different datasets: (1) a representative non-redundant dataset of different EF-hand domains to scan for structural motifs, i.e. clusters of amino acids and interactions, which are conserved across different structural families of EF-hand proteins; (2) a dataset to study effects of calcium binding on the domain conformation of double calcium binding EF-hand motifs, which would include all known pairs of EF-hand domains, whose three-dimensional structures are known with and without two bound calcium atoms; and (3) a similar dataset to study calcium binding effects on single calcium binding EF-hand motifs. For the first task of searching for conserved amino acid clusters and interactions, we had to manually choose a set of best resolution representative structures from the 11 different structural families of EF-hand proteins, given in SCOP [48]. Where a family contained only NMR structures, a most represented and complete NMR structure was chosen (fold families 7 and 11, Table 1). The search for structural motifs and creation of the first dataset involved a semi-manual pairwise and multiple global and local structural comparisons of the 11 representative structures using the Accelrys Discovery Studio molecular modeling environment (www.accelrys.com). Because of semi-manual all-against-all structural comparisons and analysis of local structural similarities, we stayed at the structural family level (11 structures), and did not go below, to the level of individual protein domains (51 structures). Beyond the 51 structures of individual protein domains, the rest of the 410 structures included either orthologs of the same proteins from the different species or structures of various mutants, which could be discarded from analysis. For the second and third tasks to study effects of calcium binding on the domain conformation, the creation of two different datasets was straightforward. We manually analyzed all 410 structures of EF-hand domains and the corresponding literature to select all double and single calcium binding EF-hand domains, whose three-dimensional structures are known with and without bound calcium and possibly other ligands.

All types of structural superpositions and RMSD calculations based on back-bone and all atoms other than hydrogen, including (1) the superposition of the entire Odd and Even EF-hand motifs of EF-hand containing proteins; and (2) local structural superpositions of clusters I and II, were done using the SuperPose superposition server [88] and the Accelrys Discovery Studio molecular modeling environment (www.accelrys.com).

Calculation of atomic contacts and interacting surface areas was done using the Contacts of Structural Units (CSU) software, which is based on the surface complementarity approach developed by Sobolev et al. [61]. The “Detailed Analysis” procedure within the CSU software was used to calculate bond and surface parameters for all amino acids of the clusters individually.

Geometric parameters to assign CH-π interactions were chosen to satisfy criteria given in Brandl et al. [62]. If X designates the center of an aromatic ring, then the C-X distance must be ≤4.5 Å; the C-H-X angle must be greater than 120°; and the dHp-X projection distance must be ≤1.2 Å [62]. We used an additional criteria, d(H-X) ≤3.5 Å, to ensure that the CH group points directly to the center of the π-ring (Table S2 in File S1).

The criteria to assign weak hydrogen bonds were taken from Derewenda et al. [63]. Firstly, the C-H-O angle ζ must be greater than 90°. Secondly, an electronegative atom must be located adjacent to the carbon atom, such that the acidity of hydrogen atoms attached to the carbon atom increases, and consequently, the carbon atom could be a hydrogen bond donor. Thirdly, the C-O distance must be ≤4.0 Å and the H-O distance must be ≤3.0 Å. Two distance criteria d(N-O) ≤3.7 Å and d(H-O) ≤2.7 Å were used for the conventional hydrogen NH-O bonds; angular criteria, as described above, were also imposed.

All geometrical calculations (i.e., angles, torsion angles and distances) were made using the Accelrys Discovery Studio molecular modeling environment. Color figures in this manuscript were produced with MOLSCRIPT [89] and Raster3D [90].

Supporting Information

File S1.

Contains Tables S1–S7. Table S1. Assignment of interatomic contacts among the residues in clusters I and II according to the Interatomic Contacts software [61]. The “+” and “−” signs show presence or absence of Stabilizing and Destabilizing contacts, respectively, as defined by the Interatomic Contacts program. NI, absence of any interactions between the amino acids. 1EXR and 1IJ5 contain two EF-hand domains within chain A, while 1K94 contains two EF-hand domains within chain A, and one additional domain, shared between chains A and B. Table S2. Values of distances and angles for the CH-π interactions shown in Figure 2 for the interactions within cluster I. All designations are as in [62]. Table S3. Values of distances and angles for the CH-O interactions shown in Figure 2 for the interactions within cluster I. All designations are as in [63]. Table S4. Summarized areas of interacting surfaces among the amino acids within cluster I and cluster II (Figure 1), and the interacting area between the two clusters (I/II), for all EF-hand domains with known apo- and holo-form structures. NI, no interaction between the clusters. Table S5. Structural similarity as RMSD values calculated for the same protein with and without bound calcium; for clusters I and II separately and cluster I and cluster II together (I/II). Table S6. Values of distances and angles for the interactions shown in Figure 3 for the interaction of cluster I with the central β sheet. All designations are as in [63]. Table S7. Five classes of EF-Hand domains, based on the characteristics and conformational changes within the Clusters I and II (local level conformational changes) and within the entire EF-hand domain (global, domain level, conformational changes).

https://doi.org/10.1371/journal.pone.0109287.s001

(DOC)

Acknowledgments

We thank Robert H. Kretsinger, University of Virginia, School of Medicine, for support, helpful discussions, and careful reading of the manuscript.

Author Contributions

Conceived and designed the experiments: KD SP AD EP MSJ. Performed the experiments: KD SP AD. Analyzed the data: KD SP AD EP MSJ. Contributed reagents/materials/analysis tools: EP MSJ. Wrote the paper: KD EP MSJ.

References

  1. 1. Krebs J, Michalak M (2007) Calcium: A matter of life or death. Paris: Elsevier.
  2. 2. Permyakov EA (2009) Metalloproteomics. A John Wiley & Sons, Inc., Hoboken, New Jersey.
  3. 3. Permyakov EA, Kretsinger RH (2011) Calcium Binding Proteins. A John Wiley & Sons, Inc., Hoboken, New Jersey.
  4. 4. Dominguez DC (2004) Calcium signalling in bacteria. Mol Microbiol 54: 291–297.
  5. 5. Clapham DE (2007) Calcium signaling. Cell 131: 1047–1058.
  6. 6. Zhou Y, Frey TK, Yang JJ (2009) Viral calciomics: interplays between Ca2+ and virus. Cell Calcium 46: 1–17.
  7. 7. Chouhan B, Denesyuk A, Heino J, Johnson MS, Denessiouk K (2011) Conservation of the human integrin-type beta-propeller domain in bacteria. PLoS One 6: e25069.
  8. 8. Marchand A, Abi-Gerges A, Saliba Y, Merlet E, Lompré AM (2012) Calcium signaling in vascular smooth muscle cells: from physiology to pathology. Adv Exp Med Biol 740: 795–810.
  9. 9. Baylor SM, Hollingworth S (2012) Intracellular calcium movements during excitation-contraction coupling in mammalian slow-twitch and fast-twitch muscle fibers. J Gen Physiol 139: 261–272.
  10. 10. Mikhaylova M, Hradsky J, Kreutz MR (2011) Between promiscuity and specificity: novel roles of EF-hand calcium sensors in neuronal Ca2+ signalling. J Neurochem 118: 695–713.
  11. 11. Woods NK, Padmanabhan J (2012) Neuronal calcium signaling and Alzheimer's disease. Adv Exp Med Biol 740: 1193–1217.
  12. 12. Machaca K (2011) Ca2+ signaling, genes and the cell cycle. Cell Calcium 49: 323–330.
  13. 13. Naranjo JR, Mellström B (2012) Ca2+-dependent transcriptional control of Ca2+ homeostasis. J Biol Chem 287: 31674–31680.
  14. 14. Santamaría-Hernando S, Krell T, Ramos-González MI (2012) Identification of a novel calcium binding motif based on the detection of sequence insertions in the animal peroxidase domain of bacterial proteins. PLoS One 7: e40698.
  15. 15. Rigden DJ, Galperin MY (2004) The DxDxDG motif for calcium binding: multiple structural contexts and implications for evolution. J Mol Biol 343: 971–984.
  16. 16. Rigden DJ, Woodhead DD, Wong PW, Galperin MY (2011) New structural and functional contexts of the Dx[DN]xDG linear motif: insights into evolution of calcium-binding proteins. PLoS One 6: e21507.
  17. 17. Gifford JL, Walsh MP, Vogel HJ (2007) Structures and metal-ion-binding properties of the Ca2+-binding helix-loop-helix EF-hand motifs. Biochem J 405: 199–221.
  18. 18. Moncrief ND, Kretsinger RH, Goodman M (1990) Evolution of EF-hand calcium-modulated proteins. I. Relationships based on amino acid sequences. J Mol Evol 30: 522–562.
  19. 19. Grabarek Z (2006) Structural basis for diversity of the EF-hand calcium-binding proteins. J Mol Biol 359: 509–525.
  20. 20. Sillitoe I, Cuff AL, Dessailly BH, Dawson NL, Furnham N, et al. (2013) New functional families (FunFams) in CATH to improve the mapping of conserved functional sites to 3D structures. Nucleic Acids Res 41(Database issue): D490–D498.
  21. 21. Strynadka NC, James MN (1989) Crystal structures of the helix-loop-helix calcium-binding proteins. Annu Rev Biochem 58: 951–998.
  22. 22. Nakayama N, Kawasaki H, Kretsinger RH (2000) Evolution of EF-hand proteins. Topics Biol Inorg Chem 3: 29–58.
  23. 23. Lewit-Bentley A, Réty S (2000) EF-hand calcium-binding proteins. Curr Opin Struct Biol 10: 637–643.
  24. 24. Zhou Y, Yang W, Kirberger M, Lee HW, Ayalasomayajula G, et al. (2006) Prediction of EF-hand calcium-binding proteins and analysis of bacterial EF-hand proteins. Proteins 65: 643–655.
  25. 25. Yanyi C, Shenghui X, Yubin Z, Jie YJ (2010) Calciomics: prediction and analysis of EF-hand calcium binding proteins by protein engineering. Sci China Chem 53: 52–60.
  26. 26. Mazumder M, Padhan N, Bhattacharya A, Gourinath S (2014) Prediction and analysis of canonical EF hand loop and qualitative estimation of Ca2+ binding affinity. PLoS One 9: e96202.
  27. 27. Braunewell KH (2005) The darker side of Ca2+ signaling by neuronal Ca2+-sensor proteins: from Alzheimer's disease to cancer. Trends Pharmacol Sci 26: 345–351.
  28. 28. Emberley ED, Murphy LC, Watson PH (2004) S100 proteins and their influence on pro-survival pathways in cancer. Biochem Cell Biol 82: 508–515.
  29. 29. Subramanian L, Polans AS (2004) Cancer-related diseases of the eye: the role of calcium and calcium-binding proteins. Biochem Biophys Res Commun 322: 1153–1165.
  30. 30. Wolf S, Haase-Kohn C, Pietzsch J (2011) S100A2 in cancerogenesis: a friend or a foe? Amino Acids 41: 849–861.
  31. 31. Seaton G, Hogg EL, Jo J, Whitcomb DJ, Cho K (2011) Sensing change: the emerging role of calcium sensors in neuronal disease. Semin Cell Dev Biol 22: 530–535.
  32. 32. Hübner CA, Jentsch TJ (2002) Ion channel diseases. Hum Mol Genet 11: 2435–2445.
  33. 33. Glaaser IW, Bankston JR, Liu H, Tateyama M, Kass RS (2006) A carboxyl-terminal hydrophobic interface is critical to sodium channel function. Relevance to inherited disorders. J Biol Chem 281: 24015–24023.
  34. 34. Schwaller B (2009) The continuing disappearance of “pure” Ca2+ buffers. Cell Mol Life Sci 66: 275–300.
  35. 35. Nelson MR, Thulin E, Fagan PA, Forsén S, Chazin WJ (2002) The EF-hand domain: a globally cooperative structural unit. Protein Sci 2: 198–205.
  36. 36. Ikura M (1996) Calcium binding and conformational response in EF-hand proteins. Trends Biochem Sci 1: 14–17.
  37. 37. Yap KL, Ames JB, Swindells MB, Ikura M (1999) Diversity of conformational states and changes within the EF-hand protein superfamily. Proteins 37: 499–507.
  38. 38. Ikura M, Ames JB (2006) Genetic polymorphism and protein conformational plasticity in the calmodulin superfamily: two ways to promote multifunctionality. Proc Natl Acad Sci USA 103: 1159–1164.
  39. 39. Stathopulos PB, Zheng L, Li GY, Plevin MJ, Ikura M (2008) Structural and mechanistic insights into STIM1-mediated initiation of store-operated calcium entry. Cell 135: 110–122.
  40. 40. Kawasaki H, Kretsinger RH (2012) Analysis of the movements of helices in EF-hands. Proteins 80: 2592–2600.
  41. 41. Permyakov SE, Bakunts AG, Denesyuk AI, Knyazeva EL, Uversky VN, et al. (2008) Apo-parvalbumin as an intrinsically disordered protein. Proteins 72: 822–836.
  42. 42. Permyakov SE, Bakunts AG, Permyakova ME, Denesyuk AI, Uversky VN, et al. (2009) Metal-controlled interdomain cooperativity in parvalbumins. Cell Calcium 46: 163–175.
  43. 43. Permyakov SE, Ismailov RG, Xue B, Denesyuk AI, Uversky VN, et al. (2011) Intrinsic disorder in S100 proteins. Mol Biosyst 7: 2164–2180.
  44. 44. Permyakov SE, Vologzhannikova AA, Emelyanenko VI, Knyazeva EL, Kazakov AS, et al. (2012) The impact of alpha-N-acetylation on structural and functional status of parvalbumin. Cell Calcium 52: 366–376.
  45. 45. Kretsinger RH, Nockolds CE (1973) Carp muscle calcium-binding protein. II. Structure determination and general description. J Biol Chem 248: 3313–3326.
  46. 46. Kawasaki H, Nakayama S, Kretsinger RH (1998) Classification and evolution of EF-hand proteins. Biometals 11: 277–295.
  47. 47. Bluhm WF, Beran B, Bi C, Dimitropoulos D, Prlic A, et al. (2011) Quality assurance for the query and distribution systems of the RCSB Protein Data Bank. Database (Oxford) 2011: bar003.
  48. 48. Andreeva A, Howorth D, Chothia C, Kulesha E, Murzin AG (2014) SCOP2 prototype: a new approach to protein structure mining. Nucleic Acids Res 42(Database issue): D310–314.
  49. 49. Andersson M, Malmendal A, Linse S, Ivarsson I, Forsén S, et al. (1997) Structural basis for the negative allostery between Ca(2+)- and Mg(2+)-binding in the intracellular Ca(2+)-receptor calbindin D9k. Protein Sci 6: 1139–1147.
  50. 50. Brodersen DE, Etzerodt M, Madsen P, Celis JE, Thøgersen HC, et al. (1998) EF-hands at atomic resolution: the structure of human psoriasin (S100A7) solved by MAD phasing. Structure 6: 477–489.
  51. 51. Verdino P, Westritschnig K, Valenta R, Keller W (2002) The cross-reactive calcium-binding pollen allergen, Phl p 7, reveals a novel dimer assembly. EMBO J 21: 5007–5016.
  52. 52. Hohenester E, Maurer P, Hohenadl C, Timpl R, Jansonius JN, et al. (1996) Structure of a novel extracellular Ca(2+)-binding module in BM-40. Nat Struct Biol 3: 67–73.
  53. 53. Declercq JP, Evrard C, Lamzin V, Parello J (1999) Crystal structure of the EF-hand parvalbumin at atomic resolution (0.91 A) and at low temperature (100 K). Evidence for conformational multistates within the hydrophobic core. Protein Sci 8: 2194–2204.
  54. 54. Wilson MA, Brunger AT (2000) The 1.0 Å crystal structure of Ca(2+)-bound calmodulin: an analysis of disorder and implications for functionally relevant plasticity. J Mol Biol 301: 1237–1256.
  55. 55. de Beer T, Hoofnagle AN, Enmon JL, Bowers RC, Yamabhai M, et al. (2000) Molecular mechanism of NPF recognition by EH domains. Nat Struct Biol 7: 1018–1022.
  56. 56. Iwasaki W, Sasaki H, Nakamura A, Kohama K, Tanokura M (2003) Metal-free and Ca2+-bound structures of a multidomain EF-hand protein, CBP40, from the lower eukaryote Physarum polycephalum. Structure 11: 75–85.
  57. 57. Jia J, Borregaard N, Lollike K, Cygler M (2001) Structure of Ca(2+)-loaded human grancalcin. Acta Crystallogr D Biol Crystallogr 57: 1843–1849.
  58. 58. Ng C, Jackson RA, Buschdorf JP, Sun Q, Guy GR, et al. (2008) Structural basis for a novel intrapeptidyl H-bond and reverse binding of c-Cbl-TKB domain substrates. EMBO J 27: 804–816.
  59. 59. Kobayashi N, Koshiba S, Inoue M, Kigawa T, Yokoyama S. Solution structure of mouse CGI-38 protein. To be Published.
  60. 60. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, et al. (2012) The Pfam protein families database. Nucleic Acids Res 40: D290–D301.
  61. 61. Sobolev V, Sorokine A, Prilusky J, Abola EE, Edelman M (1999) Automated analysis of interatomic contacts in proteins. Bioinformatics 15: 327–332.
  62. 62. Brandl M, Weiss MS, Jabs A, Sühnel J, Hilgenfeld R (2001) C-H…pi-interactions in proteins. J Mol Biol. 307: 357–377.
  63. 63. Derewenda ZS, Lee L, Derewenda U (1995) The occurrence of C-H…O hydrogen bonds in proteins. J Mol Biol 252: 248–262.
  64. 64. Jia J, Borregaard N, Lollike K, Cygler M (2001) Structure of Ca(2+)-loaded human grancalcin. Acta Crystallogr D Biol Crystallogr 57: 1843–1849.
  65. 65. Himmel DM, Mui S, O'Neall-Hennessey E, Szent-Györgyi AG, Cohen C (2009) The on-off switch in regulated myosins: different triggers but related mechanisms. J Mol Biol 394: 496–505.
  66. 66. Tanaka T, Ames JB, Harvey TS, Stryer L, Ikura M (1995) Sequestration of the membrane-targeting myristoyl group of recoverin in the calcium-free state. Nature 376: 444–447.
  67. 67. Ames JB, Ishima R, Tanaka T, Gordon JI, Stryer L, et al. (1997) Molecular mechanics of calcium-myristoyl switches. Nature 389: 198–202.
  68. 68. Spyracopoulos L, Li MX, Sia SK, Gagné SM, Chandra M, et al. (1997) Calcium-induced structural transition in the regulatory domain of human cardiac troponin C. Biochemistry. 36: 12138–12146.
  69. 69. Robertson IM, Pineda-Sanabria SE1, Holmes PC1, Sykes BD (2014) Conformation of the critical pH sensitive region of troponin depends upon a single residue in troponin I. Arch Biochem Biophys 552-553: 40–49.
  70. 70. Oleszczuk M, Robertson IM, Li MX, Sykes BD (2010) Solution structure of the regulatory domain of human cardiac troponin C in complex with the switch region of cardiac troponin I and W7: the basis of W7 as an inhibitor of cardiac muscle contraction. J Mol Cell Cardiol 48: 925–933.
  71. 71. Wang X, Li MX, Sykes BD (2002) Structure of the regulatory N-domain of human cardiac troponin C in complex with human cardiac troponin I147–163 and bepridil. J Biol Chem 277: 31124–31133.
  72. 72. Lakowski TM, Lee GM, Okon M, Reid RE, McIntosh LP (2007) Calcium-induced folding of a fragment of calmodulin composed of EF-hands 2 and 3. Protein Sci 16: 1119–1132.
  73. 73. Jiang L, Baehr W (2010) GCAP1 mutations associated with autosomal dominant cone dystrophy. Adv Exp Med Biol 664: 273–282.
  74. 74. Sokal I, Li N, Surgucheva I, Warren MJ, Payne AM, et al. (1998) GCAP1 (Y99C) mutant is constitutively active in autosomal dominant cone dystrophy. Mol Cell 2: 129–133.
  75. 75. Nishiguchi KM, Sokal I, Yang L, Roychowdhury N, Palczewski K, et al. (2004) A novel mutation (I143NT) in guanylate cyclase-activating protein 1 (GCAP1) associated with autosomal dominant cone degeneration. Invest Ophthalmol Vis Sci 45: 3863–3870.
  76. 76. Koltzscher M, Gerke V (2000) Identification of hydrophobic amino acid residues involved in the formation of S100P homodimers in vivo. Biochemistry 39: 9533–9539.
  77. 77. Kragelund BB, Jönsson M, Bifulco G, Chazin WJ, Nilsson H, et al. (1998) Hydrophobic core substitutions in calbindin D9k: effects on Ca2+ binding and dissociation. Biochemistry 37: 8926–8937.
  78. 78. Tikunova SB, Rall JA, Davis JP (2002) Effect of hydrophobic residue substitutions with glutamine on Ca2+ binding and exchange with the N-domain of troponin C. Biochemistry 41: 6697–6705.
  79. 79. Takahashi T, Suzuki H, Inuzuka T, Shibata H, Maki M (2012) Prediction of a new ligand-binding site for type 2 motif based on the crystal structure of ALG-2 by dry and wet approaches. Int J Mol Sci 13: 7532–7549.
  80. 80. da Silva AC, Reinach FC (1991) Calcium binding induces conformational changes in muscle regulatory proteins. Trends Biochem Sci 16: 53–57.
  81. 81. Nelson MR, Chazin WJ (1998) Structures of EF-hand Ca(2+)-binding proteins: diversity in the organization, packing and response to Ca2+ binding. Biometals 11: 297–318.
  82. 82. Santamaria-Kisiel L, Rintala-Dempsey AC, Shaw GS (2006) Calcium-dependent and -independent interactions of the S100 protein family. Biochem J 396: 201–214.
  83. 83. Babini E, Bertini I, Borsi V, Calderone V, Hu X, et al. (2011) Structural characterization of human S100A16, a low-affinity calcium binder. J Biol Inorg Chem 16: 243–256.
  84. 84. Blanchard H, Grochulski P, Li Y, Arthur JS, Davies PL, et al. (1997) Structure of a calpain Ca2+-binding domain reveals a novel EF-hand and Ca2+-induced conformational changes. Nat Struct Biol 4: 532–538.
  85. 85. Iwasaki W, Sasaki H, Nakamura A, Kohama K, Tanokura M (2003) Metal-free and Ca2+-bound structures of a multidomain EF-hand protein, CBP40, from the lower eukaryote Physarum polycephalum. Structure 11: 75–85.
  86. 86. Lee YC, Volk DE, Thiviyanathan V, Kleerekoper Q, Gribenko AV, et al. (2004) NMR structure of the apo-S100P protein. J Biomol NMR 29: 399–402.
  87. 87. Sigrist CJ, de Castro E, Cerutti L, Cuche BA, Hulo N, et al. (2013) New and continuing developments at PROSITE. Nucleic Acids Res 41(Database issue): D344–347.
  88. 88. Maiti R, Van Domselaar GH, Zhang H, Wishart DS (2004) SuperPose: a simple server for sophisticated structural superposition. Nucleic Acids Res 32: W590–W594.
  89. 89. Kraulis PJ (1991) MOLSCRIPT: a program to produce both detailed and schematic plots of protein structures. J Appl Cryst 24: 946–950.
  90. 90. Merritt EA, Bacon DJ (1997) Raster3D: photorealistic molecular graphics. Methods Enzymol 277: 505–524.