Calcium binding to calmodulin and its globular domains

Dynamic interactions of dimeric hub proteins underlie their diverse functions and structures: A comparative analysis of 14-3-3 and LC8

Hub proteins interact with a host of client proteins and regulate multiple cellular functions. Dynamic hubs have a single binding interface for one client at a time resulting in competition among clients with the highest affinity. Dynamic dimeric hubs with two identical sites bind either two different client proteins or two chains of the same client to form homogenous complexes and could also form heterogeneous mixtures of interconverting complexes. Here, we review the interactions of the dimeric hubs 14-3-3 and LC8. 14-3-3 is a phosphoserine/threonine binding protein involved in structuring client proteins and regulating their phosphorylation. LC8 is involved in promoting the dimerization of client peptides and the rigidification of their disordered regions. Both 14-3-3 and LC8 are essential genes, with 14-3-3 playing a crucial role in apoptosis and cell cycle regulation, while LC8 is critical for the assembly of proteins involved in transport, DNA repair, and transcription. Interestingly, both protein dimers can dissociate by phosphorylation, which results in their interactome-wide changes. Their interactions are also regulated by the phosphorylation of their clients. Both form heterogeneous complexes with various functions including phase separation, signaling, and viral hijacking where they restrict the conformational heterogeneity of their dimeric clients that bind nucleic acids. This comparative analysis highlights the importance of dynamic protein-protein interactions in the diversity of functions of 14-3-3 and LC8 and how small differences in structures of interfaces explain why 14-3-3 is primarily involved in the regulation of phosphorylation states while LC8 is primarily involved in the regulation of assembly of large dynamic complexes.

Kinetics insight into the roles of the N- and C-lobes of calmodulin in RyR1 channel regulation

Calmodulin (CaM) activates the skeletal muscle Ca2+ release channel (ryanodine receptor, RyR1) at nanomolar Ca2+ and inhibits it at micromolar Ca2+. CaM conversion from RyR1 activator to inhibitor is due to structural changes induced by Ca2+ binding at CaM's two lobes. However, it remains unclear which lobe provides the switch for this conversion. Here, we attached the environment-sensitive fluorophore acrylodan (Acr) at either lobe of intact CaM or lobe-specific Ca2+-sensitive CaM mutants, and monitored the effects of Ca2+ binding via the fluorescence change of free or RyR1-bound AcrCaM. Using steady state measurements, we found that Ca2+ binding to free CaM causes a dramatic structural change in the N-lobe, but only a slight effect on the C-lobe of the Ca2+-sensitive lobe-specific mutants, in addition to the previously known higher Ca2+ affinity at the C-lobe versus the N-lobe. Using stopped-flow measurements, we found ∼30x faster Ca2+ dissociation from the N- versus C-lobe, and ∼20x slower Ca2+ association to the N-lobe versus C-lobe. These Ca2+ binding properties hold for the CaM/RyR1 complex, and Ca2+ affinity is enhanced at the CaM C-lobe but decreased at the N-lobe by RyR1 binding. We propose that fast Ca2+-binding at the C-lobe of CaM initiates its inhibition to RyR1 at high [Ca2+], while slow Ca2+ binding to the N-lobe is necessary for timely enhancement of the inhibitory effect. The dysregulation of RyR1 by M124Q-CaM may be explained by the lower Ca2+ affinity versus WT-CaM, as suggested by both steady-state and transient kinetics results.

Assessing the predicted impact of single amino acid substitutions in calmodulin for CAGI6 challenges

Recent thermodynamic and functional studies have been conducted to evaluate the impact of amino acid substitutions on Calmodulin (CaM). The Critical Assessment of Genome Interpretation (CAGI) data provider at University of Verona (Italy) measured the melting temperature (Tm) and the percentage of unfolding (%unfold) of a set of CaM variants (CaM challenge dataset). Thermodynamic measurements for the equilibrium unfolding of CaM were obtained by monitoring far-UV Circular Dichroism as a function of temperature. These measurements were used to determine the Tm and the percentage of protein remaining unfolded at the highest temperature. The CaM challenge dataset, comprising a total of 15 single amino acid substitutions, was used to evaluate the effectiveness of computational methods in predicting the Tm and unfolding percentages associated with the variants, and categorizing them as destabilizing or not. For the sixth edition of CAGI, nine independent research groups from four continents (Asia, Australia, Europe, and North America) submitted over 52 sets of predictions, derived from various approaches. In this manuscript, we summarize the results of our assessment to highlight the potential limitations of current algorithms and provide insights into the future development of more accurate prediction tools. By evaluating the thermodynamic stability of CaM variants, this study aims to enhance our understanding of the relationship between amino acid substitutions and protein stability, ultimately contributing to more accurate predictions of the effects of genetic variants.

Fitting and comparison of calcium-calmodulin kinetic schemes to a common data set using non-linear mixed effects modelling

Calmodulin is a calcium binding protein that is essential in calcium signalling in the brain. There are many computational models of calcium-calmodulin binding that capture various calmodulin features. However, existing models have generally been fit to different data sets, with some publications not reporting their training and validation performance. Moreover, there is no model comparison using a common benchmark data set as is common practice in other modeling domains. Finally, some calmodulin models have been fit as a part of a larger kinetic scheme, which may have resulted in parameters being underdetermined. We address these three limitations of previous models by fitting the published calcium-calmodulin schemes to a common calcium-calmodulin data set comprising equilibrium data from Shifman et al. and dynamical data from Faas et al. Due to technical limitations, the amount of uncaged calcium in Faas et al. data could not be predicted with certainty. To find good parameter fits, despite this uncertainty, we used non-linear mixed effects modelling as implemented in the Pumas.jl package. The Akaike information criterion values for our reaction rate constants were significantly lower than for the published parameters, indicating that the published parameters are suboptimal. Moreover, there were significant differences in calmodulin activation, both between the schemes and between our reaction rate and those previously published. A kinetic scheme with independent lobes and unique, rather than identical, binding sites fit the data best. Our results support two hypotheses: (1) partially bound calmodulin is important in cellular signalling; and (2) calcium binding sites within a calmodulin lobe are kinetically distinct rather than identical. We conclude that more attention should be given to validation and comparison of models of individual molecules.

The Molecular Biology of Placental Transport of Calcium to the Human Foetus

From fertilisation to delivery, calcium must be transported into and within the foetoplacental unit for intracellular signalling. This requires very rapid, precisely located Ca2+ transfers. In addition, from around the eighth week of gestation, increasing amounts of calcium must be routed directly from maternal blood to the foetus for bone mineralisation through a flow-through system, which does not impact the intracellular Ca2+ concentration. These different processes are mediated by numerous membrane-sited Ca2+ channels, transporters, and exchangers. Understanding the mechanisms is essential to direct interventions to optimise foetal development and postnatal bone health and to protect the mother and foetus from pre-eclampsia. Ethical issues limit the availability of human foetal tissue for study. Our insight into the processes of placental Ca2+ handling is advancing rapidly, enabled by developing genetic, analytical, and computer technology. Because of their diverse sources, the reports of new findings are scattered. This review aims to pull the data together and to highlight areas of uncertainty. Areas needing clarification include trafficking, membrane expression, and recycling of channels and transporters in the placental microvilli; placental metabolism of vitamin D in gestational diabetes and pre-eclampsia; and the vascular effects of increased endothelial Orai expression by pregnancy-specific beta-1-glycoproteins PSG1 and PSG9.

Welding techniques and manganese concentrations in blood and brain: Results from the WELDFUMES study

This study used whole-brain mapping to investigate the effect of different welding processes on manganese (Mn) accumulation in the brain. Exposure measurements were performed at the welders' workplaces about 3 weeks before a magnetic resonance imaging (MRI) examination. The welders were categorized into three main groups based on welding method, and the T1-relaxation rate (R1) was measured using quantitative MRI (qMRI). Welders using shielded metal arc welding (SMAW) were found to have lower accumulations of total Mn in clusters encompassing white matter, thalamus, putamen, pallidum, and substantia nigra compared with welders using inert gas tungsten arc welding (GTAW) or continuous consumable electrode arc welding (CCEAW). A positive correlation was found between Mn in red blood cells (Mn-RBC) and R1 in a region encompassing pre-and post-central gyri. The results of this study show that the accumulation of free, bound, or compartmentalized Mn ions in the brain differed depending on the welding method used. These differences were predominately located in the basal ganglia but were also found in regions encompassing white matter. The level of Mn-RBC was correlated to the deposition of Mn in the left primary somatosensory and motor cortex and may therefore be linked to neurological and neurobehavioral symptoms.

Role of Calmodulin in Cardiac Disease: Insights on Genotype and Phenotype

Calmodulin, a protein critically important for the regulation of all major cardiac ion channels, is the quintessential cellular calcium sensor and plays a key role in preserving cardiac electrical stability. Its unique importance is highlighted by the presence of 3 genes in 3 different chromosomes encoding for the same protein and by their extreme conservation. Indeed, all 3 calmodulin (CALM) genes are among the most constrained genes in the human genome, that is, the observed variants are much less than expected by chance. Not surprisingly, CALM variants are poorly tolerated and accompany significant clinical phenotypes, of which the most important are those associated with increased risk for life-threatening arrhythmias. Here, we review the current knowledge about calmodulin, its specific physiological, structural, and functional characteristics, and its importance for cardiovascular disease. Given our role in the development of this knowledge, we also share some of our views about currently unanswered questions, including the rational approaches to the clinical management of the affected patients. Specifically, we present some of the most critical information emerging from the International Calmodulinopathy Registry, which we established 10 years ago. Further progress clearly requires deep phenotypic information on as many carriers as possible through international contributions to the registry, in order to expand our knowledge about Calmodulinopathies and guide clinical management.

Divergent Biochemical Properties and Disparate Impact of Arrhythmogenic Calmodulin Mutations on Zebrafish Cardiac Function

Calmodulin (CaM) is a ubiquitous, small cytosolic calcium (Ca2+)-binding sensor that plays a vital role in many cellular processes by binding and regulating the activity of over 300 protein targets. In cardiac muscle, CaM modulates directly or indirectly the activity of several proteins that play a key role in excitation-contraction coupling (ECC), such as ryanodine receptor type 2 (RyR2), l-type Ca2+ (Cav1.2), sodium (NaV1.5) and potassium (KV7.1) channels. Many recent clinical and genetic studies have reported a series of CaM mutations in patients with life-threatening arrhythmogenic syndromes, such as long QT syndrome (LQTS) and catecholaminergic polymorphic ventricular tachycardia (CPVT). We recently showed that four arrhythmogenic CaM mutations (N98I, D132E, D134H, and Q136P) significantly reduce the binding of CaM to RyR2. Herein, we investigate in vivo functional effects of these CaM mutations on the normal zebrafish embryonic heart function by microinjecting complementary RNA corresponding to CaMN98I, CaMD132E, CaMD134H, and CaMQ136P mutants. Expression of CaMD132E and CaMD134H mutants results in significant reduction of the zebrafish heart rate, mimicking a severe form of human bradycardia, whereas expression of CaMQ136P results in an increased heart rate mimicking human ventricular tachycardia. Moreover, analysis of cardiac ventricular rhythm revealed that the CaMD132E and CaMN98I zebrafish groups display an irregular pattern of heart beating and increased amplitude in comparison to the control groups. Furthermore, circular dichroism spectroscopy experiments using recombinant CaM proteins reveals a decreased structural stability of the four mutants compared to the wild-type CaM protein in the presence of Ca2+. Finally, Ca2+-binding studies indicates that all CaM mutations display reduced CaM Ca2+-binding affinities, with CaMD132E exhibiting the most prominent change. Our data suggest that CaM mutations can trigger different arrhythmogenic phenotypes through multiple and complex molecular mechanisms.

Identifying Key Binding Interactions Between the Cardiac L-Type Calcium Channel and Calmodulin Using Molecular Dynamics Simulations

Defects in the binding of the calcium sensing protein calmodulin (CaM) to the L-type calcium channel (CaV1.2) or to the ryanodine receptor type 2 (RyR2) can lead to dangerous cardiac arrhythmias with distinct phenotypes, such as long-QT syndrome (LQTS) and catecholaminergic ventricular tachycardia (CPVT). Certain CaM mutations lead to LQTS while other mutations lead to CPVT, but the mechanisms by which a specific mutation can lead to each disease phenotype are not well-understood. In this study, we use long, 2 μs molecular dynamics simulations and a multitrajectory approach to identify the key binding interactions between the IQ domain of CaV1.2 and CaM. Five key interactions are found between CaV1.2 and CaM in the C-lobe, 1 in the central linker, and 2 in the N-lobe. In addition, while 5 key interactions appear between residues 120-149 in the C-lobe of CaM when it interacts with CaV1.2, only 1 key interaction is found within this region of CaM when it interacts with the RyR2. We show that this difference in the distribution of key interactions correlates with the known distribution of CaM mutations that lead to LQTS or CPVT. This correlation suggests that a disruption of key binding interactions is a plausible mechanism that can lead to these two different disease phenotypes.

Neurogranin modulates the rate of association between calmodulin and target peptides

The best-known mode of action of calmodulin (CaM) is binding of Ca2+ to its N- and C-domains, followed by binding to target proteins. An underappreciated facet of this process is that CaM is typically bound to proteins at basal levels of free Ca2+, including the small, intrinsically disordered, neuronal IQ-motif proteins called PEP-19 and neurogranin (Ng). PEP-19 and Ng would not be effective competitive inhibitors of high-affinity Ca2+-dependent CaM targets at equilibrium because they bind to CaM with relatively low affinity, but they could influence the time course of CaM signaling by affecting the rate of association of CaM with high-affinity Ca2+-dependent targets. This mode of regulation may be domain specific because PEP-19 binds to the C-domain of CaM, whereas Ng binds to both N- and C-domains. In this report, we used a model CaM binding peptide (CKIIp) to characterize the preferred pathway of complex formation with Ca2+-CaM at low levels of free Ca2+ (0.25-1.5 μM), and how PEP-19 and Ng affect this process. We show that the dominant encounter complex involves association of CKIIp with the N-domain of CaM, even though the C-domain has a greater affinity for Ca2+. We also show that Ng greatly decreases the rate of association of Ca2+-CaM with CKIIp due to the relatively slow dissociation of Ng from CaM, and to interactions between the Gly-rich C-terminal region of Ng with the N-domain of CaM, which inhibits formation of the preferred encounter complex with CKIIp. These results provide the general mechanistic paradigms that binding CaM to targets can be driven by its N-domain, and that low-affinity regulators of CaM signaling have the potential to influence the rate of activation of high-affinity CaM targets and potentially affect the distribution of limited CaM among multiple targets during Ca2+ oscillations.

Lanmodulin's EF 2-3 Domain: Insights from Infrared Spectroscopy and Simulations

Lanmodulins are small, ∼110-residue proteins with four EF-hand motifs that demonstrate a picomolar affinity for lanthanide ions, making them efficient in the recovery and separation of these technologically important metals. In this study, we examine the thermodynamic and structural complexities of lanthanide ion binding to a 41-residue domain, EF 2-3, that constitutes the two highest-affinity metal-binding sites in the lanmodulin protein from Methylorubrum extorquens. Using a combination of circular dichroism (CD) spectroscopy, isothermal titration calorimetry (ITC), two-dimensional infrared (2D IR) spectroscopy, and molecular dynamics (MD) simulations, we characterize the metal binding capabilities of EF 2-3. ITC demonstrates that binding occurs between peptide and lanthanides with conditional dissociation constants (Kd) in the range 20-30 μM, with no significant differences in the Kd values for La3+, Eu3+, and Tb3+ at pH 7.4. In addition, CD spectroscopy suggests that only one binding site of EF 2-3 undergoes a significant conformational change in the presence of lanthanides. 2D IR spectroscopy demonstrates the presence of both mono- and bidentate binding configurations in EF 2-3 with all three lanthanides. MD simulations, supported by Eu3+ luminescence measurements, explore these results, suggesting a competition between water-lanthanide and carboxylate-lanthanide interactions in the EF 2-3 domain. These results underscore the role of the core helical bundle of the protein architecture in influencing binding affinities and communication between the metal-binding sites in the full-length protein.

Interplay of α-synuclein with Lipid Membranes: Cooperative Adsorption, Membrane Remodeling and Coaggregation

α-Synuclein is a small neuronal protein enriched at presynaptic termini. It is hypothesized to play a role in neurotransmitter release and synaptic vesicle cycling, while the formation of α-synuclein amyloid fibrils is associated with several neurodegenerative diseases, most notably Parkinson's Disease. The molecular mechanisms of both the physiological and pathological functions of α-synuclein remain to be fully understood, but in both cases, interactions with membranes play an important role. In this Perspective, we discuss several aspects of α-synuclein interactions with lipid membranes including cooperative adsorption, membrane remodeling and α-synuclein amyloid fibril formation in the presence of lipid membranes. We highlight the coupling between the different phenomena and their interplay in the context of physiological and pathological functions of α-synuclein.

Mononuclear binding and catalytic activity of europium(III) and gadolinium(III) at the active site of the model metalloenzyme phosphotriesterase

Lanthanide ions have ideal chemical properties for catalysis, such as hard Lewis acidity, fast ligand-exchange kinetics, high coordination-number preferences and low geometric requirements for coordination. As a result, many small-molecule lanthanide catalysts have been described in the literature. Yet, despite the ability of enzymes to catalyse highly stereoselective reactions under gentle conditions, very few lanthanoenzymes have been investigated. In this work, the mononuclear binding of europium(III) and gadolinium(III) to the active site of a mutant of the model enzyme phosphotriesterase are described using X-ray crystallography at 1.78 and 1.61 Å resolution, respectively. It is also shown that despite coordinating a single non-natural metal cation, the PTE-R18 mutant is still able to maintain esterase activity.

Hydrogen/Deuterium Exchange Mass Spectrometry Provides Insights into the Role of Drosophila Testis-Specific Myosin VI Light Chain AndroCaM

In Drosophila testis, myosin VI plays a special role, distinct from its motor function, by anchoring components to the unusual actin-based structures (cones) that are required for spermatid individualization. For this, the two calmodulin (CaM) light-chain molecules of myosin VI are replaced by androcam (ACaM), a related protein with 67% identity to CaM. Although ACaM has a similar bi-lobed structure to CaM, with two EF hand-type Ca2+ binding sites per lobe, only one functional Ca2+ binding site operates in the amino-terminus. To understand this light chain substitution, we used hydrogen-deuterium exchange mass spectrometry (HDX-MS) to examine dynamic changes in ACaM and CaM upon Ca2+ binding and interaction with the two CaM binding motifs of myosin VI (insert2 and IQ motif). HDX-MS reveals that binding of Ca2+ to ACaM destabilizes its N-lobe but stabilizes the entire C-lobe, whereas for CaM, Ca2+ binding induces a pattern of alternating stabilization/destabilization throughout. The conformation of this stable holo-C-lobe of ACaM seems to be a "prefigured" version of the conformation adopted by the holo-C-lobe of CaM for binding to insert2 and the IQ motif of myosin VI. Strikingly, the interaction of holo-ACaM with either peptide converts the holo-N-lobe to its Ca2+-free, more stable, form. Thus, ACaM in vivo should bind the myosin VI light chain sites in an apo-N-lobe/holo-C-lobe state that cannot fulfill the Ca2+-related functions of holo-CaM required for myosin VI motor assembly and activity. These findings indicate that inhibition of myosin VI motor activity is a precondition for transition to an anchoring function.

Energetic and structural insights behind calcium induced conformational transition in calmodulin

Calmodulin (CaM) is a key signaling protein that triggers several cellular and physiological processes inside the cell. Upon binding with calcium ion, CaM undergoes large scale conformational transition from a closed state to an open state that facilitates its interaction with various target protein and regulates their activity. This work explores the origin of the energetic and structural variation of the wild type and mutated CaM and explores the molecular origin for the structural differences between them. We first calculated the sequential calcium binding energy to CaM using the PDLD/S-LRA/β approach. This study  shows a very good correlation with experimental calcium binding energies. Next we calculated the calcium binding energies to the wild type CaM and several mutated CaM systems which were reported experimentally. On the structural aspect, it has been reported experimentally that certain mutation (Q41L-K75I) in calcium bound CaM leads to complete conformational transition from an open to a closed state. By using equilibrium molecular dynamics simulation, free energy calculation and contact frequency map analysis, we have shown that the formation of a cluster of long-range hydrophobic contacts, initiated by the Q41L-K75I CaM variant is the driving force behind its closing motion. This study unravels the energetics and structural aspects behind calcium ion induced conformational changes in wild type CaM and its variant.

Role of myristoylation in modulating PCaP1 interaction with calmodulin

Plasma membrane-associated Cation-binding Protein 1 (PCaP1) belongs to the plant-unique DREPP protein family with largely unknown biological functions but ascertained roles in plant development and calcium (Ca2+) signaling. PCaP1 is anchored to the plasma membrane via N-myristoylation and a polybasic cluster, and its N-terminal region can bind Ca2+/calmodulin (CaM). However, the molecular determinants of PCaP1-Ca2+-CaM interaction and the functional impact of myristoylation in the complex formation and Ca2+ sensitivity of CaM remained to be elucidated. Herein, we investigated the direct interaction between Arabidopsis PCaP1 (AtPCaP1) and CaM1 (AtCaM1) using both myristoylated and non-myristoylated peptides corresponding to the N-terminal region of AtPCaP1. ITC analysis showed that AtCaM1 forms a high affinity 1:1 complex with AtPCaP1 peptides and the interaction is strictly Ca2+-dependent. Spectroscopic and kinetic Ca2+ binding studies showed that the myristoylated peptide dramatically increased the Ca2+-binding affinity of AtCaM1 and slowed the Ca2+ dissociation rates from both the C- and N-lobes, thus suggesting that the myristoylation modulates the mechanism of AtPCaP1 recognition by AtCaM1. Furthermore, NMR and CD spectroscopy revealed that the structure of both the N- and C-lobes of Ca2+-AtCaM1 changes markedly in the presence of the myristoylated AtPCaP1 peptide, which assumes a helical structure in the final complex. Overall, our results indicate that AtPCaP1 biological function is strictly related to the presence of multiple ligands, i.e., the myristoyl moiety, Ca2+ ions and AtCaM1 and only a full characterization of their equilibria will allow for a complete molecular understanding of the putative role of PCaP1 as signal protein.

Assessments of Variational Autoencoder in Protein Conformation Exploration

Molecular dynamics (MD) simulations have been extensively used to study protein dynamics and subsequently functions. However, MD simulations are often insufficient to explore adequate conformational space for protein functions within reachable timescales. Accordingly, many enhanced sampling methods, including variational autoencoder (VAE) based methods, have been developed to address this issue. The purpose of this study is to evaluate the feasibility of using VAE to assist in the exploration of protein conformational landscapes. Using three modeling systems, we showed that VAE could capture high-level hidden information which distinguishes protein conformations. These models could also be used to generate new physically plausible protein conformations for direct sampling in favorable conformational spaces. We also found that VAE worked better in interpolation than extrapolation and increasing latent space dimension could lead to a trade-off between performances and complexities.

Ionic displacement of Ca2+ by Pb2+ in calmodulin is affected by arrhythmia-associated mutations

Lead is a highly toxic metal that severely perturbs physiological processes even at sub-micromolar levels, often by disrupting the Ca²⁺ signaling pathways. Recently, Pb²⁺-associated cardiac toxicity has emerged, with potential involvement of both the ubiquitous Ca²⁺ sensor protein calmodulin (CaM) and ryanodine receptors. In this work, we explored the hypothesis that Pb²⁺ contributes to the pathological phenotype of CaM variants associated with congenital arrhythmias. We performed a thorough spectroscopic and computational characterization of CaM conformational switches in the co-presence of Pb²⁺ and four missense mutations associated with congenital arrhythmias, namely N53I, N97S, E104A and F141L, and analyzed their effects on the recognition of a target peptide of RyR2. When bound to any of the CaM variants, Pb²⁺ is difficult to displace even under equimolar Ca²⁺ concentrations, thus locking all CaM variants in a specific conformation, which exhibits characteristics of coiled-coil assemblies. All arrhythmia-associated variants appear to be more susceptible to Pb²⁺ than WT CaM, as the conformational transition towards the coiled-coil conformation occurs at lower Pb²⁺, regardless of the presence of Ca²⁺, with altered cooperativity. The presence of arrhythmia-associated mutations specifically alters the cation coordination of CaM variants, in some cases involving allosteric communication between the EF-hands in the two domains. Finally, while wild type CaM increases the affinity for the RyR2 target in the presence of Pb²⁺, no specific pattern could be detected for all other variants, ruling out a synergistic effect of Pb²⁺ and mutations in the recognition process.

Investigation of the Plutonium(IV) Interactions with Two Variants of the EF-Hand Ca-Binding Site I of Calmodulin

Due to its presence in the nuclear industry and its strong radiotoxicity, plutonium is an actinide of major interest in the event of internal contamination. To improve the understanding of its mechanisms of transport and accumulation in the body, the complexation of Pu(IV) to the most common protein calcium-binding motif found in cells, the EF-hand motif of calmodulin, was investigated. Visible and X-ray absorption spectroscopies (XAS) in solution made it possible to investigate the speciation of plutonium at physiological pH (pH 7.4) and pH 6 in two variants of the calmodulin Ca-binding site I and using Pu(IV) in different media: carbonate, chloride, or nitrate solutions. Three different species of Pu were identified in the samples, with formation of 1:1 Pu(IV):calmodulin peptide complexes, Pu(IV) reduction, and formation of peptide-mediated Pu(IV) hexanuclear cluster.

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.

Biochemistry and pathophysiology of the Transient Potential Receptor Vanilloid 6 (TRPV6) calcium channel

TRPV6 is a Transient Receptor Potential Vanilloid (TRPV) cation channel with high selectivity for Ca²⁺ ions. First identified in 1999 in a search for the gene which mediates intestinal Ca²⁺ absorption, its far more extensive repertoire as a guardian of intracellular Ca²⁺ has since become apparent. Studies on TRPV6-deficient mice demonstrated additional important roles in placental Ca²⁺ transport, fetal bone development and male fertility. The first reports of inherited deficiency in newborn babies appeared in 2018, revealing its physiological importance in humans. There is currently strong evidence that TRPV6 also contributes to the pathogenesis of some common cancers. The recently reported association of TRPV6 deficiency with non-alcoholic chronic pancreatitis suggests a role in normal pancreatic function. Over time and with greater awareness of TRPV6, other disease-associations are likely to emerge. Powerful analytical tools have provided invaluable insights into the structure and operation of TRPV6. Its roles in Ca²⁺ signaling and carcinogenesis, and the use of channel inhibitors in cancer treatment are being intensively investigated. This review first briefly describes the biochemistry and physiology of the channel, and analytical methods used to investigate these. The focus subsequently shifts to the clinical disorders associated with abnormal expression and the underlying pathophysiology. The aims of this review are to increase awareness of this channel, and to draw together findings from a wide range of sources which may help to formulate new ideas for further studies.

Calcium dependence of both lobes of calmodulin is involved in binding to a cytoplasmic domain of SK channels

KCa2.1-3 Ca2+-activated K+-channels (SK) require calmodulin to gate in response to cellular Ca2+. A model for SK gating proposes that the N-terminal domain (N-lobe) of calmodulin is required for activation, but an immobile C-terminal domain (C-lobe) has constitutive, Ca2+-independent binding. Although structures support a domain-driven hypothesis of SK gate activation by calmodulin, only a partial understanding is possible without measuring both channel activity and protein binding. We measured SK2 (KCa2.2) activity using inside-out patch recordings. Currents from calmodulin-disrupted SK2 channels can be restored with exogenously applied calmodulin. We find that SK2 activity only approaches full activation with full-length calmodulin with both an N- and a C-lobe. We measured calmodulin binding to a C-terminal SK peptide (SKp) using both composition-gradient multi-angle light-scattering and tryptophan emission spectra. Isolated lobes bind to SKp with high affinity, but isolated lobes do not rescue SK2 activity. Consistent with earlier models, N-lobe binding to SKp is stronger in Ca2+, and C-lobe-binding affinity is strong independent of Ca2+. However, a native tryptophan in SKp is sensitive to Ca2+ binding to both the N- and C-lobes of calmodulin at Ca2+ concentrations that activate SK2, demonstrating that the C-lobe interaction with SKp changes with Ca2+. Our peptide-binding data and electrophysiology show that SK gating models need deeper scrutiny. We suggest that the Ca2+-dependent associations of both lobes of calmodulin to SKp are crucial events during gating. Additional investigations are necessary to complete a mechanistic gating model consistent with binding, physiology, and structure.

Binding and Functional Folding (BFF): A Physiological Framework for Studying Biomolecular Interactions and Allostery

EF-hand Ca2+-binding proteins (CBPs), such as S100 proteins (S100s) and calmodulin (CaM), are signaling proteins that undergo conformational changes upon increasing intracellular Ca2+. Upon binding Ca2+, S100 proteins and CaM interact with protein targets and induce important biological responses. The Ca2+-binding affinity of CaM and most S100s in the absence of target is weak (CaKD > 1 μM). However, upon effector protein binding, the Ca2+ affinity of these proteins increases via heterotropic allostery (CaKD < 1 μM). Because of the high number and micromolar concentrations of EF-hand CBPs in a cell, at any given time, allostery is required physiologically, allowing for (i) proper Ca2+ homeostasis and (ii) strict maintenance of Ca2+-signaling within a narrow dynamic range of free Ca2+ ion concentrations, [Ca2+]free. In this review, mechanisms of allostery are coalesced into an empirical "binding and functional folding (BFF)" physiological framework. At the molecular level, folding (F), binding and folding (BF), and BFF events include all atoms in the biomolecular complex under study. The BFF framework is introduced with two straightforward BFF types for proteins (type 1, concerted; type 2, stepwise) and considers how homologous and nonhomologous amino acid residues of CBPs and their effector protein(s) evolved to provide allosteric tightening of Ca2+ and simultaneously determine how specific and relatively promiscuous CBP-target complexes form as both are needed for proper cellular function.

Inter-Site Cooperativity of Calmodulin N-Terminal Domain and Phosphorylation Synergistically Improve the Affinity and Selectivity for Uranyl

Uranyl-protein interactions participate in uranyl trafficking or toxicity to cells. In addition to their qualitative identification, thermodynamic data are needed to predict predominant mechanisms that they mediate in vivo. We previously showed that uranyl can substitute calcium at the canonical EF-hand binding motif of calmodulin (CaM) site I. Here, we investigate thermodynamic properties of uranyl interaction with site II and with the whole CaM N-terminal domain by spectrofluorimetry and ITC. Site II has an affinity for uranyl about 10 times lower than site I. Uranyl binding at site I is exothermic with a large enthalpic contribution, while for site II, the enthalpic contribution to the Gibbs free energy of binding is about 10 times lower than the entropic term. For the N-terminal domain, macroscopic binding constants for uranyl are two to three orders of magnitude higher than for calcium. A positive cooperative process driven by entropy increases the second uranyl-binding event as compared with the first one, with ΔΔG = -2.0 ± 0.4 kJ mol-1, vs. ΔΔG = -6.1 ± 0.1 kJ mol-1 for calcium. Site I phosphorylation largely increases both site I and site II affinity for uranyl and uranyl-binding cooperativity. Combining site I phosphorylation and site II Thr7Trp mutation leads to picomolar dissociation constants Kd1 = 1.7 ± 0.3 pM and Kd2 = 196 ± 21 pM at pH 7. A structural model obtained by MD simulations suggests a structural role of site I phosphorylation in the affinity modulation.

Markov state modelling reveals heterogeneous drug-inhibition mechanism of Calmodulin

Calmodulin (CaM) is a calcium sensor which binds and regulates a wide range of target-proteins. This implicitly enables the concentration of calcium to influence many downstream physiological responses, including muscle contraction, learning and depression. The antipsychotic drug trifluoperazine (TFP) is a known CaM inhibitor. By binding to various sites, TFP prevents CaM from associating to target-proteins. However, the molecular and state-dependent mechanisms behind CaM inhibition by drugs such as TFP are largely unknown. Here, we build a Markov state model (MSM) from adaptively sampled molecular dynamics simulations and reveal the structural and dynamical features behind the inhibitory mechanism of TFP-binding to the C-terminal domain of CaM. We specifically identify three major TFP binding-modes from the MSM macrostates, and distinguish their effect on CaM conformation by using a systematic analysis protocol based on biophysical descriptors and tools from machine learning. The results show that depending on the binding orientation, TFP effectively stabilizes features of the calcium-unbound CaM, either affecting the CaM hydrophobic binding pocket, the calcium binding sites or the secondary structure content in the bound domain. The conclusions drawn from this work may in the future serve to formulate a complete model of pharmacological modulation of CaM, which furthers our understanding of how these drugs affect signaling pathways as well as associated diseases.

Calmodulin (CaM) is a calcium-sensing protein which makes other proteins dependent on the surrounding calcium concentration by binding to these proteins. Such protein-protein interactions with CaM are vital for calcium to control many physiological pathways within the cell. The antipsychotic drug trifluoperazine (TFP) inhibits CaM’s ability to bind and regulate other proteins. Here, we use molecular dynamics simulations together with Markov state modeling and machine learning to understand the structural and dynamical features by which TFP bound to the one domain of CaM prevents association to other proteins. We find that TFP encourages CaM to adopt a conformation that is like the one stabilized in absence of calcium: depending on the binding orientation of TFP, the drug indeed either affects the CaM hydrophobic binding pocket, the calcium binding sites or the secondary structure content in the domain. Understanding TFP binding is a first step towards designing better drugs targeting CaM.

Calmodulin in Paramecium: Focus on Genomic Data

Calcium (Ca²⁺) is a universal second messenger that plays a key role in cellular signaling. However, Ca²⁺ signals are transduced with the help of Ca²⁺-binding proteins, which serve as sensors, transducers, and elicitors. Among the collection of these Ca²⁺-binding proteins, calmodulin (CaM) emerged as the prototypical model in eukaryotic cells. This is a small protein that binds four Ca²⁺ ions and whose functions are multiple, controlling many essential aspects of cell physiology. CaM is universally distributed in eukaryotes, from multicellular organisms, such as human and land plants, to unicellular microorganisms, such as yeasts and ciliates. Here, we review most of the information gathered on CaM in Paramecium, a group of ciliates. We condense the information here by mentioning that mature Paramecium CaM is a 148 amino acid-long protein codified by a single gene, as in other eukaryotic microorganisms. In these ciliates, the protein is notoriously localized and regulates cilia function and can stimulate the activity of some enzymes. When Paramecium CaM is mutated, cells show flawed locomotion and/or exocytosis. We further widen this and additional information in the text, focusing on genomic data.

A Modeling and Analysis Study Reveals That CaMKII in Synaptic Plasticity Is a Dominant Affecter in CaM Systems in a T286 Phosphorylation-Dependent Manner

NMDAR-dependent synaptic plasticity in the hippocampus consists of two opposing forces: long-term potentiation (LTP), which strengthens synapses and long-term depression (LTD), which weakens synapses. LTP and LTD are associated with memory formation and loss, respectively. Synaptic plasticity is controlled at a molecular level by Ca2+-mediated protein signaling. Here, Ca2+ binds the protein, calmodulin (CaM), which modulates synaptic plasticity in both directions. This is because Ca2+-bound CaM activates both LTD-and LTP-inducing proteins. Understanding how CaM responds to Ca2+ signaling and how this translates into synaptic plasticity is therefore important to understanding synaptic plasticity induction. In this paper, CaM activation by Ca2+ and calmodulin binding to downstream proteins was mathematically modeled using differential equations. Simulations were monitored with and without theoretical knockouts and, global sensitivity analyses were performed to determine how Ca2+/CaM signaling occurred at various Ca2+ signals when CaM levels were limiting. At elevated stimulations, the total CaM pool rapidly bound to its protein binding targets which regulate both LTP and LTD. This was followed by CaM becoming redistributed from low-affinity to high-affinity binding targets. Specifically, CaM was redistributed away from LTD-inducing proteins to bind the high-affinity LTP-inducing protein, calmodulin-dependent kinase II (CaMKII). In this way, CaMKII acted as a dominant affecter and repressed activation of opposing CaM-binding protein targets. The model thereby showed a novel form of CaM signaling by which the two opposing pathways crosstalk indirectly. The model also found that CaMKII can repress cAMP production by repressing CaM-regulated proteins, which catalyze cAMP production. The model also found that at low Ca2+ stimulation levels, typical of LTD induction, CaM signaling was unstable and is therefore unlikely to alone be enough to induce synaptic depression. Overall, this paper demonstrates how limiting levels of CaM may be a fundamental aspect of Ca2+ regulated signaling which allows crosstalk among proteins without requiring directly interaction.

Interdomain Dynamics via Paramagnetic NMR on the Highly Flexible Complex Calmodulin/Munc13-1

Paramagnetic NMR constraints are very useful to study protein interdomain motion, but their interpretation is not always straightforward. On the example of the particularly flexible complex Calmodulin/Munc13-1, we present a new approach to characterize this motion with pseudocontact shifts and residual dipolar couplings. Using molecular mechanics, we sampled the conformational space of the complex and used a genetic algorithm to find ensembles that are in agreement with the data. We used the Bayesian information criterion to determine the ideal ensemble size. This way, we were able to make an accurate, unambiguous, reproducible model of the interdomain motion of Calmodulin/Munc13-1 without prior knowledge about the domain orientation from crystallography.

Time-resolved solid state NMR of biomolecular processes with millisecond time resolution

We review recent efforts to develop and apply an experimental approach to the structural characterization of transient intermediate states in biomolecular processes that involve large changes in molecular conformation or assembly state. This approach depends on solid state nuclear magnetic resonance (ssNMR) measurements that are performed at very low temperatures, typically 25-30 K, with signal enhancements from dynamic nuclear polarization (DNP). This approach also involves novel technology for initiating the process of interest, either by rapid mixing of two solutions or by a rapid inverse temperature jump, and for rapid freezing to trap intermediate states. Initiation by rapid mixing or an inverse temperature jump can be accomplished in approximately-one millisecond. Freezing can be accomplished in approximately 100 microseconds. Thus, millisecond time resolution can be achieved. Recent applications to the process by which the biologically essential calcium sensor protein calmodulin forms a complex with one of its target proteins and the process by which the bee venom peptide melittin converts from an unstructured monomeric state to a helical, tetrameric state after a rapid change in pH or temperature are described briefly. Future applications of millisecond time-resolved ssNMR are also discussed briefly.

Ca2+ signalling system initiated by endoplasmic reticulum stress stimulates PERK activation

The accumulation of unfolded proteins within the Endoplasmic Reticulum (ER) activates a signal transduction pathway termed the unfolded protein response (UPR), which attempts to restore ER homoeostasis. If this cannot be done, UPR signalling ultimately induces apoptosis. Ca²⁺ depletion in the ER is a potent inducer of ER stress. Despite the ubiquity of Ca²⁺ as an intracellular messenger, the precise mechanism(s) by which Ca²⁺ release affects the UPR remains unknown. Tethering a genetically encoded Ca²⁺ indicator (GCamP6) to the ER membrane revealed novel Ca²⁺ signalling events initiated by Ca²⁺ microdomains in human astrocytes under ER stress, induced by tunicamycin (Tm), an N-glycosylation inhibitor, as well as in a cell model deficient in all three inositol triphosphate receptor isoforms. Pharmacological and molecular studies indicate that these local events are mediated by translocons and that the Ca²⁺ microdomains impact (PKR)-like-ER kinase (PERK), an UPR sensor, activation. These findings reveal the existence of a Ca²⁺ signal mechanism by which stressor-mediated Ca²⁺ release regulates ER stress.

Role of water-bridged interactions in metal ion coupled protein allostery

Allosteric communication between distant parts of proteins controls many cellular functions, in which metal ions are widely utilized as effectors to trigger the allosteric cascade. Due to the involvement of strong coordination interactions, the energy landscape dictating the metal ion binding is intrinsically rugged. How metal ions achieve fast binding by overcoming the landscape ruggedness and thereby efficiently mediate protein allostery is elusive. By performing molecular dynamics simulations for the Ca²⁺ binding mediated allostery of the calmodulin (CaM) domains, each containing two Ca²⁺ binding helix-loop-helix motifs (EF-hands), we revealed the key role of water-bridged interactions in Ca²⁺ binding and protein allostery. The bridging water molecules between Ca²⁺ and binding residue reduces the ruggedness of ligand exchange landscape by acting as a lubricant, facilitating the Ca²⁺ coupled protein allostery. Calcium-induced rotation of the helices in the EF-hands, with the hydrophobic core serving as the pivot, leads to exposure of hydrophobic sites for target binding. Intriguingly, despite being structurally similar, the response of the two symmetrically arranged EF-hands upon Ca²⁺ binding is asymmetric. Breakage of symmetry is needed for efficient allosteric communication between the EF-hands. The key roles that water molecules play in driving allosteric transitions are likely to be general in other metal ion mediated protein allostery.

Natural proteins often utilize allostery in executing a variety of functions. Metal ions are typical cofactors to trigger the allosteric cascade. In this work, using the Ca²⁺ sensor protein calmodulin as the model system, we revealed crucial roles of water-bridged interactions in the metal ion coupled protein allostery. The coordination of the Ca²⁺ to the binding site involves an intermediate in which the water molecule bridges the Ca²⁺ and the liganding residue. The bridging water reduces the free energy barrier height of ligand exchange, therefore facilitating the ligand exchange and allosteric coupling by acting as a lubricant. We also showed that the response of the two symmetrically arranged EF-hand motifs of CaM domains upon Ca²⁺ binding is asymmetric, which is directly attributed to the differing dehydration process of the Ca²⁺ ions and is needed for efficient allosteric communication.

Calcium- and Integrin-Binding Protein 2 (CIB2) in Physiology and Disease: Bright and Dark Sides

Calcium- and integrin-binding protein 2 (CIB2) is a small EF-hand protein capable of binding Mg²⁺ and Ca²⁺ ions. While its biological function remains largely unclear, an increasing number of studies have shown that CIB2 is an essential component of the mechano-transduction machinery that operates in cochlear hair cells. Mutations in the gene encoding CIB2 have been associated with non-syndromic deafness. In addition to playing an important role in the physiology of hearing, CIB2 has been implicated in a multitude of very different processes, ranging from integrin signaling in platelets and skeletal muscle to autophagy, suggesting extensive functional plasticity. In this review, we summarize the current understanding of biochemical and biophysical properties of CIB2 and the biological roles that have been proposed for the protein in a variety of processes. We also highlight the many molecular aspects that remain unclarified and deserve further investigation.

Chelator-Based Parameterization of the 12-6-4 Lennard-Jones Molecular Mechanics Potential for More Realistic Metal Ion-Protein Interactions

Metal ions are associated with a variety of proteins and play critical roles in a wide range of biochemical processes. There are multiple ways to study and quantify protein–metal ion interactions, including molecular dynamics simulations. Recently, the AMBER molecular mechanics forcefield was modified to include a 12-6-4 Lennard-Jones potential, which allows for a better description of nonbonded terms through the additional pairwise C_ij coefficients. Here, we demonstrate a method of generating C_ij parameters that allows parametrization of specific metal ion-ligating groups in order to tune binding energies computed by thermodynamic integration. The new C_ij coefficients were tested on a series of chelators: ethylenediaminetetraacetic acid, nitrilotriacetic acid, egtazic acid, and the EF1 loop peptides from the proteins lanmodulin and calmodulin. The new parameters show significant improvements in computed binding energies relative to existing force fields and produce coordination numbers and ion-oxygen distances that are in good agreement with experimental values. This parametrization method should be extensible to a range of other systems and could be readily adapted to tune properties other than binding energies.

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.

Rapid kinetics of calcium dissociation from plant calmodulin and calmodulin-like proteins and effect of target peptides

Calcium (Ca²⁺) signaling represents a universal information code in plants, playing crucial roles spanning developmental processes to stress responses. Ca²⁺ signals are decoded into defined plant adaptive responses by different Ca²⁺ sensing proteins, including calmodulin (CaM) and calmodulin-like (CML) proteins. Although major advances have been achieved in describing how these Ca²⁺ decoding proteins interact and regulate downstream target effectors, the molecular details of these processes remain largely unknown. Herein, the kinetics of Ca²⁺ dissociation from a conserved CaM and two CML isoforms from A. thaliana has been studied by fluorescence stopped-flow spectroscopy. Kinetic data were obtained for the isolated Ca²⁺-bound proteins as well as for the proteins complexed with different target peptides. Moreover, the lobe specific interactions between the Ca²⁺ sensing proteins and their targets were characterized by using a panel of protein mutants deficient in Ca²⁺ binding at the N-lobe or C-lobe. Results were analyzed and discussed in the context of the Ca²⁺-decoding and Ca²⁺-controlled target binding mechanisms in plants.

CPVT-associated calmodulin variants N53I and A102V dysregulate Ca2+ signalling via different mechanisms

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is an inherited condition that can cause fatal cardiac arrhythmia. Human mutations in the Ca2+ sensor calmodulin (CaM) have been associated with CPVT susceptibility, suggesting that CaM dysfunction is a key driver of the disease. However, the detailed molecular mechanism remains unclear. Focusing on the interaction with the cardiac ryanodine receptor (RyR2), we determined the effect of CPVT-associated variants N53I and A102V on the structural characteristics of CaM and on Ca2+ fluxes in live cells. We provide novel data showing that interaction of both Ca2+/CaM-N53I and Ca2+/CaM-A102V with the RyR2 binding domain is decreased. Ca2+/CaM-RyR23583-3603 high-resolution crystal structures highlight subtle conformational changes for the N53I variant, with A102V being similar to wild type (WT). We show that co-expression of CaM-N53I or CaM-A102V with RyR2 in HEK293 cells significantly increased the duration of Ca2+ events; CaM-A102V exhibited a lower frequency of Ca2+ oscillations. In addition, we show that CaMKIIδ (also known as CAMK2D) phosphorylation activity is increased for A102V, compared to CaM-WT. This paper provides novel insight into the molecular mechanisms of CPVT-associated CaM variants and will facilitate the development of strategies for future therapies.

Mass Spectrometry-Based Structural Proteomics for Metal Ion/Protein Binding Studies

Metal ions are critical for the biological and physiological functions of many proteins. Mass spectrometry (MS)-based structural proteomics is an ever-growing field that has been adopted to study protein and metal ion interactions. Native MS offers information on metal binding and its stoichiometry. Footprinting approaches coupled with MS, including hydrogen/deuterium exchange (HDX), "fast photochemical oxidation of proteins" (FPOP) and targeted amino-acid labeling, identify binding sites and regions undergoing conformational changes. MS-based titration methods, including "protein-ligand interactions by mass spectrometry, titration and HD exchange" (PLIMSTEX) and "ligand titration, fast photochemical oxidation of proteins and mass spectrometry" (LITPOMS), afford binding stoichiometry, binding affinity, and binding order. These MS-based structural proteomics approaches, their applications to answer questions regarding metal ion protein interactions, their limitations, and recent and potential improvements are discussed here. This review serves as a demonstration of the capabilities of these tools and as an introduction to wider applications to solve other questions.

Harnessing Thorpe-Ingold Dialkylation to Access High-Hill-Percentage pH Probes

Sensitivity is an important parameter for a molecular probe. Hill-type pH probes exhibit improved detection sensitivity compared to the traditional pH probes following the Henderson-Hasselbalch equation. Exploiting positive cooperativity, we recently devised a novel molecular scaffold (PHX) to offer such an unconventional Hill-type pH titration profile. We previously confirmed that PHX is not a pure Hill-type probe yet. Only 64% of its absorbance/fluorescence turn-on is the result of a Hill-type pathway. The remaining 36% is from an undesired Henderson-Hasselbalch-type pathway (HH pathway). In this work, the Thorpe-Ingold dialkylation was harnessed to further suppress the percent contribution of the HH pathway down to 16%. We also propose that PHX is a viable molecular model for assessing the efficacy of the steric compressing effect induced by different Thorpe-Ingold dialkylations.

Analysis of protein-ligand interactions from titrations and nuclear magnetic resonance relaxation dispersions

We present PLIS, a publicly available, open‐source software for the determination of protein–ligand dissociation constants that can be used to characterize biological processes or to shed light on biophysical aspects of interactions. PLIS can analyze data from titration experiments monitored by for instance fluorescence spectroscopy or from nuclear magnetic resonance relaxation dispersion experiments. In addition to analysis of experimental data, PLIS includes functionality for generation of synthetic data, useful for understanding how different parameters effect the data in order to better analyze experiments.

Structural Basis for the Functional Diversity of Centrins: A Focus on Calcium Sensing Properties and Target Recognition

Centrins are a family of small, EF hand-containing proteins that are found in all eukaryotes and are often complexed with centrosome-related structures. Since their discovery, centrins have attracted increasing interest due to their multiple, diverse cellular functions. Centrins are similar to calmodulin (CaM) in size, structure and domain organization, although in contrast to CaM, the majority of centrins possess at least one calcium (Ca²⁺) binding site that is non-functional, thus displaying large variance in Ca²⁺ sensing abilities that could support their functional versatility. In this review, we summarize current knowledge on centrins from both biophysical and structural perspectives with an emphasis on centrin-target interactions. In-depth analysis of the Ca²⁺ sensing properties of centrins and structures of centrins complexed with target proteins can provide useful insight into the mechanisms of the different functions of centrins and how these proteins contribute to the complexity of the Ca²⁺ signaling cascade. Moreover, it can help to better understand the functional redundancy of centrin isoforms and centrin-binding proteins.

Physiologically Relevant Free Ca2+ Ion Concentrations Regulate STRA6-Calmodulin Complex Formation via the BP2 Region of STRA6

The interaction of calmodulin (CaM) with the receptor for retinol uptake, STRA6, involves an α-helix termed BP2 that is located on the intracellular side of this homodimeric transporter (Chen et al., 2016 [1]). In the absence of Ca2+, NMR data showed that a peptide derived from BP2 bound to the C-terminal lobe (C-lobe) of Mg2+-bound CaM (MgCaM). Upon titration of Ca2+ into MgCaM-BP2, NMR chemical shift perturbations (CSPs) were observed for residues in the C-lobe, including those in the EF-hand Ca2+-binding domains, EF3 and EF4 (CaKD = 60 ± 7 nM). As higher concentrations of free Ca2+ were achieved, CSPs occurred for residues in the N-terminal lobe (N-lobe) including those in EF1 and EF2 (CaKD = 1000 ± 160 nM). Thermodynamic and kinetic Ca2+ binding studies showed that BP2 addition increased the Ca2+-binding affinity of CaM and slowed its Ca2+ dissociation rates (koff) in both the C- and N-lobe EF-hand domains, respectively. These data are consistent with BP2 binding to the C-lobe of CaM at low free Ca2+ concentrations (<100 nM) like those found at resting intracellular levels. As free Ca2+ levels approach 1000 nM, which is typical inside a cell upon an intracellular Ca2+-signaling event, BP2 is shown here to interact with both the N- and C-lobes of Ca2+-loaded CaM (CaCaM-BP2). Because this structural rearrangement observed for the CaCaM-BP2 complex occurs as intracellular free Ca2+ concentrations approach those typical of a Ca2+-signaling event (CaKD = 1000 ± 160 nM), this conformational change could be relevant to vitamin A transport by full-length CaCaM-STRA6.

Calmodulinopathy in inherited arrhythmia syndromes

Calmodulin (CaM) is a ubiquitous intracellular calcium sensor that controls and regulates key cellular functions. In all vertebrates, three CaM genes located on separate chromosomes encode an identical 149 amino acid protein, implying an extraordinarily high level of evolutionary importance and suggesting that CaM mutations would be possibly fatal. Inherited arrhythmia syndromes comprise a spectrum of primary electrical disorders caused by mutations in genes encoding ion channels or associated proteins leading to various cardiac arrhythmias, unexplained syncope, and sudden cardiac death. CaM mutations have emerged as an independent entity among inherited arrhythmia syndromes, referred to as calmodulinopathies. The most common clinical presentation associated with calmodulinopathy is congenital long QT syndrome, followed by catecholaminergic polymorphic ventricular tachycardia, both of which significantly increase the possibility of repeated syncope, lethal arrhythmic events, and sudden cardiac death, especially in young individuals. Here, we aim to give an overview of biochemical and structural characteristics of CaM and progress toward updating current known CaM mutations and associated clinical phenotypes. We also review the possible mechanisms underlying calmodulinopathy, based on several key in vitro studies. We expect that further experimental studies are needed to explore the complexity of calmodulinopathy.

High-affinity tamoxifen analogues retain extensive positional disorder when bound to calmodulin

Using a combination of NMR and fluorescence measurements, we have investigated the structure and dynamics of the complexes formed between calcium-loaded calmodulin (CaM) and the potent breast cancer inhibitor idoxifene, a derivative of tamoxifen. High-affinity binding (K d ∼ 300  nM) saturates with a 2 : 1 idoxifene : CaM complex. The complex is an ensemble where each idoxifene molecule is predominantly in the vicinity of one of the two hydrophobic patches of CaM but, in contrast with the lower-affinity antagonists TFP, J-8, and W-7, does not substantially occupy the hydrophobic pocket. At least four idoxifene orientations per domain of CaM are necessary to satisfy the intermolecular nuclear Overhauser effect (NOE) restraints, and this requires that the idoxifene molecules switch rapidly between positions. The CaM molecule is predominantly in the form where the N and C-terminal domains are in close proximity, allowing for the idoxifene molecules to contact both domains simultaneously. Hence, the 2 : 1 idoxifene : CaM complex illustrates how high-affinity binding occurs without the loss of extensive positional dynamics.

Distinctive Properties and Powerful Neuromodulation of Nav1.6 Sodium Channels Regulates Neuronal Excitability

Voltage-gated sodium channels (Navs) are critical determinants of cellular excitability. These ion channels exist as large heteromultimeric structures and their activity is tightly controlled. In neurons, the isoform Nav1.6 is highly enriched at the axon initial segment and nodes, making it critical for the initiation and propagation of neuronal impulses. Changes in Nav1.6 expression and function profoundly impact the input-output properties of neurons in normal and pathological conditions. While mutations in Nav1.6 may cause channel dysfunction, aberrant changes may also be the result of complex modes of regulation, including various protein-protein interactions and post-translational modifications, which can alter membrane excitability and neuronal firing properties. Despite decades of research, the complexities of Nav1.6 modulation in health and disease are still being determined. While some modulatory mechanisms have similar effects on other Nav isoforms, others are isoform-specific. Additionally, considerable progress has been made toward understanding how individual protein interactions and/or modifications affect Nav1.6 function. However, there is still more to be learned about how these different modes of modulation interact. Here, we examine the role of Nav1.6 in neuronal function and provide a thorough review of this channel's complex regulatory mechanisms and how they may contribute to neuromodulation.

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.

The leucine zipper EF-hand containing transmembrane protein-1 EF-hand is a tripartite calcium, temperature, and pH sensor

Leucine Zipper EF-hand containing transmembrane protein-1 (LETM1) is an inner mitochondrial membrane protein that mediates mitochondrial calcium (Ca²⁺ )/proton exchange. The matrix residing carboxyl (C)-terminal domain contains a sequence identifiable EF-hand motif (EF1) that is highly conserved among orthologues. Deletion of EF1 abrogates LETM1 mediated mitochondrial Ca²⁺ flux, highlighting the requirement of EF1 for LETM1 function. To understand the mechanistic role of this EF-hand in LETM1 function, we characterized the biophysical properties of EF1 in isolation. Our data show that EF1 exhibits α-helical secondary structure that is augmented in the presence of Ca²⁺ . Unexpectedly, EF1 features a weak (~mM), but specific, apparent Ca²⁺ -binding affinity, consistent with the canonical Ca²⁺ coordination geometry, suggested by our solution NMR. The low affinity is, at least in part, due to an Asp at position 12 of the binding loop, where mutation to Glu increases the affinity by ~4-fold. Further, the binding affinity is sensitive to pH changes within the physiological range experienced by mitochondria. Remarkably, EF1 unfolds at high and low temperatures. Despite these unique EF-hand properties, Ca²⁺ binding increases the exposure of hydrophobic regions, typical of EF-hands; however, this Ca²⁺ -induced conformational change shifts EF1 from a monomer to higher order oligomers. Finally, we showed that a second, putative EF-hand within LETM1 is unreactive to Ca²⁺ either in isolation or tandem with EF1. Collectively, our data reveal that EF1 is structurally and biophysically responsive to pH, Ca²⁺ and temperature, suggesting a role as a multipartite environmental sensor within LETM1.

The S2-S3 Loop of Kv7.4 Channels Is Essential for Calmodulin Regulation of Channel Activation

Kv7.4 (KCNQ4) voltage-gated potassium channels control excitability in the inner ear and the central auditory pathway. Mutations in Kv7.4 channels result in inherited progressive deafness in humans. Calmodulin (CaM) is crucial for regulating Kv7 channels, but how CaM affects Kv7 activity has remained unclear. Here, based on electrophysiological recordings, we report that the third EF hand (EF3) of CaM controls the calcium-dependent regulation of Kv7.4 activation and that the S2–S3 loop of Kv7.4 is essential for the regulation mediated by CaM. Overexpression of the mutant CaM₁₂₃₄, which loses the calcium binding ability of all four EF hands, facilitates Kv7.4 activation by accelerating activation kinetics and shifting the voltage dependence of activation leftwards. The single mutant CaM₃, which loses the calcium binding ability of the EF3, phenocopies facilitating effects of CaM₁₂₃₄ on Kv7.4 activation. Kv7.4 channels co-expressed with wild-type (WT) CaM show inhibited activation when intracellular calcium levels increase, while Kv7.4 channels co-expressed with CaM₁₂₃₄ or CaM₃ are insensitive to calcium. Mutations C156A, C157A, C158V, R159, and R161A, which are located within the Kv7.4 S2–S3 loop, dramatically facilitate activation of Kv7.4 channels co-expressed with WT CaM but have no effect on activation of Kv7.4 channels co-expressed with CaM₃, indicating that these five mutations decrease the inhibitory effect of Ca²⁺/CaM. The double mutation C156A/R159A decreases Ca²⁺/CaM binding and completely abolishes CaM-mediated calcium-dependent regulation of Kv7.4 activation. Taken together, our results provide mechanistic insights into CaM regulation of Kv7.4 activation and highlight the crucial role of the Kv7.4 S2–S3 loop in CaM regulation.

More Than Just Simple Interaction between STIM and Orai Proteins: CRAC Channel Function Enabled by a Network of Interactions with Regulatory Proteins

The calcium-release-activated calcium (CRAC) channel, activated by the release of Ca2+ from the endoplasmic reticulum (ER), is critical for Ca2+ homeostasis and active signal transduction in a plethora of cell types. Spurred by the long-sought decryption of the molecular nature of the CRAC channel, considerable scientific effort has been devoted to gaining insights into functional and structural mechanisms underlying this signalling cascade. Key players in CRAC channel function are the Stromal interaction molecule 1 (STIM1) and Orai1. STIM1 proteins span through the membrane of the ER, are competent in sensing luminal Ca2+ concentration, and in turn, are responsible for relaying the signal of Ca2+ store-depletion to pore-forming Orai1 proteins in the plasma membrane. A direct interaction of STIM1 and Orai1 allows for the re-entry of Ca2+ from the extracellular space. Although much is already known about the structure, function, and interaction of STIM1 and Orai1, there is growing evidence that CRAC under physiological conditions is dependent on additional proteins to function properly. Several auxiliary proteins have been shown to regulate CRAC channel activity by means of direct interactions with STIM1 and/or Orai1, promoting or hindering Ca2+ influx in a mechanistically diverse manner. Various proteins have also been identified to exert a modulatory role on the CRAC signalling cascade although inherently lacking an affinity for both STIM1 and Orai1. Apart from ubiquitously expressed representatives, a subset of such regulatory mechanisms seems to allow for a cell-type-specific control of CRAC channel function, considering the rather restricted expression patterns of the specific proteins. Given the high functional and clinical relevance of both generic and cell-type-specific interacting networks, the following review shall provide a comprehensive summary of regulators of the multilayered CRAC channel signalling cascade. It also includes proteins expressed in a narrow spectrum of cells and tissues that are often disregarded in other reviews of similar topics.

Millisecond Time-Resolved Solid-State NMR Reveals a Two-Stage Molecular Mechanism for Formation of Complexes between Calmodulin and a Target Peptide from Myosin Light Chain Kinase

Calmodulin (CaM) mediates a wide range of biological responses to changes in intracellular Ca²⁺ concentrations through its calcium-dependent binding affinities to numerous target proteins. Binding of two Ca²⁺ ions to each of the two four-helix-bundle domains of CaM results in major conformational changes that create a potential binding site for the CaM binding domain of a target protein, which also undergoes major conformational changes to form the complex with CaM. Details of the molecular mechanism of complex formation are not well established, despite numerous structural, spectroscopic, thermodynamic, and kinetic studies. Here, we report a study of the process by which the 26-residue peptide M13, which represents the CaM binding domain of skeletal muscle myosin light chain kinase, forms a complex with CaM in the presence of excess Ca²⁺ on the millisecond time scale. Our experiments use a combination of selective ¹³C labeling of CaM and M13, rapid mixing of CaM solutions with M13/Ca²⁺ solutions, rapid freeze-quenching of the mixed solutions, and low-temperature solid state nuclear magnetic resonance (ssNMR) enhanced by dynamic nuclear polarization. From measurements of the dependence of 2D ¹³C-¹³C ssNMR spectra on the time between mixing and freezing, we find that the N-terminal portion of M13 converts from a conformationally disordered state to an α-helix and develops contacts with the C-terminal domain of CaM in about 2 ms. The C-terminal portion of M13 becomes α-helical and develops contacts with the N-terminal domain of CaM more slowly, in about 8 ms. The level of structural order in the CaM/M13/Ca²⁺ complexes, indicated by ¹³C ssNMR line widths, continues to increase beyond 27 ms.

Benefits and constrains of covalency: the role of loop length in protein stability and ligand binding

Protein folding is governed by non-covalent interactions under the benefits and constraints of the covalent linkage of the backbone chain. In the current work we investigate the influence of loop length variation on the free energies of folding and ligand binding in a small globular single-domain protein containing two EF-hand subdomains—calbindin D_9k. We introduce a linker extension between the subdomains and vary its length between 1 to 16 glycine residues. We find a close to linear relationship between the linker length and the free energy of folding of the Ca²⁺-free protein. In contrast, the linker length has only a marginal effect on the Ca²⁺ affinity and cooperativity. The variant with a single-glycine extension displays slightly increased Ca²⁺ affinity, suggesting that the slightly extended linker allows optimized packing of the Ca²⁺-bound state. For the extreme case of disconnected subdomains, Ca²⁺ binding becomes coupled to folding and assembly. Still, a high affinity between the EF-hands causes the non-covalent pair to retain a relatively high apparent Ca²⁺ affinity. Our results imply that loop length variation could be an evolutionary option for modulating properties such as protein stability and turnover without compromising the energetics of the specific function of the protein.