Conserved properties of individual Ca2+-binding sites in calmodulin

Zinc mitigates cadmium-induced sperm dysfunction through regulating Ca2+ and metallothionein expression in the freshwater crab Sinopotamon henanense

Cadmium (Cd) is a highly toxic heavy metal element that might adversely affect sperm function such as the acrosome reaction (AR). Although it is widely recognized that zinc (Zn) plays a crucial role in sperm quality, the complete elucidation of how Zn ameliorates Cd-induced sperm dysfunction is still unclear. In this study, we aimed to explore the protective effects of Zn against the sperm dysfunction induced by Cd in the freshwater crab Sinopotamon henanense. The results demonstrated that Cd exposure not only impaired the sperm ultrastructure, but also caused sperm dysfunction by decreasing the AR induction rate, acrosome enzyme activity, and Ca2+ content in sperm while elevating the activity and transcription expression of key Ca2+ signaling pathway-related proteins Calmodulin (CAM) and Ca2+-ATPase. However, the administration of Zn was found to alleviate Cd-induced sperm morphological and functional disorders by increasing the activity and transcription levels of CaM and Ca2+-ATPase, thereby regulating intracellular Ca2+ homeostasis and reversing the decrease in Ca2+ contents caused by Cd. Furthermore, this study was the first to investigate the distribution of metallothionein (MT) in the AR of S. henanense, and it was found that Zn can reduce the elevated levels of MT in crabs caused by Cd, demonstrating the significance of Zn in inducing MT to participate in the AR process and in metal detoxification in S. henanense. These findings offer novel perspectives and substantiation regarding the utilization of Zn as a protective agent against Cd-induced toxicity and hold significant practical implications for mitigating Cd-induced sperm dysfunction.

Real-time single-molecule imaging of CaMKII-calmodulin interactions

The binding of calcium/calmodulin (CAM) to calcium/calmodulin-dependent protein kinase II (CaMKII) initiates an ATP-driven cascade that triggers CaMKII autophosphorylation. The autophosphorylation in turn increases the CaMKII affinity for CAM. Here, we studied the ATP dependence of CAM association with the actin-binding CaMKIIβ isoform using single-molecule total internal reflection fluorescence microscopy. Rhodamine-CAM associations/dissociations to surface-immobilized Venus-CaMKIIβ were resolved with 0.5 s resolution from video records, batch-processed with a custom algorithm. CAM occupancy was determined simultaneously with spot-photobleaching measurement of CaMKII holoenzyme stoichiometry. We show the ATP-dependent increase of the CAM association requires dimer formation for both the α and β isoforms. The study of mutant β holoenzymes revealed that the ATP-dependent increase in CAM affinity results in two distinct states. The phosphorylation-defective (T287.306-307A) holoenzyme resides only in the low-affinity state. CAM association is further reduced in the T287A holoenzyme relative to T287.306-307A. In the absence of ATP, the affinity of CAM for the T287.306-307A mutant and the wild-type monomer are comparable. The affinity of the ATP-binding impaired (K43R) mutant is even weaker. In ATP, the K43R holoenzyme resides in the low-affinity state. The phosphomimetic mutant (T287D) resides only in a 1000-fold higher-affinity state, with mean CAM occupancy of more than half of the 14-mer holoenzyme stoichiometry in picomolar CAM. ATP promotes T287D holoenzyme disassembly but does not elevate CAM occupancy. Single Poisson distributions characterized the ATP-dependent CAM occupancy of mutant holoenzymes. In contrast, the CAM occupancy of the wild-type population had a two-state distribution with both low- and high-affinity states represented. The low-affinity state was the dominant state, a result different from published in vitro assays. Differences in assay conditions can alter the balance between activating and inhibitory autophosphorylation. Bound ATP could be sufficient for CaMKII structural function, while antagonistic autophosphorylations may tune CaMKII kinase-regulated action-potential frequency decoding in vivo.

Bioseparation of rare earth elements and high value-added biomaterials applications

Rare earth elements (REEs) are a group of critical minerals and extensively employed in new material manufacturing. However, separation of lanthanides is difficult because of their similar chemical natures. Current lanthanide leaching and separation methods require hazardous compounds, resulting in severe environmental concerns. Bioprocessing of lanthanides offers an emerging class of tools for REE separation due to mild leaching conditions and highly selective separation scenarios. In the course of biopreparation, engineered microbes not only dissolve REEs from ores but also allow for selective separation of the lanthanides. In this review, we present an overview of recent advances in microbes and proteins used for the biomanufacturing of lanthanides and discuss high value-added applications of REE-derived biomaterials. We begin by introducing the fundamental interactions between natural microbes and REEs. Then we discuss the rational design of chassis microbes for bioleaching and biosorption. We also highlight the investigations on REE binding proteins and their applications in the synthesis of high value-added biomaterials. Finally, future opportunities and challenges for the development of next generation lanthanide-binding biological systems are discussed.

Structure and Function of Calcium-Activated Chloride Channels and Phospholipid Scramblases in the TMEM16 Family

The transmembrane protein 16 (TMEM16) family consists of Ca2+-activated chloride channels and phospholipid scramblases. Ten mammalian TMEM16 proteins, TMEM16A-K (with no TMEM16I), and several non-mammalian TMEM16 proteins, such as afTMEM16 and nhTMEM16, have been discovered. All known TMEM16 proteins are homodimeric proteins containing two subunits. Each subunit consists of ten transmembrane helices with Ca2+-binding sites and a single ion-permeation/phospholipid transport pathway. The ion-permeation pathway and the phospholipid transport pathway of TMEM16 proteins have a wide intracellular vestibule, a narrow neck, and a smaller extracellular vestibule. Interestingly, the lining wall of the ion-permeation/phospholipid transport pathway may be formed, at least partially, by membrane phospholipids, though the degree of pore-wall forming by phospholipids likely varies among TMEM16 proteins. Thus, the biophysical properties and activation mechanisms of TMEM16 proteins could differ from each other accordingly. Here we review the current understanding of the structure and function of TMEM16 molecules.

Interaction of Calmodulin with TRPM: An Initiator of Channel Modulation

Transient receptor potential melastatin (TRPM) channels, a subfamily of the TRP superfamily, constitute a diverse group of ion channels involved in mediating crucial cellular processes like calcium homeostasis. These channels exhibit complex regulation, and one of the key regulatory mechanisms involves their interaction with calmodulin (CaM), a cytosol ubiquitous calcium-binding protein. The association between TRPM channels and CaM relies on the presence of specific CaM-binding domains in the channel structure. Upon CaM binding, the channel undergoes direct and/or allosteric structural changes and triggers down- or up-stream signaling pathways. According to current knowledge, ion channel members TRPM2, TRPM3, TRPM4, and TRPM6 are directly modulated by CaM, resulting in their activation or inhibition. This review specifically focuses on the interplay between TRPM channels and CaM and summarizes the current known effects of CaM interactions and modulations on TRPM channels in cellular physiology.

Calmodulin, Calcium/Calmodulin-Dependent Kinases-1 and 2 Regulate Expression of the Heat Shock Proteins for Heat Shock Tolerance and Pheromone Signaling Genes for Sexual Development in Neurospora crassa

Calmodulin (CaM) is a primary Ca2+ sensor that binds and activates numerous target proteins and regulates several cellular processes in eukaryotes. CaM is essential in Neurospora crassa; therefore, we generated a CaM mutant using repeat-induced point (RIP) mutation and investigated the cmdRIP mutant phenotypes. We also studied knockout mutants of four Ca2+/CaM kinases (camk-1, 2, 3, and 4) for their role during stress conditions and sexual development. The cmdRIP, ∆camk-1, and ∆camk-2 mutants showed reduced survival and growth rates under heat stress, oxidative stress, pH, and ER stress conditions. In addition, under the heat stress conditions, expression of the heat shock protein genes hsp70 and hsp80 was reduced in the cmdRIP, ∆camk-1, and ∆camk-2 mutants. The cmdRIP mutant was also defective in cell fusion, its vegetative hyphae could not support the fertilized wild type perithecia graft, and female sterile. Furthermore, the expression of pheromone signaling genes pre-1, pre-2, ccg-4, mfa-1, and fmf-1 was reduced in the cmdRIP, ∆camk-1, and ∆camk-2 mutants. Therefore, CaM, Ca2+/CaMK-1 and 2 are involved in the tolerance to heat stress conditions and sexual development by regulating the heat shock and pheromone response pathways, respectively, in N. crassa.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-023-01091-8.

Machine learning-based modulation of Ca2+-binding affinity in EF-hand proteins and comparative structural insights into site-specific cooperative binding

Ca2+-binding proteins are present in almost all living organisms and different types display different levels of binding affinities for the cation. Here, we report two new scoring schemes enabling the user to estimate and manipulate the calcium binding affinities in EF hand containing proteins. To validate this, we designed a unique EF-hand loop capable of binding calcium with high affinity by altering five residues. The N-terminal domain of Entamoeba histolytica calcium-binding protein1 (NtEhCaBP1) is used for site-directed mutagenesis to incorporate the designed loop sequence into the second EF-hand motif of this protein, referred as Nt-EhCaBP1-EF2 mutant. The binding isotherms calculated using ITC calorimetry showed that Nt-EhCaBP1-EF2 mutant site binds Ca2+ with higher affinity than Wt-Nt-EhCaBP1, by ∼600 times. The crystal structure of the mutant displayed more compact Ca2+-coordination spheres in both of its EF loops than the structure of the wildtype protein. The compact coordination sphere of EF-2 causes the bend in the helix-3, which leads to the formation of unexpected hexamer of NtEhCaBP1-EF2 mutant structure. Further dynamic correlation analysis revealed that the mutation in the second EF loop changed the entire residue network of the monomer, resulting in stronger coordination of Ca2+ even in another EF-hand loop.

Long QT syndrome-associated calmodulin variants disrupt the activity of the slowly activating delayed rectifier potassium channel

Calmodulin (CaM) is a highly conserved mediator of calcium (Ca2+ )-dependent signalling and modulates various cardiac ion channels. Genotyping has revealed several CaM mutations associated with long QT syndrome (LQTS). LQTS patients display prolonged ventricular recovery times (QT interval), increasing their risk of incurring life-threatening arrhythmic events. Loss-of-function mutations to Kv7.1 (which drives the slow delayed rectifier potassium current, IKs, a key ventricular repolarising current) are the largest contributor to congenital LQTS (>50% of cases). CaM modulates Kv7.1 to produce a Ca2+ -sensitive IKs, but little is known about the consequences of LQTS-associated CaM mutations on Kv7.1 function. Here, we present novel data characterising the biophysical and modulatory properties of three LQTS-associated CaM variants (D95V, N97I and D131H). We showed that mutations induced structural alterations in CaM and reduced affinity for Kv7.1, when compared with wild-type (WT). Using HEK293T cells expressing Kv7.1 channel subunits (KCNQ1/KCNE1) and patch-clamp electrophysiology, we demonstrated that LQTS-associated CaM variants reduced current density at systolic Ca2+ concentrations (1 μm), revealing a direct QT-prolonging modulatory effect. Our data highlight for the first time that LQTS-associated perturbations to CaM's structure impede complex formation with Kv7.1 and subsequently result in reduced IKs. This provides a novel mechanistic insight into how the perturbed structure-function relationship of CaM variants contributes to the LQTS phenotype. KEY POINTS: Calmodulin (CaM) is a ubiquitous, highly conserved calcium (Ca2+ ) sensor playing a key role in cardiac muscle contraction. Genotyping has revealed several CaM mutations associated with long QT syndrome (LQTS), a life-threatening cardiac arrhythmia syndrome. LQTS-associated CaM variants (D95V, N97I and D131H) induced structural alterations, altered binding to Kv7.1 and reduced IKs. Our data provide a novel mechanistic insight into how the perturbed structure-function relationship of CaM variants contributes to the LQTS phenotype.

K-Ras Binds Calmodulin-Related Centrin1 with Potential Implications for K-Ras Driven Cancer Cell Stemness

Recent data suggest that K-Ras4B (hereafter K-Ras) can drive cancer cell stemness via calmodulin (CaM)-dependent, non-canonical Wnt-signalling. Here we examined whether another Ca²⁺-binding protein, the CaM-related centrin1, binds to K-Ras and could mediate some K-Ras functions that were previously ascribed to CaM. While CaM and centrin1 appear to distinguish between peptides that were derived from their classical targets, they both bind to K-Ras in cells. Cellular BRET- and immunoprecipitation data suggest that CaM engages more with K-Ras than centrin1 and that the interaction with the C-terminal membrane anchor of K-Ras is sufficient for this. Surprisingly, binding of neither K-Ras nor its membrane anchor alone to CaM or centrin1 is sensitive to inhibition of prenylation. In support of an involvement of the G-domain of K-Ras in cellular complexes with these Ca²⁺-binding proteins, we find that oncogenic K-RasG12V displays increased engagement with both CaM and centrin1. This is abrogated by addition of the D38A effector-site mutation, suggesting that K-RasG12V is held together with CaM or centrin1 in complexes with effectors. When treated with CaM inhibitors, the BRET-interaction of K-RasG12V with centrin1 was also disrupted in the low micromolar range, comparable to that with CaM. While CaM predominates in regulating functional membrane anchorage of K-Ras, it has a very similar co-distribution with centrin1 on mitotic organelles. Given these results, a significant overlap of the CaM- and centrin1-dependent functions of K-Ras is suggested.

Periodicity of the Affinity of Lanmodulin for Trivalent Lanthanides and Actinides: Structural and Electronic Insights from Quantum Chemical Calculations

Lanmodulin (LanM) is the first identified macrochelator that has naturally evolved to sequester ions of rare earth elements (REEs) such as Y and all lanthanides, reversibly. This natural protein showed a 10⁶ times better affinity for lanthanide cations than for Ca, which is a naturally abundant and biologically relevant element. Recent experiments have shown that its metal ion binding activity can be further extended to some actinides, like Np, Pu, and Am. For this reason, it was thought that LanM could be adopted for the separation of REE ions and actinides, thus increasing the interest in its potential use for industry-oriented applications. In this work, a systematic study of the affinity of LanM for lanthanides and actinides has been carried out, taking into account all trivalent ions belonging to the 4f (from La to Lu) and 5f (from Ac to Lr) series, starting from their chemistry in solution. On the basis of a recently published nuclear magnetic resonance structure, a model of the LanM-binding site was built and a detailed structural and electronic description of initial aquo- and LanM-metal ion complexes was provided. The obtained binding energies are in agreement with the available experimental data. A possible reason that could explain the origin of the affinity of LanM for these metal ions is also discussed.

Clickable polymer scaffolds enable Ce recovery with peptide ligands

Rare earth elements (REEs) are a vital part of many technologies with particular importance to the renewable energy sector and there is a pressing need for environmentally friendly and sustainable processes to recover and recycle them from waste streams. Functionalized polymer scaffolds are a promising means to recover REEs due to the ability to engineer both transport properties of the porous material and specificity for target ions. In this work, REE adsorbing polymer scaffolds were synthesized by first introducing poly(glycidyl methacrylate) (GMA) brushes onto porous polyvinylidene fluoride (PVDF) surface through activator generated electron transfer atom transfer radical polymerization (AGET ATRP). Azide moieties were then introduced through a ring opening reaction of GMA. Subsequently, REE-binding peptides were conjugated to the polymer surface through copper catalyzed azide alkyne cycloaddition (CuAAC) click chemistry. The presence of GMA, azide, and peptide was confirmed through Fourier transform infrared spectroscopy. Polymer scaffolds functionalized with the REE-binding peptide bound cerium, while polymer scaffolds functionalized with a scrambled control peptide bound significantly less cerium. Importantly, this study shows that the REE binding peptide retains its functionality when bound to a polymer surface. The conjugation strategy employed in this work can be used to introduce peptides onto other polymeric surfaces and tailor surface specificity for a wide variety of ions and small molecules.

Ca2+-calmodulin signalling at the host-pathogen interface

Multiple eukaryotic cell processes are modulated by calcium ions (Ca²⁺). As such, Ca²⁺ is emerging as a crucial regulator of innate immunity in multicellular organisms. In particular, recent studies have identified roles of Ca²⁺ signalling at the host-bacteria interface. Following microbial exposure, Ca²⁺ signals mobilised from the extracellular milieu or intracellular stores are transduced into cell physiological responses. However, during infection with host-adapted pathogens, Ca²⁺ signals are often atypical, due to the activities of virulence factors, with varied consequences for both the pathogen and the host cell. In this review, we describe the Ca²⁺-dependent host factors regulating antibacterial immunity, in addition to bacterial effectors that promote, inhibit, or co-opt Ca²⁺-calmodulin signalling to promote infection.

Regulation of Cardiac Cav1.2 Channels by Calmodulin

Cav1.2 Ca2+ channels, a type of voltage-gated L-type Ca2+ channel, are ubiquitously expressed, and the predominant Ca2+ channel type, in working cardiac myocytes. Cav1.2 channels are regulated by the direct interactions with calmodulin (CaM), a Ca2+-binding protein that causes Ca2+-dependent facilitation (CDF) and inactivation (CDI). Ca2+-free CaM (apoCaM) also contributes to the regulation of Cav1.2 channels. Furthermore, CaM indirectly affects channel activity by activating CaM-dependent enzymes, such as CaM-dependent protein kinase II and calcineurin (a CaM-dependent protein phosphatase). In this article, we review the recent progress in identifying the role of apoCaM in the channel ‘rundown’ phenomena and related repriming of channels, and CDF, as well as the role of Ca2+/CaM in CDI. In addition, the role of CaM in channel clustering is reviewed.

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.

Calcium Signalling in Heart and Vessels: Role of Calmodulin and Downstream Calmodulin-Dependent Protein Kinases

Cardiovascular disease is the major cause of death worldwide. The success of medication and other preventive measures introduced in the last century have not yet halted the epidemic of cardiovascular disease. Although the molecular mechanisms of the pathophysiology of the heart and vessels have been extensively studied, the burden of ischemic cardiovascular conditions has risen to become a top cause of morbidity and mortality. Calcium has important functions in the cardiovascular system. Calcium is involved in the mechanism of excitation-contraction coupling that regulates numerous events, ranging from the production of action potentials to the contraction of cardiomyocytes and vascular smooth muscle cells. Both in the heart and vessels, the rise of intracellular calcium is sensed by calmodulin, a protein that regulates and activates downstream kinases involved in regulating calcium signalling. Among them is the calcium calmodulin kinase family, which is involved in the regulation of cardiac functions. In this review, we present the current literature regarding the role of calcium/calmodulin pathways in the heart and vessels with the aim to summarize our mechanistic understanding of this process and to open novel avenues for research.

Effects of pre-partum dietary crude protein level on colostrum fat globule membrane proteins and the performance of Hu ewes and their offspring

Dietary proteins play important roles in the growth and reproduction of sheep, and the ewe's demand for proteins increases dramatically during late pregnancy. This research aimed to investigate the effect of dietary crude protein (CP) levels during late pregnancy on colostrum fat globule membrane (MFGM) protein and the growth performance of Hu sheep and their offspring, and provide a reference for the protein intake of ewes during late pregnancy. A total of 108 multiparous Hu sheep (45.6 ± 1.18 kg) were selected for this study, then 60 pregnant ewes confirmed by B-scan ultrasonography were randomly divided into three treatments (20 ewes/treatment) and fed by total mixed ration pellet with CP levels at 9.00% (LP), 12.0% (MP), and 15.0% (HP) during late pregnancy, respectively. The weight and dry matter intake of ewes during late pregnancy were recorded to calculate the average daily gain (ADG) and feed conversion ratio (FCR). Twin lambs were weighed on days 0, 7, 14, 30, 60, and 180 after birth to calculate ADG. Meanwhile, the colostrum of ewes was collected within 12 h after delivery. The colostrum MFGM proteins were identified and quantified by the isobaric tag for relative and absolute quantification (iTRAQ) coupled with liquid chromatography-tandem mass spectrometry methods. In addition, biological functions of differentially expressed proteins (DEPs) were annotated by Gene Ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. The results revealed that a 15.0% CP level had significant effects on the BW of lambs on days 0, 7, and 30 (P < 0.05). Notably, a total of 1,529 MFGM proteins were identified and 286 DEPs were found among three treatments. Functional analysis showed that DEPs were mainly involved in cell growth, differentiation, and tissue repair, and involved in metabolic pathways, such as the porphyrin and chlorophyll metabolism pathways. In this study, lambs in HP treatment had better growth performance; moreover, dietary 15.0% CP level also affected the colostrum MFGM proteins composition of Hu ewes. These observations can facilitate future studies on the feeding regimen of ewes during late pregnancy.

Threonine Phosphorylation of an Electrochemical Peptide-Based Sensor to Achieve Improved Uranyl Ion Binding Affinity

We have successfully designed a uranyl ion (U(VI)-specific peptide and used it in the fabrication of an electrochemical sensor. The 12-amino acid peptide sequence, (n) DKDGDGYIpTAAE (c), originates from calmodulin, a Ca(II)-binding protein, and contains a phosphothreonine that enhances the sequence's affinity for U(VI) over Ca(II). The sensing mechanism of this U(VI) sensor is similar to other electrochemical peptide-based sensors, which relies on the change in the flexibility of the peptide probe upon interacting with the target. The sensor was systematically characterized using alternating current voltammetry (ACV) and cyclic voltammetry. Its limit of detection was 50 nM, which is lower than the United States Environmental Protection Agency maximum contaminant level for uranium. The signal saturation time was ~40 min. In addition, it showed minimal cross-reactivity when tested against nine different metal ions, including Ca(II), Mg(II), Pb(II), Hg(II), Cu(II), Fe(II), Zn(II), Cd(II), and Cr(VI). Its reusability and ability to function in diluted aquifer and drinking water samples were further confirmed and validated. The response of the sensor fabricated with the same peptide sequence but with a nonphosphorylated threonine was also analyzed, substantiating the positive effects of threonine phosphorylation on U(VI) binding. This study places emphasis on strategic utilization of non-standard amino acids in the design of metal ion-chelating peptides, which will further diversify the types of peptide recognition elements available for metal ion sensing applications.

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's Interdomain Linker Is Optimized for Dynamics Signal Transmission and Calcium Binding

Linkers are ubiquitous in multidomain proteins. These linkers are integral to protein functions, and accumulating evidence suggests that the linkers' versatile roles are encoded in their sequences. However, a molecular picture of how amino acid differences in the linker influence protein function is still lacking. By using extensive Gaussian-accelerated MD coupled with dynamic network analysis, we reveal the molecular bases underlying the linker's role in Calmodulin (CaM), a highly conserved Ca²⁺-signaling hub in eukaryotes. Three CaM constructs comprising a wild-type linker, a flexible linker (four glycines at position D78-S81), and a rigid linker (four prolines at position D78-S81) were simulated. We show that the flexible linker resembles the wild type in allowing CaM to sample a large ensemble of conformations while the rigid linker confines the sampling. Our simulations recapture experimental observations that target binding enhances the Ca²⁺ affinity to CaM's EF-hand sites at the N-domain. However, only the wild-type linker can both correctly capture the Ca²⁺ binding order and maintain the α-helical structure of the domain. The other two constructs either bind Ca²⁺ in an incorrect order or exhibit unfolding of an N-domain helix. We demonstrate that the wild-type linker achieves these outcomes by transmitting interdomain dynamics efficiently. This was evidenced by stronger (anti)correlations among the linker residues, decoupling of the hydrogen bonds between A1-A15 and V35-E45, and structuring of the N-domain for Ca²⁺ binding. This decoupling was not evident for the other two constructs. Lastly, we show that the wild-type linker's optimal transmission stems from its thermodynamically favorable strain and solvation relative to the other two constructs. Our results show how the linker sequence tunes CaM function, suggesting possible mechanisms for changes in linker properties such as mutations or post-translational modifications to modulate protein/substrate binding.

Fibrillarin RNA methylase is an interacting protein of Cryptosporidium parvum calmodulin-like protein (CpCML)

Cryptosporidium parvum is an obligate protozoan parasite invading epithelial cells of small intestine of human and animals, and causing diarrheal disease. In apicomplexan parasites, calcium signaling can regulate many essential biological processes such as invasion and migration. As the main intracellular receptor for calcium ions, calmodulins control the activities of hundreds of enzymes and proteins. Calmodulin-like protein (CML) is an important member of the calmodulin family and may play a key role in C. parvum, however, the actual situation is still not clear. The present study aimed to identify the parasite interaction partner proteins of C. parvum calmodulin-like protein (CpCML). By constructing the cpcml bait plasmid, 5 potential CpCML - interacting proteins in C. parvum oocyst were screened by yeast-two-hybrid system (Y2H). Bimolecular fluorescence complementation (BiFC) and Co-immunoprecipitation (Co-IP) were performed as subsequent validations. Fibrillarin RNA methylase (FBL) was identified via this screening method as CpCML interacting protein in C. parvum. The identification of this interaction made it possible to get a further understanding of the function of CpCML and its contribution to the pathogenicity of C. parvum.

The Dually Localized EF-Hand Domain-Containing Protein TgEFP1 Regulates the Lytic Cycle of Toxoplasma gondii

The propagation of the obligate intracellular parasite Toxoplasma gondii is tightly regulated by calcium signaling. However, the mechanisms by which calcium homeostasis and fluxes are regulated in this human pathogen are not fully understood. To identify Toxoplasma's calcium homeostasis network, we have characterized a novel EF-hand domain-containing protein, which we have named TgEFP1. We have determined that TgEFP1 localizes to a previously described compartment known as the plant-like vacuole or the endosomal-like compartment (PLV/ELC), which harbors several proteins related to ionic regulation. Interestingly, partial permeabilization techniques showed that TgEFP1 is also secreted into the parasitophorous vacuole (PV), within which the parasite divides. Ultrastructure expansion microscopy confirmed the unusual dual localization of TgEFP1 at the PLV/ELC and the PV. Furthermore, we determined that the localization of TgEFP1 to the PV, but not to the PLV/ELC, is affected by disruption of Golgi-dependent transport with Brefeldin A. Knockout of TgEFP1 results in faster propagation in tissue culture, hypersensitivity to calcium ionophore-induced egress, and premature natural egress. Thus, our work has revealed an interplay between the PV and the PLV/ELC and a role for TgEFP1 in the regulation of calcium-dependent events.

Engineered CaM2 modulates nuclear calcium oscillation and enhances legume root nodule symbiosis

Oscillations in intracellular calcium concentration play an essential role in the regulation of multiple cellular processes. In plants capable of root endosymbiosis with nitrogen-fixing bacteria and/or arbuscular mycorrhizal fungi, nuclear localized calcium oscillations are essential to transduce the microbial signal. Although the ion channels required to generate the nuclear localized calcium oscillations have been identified, their mechanisms of regulation are unknown. Here, we combined proteomics and engineering approaches to demonstrate that the calcium-bound form of the calmodulin 2 (CaM2) associates with CYCLIC NUCLEOTIDE GATED CHANNEL 15 (CNGC15s), closing the channels and providing the negative feedback to sustain the oscillatory mechanism. We further unraveled that the engineered CaM2 accelerates early endosymbioses and enhanced root nodule symbiosis but not arbuscular mycorrhization.

The key physiological event essential to the establishment of nitrogen-fixing bacteria and phosphate-delivering arbuscular mycorrhizal symbioses is the induction of nuclear calcium oscillations that are required for endosymbioses. These regular fluctuations in nucleoplasmic calcium concentrations are generated by ion channels and a pump located at the nuclear envelope, including the CYCLIC NUCLEOTIDE GATED CHANNEL 15 (CNGC15). However, how the CNGC15s are regulated in planta to sustain a calcium oscillatory mechanism remains unknown. Here, we demonstrate that the CNGC15s are regulated by the calcium-bound form of the calmodulin 2 (holo-CaM2), which, upon release of calcium, provides negative feedback to close the CNGC15s. Combining structural and evolutionary analyses of CaM residues with bioinformatic analysis, we engineered a holo-CaM2 with an increased affinity for CNGC15s. In planta, the expression of the engineered holo-CaM2 accelerates the calcium oscillation frequency, early endosymbioses signaling and is sufficient to sustain over time an enhanced root nodule symbiosis but not an increased arbuscular mycorrhization. Together, these results reveal that holo-CaM2 is a component of endosymbiosis signaling required to modulate CNGC15s activity and the downstream root nodule symbiosis pathway.

Overcoming universal restrictions on metal selectivity by protein design

Selective metal coordination is central to the functions of metalloproteins:1,2 each metalloprotein must pair with its cognate metallocofactor to fulfil its biological role3. However, achieving metal selectivity solely through a three-dimensional protein structure is a great challenge, because there is a limited set of metal-coordinating amino acid functionalities and proteins are inherently flexible, which impedes steric selection of metals3,4. Metal-binding affinities of natural proteins are primarily dictated by the electronic properties of metal ions and follow the Irving-Williams series5 (Mn2+ < Fe2+ < Co2+ < Ni2+ < Cu2+ > Zn2+) with few exceptions6,7. Accordingly, metalloproteins overwhelmingly bind Cu2+ and Zn2+ in isolation, regardless of the nature of their active sites and their cognate metal ions1,3,8. This led organisms to evolve complex homeostatic machinery and non-equilibrium strategies to achieve correct metal speciation1,3,8-10. Here we report an artificial dimeric protein, (AB)2, that thermodynamically overcomes the Irving-Williams restrictions in vitro and in cells, favouring the binding of lower-Irving-Williams transition metals over Cu2+, the most dominant ion in the Irving-Williams series. Counter to the convention in molecular design of achieving specificity through structural preorganization, (AB)2 was deliberately designed to be flexible. This flexibility enabled (AB)2 to adopt mutually exclusive, metal-dependent conformational states, which led to the discovery of structurally coupled coordination sites that disfavour Cu2+ ions by enforcing an unfavourable coordination geometry. Aside from highlighting flexibility as a valuable element in protein design, our results illustrate design principles for constructing selective metal sequestration agents.

Calmodulin and calmodulin-like gene family in barley: Identification, characterization and expression analyses

Calmodulin (CaM) and calmodulin-like (CML) proteins are Ca²⁺ relays and play diverse and multiple roles in plant growth, development and stress responses. However, CaM/CML gene family has not been identified in barley (Hordeum vulgare). In the present study, 5 HvCaMs and 80 HvCMLs were identified through a genome-wide analysis. All HvCaM proteins possessed 4 EF-hand motifs, whereas HvCMLs contained 1 to 4 EF-hand motifs. HvCaM2, HvCaM3 and HvCaM5 coded the same polypeptide although they differed in nucleotide sequence, which was identical to the polypeptides coded by OsCaM1-1, OsCaM1-2 and OsCaM1-3. HvCaMs/CMLs were unevenly distributed over barley 7 chromosomes, and could be phylogenetically classified into 8 groups. HvCaMs/CMLs differed in gene structure, cis-acting elements and tissue expression patterns. Segmental and tandem duplication were observed among HvCaMs/CMLs during evolution. HvCML16, HvCML18, HvCML50 and HvCML78 were dispensable genes and the others were core genes in barley pan-genome. In addition, 14 HvCaM/CML genes were selected to examine their responses to salt, osmotic and low potassium stresses by qRT-PCR, and their expression were stress-and time-dependent. These results facilitate our understanding and further functional identification of HvCaMs/CMLs.

Engineering lanmodulin's selectivity for actinides over lanthanides by controlling solvent coordination and second-sphere interactions

Developing chelators that combine high affinity and selectivity for lanthanides and/or actinides is paramount for numerous industries, including rare earths mining, nuclear waste management, and cancer medicine. In particular, achieving selectivity between actinides and lanthanides is notoriously difficult. The protein lanmodulin (LanM) is one of Nature's most selective chelators for trivalent actinides and lanthanides. However, mechanistic understanding of LanM's affinity and selectivity for f-elements remains limited. In order to decipher, and possibly improve, the features of LanM's metal-binding sites that contribute to this actinide/lanthanide selectivity, we characterized five LanM variants, substituting the aspartate residue at the 9^th position of each metal-binding site with asparagine, histidine, alanine, methionine, and selenomethionine. Spectroscopic measurements with lanthanides (Nd³⁺ and Eu³⁺) and actinides (²⁴³Am³⁺ and ²⁴⁸Cm³⁺) reveal that, contrary to the behavior of small chelator complexes, metal-coordinated water molecules enhance LanM's affinity for f-elements and pH-stability of its complexes. Furthermore, the results show that the native aspartate does not coordinate the metal directly but rather hydrogen bonds to coordinated solvent. By tuning this first-sphere/second-sphere interaction, the asparagine variant nearly doubles LanM's selectivity for actinides versus lanthanides. This study not only clarifies the essential role of coordinated solvent for LanM's physiological function and separation applications, but it also demonstrates that LanM's preference for actinides over lanthanides can be further improved. More broadly, it demonstrates how biomolecular scaffolds possess an expanded repertoire of tunable interactions compared to most small-molecule ligands – providing an avenue for high-performance LanM-based actinide/lanthanide separation methods and bio-engineered chelators optimized for specific medical isotopes.

Nature’s most potent protein for f-elements, lanmodulin, relies on subtle first-sphere/second-sphere interactions to bind metal ions. Dissecting lanmodulin’s binding mechanism yielded variants with enhanced actinide/lanthanide selectivity.

Lanmodulin peptides – unravelling the binding of the EF-Hand loop sequences stripped from the structural corset

Lanmodulin (LanM), a naturally lanthanide (Ln)-binding protein with a remarkable selectivity for Lns over Ca(ii) and affinities in the picomolar range, is an attractive target to address challenges in Ln separation. Why LanM has such a high selectivity is currently not entirely understood; both specific amino acid sequences of the EF-Hand loops and cooperativity effects have been suggested. Here, we removed the effect of cooperativity and synthesised all four 12-amino acid EF-Hand loop peptides, and investigated their affinity for two Lns (Eu(iii) and Tb(iii)), the actinide Cm(iii) and Ca(ii). Using isothermal titration calorimetry and time-resolved laser fluorescence spectroscopy (TRLFS) combined with parallel factor analysis, we show that the four short peptides behave very similarly, having affinities in the micromolar range for Eu(iii) and Tb(iii). Ca(ii) was shown not to bind to the peptides, which was verified with circular dichroism spectroscopy. This technique also revealed an increase in structural organisation upon Eu(iii) addition, which was supported by molecular dynamics simulations. Lastly, we put Eu(iii) and Cm(iii) in direct competition using TRLFS. Remarkably, a slightly higher affinity for Cm(iii) was found. Our results demonstrate that the picomolar affinities in LanM are largely an effect of pre-structuring and therefore a reduction of flexibility in combination with cooperative effects, and that all EF-Hand loops possess similar affinities when detached from the protein backbone, albeit still retaining the high selectivity for lanthanides and actinides over calcium.

Taking a closer look at Lanmodulin’s remarkable selectivity for lanthanides (Ln) over Ca(ii) and high Ln/actinide affinities on the amino acid level by investigating the four binding-loops as peptides with Ca(ii), Eu(iii), Tb(iii) and Cm(iii).

Genome-Wide Analysis of the Soybean Calmodulin-Binding Protein 60 Family and Identification of GmCBP60A-1 Responses to Drought and Salt Stresses

Calmodulin-binding protein 60 (CBP60) members constitute a plant-specific protein family that plays an important role in plant growth and development. In the soybean genome, nineteen CBP60 members were identified and analyzed for their corresponding sequences and structures to explore their functions. Among GmCBP60A-1, which primarily locates in the cytomembrane, was significantly induced by drought and salt stresses. The overexpression of GmCBP60A-1 enhanced drought and salt tolerance in Arabidopsis, which showed better state in the germination of seeds and the root growth of seedlings. In the soybean hairy roots experiment, the overexpression of GmCBP60A-1 increased proline content, lowered water loss rate and malondialdehyde (MDA) content, all of which likely enhanced the drought and salt tolerance of soybean seedlings. Under stress conditions, drought and salt response-related genes showed significant differences in expression in hairy root soybean plants of GmCBP60A-1-overexpressing and hairy root soybean plants of RNAi. The present study identified GmCBP60A-1 as an important gene in response to salt and drought stresses based on the functional analysis of this gene and its potential underlying mechanisms in soybean stress-tolerance.

A turquoise fluorescence lifetime-based biosensor for quantitative imaging of intracellular calcium

The most successful genetically encoded calcium indicators (GECIs) employ an intensity or ratiometric readout. Despite a large calcium-dependent change in fluorescence intensity, the quantification of calcium concentrations with GECIs is problematic, which is further complicated by the sensitivity of all GECIs to changes in the pH in the biological range. Here, we report on a sensing strategy in which a conformational change directly modifies the fluorescence quantum yield and fluorescence lifetime of a circular permutated turquoise fluorescent protein. The fluorescence lifetime is an absolute parameter that enables straightforward quantification, eliminating intensity-related artifacts. An engineering strategy that optimizes lifetime contrast led to a biosensor that shows a 3-fold change in the calcium-dependent quantum yield and a fluorescence lifetime change of 1.3 ns. We dub the biosensor Turquoise Calcium Fluorescence LIfeTime Sensor (Tq-Ca-FLITS). The response of the calcium sensor is insensitive to pH between 6.2-9. As a result, Tq-Ca-FLITS enables robust measurements of intracellular calcium concentrations by fluorescence lifetime imaging. We demonstrate quantitative imaging of calcium concentrations with the turquoise GECI in single endothelial cells and human-derived organoids.

Ca2+-dependent modulation of voltage-gated myocyte sodium channels

Voltage-dependent Na⁺ channel activation underlies action potential generation fundamental to cellular excitability. In skeletal and cardiac muscle this triggers contraction via ryanodine-receptor (RyR)-mediated sarcoplasmic reticular (SR) Ca²⁺ release. We here review potential feedback actions of intracellular [Ca²⁺] ([Ca²⁺]_i) on Na⁺ channel activity, surveying their structural, genetic and cellular and functional implications, translating these to their possible clinical importance. In addition to phosphorylation sites, both Nav1.4 and Nav1.5 possess potentially regulatory binding sites for Ca²⁺ and/or the Ca^2+-sensor calmodulin in their inactivating III–IV linker and C-terminal domains (CTD), where mutations are associated with a range of skeletal and cardiac muscle diseases. We summarize in vitro cell-attached patch clamp studies reporting correspondingly diverse, direct and indirect, Ca²⁺ effects upon maximal Nav1.4 and Nav1.5 currents (I_max) and their half-maximal voltages (V_1/2) characterizing channel gating, in cellular expression systems and isolated myocytes. Interventions increasing cytoplasmic [Ca²⁺]_i down-regulated I_max leaving V_1/2 constant in native loose patch clamped, wild-type murine skeletal and cardiac myocytes. They correspondingly reduced action potential upstroke rates and conduction velocities, causing pro-arrhythmic effects in intact perfused hearts. Genetically modified murine RyR2-P2328S hearts modelling catecholaminergic polymorphic ventricular tachycardia (CPVT), recapitulated clinical ventricular and atrial pro-arrhythmic phenotypes following catecholaminergic challenge. These accompanied reductions in action potential conduction velocities. The latter were reversed by flecainide at RyR-blocking concentrations specifically in RyR2-P2328S as opposed to wild-type hearts, suggesting a basis for its recent therapeutic application in CPVT. We finally explore the relevance of these mechanisms in further genetic paradigms for commoner metabolic and structural cardiac disease.

The mechanism of complex formation between calmodulin and voltage gated calcium channels revealed by molecular dynamics

Calmodulin, a ubiquitous eukaryotic calcium sensor responsible for the regulation of many fundamental cellular processes, is a highly flexible protein and exhibits an unusually wide range of conformations. Furthermore, CaM is known to interact with more than 300 cellular targets. Molecular dynamics (MD) simulation trajectories suggest that EF-hand loops show different magnitudes of flexibility. Therefore, the four EF-hand motifs have different affinities for Ca2+ ions, which enables CaM to function on wide range of Ca2+ ion concentrations. EF-hand loops are 2-3 times more flexible in apo CaM whereas least flexible in Ca2+/CaM-IQ motif complexes. We report a unique intermediate conformation of Ca2+/CaM while transitioning from extended to compact form. We also report the complex formation process between Ca2+/CaM and IQ CaM-binding motifs. Our results showed how IQ motif recognise its binding site on the CaM and how CaM transforms from extended to compact form upon binding to IQ motif.

Lanthanide-dependent coordination interactions in lanmodulin: a 2D IR and molecular dynamics simulations study

The biological importance of lanthanides, and the early lanthanides (La³⁺-Nd³⁺) in particular, has only recently been recognized, and the structural principles underlying selective binding of lanthanide ions in biology are not yet well established. Lanmodulin (LanM) is a novel protein that displays unprecedented affinity and selectivity for lanthanides over most other metal ions, with an uncommon preference for the early lanthanides. Its utilization of EF-hand motifs to bind lanthanides, rather than the Ca²⁺ typically recognized by these motifs in other proteins, has led it to be used as a model system to understand selective lanthanide recognition. Two-dimensional infrared (2D IR) spectroscopy combined with molecular dynamics simulations were used to investigate LanM's selectivity mechanisms by characterizing local binding site geometries upon coordination of early and late lanthanides as well as calcium. These studies focused on the protein's uniquely conserved proline residues in the second position of each EF-hand binding loop. We found that these prolines constrain the EF-hands for strong coordination of early lanthanides. Substitution of this proline results in a more flexible binding site to accommodate a larger range of ions but also results in less compact coordination geometries and greater disorder within the binding site. Finally, we identify the conserved glycine in the sixth position of each EF-hand as a mediator of local binding site conformation and global secondary structure. Uncovering fundamental structure-function relationships in LanM informs the development of synthetic biology technologies targeting lanthanides in industrial applications.

Conformational scanning of individual EF-hand motifs of calcium sensor protein centrin-1

Centrin-1, a Ca²⁺ sensor protein of the centrin family is a crucial player for cell division in eukaryotes and plays a key role in the microtubule organising centre. Despite being regarded as a calcium sensor with a matched structure to calmodulin/troponin C, the protein undergoes mild changes in conformation and binds Ca²⁺ with moderate affinity. We present an in-depth analysis of the Ca²⁺ sensing by individual EF-hand motifs of centrin-1 and address unsolved questions of the rationales for moderate affinity and conformational transitions of the protein. Employing the more sensitive approach of Trp scanning of individual EF-hand motif, we have undertaken an exhaustive investigation of Ca²⁺ binding to individual EF-hand motifs, named EF1 to EF4. All four EF-hand motifs of centrin-1 are structural as all of them bind both Ca²⁺ and Mg²⁺. EF1 and EF4 are the most flexible sites as they undergo drastic conformational changes following Ca²⁺ binding, whereas EF3 responds to Ca²⁺ minimally. On the other hand, EF2 moves towards the protein surface upon binding Ca²⁺. The independent filling mode of Ca²⁺ to EF-hand motifs and lack of intermotif communication explain the lack of cooperativity of binding, thus constraining centrin-1 to a moderate affinity binding protein. Thus, centrin-1 is distinct from other calcium sensors such as calmodulin.

Influence of Subcellular Localization and Functional State on Protein Turnover

Protein homeostasis is an equilibrium of paramount importance that maintains cellular performance by preserving an efficient proteome. This equilibrium avoids the accumulation of potentially toxic proteins, which could lead to cellular stress and death. While the regulators of proteostasis are the machineries controlling protein production, folding and degradation, several other factors can influence this process. Here, we have considered two factors influencing protein turnover: the subcellular localization of a protein and its functional state. For this purpose, we used an imaging approach based on the pulse-labeling of 17 representative SNAP-tag constructs for measuring protein lifetimes. With this approach, we obtained precise measurements of protein turnover rates in several subcellular compartments. We also tested a selection of mutants modulating the function of three extensively studied proteins, the Ca²⁺ sensor calmodulin, the small GTPase Rab5a and the brain creatine kinase (CKB). Finally, we followed up on the increased lifetime observed for the constitutively active Rab5a (Q79L), and we found that its stabilization correlates with enlarged endosomes and increased interaction with membranes. Overall, our data reveal that both changes in protein localization and functional state are key modulators of protein turnover, and protein lifetime fluctuations can be considered to infer changes in cellular behavior.

A Covalent Calmodulin Inhibitor as a Tool to Study Cellular Mechanisms of K-Ras-Driven Stemness

Recently, the highly mutated oncoprotein K-Ras4B (hereafter K-Ras) was shown to drive cancer cell stemness in conjunction with calmodulin (CaM). We previously showed that the covalent CaM inhibitor ophiobolin A (OphA) can potently inhibit K-Ras stemness activity. However, OphA, a fungus-derived natural product, exhibits an unspecific, broad toxicity across all phyla. Here we identified a less toxic, functional analog of OphA that can efficiently inactivate CaM by covalent inhibition. We analyzed a small series of benzazulenones, which bear some structural similarity to OphA and can be synthesized in only six steps. We identified the formyl aminobenzazulenone 1, here named Calmirasone1, as a novel and potent covalent CaM inhibitor. Calmirasone1 has a 4-fold increased affinity for CaM as compared to OphA and was active against K-Ras in cells within minutes, as compared to hours required by OphA. Calmirasone1 displayed a 2.5–4.5-fold higher selectivity for KRAS over BRAF mutant 3D spheroid growth than OphA, suggesting improved relative on-target activity. Importantly, Calmirasone1 has a 40–260-fold lower unspecific toxic effect on HRAS mutant cells, while it reaches almost 50% of the activity of novel K-RasG12C specific inhibitors in 3D spheroid assays. Our results suggest that Calmirasone1 can serve as a new tool compound to further investigate the cancer cell biology of the K-Ras and CaM associated stemness activities.

CH Mode Mixing Determines the Band Shape of the Carboxylate Symmetric Stretch in Apo-EDTA, Ca2+-EDTA, and Mg2+-EDTA

The infrared spectra of EDTA complexed with Ca²⁺ and Mg²⁺ contain, to date, unidentified vibrational bands. This study assigns the peaks in the linear and two-dimensional infrared spectra of EDTA, with and without either Ca²⁺ or Mg²⁺ ions. Two-dimensional infrared spectroscopy and DFT calculations reveal that, in both the presence and absence of ions, the carboxylate symmetric stretch and the terminal CH bending vibrations mix. We introduce a method to calculate participation coefficients that quantify the contribution of the carboxylate symmetric stretch, CH wag, CH twist, and CH scissor in the 1400-1550 cm^-1 region. With the help of participation coefficients, we assign the 1400-1430 cm^-1 region to the carboxylate symmetric stretch, which can mix with CH modes. We assign the 1000-1380 cm^-1 region to CH twist modes, the 1380-1430 cm^-1 region to wag modes, and the 1420-1650 cm^-1 region to scissor modes. The difference in binding geometry between the carboxylate-Ca²⁺ and carboxylate-Mg²⁺ complex manifests as new diagonal and cross-peaks between the mixed modes in the two complexes. The small Mg²⁺ ion binds EDTA tighter than the Ca²⁺ ion, which causes a redshift of the COO symmetric stretches of the sagittal carboxylates. Energy decomposition analysis further characterizes the importance of electrostatics and deformation energy in the bound complexes.

Infrared spectroscopy probes ion binding geometries

Infrared (IR) spectroscopy is a well-established technique for probing the structure, behavior, and surroundings of molecules in their native environments. Its characteristics-most specifically high structural sensitivity, ready applicability to aqueous samples, and broad availability-make it a valuable enzymological technique, particularly for the interrogation of ion binding sites. While IR spectroscopy of the "garden variety" (steady state at room temperature with wild-type proteins) is versatile and powerful in its own right, the combination of IR spectroscopy with specialized experimental schemes for leveraging ultrafast time resolution, protein labeling, and other enhancements further extends this utility. This book chapter provides the fundamental physical background and literature context essential for harnessing IR spectroscopy in the general context of enzymology with specific focus on interrogation of ion binding. Studies of lanthanide ions binding to calmodulin are highlighted as illustrative examples of this process. Appropriate sample preparation, data collection, and spectral interpretation are discussed from a detail-oriented and practical perspective with the goal of facilitating the reader's rapid progression from reading words in a book to collecting and analyzing their own data in the lab.

Genome-Wide Identification and Characterization of Calcium Metabolism Related Gene Families in Arabidopsis thaliana and Their Regulation by Bacillus amyloliquefaciens Under High Calcium Stress

Several gene families involved in calcium signaling have been detected in plants, including calmodulin (CaM), calcium dependent protein kinases (CDPK), calcineurin B-like (CBL) and cyclic nucleotide-gated channels (CNGCs). In our previous study, we demonstrated that Bacillus amyloliquefaciens LZ04 (B. amyloliquefaciens LZ04) regulate genes involved in calcium stress in Arabidopsis thaliana (A. thaliana). Here, we aimed to explore the potential involvement of calcium-related gene families in the response of A. thaliana to calcium stress and the potential regulatory effects of B. amyloliquefaciens LZ04 on these genes. The structure, duplication, synteny, and expression profiles of 102 genes in calcium-related gene families in A. thaliana were investigated. Hidden Markov Models (HMMs) and BLASTP were used to predict candidate genes and conserved domains of the candidate genes were confirmed in SMART and NCBI CDD databases. Gene duplications and synteny were uncovered by BLASTP and phylogenetic analysis. The transcriptome expression profiles of candidate genes were investigated by strand-specific sequencing. Cluster analysis was used to find the expression profiles of calcium-related genes families under different treatment conditions. A total of 102 genes in calcium-related gene families were detected in A. thaliana genome, including 34 CDPK genes, 20 CNGC genes, 18 CIPK genes, 22 IQD genes, and 10 CBP genes. Additionally, of the 102 genes, 33 duplications (32.35%) and 26 gene pairs including 48 genes (47.06%) were detected. Treatment with B. amyloliquefaciens LZ04 enhanced the resistance of A. thaliana under high calcium stress by regulating some of the genes in the calcium-related gene families. Functional enrichment analysis revealed that the genes clustered in the 42nd expression profile which may be B. amyloliquefaciens-responsive genes under calcium stress were enriched in protein phosphorylation and protein modification process. Transcriptome data was validated by RT-PCR and the results generally corroborated the transcriptome sequencing results. These results may be useful for agricultural improvement in high calcium stress regions.

Structural Aspects and Prediction of Calmodulin-Binding Proteins

Calmodulin (CaM) is an important intracellular protein that binds Ca²⁺ and functions as a critical second messenger involved in numerous biological activities through extensive interactions with proteins and peptides. CaM’s ability to adapt to binding targets with different structures is related to the flexible central helix separating the N- and C-terminal lobes, which allows for conformational changes between extended and collapsed forms of the protein. CaM-binding targets are most often identified using prediction algorithms that utilize sequence and structural data to predict regions of peptides and proteins that can interact with CaM. In this review, we provide an overview of different CaM-binding proteins, the motifs through which they interact with CaM, and shared properties that make them good binding partners for CaM. Additionally, we discuss the historical and current methods for predicting CaM binding, and the similarities and differences between these methods and their relative success at prediction. As new CaM-binding proteins are identified and classified, we will gain a broader understanding of the biological processes regulated through changes in Ca²⁺ concentration through interactions with CaM.

The Role of Calcium Signaling in Regulation of Epithelial-Mesenchymal Transition

Despite substantial advances in the field of cancer therapeutics, metastasis is a significant challenge for a favorable clinical outcome. Epithelial to mesenchymal transition (EMT) is a process of acquiring increased motility, invasiveness, and therapeutic resistance by cancer cells for their sustained growth and survival. A plethora of intrinsic mechanisms and extrinsic microenvironmental factors drive the process of cancer metastasis. Calcium (Ca2+) signaling plays a critical role in dictating the adaptive metastatic cell behavior comprising of cell migration, invasion, angiogenesis, and intravasation. By modulating EMT, Ca2+ signaling can regulate the complexity and dynamics of events leading to metastasis. This review summarizes the role of Ca2+ signal remodeling in the regulation of EMT and metastasis in cancer.

Calmodulin mutations and life-threatening cardiac arrhythmias: insights from the International Calmodulinopathy Registry

AIMS

Calmodulinopathies are rare life-threatening arrhythmia syndromes which affect mostly young individuals and are, caused by mutations in any of the three genes (CALM 1-3) that encode identical calmodulin proteins. We established the International Calmodulinopathy Registry (ICalmR) to understand the natural history, clinical features, and response to therapy of patients with a CALM-mediated arrhythmia syndrome.

METHODS AND RESULTS

A dedicated Case Report File was created to collect demographic, clinical, and genetic information. ICalmR has enrolled 74 subjects, with a variant in the CALM1 (n = 36), CALM2 (n = 23), or CALM3 (n = 15) genes. Sixty-four (86.5%) were symptomatic and the 10-year cumulative mortality was 27%. The two prevalent phenotypes are long QT syndrome (LQTS; CALM-LQTS, n = 36, 49%) and catecholaminergic polymorphic ventricular tachycardia (CPVT; CALM-CPVT, n = 21, 28%). CALM-LQTS patients have extremely prolonged QTc intervals (594 ± 73 ms), high prevalence (78%) of life-threatening arrhythmias with median age at onset of 1.5 years [interquartile range (IQR) 0.1-5.5 years] and poor response to therapies. Most electrocardiograms (ECGs) show late onset peaked T waves. All CALM-CPVT patients were symptomatic with median age of onset of 6.0 years (IQR 3.0-8.5 years). Basal ECG frequently shows prominent U waves. Other CALM-related phenotypes are idiopathic ventricular fibrillation (IVF, n = 7), sudden unexplained death (SUD, n = 4), overlapping features of CPVT/LQTS (n = 3), and predominant neurological phenotype (n = 1). Cardiac structural abnormalities and neurological features were present in 18 and 13 patients, respectively.

CONCLUSION

Calmodulinopathies are largely characterized by adrenergically-induced life-threatening arrhythmias. Available therapies are disquietingly insufficient, especially in CALM-LQTS. Combination therapy with drugs, sympathectomy, and devices should be considered.

Regulation of Connexin Gap Junctions and Hemichannels by Calcium and Calcium Binding Protein Calmodulin

Connexins are the structural components of gap junctions and hemichannels that mediate the communication and exchange of small molecules between cells, and between the intracellular and extracellular environment, respectively. Connexin (Cx) 46 is predominately expressed in lens fiber cells, where they function in maintaining the homeostasis and transparency of the lens. Cx46 mutations are associated with impairment of channel function, which results in the development of congenital cataracts. Cx46 gap junctions and hemichannels are closely regulated by multiple mechanisms. Key regulators of Cx46 channel function include Ca2+ and calmodulin (CaM). Ca2+ plays an essential role in lens homeostasis, and its dysregulation causes cataracts. Ca2+ associated CaM is a well-established inhibitor of gap junction coupling. Recent studies suggest that elevated intracellular Ca2+ activates Cx hemichannels in lens fiber cells and Cx46 directly interacts with CaM. A Cx46 site mutation (Cx46-G143R), which is associated with congenital Coppock cataracts, shows an increased Cx46-CaM interaction and this interaction is insensitive to Ca2+, given that depletion of Ca2+ reduces the interaction between CaM and wild-type Cx46. Moreover, inhibition of CaM function greatly reduces the hemichannel activity in the Cx46 G143R mutant. These research findings suggest a new regulatory mechanism by which enhanced association of Cx46 with CaM leads to the increase in hemichannel activity and dysregulation may lead to cataract development. In this review, we will first discuss the involvement of Ca2+/CaM in lens homeostasis and pathology, and follow by providing a general overview of Ca2+/CaM in the regulation of Cx46 gap junctions. We discuss the most recent studies concerning the molecular mechanism of Ca2+/CaM in regulating Cx46 hemichannels. Finally, we will offer perspectives of the impacts of Ca2+/CaM and dysregulation on Cx46 channels and vice versa.

The biochemistry of lanthanide acquisition, trafficking, and utilization

Lanthanides are relative newcomers to the field of cell biology of metals; their specific incorporation into enzymes was only demonstrated in 2011, with the isolation of a bacterial lanthanide- and pyrroloquinoline quinone-dependent methanol dehydrogenase. Since that discovery, the efforts of many investigators have revealed that lanthanide utilization is widespread in environmentally important bacteria, and parallel efforts have focused on elucidating the molecular details involved in selective recognition and utilization of these metals. In this review, we discuss the particular chemical challenges and advantages associated with biology's use of lanthanides, as well as the currently known lanthano-enzymes and -proteins (the lanthanome). We also review the emerging understanding of the coordination chemistry and biology of lanthanide acquisition, trafficking, and regulatory pathways. These studies have revealed significant parallels with pathways for utilization of other metals in biology. Finally, we discuss some of the many unresolved questions in this burgeoning field and their potentially far-reaching applications.

Gene polymorphisms in calcium-calmodulin pathway: Focus on cardiovascular disease

Cardiovascular disease is the leading cause of death in industrialized countries and affects an increasing number of people. Several risk factors play an important role in the etiology of this disease, such as an unhealthy lifestyle. It is increasingly clear that genetic factors influencing the molecular basis of excitation-contraction mechanisms in the heart could contribute to modify the individual's risk. Thanks to the progress that has been made in understanding calcium signaling in the heart, it is assumed that calmodulin can play a crucial role in the excitation-contraction coupling. In fact, calmodulin (CaM) binds calcium and consequently regulates calcium channels. Several works show how some polymorphic variants can be considered predisposing factors to complex pathologies. Therefore, we hypothesize that the identification of polymorphic variants of proteins involved in the CaM pathway could be important for understanding how genetic traits can influence predisposition to myocardial infarction. This review considers each pathway of the three different isoforms of calmodulin (CaM1; CaM2; CaM3) and focuses on some common proteins involved in the three pathways, with the aim of analyzing the polymorphisms studied in the literature and understanding if they are associated with cardiovascular disease.

Calcium signaling is involved in diverse cellular processes in fungi

Calcium (Ca²⁺) is a universal signalling molecule of life. The Ca²⁺ signalling is an evolutionarily conserved process from prokaryotes to eukaryotes. Ca²⁺ at high concentration is deleterious to the cell; therefore, cell maintains a low resting level of intracellular free Ca²⁺ concentration ([Ca²⁺]_c). The resting [Ca²⁺]_c is tightly regulated, and a transient increase of the [Ca²⁺]_c initiates a signalling cascade in the cell. Ca²⁺ signalling plays an essential role in various processes, including growth, development, reproduction, tolerance to stress conditions, and virulence in fungi. In this review, we describe the evolutionary aspects of Ca²⁺ signalling and cell functions of major Ca²⁺ signalling proteins in different fungi.

PIRT the TRP Channel Regulating Protein Binds Calmodulin and Cholesterol-Like Ligands

Transient receptor potential (TRP) ion channels are polymodal receptors that have been implicated in a variety of pathophysiologies, including pain, obesity, and cancer. The capsaicin and heat sensor TRPV1, and the menthol and cold sensor TRPM8, have been shown to be modulated by the membrane protein PIRT (Phosphoinositide-interacting regulator of TRP). The emerging mechanism of PIRT-dependent TRPM8 regulation involves a competitive interaction between PIRT and TRPM8 for the activating phosphatidylinositol 4,5-bisphosphate (PIP2) lipid. As many PIP2 modulated ion channels also interact with calmodulin, we investigated the possible interaction between PIRT and calmodulin. Using microscale thermophoresis (MST), we show that calmodulin binds to the PIRT C-terminal α-helix, which we corroborate with a pull-down experiment, nuclear magnetic resonance-detected binding study, and Rosetta-based computational studies. Furthermore, we identify a cholesterol-recognition amino acid consensus (CRAC) domain in the outer leaflet of the first transmembrane helix of PIRT, and with MST, show that PIRT specifically binds to a number of cholesterol-derivatives. Additional studies identified that PIRT binds to cholecalciferol and oxytocin, which has mechanistic implications for the role of PIRT regulation of additional ion channels. This is the first study to show that PIRT specifically binds to a variety of ligands beyond TRP channels and PIP2.

A native prokaryotic voltage-dependent calcium channel with a novel selectivity filter sequence

Voltage-dependent Ca²⁺ channels (Cavs) are indispensable for coupling action potentials with Ca²⁺ signaling in living organisms. The structure of Cavs is similar to that of voltage-dependent Na⁺ channels (Navs). It is known that prokaryotic Navs can obtain Ca²⁺ selectivity by negative charge mutations of the selectivity filter, but native prokaryotic Cavs had not yet been identified. We report the first identification of a native prokaryotic Cav, CavMr, whose selectivity filter contains a smaller number of negatively charged residues than that of artificial prokaryotic Cavs. A relative mutant whose selectivity filter was replaced with that of CavMr exhibits high Ca²⁺ selectivity. Mutational analyses revealed that the glycine residue of the CavMr selectivity filter is a determinant for Ca²⁺ selectivity. This glycine residue is well conserved among subdomains I and III of eukaryotic Cavs. These findings provide new insight into the Ca²⁺ selectivity mechanism that is conserved from prokaryotes to eukaryotes.

Electrical signals in the brain and muscles allow animals – including humans – to think, make memories and move around. Cells generate these signals by enabling charged particles known as ions to pass through the physical barrier that surrounds all cells, the cell membrane, at certain times and in certain locations.

The ions pass through pores made by various channel proteins, which generally have so-called “selectivity filters” that only allow particular types of ions to fit through. For example, the selectivity filters of a family of channels in mammals known as the Cavs only allow calcium ions to pass through. Another family of ion channels in mammals are similar in structure to the Cavs but their selectivity filters only allow sodium ions to pass through instead of calcium ions.

Ion channels are found in all living cells including in bacteria. It is thought that the Cavs and sodium-selective channels may have both evolved from Cav-like channels in an ancient lifeform that was the common ancestor of modern bacteria and animals. Previous studies in bacteria found that modifying the selectivity filters of some sodium-selective channels known as BacNavs allowed calcium ions to pass through the mutant channels instead of sodium ions. However, no Cav channels had been identified in bacteria so far, representing a missing link in the evolutionary history of ion channels.

Shimomura et al. have now found a Cav-like channel in a bacterium known as Meiothermus ruber. Like all proteins, ion channels are made from amino acids and comparing the selectivity filter of the M. ruber Cav with those of mammalian Cavs and the calcium-selective BacNav mutants from previous studies revealed one amino acid that plays a particularly important role. This amino acid is a glycine that helps select which ions may pass through the pore and is also present in the selectivity filters of many Cavs in mammals.

Together these findings suggest that the Cav channel from M. ruber is similar to the mammal Cav channels and may more closely resemble the Cav-like channels thought to have existed in the common ancestor of bacteria and animals. Since other channel proteins from bacteria are useful genetic tools for studies in human and other animal cells, the Cav channel from M. ruber has the potential to be used to stimulate calcium signaling in experiments.

Calmodulin and Calmodulin Binding Proteins in Dictyostelium: A Primer

Dictyostelium discoideum is gaining increasing attention as a model organism for the study of calcium binding and calmodulin function in basic biological events as well as human diseases. After a short overview of calcium-binding proteins, the structure of Dictyostelium calmodulin and the conformational changes effected by calcium ion binding to its four EF hands are compared to its human counterpart, emphasizing the highly conserved nature of this central regulatory protein. The calcium-dependent and -independent motifs involved in calmodulin binding to target proteins are discussed with examples of the diversity of calmodulin binding proteins that have been studied in this amoebozoan. The methods used to identify and characterize calmodulin binding proteins is covered followed by the ways Dictyostelium is currently being used as a system to study several neurodegenerative diseases and how it could serve as a model for studying calmodulinopathies such as those associated with specific types of heart arrythmia. Because of its rapid developmental cycles, its genetic tractability, and a richly endowed stock center, Dictyostelium is in a position to become a leader in the field of calmodulin research.

Multi-state design of flexible proteins predicts sequences optimal for conformational change

Computational protein design of an ensemble of conformations for one protein–i.e., multi-state design–determines the side chain identity by optimizing the energetic contributions of that side chain in each of the backbone conformations. Sampling the resulting large sequence-structure search space limits the number of conformations and the size of proteins in multi-state design algorithms. Here, we demonstrated that the REstrained CONvergence (RECON) algorithm can simultaneously evaluate the sequence of large proteins that undergo substantial conformational changes. Simultaneous optimization of side chain conformations across all conformations increased sequence conservation when compared to single-state designs in all cases. More importantly, the sequence space sampled by RECON MSD resembled the evolutionary sequence space of flexible proteins, particularly when confined to predicting the mutational preferences of limited common ancestral descent, such as in the case of influenza type A hemagglutinin. Additionally, we found that sequence positions which require substantial changes in their local environment across an ensemble of conformations are more likely to be conserved. These increased conservation rates are better captured by RECON MSD over multiple conformations and thus multiple local residue environments during design. To quantify this rewiring of contacts at a certain position in sequence and structure, we introduced a new metric designated ‘contact proximity deviation’ that enumerates contact map changes. This measure allows mapping of global conformational changes into local side chain proximity adjustments, a property not captured by traditional global similarity metrics such as RMSD or local similarity metrics such as changes in φ and ψ angles.

Multi-state design can be used to engineer proteins that need to exist in multiple conformations or that bind to multiple partner molecules. In essence, multi-state design selects a compromise of protein sequences that allow for an ensemble of protein conformations, or states, associated with a particular biological function. In this paper, we used the REstrained CONvergence (RECON) algorithm with Rosetta to show that multi-state design of flexible proteins predicts sequences optimal for conformational change, mimicking mutation preferences sampled in evolution. Modeling optimal local side chain physicochemical environments within an ensemble selected significantly more native-like sequences than selections performed when all conformations states are designed independently. This outcome was particularly true for amino acids whose local side chain environment change between conformations. To quantify such contact map changes, we introduced a novel metric to show that sequence conservation is dependent on protein flexibility, i.e., changes in local side chain environments between stated limit the space of tolerated mutations. Additionally, such positions in sequence and structure are more likely to be energetically frustrated, at least in some states. Importantly, we showed that multi-state design over an ensemble of conformations (space) can explore evolutionary tolerated sequence space (time), thus enabling RECON to not only design proteins that require multiple states for function but also predict mutations that might be tolerated in native proteins but have not yet been explored by evolution. The latter aspect can be important to anticipate escape mutations, for example in pathogens or oncoproteins.

Molecular Insights from Conformational Ensembles via Machine Learning

Biomolecular simulations are intrinsically high dimensional and generate noisy data sets of ever-increasing size. Extracting important features from the data is crucial for understanding the biophysical properties of molecular processes, but remains a big challenge. Machine learning (ML) provides powerful dimensionality reduction tools. However, such methods are often criticized as resembling black boxes with limited human-interpretable insight. We use methods from supervised and unsupervised ML to efficiently create interpretable maps of important features from molecular simulations. We benchmark the performance of several methods, including neural networks, random forests, and principal component analysis, using a toy model with properties reminiscent of macromolecular behavior. We then analyze three diverse biological processes: conformational changes within the soluble protein calmodulin, ligand binding to a G protein-coupled receptor, and activation of an ion channel voltage-sensor domain, unraveling features critical for signal transduction, ligand binding, and voltage sensing. This work demonstrates the usefulness of ML in understanding biomolecular states and demystifying complex simulations.

Two SnRK2-Interacting Calcium Sensor Isoforms Negatively Regulate SnRK2 Activity by Different Mechanisms

SNF1-related protein kinases 2 (SnRK2s) are key signaling elements regulating abscisic acid-dependent plant development and responses to environmental stresses. Our previous data showed that the SnRK2-interacting Calcium Sensor (SCS) inhibits SnRK2 activity. Use of alternative transcription start sites located within the Arabidopsis (Arabidopsis thaliana) AtSCS gene results in two in-frame transcripts and subsequently two proteins, that differ only by the sequence position of the N terminus. We previously described the longer AtSCS-A, and now describe the shorter AtSCS-B and compare the two isoforms. The two isoforms differ substantially in their expression profiles in plant organs and in response to environmental stresses, in their calcium binding properties, and in their conformational dynamics in the presence and absence of Ca²⁺ Only AtSCS-A has the features of a calcium sensor. Both forms inhibit SnRK2 activity, but while AtSCS-A requires calcium for inhibition, AtSCS-B does not. Analysis of Arabidopsis plants stably expressing 35S::AtSCS-A-c-myc or 35S::AtSCS-B-c-myc in the scs-1 knockout mutant background revealed that, in planta, both forms are negative regulators of abscisic acid-induced SnRK2 activity and regulate plant resistance against water deficit. Moreover, the data highlight biochemical, biophysical, and functional properties of EF-hand-like motifs in plant proteins.