Calcium-dependent stabilization of the central sequence between Met(76) and Ser(81) in vertebrate calmodulin

Coarse-Grained Modeling and Molecular Dynamics Simulations of Ca2+-Calmodulin

Calmodulin (CaM) is a calcium-binding protein that transduces signals to downstream proteins through target binding upon calcium binding in a time-dependent manner. Understanding the target binding process that tunes CaM’s affinity for the calcium ions (Ca²⁺), or vice versa, may provide insight into how Ca²⁺-CaM selects its target binding proteins. However, modeling of Ca²⁺-CaM in molecular simulations is challenging because of the gross structural changes in its central linker regions while the two lobes are relatively rigid due to tight binding of the Ca²⁺ to the calcium-binding loops where the loop forms a pentagonal bipyramidal coordination geometry with Ca²⁺. This feature that underlies the reciprocal relation between Ca²⁺ binding and target binding of CaM, however, has yet to be considered in the structural modeling. Here, we presented a coarse-grained model based on the Associative memory, Water mediated, Structure, and Energy Model (AWSEM) protein force field, to investigate the salient features of CaM. Particularly, we optimized the force field of CaM and that of Ca²⁺ ions by using its coordination chemistry in the calcium-binding loops to match with experimental observations. We presented a “community model” of CaM that is capable of sampling various conformations of CaM, incorporating various calcium-binding states, and carrying the memory of binding with various targets, which sets the foundation of the reciprocal relation of target binding and Ca²⁺ binding in future studies.

Resolved Structural States of Calmodulin in Regulation of Skeletal Muscle Calcium Release

Calmodulin (CaM) is proposed to modulate activity of the skeletal muscle sarcoplasmic reticulum (SR) calcium release channel (ryanodine receptor, RyR1 isoform) via a mechanism dependent on the conformation of RyR1-bound CaM. However, the correlation between CaM structure and functional regulation of RyR in physiologically relevant conditions is largely unknown. Here, we have used time-resolved fluorescence resonance energy transfer (TR-FRET) to study structural changes in CaM that may play a role in the regulation of RyR1. We covalently labeled each lobe of CaM (N and C) with fluorescent probes and used intramolecular TR-FRET to assess interlobe distances when CaM is bound to RyR1 in SR membranes, purified RyR1, or a peptide corresponding to the CaM-binding domain of RyR (RyRp). TR-FRET resolved an equilibrium between two distinct structural states (conformations) of CaM, each characterized by an interlobe distance and Gaussian distribution width (disorder). In isolated CaM, at low Ca²⁺, the two conformations of CaM are resolved, centered at 5 nm (closed) and 7 nm (open). At high Ca²⁺, the equilibrium shifts to favor the open conformation. In the presence of RyRp at high Ca²⁺, the closed conformation shifts to a more compact conformation and is the major component. When CaM is bound to full-length RyR1, either purified or in SR membranes, strikingly different results were obtained: 1) the two conformations are resolved and more ordered, 2) the open state is the major component, and 3) Ca²⁺ stabilized the closed conformation by a factor of two. We conclude that the Ca²⁺-dependent structural distribution of CaM bound to RyR1 is distinct from that of CaM bound to RyRp. We propose that the function of RyR1 is tuned to the Ca²⁺-dependent structural dynamics of bound CaM.

Electron paramagnetic resonance spectroscopy of nitroxide-labeled calmodulin

Calmodulin (CaM) is a highly conserved calcium-binding protein consisting of two homologous domains, each of which contains two EF-hands, that is known to bind well over 300 proteins and peptides. In most cases the (Ca(2+))(4-)form of CaM leads to the activation of a key regulatory enzyme or protein in a myriad of biological processes. Using the nitroxide spin-labeling reagent, 3-(2-iodoacetamido)-2,2,5,5-tetramethyl-1-pyrrolidinyl oxyl, bovine brain CaM was modified at 2-3 methionines with retention of activity as judged by the activation of cyclic nucleotide phosphodiesterase. X-band electron paramagnetic resonance (EPR) spectroscopy was used to measure the spectral changes upon addition of Ca(2+) to the apo-form of spin-labeled protein. A significant loss of spectral intensity, arising primarily from reductions in the heights of the low, intermediate, and high field peaks, accompanied Ca(2+) binding. The midpoint of the Ca(2+)-mediated transition determined by EPR occurred at a higher Ca(2+) concentration than that measured with circular dichroic spectroscopy and enzyme activation. Recent data have indicated that the transition from the apo-state of CaM to the fully saturated form, [(Ca(2+))(4-)CaM], contains a compact intermediate corresponding to [(Ca(2+))(2-)CaM], and the present results suggest that the spin probes are reporting on Ca(2+) binding to the last two sites in the N-terminal domain, i.e. for the [(Ca(2+))(2)-CaM] → [(Ca(2+))(4-)CaM] transition in which the compact structure becomes more extended. EPR of CaM, spin-labeled at methionines, offers a different approach for studying Ca(2+)-mediated conformational changes and may emerge as a useful technique for monitoring interactions with target proteins.

Site-directed mutagenesis of rat α-parvalbumin: replacement of canonical CD-site residues with their non-consensus counterparts from rat β-parvalbumin

Rat β-parvalbumin (β-PV) displays low divalent-ion affinity. Its CD site is distinguished by six non-consensus residues--the "CD-loop residues"--at positions 49, 50, 57-60. Additionally, leucine occupies position 85, rather than phenylalanine, the β-lineage-consensus residue. Replacement of the CD-loop residues in rat β with the canonical residues was previously found to have little effect on divalent-ion affinity, unless L85 is replaced by phenylalanine. Herein, we replace the canonical CD-loop residues in rat α-PV with their rat β-PV counterparts. Although the mutations have a generally modest impact on affinity, E59D confers Ca(2+)-specificity on the CD site, in the presence or absence of the other mutations. Despite their minimal impact on ΔG, several CD-loop mutations markedly alter ΔH, evidently by perturbing the apo-protein conformation. The L85F mutation was also examined. In wild-type rat α, L85F increases EF-site Ca(2+) affinity. In the CD-loop variants, the mutation leaves the ΔG for Ca(2+)-binding largely unaffected. However, several variants display highly exothermic binding enthalpies, indicative of ligation-linked protein-folding. Consistent with that idea, scanning-calorimetry data confirm that L85F has significantly destabilized those proteins.

Impact of methionine oxidation on calmodulin structural dynamics

We have used electron paramagnetic resonance (EPR) to examine the structural impact of oxidizing specific methionine (M) side chains in calmodulin (CaM). It has been shown that oxidation of either M109 or M124 in CaM diminishes CaM regulation of the muscle calcium release channel, the ryanodine receptor (RyR), and that mutation of M to Q (glutamine) in either case produces functional effects identical to those of oxidation. Here we have used site-directed spin labeling and double electron-electron resonance (DEER), a pulsed EPR technique that measures distances between spin labels, to characterize the structural changes resulting from these mutations. Spin labels were attached to a pair of introduced cysteine residues, one in the C-lobe (T117C) and one in the N-lobe (T34C) of CaM, and DEER was used to determine the distribution of interspin distances. Ca binding induced a large increase in the mean distance, in concert with previous X-ray crystallography and NMR data, showing a closed structure in the absence of Ca and an open structure in the presence of Ca. DEER revealed additional information about CaM's structural heterogeneity in solution: in both the presence and absence of Ca, CaM populates both structural states, one with probes separated by ∼4nm (closed) and another at ∼6nm (open). Ca shifts the structural equilibrium constant toward the open state by a factor of 13. DEER reveals the distribution of interprobe distances, showing that each of these states is itself partially disordered, with the width of each population ranging from 1 to 3nm. Both mutations (M109Q and M124Q) decrease the effect of Ca on the structure of CaM, primarily by decreasing the closed-to-open equilibrium constant in the presence of Ca. We propose that Met oxidation alters CaM's functional interaction with its target proteins by perturbing this Ca-dependent structural shift.

Distinguishing unfolding and functional conformational transitions of calmodulin using ultraviolet resonance Raman spectroscopy

Calmodulin (CaM) is a ubiquitous moderator protein for calcium signaling in all eukaryotic cells. This small calcium-binding protein exhibits a broad range of structural transitions, including domain opening and folding-unfolding, that allow it to recognize a wide variety of binding partners in vivo. While the static structures of CaM associated with its various binding activities are fairly well-known, it has been challenging to examine the dynamics of transition between these structures in real-time, due to a lack of suitable spectroscopic probes of CaM structure. In this article, we examine the potential of ultraviolet resonance Raman (UVRR) spectroscopy for clarifying the nature of structural transitions in CaM. We find that the UVRR spectral change (with 229 nm excitation) due to thermal unfolding of CaM is qualitatively different from that associated with opening of the C-terminal domain in response to Ca(2+) binding. This spectral difference is entirely due to differences in tertiary contacts at the interdomain tyrosine residue Tyr138, toward which other spectroscopic methods are not sensitive. We conclude that UVRR is ideally suited to identifying the different types of structural transitions in CaM and other proteins with conformation-sensitive tyrosine residues, opening a path to time-resolved studies of CaM dynamics using Raman spectroscopy.

Stereospecific oxidation of calmodulin by methionine sulfoxide reductase A

Methionine sulfoxide reductase A has long been known to reduce S-methionine sulfoxide, both as a free amino acid and within proteins. Recently the enzyme was shown to be bidirectional, capable of oxidizing free methionine and methionine in proteins to S-methionine sulfoxide. A feasible mechanism for controlling the directionality has been proposed, raising the possibility that reversible oxidation and reduction of methionine residues within proteins is a redox-based mechanism for cellular regulation. We undertook studies aimed at identifying proteins that are subject to site-specific, stereospecific oxidation and reduction of methionine residues. We found that calmodulin, which has nine methionine residues, is such a substrate for methionine sulfoxide reductase A. When calmodulin is in its calcium-bound form, Met77 is oxidized to S-methionine sulfoxide by methionine sulfoxide reductase A. When methionine sulfoxide reductase A operates in the reducing direction, the oxidized calmodulin is fully reduced back to its native form. We conclude that reversible covalent modification of Met77 may regulate the interaction of calmodulin with one or more of its many targets.

Retention of conformational entropy upon calmodulin binding to target peptides is driven by transient salt bridges

Calmodulin (CaM) is a highly flexible calcium-binding protein that mediates signal transduction through an ability to differentially bind to highly variable binding sequences in target proteins. To identify how binding affects CaM motions, and its relationship to conformational entropy and target peptide sequence, we have employed fully atomistic, explicit solvent molecular dynamics simulations of unbound CaM and CaM bound to five different target peptides. The calculated CaM conformational binding entropies correlate with experimentally derived conformational entropies with a correlation coefficient R(2) of 0.95. Selected side-chain interactions with target peptides restrain interhelical loop motions, acting to tune the conformational entropy of the bound complex via widely distributed CaM motions. In the complex with the most conformational entropy retention (CaM in complex with the neuronal nitric oxide synthase binding sequence), Lys-148 at the C-terminus of CaM forms transient salt bridges alternating between Glu side chains in the N-domain, the central linker, and the binding target. Additional analyses of CaM structures, fluctuations, and CaM-target interactions illuminate the interplay between electrostatic, side chain, and backbone properties in the ability of CaM to recognize and discriminate against targets by tuning its conformational entropy, and suggest a need to consider conformational dynamics in optimizing binding affinities.

Expression, purification, and characterization of proteins from high-quality combinatorial libraries of the mammalian calmodulin central linker

Combinatorial libraries offer an attractive approach towards exploring protein sequence, structure and function. Although several strategies introduce sequence diversity, the likelihood of identifying proteins with novel functions is increased when the library of genes encodes for folded and soluble structures. Here we present the first application of the binary patterning approach of combinatorial protein library design to the unique central linker region of the highly-conserved protein, calmodulin (CaM). We show that this high-quality approach translates very well to the CaM protein scaffold: all library members over-express and are functionally diverse, having a range of conformations in the presence and absence of calcium as determined by circular dichroism spectroscopy. Collectively, these data support that the binary patterning approach, when applied to the highly-conserved protein fold, can yield large collections of folded, soluble and highly-expressible proteins.

Global and local mobility of apocalmodulin monitored through fast-field cycling relaxometry

Calmodulin (CaM) is a ubiquitous eukaryotic protein with two conformationally independent domains that can bind up to two calcium ions each. In the calcium-bound state, CaM is able to regulate a vast number of cellular activities by binding to a multiplicity of target proteins in different modes. Its versatility has been ascribed to its anomalously high flexibility. The calcium-free form (apoCaM), which is the resting state of CaM in cells, is also able to functionally bind a number of protein targets, but its dynamics has received less attention. At variance with the calcium-bound form, the crystal structure of apoCaM shows a compact organization of the two domains, but NMR measurements could not detect any contact between them, thus indicating the presence of mobility in solution. The mobility of apoCaM is here investigated through protein proton relaxation rate measurements performed with a high-sensitivity fast-field cycling relaxometer. Such measurements provide direct access to the spectral density function and show that 1), the reorientation time is in agreement with a closed form of the protein; but 2), the collective order parameter is much smaller than for other well folded compact proteins, indicating that a remarkably large side-chain mobility must be present.

The extra C-terminal tail is involved in the conformation, stability changes and the N/C-domain interactions of the calmodulin-like protein from pearl oyster Pinctada fucata

Pearl oyster Pinctada fucata calmodulin-like protein (PfCaLP), containing an extra tail (D150-K161) at the C-terminal, is a novel protein involved in the regulation of oyster calcium metabolism. The purpose of this study is to gain insight into the conformational characteristics of the N/C-domain of PfCaLP, especially the detailed contribution of the extra tail to the Ca(2+)/Mg(2+)-induced conformational changes, the stability of the intact PfCaLP molecule and its C-domain, as well as to the interdomain communications in PfCaLP. Our results demonstrate that a strong interaction exists between the hydrophilic tail and the C-domain of PfCaLP. The extra tail, through affecting the C-domain conformational changes, further influences the migration rate, conformational changes, N/C-domain interactions and exposure of the hydrophobic patches of the intact PfCaLP molecule. Furthermore, the tail could actively regulate the stability of PfCaLP and its C-domain. Our studies are helpful to explain our previous finding that the tail plays important roles in PfCaLP-target interaction in the oyster calcium metabolism.

Conformational changes of calmodulin upon Ca2+ binding studied with a microfluidic mixer

A microfluidic mixer is applied to study the kinetics of calmodulin conformational changes upon Ca2+ binding. The device facilitates rapid, uniform mixing by decoupling hydrodynamic focusing from diffusive mixing and accesses time scales of tens of microseconds. The mixer is used in conjunction with multiphoton microscopy to examine the fast Ca2+-induced transitions of acrylodan-labeled calmodulin. We find that the kinetic rates of the conformational changes in two homologous globular domains differ by more than an order of magnitude. The characteristic time constants are approximately 490 micros for the transitions in the C-terminal domain and approximately 20 ms for those in the N-terminal domain of the protein. We discuss possible mechanisms for the two distinct events and the biological role of the stable intermediate, half-saturated calmodulin.

2,2,2-Trifluoroethanol changes the transition kinetics and subunit interactions in the small bacterial mechanosensitive channel MscS

2,2,2-Trifluoroethanol (TFE), a low-dielectric solvent, has recently been used as a promising tool to probe the strength of intersubunit interactions in membrane proteins. An analysis of inner membrane proteins of Escherichia coli has identified several SDS-resistant protein complexes that separate into subunits upon exposure to TFE. One of these was the homo-heptameric stretch-activated mechanosensitive channel of small conductance (MscS), a ubiquitous component of the bacterial turgor-regulation system. Here we show that a substantial fraction of MscS retains its oligomeric state in cold lithium-dodecyl-sulfate gel electrophoresis. Exposure of MscS complexes to 10-15 vol % TFE in native membranes or nonionic detergent micelles before lithium-dodecyl-sulfate electrophoresis results in a complete dissociation into monomers, suggesting that at these concentrations TFE by itself disrupts or critically compromises intersubunit interactions. Patch-clamp analysis of giant E. coli spheroplasts expressing MscS shows that exposure to TFE in lower concentrations (0.5-5.0 vol %) causes leftward shifts of the dose-response curves when applied extracellularly, and rightward shifts when added from the cytoplasmic side. In the latter case, TFE increases the rate of tension-dependent inactivation and lengthens the process of recovery to the resting state. MscS responses to pressure ramps of different speeds indicate that in the presence of TFE most channels reside in the resting state and only at tensions near the activation threshold does TFE dramatically speed up inactivation. The effect of TFE is reversible as normal channel activity returns 15-30 min after a TFE washout. We interpret the observed midpoint shifts in terms of asymmetric partitioning of TFE into the membrane and distortion of the bilayer lateral pressure profile. We also relate the increased rate of inactivation and subunit separation with the capacity of TFE to perturb buried interhelical contacts in proteins and discuss these effects in the framework of the proposed gating mechanism of MscS.

Conformational changes upon calcium binding and phosphorylation in a synthetic fragment of calmodulin

We have recently investigated by far-UV circular dichroism (CD) the effects of Ca(2+) binding and the phosphorylation of Ser 81 for the synthetic peptide CaM [54-106] encompassing the Ca(2+)-binding loops II and III and the central alpha helix of calmodulin (CaM) (Arrigoni et al., Biochemistry 2004, 43, 12788-12798). Using computational methods, we studied the changes in the secondary structure implied by these spectra with the aim to investigate the effect of Ca(2+) binding and the functional role of the phosphorylation of Ser 81 in the action of the full-length CaM. Ca(2+) binding induces the nucleation of helical structure by inducing side chain stacking of hydrophobic residues. We further investigated the effect of Ca(2+) binding by using near-UV CD spectroscopy. Molecular dynamics simulations of different fragments containing the central alpha-helix of CaM using various experimentally determined structures of CaM with bound Ca(2+) disclose the structural effects provided by the phosphorylation of Ser 81. This post-translational modification is predicted to alter the secondary structure in its surrounding and also to hinder the physiological bending of the central helix of CaM through an alteration of the hydrogen bond network established by the side chain of residue 81. Using quantum mechanical methods to predict the CD spectra for the frames obtained during the MD simulations, we are able to reproduce the relative experimental intensities in the far-UV CD spectra for our peptides. Similar conformational changes that take place in CaM [54-106] upon Ca(2+) binding and phosphorylation may occur in the full-length CaM.

Unwinding the helical linker of calcium-loaded calmodulin: A molecular dynamics study

The fold of calmodulin (CaM) consists of two globular domains connected by a helical segment (the linker), whose conformational properties play a crucial role for the protein's molecular recognition processes. Here we investigate the structural properties of the linker by performing a 11.5 ns molecular dynamics (MD) simulation of calcium-loaded human CaM in aqueous solution. The calculations are based on the AMBER force field. The calculated S2 order parameters are in good accord with NMR data: The structure of the linker in our simulations is much more flexible than that emerging from the Homo sapiens X-ray structure, consistently with the helix unwinding observed experimentally in solution. This process occurs spontaneously in a nanosecond timescale, as observed also in a very recent simulation based on the GROMOS force field. A detailed description of the mechanism that determines the linker unwinding is provided, in which electrostatic contacts between the two globular domains play a critical role. The orientation of the domains emerging from our MD calculations is consistent both with former X-ray scattering data and a recent NMR work. Based on our findings, a rationale for the experimentally measured entropy cost associated to binding to the protein's cellular partners is also given.

Rapid method for quantifying the extent of methionine oxidation in intact calmodulin

We have developed a method for rapidly quantifying the extent to which the functionally important Met144 and Met145 residues near the C-terminus of calmodulin (CaM) are converted to the corresponding sulfoxides, Met(O). The method utilizes a whole protein collision-induced dissociation (CID) approach on an electrospray ionization quadrupole time-of-flight (ESI-Q-TOF) mass spectrometer. Using standards of CaM oxidized by hydrogen peroxide (H2O2) or peroxynitrite (ONOO-), we demonstrated that CID fragmentation of the protein ions resulted in a series of C-terminal singly charged y1-y15 ions. Fragments larger than y4 exhibited mass shifts of +16 or +32 Da, corresponding to oxidation of one or two methionines, respectively. To assess the extent of oxidative modification for Met144 and Met145 to Met(O), we averaged the ratio of intensities for yn, yn+16, and yn+32 ions, where n=6-9. By alternating MS and CID scans at low and high collision energies, this technique allowed us to rapidly determine both the distribution of intact CaM oxiforms and the extent of oxidative modification in the C-terminal region of the protein in a single run. We have applied the method to studies of the repair of fully oxidized CaM by methionine sulfoxide reductases (MsrA and MsrB), which normally function in concert to reduce the S and R stereoisomers of methionine sulfoxide. We found that repair of Met(O)144 and Met(O)145 did not go to completion, but was more efficient than average Met repair. Absence of complete repair is consistent with previous studies showing that accumulation of methionine sulfoxide in CaM can occur during aging (Gao, J.; Yin, D.; Yao, Y.; Williams, T. D.; Squier, T. C. Biochemistry1998, 37, 9536-9548).

The conformationally dynamic C helix of the RIalpha subunit of protein kinase A mediates isoform-specific domain reorganization upon C subunit binding

Different isoforms of the full-length protein kinase A (PKA) regulatory subunit homodimer (R2) and the catalytic (C) subunit-bound holoenzyme (R2C2) have very different global structures despite similar molecular weights and domain organization within their primary sequences. To date, it has been the linker sequence between the R subunit dimerization/docking domain and cAMP-binding domain A that has been implicated in modulating domain interactions to give rise to these differences in global structure. The small angle solution scattering data presented here for three different isoforms of PKA heterodimer (deltaR-C) complexes reveal a role for another conformationally dynamic sequence in modulating inter-subunit and domain interactions, the C helix that connects the cAMP-binding domains A and B of the R subunit. The deltaR-C heterodimer complexes studied here were each formed with a monomeric N-terminal deletion mutant of the R subunit (deltaR) that contains the inhibitor sequence and both cAMP-binding domains. The scattering data show that type IIalpha and type IIbeta deltaR-C heterodimers are relatively compact and globular, with the C subunit contacting the inhibitor sequence and both cAMP-binding domains. In contrast, the type Ialpha heterodimer is significantly more extended, with the C subunit interacting with the inhibitor sequence and cAMP-binding domain A, whereas domain B extends out such that its surface is almost completely solvent exposed. These data implicate the C helix of RIalpha in modulating isoform-specific interdomain communication in the PKA holoenzyme, adding another layer of structural complexity to our understanding of signaling dynamics in this multisubunit, multidomain protein kinase.

A molecular dynamics study of Ca(2+)-calmodulin: evidence of interdomain coupling and structural collapse on the nanosecond timescale

A 20-ns molecular dynamics simulation of Ca(2+)-calmodulin (CaM) in explicit solvent is described. Within 5 ns, the extended crystal structure adopts a compact shape similar in dimension to complexes of CaM and target peptides but with a substantially different orientation between the N- and C-terminal domains. Significant interactions are observed between the terminal domains in this compact state, which are mediated through the same regions of CaM that bind to target peptides derived from protein kinases and most other target proteins. The process of compaction is driven by the loss of helical structure in two separate regions between residues 75-79 and 82-86, the latter being driven by unfavorable electrostatic interactions between acidic residues. In the first 5 ns of the simulation, a substantial number of contacts are observed between the first helix of the N-terminal domain and residues 74-77 of the central linker. These contacts are correlated with the closing of the second EF-hand, indicating a mechanism by which they can lower calcium affinity in the N-terminal domain.

Structural Transitions as Determinants of the Action of the Calcium-Dependent Antibiotic Daptomycin

Daptomycin is a cyclic anionic lipopeptide antibiotic recently approved for the treatment of complicated skin infections (Cubicin). Its function is dependent on calcium (as Ca2+). Circular dichroism spectroscopy indicated that daptomycin experienced two structural transitions: a transition upon interaction of daptomycin with Ca2+, and a further transition upon interaction with Ca2+ and the bacterial acidic phospholipid, phosphatidyl glycerol. The Ca2+-dependent insertion of daptomycin into model membranes promoted mild and more pronounced perturbations as assessed by the increase of lipid flip-flop and membrane leakage, respectively. The NMR structure of daptomycin indicated that Ca2+ induced a conformational change in daptomycin that increased its amphipathicity. These results are consistent with the hypothesis that the association of Ca2+ with daptomycin permits it to interact with bacterial membranes with effects that are similar to those of the cationic antimicrobial peptides.