Elucidating the mechanisms of cooperative calcium-calmodulin interactions: a structural systems biology approach

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 Ca2+ dissociation and cooperatively effected dissociation constants, and variations in the number of Ca2+ ions required to bind CaM for target activation

A Dynamic Model of Interactions of Ca^(2+), Calmodulin, and Catalytic Subunits of Ca^(2+)/Calmodulin-Dependent Protein Kinase II

During the acquisition of memories, influx of Ca^(2+) into the postsynaptic spine through the pores of activated N-methyl-D-aspartate-type glutamate receptors triggers processes that change the strength of excitatory synapses. The pattern of Ca^(2+) influx during the first few seconds of activity is interpreted within the Ca^(2+)-dependent signaling network such that synaptic strength is eventually either potentiated or depressed. Many of the critical signaling enzymes that control synaptic plasticity, including Ca^(2+)/calmodulin-dependent protein kinase II (CaMKII), are regulated by calmodulin, a small protein that can bind up to 4 Ca^(2+) ions. As a first step toward clarifying how the Ca^(2+)-signaling network decides between potentiation or depression, we have created a kinetic model of the interactions of Ca^(2+), calmodulin, and CaMKII that represents our best understanding of the dynamics of these interactions under conditions that resemble those in a postsynaptic spine. We constrained parameters of the model from data in the literature, or from our own measurements, and then predicted time courses of activation and autophosphorylation of CaMKII under a variety of conditions. Simulations showed that species of calmodulin with fewer than four bound Ca^(2+) play a significant role in activation of CaMKII in the physiological regime, supporting the notion that processing ofCa^(2+) signals in a spine involves competition among target enzymes for binding to unsaturated species of CaM in an environment in which the concentration of Ca^(2+) is fluctuating rapidly. Indeed, we showed that dependence of activation on the frequency of Ca^(2+) transients arises from the kinetics of interaction of fluctuating Ca^(2+) with calmodulin/CaMKII complexes. We used parameter sensitivity analysis to identify which parameters will be most beneficial to measure more carefully to improve the accuracy of predictions. This model provides a quantitative base from which to build more complex dynamic models of postsynaptic signal transduction during learning

Site-Directed Mutagenesis and Site-Specific Binding Analysis of Calmodulin (CaM)

Calcium signaling is a major regulatory system in cells and a crucial part of cell biology. An important element in the decoding of intracellular calcium concentration into downstream processes is the ubiquitous and highly conserved calcium binding protein calmodulin (CaM) which can bind to and modulate the function of hundreds of different target proteins, regulating such processes as synaptic plasticity, gene expression and electrical signaling. The biophysical characterization of binding affinity and cooperative interactions between each of calmodulin’s four EF-hand calcium binding sites is essential for understanding calcium signaling. Highly conserved amino acid sequence differences in the ion binding loops of the EF-hands give each site unique affinity for calcium. EF-hands are almost always found in pairs, where binding to one of the sites affects the affinity of the paired site. We have used spectroscopy to measure site-specific binding in each of the paired binding sites in the CaM N-lobe, along with site-directed mutagenesis, to study the contributions of individual amino acids to the ion binding affinity in the mutated site (cis effects) and in the neighboring site (trans effects). Of the twelve amino acids in the binding loops, five are different between Site 1 and Site 2. We constructed proteins with substituted individual residues from Site 1 to Site 2. CaM with the full Site 1 sequence in both Site 1 and Site 2 shows significant changes in affinity and binding characteristics in both sites. To investigate the contributions of the individual amino acid differences, we made intermediate mutants containing individual amino acid changes in Site 2. The cis-effects of the intermediate mutations on the mutated site, Site 2, seem to be independent and additive, whereas the trans-effects on the non-mutated Site 1 showed unexpected dependence on combinations of amino acid changes in Site 2.Neuroscienc

Evolutionary Covariant Positions within Calmodulin EF-hand Sequences Promote Ligand Binding

Intracellular calcium signaling is an essential regulatory mechanism through calcium-mediated signal transduction pathways involved in many cell processes, such as exocytosis, motility, apoptosis, excitability, transcription, and muscle contraction. The calcium-binding, ubiquitous, and highly conserved protein calmodulin (CaM) is an important regulator of hundreds of target proteins involved in cellular calcium signaling. CaM comprises of two pairs of EF-hand calcium-binding domains and these structural regions of the protein are highly conserved. Studying the molecular mechanisms underlying the binding of calcium to the EF-hands of CaM is critical in understanding the calcium-mediated cellular processes and how improper binding of calcium can lead to various human pathologies. Previous site-specific binding measurements indicate that each of the four EF-hands of CaM have distinct affinities for calcium. In this study, we have utilized covariance patterns and site-specific mutagenesis to analyze calcium affinity in the two EF-hands of the N-lobe of CaM in order to determine the specific amino acids that are evolutionarily conserved to coordinate calcium. The specific amino acids in CaM that we studied are theorized to coevolve, which means that in their protein coding genes, when a mutation occurs, a compensatory mutation is likely to follow to conserve structure and function of CaM. Since CaM is a highly conserved protein with a known structure, covariance analyses will help in understanding which amino acid contacts are most important for the coordination of calcium in the EF-hands of CaM and to determine which amino acids are under evolutionary constraint. Covariance algorithms, multiple sequence analyses and accompanied protein structure analyses were used to identify the two high scoring amino acid pairs in the N-lobe EF-hands: positions 22 and 24 in EF-hand site 1 and positions 58 and 60 in EF-hand site 2. The amino acids in these locations were mutated and accompanied calcium binding was measured to better understand the effects of the mutations on calcium binding. We have found that both the D24N mutation in site 1 and the D58N mutation in site 2 disrupt binding likely due to the removal of a necessary aspartate in the binding site. However, the combined D58N and N60D mutations restore binding in site 2 by providing the necessary aspartate in the covariant location. The N60D mutation by itself has little impact on calcium binding in site 2. Therefore, it is evident that evolution conserves at least one aspartate in the covariant positions of the binding site and the presence of two aspartates in the covariant positions of the binding site has little affect on calcium binding. We are currently studying the covariant positions in site 1 and future work includes structurally analyzing the covariant positions in the C-lobe of CaM and studying covariance patterns of other calcium-binding proteins with EF-hand binding domains.Biochemistr

Purification and characterization of a calmodulin-dependent protein kinase that is highly concentrated in brain

A calcium and calmodulin-dependent protein kinase has been purified from rat brain. It was monitored during the purification by its ability to phosphorylate the synaptic vesicle-associated protein, synapsin I. A 300-fold purification was sufficient to produce kinase that is 90-95% pure as determined by scans of stained sodium dodecyl sulfate-polyacrylamide gels and has a specific activity of 2.9 mumol of 32P transferred per min/mg of protein. Thus, the kinase is a relatively abundant brain enzyme, perhaps comprising as much as 0.3% of the total brain protein. The Stokes radius (95 A) and sedimentation coefficient (16.4 S) of the kinase indicate a holoenzyme molecular weight of approximately 650,000. The holoenzyme is composed of three subunits as judged by their co-migration with kinase activity during the purification steps and co-precipitation with kinase activity by a specific anti-kinase monoclonal antibody. The three subunits have molecular weights of 50,000, 58,000, and 60,000, and have been termed alpha, beta', and beta, respectively. The alpha- and beta-subunits are distinct peptides, however, beta' may have been generated from beta by proteolysis. All three of these subunits bind calmodulin in the presence of calcium and are autophosphorylated under conditions in which the kinase is active. The subunits are present in a ratio of about 3 alpha-subunits to 1 beta/beta'-subunit. We therefore postulate that the 650,000-Da holoenzyme consists of approximately 9 alpha-subunits and 3 beta/beta'-subunits. The abundance of this calmodulin-dependent protein kinase indicates that its activation is likely to be an important biochemical response to increases in calcium ion concentration in neuronal tissue

Structural Analysis and Stochastic Modelling Suggest a Mechanism for Calmodulin Trapping by CaMKII

Activation of CaMKII by calmodulin and the subsequent maintenance of constitutive activity through autophosphorylation at threonine residue 286 (Thr286) are thought to play a major role in synaptic plasticity. One of the effects of autophosphorylation at Thr286 is to increase the apparent affinity of CaMKII for calmodulin, a phenomenon known as “calmodulin trapping”. It has previously been suggested that two binding sites for calmodulin exist on CaMKII, with high and low affinities, respectively. We built structural models of calmodulin bound to both of these sites. Molecular dynamics simulation showed that while binding of calmodulin to the supposed low-affinity binding site on CaMKII is compatible with closing (and hence, inactivation) of the kinase, and could even favour it, binding to the high-affinity site is not. Stochastic simulations of a biochemical model showed that the existence of two such binding sites, one of them accessible only in the active, open conformation, would be sufficient to explain calmodulin trapping by CaMKII. We can explain the effect of CaMKII autophosphorylation at Thr286 on calmodulin trapping: It stabilises the active state and therefore makes the high-affinity binding site accessible. Crucially, a model with only one binding site where calmodulin binding and CaMKII inactivation are strictly mutually exclusive cannot reproduce calmodulin trapping. One of the predictions of our study is that calmodulin binding in itself is not sufficient for CaMKII activation, although high-affinity binding of calmodulin is

Modulation of the asymmetry of sea urchin sperm flagellar bending by calmodulin

Sea urchin spermatozoa demembranated with Triton X-100 in the presence of EGTA, termed potentially asymmetric, generate asymmetric bending waves in reactivation solutions containing EGTA. After they are converted to the potentially symmetric condition by extraction with Triton and millimolar Ca++, they generate symmetric bending waves in reactivation solutions containing EGTA. In the presence of EGTA, their asymmetry can be restored by addition of brain calmodulin or the concentrated supernatant obtained from extraction with Triton and millimolar Ca++. These extracts contain calmodulin, as assayed by gel electrophoresis, radioimmunoassay, activation of brain phosphodiesterase, and Ca++-dependent binding of asymmetry-restoring activity to a trifluorophenothiazine-affinity resin. Conversion to the potentially symmetric condition can also be achieved with trifluoperazine substituted for Triton during the exposure to millimolar Ca++, which suggests that the calmodulin-binding activity of Triton is important for this conversion. These observations suggest that the conversion to the potentially symmetric condition is the result of removal of some of the axonemal calmodulin and provide additional evidence for axonemal calmodulin as a mediator of the effect of Ca++ on the asymmetry of flagellar bending

Phosphatidylinositol (4,5)-bisphosphate turnover by INP51 regulates the cell wall integrity pathway in "Saccharomyces cerevisiae"

Signal transduction pathways are important for the cell to transduce external or internal stimuli where second messengers play an important role as mediators of the stimuli. One important group of second messengers are the phosphoinositide family present in organisms ranging from yeast to mammals. The dephosphorylation and phosphorylation cycle of the phosphatidylinositol species are thought to be important in signaling for recruitment or activation of proteins involved in vesicular transport and/or to control the organization of the actin cytoskeleton. In mammals, phosphatidylinositol (4,5)bisphosphate (PI(4,5)P2) signaling is essential and regulated by various kinases and phosphatases. In the model organism Saccharomyces cerevisiae PI(4,5)P2 signaling is also essential but the regulation remains unclear. My dissertation focuses on the regulation of PI(4,5)P2 signaling in Saccharomyces cerevisiae. The organization of the actin cytoskeleton in Saccharomyces cerevisiae is regulated by different proteins such as calmodulin, CMD1, and here I present data that CMD1 plays a role in the regulation of the only phosphatidylinositol 4-phosphate 5-kinase, MSS4, in Saccharomyces cerevisiae. CMD1 regulates MSS4 activity through an unknown mechanism and thereby controls the organization of the actin cytoskeleton. MSS4 and CMD1 do not physically interact but MSS4 seems to be part of a large molecular weight complex as shown by gel filtration chromatography. This complex could contain regulators of the MSS4 activity. The complex is not caused by dimerization of MSS4 since MSS4 does not interact with itself. Two pathways, the cell wall integrity pathway and TORC2 (target of rapamycin complex 2) signaling cascade are important for the organization of the actin cytoskeleton. Loss of TOR2 function results in a growth defect that can be suppressed by MSS4 overexpression. To further characterize the link between MSS4 and the TORC2 signaling pathway and the cell wall integrity pathway we looked for targets of PI(4,5)P2. The TORC2 pathway and the cell wall integrity pathway signal to the GEF ROM2, an activator of the small GTPase RHO1. In our study we identified ROM2 as a target of PI(4,5)P2 signaling. We observed that the ROM2 localization changes in an mss4 conditional mutant. This suggests that the proper localization needs PI(4,5)P2. This could be mediated by the putative PI(4,5)P2 binding pleckstrin homology (PH) domain of ROM2. To better understand the regulation of PI(4,5)P2 levels in Saccharomyces cerevisiae we focused on one of the PI(4,5)P2 5-phosphatases, INP51. Here we present evidence that INP51 is a new negative regulator of the cell wall integrity pathway as well as the TORC2 pathway. INP51 probably regulates these two pathways by the turnover of PI(4,5)P2 thereby inactivating the effector/s. The deletion of INP51 does not result in any phenotype, but when combined with mutations of the cell wall integrity pathway we observe synthetic interaction. INP51 together with the GTPase activating protein (GAP) SAC7, responsible for the negative regulation of RHO1, negatively regulates the cell wall integrity pathway during vegetative growth. One of the targets of cell wall integrity pathway, the cell wall component chitin, which is normally deposited at the bud end, bud neck and forms bud scars, is delocalized in the mother cell in the sac7 inp51 double deletion mutant. In addition, another downstream component of the cell wall integrity pathway, the MAP kinase MPK1, has increased phosphorylation and protein level in the sac7 inp51 double deletion mutant. This suggests that INP51 is important for the negative regulation of the cell wall integrity pathway. Furthermore, we show evidence that INP51 forms a complex with TAX4 or IRS4, with two EH-domain containing proteins, that positively regulates the activity of INP51 and in this manner negatively regulate the cell wall integrity pathway. The EH-domain is known to bind the NPF-motif. This motif is present in INP51 and is important for INP51 interaction with TAX4 or IRS4. The EH-NPF interaction is a conserved mechanism to build up protein networks. The interaction between an EH-domain containing protein and a PI(4,5)P2 5-phosphatase is conserved. This is demonstrated by the epidermal growth factor substrate EPS15 (EH) interaction with the PI(4,5)P2 5-phosphatase synaptojanin the mammalian orthologue of the Saccharomyces cerevisiae INP proteins. In summary, INP51 together with TAX4 and IRS4, forms complexes important for regulation of PI(4,5)P2 levels. The complexes are linked to the TORC2 signaling pathway and the cell wall integrity pathway, specifically regulating MPK1 activation and chitin biosynthesis. The work presented in this dissertation facilitates the development of a model of the complex regulation of PI(4,5)P2 signaling in Saccharomyces cerevisiae