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

Intracellular acidity impedes KCa3.1 activation by Riluzole and SKA-31

Background

The unique microenvironment in tumors inhibits the normal functioning of tumor-infiltrating lymphocytes, leading to immune evasion and cancer progression. Over-activation of KCa3.1 using positive modulators has been proposed to rescue the anti-tumor response. One of the key characteristics of the tumor microenvironment is extracellular acidity. Herein, we analyzed how intra- and extracellular pH affects K+ currents through KCa3.1 and if the potency of two of its positive modulators, Riluzole and SKA-31, is pH sensitive.

Methods

Whole-cell patch-clamp was used to measure KCa3.1 currents either in activated human peripheral lymphocytes or in CHO cells transiently transfected with either the H192A mutant or wild-type hKCa3.1 in combination with T79D-Calmodulin, or with KCa2.2.

Results

We found that changes in the intra- and extracellular pH minimally influenced the KCa3.1-mediated K+ current. Extracellular pH, in the range of 6.0-8.0, does not interfere with the capacity of Riluzole and SKA-31 to robustly activate the K+ currents through KCa3.1. Contrariwise, an acidic intracellular solution causes a slow, but irreversible loss of potency of both the activators. Using different protocols of perfusion and depolarization we demonstrated that the loss of potency is strictly time and pH-dependent and that this peculiar effect can be observed with a structurally similar channel KCa2.2. While two different point mutations of both KCa3.1 (H192A) and its associated protein Calmodulin (T79D) do not limit the effect of acidity, increasing the cytosolic Ca2+ concentration to saturating levels eliminated the loss-of-potency phenotype.

Conclusion

Based on our data we conclude that KCa3.1 currents are not sensitive the either the intracellular or the extracellular pH in the physiological and pathophysiological range. However, intracellular acidosis in T cells residing in the tumor microenvironment could hinder the potentiating effect of KCa3.1 positive modulators administered to boost their activity. Further research is warranted both to clarify the molecular interactions between the modulators and KCa3.1 at different intracellular pH conditions and to define whether this loss of potency can be observed in cancer models as well.

Human disease-associated calmodulin mutations alter calcineurin function through multiple mechanisms

Calmodulin (CaM) is a ubiquitous, calcium-sensing protein that regulates a multitude of processes throughout the body. In response to changes in [Ca2+], CaM modifies, activates, and deactivates enzymes and ion channels, as well as many other cellular processes. The importance of CaM is highlighted by the conservation of an identical amino acid sequence in all mammals. Alterations to CaM amino acid sequence were once thought to be incompatible with life. During the last decade modifications to the CaM protein sequence have been observed in patients suffering from life-threatening heart disease (calmodulinopathy). Thus far, inadequate or untimely interaction between mutant CaM and several proteins (LTCC, RyR2, and CaMKII) have been identified as mechanisms underlying calmodulinopathy. Given the extensive number of CaM interactions in the body, there are likely many consequences for altering CaM protein sequence. Here, we demonstrate that disease-associated CaM mutations alter the sensitivity and activity of the Ca2+-CaM-enhanced serine/threonine phosphatase calcineurin (CaN). Biophysical characterization by circular dichroism, solution NMR spectroscopy, stopped-flow kinetic measurements, and MD simulations provide mechanistic insight into mutation dysfunction as well as highlight important aspects of CaM Ca2+ signal transduction. We find that individual CaM point mutations (N53I, F89L, D129G, and F141L) impair CaN function, however, the mechanisms are not the same. Specifically, individual point mutations can influence or modify the following properties: CaM binding, Ca2+ binding, and/or Ca2+kinetics. Moreover, structural aspects of the CaNCaM complex can be altered in manners that indicate changes to allosteric transmission of CaM binding to the enzyme active site. Given that loss of CaN function can be fatal, as well as evidence that CaN modifies ion channels already associated with calmodulinopathy, our results raise the possibility that altered CaN function contributes to calmodulinopathy.

KCa2 and KCa3.1 Channels in the Airways: A New Therapeutic Target

K⁺ channels are involved in many critical functions in lung physiology. Recently, the family of Ca²⁺-activated K⁺ channels (K_Ca) has received more attention, and a massive amount of effort has been devoted to developing selective medications targeting these channels. Within the family of K_Ca channels, three small-conductance Ca²⁺-activated K⁺ (K_Ca2) channel subtypes, together with the intermediate-conductance K_Ca3.1 channel, are voltage-independent K⁺ channels, and they mediate Ca²⁺-induced membrane hyperpolarization. Many K_Ca2 channel members are involved in crucial roles in physiological and pathological systems throughout the body. In this article, different subtypes of K_Ca2 and K_Ca3.1 channels and their functions in respiratory diseases are discussed. Additionally, the pharmacology of the K_Ca2 and K_Ca3.1 channels and the link between these channels and respiratory ciliary regulations will be explained in more detail. In the future, specific modulators for small or intermediate Ca²⁺-activated K⁺ channels may offer a unique therapeutic opportunity to treat muco-obstructive lung diseases.

Calmodulin variants associated with congenital arrhythmia impair selectivity for ryanodine receptors

Among its many molecular targets, the ubiquitous calcium sensor protein calmodulin (CaM) recognizes and regulates the activity of ryanodine receptors type 1 (RyR1) and 2 (RyR2), mainly expressed in skeletal and cardiac muscle, respectively. Such regulation is essential to achieve controlled contraction of muscle cells. To unravel the molecular mechanisms underlying the target recognition process, we conducted a comprehensive biophysical investigation of the interaction between two calmodulin variants associated with congenital arrhythmia, namely N97I and Q135P, and a highly conserved calmodulin-binding region in RyR1 and RyR2. The structural, thermodynamic, and kinetic properties of protein-peptide interactions were assessed together with an in-depth structural and topological investigation based on molecular dynamics simulations. This integrated approach allowed us to identify amino acids that are crucial in mediating allosteric processes, which enable high selectivity in molecular target recognition. Our results suggest that the ability of calmodulin to discriminate between RyR1 an RyR2 targets depends on kinetic discrimination and robust allosteric communication between Ca²⁺-binding sites (EF1-EF3 and EF3-EF4 pairs), which is perturbed in both N97I and Q135P arrhythmia-associated variants.

The sequestration mechanism as a generalizable approach to improve the sensitivity of biosensors and bioassays

Biosensors and bioassays, both of which employ proteins and nucleic acids to detect specific molecular targets, have seen significant applications in both biomedical research and clinical practice. This success is largely due to the extraordinary versatility, affinity, and specificity of biomolecular recognition. Nevertheless, these receptors suffer from an inherent limitation: single, saturable binding sites exhibit a hyperbolic relationship (the "Langmuir isotherm") between target concentration and receptor occupancy, which in turn limits the sensitivity of these technologies to small variations in target concentration. To overcome this and generate more responsive biosensors and bioassays, here we have used the sequestration mechanism to improve the steepness of the input/output curves of several bioanalytical methods. As our test bed for this we employed sensors and assays against neutrophil gelatinase-associated lipocalin (NGAL), a kidney biomarker for which enhanced sensitivity will improve the monitoring of kidney injury. Specifically, by introducing sequestration we have improved the responsiveness of an electrochemical aptamer based (EAB) biosensor, and two bioassays, a paper-based "dipstick" assay and an enzyme-linked immunosorbent assay (ELISA). Doing so we have narrowed the dynamic range of these sensors and assays several-fold, thus enhancing their ability to measure small changes in target concentration. Given that introducing sequestration requires only the addition of the appropriate concentration of a high-affinity "depletant," the mechanism appears simple and easily adaptable to tuning the binding properties of the receptors employed in a wide range of biosensors and bioassays.

Alterations in calmodulin-cardiac ryanodine receptor molecular recognition in congenital arrhythmias

Calmodulin (CaM), a ubiquitous and highly conserved Ca²⁺-sensor protein involved in the regulation of over 300 molecular targets, has been recently associated with severe forms of lethal arrhythmia. Here, we investigated how arrhythmia-associated mutations in CaM localized at the C-terminal lobe alter the molecular recognition with Ryanodine receptor 2 (RyR2), specifically expressed in cardiomyocytes. We performed an extensive structural, thermodynamic, and kinetic characterization of the variants D95V/H in the EF3 Ca²⁺-binding motif and of the D129V and D131H/E variants in the EF4 motif, and probed their interaction with RyR2. Our results show that the specific structural changes observed for individual CaM variants do not extend to the complex with the RyR2 target. Indeed, some common alterations emerge at the protein–protein interaction level, suggesting the existence of general features shared by the arrhythmia-associated variants. All mutants showed a faster rate of dissociation from the target peptide than wild-type CaM. Integration of spectroscopic data with exhaustive molecular dynamics simulations suggests that, in the presence of Ca²⁺, functional recognition involves allosteric interactions initiated by the N-terminal lobe of CaM, which shows a lower affinity for Ca²⁺ compared to the C-terminal lobe in the isolated protein.

Simultaneous detection of reciprocal interactions between calmodulin, Ca2+ and molecular targets: a focus on the calmodulin-RyR2 complex

In a recent issue of Biochemical Journal, Brohus et al. (Biochem. J.476, 193-209) investigated the interaction between the ubiquitous intracellular Ca2+-sensor calmodulin (CaM) and peptides that mimic different structural regions of the cardiac ryanodine receptor (RyR2) at different Ca2+ concentrations. For the purpose, a novel bidimensional titration assay based on changes in fluorescence anisotropy was designed. The study identified the CaM domains that selectively bind to a specific CaM-binding domain in RyR2 and demonstrated that the interaction occurs essentially under Ca2+-saturating conditions. This study provides an elegant and experimentally accessible framework for detailed molecular investigations of the emerging life-threatening arrhythmia diseases associated with mutations in the genes encoding CaM. Furthermore, by allowing the measurement of the equilibrium dissociation constant in a protein-protein complex as a function of [Ca2+], the methodology presented by Brohus et al. may have broad applicability to the study of Ca2+ signalling.

Met125 is essential for maintaining the structural integrity of calmodulin's C-terminal domain

We have used NMR and circular dichroism spectroscopy to investigate the structural and dynamic effects of oxidation on calmodulin (CaM), using peroxide and the Met to Gln oximimetic mutations. CaM is a Ca²⁺-sensitive regulatory protein that interacts with numerous targets. Due to its high methionine content, CaM is highly susceptible to oxidation by reactive oxygen species under conditions of cell stress and age-related muscle degeneration. CaM oxidation alters regulation of a host of CaM's protein targets, emphasizing the importance of understanding the mechanism of CaM oxidation in muscle degeneration and overall physiology. It has been shown that the M125Q CaM mutant can mimic the functional effects of methionine oxidation on CaM's regulation of the calcium release channel, ryanodine receptor (RyR). We report here that the M125Q mutation causes a localized unfolding of the C-terminal lobe of CaM, preventing the formation of a hydrophobic cluster of residues near the EF-hand Ca²⁺ binding sites. NMR analysis of CaM oxidation by peroxide offers further insights into the susceptibility of CaM's Met residues to oxidation and the resulting structural effects. These results further resolve oxidation-driven structural perturbation of CaM, with implications for RyR regulation and the decay of muscle function in aging.

Design Principle for Decoding Calcium Signals to Generate Specific Gene Expression Via Transcription

The second messenger calcium plays a key role in conveying specificity of signaling pathways in plant cells. Specific calcium signatures are decoded to generate correct gene expression responses and amplification of calcium signatures is vital to this process. (1) It is not known if this amplification is an intrinsic property of all calcium-regulated gene expression responses and whether all calcium signatures have the potential to be amplified, or (2) how a given calcium signature maintains specificity in cells containing a great number of transcription factors (TFs) and other proteins with the potential to be calcium-regulated. The work presented here uncovers the design principle by which it is possible to decode calcium signals into specific changes in gene transcription in plant cells. Regarding the first question, we found that the binding mechanism between protein components possesses an intrinsic property that will nonlinearly amplify any calcium signal. This nonlinear amplification allows plant cells to effectively distinguish the kinetics of different calcium signatures to produce specific and appropriate changes in gene expression. Regarding the second question, we found that the large number of calmodulin (CaM)-binding TFs or proteins in plant cells form a buffering system such that the concentration of an active CaM-binding TF is insensitive to the concentration of any other CaM-binding protein, thus maintaining specificity. The design principle revealed by this work can be used to explain how any CaM-binding TF decodes calcium signals to generate specific gene expression responses in plant cells via transcription.

Calmodulin Mutations Associated with Heart Arrhythmia: A Status Report

Calmodulin (CaM) is a ubiquitous intracellular Ca²⁺ sensing protein that modifies gating of numerous ion channels. CaM has an extraordinarily high level of evolutionary conservation, which led to the fundamental assumption that mutation would be lethal. However, in 2012, complete exome sequencing of infants suffering from recurrent cardiac arrest revealed de novo mutations in the three human CALM genes. The correlation between mutations and pathophysiology suggests defects in CaM-dependent ion channel functions. Here, we review the current state of the field for all reported CaM mutations associated with cardiac arrhythmias, including knowledge of their biochemical and structural characteristics, and progress towards understanding how these mutations affect cardiac ion channel function.

Computational Modeling Reveals Frequency Modulation of Calcium-cAMP/PKA Pathway in Dendritic Spines

Dendritic spines are the primary excitatory postsynaptic sites that act as subcompartments of signaling. Ca²⁺ is often the first and most rapid signal in spines. Downstream of calcium, the cyclic adenosine monophosphate (cAMP)/protein kinase A (PKA) pathway plays a critical role in the regulation of spine formation, morphological modifications, and ultimately, learning and memory. Although the dynamics of calcium are reasonably well-studied, calcium-induced cAMP/PKA dynamics, particularly with respect to frequency modulation, are not fully explored. In this study, we present a well-mixed model for the dynamics of calcium-induced cAMP/PKA dynamics in dendritic spines. The model is constrained using experimental observations in the literature. Further, we measured the calcium oscillation frequency in dendritic spines of cultured hippocampal CA1 neurons and used these dynamics as model inputs. Our model predicts that the various steps in this pathway act as frequency modulators for calcium, and the high frequency of calcium input is filtered by adenylyl cyclase 1 and phosphodiesterases in this pathway such that cAMP/PKA only responds to lower frequencies. This prediction has important implications for noise filtering and long-timescale signal transduction in dendritic spines. A companion manuscript presents a three-dimensional spatial model for the same pathway.

An Information-Theoretical Approach for Calcium Signaling Specificity

Calcium is a key signaling agent in animals and plants. Its involvement in the regulation of a wide range of processes has led to the question of how calcium signals can activate stimulus-specific responses. We introduce a computational framework for studying intracellular calcium signaling using elements of information theory. We use mutual information to quantify the differential activation of proteins in response to different calcium signals to provide an operational definition of specificity. Using optimization procedures this framework allows us to explore the biochemical determinants of calcium decoding. We explore simple toy models and general binding kinetics approaches to demonstrate the utility and limitations of the proposed framework. Unravelling signaling specificity is key for understanding information processing within cells and for the future design of synthetic nanodevices for molecular communications.

Calmodulinopathy: Functional Effects of CALM Mutations and Their Relationship With Clinical Phenotypes

In spite of the widespread role of calmodulin (CaM) in cellular signaling, CaM mutations lead specifically to cardiac manifestations, characterized by remarkable electrical instability and a high incidence of sudden death at young age. Penetrance of the mutations is surprisingly high, thus postulating a high degree of functional dominance. According to the clinical patterns, arrhythmogenesis in CaM mutations can be attributed, in the majority of cases, to either prolonged repolarization (as in long-QT syndrome, LQTS phenotype), or to instability of the intracellular Ca2+ store (as in catecholamine-induced tachycardias, CPVT phenotype). This review discusses how mutations affect CaM signaling function and how this may relate to the distinct arrhythmia phenotypes/mechanisms observed in patients; this involves mechanistic interpretation of negative dominance and mutation-specific CaM-target interactions. Knowledge of the mechanisms involved may allow critical approach to clinical manifestations and aid in the development of therapeutic strategies for "calmodulinopathies," a recently identified nosological entity.

Calmodulinopathy: A Novel, Life-Threatening Clinical Entity Affecting the Young

Sudden cardiac death (SCD) in the young may often be the first manifestation of a genetic arrythmogenic disease that had remained undiagnosed. Despite the significant discoveries of the genetic bases of inherited arrhythmia syndromes, there remains a measurable fraction of cases where in-depth clinical and genetic investigations fail to identify the underlying SCD etiology. A few years ago, 2 cases of infants with recurrent cardiac arrest episodes, due to what appeared to be as a severe form of long QT syndrome (LQTS), came to our attention. These prompted a number of clinical and genetic research investigations that allowed us to identify a novel, closely associated to LQTS but nevertheless distinct, clinical entity that is now known as calmodulinopathy. Calmodulinopathy is a life-threatening arrhythmia syndrome, affecting mostly young individuals, caused by mutations in any of the 3 genes encoding calmodulin (CaM). Calmodulin is a ubiquitously expressed Ca²⁺ signaling protein that, in the heart, modulates several ion channels and participates in a plethora of cellular processes. We will hereby provide an overview of CaM's structure and function under normal and disease states, highlighting the genetic etiology of calmodulinopathy and the related disease mechanisms. We will also discuss the phenotypic spectrum of patients with calmodulinopathy and present state-of-the art approaches with patient-derived induced pluripotent stem cells that have been thus far adopted in order to accurately model calmodulinopathy in vitro, decipher disease mechanisms and identify novel therapies.

Structural Basis for the Recognition of Eukaryotic Elongation Factor 2 Kinase by Calmodulin

Binding of Ca(2+)-loaded calmodulin (CaM) activates eukaryotic elongation factor 2 kinase (eEF-2K) that phosphorylates eEF-2, its only known cellular target, leading to a decrease in global protein synthesis. Here, using an eEF-2K-derived peptide (eEF-2KCBD) that encodes the region necessary for its CaM-mediated activation, we provide a structural basis for their interaction. The striking feature of this association is the absence of Ca(2+) from the CaM C-lobe sites, even under high Ca(2+) conditions. eEF-2KCBD engages CaM largely through the C lobe of the latter in an anti-parallel 1-5-8 hydrophobic mode reinforced by a pair of unique electrostatic contacts. Sparse interactions of eEF-2KCBD with the CaM N lobe results in persisting inter-lobe mobility. A conserved eEF-2K residue (W85) anchors it to CaM by inserting into a deep hydrophobic cavity within the CaM C lobe. Mutation of this residue (W85S) substantially weakens interactions between full-length eEF-2K and CaM in vitro and reduces eEF-2 phosphorylation in cells.

Kinetic regulation of multi-ligand binding proteins

BACKGROUND

Second messengers, such as calcium, regulate the activity of multisite binding proteins in a concentration-dependent manner. For example, calcium binding has been shown to induce conformational transitions in the calcium-dependent protein calmodulin, under steady state conditions. However, intracellular concentrations of these second messengers are often subject to rapid change. The mechanisms underlying dynamic ligand-dependent regulation of multisite proteins require further elucidation.

RESULTS

In this study, a computational analysis of multisite protein kinetics in response to rapid changes in ligand concentrations is presented. Two major physiological scenarios are investigated: i) Ligand concentration is abundant and the ligand-multisite protein binding does not affect free ligand concentration, ii) Ligand concentration is of the same order of magnitude as the interacting multisite protein concentration and does not change. Therefore, buffering effects significantly influence the amounts of free ligands. For each of these scenarios the influence of the number of binding sites, the temporal effects on intermediate apo- and fully saturated conformations and the multisite regulatory effects on target proteins are investigated.

CONCLUSIONS

The developed models allow for a novel and accurate interpretation of concentration and pressure jump-dependent kinetic experiments. The presented model makes predictions for the temporal distribution of multisite protein conformations in complex with variable numbers of ligands. Furthermore, it derives the characteristic time and the dynamics for the kinetic responses elicited by a ligand concentration change as a function of ligand concentration and the number of ligand binding sites. Effector proteins regulated by multisite ligand binding are shown to depend on ligand concentration in a highly nonlinear fashion.

Engineering PQQ-glucose dehydrogenase into an allosteric electrochemical Ca2+ sensor

Electrochemical biosensors convert biological events to an electrical current.

Calmidazolium evokes high calcium fluctuations in Plasmodium falciparum

Calcium and calmodulin (CaM) are important players in eukaryote cell signaling. In the present study, by using a knockin approach, we demonstrated the expression and localization of CaM in all erythrocytic stages of Plasmodium falciparum. Under extracellular Ca(2+)-free conditions, calmidazolium (CZ), a potent CaM inhibitor, promoted a transient cytosolic calcium ([Ca(2+)]cyt) increase in isolated trophozoites, indicating that CZ mobilizes intracellular sources of calcium. In the same extracellular Ca(2+)-free conditions, the [Ca(2+)]cyt rise elicited by CZ treatment was ~3.5 fold higher when the endoplasmic reticulum (ER) calcium store was previously depleted ruling out the mobilization of calcium from the ER by CZ. The effects of the Ca(2+)/H(+) ionophore ionomycin (ION) and the Na(+)/H(+) ionophore monensin (MON) suggest that the [Ca(2+)]cyt-increasing effect of CZ is driven by the removal of Ca(2+) from at least one Ca(2+)-CaM-related (CaMR) protein as well as by the mobilization of Ca(2+) from intracellular acidic calcium stores. Moreover, we showed that the mitochondrion participates in the sequestration of the cytosolic Ca(2+) elicited by CZ. Finally, the modulation of membrane Ca(2+) channels by CZ and thapsigargin (THG) was demonstrated. The opened channels were blocked by the unspecific calcium channel blocker Co(2+) but not by 2-APB (capacitative calcium entry inhibitor) or nifedipine (L-type Ca(2+) channel inhibitor). Taken together, the results suggested that one CaMR protein is an important modulator of calcium signaling and homeostasis during the Plasmodium intraerythrocytic cell cycle, working as a relevant intracellular Ca(2+) reservoir in the parasite.

Combining modelling and experimental approaches to explain how calcium signatures are decoded by calmodulin-binding transcription activators (CAMTAs) to produce specific gene expression responses

Experimental data show that Arabidopsis thaliana is able to decode different calcium signatures to produce specific gene expression responses. It is also known that calmodulin-binding transcription activators (CAMTAs) have calmodulin (CaM)-binding domains. Therefore, the gene expression responses regulated by CAMTAs respond to calcium signals. However, little is known about how different calcium signatures are decoded by CAMTAs to produce specific gene expression responses. A dynamic model of Ca(2+) -CaM-CAMTA binding and gene expression responses is developed following thermodynamic and kinetic principles. The model is parameterized using experimental data. Then it is used to analyse how different calcium signatures are decoded by CAMTAs to produce specific gene expression responses. Modelling analysis reveals that: calcium signals in the form of cytosolic calcium concentration elevations are nonlinearly amplified by binding of Ca(2+) , CaM and CAMTAs; amplification of Ca(2+) signals enables calcium signatures to be decoded to give specific CAMTA-regulated gene expression responses; gene expression responses to a calcium signature depend upon its history and accumulate all the information during the lifetime of the calcium signature. Information flow from calcium signatures to CAMTA-regulated gene expression responses has been established by combining experimental data with mathematical modelling.

The Actin Nucleator Cobl Is Controlled by Calcium and Calmodulin

Actin nucleation triggers the formation of new actin filaments and has the power to shape cells but requires tight control in order to bring about proper morphologies. The regulation of the members of the novel class of WASP Homology 2 (WH2) domain-based actin nucleators, however, thus far has largely remained elusive. Our study reveals signal cascades and mechanisms regulating Cordon-Bleu (Cobl). Cobl plays some, albeit not fully understood, role in early arborization of neurons and nucleates actin by a mechanism that requires a combination of all three of its actin monomer–binding WH2 domains. Our experiments reveal that Cobl is regulated by Ca²⁺ and multiple, direct associations of the Ca²⁺ sensor Calmodulin (CaM). Overexpression analyses and rescue experiments of Cobl loss-of-function phenotypes with Cobl mutants in primary neurons and in tissue slices demonstrated the importance of CaM binding for Cobl’s functions. Cobl-induced dendritic branch initiation was preceded by Ca²⁺ signals and coincided with local F-actin and CaM accumulations. CaM inhibitor studies showed that Cobl-mediated branching is strictly dependent on CaM activity. Mechanistic studies revealed that Ca²⁺/CaM modulates Cobl’s actin binding properties and furthermore promotes Cobl’s previously identified interactions with the membrane-shaping F-BAR protein syndapin I, which accumulated with Cobl at nascent dendritic protrusion sites. The findings of our study demonstrate a direct regulation of an actin nucleator by Ca²⁺/CaM and reveal that the Ca²⁺/CaM-controlled molecular mechanisms we discovered are crucial for Cobl’s cellular functions. By unveiling the means of Cobl regulation and the mechanisms, by which Ca²⁺/CaM signals directly converge on a cellular effector promoting actin filament formation, our work furthermore sheds light on how local Ca²⁺ signals steer and power branch initiation during early arborization of nerve cells—a key process in neuronal network formation.

The calcium sensor calmodulin directly regulates the actin filament-promoting factor Cobl to help shape the complex architecture of neurons underlying neuronal network formation.

The organization and the formation of new actin filaments by polymerization of actin monomers has the power to shape cells. The rate-limiting step in actin polymerization is “nucleation”—a process during which the first actin monomers are assembled with the help of actin nucleators. This nucleation step requires tight temporal and spatial control in order to achieve proper cell morphologies. Here, we analyse signaling cascades and mechanisms regulating the actin nucleator Cobl, which is crucial for the formation of dendritic arbors of nerve cells—a key process in neuronal network formation. We show that the calcium (Ca²⁺)-binding signaling component calmodulin (CaM) binds to Cobl and regulates its functions. Using 3-D time-lapse analyses of developing neurons, we visualized how Cobl works. We observed local accumulation of CaM, Cobl, actin, and syndapin I—a membrane-shaping protein—at dendritic branch initiation sites. We find that Ca²⁺/CaM modulates Cobl’s actin-binding properties and promotes its interactions with syndapin I, which then serves as a membrane anchor for Cobl. In summary, we i) show a direct regulation of the actin nucleator Cobl by Ca²⁺/CaM, ii) demonstrate that the molecular mechanisms we discovered are crucial for shaping nerve cells, and iii) underscore how local Ca²⁺ signals steer and power branch initiation during early arborization of neurons.

Structures and energies of the transition between two conformations of the alternate frame folding calbindin-D9k protein: a theoretical study

We carried out molecular dynamics simulations and energy calculations for the two states of the alternate frame folding (AFF) calbindin-D_9k protein and their conformational transition in Ca²⁺-free form to address their dynamical transition mechanism.

Ultrasensitive response motifs: basic amplifiers in molecular signalling networks

Multi-component signal transduction pathways and gene regulatory circuits underpin integrated cellular responses to perturbations. A recurring set of network motifs serve as the basic building blocks of these molecular signalling networks. This review focuses on ultrasensitive response motifs (URMs) that amplify small percentage changes in the input signal into larger percentage changes in the output response. URMs generally possess a sigmoid input–output relationship that is steeper than the Michaelis–Menten type of response and is often approximated by the Hill function. Six types of URMs can be commonly found in intracellular molecular networks and each has a distinct kinetic mechanism for signal amplification. These URMs are: (i) positive cooperative binding, (ii) homo-multimerization, (iii) multistep signalling, (iv) molecular titration, (v) zero-order covalent modification cycle and (vi) positive feedback. Multiple URMs can be combined to generate highly switch-like responses. Serving as basic signal amplifiers, these URMs are essential for molecular circuits to produce complex nonlinear dynamics, including multistability, robust adaptation and oscillation. These dynamic properties are in turn responsible for higher-level cellular behaviours, such as cell fate determination, homeostasis and biological rhythm.

Surface dynamics in allosteric regulation of protein-protein interactions: modulation of calmodulin functions by Ca2+

Knowledge of the structural basis of protein-protein interactions (PPI) is of fundamental importance for understanding the organization and functioning of biological networks and advancing the design of therapeutics which target PPI. Allosteric modulators play an important role in regulating such interactions by binding at site(s) orthogonal to the complex interface and altering the protein's propensity for complex formation. In this work, we apply an approach recently developed by us for analyzing protein surfaces based on steered molecular dynamics simulation (SMD) to the study of the dynamic properties of functionally distinct conformations of a model protein, calmodulin (CaM), whose ability to interact with target proteins is regulated by the presence of the allosteric modulator Ca(2+). Calmodulin is a regulatory protein that acts as an intracellular Ca(2+) sensor to control a wide variety of cellular processes. We demonstrate that SMD analysis is capable of pinpointing CaM surfaces implicated in the recognition of both the allosteric modulator Ca(2+) and target proteins. Our analysis of changes in the dynamic properties of the CaM backbone elicited by Ca(2+) binding yielded new insights into the molecular mechanism of allosteric regulation of CaM-target interactions.

Genetically encoded Ca(2+) indicators: properties and evaluation

Genetically encoded calcium ion (Ca(2+)) indicators have become very useful and widely used tools for Ca(2+) imaging, not only in cellular models, but also in living organisms. However, the in vivo and in situ characterization of these indicators is tedious and time consuming, and it does not provide information regarding the suitability of an indicator for particular experimental environments. Thus, initial in vitro evaluation of these tools is typically performed to determine their properties. In this review, we examined the properties of dynamic range, affinity, selectivity, and kinetics for Ca(2+) indicators. Commonly used strategies for evaluating these properties are presented. This article is part of a Special Issue entitled: 12th European Symposium on Calcium.

Analysis and elimination of a bias in targeted molecular dynamics simulations of conformational transitions: application to calmodulin

The popular targeted molecular dynamics (TMD) method for generating transition paths in complex biomolecular systems is revisited. In a typical TMD transition path, the large-scale changes occur early and the small-scale changes tend to occur later. As a result, the order of events in the computed paths depends on the direction in which the simulations are performed. To identify the origin of this bias, and to propose a method in which the bias is absent, variants of TMD in the restraint formulation are introduced and applied to the complex open ↔ closed transition in the protein calmodulin. Due to the global best-fit rotation that is typically part of the TMD method, the simulated system is guided implicitly along the lowest-frequency normal modes, until the large spatial scales associated with these modes are near the target conformation. The remaining portion of the transition is described progressively by higher-frequency modes, which correspond to smaller-scale rearrangements. A straightforward modification of TMD that avoids the global best-fit rotation is the locally restrained TMD (LRTMD) method, in which the biasing potential is constructed from a number of TMD potentials, each acting on a small connected portion of the protein sequence. With a uniform distribution of these elements, transition paths that lack the length-scale bias are obtained. Trajectories generated by steered MD in dihedral angle space (DSMD), a method that avoids best-fit rotations altogether, also lack the length-scale bias. To examine the importance of the paths generated by TMD, LRTMD, and DSMD in the actual transition, we use the finite-temperature string method to compute the free energy profile associated with a transition tube around a path generated by each algorithm. The free energy barriers associated with the paths are comparable, suggesting that transitions can occur along each route with similar probabilities. This result indicates that a broad ensemble of paths needs to be calculated to obtain a full description of conformational changes in biomolecules. The breadth of the contributing ensemble suggests that energetic barriers for conformational transitions in proteins are offset by entropic contributions that arise from a large number of possible paths.

Computational modelling elucidates the mechanism of ciliary regulation in health and disease

Background

Ciliary dysfunction leads to a number of human pathologies, including primary ciliary dyskinesia, nephronophthisis, situs inversus pathology or infertility. The mechanism of cilia beating regulation is complex and despite extensive experimental characterization remains poorly understood. We develop a detailed systems model for calcium, membrane potential and cyclic nucleotide-dependent ciliary motility regulation.

Results

The model describes the intimate relationship between calcium and potassium ionic concentrations inside and outside of cilia with membrane voltage and, for the first time, describes a novel type of ciliary excitability which plays the major role in ciliary movement regulation. Our model describes a mechanism that allows ciliary excitation to be robust over a wide physiological range of extracellular ionic concentrations. The model predicts the existence of several dynamic modes of ciliary regulation, such as the generation of intraciliary Ca²⁺ spike with amplitude proportional to the degree of membrane depolarization, the ability to maintain stable oscillations, monostable multivibrator regimes, all of which are initiated by variability in ionic concentrations that translate into altered membrane voltage.

Conclusions

Computational investigation of the model offers several new insights into the underlying molecular mechanisms of ciliary pathologies. According to our analysis, the reported dynamic regulatory modes can be a physiological reaction to alterations in the extracellular environment. However, modification of the dynamic modes, as a result of genetic mutations or environmental conditions, can cause a life threatening pathology.

Quantitative measurement of Ca(2+)-dependent calmodulin-target binding by Fura-2 and CFP and YFP FRET imaging in living cells

Calcium dynamics and its linked molecular interactions cause a variety of biological responses; thus, exploiting techniques for detecting both concurrently is essential. Here we describe a method for measuring the cytosolic Ca(2+) concentration ([Ca(2+)](i)) and protein-protein interactions within the same cell, using Fura-2 and superenhanced cyan and yellow fluorescence protein (seCFP and seYFP, respectively) FRET imaging techniques. Concentration-independent corrections for bleed-through of Fura-2 into FRET cubes across different time points and [Ca(2+)](i) values allowed for an effective separation of Fura-2 cross-talk signals and seCFP and seYFP cross-talk signals, permitting calculation of [Ca(2+)](i) and FRET with high fidelity. This correction approach was particularly effective at lower [Ca(2+)](i) levels, eliminating bleed-through signals that resulted in an artificial enhancement of FRET. By adopting this correction approach combined with stepwise [Ca(2+)](i) increases produced in living cells, we successfully elucidated steady-state relationships between [Ca(2+)](i) and FRET derived from the interaction of seCFP-tagged calmodulin (CaM) and the seYFP-fused CaM binding domain of myosin light chain kinase. The [Ca(2+)](i) versus FRET relationship for voltage-gated sodium, calcium, and TRPC6 channel CaM binding domains (IQ domain or CBD) revealed distinct sensitivities for [Ca(2+)](i). Moreover, the CaM binding strength at basal or subbasal [Ca(2+)](i) levels provided evidence of CaM tethering or apoCaM binding in living cells. Of the ion channel studies, apoCaM binding was weakest for the TRPC6 channel, suggesting that more global Ca(2+) and CaM changes rather than the local CaM-channel interface domain may be involved in Ca(2+)CaM-mediated regulation of this channel. This simultaneous Fura-2 and CFP- and YFP-based FRET imaging system will thus serve as a simple but powerful means of quantitatively elucidating cellular events associated with Ca(2+)-dependent functions.

Activation of the edema factor of Bacillus anthracis by calmodulin: evidence of an interplay between the EF-calmodulin interaction and calcium binding

Calmodulin (CaM) is a remarkably flexible protein which can bind multiple targets in response to changes in intracellular calcium concentration. It contains four calcium-binding sites, arranged in two globular domains. The calcium affinity of CaM N-terminal domain (N-CaM) is dramatically reduced when the complex with the edema factor (EF) of Bacillus anthracis is formed. Here, an atomic explanation for this reduced affinity is proposed through molecular dynamics simulations and free energy perturbation calculations of the EF-CaM complex starting from different crystallographic models. The simulations show that electrostatic interactions between CaM and EF disfavor the opening of N-CaM domains usually induced by calcium binding. Relative calcium affinities of the N-CaM binding sites are probed by free energy perturbation, and dissociation probabilities are evaluated with locally enhanced sampling simulations. We show that EF impairs calcium binding on N-CaM through a direct conformational restraint on Site 1, by an indirect destabilization of Site 2, and by reducing the cooperativity between the two sites.

Systems biology of cell behavior

Systems Biology approaches to drug discovery largely focus on the increasing understanding of intracellular and cellular circuits, by computational representation of a molecular system followed by parameter validation against experimental data. This chapter outlines a universal approach to systems biology that allows the linking of intracellular molecular machinery and cellular activity. This procedure is achieved by applying mathematical modeling to molecular modules of a cell in the light of systems biology techniques.

Computational modelling suggests dynamic interactions between Ca2+, IP3 and G protein-coupled modules are key to robust Dictyostelium aggregation

Under conditions of starvation, Dictyostelium cells begin a programme of development during which they aggregate to form a multicellular structure by chemotaxis, guided by propagating waves of cyclic AMP that are relayed robustly from cell to cell. In this paper, we develop and analyse a new model for the intracellular and extracellular cAMP dependent processes that regulate Dictyostelium migration. The model allows, for the first time, a quantitative analysis of the dynamic interactions between calcium, IP(3) and G protein-dependent modules that are shown to be key to the generation of robust cAMP oscillations in Dictyostelium cells. The model provides a mechanistic explanation for the transient increase in cytosolic free Ca(2+) concentration seen in recent experiments with the application of the calmodulin inhibitor calmidazolium (R24571) to Dictyostelium cells, and also allows elucidation of the effects of varying both the conductivity of stretch-activated channels and the concentration of external phosphodiesterase on the oscillatory regime of an individual cell. A rigorous analysis of the robustness of the new model shows that interactions between the different modules significantly reduce the sensitivity of the resulting cAMP oscillations to variations in the kinetics of different Dictyostelium cells, an essential requirement for the generation of the spatially and temporally synchronised chemoattractant cAMP waves that guide Dictyostelium aggregation.