Interactions between domains of apo calmodulin alter calcium binding and stability

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

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

Redox regulation of KV7 channels through EF3 hand of calmodulin

Neuronal KV7 channels, important regulators of cell excitability, are among the most sensitive proteins to reactive oxygen species. The S2S3 linker of the voltage sensor was reported as a site-mediating redox modulation of the channels. Recent structural insights reveal potential interactions between this linker and the Ca2+-binding loop of the third EF-hand of calmodulin (CaM), which embraces an antiparallel fork formed by the C-terminal helices A and B, constituting the calcium responsive domain (CRD). We found that precluding Ca2+ binding to the EF3 hand, but not to EF1, EF2, or EF4 hands, abolishes oxidation-induced enhancement of KV7.4 currents. Monitoring FRET (Fluorescence Resonance Energy Transfer) between helices A and B using purified CRDs tagged with fluorescent proteins, we observed that S2S3 peptides cause a reversal of the signal in the presence of Ca2+ but have no effect in the absence of this cation or if the peptide is oxidized. The capacity of loading EF3 with Ca2+ is essential for this reversal of the FRET signal, whereas the consequences of obliterating Ca2+ binding to EF1, EF2, or EF4 are negligible. Furthermore, we show that EF3 is critical for translating Ca2+ signals to reorient the AB fork. Our data are consistent with the proposal that oxidation of cysteine residues in the S2S3 loop relieves KV7 channels from a constitutive inhibition imposed by interactions between the EF3 hand of CaM which is crucial for this signaling.

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.

Novel CALM3 Variant Causing Calmodulinopathy With Variable Expressivity in a 4-Generation Family

BACKGROUND

CaM (calmodulin), encoded by 3 separate genes (CALM1, CALM2, and CALM3), is a multifunctional Ca²⁺-binding protein involved in many signal transduction events including ion channel regulation. CaM variants may present with early-onset long QT syndrome (LQTS), catecholaminergic polymorphic ventricular tachycardia, or sudden cardiac death. Most reported variants occurred de novo. We identified a novel CALM3 variant, p.Asn138Lys (N138K), in a 4-generation family segregating with LQTS. The aim of this study was to elucidate its pathogenicity and to compare it with that of p.D130G-CaM-a variant associated with a severe LQTS phenotype.

METHODS

We performed whole exome sequencing for a large, 4-generation family affected by LQTS. To assess the effect of the detected CALM3 variant, the intrinsic Ca²⁺-binding affinity was measured by stoichiometric Ca²⁺ titrations and equilibrium titrations. L-type Ca²⁺ and slow delayed rectifier potassium currents (I_CaL and I_Ks) were recorded by whole-cell patch-clamp. Cav1.2 and Kv7.1 membrane expression were determined by optical fluorescence assays.

RESULTS

We identified 14 p.N138K-CaM carriers in a family where 2 sudden deaths occurred in children. Several members were only mildly affected compared with CaM-LQTS patients to date described in literature. The intrinsic Ca²⁺-binding affinity of the CaM C-terminal domain was 10-fold lower for p.N138K-CaM compared with wild-type-CaM. I_CaL inactivation was slowed in cells expressing p.N138K-CaM but less than in p.D130G-CaM cells. Unexpectedly, a larger I_Ks current density was observed in cells expressing p.N138K-CaM, but not for p.D130G-CaM, compared with wild-type-CaM.

CONCLUSIONS

The p.N138K CALM3 variant impairs Ca²⁺-binding affinity of CaM and I_CaL inactivation but potentiates I_Ks. The variably expressed phenotype of this variant compared with previously published de novo LQTS-CaM variants is likely explained by a milder impairment of I_CaL inactivation combined with I_Ks augmentation.

NaV1.2 EFL domain allosterically enhances Ca2+ binding to sites I and II of WT and pathogenic calmodulin mutants bound to the channel CTD

Neuronal voltage-gated sodium channel NaV1.2 C-terminal domain (CTD) binds calmodulin (CaM) constitutively at its IQ motif. A solution structure (6BUT) and other NMR evidence showed that the CaM N domain (CaMN) is structurally independent of the C-domain (CaMC) whether CaM is bound to the NaV1.2IQp (1,901-1,927) or NaV1.2CTD (1,777-1,937) with or without calcium. However, in the CaM + NaV1.2CTD complex, the Ca2+ affinity of CaMN was more favorable than in free CaM, while Ca2+ affinity for CaMC was weaker than in the CaM + NaV1.2IQp complex. The CTD EF-like (EFL) domain allosterically widened the energetic gap between CaM domains. Cardiomyopathy-associated CaM mutants (N53I(N54I), D95V(D96V), A102V(A103V), E104A(E105A), D129G(D130G), and F141L(F142L)) all bound the NaV1.2 IQ motif favorably under resting (apo) conditions and bound calcium normally at CaMN sites. However, only N53I and A102V bound calcium at CaMC sites at [Ca2+] < 100 μM. Thus, they are expected to respond like wild-type CaM to Ca2+ spikes in excitable cells.

Ca2+-saturated calmodulin binds tightly to the N-terminal domain of A-type fibroblast growth factor homologous factors

Voltage-gated sodium channels (Na_vs) are tightly regulated by multiple conserved auxiliary proteins, including the four fibroblast growth factor homologous factors (FGFs), which bind the Na_v EF-hand like domain (EFL), and calmodulin (CaM), a multifunctional messenger protein that binds the Na_V IQ motif. The EFL domain and IQ motif are contiguous regions of Na_V cytosolic C-terminal domains (CTD), placing CaM and FGF in close proximity. However, whether the FGFs and CaM act independently, directly associate, or operate through allosteric interactions to regulate channel function is unknown. Titrations monitored by steady-state fluorescence spectroscopy, structural studies with solution NMR, and computational modeling demonstrated for the first time that both domains of (Ca²⁺)₄-CaM (but not apo CaM) directly bind two sites in the N-terminal domain (NTD) of A-type FGF splice variants (FGF11A, FGF12A, FGF13A, and FGF14A) with high affinity. The weaker of the (Ca²⁺)₄-CaM-binding sites was known via electrophysiology to have a role in long-term inactivation of the channel but not known to bind CaM. FGF12A binding to a complex of CaM associated with a fragment of the Na_V1.2 CTD increased the Ca²⁺-binding affinity of both CaM domains, consistent with (Ca²⁺)₄-CaM interacting preferentially with its higher-affinity site in the FGF12A NTD. Thus, A-type FGFs can compete with Na_V IQ motifs for (Ca²⁺)₄-CaM. During spikes in the cytosolic Ca²⁺ concentration that accompany an action potential, CaM may translocate from the Na_V IQ motif to the FGF NTD, or the A-type FGF NTD may recruit a second molecule of CaM to the channel.

The arrhythmogenic N53I variant subtly changes the structure and dynamics in the calmodulin N-terminal domain, altering its interaction with the cardiac ryanodine receptor

Mutations in the genes encoding the highly conserved Ca²⁺-sensing protein calmodulin (CaM) cause severe cardiac arrhythmias, including catecholaminergic polymorphic ventricular tachycardia or long QT syndrome and sudden cardiac death. Most of the identified arrhythmogenic mutations reside in the C-terminal domain of CaM and mostly affect Ca²⁺-coordinating residues. One exception is the catecholaminergic polymorphic ventricular tachycardia-causing N53I substitution, which resides in the N-terminal domain (N-domain). It does not affect Ca²⁺ coordination and has only a minor impact on binding affinity toward Ca²⁺ and on other biophysical properties. Nevertheless, the N53I substitution dramatically affects CaM's ability to reduce the open probability of the cardiac ryanodine receptor (RyR2) while having no effect on the regulation of the plasmalemmal voltage-gated Ca²⁺ channel, Ca_v1.2. To gain more insight into the molecular disease mechanism of this mutant, we used NMR to investigate the structures and dynamics of both apo- and Ca²⁺-bound CaM-N53I in solution. We also solved the crystal structures of WT and N53I CaM in complex with the primary calmodulin-binding domain (CaMBD2) from RyR2 at 1.84-2.13 Å resolutions. We found that all structures of the arrhythmogenic CaM-N53I variant are highly similar to those of WT CaM. However, we noted that the N53I substitution exposes an additional hydrophobic surface and that the intramolecular dynamics of the protein are significantly altered such that they destabilize the CaM N-domain. We conclude that the N53I-induced changes alter the interaction of the CaM N-domain with RyR2 and thereby likely cause the arrhythmogenic phenotype of this mutation.

Calmodulin Mutations Associated with Heart Arrhythmia: A Status Report

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

Arrhythmia mutations in calmodulin can disrupt cooperativity of Ca2+ binding and cause misfolding

KEY POINTS

Mutations in the calmodulin protein (CaM) are associated with arrhythmia syndromes. This study focuses on understanding the structural characteristics of CaM disease mutants and their interactions with the voltage-gated calcium channel Ca_V 1.2. Arrhythmia mutations in CaM can lead to loss of Ca²⁺ binding, uncoupling of Ca²⁺ binding cooperativity, misfolding of the EF-hands and altered affinity for the calcium channel. These results help us to understand how different CaM mutants have distinct effects on structure and interactions with protein targets to cause disease.

ABSTRACT

Calmodulinopathies are life-threatening arrhythmia syndromes that arise from mutations in calmodulin (CaM), a calcium sensing protein whose sequence is completely conserved across all vertebrates. These mutations have been shown to interfere with the function of cardiac ion channels, including the voltage-gated Ca²⁺ channel Ca_V 1.2 and the ryanodine receptor (RyR2), in a mutation-specific manner. The ability of different CaM disease mutations to discriminate between these channels has been enigmatic. We present crystal structures of several C-terminal lobe mutants and an N-terminal lobe mutant in complex with the Ca_V 1.2 IQ domain, in conjunction with binding assays and complementary structural biology techniques. One mutation (D130G) causes a pathological conformation, with complete separation of EF-hands within the C-lobe and loss of Ca²⁺ binding in EF-hand 4. Another variant (Q136P) has severely reduced affinity for the IQ domain, and shows changes in the CD spectra under Ca²⁺ -saturating conditions when unbound to the IQ domain. Ca²⁺ binding to a pair of EF-hands normally proceeds with very high cooperativity, but we found that N98S CaM can adopt different conformations with either one or two Ca²⁺ ions bound to the C-lobe, possibly disrupting the cooperativity. An N-lobe variant (N54I), which causes severe stress-induced arrhythmia, does not show any major changes in complex with the IQ domain, providing a structural basis for why this mutant does not affect function of Ca_V 1.2. These findings show that different CaM mutants have distinct effects on both the CaM structure and interactions with protein targets, and act via distinct pathological mechanisms to cause disease.

Composite Conformational Changes of Signaling Proteins upon Ligand Binding Revealed by a Single Approach: Calcium-Calmodulin Study

Signaling proteins exemplified by calmodulin usually bind cooperatively to multiple ligands. Intermediate states and allosteric behavior are difficult to characterize. Here we extend a recently reported mass spectrometry (MS)-based method named LITPOMS (ligand titration, fast photochemical oxidation of proteins and mass spectrometry) that characterizes complex binding systems typically found as signaling proteins. As reported previously, calmodulin's response to binding four Ca²⁺ can be determined by LITPOMS to reveal binding sites, binding order, and most importantly composite binding behavior. Modeling this behavior provides site-specific binding affinities. In this article, we dissect the composite, peptide-level conformational changes at several regions either by digestion with a different protease or by tandem MS of LITPOMS behavior at the amino-acid residue level. Such dissection greatly elevates spatial resolution and increases the confidence of binding-order assignment. These complementary views of complex protein conformational change recapitulate the cumulative understanding via a single approach, providing new insights on poorly understood yet important allostery and underpin an approach applicable for exploring other signaling systems.

The CaMKII inhibitor KN93-calmodulin interaction and implications for calmodulin tuning of NaV1.5 and RyR2 function

Here we report the structure of the widely utilized calmodulin (CaM)-dependent protein kinase II (CaMKII) inhibitor KN93 bound to the Ca²⁺-sensing protein CaM. KN93 is widely believed to inhibit CaMKII by binding to the kinase. The CaM-KN93 interaction is significant as it can interfere with the interaction between CaM and it's physiological targets, thereby raising the possibility of ascribing modified protein function to CaMKII phosphorylation while concealing a CaM-protein interaction. NMR spectroscopy, stopped-flow kinetic measurements, and x-ray crystallography were used to characterize the structure and biophysical properties of the CaM-KN93 interaction. We then investigated the functional properties of the cardiac Na⁺ channel (Na_V1.5) and ryanodine receptor (RyR2). We find that KN93 disrupts a high affinity CaM-Na_V1.5 interaction and alters channel function independent of CaMKII. Moreover, KN93 increases RyR2 Ca²⁺ release in cardiomyocytes independent of CaMKII. Therefore, when interpreting KN93 data, targets other than CaMKII need to be considered.

Structural and dynamic responses of calcium ion binding loop residues in metallo-proteins

Conformational changes in bio-molecular systems are fundamental to several biological processes. It is important to study changes in responses of underlying microscopic variables, like dihedral angles as conformational change takes place. We perform all-atom simulations and modelling via Langevin equation to illustrate the changes in structural and dynamic responses of dihedral angles of calcium ion binding residues of different proteins in metal ion free (apo) and bound (holo) states. The equilibrium distributions of dihedral angles in apo- and holo-states represent structural response. Our studies show the presence of dihedrals with multiple peaks (isomeric states) separated by barrier heights is more frequent in apo- than in holo-state. The relaxation time-scale of dihedral fluctuations is found to increase linearly with decreasing barrier height due to more frequent barrier re-crossing events. The slow kinetic response of the dihedrals also contributes to slowing down of macro-scale fluctuations, which may be useful to understand kinetics of various bio-molecular processes.

A Single Approach Reveals the Composite Conformational Changes, Order of Binding, and Affinities for Calcium Binding to Calmodulin

We found that a newly developed method named LITPOMS (ligand titration, fast photochemical oxidation of proteins and mass spectrometry) can characterize section-by-section of a protein the conformational changes induced by metal-ion binding. Peptide-level LITPOMS applied to Ca²⁺ binding to calmodulin reveals binding order and site-specific affinity, providing new insights on the behavior of proteins upon binding Ca²⁺. We established that EF hand-4 (EF-4) binds calcium first, followed by EF-3, EF-2, and EF-1 and determined the four affinity constants by modeling the extent-of-modification curves. We also found positive cooperativity between EF-4, EF-3 and EF-2, EF-1 and allostery involving the four EF-hands. LITPOMS recapitulates via one approach the calcium-calmodulin binding that required decades of sophisticated development to afford versatility, comprehensiveness, and outstanding spatial resolution.

Ca2+-dependent calmodulin binding to cardiac ryanodine receptor (RyR2) calmodulin-binding domains

The Ca²⁺ sensor calmodulin (CaM) regulates cardiac ryanodine receptor (RyR2)-mediated Ca²⁺ release from the sarcoplasmic reticulum. CaM inhibits RyR2 in a Ca²⁺-dependent manner and aberrant CaM-dependent inhibition results in life-threatening cardiac arrhythmias. However, the molecular details of the CaM–RyR2 interaction remain unclear. Four CaM-binding domains (CaMBD1a, -1b, -2, and -3) in RyR2 have been proposed. Here, we investigated the Ca²⁺-dependent interactions between CaM and these CaMBDs by monitoring changes in the fluorescence anisotropy of carboxytetramethylrhodamine (TAMRA)-labeled CaMBD peptides during titration with CaM at a wide range of Ca²⁺ concentrations. We showed that CaM bound to all four CaMBDs with affinities that increased with Ca²⁺ concentration. CaM bound to CaMBD2 and -3 with high affinities across all Ca²⁺ concentrations tested, but bound to CaMBD1a and -1b only at Ca²⁺ concentrations above 0.2 µM. Binding experiments using individual CaM domains revealed that the CaM C-domain preferentially bound to CaMBD2, and the N-domain to CaMBD3. Moreover, the Ca²⁺ affinity of the CaM C-domain in complex with CaMBD2 or -3 was so high that these complexes are essentially Ca²⁺ saturated under resting Ca²⁺ conditions. Conversely, the N-domain senses Ca²⁺ exactly in the transition from resting to activating Ca²⁺ when complexed to either CaMBD2 or -3. Altogether, our results support a binding model where the CaM C-domain is anchored to RyR2 CaMBD2 and saturated with Ca²⁺ during Ca²⁺ oscillations, while the CaM N-domain functions as a dynamic Ca²⁺ sensor that can bridge noncontiguous regions of RyR2 or clamp down onto CaMBD2.

Backbone resonance assignments of complexes of apo human calmodulin bound to IQ motif peptides of voltage-dependent sodium channels NaV1.1, NaV1.4 and NaV1.7

Human voltage-gated sodium (Na_V) channels are critical for initiating and propagating action potentials in excitable cells. Nine isoforms have different roles but similar topologies, with a pore-forming α-subunit and auxiliary transmembrane β-subunits. Na_V pathologies lead to debilitating conditions including epilepsy, chronic pain, cardiac arrhythmias, and skeletal muscle paralysis. The ubiquitous calcium sensor calmodulin (CaM) binds to an IQ motif in the C-terminal tail of the α-subunit of all Na_V isoforms, and contributes to calcium-dependent pore-gating in some channels. Previous structural studies of calcium-free (apo) CaM bound to the IQ motifs of Na_V1.2, Na_V1.5, and Na_V1.6 showed that CaM binding was mediated by the C-domain of CaM (CaM_C), while the N-domain (CaM_N) made no detectable contacts. To determine whether this domain-specific recognition mechanism is conserved in other Na_V isoforms, we used solution NMR spectroscopy to assign the backbone resonances of complexes of apo CaM bound to peptides of IQ motifs of Na_V1.1, Na_V1.4, and Na_V1.7. Analysis of chemical shift differences showed that peptide binding only perturbed resonances in CaM_C; resonances of CaM_N were identical to free CaM. Thus, CaM_C residues contribute to the interface with the IQ motif, while CaM_N is available to interact elsewhere on the channel.

Defining potential roles of Pb(2+) in neurotoxicity from a calciomics approach

Metal ions play crucial roles in numerous biological processes, facilitating biochemical reactions by binding to various proteins. An increasing body of evidence suggests that neurotoxicity associated with exposure to nonessential metals (e.g., Pb(2+)) involves disruption of synaptic activity, and these observed effects are associated with the ability of Pb(2+) to interfere with Zn(2+) and Ca(2+)-dependent functions. However, the molecular mechanism behind Pb(2+) toxicity remains a topic of debate. In this review, we first discuss potential neuronal Ca(2+) binding protein (CaBP) targets for Pb(2+) such as calmodulin (CaM), synaptotagmin, neuronal calcium sensor-1 (NCS-1), N-methyl-d-aspartate receptor (NMDAR) and family C of G-protein coupled receptors (cGPCRs), and their involvement in Ca(2+)-signalling pathways. We then compare metal binding properties between Ca(2+) and Pb(2+) to understand the structural implications of Pb(2+) binding to CaBPs. Statistical and biophysical studies (e.g., NMR and fluorescence spectroscopy) of Pb(2+) binding are discussed to investigate the molecular mechanism behind Pb(2+) toxicity. These studies identify an opportunistic, allosteric binding of Pb(2+) to CaM, which is distinct from ionic displacement. Together, these data suggest three potential modes of Pb(2+) activity related to molecular and/or neural toxicity: (i) Pb(2+) can occupy Ca(2+)-binding sites, inhibiting the activity of the protein by structural modulation, (ii) Pb(2+) can mimic Ca(2+) in the binding sites, falsely activating the protein and perturbing downstream activities, or (iii) Pb(2+) can bind outside of the Ca(2+)-binding sites, resulting in the allosteric modulation of the protein activity. Moreover, the data further suggest that even low concentrations of Pb(2+) can interfere at multiple points within the neuronal Ca(2+) signalling pathways to cause neurotoxicity.

A cryptophane-based "turn-on" 129Xe NMR biosensor for monitoring calmodulin

We present the first cryptophane-based "turn-on" 129Xe NMR biosensor, employing a peptide-functionalized cryptophane to monitor the activation of calmodulin (CaM) protein in solution. In the absence of CaM binding, interaction between the peptide and cryptophane completely suppresses the hyperpolarized 129Xe-cryptophane NMR signal. Biosensor binding to Ca2+-activated CaM produces the expected 129Xe-cryptophane NMR signal.

Backbone resonance assignments of complexes of human voltage-dependent sodium channel NaV1.2 IQ motif peptide bound to apo calmodulin and to the C-domain fragment of apo calmodulin

Human voltage-gated sodium channel Na_V1.2 has a single pore-forming α-subunit and two transmembrane β-subunits. Expressed primarily in the brain, Na_V1.2 is critical for initiation and propagation of action potentials. Milliseconds after the pore opens, sodium influx is terminated by inactivation processes mediated by regulatory proteins including calmodulin (CaM). Both calcium-free (apo) CaM and calcium-saturated CaM bind tightly to an IQ motif in the C-terminal tail of the α-subunit. Our thermodynamic studies and solution structure (2KXW) of a C-domain fragment of apo ¹³C,¹⁵N- CaM (CaM_C) bound to an unlabeled peptide with the sequence of rat Na_V1.2 IQ motif showed that apo CaM_C (a) was necessary and sufficient for binding, and (b) bound more favorably than calcium-saturated CaM_C. However, we could not monitor the Na_V1.2 residues directly, and no structure of full-length CaM (including the N-domain of CaM (CaM_N)) was determined. To distinguish contributions of CaM_N and CaM_C, we used solution NMR spectroscopy to assign the backbone resonances of a complex containing a ¹³C,¹⁵N-labeled peptide with the sequence of human Na_V1.2 IQ motif (Na_V1.2_IQp) bound to apo ¹³C,¹⁵N-CaM or apo ¹³C,¹⁵N-CaM_C. Comparing the assignments of apo CaM in complex with Na_V1.2_IQp to those of free apo CaM showed that residues within CaM_C were significantly perturbed, while residues within CaM_N were essentially unchanged. The chemical shifts of residues in Na_V1.2_IQp and in the C-domain of CaM were nearly identical regardless of whether CaM_N was covalently linked to CaM_C. This suggests that CaM_N does not influence apo CaM binding to Na_V1.2_IQp.

The Recognition of Calmodulin to the Target Sequence of Calcineurin-A Novel Binding Mode

Calcineurin (CaN) is a Ca2+/calmodulin-dependent Ser/Thr protein phosphatase, which plays essential roles in many cellular and developmental processes. CaN comprises two subunits, a catalytic subunit (CaN-A, 60 kDa) and a regulatory subunit (CaN-B, 19 kDa). CaN-A tightly binds to CaN-B in the presence of minimal levels of Ca2+, but the enzyme is inactive until activated by CaM. Upon binding to CaM, CaN then undergoes a conformational rearrangement, the auto inhibitory domain is displaced and thus allows for full activity. In order to elucidate the regulatory role of CaM in the activation processes of CaN, we used NMR spectroscopy to determine the structure of the complex of CaM and the target peptide of CaN (CaNp). The CaM/CaNp complex shows a compact ellipsoidal shape with 8 α-helices of CaM wrapping around the CaNp helix. The RMSD of backbone and heavy atoms of twenty lowest energy structures of CaM/CaNp complex are 0.66 and 1.14 Å, respectively. The structure of CaM/CaNp complex can be classified as a novel binding mode family 1-18 with major anchor residues Ile396 and Leu413 to allocate the largest space between two domains of CaM. The relative orientation of CaNp to CaM is similar to the CaMKK peptide in the 1-16 binding mode with N- and C-terminal hydrophobic anchors of target sequence engulfed in the hydrophobic pockets of the N- and C-domain of CaM, respectively. In the light of the structural model of CaM/CaNp complex reported here, we provide new insight in the activation processes of CaN by CaM. We propose that the hydrophobic interactions between the Ca2+-saturated C-domain and C-terminal half of the target sequence provide driving forces for the initial recognition. Subsequent folding in the target sequence and structural readjustments in CaM enhance the formation of the complex and affinity to calcium. The electrostatic repulsion between CaM/CaNp complex and AID may result in the displacement of AID from active site for full activity.

Calcium triggers reversal of calmodulin on nested anti-parallel sites in the IQ motif of the neuronal voltage-dependent sodium channel NaV1.2

Several members of the voltage-gated sodium channel family are regulated by calmodulin (CaM) and ionic calcium. The neuronal voltage-gated sodium channel NaV1.2 contains binding sites for both apo (calcium-depleted) and calcium-saturated CaM. We have determined equilibrium dissociation constants for rat NaV1.2 IQ motif [IQRAYRRYLLK] binding to apo CaM (~3nM) and (Ca2+)4-CaM (~85nM), showing that apo CaM binding is favored by 30-fold. For both apo and (Ca2+)4-CaM, NMR demonstrated that NaV1.2 IQ motif peptide (NaV1.2IQp) exclusively made contacts with C-domain residues of CaM (CaMC). To understand how calcium triggers conformational change at the CaM-IQ interface, we determined a solution structure (2M5E.pdb) of (Ca2+)2-CaMC bound to NaV1.2IQp. The polarity of (Ca2+)2-CaMC relative to the IQ motif was opposite to that seen in apo CaMC-Nav1.2IQp (2KXW), revealing that CaMC recognizes nested, anti-parallel sites in Nav1.2IQp. Reversal of CaM may require transient release from the IQ motif during calcium binding, and facilitate a re-orientation of CaMN allowing interactions with non-IQ NaV1.2 residues or auxiliary regulatory proteins interacting in the vicinity of the IQ motif.

The effects of thioamide backbone substitution on protein stability: a study in α-helical, β-sheet, and polyproline II helical contexts

Thioamides are single atom substitutions of the peptide bond that serve as versatile probes of protein structure. Effective use of thioamides requires a robust understanding of the impact that the substitution has on a protein of interest. However, the thermodynamic effects of thioamide incorporation have only been studied in small structural motifs, and their influence on secondary structure in the context of full-length proteins is not known. Here we describe a comprehensive survey of thioamide substitutions in three benchmark protein systems (calmodulin, the B1 domain of protein G, and collagen) featuring the most prevalent secondary structure motifs: α-helix, β-sheet, and polyproline type II helix. We find that in most cases, effects on thermostability can be understood in terms of the positioning and local environment of the thioamide relative to proximal structural elements and hydrogen bonding networks. These observations set the stage for the rational design of thioamide substituted proteins with predictable stabilities.

Comparing allosteric transitions in the domains of calmodulin through coarse-grained simulations

Calmodulin (CaM) is a ubiquitous Ca(2+)-binding protein consisting of two structurally similar domains with distinct stabilities, binding affinities, and flexibilities. We present coarse grained simulations that suggest that the mechanism for the domain's allosteric transitions between the open and closed conformations depends on subtle differences in the folded state topology of the two domains. Throughout a wide temperature range, the simulated transition mechanism of the N-terminal domain (nCaM) follows a two-state transition mechanism while domain opening in the C-terminal domain (cCaM) involves unfolding and refolding of the tertiary structure. The appearance of the unfolded intermediate occurs at a higher temperature in nCaM than it does in cCaM consistent with nCaM's higher thermal stability. Under approximate physiological conditions, the simulated unfolded state population of cCaM accounts for 10% of the population with nearly all of the sampled transitions (approximately 95%) unfolding and refolding during the conformational change. Transient unfolding significantly slows the domain opening and closing rates of cCaM, which can potentially influence its Ca(2+)-binding mechanism.

Relative Cosolute Size Influences the Kinetics of Protein-Protein Interactions

Protein signaling occurs in crowded intracellular environments, and while high concentrations of macromolecules are postulated to modulate protein-protein interactions, analysis of their impact at each step of the reaction pathway has not been systematically addressed. Potential cosolute-induced alterations in target association are particularly important for a signaling molecule like calmodulin (CaM), where competition among >300 targets governs which pathways are selectively activated. To explore how high concentrations of cosolutes influence CaM-target affinity and kinetics, we methodically investigated each step of the CaM-target binding mechanism under crowded or osmolyte-rich environments mimicked by ficoll-70, dextran-10, and sucrose. All cosolutes stabilized compact conformers of CaM and modulated association kinetics by affecting diffusion and rates of conformational change; however, the results showed that differently sized molecules had variable effects to enhance or impede unique steps of the association pathway. On- and off-rates were modulated by all cosolutes in a compensatory fashion, producing little change in steady-state affinity. From this work insights were gained on how high concentrations of inert crowding agents and osmolytes fit into a kinetic framework to describe protein-protein interactions relevant for cellular signaling.

Biochemical and biophysical characterization of a plant calmodulin: Role of the N- and C-lobes in calcium binding, conformational change, and target interaction

In plants, transient elevation of intracellular Ca(2+) concentration in response to abiotic stress is responsible for glutamate decarboxylase (GAD) activation via association with calmodulin (CaM), an EF-hand protein consisting of two homologous domains (N and C). An unusual 1:2 binding mode of CaM to CaM-binding domains of GAD has long been known, however the contribution of the two CaM domains in target recognition and activation remains to be clarified. Here, we explored the coupling between physicochemical properties of Arabidopsis CaM1 (AtCaM1) and Arabidopsis GAD1 activation, focusing on each AtCaM1 lobe. We found that the four EF-loops of AtCaM1 differently contribute to the ~20 μM apparent affinity for Ca(2+) and the C-lobe shows a ~6-fold higher affinity than N-lobe (Kd(app) 5.6 μM and 32 μM for C- and N-lobes, respectively). AtCaM1 responds structurally to Ca(2+) in a manner similar to vertebrate CaM based on comparison of Ca(2+)-induced changes in hydrophobicity exposure, secondary structure, and hydrodynamic behavior. Molecular dynamics simulations of AtCaM1 apo and Ca(2+)-bound reveal that the latter state is significantly less flexible, although regions of the N-lobe remain quite flexible; this suggests the importance of N-lobe for completing the transition to the extended structure of holoprotein, consistent with data from ANS fluorescence, CD spectroscopy, and SEC analysis. Moreover, enzymatic analysis reveal that mutations in the two lobes affect GAD1 activation in similar ways and only intact AtCaM1 can fully activate GAD1. Taken together, our data provide new insights into the CaM lobes role in interactions between CaM and plant GAD.

Insights into molecular interactions between CaM and its inhibitors from molecular dynamics simulations and experimental data

In order to contribute to the structural basis for rational design of calmodulin (CaM) inhibitors, we analyzed the interaction of CaM with 14 classic antagonists and two compounds that do not affect CaM, using docking and molecular dynamics (MD) simulations, and the data were compared to available experimental data. The Ca(2+)-CaM-Ligands complexes were simulated 20 ns, with CaM starting in the "open" and "closed" conformations. The analysis of the MD simulations provided insight into the conformational changes undergone by CaM during its interaction with these ligands. These simulations were used to predict the binding free energies (ΔG) from contributions ΔH and ΔS, giving useful information about CaM ligand binding thermodynamics. The ΔG predicted for the CaM's inhibitors correlated well with available experimental data as the r(2) obtained was 0.76 and 0.82 for the group of xanthones. Additionally, valuable information is presented here: I) CaM has two preferred ligand binding sites in the open conformation known as site 1 and 4, II) CaM can bind ligands of diverse structural nature, III) the flexibility of CaM is reduced by the union of its ligands, leading to a reduction in the Ca(2+)-CaM entropy, IV) enthalpy dominates the molecular recognition process in the system Ca(2+)-CaM-Ligand, and V) the ligands making more extensive contact with the protein have higher affinity for Ca(2+)-CaM. Despite their limitations, docking and MD simulations in combination with experimental data continue to be excellent tools for research in pharmacology, toward a rational design of new drugs.

Opposing orientations of the anti-psychotic drug trifluoperazine selected by alternate conformations of M144 in calmodulin

The anti-psychotic drug trifluoperazine (TFP) is an antagonist observed to bind to calcium-saturated calmodulin ((Ca(2+) )4 -CaM) at ratios of 1:1 (1CTR), 2:1 (1A29), and 4:1 (1LIN). Each structure contains one TFP bound in the hydrophobic cleft of the C-domain of CaM. However, the orientation of the trifluoromethyl (CF3 ) moiety differs among them: it is buried in the C-domain cleft of 1A29 and 1LIN, but protrudes from 1CTR. We report a 2.0 Å resolution crystallographic structure (4RJD) of TFP bound to the (Ca(2+) )-saturated C-domain of CaM (CaMC ). The asymmetric unit contains two molecules of (Ca(2+) )2 -CaMC . Chain backbones were nearly identical, but the orientation of TFP in the cleft of Chain A matched 1A29/1LIN, while TFP bound to Chain B matched 1CTR. This was accommodated by a flip of the M144 sidechain and small changes in sidechains of M109 and M145. Docking simulations suggested that the rotamer conformation of M144 determined the orientation of TFP within the cleft of (Ca(2+) )2 -CaMC . Chains A and B show that the open cleft of (Ca(2+) )2 -CaMC is promiscuous in accepting TFP in reversed directions under the same crystallization conditions. Observing multiple orientations of an antagonist bound to a single protein highlights the challenge of designing highly specific pharmaceuticals, and may have importance for QSAR of other CF3 -containing drugs such as fluoxetine (anti-depressant) or efavirenz (reverse transcriptase inhibitor). This study emphasizes that a single structure of a complex represents an energetically accessible state, but does not necessarily show the full range of energetically equivalent states.

Calmodulin mutations causing catecholaminergic polymorphic ventricular tachycardia confer opposing functional and biophysical molecular changes

Calmodulin (CaM) is the central mediator of intracellular Ca(2+) signalling in cardiomyocytes, where it conveys the intricate Ca(2+) transients to the proteins controlling cardiac contraction. We recently linked two separate mutations in CaM (N53I and N97S) to dominantly inherited catecholaminergic polymorphic ventricular tachycardia (CPVT), an arrhythmic disorder in which exercise or acute emotion can lead to syncope and sudden cardiac death. Given the ubiquitous presence of CaM in all eukaryote cells, it is particular intriguing that carriers of either mutation show no additional symptoms. Here, we investigated the effects of the CaM CPVT mutations in a zebrafish animal model. Three-day-old embryos injected with either CaM mRNA showed no detectable pathologies or developmental abnormalities. However, embryos injected with CPVT CaM mRNA displayed increased heart rate compared to wild-type CaM mRNA under β-adrenergic stimulation, demonstrating a conserved dominant cardiac specific effect between zebrafish and human carriers of these mutations. Motivated by the highly similar physiological phenotypes, we compared the effects of the N53I and N97S mutations on the biophysical and functional properties of CaM. Surprisingly, the mutations have opposing effects on CaM C-lobe Ca(2+) binding affinity and kinetics, and changes to the CaM N-lobe Ca(2+) binding are minor and specific to the N53I mutation. Furthermore, both mutations induce differential perturbations to structure and stability towards unfolding. Our results suggest different molecular disease mechanisms for the CPVT (N53I and N97S mutations) and strongly support that cardiac contraction is the physiological process most sensitive to CaM integrity.

Calcium-dependent energetics of calmodulin domain interactions with regulatory regions of the Ryanodine Receptor Type 1 (RyR1)

Calmodulin (CaM) allosterically regulates the homo-tetrameric human Ryanodine Receptor Type 1 (hRyR1): apo CaM activates the channel, while (Ca(2+))4-CaM inhibits it. CaM-binding RyR1 residues 1975-1999 and 3614-3643 were proposed to allow CaM to bridge adjacent RyR1 subunits. Fluorescence anisotropy titrations monitored the binding of CaM and its domains to peptides encompassing hRyR(11975-1999) or hRyR1(3614-3643). Both CaM and its C-domain associated in a calcium-independent manner with hRyR1(3614-3643) while N-domain required calcium and bound ~250-fold more weakly. Association with hRyR1(11975-1999) was weak. Both hRyR1 peptides increased the calcium-binding affinity of both CaM domains, while maintaining differences between them. These energetics support the CaM C-domain association with hRyR1(3614-3643) at low calcium, positioning CaM to respond to calcium efflux. However, the CaM N-domain affinity for hRyR(11975-1999) alone was insufficient to support CaM bridging adjacent RyR1 subunits. Other proteins or elements of the hRyR1 structure must contribute to the energetics of CaM-mediated regulation.

Altered RyR2 regulation by the calmodulin F90L mutation associated with idiopathic ventricular fibrillation and early sudden cardiac death

Calmodulin (CaM) association with the cardiac muscle ryanodine receptor (RyR2) regulates excitation-contraction coupling. Defective CaM-RyR2 interaction is associated with heart failure. A novel CaM mutation (CaM(F90L)) was recently identified in a family with idiopathic ventricular fibrillation (IVF) and early onset sudden cardiac death. We report the first biochemical characterization of CaM(F90L). F90L confers a deleterious effect on protein stability. Ca(2+)-binding studies reveal reduced Ca(2+)-binding affinity and a loss of co-operativity. Moreover, CaM(F90L) displays reduced RyR2 interaction and defective modulation of [(3)H]ryanodine binding. Hence, dysregulation of RyR2-mediated Ca(2+) release via aberrant CaM(F90L)-RyR2 interaction is a potential mechanism that underlies familial IVF.

Novel calmodulin mutations associated with congenital arrhythmia susceptibility

BACKGROUND

Genetic predisposition to life-threatening cardiac arrhythmias such as congenital long-QT syndrome (LQTS) and catecholaminergic polymorphic ventricular tachycardia (CPVT) represent treatable causes of sudden cardiac death in young adults and children. Recently, mutations in calmodulin (CALM1, CALM2) have been associated with severe forms of LQTS and CPVT, with life-threatening arrhythmias occurring very early in life. Additional mutation-positive cases are needed to discern genotype-phenotype correlations associated with calmodulin mutations.

METHODS AND RESULTS

We used conventional and next-generation sequencing approaches, including exome analysis, in genotype-negative LQTS probands. We identified 5 novel de novo missense mutations in CALM2 in 3 subjects with LQTS (p.N98S, p.N98I, p.D134H) and 2 subjects with clinical features of both LQTS and CPVT (p.D132E, p.Q136P). Age of onset of major symptoms (syncope or cardiac arrest) ranged from 1 to 9 years. Three of 5 probands had cardiac arrest and 1 of these subjects did not survive. The clinical severity among subjects in this series was generally less than that originally reported for CALM1 and CALM2 associated with recurrent cardiac arrest during infancy. Four of 5 probands responded to β-blocker therapy, whereas 1 subject with mutation p.Q136P died suddenly during exertion despite this treatment. Mutations affect conserved residues located within Ca(2+)-binding loops III (p.N98S, p.N98I) or IV (p.D132E, p.D134H, p.Q136P) and caused reduced Ca(2+)-binding affinity.

CONCLUSIONS

CALM2 mutations can be associated with LQTS and with overlapping features of LQTS and CPVT.

Neurogranin alters the structure and calcium binding properties of calmodulin

Neurogranin (Ng) is a member of the IQ motif class of calmodulin (CaM)-binding proteins, and interactions with CaM are its only known biological function. In this report we demonstrate that the binding affinity of Ng for CaM is weakened by Ca(2+) but to a lesser extent (2-3-fold) than that previously suggested from qualitative observations. We also show that Ng induced a >10-fold decrease in the affinity of Ca(2+) binding to the C-terminal domain of CaM with an associated increase in the Ca(2+) dissociation rate. We also discovered a modest, but potentially important, increase in the cooperativity in Ca(2+) binding to the C-lobe of CaM in the presence of Ng, thus sharpening the threshold for the C-domain to become Ca(2+)-saturated. Domain mapping using synthetic peptides indicated that the IQ motif of Ng is a poor mimetic of the intact protein and that the acidic sequence just N-terminal to the IQ motif plays an important role in reproducing Ng-mediated decreases in the Ca(2+) binding affinity of CaM. Using NMR, full-length Ng was shown to make contacts largely with residues in the C-domain of CaM, although contacts were also detected in residues in the N-terminal domain. Together, our results can be consolidated into a model where Ng contacts residues in the N- and C-lobes of both apo- and Ca(2+)-bound CaM and that although Ca(2+) binding weakens Ng interactions with CaM, the most dramatic biochemical effect is the impact of Ng on Ca(2+) binding to the C-terminal lobe of CaM.

Importance of the interaction protein-protein of the CaM-PDE1A and CaM-MLCK complexes in the development of new anti-CaM drugs

Protein-protein interactions play central roles in physiological and pathological processes. The bases of the mechanisms of drug action are relevant to the discovery of new therapeutic targets. This work focuses on understanding the interactions in protein-protein-ligands complexes, using proteins calmodulin (CaM), human calcium/calmodulin-dependent 3',5'-cyclic nucleotide phosphodiesterase 1A active human (PDE1A), and myosin light chain kinase (MLCK) and ligands αII-spectrin peptide (αII-spec), and two inhibitors of CaM (chlorpromazine (CPZ) and malbrancheamide (MBC)). The interaction was monitored with a fluorescent biosensor of CaM (hCaM M124C-mBBr). The results showed changes in the affinity of CPZ and MBC depending on the CaM-protein complex under analysis. For the Ca(2+) -CaM, Ca(2+) -CaM-PDE1A, and Ca(2+) -CaM-MLCK complexes, CPZ apparent dissociation constants (Kds ) were 1.11, 0.28, and 0.55 μM, respectively; and for MBC Kds were 1.43, 1.10, and 0.61 μM, respectively. In competition experiments the addition of calmodulin binding peptide 1 (αII-spec) to Ca(2+) -hCaM M124C-mBBr quenched the fluorescence (Kd  = 2.55 ± 1.75 pM) and the later addition of MBC (up to 16 μM) did not affect the fluorescent signal. Instead, the additions of αII-spec to a preformed Ca(2+) -hCaM M124C-mBBr-MBC complex modified the fluorescent signal. However, MBC was able to displace the PDE1A and MLCK from its complex with Ca(2+) -CaM. In addition, docking studies were performed for all complexes with both ligands showing an excellent correlation with experimental data. These experiments may help to explain why in vivo many CaM drugs target prefer only a subset of the Ca(2+) -CaM regulated proteins and adds to the understanding of molecular interactions between protein complexes and small ligands.

Binding free-energy calculation of an ion-peptide complex by constrained dynamics

Binding free energy is the most important physical parameter that describes the binding affinity of a receptor-ligand complex. Conventionally, it was obtained based on the thermodynamic cycle or alchemical reaction. These strategies have been widely used, but they would be problematic if the receptors and/or ligands have large conformational changes during the binding processes. In this paper, we present a way to calculate the binding free energy: constrained dynamics along a fragmental and high-dimensional transition path. This method directly considers unbound states in the simulation. The application to the calmodulin loop-calcium complexes shows that it is practical and the calculated relative binding affinities are in good agreement with experimental results.

Metal toxicity and opportunistic binding of Pb(2+) in proteins

Lead toxicity is associated with various human diseases. While Ca(2+) binding proteins such as calmodulin (CaM) are often reported to be molecular targets for Pb(2+)-binding and lead toxicity, the effect of Pb(2+) on the Ca(2+)/CaM regulated biological activities cannot be described by the primary mechanism of ionic displacement (e.g., ionic mimicry). The focus of this study was to investigate the mechanism of lead toxicity through binding differences between Ca(2+) and Pb(2+) for CaM, an essential intracellular trigger protein with two EF-Hand Ca(2+)-binding sites in each of its two domains that regulates many molecular targets via Ca(2+)-induced conformational change. Fluorescence changes in phenylalanine indicated that Pb(2+) binds with 8-fold higher affinity than Ca(2+) in the N-terminal domain. Additionally, NMR chemical shift changes and an unusual biphasic response observed in tyrosine fluorescence associated with C-terminal domain sites EF-III and EF-IV suggest a single higher affinity Pb(2+)-binding site with a 3-fold higher affinity than Ca(2+), coupled with a second site exhibiting affinity nearly equivalent to that of the N-terminal domain sites. Our results further indicate that Pb(2+) displaces Ca(2+) only in the N-terminal domain, with minimal perturbation of the C-terminal domain, however significant structural/dynamic changes are observed in the trans-domain linker region which appear to be due to Pb(2+)-binding outside of the known calcium-binding sites. These data suggest that opportunistic Pb(2+)-binding in Ca(2+)/CaM has a profound impact on the conformation and dynamics of the essential molecular recognition sites of the central helix, and provides insight into the molecular toxicity of non-essential metal ions.

Relative stability of human centrins and its relationship to calcium binding

Centrins are calcium binding proteins that belong to the EF-hand superfamily with diverse biological functions. Herein we present the first systematic study that establishes the relative stability of related centrins via complementary biophysical techniques. Our results define the stepwise molecular behavior of human centrins by two-dimensional infrared (2D IR) correlation spectroscopy, the change in heat capacity and enthalpy of denaturation by differential scanning calorimetry, and the relative stability of the helical regions of centrins by circular dichroism. More importantly, 2D IR correlation spectroscopy provides unique information about the similarities and differences in dynamics between these related proteins. The thermally induced molecular behavior of human centrins can be used to predict biological target interactions that have a relative dependence on calcium affinity. This information is essential for understanding why certain isoforms may be used to rescue a phenotype and therefore also for explaining the different functions these proteins may have in vivo. Furthermore, this comparative approach can be applied to the study of recombinant therapeutic protein candidates for the treatment of disease states.

Calcium-dependent conformational transition of calmodulin determined by Fourier transform infrared spectroscopy

The Ca(2+)-induced conformational changes in calmodulin (CaM) were monitored by Fourier transform infrared spectroscopy (FT-IR) at different molar ratios of Ca(2+) to CaM. The results show that these changes occur in two distinctive transitions. The first transition involves significant changes in the overall secondary structure with a small gain in solvent accessibility, and is completed after the second Ca(2+) binds to both EF-hands of its C-terminal domain. The second transition is accompanied by CaM folding into a tighter, less hydrogen-exchangeable structure, and is completed by the addition of the fourth Ca(2+) to have four Ca(2+) per molecule. Particularly, α-helices in CaM-nCa(2+)(n=0, 1, 2) are less stable than those in CaM-nCa(2+)(n=3, 4).

Utilization of a calmodulin lysine methyltransferase co-expression system for the generation of a combinatorial library of post-translationally modified proteins

By successfully incorporating sequence diversity into proteins, combinatorial libraries have been a staple technology used in protein engineering, directed evolution, and synthetic biology for generating proteins with novel specificities and activities. However, these approaches mostly overlook the incorporations of post-translational modifications, which nature extensively uses for modulating protein activities in vivo. As an initial step of incorporating post-translational modifications into combinatorial libraries, we present a bacterial co-expression system, utilizing a recently characterized calmodulin methyltransferase (CaM KMT), to trimethylate a combinatorial library of the calmodulin central linker region. We show that this system is robust, with the successful over-expression and post-translational modification performed in Escherichia coli. Furthermore we show that trimethylation differentially affected the conformational dynamics of the protein upon the binding of calcium, and the thermal stability of the apoprotein. Collectively, these data support that when applied to an appropriately designed protein library scaffold, CaM KMT is able to produce a post-translationally modified library of protein sequences, thus providing a powerful tool for future protein library designs and constructions.

Binding orientation and specificity of calmodulin to rat olfactory cyclic nucleotide-gated ion channel

Calmodulin (CaM), the primary intracellular Ca(2+) receptor, regulates a large number of key enzymes and controls a wide spectrum of important biological responses. Recognition between CaM and its target sequence in rat olfactory cyclic nucleotide-gated ion channel (OLFp) was investigated by circular dichroism (CD), fluorescence, and NMR spectroscopy. Fluorescence data showed the OLFp tightly bound to CaM with a dissociation constant of 12 nM in a 1:1 stoichiometry. Far-UV CD data showed that approximately 60% of OLFp residues formed α-helical structures when associated with CaM. NMR data showed that most of the (15)N-(1)H HSQC cross-peaks of the (15)N-labeled CaM not only shifted but also split into two sets of peaks upon association with the OLFp. Our data indicated that the two distinct CaM/OLFp complexes existed simultaneously with stable structures that were not interexchangeable within the NMR time scale. In light of the palindromic sequence of OLFp (FQRIVRLVGVIRDW) for CaM targeting, we proposed that the helical OLFp with C2 symmetry may bind to CaM in two orientations. This hypothesis is supported by the observation that only one set of (15)N-(1)H HSQC cross-peaks of the (15)N-labeled CaM was detected upon association with OLFp-M13 chimeric peptide (OLFMp), a mutated OLFp lacking the palindromic feature. The binding specificity of OLFMp to CaM was restored when the palindromic feature was destroyed. Binding modes of CaM/OLFp and CaM/OLFMp simulated by molecular docking were in accord with their distinct patterns observed in HSQC spectra. Our studies suggest that the palindromic residues in OLFp are crucial for the orientation-specific recognition by CaM.

The effect of macromolecular crowding, ionic strength and calcium binding on calmodulin dynamics

The flexibility in the structure of calmodulin (CaM) allows its binding to over 300 target proteins in the cell. To investigate the structure-function relationship of CaM, we combined methods of computer simulation and experiments based on circular dichroism (CD) to investigate the structural characteristics of CaM that influence its target recognition in crowded cell-like conditions. We developed a unique multiscale solution of charges computed from quantum chemistry, together with protein reconstruction, coarse-grained molecular simulations, and statistical physics, to represent the charge distribution in the transition from apoCaM to holoCaM upon calcium binding. Computationally, we found that increased levels of macromolecular crowding, in addition to calcium binding and ionic strength typical of that found inside cells, can impact the conformation, helicity and the EF hand orientation of CaM. Because EF hand orientation impacts the affinity of calcium binding and the specificity of CaM's target selection, our results may provide unique insight into understanding the promiscuous behavior of calmodulin in target selection inside cells.

Proteins are workhorses for driving biological functions inside cells. Calmodulin (CaM) is a protein that can carry cellular signals by triggered conformational changes due to calcium binding that alters target binding. Interestingly, CaM is able to bind over 300 targets. One of the challenges in characterizing CaM's ability to bind multiple targets lies in that CaM is a flexible protein and its structure is easily modulated by the physicochemical changes in its surroundings, particularly inside a complex cellular milieu. In order to determine structure-function relationships of CaM, we employed a combined approach of experiments, computer simulations and statistical physics in the investigation of the effect of calcium-binding, salt concentration, and macromolecular crowding on CaM. The results revealed unique folding energy landscapes of CaM in the absence and presence of calcium ions and the structural implications of CaM are interpreted under cell-like conditions. Further, a large conformational change in CaM in response to environmental impacts, dictates the packing of local helices that may be critical to its function of target binding and recognition among vast target selections.

NMR, biophysical, and biochemical studies reveal the minimal Calmodulin binding domain of the HIV-1 matrix protein

Subcellular distribution of Calmodulin (CaM) in human immunodeficiency virus type-1 (HIV-1)-infected cells is distinct from that observed in uninfected cells. CaM has been shown to interact and co-localize with the HIV-1 Gag protein in infected cells. However, the precise molecular mechanism of this interaction is not known. Binding of Gag to CaM is dependent on calcium and is mediated by the N-terminal-myristoylated matrix (myr(+)MA) domain. We have recently shown that CaM binding induces a conformational change in the MA protein, triggering exposure of the myristate group. To unravel the molecular mechanism of CaM-MA interaction and to identify the minimal CaM binding domain of MA, we devised multiple approaches utilizing NMR, biochemical, and biophysical methods. Short peptides derived from the MA protein have been examined. Our data revealed that whereas peptides spanning residues 11-28 (MA-(11-28)) and 31-46 (MA-(31-46)) appear to bind preferentially to the C-terminal lobe of CaM, a peptide comprising residues 11-46 (MA-(11-46)) appears to engage both domains of CaM. Limited proteolysis data conducted on the MA-CaM complex yielded a MA peptide (residues 8-43) that is protected by CaM and resistant to proteolysis. MA-(8-43) binds to CaM with a very high affinity (dissociation constant = 25 nm) and in a manner that is similar to that observed for the full-length MA protein. The present findings provide new insights on how MA interacts with CaM that may ultimately help in identification of the functional role of CaM-Gag interactions in the HIV replication cycle.

Thermodynamic linkage between calmodulin domains binding calcium and contiguous sites in the C-terminal tail of Ca(V)1.2

Calmodulin (CaM) binding to the intracellular C-terminal tail (CTT) of the cardiac L-type Ca(2+) channel (Ca(V)1.2) regulates Ca(2+) entry by recognizing sites that contribute to negative feedback mechanisms for channel closing. CaM associates with Ca(V)1.2 under low resting [Ca(2+)], but is poised to change conformation and position when intracellular [Ca(2+)] rises. CaM binding Ca(2+), and the domains of CaM binding the CTT are linked thermodynamic functions. To better understand regulation, we determined the energetics of CaM domains binding to peptides representing pre-IQ sites A(1588), and C(1614) and the IQ motif studied as overlapping peptides IQ(1644) and IQ'(1650) as well as their effect on calcium binding. (Ca(2+))(4)-CaM bound to all four peptides very favorably (K(d)≤2 nM). Linkage analysis showed that IQ(1644-1670) bound with a K(d)~1 pM. In the pre-IQ region, (Ca(2+))(2)-N-domain bound preferentially to A(1588), while (Ca(2+))(2)-C-domain preferred C(1614). When bound to C(1614), calcium binding in the N-domain affected the tertiary conformation of the C-domain. Based on the thermodynamics, we propose a structural mechanism for calcium-dependent conformational change in which the linker between CTT sites A and C buckles to form an A-C hairpin that is bridged by calcium-saturated CaM.

Conformation of the calmodulin-binding domain of metabotropic glutamate receptor subtype 7 and its interaction with calmodulin

Calmodulin (CaM), a Ca(2+)-binding protein, is a well-known regulator of various cellular functions. One of the targets of CaM is metabotropic glutamate receptor 7 (mGluR7), which serves as a low-pass filter for glutamate in the pre-synaptic terminal to regulate neurotransmission. Surface plasmon resonance (SPR), circular dichroism (CD) spectroscopy and nuclear magnetic spectroscopy (NMR) were performed to study the structure of the peptides corresponding to the CaM-binding domain of mGluR7 and their interaction with CaM. Unlike well-known CaM-binding peptides, mGluR7 has a random coil structure even in the presence of trifluoroethanol. Moreover, NMR data suggested that the complex between Ca(2+)/CaM and the mGluR7 peptide has multiple conformations. The mGluR7 peptide has been found to interact with CaM even in the absence of Ca(2+), and the binding is directed toward the C-domain of apo-CaM rather than the N-domain. We propose a possible mechanism for the activation of mGluR7 by CaM. A pre-binding occurs between apo-CaM and mGluR7 in the resting state of cells. Then, the Ca(2+)/CaM-mGluR7 complex is formed once Ca(2+) influx occurs. The weak interaction at lower Ca(2+) concentrations is likely to bind CaM to mGluR7 for the fast complex formation in response to the elevation of Ca(2+) concentration.

Binding of calmodulin to the HIV-1 matrix protein triggers myristate exposure

Steady progress has been made in defining both the viral and cellular determinants of retroviral assembly and release. Although it is widely accepted that targeting of the Gag polypeptide to the plasma membrane is critical for proper assembly of HIV-1, the intracellular interactions and trafficking of Gag to its assembly sites in the infected cell are poorly understood. HIV-1 Gag was shown to interact and co-localize with calmodulin (CaM), a ubiquitous and highly conserved Ca(2+)-binding protein expressed in all eukaryotic cells, and is implicated in a variety of cellular functions. Binding of HIV-1 Gag to CaM is dependent on calcium and is mediated by the N-terminally myristoylated matrix (myr(+)MA) domain. Herein, we demonstrate that CaM binds to myr(+)MA with a dissociation constant (K(d)) of ∼2 μm and 1:1 stoichiometry. Strikingly, our data revealed that CaM binding to MA induces the extrusion of the myr group. However, in contrast to all known examples of CaM-binding myristoylated proteins, our data show that the myr group is exposed to solvent and not involved in CaM binding. The interactions between CaM and myr(+)MA are endothermic and entropically driven, suggesting that hydrophobic contacts are critical for binding. As revealed by NMR data, both CaM and MA appear to engage substantial regions and/or undergo significant conformational changes upon binding. We believe that our findings will provide new insights on how Gag may interact with CaM during the HIV replication cycle.

Recognition of β-calcineurin by the domains of calmodulin: thermodynamic and structural evidence for distinct roles

Calcineurin (CaN, PP2B, PPP3), a heterodimeric Ca(2+)-calmodulin-dependent Ser/Thr phosphatase, regulates swimming in Paramecia, stress responses in yeast, and T-cell activation and cardiac hypertrophy in humans. Calcium binding to CaN(B) (the regulatory subunit) triggers conformational change in CaN(A) (the catalytic subunit). Two isoforms of CaN(A) (α, β) are both abundant in brain and heart and activated by calcium-saturated calmodulin (CaM). The individual contribution of each domain of CaM to regulation of calcineurin is not known. Hydrodynamic analyses of (Ca(2+))₄-CaM(1-148) bound to βCaNp, a peptide representing its CaM-binding domain, indicated a 1:1 stoichiometry. βCaNp binding to CaM increased the affinity of calcium for the N- and C-domains equally, thus preserving intrinsic domain differences, and the preference of calcium for sites III and IV. The equilibrium constants for individual calcium-saturated CaM domains dissociating from βCaNp were ∼1 μM. A limiting K(d) ≤ 1 nM was measured directly for full-length CaM, while thermodynamic linkage analysis indicated that it was approximately 1 pM. βCaNp binding to ¹⁵N-(Ca(2+))₄-CaM(1-148) monitored by ¹⁵N/¹HN HSQC NMR showed that association perturbed the N-domain of CaM more than its C-domain. NMR resonance assignments of CaM and βCaNp, and interpretation of intermolecular NOEs observed in the ¹³C-edited and ¹²C-¹⁴N-filtered 3D NOESY spectrum indicated anti-parallel binding. The sole aromatic residue (Phe) located near the βCaNp C-terminus was in close contact with several residues of the N-domain of CaM outside the hydrophobic cleft. These structural and thermodynamic properties would permit the domains of CaM to have distinct physiological roles in regulating activation of βCaN.

Allosteric effects of the antipsychotic drug trifluoperazine on the energetics of calcium binding by calmodulin

Trifluoperazine (TFP; Stelazine) is an antagonist of calmodulin (CaM), an essential regulator of calcium-dependent signal transduction. Reports differ regarding whether, or where, TFP binds to apo CaM. Three crystallographic structures (1CTR, 1A29, and 1LIN) show TFP bound to (Ca(2+))(4)-CaM in ratios of 1, 2, or 4 TFP per CaM. In all of these, CaM domains adopt the "open" conformation seen in CaM-kinase complexes having increased calcium affinity. Most reports suggest TFP also increases calcium affinity of CaM. To compare TFP binding to apo CaM and (Ca(2+))(4)-CaM and explore differential effects on the N- and C-domains of CaM, stoichiometric TFP titrations of CaM were monitored by (15)N-HSQC NMR. Two TFP bound to apo CaM, whereas four bound to (Ca(2+))(4)-CaM. In both cases, the preferred site was in the C-domain. During the titrations, biphasic responses for some resonances suggested intersite interactions. TFP-binding sites in apo CaM appeared distinct from those in (Ca(2+))(4)-CaM. In equilibrium calcium titrations at defined ratios of TFP:CaM, TFP reduced calcium affinity at most levels tested; this is similar to the effect of many IQ-motifs on CaM. However, at the highest level tested, TFP raised the calcium affinity of the N-domain of CaM. A model of conformational switching is proposed to explain how TFP can exert opposing allosteric effects on calcium affinity by binding to different sites in the "closed," "semi-open," and "open" domains of CaM. In physiological processes, apo CaM, as well as (Ca(2+))(4)-CaM, needs to be considered a potential target of drug action.

Thermodynamics and conformational change governing domain-domain interactions of calmodulin

Calmodulin (CaM) is a small (148 amino acid), ubiquitously expressed eukaryotic protein essential for Ca(2+) regulation and signaling. This highly acidic polypeptide (pI<4) has two homologous domains (N and C), each consisting of two EF-hand Ca(2+)-binding sites. Despite significant homology, the domains have intrinsic differences in their Ca(2+)-binding properties and separable roles in regulating physiological targets such as kinases and ion channels. In mammalian full-length CaM, sites III and IV in the C-domain bind Ca(2+) cooperatively with ~10-fold higher affinity than sites I and II in the N-domain. However, the difference is only twofold when CaM is severed at residue 75, indicating that anticooperative interactions occur in full-length CaM. The Ca(2+)-binding properties of sites I and II are regulated by several factors including the interplay of interdomain linker residues far from the binding sites. Our prior thermodynamic studies showed that these residues inhibit thermal denaturation and decrease calcium affinity. Based on high-resolution structures and NMR spectra, there appear to be interactions between charged residues in the sequence 75-80 and those near the amino terminus of CaM. To explore electrostatic contributions to interdomain interactions in CaM, KCl was used to perturb the Ca(2+)-binding affinity, thermal stability, and hydrodynamic size of a nested set of recombinant mammalian CaM (rCaM) fragments terminating at residues 75, 80, 85, or 90. Potassium chloride is known to decrease Ca(2+)-binding affinity of full-length CaM. It may act directly by competition with acidic side chains that chelate Ca(2+) in the binding sites, and indirectly elsewhere in the molecule by changing tertiary constraints and conformation. In all proteins studied, KCl decreased Ca(2+)-affinity, decreased Stokes radius, and increased thermal stability, but not monotonically. Crystallographic structures of Ca(2+)-saturated rCaM(1-75) (3B32.pdb) and rCaM(1-90) (3IFK.pdb) were determined, offering cautionary notes about the effect of packing interactions on flexible linkers. This chapter describes an array of methods for characterizing system-specific thermodynamic properties that in concert govern structure and function.