Elucidating the mechanisms of cooperative calcium-calmodulin interactions: a structural systems biology approach
- Najl V Valeyev^{1}Email author,
- Declan G Bates^{1},
- Pat Heslop-Harrison^{1, 2},
- Ian Postlethwaite^{1} and
- Nikolay V Kotov^{3}
https://doi.org/10.1186/1752-0509-2-48
© Valeyev et al; licensee BioMed Central Ltd. 2008
Received: 30 July 2007
Accepted: 02 June 2008
Published: 02 June 2008
Abstract
Background
Calmodulin is an important multifunctional molecule that regulates the activities of a large number of proteins in the cell. Calcium binding induces conformational transitions in calmodulin that make it specifically active to particular target proteins. The precise mechanisms underlying calcium binding to calmodulin are still, however, quite poorly understood.
Results
In this study, we adopt a structural systems biology approach and develop a mathematical model to investigate various types of cooperative calcium-calmodulin interactions. We compare the predictions of our analysis with physiological dose-response curves taken from the literature, in order to provide a quantitative comparison of the effects of different mechanisms of cooperativity on calcium-calmodulin interactions. The results of our analysis reduce the gap between current understanding of intracellular calmodulin function at the structural level and physiological calcium-dependent calmodulin target activation experiments.
Conclusion
Our model predicts that the specificity and selectivity of CaM target regulation is likely to be due to the following factors: variations in the target-specific Ca^{2+} dissociation and cooperatively effected dissociation constants, and variations in the number of Ca^{2+} ions required to bind CaM for target activation.
Background
Calmodulin (CaM) is a multisite and multifunctional protein that contains four EF-hand Ca^{2+} binding sites [1], and is involved in a wide variety of cellular functions [2]. For example, it regulates the concentration of intracellular cAMP concentration in a very complex manner by regulating activities of cAMP producing adenylate cyclases (AC) and cAMP hydrolysing enzyme phosphodiesterase (PDE). CaM also regulates a large number of kinases and phosphatases as well as other enzymes with opposing cellular effects. Despite a large number of experimental studies [3–15], the detailed mechanisms underlying CaM-dependent intracellular regulation of such a large variety of target proteins are still not fully understood.
Recent experimental studies have attempted to elucidate the mechanisms underlying Ca^{2+}/CaM-dependent target regulation by measuring the kinetics and steady-state levels of CaM-target binding [3, 11, 31–33] as well as by analysing the mechanisms of Ca^{2+}-CaM interactions [9, 12, 34–37]. Ca^{2+} ion binding to EF-hand sites was shown to lead to CaM conformational alterations [1, 38–41]. In the modified conformational state, CaM is likely to alter its affinity to different targets by increasing and decreasing its affinity to certain proteins. Ca^{2+} ion binding to CaM is also argued to positively modulate the affinity of other Ca^{2+} binding sites of the molecule. There is still, however, an ongoing debate about the existence, the mechanisms and the degree of cooperativity in Ca^{2+}-CaM interactions. In some studies, Ca^{2+} binding to CaM has been reported to be independent [15, 37]. On the other hand, other studies have reported cooperative interactions between the neighbouring EF-hand binding sites [12, 35] or cooperativity linking all sites of the CaM molecule [9].
A number of previous studies have attempted to use mathematical modelling to obtain a quantitative understanding of the mechanisms involved in Ca^{2+}-CaM interactions. Different mathematical models, including the well known Hill [42], Adair [43] and Monod-Wyman-Changeux (MWC) [44] models, have been used in the literature to describe the cooperativity of ligand binding to a multisite protein. The Hill equation is frequently used to qualitatively measure the degree of cooperativity in multisite binding. It describes the simultaneous binding of n ligand molecules to a protein where the parameter n can be interpreted as the number of bound molecules. The Adair model represents ligand-protein interactions in terms of successive binding steps. The MWC model is based on two conformations that are in equilibrium and have different affinities for a ligand. To date, the Hill and Adair models have been most frequently used to investigate Ca^{2+} binding to CaM [4, 9–12, 14, 33, 45]. While these studies have provided much useful information, the use of the classical models mentioned above also introduces some limitations in the analysis – see [46] for a full discussion of this issue. In particular, the detailed analysis of ligand-protein interactions which are unique to CaM requires the development of a model that captures multiple functionally important intermediate conformations of the protein. A steady-state solution to the cooperativity problem for Ca2+ binding has been analysed in [47]. In this paper, we develop a new model, based on the assumption that the specificity in CaM target regulation arises from the Ca^{2+}-CaM complex specific target interactions with variable numbers of bound Ca^{2+} ions. In this approach, Ca^{2+} binding to each EF-hand sites causes conformational transitions in the CaM molecule leading to a model that has multiple conformational states in complex with variable numbers of Ca^{2+} ions. In the proposed model, CaM may regulate its targets with one, two or three Ca^{2+} ions as well as in the apo- or fully bound states. In particular, we address the Ca^{2+}-CaM interaction in significant detail, although we do not incorporate detailed Ca^{2+}-CaM species interactions with target proteins. This approach is in agreement with recent experimental evidence that the concentration-dependent profiles for several Ca^{2+}-CaM-dependent protein targets exhibit quite a diverse range of behaviour. PMCA and PDE protein concentrations in the active state, for example, reveal "Hill-shape"-like curves, whereas the ACII isoform is inhibited by increasing Ca^{2+} concentration. The ACVI isoform exhibits inhibition with an interesting plateau feature on the Ca^{2+}-dependent profile. Yet ACI isoforms have bell-shaped concentration-dependent profiles [48]. It has also been shown that CaMPKII [49] as well as the K^{+} channel from Paramecium [50] are activated by CaM with two bound Ca^{2+} ions.
The structural systems biology approach, [51], employed in this paper provides new insights into the Ca^{2+}-CaM-target binding dose-response curves which have been derived experimentally, and allows us to advance testable hypotheses about the nature of cooperative mechanisms unique to calcium-CaM interactions. The resulting analysis further bridges the gap between our understanding of CaM structural properties and intracellular Ca^{2+}-CaM-dependent target regulation.
Results
Multisite binding of cooperatively linked binding centres
In order to investigate the dependence of one Ca^{2+} binding site on another in both the N- and C-terminal domains, we assumed the alteration of a dissociation constant when a neighboring site is occupied as illustrated schematically in Figure 3 (compare the Model 1 and Model 2 descriptions in the Materials and Methods section). This approximation allows the derivation of a model that has an analytical solution in the form of conditional probabilities (Equation 12). It provides a quantitative comparison of the concentration of Ca^{2+} bound to CaM in the presence and absence of Ca^{2+} binding site cooperative interactions (Equations 13–14). Figure 5A shows the model predictions in the case where CaM is assumed to have two pairs of independent EF-hand globular domains. Within these domains, one Ca^{2+} binding site influences the other. In the N-terminal domain, the affinity of the second site depends on the state of the first and changes from K_{2} = 0.9 to ${K}_{2}^{coop}$ = 0.2 (μM l^{-1}) when a Ca^{2+} ion occupies the first centre. In the C-terminal, the affinity of the fourth site depends on the state of the third and changes from K_{4} = 0.8 to ${K}_{4}^{coop}$ = 0.1 (μM l^{-1}) when a Ca^{2+} ion occupies the third centre. The model (Model 2 in Materials and Methods) predicts that such changes will mainly influence the "amplitudes" of the intermediate conformations of the concentration-dependent profiles while leaving the ligand concentrations that produce their maximum values largely unchanged. Figure 5B shows the difference in the total amount of bound ligand with and without the type of cooperativity described above. Since the amount of bound Ca^{2+} is frequently measured in Ca^{2+}-CaM or other ligand-multisite protein interaction experiments, these results allow a direct quantitative comparison to be made of the binding reactions with and without the presence of cooperative binding.
N- and C-terminal domains reveal unique cooperative properties
A more realistic description for EF-hand Ca^{2+} binding sites would involve the incorporation of the influence of both Ca^{2+} binding sites on each other within the CaM globular domains as schematically illustrated on Figure 5 (Model 3 in Materials and Methods). The resulting system of differential equations describing the Ca^{2+}-CaM interactions is given by (Equation 15) in Materials and Methods. The model based on these assumptions leads to some interesting predictions regarding Ca^{2+}-dependent CaM interactions with various CaM target protein peptides [33, 49] and reveals a complex story of specificity in CaM regulation. While it is well established that Ca^{2+} ions are required to modulate the CaM-target protein interactions, the mechanism of Ca^{2+}-induced CaM conformational transitions that allow selective interactions with a particular target protein is still unclear. The presented model for Ca^{2+} binding to CaM provides new insights into how the cooperative interactions between EF-hand binding sites contribute to the mechanism of selective target regulation by CaM, as described below.
Calmodulin-target interactions reveal target-specific cooperativity in Ca^{2+} binding
The results of the present study, when combined with previous experimental data from the literature, suggest that CaM interacts with phosphorylase kinase, CaATPase and skMLCK in the apo state, but activates these proteins only when Ca^{2+} ions bind to a CaM-target protein complex. CaMKII kinase, on the other hand, binds to the Ca^{2+}-CaM complex rather than apo-CaM. Each kinase has a unique combination of K_{1} and K_{2} dissociation constants. The half length peptides of CaATPase and skMLCK appear to have an even more complex mechanism of binding. As mentioned earlier, each Ca^{2+}-CaM complex (with variable numbers of Ca^{2+} ions bound) may have a unique set of K_{1} and K_{2} constants for a target protein or peptide. For simplicity, all possible CaM species were divided into two groups: i) those with less than 3 Ca^{2+} ions bound and ii) those with three and four Ca^{2+} ions bound to CaM. A comparison of the model predictions with the Ca^{2+}-CaM-parts of skMLCK and CaATPase peptide binding data (Figure 7B and 7C) suggests that both the K_{1} and the cooperatively influenced K_{2} dissociation constants are different when peptides are bound to CaM species with less than 3 Ca^{2+} ions or to CaM species with 3 or 4 ions. The conclusion from this observation is that the specificity in Ca^{2+}-CaM-dependent regulation arises from a combination of the target specific affinity between Ca^{2+} and CaM, target specific cooperative constants, the order of the Ca^{2+}-CaM-target complex assembly, as well as the number of Ca^{2+} ions bound to CaM. All these factors contribute to the mechanism of selective Ca^{2+}-CaM dependent regulation in addition to the diversity of CaM-target interfaces [39].
pH dependence of cooperative Ca^{2+} binding to CaM
Ca^{2+}-CaM dissociation constants derived by the different mathematical models.
Proteins | Hill | Adair | present model | Reference |
---|---|---|---|---|
PhK5 | K_{D 1}= 0.24, K_{D 2}= 13 | K_{1} = 1, Kc_{1} = 1, K_{2} = 1, Kc_{2} = 1 | [33] | |
skMLCK | K_{D 1}= 0.02, K_{D 2}= 0.08 | K_{1} = 0.04, Kc_{1} = 0.02, K_{2} = 0.04, Kc_{2} = 0.02 | ||
sk-N11 | K_{D 1}= 0.26, K_{D 2}= 6 | K_{1} = 1.2, K_{2} = 0.5 K_{1} = 5, K_{2} = 1 | ||
sk-C10 | K_{D 1}= 3.4, K_{D 2}= 4 | K_{1} = 0.06, Kc_{1} = 0.02, K_{2} = 0.06, Kc_{2} = 0.02 | ||
CaATPase | K_{D 1}= 0.09, K_{D 2}= 0.2 | K_{1} = 0.15, Kc_{1} = 0.05, K_{2} = 0.15, Kc_{2} = 0.05 | ||
ATPase-N18 | K_{D 1}= 0.12, K_{D 2}= 3.9 | K_{1} = 2, Kc_{1} = 1, K_{2} = 2, Kc_{2} = 1 | ||
ATPase-C17 | K_{D 1}= 0.66, K_{D 2}= 2.4 | K_{1} = 0.4, Kc_{1} = 0.2 K_{2} = 2, Kc_{2} = 1 | ||
CaMKII-cbp | K_{1} = 0.5, Kc_{1} = 0.5, K_{2} = 0.5, Kc_{2} = 0.5 | [49] | ||
CaMKII | K_{1} = 5, Kc_{1} = 5, K_{2} = 5, Kc_{2} = 5 | |||
CaM pH = 7.2 | K_{1} = 0.34, K_{2} = 0.36, K_{3} = 0.13, K_{4} = 0.06 | K_{1} = 17, Kc_{1} = 7, K_{2} = 20, Kc_{2} = 0.5 | [12] | |
F12 | K_{1} = 0.142, K_{2} = 0.062 | K_{1} = 17, Kc_{1} = 7, K_{2} = 17, Kc_{2} = 7 | ||
F34 | K_{3} = 0.0543, K_{4} = 1.82 | K_{1} = 20, Kc_{1} = 0.5, K_{2} = 20, Kc_{2} = 0.5 | ||
CaM pH = 6 | K_{1} = 10, Kc_{1} = 5, K_{2} = 10, Kc_{2} = 5 | [9] | ||
CaM pH = 10.1 | K_{1} = 2, Kc_{1} = 1.8, K_{2} = 2, Kc_{2} = 1.8 |
Discussion
In this paper a structural model of CaM interactions with and without cooperativity has been used to elucidate the mechanisms of Ca^{2+}-CaM-target complex assembly. The differences seen in dose-response curves for proteins activated by Ca^{2+}-CaM pairs were explained in terms of cooperative interactions between the EF-hand pairs [38, 40, 41] in both CaM domains. This study predicts that the specific interaction interface between CaM and CaM-regulated proteins [8, 10, 39] is complemented by a number of additional factors influencing the Ca^{2+}-CaM-target complex assembly. By comparing our model predictions with experimentally measured dose-response curves from the literature, we propose that some proteins bind CaM without Ca^{2+} ions and only become activated when Ca^{2+} ions interact with the CaM-target complex, whereas others are activated by CaM molecules with already bound Ca^{2+} ions. The Ca^{2+}-CaM interaction properties are tuned by the target proteins and characterized not only by the macroscopic dissociation constant set for Ca^{2+} sites, but also by the macroscopic cooperatively altered dissociation constants that are also unique to the CaM binding proteins. In other words, the order of Ca^{2+}-CaM-target complex assembly, the number of bound Ca^{2+} ions, target specific Ca^{2+}-CaM cooperative affinities, in addition to unique CaM-target interaction interfaces, all allow CaM to achieve its highly versatile intracellular multifunctionality. This proposition also explains the effects of pH on the considered dose-response curves by allowing for the modulation of the cooperatively effected dissociation constants. We would also like to point out that while we addressed the Ca^{2+}-CaM interactions in great detail, the model could still be developed further by incorporating detailed dissociation constants between the intermediate Ca^{2+}-CaM complexes and target proteins in a similar way to how it has been done for the Ca^{2+}-CaM interactions.
Although the presented model predicts similar curves to the ones already used to approximate the experimental Ca^{2+}-CaM-target dose-response data using the Hill and Adair models, it allows a far more detailed interpretation of the Ca^{2+}-CaM dependent interactions involved. In particular, it distinguishes the structure-dependent properties of CaM molecules and suggests potential scenarios for Ca^{2+}-CaM-target complex assembly. Importantly, it reveals the CaM specific type of cooperativity involved in this process and helps to explain the contribution of cooperativity in the specificity of CaM-dependent regulation.
Conclusion
- (1)
Mathematical models for protein interactions are usually derived according to a number of assumptions which will inevitably be more or less applicable to each particular protein. The structure of CaM suggests that this molecule is very likely to have non-sequential Ca^{2+} access to EF-hand binding sites. The results of our analysis support the theory of non-sequential cooperative access of Ca^{2+} to CaM binding sites, and also allow the derivation of cooperatively effected dissociation constants, thus providing a more realistic tool for fitting experimental dose-response curves.
- (2)
Our model suggests that the structural data alone cannot provide the required level of information and comparisons with the dose-response data are required. Predictions from the mathematical model used in this study were compared with the dose-response curves for Ca^{2+} binding to CaM and Ca^{2+}-CaM-target peptides. This analysis allowed us to distinguish between (a) proteins that form a complex with CaM in its Ca^{2+} free state and then interact with Ca^{2+} ions and (b) other proteins which interact with Ca^{2+} bound CaM with variable numbers of Ca^{2+} ions. However, the transient kinetics has not been addressed in this study.
- (3)
In addition to the diversity of interaction interfaces, the specificity and selectivity of CaM target activation may be achieved by variations in the target-specific dissociation and cooperatively effected dissociation constants, the order of Ca^{2+}-CaM-target complex assembly and the number of Ca^{2+} ions required to bind CaM for target activation.
Methods
A structural mathematical model for Ca^{2+}-CaM activation
Here, we present the mathematical equations used to describe Ca^{2+}-CaM interactions in our study. To clarify how various factors contribute to the CaM-dependent regulation, we describe three models for CaM starting from a very basic approximation of completely independent Ca^{2+} binding sites and gradually progressing to more realistic models that take account of cooperativity mechanisms and the Ca^{2+}-CaM-target peptide complex assembly.
Model 1. CaM with independent Ca^{2+} binding sites
where cm 0 is the total concentration of CaM, u is normalised Ca^{2+} concentration,${p}_{i}^{0}(u)$ is the probability of binding site i not being occupied and ${p}_{i}^{1}(u)$ is the probability of binding site i being bound. c_{ i }equals 1 if a binding site is occupied and 0 if it is not. The probability of CaM being in a particular bound state is equal to the product of the probabilities of each individual binding site.
where the K and u are the microscopic equilibrium dissociation constant and the ligand concentration, respectively. Effectively these are Michaelis-Menten equations for a protein in a complex with and without a ligand molecule, but normalized by the total protein concentration.
where the K and u are the equilibrium dissociation constant and the ligand concentration, respectively.
A complete mathematical description of the relationship between macroscopic constants derived from the Adair equation and the proposed model is provided in the Supplementary Materials Section (Additional File 1).
where L 0 is the total protein concentration and L_{ i }is the concentration of conformation i, K_{ j }are the equilibrium dissociation constants of each binding site, and u is the ligand concentration.
The conformations L_{1}(u), ..., L_{n-1}(u) of a multisite protein are all bell-shaped curves, the conformation L_{0}(u) is the apo state of a multisite protein, whereas the L_{ n }(u) is the fully bound multisite protein. If equation (5) is divided by L 0, then instead of predicting protein concentrations in specific ligand-bound conformations, it predicts the probability of a particular conformation to be in that state as a function of ligand concentration.
where cm_{ i }(u) is substituted from equation 5, and K_{d} is the dissociation constant for CaM-target interactions, N 0 is the total concentration of target protein.
Model 2. Cooperative Ca^{2+}-CaM interactions
While the previous model provides predictions for the number of Ca^{2+}-CaM complexes as a function of Ca^{2+} concentration with a reasonable accuracy, it does not capture the effects of the cooperative influence of Ca^{2+} binding sites. There are several possible ways to incorporate these cooperative mechanisms into the model. In order to derive a model that illustrates what contribution cooperativity makes to the distribution of concentration profiles of Ca^{2+}-CaM complexes, we assume that in the N-terminal domain, the first centre is cooperatively bound to the second, and in the C-terminal, the third is cooperatively bound to the fourth. In this case, we will define the dissociation constants as K_{1}, K_{2}, ${K}_{2}^{coop}$, K_{3}, K_{4}, ${K}_{4}^{coop}$, where ${K}_{2}^{coop}$ and ${K}_{4}^{coop}$ are the cooperatively influenced dissociation constants for the second and fourth centres when the ligand is bound to the first and third binding sites, correspondingly.
where K_{1} and u are the equilibrium dissociation constant for the first centre and the ligand concentration, respectively.
where ${p}_{2}^{00}$ is the probability for both the first and the second centres to be free, ${p}_{2}^{01}$ is the probability for the first site to be free and the second to be occupied, ${p}_{2}^{10}$ is the probability for the first site to be bound and the second to be free, and ${p}_{2}^{11}$ is the probability for both sites to be ligand bound. K_{2} and ${K}_{2}^{coop}$ are the dissociation and the cooperatively modified dissociation constants for the second Ca^{2+} binding site.
where K_{3} is a dissociation of the third binding site.
where ${p}_{4}^{00}$ is the probability for both the third and the fourth centres to be free, ${p}_{4}^{01}$ is the probability for the third site to be free and the fourth to be occupied, ${p}_{3}^{10}$ is the probability for the third site to be bound and the fourth to be free, and ${p}_{4}^{11}$ is the probability for both sites to be ligand bound. K_{4} and ${K}_{4}^{coop}$ are the dissociation and the cooperatively modified dissociation constants for the second Ca^{2+} binding site.
where k_{i} = 0, if the binding site i is not occupied and k_{i} = 1 if the centre i is occupied by a ligand molecule.
The distribution of intermediate Ca^{2+}-CaM complexes with 1, 2 and 3 Ca^{2+} ions with and without cooperativity is shown in Figure 5A.
Figure 5B shows the total amount of ligand bound to CaM in the presence (1) and absence (2) of cooperative binding. The line (3) shows the difference in the level of bound ligand between the two types of binding mechanisms.
Model 3. The Ca^{2+}-CaM-target protein complex assembly
where cm_{00}, cm_{01}, cm_{10}, cm_{11} are CaM molecules without Ca^{2+} ions, with one Ca^{2+} ion bound to the N-terminal domain, with one Ca^{2+} ion bound to the N-terminal domain, and CaM species with two bound Ca^{2+} ions at each terminal domain, respectively. k_{1} and k_{2} are the association constants and kc_{1} and kc_{2} are the cooperatively modified association constants for the N- terminal binding sites of CaM, respectively. Similarly, d_{1}, d_{2}, dc_{1}, dc_{2} are the dissociation and cooperatively modified dissociation constants for the N- and terminal binding sites, respectively. Note that a similar system of differential equations can be developed for the C-terminal.
The conservation law gives:
cm_{00} +cm_{10} + cm_{01} + cm_{11} = 1
In steady-state, the matrix for the system (15) with the last equation substituted by (16) is given by:
where ${K}_{1}=\frac{{d}_{1}}{{k}_{1}}$, $K{c}_{1}=\frac{d{c}_{1}}{kc}$. k_{1} = k_{2} = kc_{1} = kc_{2} is a simplifying assumption that allows analytical solution of the current system.
where N_{00}, N_{01}, N_{10} and N_{11} correspond to the N-terminal domain without any Ca^{2+} ions bound, with one Ca^{2+} ion bound to one or another binding site and to the fully bound state, respectively.
where ${p}_{{N}_{00}}(u)$ is the probability for the N-terminal to be in a Ca^{2+} state, ${p}_{{N}_{01}}(u)$ and ${p}_{{N}_{10}}(u)$ describe the probabilities for a complex with one Ca^{2+} ion, and ${p}_{{N}_{11}}(u)$ is the probability function for the N-terminal to be occupied by Ca^{2+}.
Where ${K}_{2}=\frac{{d}_{2}}{{k}_{2}},K{c}_{2}=\frac{d{c}_{2}}{k{c}_{2}},{k}_{2}=k{c}_{2}$
where ${p}_{{N}_{00}}(u)$ is the probability for the N-terminal to be in a Ca^{2+} state, ${p}_{{N}_{01}}(u)$ and ${p}_{{N}_{10}}(u)$ describe the probabilities for a complex with one Ca^{2+} ion, and ${p}_{{N}_{11}}(u)$ is the probability function for the N-terminal occupied by Ca^{2+}.
where $Ca{M}_{00}^{00}$ are CaM species in apo state, $Ca{M}_{11}^{11}$ is fully bound CaM, and $Ca{M}_{00}^{11}$ and $Ca{M}_{11}^{00}$ are CaM species with fully bound N- and C-terminals, respectively. $Ca{M}_{00}^{01}$ and $Ca{M}_{01}^{00}$ are CaM species with one Ca^{2+} ion bound to N- and C-terminals, respectively. Comparison of numerical solutions of this system with the available experimental data [9, 15, 33] allows us to propose that CaM molecules can be represented as a pair of two independent N- and C-terminal globular domains, each containing two symmetrical and cooperatively bound EF-hand Ca^{2+} binding sites.
The equations (28) are specific to Ca^{2+}-CaM interactions and incorporate the pairwise cooperative interactions between the EF-hand binding sites within the N- and C-terminals, whereas the N- and C-terminal domains are considered to be independent of each other. Note that the equations developed here for the Ca^{2+}-CaM complexes are essentially different from the Michelis-Menten, Hill and Adair models, and also differ from models with independent binding sites (5) or with limited amounts of cooperativity (11) and (12).
where T 0 is the total concentration of a target protein or a peptide, K_{ d }is the equilibrium dissociation constant between CaM and a target protein, and $\sum Ca{M}_{m,l}^{i,j}(u)$ is a single Ca^{2+} complex or the sum of several CaM complexes.
The combinations of Ca^{2+}-CaM complexes have been varied simultaneously with the dissociation constants to fit the experimental data. This analysis allows us to predict the Ca^{2+}-CaM complexes required for activation of specific protein targets. The fitting of dissociation constants of Ca^{2+} binding sites on CaM molecules to the experimental dose-response curves reveals the impact of the target protein on the Ca^{2+}-CaM interactions. The dissociation constants calculated based on the cooperative Model 3, which also takes into account the impact of target proteins, are compared with the dissociation constants calculated using the Hill and modified Adair equation in the original experimental publications for the same data in Table 1.
Declarations
Acknowledgements
This work was carried out under EPSRC platform grant (EP/D029937/1), BBSRC grant (BB/D015340/1) and by RFBR grant.
Authors’ Affiliations
References
- Yap KL, Ames JB, Swindells MB, Ikura M: Diversity of conformational states and changes within the EF-hand protein superfamily. Proteins. 1999, 37: 499-507. 10.1002/(SICI)1097-0134(19991115)37:3<499::AID-PROT17>3.0.CO;2-YView ArticlePubMedGoogle Scholar
- Berridge MJ, Bootman MD, Lipp P: Calcium–a life and death signal. Nature. 1998, 395: 645-648. 10.1038/27094View ArticlePubMedGoogle Scholar
- Guo Q, Shen Y, Lee YS, Gibbs CS, Mrksich M, Tang WJ: Structural basis for the interaction of Bordetella pertussis adenylyl cyclase toxin with calmodulin. Embo J. 2005, 24: 3190-3201. 10.1038/sj.emboj.7600800PubMed CentralView ArticlePubMedGoogle Scholar
- Grabarek Z: Structure of a trapped intermediate of calmodulin: calcium regulation of EF-hand proteins from a new perspective. J Mol Biol. 2005, 346: 1351-1366. 10.1016/j.jmb.2005.01.004View ArticlePubMedGoogle Scholar
- Elshorst B, Hennig M, Forsterling H, Diener A, Maurer M, Schulte P, Schwalbe H, Griesinger C, Krebs J, Schmid H, et al.: NMR solution structure of a complex of calmodulin with a binding peptide of the Ca2+ pump. Biochemistry. 1999, 38: 12320-12332. 10.1021/bi9908235View ArticlePubMedGoogle Scholar
- Drum CL, Yan SZ, Bard J, Shen YQ, Lu D, Soelaiman S, Grabarek Z, Bohm A, Tang WJ: Structural basis for the activation of anthrax adenylyl cyclase exotoxin by calmodulin. Nature. 2002, 415: 396-402. 10.1038/415396aView ArticlePubMedGoogle Scholar
- Haiech J, Klee CB, Demaille JG: Effects of cations on affinity of calmodulin for calcium: ordered binding of calcium ions allows the specific activation of calmodulin-stimulated enzymes. Biochemistry. 1981, 20: 3890-3897. 10.1021/bi00516a035View ArticlePubMedGoogle Scholar
- Hoeflich KP, Ikura M: Calmodulin in action: diversity in target recognition and activation mechanisms. Cell. 2002, 108: 739-742. 10.1016/S0092-8674(02)00682-7View ArticlePubMedGoogle Scholar
- Iida S, Potter JD: Calcium binding to calmodulin. Cooperativity of the calcium-binding sites. J Biochem. 1986, 99 (6): 1765-1772.PubMedGoogle Scholar
- Ikura M: Calcium binding and conformational response in EF-hand proteins. Trends Biochem Sci. 1996, 21: 14-17.View ArticlePubMedGoogle Scholar
- Mirzoeva S, Weigand S, Lukas TJ, Shuvalova L, Anderson WF, Watterson DM: Analysis of the functional coupling between Calmodulin's calcium binding and peptide recognition properties. Biochemistry. 1999, 38: 14117-14118. 10.1021/bi9950894View ArticlePubMedGoogle Scholar
- Minowa O, Yagi K: Calcium binding to tryptic fragments of calmodulin. J Biochem. 1984, 96 (4): 1175-1182.PubMedGoogle Scholar
- Meador WE, Means AR, Quiocho FA: Target enzyme recognition by calmodulin: 2.4 A structure of a calmodulin-peptide complex. Science. 1992, 257: 1251-1255. 10.1126/science.1519061View ArticlePubMedGoogle Scholar
- Maune JF, Klee CB, Beckingham K: Ca2+ binding and conformational change in two series of point mutations to the individual Ca(2+)-binding sites of calmodulin. J Biol Chem. 1992, 267: 5286-5295.PubMedGoogle Scholar
- Ogawa Y, Tanokura M: Calcium binding to calmodulin: effects of ionic strength, Mg2+, pH and temperature. J Biochem. 1984, 95 (1): 19-28.PubMedGoogle Scholar
- Andre I, Kesvatera T, Jonsson B, Akerfeldt KS, Linse S: The role of electrostatic interactions in calmodulin-peptide complex formation. Biophys J. 2004, 87: 1929-1938. 10.1529/biophysj.104.040998PubMed CentralView ArticlePubMedGoogle Scholar
- Andre I, Kesvatera T, Jonsson B, Linse S: Salt enhances calmodulin-target interaction. Biophys J. 2006, 90: 2903-2910. 10.1529/biophysj.105.068718PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang M, Tanaka T, Ikura M: Calcium-induced conformational transition revealed by the solution structure of apo calmodulin. Nat Struct Biol. 1995, 2: 758-767. 10.1038/nsb0995-758View ArticlePubMedGoogle Scholar
- Chattopadhyaya R, Meador WE, Means AR, Quiocho FA: Calmodulin structure refined at 1.7 A resolution. J Mol Biol. 1992, 228: 1177-1192. 10.1016/0022-2836(92)90324-DView ArticlePubMedGoogle Scholar
- Ohki S, Iwamoto U, Aimoto S, Yazawa M, Hikichi K: Mg2+ inhibits formation of 4Ca(2+)-calmodulin-enzyme complex at lower Ca2+ concentration. 1H and 113Cd NMR studies. J Biol Chem. 1993, 268: 12388-12392.PubMedGoogle Scholar
- Ohki S, Ikura M, Zhang M: Identification of Mg2+-binding sites and the role of Mg2+ on target recognition by calmodulin. Biochemistry. 1997, 36: 4309-4316. 10.1021/bi962759mView ArticlePubMedGoogle Scholar
- Osawa M, Tokumitsu H, Swindells MB, Kurihara H, Orita M, Shibanuma T, Furuya T, Ikura M: A novel target recognition revealed by calmodulin in complex with Ca2+-calmodulin-dependent kinase kinase. Nat Struct Biol. 1999, 6: 819-824. 10.1038/12271View ArticlePubMedGoogle Scholar
- Matsubara M, Nakatsu T, Kato H, Taniguchi H: Crystal structure of a myristoylated CAP-23/NAP-22 N-terminal domain complexed with Ca2+/calmodulin. Embo J. 2004, 23: 712-718. 10.1038/sj.emboj.7600093PubMed CentralView ArticlePubMedGoogle Scholar
- Houdusse A, Gaucher JF, Krementsova E, Mui S, Trybus KM, Cohen C: Crystal structure of apo-calmodulin bound to the first two IQ motifs of myosin V reveals essential recognition features. Proc Natl Acad Sci USA. 2006, 103: 19326-19331. 10.1073/pnas.0609436103PubMed CentralView ArticlePubMedGoogle Scholar
- Bahler M, Rhoads A: Calmodulin signaling via the IQ motif. FEBS Lett. 2002, 513: 107-113. 10.1016/S0014-5793(01)03239-2View ArticlePubMedGoogle Scholar
- Schumacher MA, Crum M, Miller MC: Crystal structures of apocalmodulin and an apocalmodulin/SK potassium channel gating domain complex. Structure. 2004, 12: 849-860. 10.1016/j.str.2004.03.017View ArticlePubMedGoogle Scholar
- Schumacher MA, Rivard AF, Bachinger HP, Adelman JP: Structure of the gating domain of a Ca2+-activated K+ channel complexed with Ca2+/calmodulin. Nature. 2001, 410: 1120-1124. 10.1038/35074145View ArticlePubMedGoogle Scholar
- Yuan T, Vogel HJ: Substitution of the methionine residues of calmodulin with the unnatural amino acid analogs ethionine and norleucine: biochemical and spectroscopic studies. Protein Sci. 1999, 8: 113-121.PubMed CentralView ArticlePubMedGoogle Scholar
- Corneliussen B, Holm M, Waltersson Y, Onions J, Hallberg B, Thornell A, Grundstrom T: Calcium/calmodulin inhibition of basic-helix-loop-helix transcription factor domains. Nature. 1994, 368: 760-764. 10.1038/368760a0View ArticlePubMedGoogle Scholar
- Larsson G, Schleucher J, Onions J, Hermann S, Grundstrom T, Wijmenga SS: A novel target recognition revealed by calmodulin in complex with the basic helix–loop–helix transcription factor SEF2-1/E2-2. Protein Sci. 2001, 10: 169-186. 10.1110/ps.28401PubMed CentralView ArticlePubMedGoogle Scholar
- Bayley PM, Findlay WA, Martin SR: Target recognition by calmodulin: dissecting the kinetics and affinity of interaction using short peptide sequences. Protein Sci. 1996, 5: 1215-1228.PubMed CentralView ArticlePubMedGoogle Scholar
- Green DF, Dennis AT, Fam PS, Tidor B, Jasanoff A: Rational design of new binding specificity by simultaneous mutagenesis of calmodulin and a target peptide. Biochemistry. 2006, 45: 12547-12559. 10.1021/bi060857uPubMed CentralView ArticlePubMedGoogle Scholar
- Peersen OB, Madsen TS, Falke JJ: Intermolecular tuning of calmodulin by target peptides and proteins: differential effects on Ca2+ binding and implications for kinase activation. Protein Sci. 1997, 6: 794-807.PubMed CentralView ArticlePubMedGoogle Scholar
- Evenas J, Malmendal A, Thulin E, Carlstrom G, Forsen S: Ca2+ binding and conformational changes in a calmodulin domain. Biochemistry. 1998, 37: 13744-13754. 10.1021/bi9806448View ArticlePubMedGoogle Scholar
- Linse S, Helmersson A, Forsen S: Calcium binding to calmodulin and its globular domains. J Biol Chem. 1991, 266: 8050-8054.PubMedGoogle Scholar
- Tan RY, Mabuchi Y, Grabarek Z: Blocking the Ca2+-induced conformational transitions in calmodulin with disulfide bonds. J Biol Chem. 1996, 271: 7479-7483. 10.1074/jbc.271.13.7479View ArticlePubMedGoogle Scholar
- Wang CL: A note on Ca2+ binding to calmodulin. Biochem Biophys Res Commun. 1985, 130: 426-430. 10.1016/0006-291X(85)90434-6View ArticlePubMedGoogle Scholar
- Gifford JL, Walsh MP, Vogel HJ: Structures and metal-ion-binding properties of the Ca2+-binding helix-loop-helix EF-hand motifs. Biochem J. 2007, 405: 199-221. 10.1042/BJ20070255View ArticlePubMedGoogle Scholar
- Bhattacharya S, Bunick CG, Chazin WJ: Target selectivity in EF-hand calcium binding proteins. Biochim Biophys Acta. 2004, 1742: 69-79. 10.1016/j.bbamcr.2004.09.002View ArticlePubMedGoogle Scholar
- Nelson MR, Thulin E, Fagan PA, Forsen S, Chazin WJ: The EF-hand domain: a globally cooperative structural unit. Protein Sci. 2002, 11: 198-205. 10.1110/ps.33302PubMed CentralView ArticlePubMedGoogle Scholar
- Nelson MR, Chazin WJ: An interaction-based analysis of calcium-induced conformational changes in Ca2+ sensor proteins. Protein Sci. 1998, 7: 270-282.PubMed CentralView ArticlePubMedGoogle Scholar
- Hill A: The combinations of haemoglobin with oxygen and with carbon monoxide. J Physiology. 1910, 40: 4-7.Google Scholar
- Adair G: The hemolglobin system. The oxygen dissociation curve of hemoglobin. J Biol Chem. 1925, 63: 529-545.Google Scholar
- Monod J, Wyman J, Changeux JP: On the Nature of Allosteric Transitions: A Plausible Model. J Mol Biol. 1965, 12: 88-118.View ArticlePubMedGoogle Scholar
- Porumb T: Determination of calcium-binding constants by flow dialysis. Anal Biochem. 1994, 220: 227-237. 10.1006/abio.1994.1332View ArticlePubMedGoogle Scholar
- Weiss JN: The Hill equation revisited: uses and misuses. Faseb J. 1997, 11: 835-841.PubMedGoogle Scholar
- Haiech J, Kilhoffer MC: Deconvolution of calcium-binding curves. Facts and fantasies. Methods Mol Biol. 2002, 173: 25-42.PubMedGoogle Scholar
- Guillou JL, Nakata H, Cooper DM: Inhibition by calcium of mammalian adenylyl cyclases. J Biol Chem. 1999, 274: 35539-35545. 10.1074/jbc.274.50.35539View ArticlePubMedGoogle Scholar
- Shifman JM, Choi MH, Mihalas S, Mayo SL, Kennedy MB: Ca2+/calmodulin-dependent protein kinase II (CaMKII) is activated by calmodulin with two bound calciums. Proc Natl Acad Sci USA. 2006, 103: 13968-13973. 10.1073/pnas.0606433103PubMed CentralView ArticlePubMedGoogle Scholar
- Kung C, Preston RR, Maley ME, Ling KY, Kanabrocki JA, Seavey BR, Saimi Y: In vivo Paramecium mutants show that calmodulin orchestrates membrane responses to stimuli. Cell Calcium. 1992, 13: 413-425. 10.1016/0143-4160(92)90054-VView ArticlePubMedGoogle Scholar
- Aloy P, Russell RB: Structural systems biology: modelling protein interactions. Nat Rev Mol Cell Biol. 2006, 7: 188-197. 10.1038/nrm1859View ArticlePubMedGoogle Scholar
- Stemmer PM, Klee CB: Dual calcium ion regulation of calcineurin by calmodulin and calcineurin B. Biochemistry. 1994, 33: 6859-6866. 10.1021/bi00188a015View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Comments
View archived comments (1)