Target enzyme recognition by calmodulin: 2.4 A structure of a calmodulin-peptide complex

The structure of Ca2+ sensor Case16 reveals the mechanism of reaction to low Ca2+ concentrations

Here we report the first crystal structure of a high-contrast genetically encoded circularly permuted green fluorescent protein (cpGFP)-based Ca²⁺ sensor, Case16, in the presence of a low Ca²⁺ concentration. The structure reveals the positioning of the chromophore within Case16 at the first stage of the Ca²⁺-dependent response when only two out of four Ca²⁺-binding pockets of calmodulin (CaM) are occupied with Ca²⁺ ions. In such a “half Ca²⁺-bound state”, Case16 is characterized by an incomplete interaction between its CaM-/M13-domains. We also report the crystal structure of the related Ca²⁺ sensor Case12 at saturating Ca²⁺ concentration. Based on this structure, we postulate that cpGFP-based Ca²⁺ sensors can form non-functional homodimers where the CaM-domain of one sensor molecule binds symmetrically to the M13-peptide of the partner sensor molecule. Case12 and Case16 behavior upon addition of high concentrations of free CaM or M13-peptide reveals that the latter effectively blocks the fluorescent response of the sensor. We speculate that the demonstrated intermolecular interaction with endogenous substrates and homodimerization can impede proper functioning of this type of Ca²⁺ sensors in living cells.

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.

Calmodulin disrupts the structure of the HIV-1 MA protein

The MA protein from HIV-1 is a small, multifunctional protein responsible for regulating various stages of the viral replication cycle. To achieve its diverse tasks, MA interacts with host cell proteins and it has been reported that one of these is the ubiquitous calcium-sensing calmodulin (CaM), which is up-regulated upon HIV-1 infection. The nature of the CaM-MA interaction has been the subject of structural studies, using peptides based on the MA sequence, that have led to conflicting conclusions. The results presented here show that CaM binds intact MA with 1:1 stoichiometry in a Ca(2+)-dependent manner and that the complex adopts a highly extended conformation in solution as revealed by small-angle X-ray scattering. Alterations in tryptophan fluorescence suggest that the two buried tryptophans (W16 and W36) located in the first two alpha-helices of MA mediate the CaM interaction. Major chemical shift changes occur in the NMR spectrum of MA upon complex formation, whereas chemical shift changes in the CaM spectrum are quite modest and are assigned to residues within the normal target protein-binding hydrophobic clefts of CaM. The NMR data indicate that CaM binds MA via its N- and C-terminal lobes and induces a dramatic conformational change involving a significant loss of secondary and tertiary structure within MA. Circular dichroism experiments suggest that MA loses approximately 20% of its alpha-helical content upon CaM binding. Thus, CaM binding is expected to impact upon the accessibility of interaction sites within MA that are involved in its various functions.

Structure and mechanism of calmodulin binding to a signaling sphingolipid reveal new aspects of lipid-protein interactions

Lipid-protein interactions are rarely characterized at a structural molecular level due to technical difficulties; however, the biological significance of understanding the mechanism of these interactions is outstanding. In this report, we provide mechanistic insight into the inhibitory complex formation of the lipid mediator sphingosylphosphorylcholine with calmodulin, the most central and ubiquitous regulator protein in calcium signaling. We applied crystallographic, thermodynamic, kinetic, and spectroscopic approaches using purified bovine calmodulin and bovine cerebral microsomal fraction to arrive at our conclusions. Here we present 1) a 1.6-Å resolution crystal structure of their complex, in which the sphingolipid occupies the conventional hydrophobic binding site on calmodulin; 2) a peculiar stoichiometry-dependent binding process: at low or high protein-to-lipid ratio calmodulin binds lipid micelles or a few lipid molecules in a compact globular conformation, respectively, and 3) evidence that the sphingolipid displaces calmodulin from its targets on cerebral microsomes. We have ascertained the specificity of the interaction using structurally related lipids as controls. Our observations reveal the structural basis of selective calmodulin inhibition by the sphingolipid. On the basis of the crystallographic and biophysical characterization of the calmodulin-sphingosylphosphorylcholine interaction, we propose a novel lipid-protein binding model, which might be applicable to other interactions as well.

An efficient algorithm for multistate protein design based on FASTER

Most of the methods that have been developed for computational protein design involve the selection of side-chain conformations in the context of a single, fixed main-chain structure. In contrast, multistate design (MSD) methods allow sequence selection to be driven by the energetic contributions of multiple structural or chemical states simultaneously. This methodology is expected to be useful when the design target is an ensemble of related states rather than a single structure, or when a protein sequence must assume several distinct conformations to function. MSD can also be used with explicit negative design to suggest sequences with altered structural, binding, or catalytic specificity. We report implementation details of an efficient multistate design optimization algorithm based on FASTER (MSD-FASTER). We subjected the algorithm to a battery of computational tests and found it to be generally applicable to various multistate design problems; designs with a large number of states and many designed positions are completely feasible. A direct comparison of MSD-FASTER and multistate design Monte Carlo indicated that MSD-FASTER discovers low-energy sequences much more consistently. MSD-FASTER likely performs better because amino acid substitutions are chosen on an energetic basis rather than randomly, and because multiple substitutions are applied together. Through its greater efficiency, MSD-FASTER should allow protein designers to test experimentally better-scoring sequences, and thus accelerate progress in the development of improved scoring functions and models for computational protein design.

Understanding protein non-folding

This review describes the family of intrinsically disordered proteins, members of which fail to form rigid 3-D structures under physiological conditions, either along their entire lengths or only in localized regions. Instead, these intriguing proteins/regions exist as dynamic ensembles within which atom positions and backbone Ramachandran angles exhibit extreme temporal fluctuations without specific equilibrium values. Many of these intrinsically disordered proteins are known to carry out important biological functions which, in fact, depend on the absence of a specific 3-D structure. The existence of such proteins does not fit the prevailing structure-function paradigm, which states that a unique 3-D structure is a prerequisite to function. Thus, the protein structure-function paradigm has to be expanded to include intrinsically disordered proteins and alternative relationships among protein sequence, structure, and function. This shift in the paradigm represents a major breakthrough for biochemistry, biophysics and molecular biology, as it opens new levels of understanding with regard to the complex life of proteins. This review will try to answer the following questions: how were intrinsically disordered proteins discovered? Why don't these proteins fold? What is so special about intrinsic disorder? What are the functional advantages of disordered proteins/regions? What is the functional repertoire of these proteins? What are the relationships between intrinsically disordered proteins and human diseases?

The role of conformational entropy in molecular recognition by calmodulin

The physical basis for high-affinity interactions involving proteins is complex and potentially involves a range of energetic contributions. Among these are changes in protein conformational entropy, which cannot yet be reliably computed from molecular structures. We have recently used changes in conformational dynamics as a proxy for changes in conformational entropy of calmodulin upon association with domains from regulated proteins. The apparent change in conformational entropy was linearly related to the overall binding entropy. This view warrants a more quantitative foundation. Here we calibrate an 'entropy meter' using an experimental dynamical proxy based on NMR relaxation and show that changes in the conformational entropy of calmodulin are a significant component of the energetics of binding. Furthermore, the distribution of motion at the interface between the target domain and calmodulin is surprisingly noncomplementary. These observations promote modification of our understanding of the energetics of protein-ligand interactions.

APIS: accurate prediction of hot spots in protein interfaces by combining protrusion index with solvent accessibility

BACKGROUND

It is well known that most of the binding free energy of protein interaction is contributed by a few key hot spot residues. These residues are crucial for understanding the function of proteins and studying their interactions. Experimental hot spots detection methods such as alanine scanning mutagenesis are not applicable on a large scale since they are time consuming and expensive. Therefore, reliable and efficient computational methods for identifying hot spots are greatly desired and urgently required.

RESULTS

In this work, we introduce an efficient approach that uses support vector machine (SVM) to predict hot spot residues in protein interfaces. We systematically investigate a wide variety of 62 features from a combination of protein sequence and structure information. Then, to remove redundant and irrelevant features and improve the prediction performance, feature selection is employed using the F-score method. Based on the selected features, nine individual-feature based predictors are developed to identify hot spots using SVMs. Furthermore, a new ensemble classifier, namely APIS (A combined model based on Protrusion Index and Solvent accessibility), is developed to further improve the prediction accuracy. The results on two benchmark datasets, ASEdb and BID, show that this proposed method yields significantly better prediction accuracy than those previously published in the literature. In addition, we also demonstrate the predictive power of our proposed method by modelling two protein complexes: the calmodulin/myosin light chain kinase complex and the heat shock locus gene products U and V complex, which indicate that our method can identify more hot spots in these two complexes compared with other state-of-the-art methods.

CONCLUSION

We have developed an accurate prediction model for hot spot residues, given the structure of a protein complex. A major contribution of this study is to propose several new features based on the protrusion index of amino acid residues, which has been shown to significantly improve the prediction performance of hot spots. Moreover, we identify a compact and useful feature subset that has an important implication for identifying hot spot residues. Our results indicate that these features are more effective than the conventional evolutionary conservation, pairwise residue potentials and other traditional features considered previously, and that the combination of our and traditional features may support the creation of a discriminative feature set for efficient prediction of hot spot residues. The data and source code are available on web site http://home.ustc.edu.cn/~jfxia/hotspot.html.

Design of multispecific protein sequences using probabilistic graphical modeling

In nature, proteins partake in numerous protein- protein interactions that mediate their functions. Moreover, proteins have been shown to be physically stable in multiple structures, induced by cellular conditions, small ligands, or covalent modifications. Understanding how protein sequences achieve this structural promiscuity at the atomic level is a fundamental step in the drug design pipeline and a critical question in protein physics. One way to investigate this subject is to computationally predict protein sequences that are compatible with multiple states, i.e., multiple target structures or binding to distinct partners. The goal of engineering such proteins has been termed multispecific protein design. We develop a novel computational framework to efficiently and accurately perform multispecific protein design. This framework utilizes recent advances in probabilistic graphical modeling to predict sequences with low energies in multiple target states. Furthermore, it is also geared to specifically yield positional amino acid probability profiles compatible with these target states. Such profiles can be used as input to randomly bias high-throughput experimental sequence screening techniques, such as phage display, thus providing an alternative avenue for elucidating the multispecificity of natural proteins and the synthesis of novel proteins with specific functionalities. We prove the utility of such multispecific design techniques in better recovering amino acid sequence diversities similar to those resulting from millions of years of evolution. We then compare the approaches of prediction of low energy ensembles and of amino acid profiles and demonstrate their complementarity in providing more robust predictions for protein design.

The effect of calcium ions and peptide ligands on the relative stabilities of the calmodulin dumbbell and compact structures

A combination of ion mobility and mass spectrometry methods was used to characterize the molecular shape of the protein calmodulin (CaM) and its complexes with calcium and a number of peptide ligands. CaM, a calcium-binding protein composed of 148 amino acid residues, was found by X-ray crystallography to occur both in a globular shape and in the shape of an extended dumbbell. Here, it was found, as solutions of CaM and CaM complexes were sprayed into the solvent-free environment of the mass spectrometer, that major structural features of the molecule and the stoichiometry of the units constituting a complex in solution were preserved in the desolvation process. Two types of CaM structures were observed in our experiments: a compact and an extended form of CaM with measured cross sections in near-perfect agreement with those calculated for the known globular and extended dumbbell X-ray geometries. Calcium-free solutions yielded predominantly an extended CaM conformation. Ca(n)(2+)-CaM complexes were observed in calcium-containing solutions, n = 0-4, with the population of the compact conformation increasing relative to the elongated conformation as n increases. For n = 4, a predominantly compact globular conformation was observed. Solutions containing the peptide CaMKII(290-309), the CaM target domain of the Ca(2+)/calmodulin-dependent protein kinase II (CaMKII) enzyme, yielded predominantly globular Ca(4)(2+)-CaM-CaMKII(290-309) complexes. Similar results were obtained with the 26-residue peptide melittin. For the 14-residue C-terminal melittin fragment, on the other hand, formation of both a 1:1 and a 1:2 CaM-peptide complex was detected. On the basis of the entirety of our results, we conclude that the collapse of extended (dumbbell-like) CaM structures into more compact globular structures occurs upon specific binding of four calcium ions. Furthermore, this calcium-induced structural collapse of CaM appears to be a prerequisite for formation of a particularly stable CaM-peptide complex involving peptides long enough to be engaged in interactions with both lobes of CaM.

The solution structure of the Mg2+ form of soybean calmodulin isoform 4 reveals unique features of plant calmodulins in resting cells

Soybean calmodulin isoform 4 (sCaM4) is a plant calcium-binding protein, regulating cellular responses to the second messenger Ca(2+). We have found that the metal ion free (apo-) form of sCaM4 possesses a half unfolded structure, with the N-terminal domain unfolded and the C-terminal domain folded. This result was unexpected as the apo-forms of both soybean calmodulin isoform 1 (sCaM1) and mammalian CaM (mCaM) are fully folded. Because of the fact that free Mg(2+) ions are always present at high concentrations in cells (0.5-2 mM), we suggest that Mg(2+) should be bound to sCaM4 in nonactivated cells. CD studies revealed that in the presence of Mg(2+) the initially unfolded N-terminal domain of sCaM4 folds into an alpha-helix-rich structure, similar to the Ca(2+) form. We have used the NMR backbone residual dipolar coupling restraints (1)D(NH), (1)D(C alpha H alpha), and (1)D(C'C alpha) to determine the solution structure of the N-terminal domain of Mg(2+)-sCaM4 (Mg(2+)-sCaM4-NT). Compared with the known structure of Ca(2+)-sCaM4, the structure of the Mg(2+)-sCaM4-NT does not fully open the hydrophobic pocket, which was further confirmed by the use of the fluorescent probe ANS. Tryptophan fluorescence experiments were used to study the interactions between Mg(2+)-sCaM4 and CaM-binding peptides derived from smooth muscle myosin light chain kinase and plant glutamate decarboxylase. These results suggest that Mg(2+)-sCaM4 does not bind to Ca(2+)-CaM target peptides and therefore is functionally similar to apo-mCaM. The Mg(2+)- and apo-structures of the sCaM4-NT provide unique insights into the structure and function of some plant calmodulins in resting cells.

Expression, purification, crystallization and preliminary X-ray analysis of calmodulin in complex with the regulatory domain of the plasma-membrane Ca2+-ATPase ACA8

Plasma-membrane Ca(2+)-ATPases (PMCAs) are calcium pumps that expel Ca(2+) from eukaryotic cells to maintain overall Ca(2+) homoeostasis and to provide local control of intracellular Ca(2+) signalling. They are of major physiological importance, with different isoforms being essential, for example, for presynaptic and postsynaptic Ca(2+) regulation in neurons, feedback signalling in the heart and sperm motility. In the resting state, PMCAs are autoinhibited by binding of their C-terminal (in mammals) or N-terminal (in plants) tail to two major intracellular loops. Activation requires the binding of calcium-bound calmodulin (Ca(2+)-CaM) to this tail and a conformational change that displaces the autoinhibitory tail from the catalytic domain. The complex between calmodulin and the regulatory domain of the plasma-membrane Ca(2+)-ATPase ACA8 from Arabidopsis thaliana has been crystallized. The crystals belonged to space group C2, with unit-cell parameters a = 176.8, b = 70.0, c = 69.8 A, beta = 113.2 degrees. A complete data set was collected to 3.0 A resolution and structure determination is in progress in order to elucidate the mechanism of PMCA activation by calmodulin.

A dynamic model of interactions of Ca2+, calmodulin, and catalytic subunits of Ca2+/calmodulin-dependent protein kinase II

During the acquisition of memories, influx of Ca2+ into the postsynaptic spine through the pores of activated N-methyl-D-aspartate-type glutamate receptors triggers processes that change the strength of excitatory synapses. The pattern of Ca2+influx during the first few seconds of activity is interpreted within the Ca2+-dependent signaling network such that synaptic strength is eventually either potentiated or depressed. Many of the critical signaling enzymes that control synaptic plasticity,including Ca2+/calmodulin-dependent protein kinase II (CaMKII), are regulated by calmodulin, a small protein that can bindup to 4 Ca2+ ions. As a first step toward clarifying how the Ca2+-signaling network decides between potentiation or depression, we have created a kinetic model of the interactions of Ca2+, calmodulin, and CaMKII that represents our best understanding of the dynamics of these interactions under conditions that resemble those in a postsynaptic spine. We constrained parameters of the model from data in the literature, or from our own measurements, and then predicted time courses of activation and autophosphorylation of CaMKII under a variety of conditions. Simulations showed that species of calmodulin with fewer than four bound Ca2+ play a significant role in activation of CaMKII in the physiological regime,supporting the notion that processing of Ca2+ signals in a spine involves competition among target enzymes for binding to unsaturated species of CaM in an environment in which the concentration of Ca2+ is fluctuating rapidly. Indeed, we showed that dependence of activation on the frequency of Ca2+ transients arises from the kinetics of interaction of fluctuating Ca2+with calmodulin/CaMKII complexes. We used parameter sensitivity analysis to identify which parameters will be most beneficial to measure more carefully to improve the accuracy of predictions. This model provides a quantitative base from which to build more complex dynamic models of postsynaptic signal transduction during learning.

Integrated protein array screening and high throughput validation of 70 novel neural calmodulin-binding proteins

Calmodulin is an essential regulator of intracellular processes in response to extracellular stimuli mediated by a rise in Ca(2+) ion concentration. To profile protein-protein interactions of calmodulin in human brain, we probed a high content human protein array with fluorophore-labeled calmodulin in the presence of Ca(2+). This protein array contains 37,200 redundant proteins, incorporating over 10,000 unique human neural proteins from a human brain cDNA library. We designed a screen to find high affinity (K(D) < or = 1 microm) binding partners of calmodulin and identified 76 human proteins from all intracellular compartments of which 72 are novel. We measured the binding kinetics of 74 targets with calmodulin using a high throughput surface plasmon resonance assay. Most of the novel calmodulin-target complexes identified have low dissociation rates (k(off) < or = 10(-3) s(-1)) and high affinity (K(D) </= 1 mum), consistent with the design of the screen. Many of the identified proteins are known to assemble in neural tissue, forming assemblies such as the spectrin scaffold and the postsynaptic density. We developed a microarray of the identified target proteins with which we can characterize the biochemistry of calmodulin for all targets in parallel. Four novel targets were verified in neural cells by co-immunoprecipitation, and four were selected for exploration of the calmodulin-binding regions. Using synthetic peptides and isothermal titration calorimetry, calmodulin binding motifs were identified in the potassium voltage-gated channel Kv6.1 (residues 474-493), calmodulin kinase-like vesicle-associated protein (residues 302-316), EF-hand domain family member A2 (residues 202-216), and phosphatidylinositol-4-phosphate 5-kinase, type I, gamma (residues 400-415).

N-myristoylated proteins, key components in intracellular signal transduction systems enabling rapid and flexible cell responses

N-myristoylation, one of the co- or post-translational modifications of proteins, has so far been regarded as necessary for anchoring of proteins to membranes. Recently, we have revealed that N(alpha)-myristoylation of several brain proteins unambiguously regulates certain protein-protein interactions that may affect signaling pathways in brain. Comparison of the amino acid sequences of myristoylated proteins including those in other organs suggests that this regulation is involved in signaling pathways not only in brain but also in other organs. Thus, it has been shown that myristoylated proteins in cells regulate the signal transduction between membranes and cytoplasmic fractions. An algorithm we have developed to identify myristoylated proteins in cells predicts the presence of hundreds of myristoylated proteins. Interestingly, a large portion of the myristoylated proteins thought to take part in signal transduction between membranes and cytoplasmic fractions are included in the predicted myristoylated proteins. If the proteins functionally regulated by myristoylation, a posttranslational protein modification, were understood as cross-talk points within the intracellular signal transduction system, known signaling pathways could thus be linked to each other, and a novel map of this intracellular network could be constructed. On the basis of our recent results, this review will highlight the multifunctional aspects of protein N-myristoylation in brain.

Tradeoff between stability and multispecificity in the design of promiscuous proteins

Natural proteins often partake in several highly specific protein-protein interactions. They are thus subject to multiple opposing forces during evolutionary selection. To be functional, such multispecific proteins need to be stable in complex with each interaction partner, and, at the same time, to maintain affinity toward all partners. How is this multispecificity acquired through natural evolution? To answer this compelling question, we study a prototypical multispecific protein, calmodulin (CaM), which has evolved to interact with hundreds of target proteins. Starting from high-resolution structures of sixteen CaM-target complexes, we employ state-of-the-art computational methods to predict a hundred CaM sequences best suited for interaction with each individual CaM target. Then, we design CaM sequences most compatible with each possible combination of two, three, and all sixteen targets simultaneously, producing almost 70,000 low energy CaM sequences. By comparing these sequences and their energies, we gain insight into how nature has managed to find the compromise between the need for favorable interaction energies and the need for multispecificity. We observe that designing for more partners simultaneously yields CaM sequences that better match natural sequence profiles, thus emphasizing the importance of such strategies in nature. Furthermore, we show that the CaM binding interface can be nicely partitioned into positions that are critical for the affinity of all CaM-target complexes and those that are molded to provide interaction specificity. We reveal several basic categories of sequence-level tradeoffs that enable the compromise necessary for the promiscuity of this protein. We also thoroughly quantify the tradeoff between interaction energetics and multispecificity and find that facilitating seemingly competing interactions requires only a small deviation from optimal energies. We conclude that multispecific proteins have been subjected to a rigorous optimization process that has fine-tuned their sequences for interactions with a precise set of targets, thus conferring their multiple cellular functions.

In nature, some proteins are more social than others, interacting with a large number of partners. These “promiscuous” proteins play key roles in cellular signaling pathways whose disruption may lead to diseases such as cancer. The amino acid sequences of such proteins must have evolved to be optimal for combined interactions with all natural partners. However, the evolutionary process leading to this promiscuity is not fully understood. We address this subject by predicting amino acid sequences that would be most compatible for interaction with each partner on its own and those most compatible for binding multiple proteins. We find that these two types of sequences are substantially different, the latter more closely resembling the natural sequences of promiscuous proteins. We also find that promiscuous proteins contain certain regions that are necessary for interfacing with all of their partners, while other regions convey specific interactions with each particular target protein. We analyze the tradeoffs required for such proteins to bind multiple partners and find that only some degree of compromise is typically needed in order to permit interactions that are seemingly antagonistic. We conclude that the simulations reported here mimic well the natural evolution of proteins that associate with multiple partners.

Calmodulin wraps around its binding domain in the plasma membrane Ca2+ pump anchored by a novel 18-1 motif

Using solution NMR spectroscopy, we obtained the structure of Ca(2+)-calmodulin (holoCaM) in complex with peptide C28 from the binding domain of the plasma membrane Ca(2+)-ATPase (PMCA) pump isoform 4b. This provides the first atomic resolution insight into the binding mode of holoCaM to the full-length binding domain of PMCA. Structural comparison of the previously determined holoCaM.C20 complex with this holoCaM.C28 complex supports the idea that the initial binding step is represented by (holoCaM.C20) and the final bound complex by (holoCaM.C28). This affirms the existing multi-step kinetic model of PMCA4b activation by CaM. The complex exhibits a new binding motif in which holoCaM is wrapped around helical C28 peptide using two anchoring residues from the peptide at relative positions 18 and 1. The anchors correspond to Phe-1110 and Trp-1093, respectively, in full-length PMCA4b, and the peptide and CaM are oriented in an anti-parallel manner. This is a greater sequence distance between anchors than in any of the known holoCaM complexes with a helical peptide. Analysis of the geometry of holoCaM-peptide binding for the cases where the target peptide adopts an alpha(D)-helix with its anchors buried in the main hydrophobic pockets of the two CaM lobes establishes that only relative sequential positions of 10, 14, 17, and 18 are allowed for the second anchor.

Dissociation of calmodulin-target peptide complexes by the lipid mediator sphingosylphosphorylcholine: implications in calcium signaling

Previously we have identified the lipid mediator sphingosylphosphorylcholine (SPC) as the first potentially endogenous inhibitor of the ubiquitous Ca2+ sensor calmodulin (CaM) (Kovacs, E., and Liliom, K. (2008) Biochem. J. 410, 427-437). Here we give mechanistic insight into CaM inhibition by SPC, based on fluorescence stopped-flow studies with the model CaM-binding domain melittin. We demonstrate that both the peptide and SPC micelles bind to CaM in a rapid and reversible manner with comparable affinities. Furthermore, we present kinetic evidence that both species compete for the same target site on CaM, and thus SPC can be considered as a competitive inhibitor of CaM-target peptide interactions. We also show that SPC disrupts the complex of CaM and the CaM-binding domain of ryanodine receptor type 1, inositol 1,4,5-trisphosphate receptor type 1, and the plasma membrane Ca2+ pump. By interfering with these interactions, thus inhibiting the negative feedback that CaM has on Ca2+ signaling, we hypothesize that SPC could lead to Ca2+ mobilization in vivo. Hence, we suggest that the action of the sphingolipid on CaM might explain the previously recognized phenomenon that SPC liberates Ca2+ from intracellular stores. Moreover, we demonstrate that unlike traditional synthetic CaM inhibitors, SPC disrupts the complex between not only the Ca2+-saturated but also the apo form of the protein and the target peptide, suggesting a completely novel regulation for target proteins that constitutively bind CaM, such as ryanodine receptors.

Structurally delineating stromal interaction molecules as the endoplasmic reticulum calcium sensors and regulators of calcium release-activated calcium entry

The endoplasmic reticulum (ER) lumen stores a crucial source of calcium (Ca2+) maintained orders of magnitude higher than the cytosol for the activation of a plethora of cellular responses transmitted in health and disease by a mutually efficient and communicative exchange of Ca2+ between compartments. A coordination of the Ca2+ signal is evident in the development of Ca2+ release-activated Ca2+ (CRAC) entry, vital to lymphocyte activation and replenishing of the ER Ca2+ stores, where modest decreases in ER luminal Ca2+ induce sustained increases in cytosolic Ca2+ sourced from steadfast extracellular Ca2+ supplies. While protein sensors that transduce Ca2+ signals in the cytosol such as calmodulin are succinctly understood, comparative data on the ER luminal Ca2+ sensors is only recently coming to light with the discovery that stromal interaction molecules (STIMs) sense variations in ER stored Ca2+ levels in the functional regulation of plasma membrane Orai proteins, the major component of CRAC channel pores. Drawing from data on the role of STIMs in the modulation of CRAC entry, this review illustrates the structural features that delimit the functional characteristics of ER Ca2+ sensors relative to well known cytoplasmic Ca2+ sensors.

Accurate prediction for atomic-level protein design and its application in diversifying the near-optimal sequence space

The task of engineering a protein to assume a target three-dimensional structure is known as protein design. Computational search algorithms are devised to predict a minimal energy amino acid sequence for a particular structure. In practice, however, an ensemble of low-energy sequences is often sought. Primarily, this is performed because an individual predicted low-energy sequence may not necessarily fold to the target structure because of both inaccuracies in modeling protein energetics and the nonoptimal nature of search algorithms employed. Additionally, some low-energy sequences may be overly stable and thus lack the dynamic flexibility required for biological functionality. Furthermore, the investigation of low-energy sequence ensembles will provide crucial insights into the pseudo-physical energy force fields that have been derived to describe structural energetics for protein design. Significantly, numerous studies have predicted low-energy sequences, which were subsequently synthesized and demonstrated to fold to desired structures. However, the characterization of the sequence space defined by such energy functions as compatible with a target structure has not been performed in full detail. This issue is critical for protein design scientists to successfully continue using these force fields at an ever-increasing pace and scale. In this paper, we present a conceptually novel algorithm that rapidly predicts the set of lowest energy sequences for a given structure. Based on the theory of probabilistic graphical models, it performs efficient inspection and partitioning of the near-optimal sequence space, without making any assumptions of positional independence. We benchmark its performance on a diverse set of relevant protein design examples and show that it consistently yields sequences of lower energy than those derived from state-of-the-art techniques. Thus, we find that previously presented search techniques do not fully depict the low-energy space as precisely. Examination of the predicted ensembles indicates that, for each structure, the amino acid identity at a majority of positions must be chosen extremely selectively so as to not incur significant energetic penalties. We investigate this high degree of similarity and demonstrate how more diverse near-optimal sequences can be predicted in order to systematically overcome this bottleneck for computational design. Finally, we exploit this in-depth analysis of a collection of the lowest energy sequences to suggest an explanation for previously observed experimental design results. The novel methodologies introduced here accurately portray the sequence space compatible with a protein structure and further supply a scheme to yield heterogeneous low-energy sequences, thus providing a powerful instrument for future work on protein design.

Response of small-scale, methyl rotors to protein-ligand association: a simulation analysis of calmodulin-peptide binding

Changes in the free energy barrier (DeltaE), entropy, and motional parameters associated with the rotation of methyl groups in a protein (calmodulin (CaM)) on binding a ligand (the calmodulin-binding domain of smooth-muscle myosin (smMLCKp)) are investigated using molecular dynamics simulation. In both the bound and uncomplexed forms of CaM, the methyl rotational free energy barriers follow skewed-Gaussian distributions that are not altered significantly upon ligand binding. However, site-specific perturbations are found. Around 11% of the methyl groups in CaM exhibit changes in DeltaE greater than 0.7 kcal/mol on binding. The rotational entropies of the methyl groups exhibit a nonlinear dependence on DeltaE. The relations are examined between motional parameters (the methyl rotational NMR order parameter and the relaxation time) and DeltaE. Low-barrier methyl group rotational order parameters deviate from ideal tetrahedrality by up to approximately 20%. There is a correlation between rotational barrier changes and proximity to the protein-peptide binding interface. Methyl groups that exhibit large changes in DeltaE are found to report on elements in the protein undergoing structural change on binding.

Structural studies of soybean calmodulin isoform 4 bound to the calmodulin-binding domain of tobacco mitogen-activated protein kinase phosphatase-1 provide insights into a sequential target binding mode

The calcium regulatory protein calmodulin (CaM) binds in a calcium-dependent manner to numerous target proteins. The calmodulin-binding domain (CaMBD) region of Nicotiana tabacum MAPK phosphatase has an amino acid sequence that does not resemble the CaMBD of any other known Ca(2+)-CaM-binding proteins. Using a unique fusion protein strategy, we have been able to obtain a high resolution solution structure of the complex of soybean Ca(2+)-CaM4 (SCaM4) and this CaMBD. Complete isotope labeling of both parts of the complex in the fusion protein greatly facilitated the structure determination by NMR. The 12-residue CaMBD region was found to bind exclusively to the C-lobe of SCaM4. A specific Trp and Leu side chain are utilized to facilitate strong binding through a novel "double anchor" motif. Moreover, the orientation of the helical peptide on the surface of Ca(2+)-SCaM4 is distinct from other known complexes. The N-lobe of Ca(2+)-SCaM4 in the complex remains free for additional interactions and could possibly act as a calcium-dependent adapter protein. Signaling through the MAPK pathway and increases in intracellular Ca(2+) are both hallmarks of the plant stress response, and our data support the notion that coordination of these responses may occur through the formation of a unique CaM-MAPK phosphatase multiprotein complex.

A targeted releasable affinity probe (TRAP) for in vivo photocrosslinking

Protein crosslinking, especially coupled to mass-spectrometric identification, is increasingly used to determine protein binding partners and protein-protein interfaces for isolated protein complexes. The modification of crosslinkers to permit their targeted use in living cells is of considerable importance for studying protein-interaction networks, which are commonly modulated through weak interactions that are formed transiently to permit rapid cellular response to environmental changes. We have therefore synthesized a targeted and releasable affinity probe (TRAP) consisting of a biarsenical fluorescein linked to benzophenone that binds to a tetracysteine sequence in a protein engineered for specific labeling. Here, the utility of TRAP for capturing protein binding partners upon photoactivation of the benzophenone moiety has been demonstrated in living bacteria and mammalian cells. In addition, ligand exchange of the arsenic-sulfur bonds between TRAP and the tetracysteine sequence to added dithiols results in fluorophore transfer to the crosslinked binding partner. In isolated protein complexes, this release from the original binding site permits the identification of the proximal binding interface through mass spectrometric fragmentation and computational sequence identification.

Regulation of interdomain interactions by calmodulin in inducible nitric-oxide synthase

Nitric-oxide synthases (NOSs) catalyze the conversion of l-arginine to nitric oxide and citrulline. There are three NOS isozymes, each with a different physiological role: neuronal NOS, endothelial NOS, and inducible NOS (iNOS). NOSs consist of an N-terminal oxygenase domain and a C-terminal reductase domain, linked by a calmodulin (CaM)-binding region. CaM is required for NO production, but unlike other NOS isozymes, iNOS binds CaM independently of the exogenous Ca(2+) concentration. We have co-expressed CaM and the FMN domain of human iNOS, which includes the CaM-binding region. The Ca(2+)-bound protein complex (CaCaMxFMN) forms an air-stable semiquinone when reduced with NADPH and reduces cytochrome c when reconstituted with the iNOS FAD/NADPH domain. We have solved the crystal structure of the CaCaMxFMN complex in four different conformations, each with a different relative orientation, between the FMN domain and the bound CaM. The CaM-binding region together with bound CaM forms a hinge, pivots on the conserved Arg(536), and regulates electron transfer from FAD to FMN and from FMN to heme by adjusting the relative orientation and distance among the three cofactors. In addition, the relative orientations of the N- and C-terminal lobes of CaM are also different among the four conformations, suggesting that the flexibility between the two halves of CaM also contributes to the fine tuning of the orientation/distance between the redox centers. The data demonstrate a possible mode for precise control of electron transfer by altering the distance and orientation of redox centers in a protein displaying domain movement.

Novel shuttling domain in a regulator (RSC1A1) of transporter SGLT1 steers cell cycle-dependent nuclear location

The gene product of RSC1A1, RS1, participates in the regulation of the Na(+)-D-glucose cotransporter SGLT1. RS1 inhibits release of SGLT1 from the trans Golgi network. In subconfluent LLC-PK(1) cells, RS1 migrates into the nucleus and modulates transcription of SGLT1, whereas most confluent cells do not contain RS1 in the nuclei. We showed that confluence-dependent nuclear location of RS1 is because of different phases of the cell cycle and identified a RS1 nuclear shuttling domain (RNS) with an associated protein kinase C (PKC) phosphorylation site (RNS-PKC) that mediates cell cycle-dependent nuclear location. RNS-PKC contains a novel non-conventional nuclear localization signal interacting with importin beta1, a nuclear export signal mediating export via protein CRM1 and a Ca(2+)-dependent calmodulin binding site. PKC and calmodulin compete for binding to RNS-PKC. Mutagenesis experiments and analyses of the phosphorylation status suggest the following sequences of events. Subconfluent cells without and with synchronization to the G2/M phase contain non-phosphorylated RNS-PKC that mediates nuclear import of RS1 but not its export. During confluence or synchronization of subconfluent cells to the G2/M phase, phosphorylation of RNS-PKC mediates rapid nuclear export of RS1.

StoneHinge: hinge prediction by network analysis of individual protein structures

Hinge motions are important for molecular recognition, and knowledge of their location can guide the sampling of protein conformations for docking. Predicting domains and intervening hinges is also important for identifying structurally self-determinate units and anticipating the influence of mutations on protein flexibility and stability. Here we present StoneHinge, a novel approach for predicting hinges between domains using input from two complementary analyses of noncovalent bond networks: StoneHingeP, which identifies domain-hinge-domain signatures in ProFlex constraint counting results, and StoneHingeD, which does the same for DomDecomp Gaussian network analyses. Predictions for the two methods are compared to hinges defined in the literature and by visual inspection of interpolated motions between conformations in a series of proteins. For StoneHingeP, all the predicted hinges agree with hinge sites reported in the literature or observed visually, although some predictions include extra residues. Furthermore, no hinges are predicted in six hinge-free proteins. On the other hand, StoneHingeD tends to overpredict the number of hinges, while accurately pinpointing hinge locations. By determining the consensus of their results, StoneHinge improves the specificity, predicting 11 of 13 hinges found both visually and in the literature for nine different open protein structures, and making no false-positive predictions. By comparison, a popular hinge detection method that requires knowledge of both the open and closed conformations finds 10 of the 13 known hinges, while predicting four additional, false hinges.

Time-dependent autoinactivation of phospho-Thr286-alphaCa2+/calmodulin-dependent protein kinase II

Ca²⁺/calmodulin-dependent protein kinase II (αCaMKII) is thought to exert its role in memory formation by autonomous Ca²⁺-independent persistent activity conferred by Thr²⁸⁶ autophosphorylation, allowing the enzyme to remain active even when intracellular [Ca²⁺] has returned to resting levels. Ca²⁺ sequestration-induced inhibition, caused by a burst of Thr^305/306 autophosphorylation via calmodulin (CaM) dissociation from the Thr^305/306 sites, is in conflict with this view. The processes of CaM binding, autophosphorylation, and inactivation are dissected to resolve this conflict. Upon Ca²⁺ withdrawal, CaM sequential domain dissociation is observed, starting with the rapid release of the first (presumed N-terminal) CaM lobe, thought to be bound at the Thr^305/306 sites. The time courses of Thr^305/306 autophosphorylation and inactivation, however, correlate with the slow dissociation of the second (presumed C-terminal) CaM lobe. Exposure of the Thr^305/306 sites is thus not sufficient for their autophosphorylation. Moreover, Thr^305/306 autophosphorylation and autoinactivation are shown to occur in the continuous presence of Ca²⁺ and bound Ca²⁺/CaM by time courses similar to those seen following Ca²⁺ sequestration. Our investigation of the activity and mechanisms of phospho-Thr₂₈₆-αCaMKII thus shows time-dependent autoinactivation, irrespective of the continued presence of Ca²⁺ and CaM, allowing a very short, if any, time window for Ca²⁺/CaM-free phospho-Thr²⁸⁶-αCaMKII activity. Physiologically, the time-dependent autoinactivation mechanisms of phospho-Thr²⁸⁶-αCaMKII (t½ of ∼50 s at 37 °C) suggest a transient kinase activity of ∼1 min duration in the induction of long term potentiation and thus memory formation.

Structural insights into the calmodulin-Munc13 interaction obtained by cross-linking and mass spectrometry

Munc13 proteins are essential regulators of synaptic vesicle priming and play a key role in adaptive synaptic plasticity phenomena. We recently identified and characterized the Ca(2+)-dependent interaction of Munc13 and calmodulin (CaM) as the molecular mechanism linking changes in residual Ca(2+) concentrations to presynaptic vesicle priming and short-term plasticity. Here, we used peptidic photoprobes covering the established CaM-binding motif of Munc13 for photoaffinity labeling (PAL) of CaM, followed by structural characterization of the covalent photoadducts. Our innovative analytical workflow based on isotopically labeled CaM and mass spectrometry revealed that, in the bound state, the hydrophobic anchor residue of the CaM-binding motif in Munc13s contacts two distinct methionine residues in the C-terminal domain of CaM. To address the orientation of the peptide during binding, we obtained additional distance constraints from the mass spectrometric analysis of chemically cross-linked CaM-Munc13 peptide adducts. The constraints from both complementary cross-linking approaches were integrated into low-resolution three-dimensional structure models of the CaM-Munc13 peptide complexes. Our experimental data are best compatible with the structure of the complex formed by CaM and a CaM-binding peptide derived from neuronal NO synthase and show that Munc13-1 and ubMunc13-2 bind to CaM in an antiparallel orientation through a 1-5-8 motif. The structural information about the CaM-Munc13 peptide complexes will facilitate the design of Munc13 variants with altered CaM affinity and thereby advance the detailed functional analysis of the role of Munc13 proteins in synaptic transmission and plasticity.

Functional aspects of protein flexibility

Proteins are dynamic entities, and they possess an inherent flexibility that allows them to function through molecular interactions within the cell, among cells and even between organisms. Appreciation of the non-static nature of proteins is emerging, but to describe and incorporate this into an intuitive perception of protein function is challenging. Flexibility is of overwhelming importance for protein function, and the changes in protein structure during interactions with binding partners can be dramatic. The present review addresses protein flexibility, focusing on protein-ligand interactions. The thermodynamics involved are reviewed, and examples of structure-function studies involving experimentally determined flexibility descriptions are presented. While much remains to be understood about protein flexibility, it is clear that it is encoded within their amino acid sequence and should be viewed as an integral part of their structure.

Energetics of calmodulin domain interactions with the calmodulin binding domain of CaMKII

Calmodulin (CaM) is an essential eukaryotic calcium receptor that regulates many kinases, including CaMKII. Calcium-depleted CaM does not bind to CaMKII under physiological conditions. However, binding of (Ca(2+))(4)-CaM to a basic amphipathic helix in CaMKII releases auto-inhibition of the kinase. The crystal structure of CaM bound to CaMKIIp, a peptide representing the CaM-binding domain (CaMBD) of CaMKII, shows an antiparallel interface: the C-domain of CaM primarily contacts the N-terminal half of the CaMBD. The two domains of calcium-saturated CaM are believed to play distinct roles in releasing auto-inhibition. To investigate the underlying mechanism of activation, calcium-dependent titrations of isolated domains of CaM binding to CaMKIIp were monitored using fluorescence anisotropy. The binding affinity of CaMKIIp for the domains of CaM increased upon saturation with calcium, with the C-domain having a 35-fold greater affinity than the N-domain. Because the interdomain linker of CaM regulates calcium-binding affinity and contribute to conformational change, the role of each CaM domain was explored further by investigating effects of CaMKIIp on site-knockout mutants affecting the calcium-binding sites of a single domain. Investigation of the thermodynamic linkage between saturation of individual calcium-binding sites and CaM-domain binding to CaMKIIp showed that calcium binding to Sites III and IV was sufficient to recapitulate the behavior of (Ca(2+))(4)-CaM. The magnitude of favorable interdomain cooperativity varied depending on which of the four calcium-binding sites were mutated, emphasizing differential regulatory roles for the domains of CaM, despite the high degree of homology among the four EF-hands of CaM.

Two synthetic peptides corresponding to the proximal heme-binding domain and CD1 domain of human endothelial nitric-oxide synthase inhibit the oxygenase activity by interacting with CaM

Human endothelial nitric-oxide synthase (eNOS) is a complex enzyme, requiring binding of calmodulin (CaM) for electron transfer. The prevailing view is that calcium-activated CaM binds eNOS at the canonical binding site located at residues 493-510, which induces a conformational change to facilitate electron transfer. Here we demonstrated that the CaM enhances the rate of electron transfer from NADPH to FAD on a truncated eNOS FAD subdomain (residues 682-1204) purified from baculovirus-infected Sf9 cells, suggesting more complicated regulatory mechanism of CaM on eNOS. Metabolically (35)S-labeled CaM overlay on fusion proteins spanning the entire linear sequence of eNOS revealed three positive (35)S-CaM binding fragments: sequence 66-205, sequence 460-592, and sequence 505-759. Synthetic peptides derived from these fragments are tested for their effects on CaM binding and eNOS catalytic activities. Peptides corresponding to the proximal heme-binding site (E1, residues 174-193) and the CD1 linker connecting FAD/FMN subdomains (E4, residues 729-757) bind CaM at both high Ca(2+) (Ca(2+)CaM) and low Ca(2+) (apoCaM) concentrations, whereas peptide of the canonical CaM-binding helix (E2, residues 493-510) binds only Ca(2+)CaM. All three peptides E1, E2 and E4 significantly inhibit oxygenase activity in a concentration-dependent manner, but only E2 effectively inhibits reductase activity. Concurrent experiments with human iNOS showed major differences in the CaM binding properties between eNOS and iNOS. The results suggest that multiple regions of eNOS might interact with CaM with differential Ca(2+) sensitivity in vivo. A possible mechanism in regulating eNOS activation and deactivation is proposed.

Identification and characterization of wolframin, the product of the wolfram syndrome gene (WFS1), as a novel calmodulin-binding protein

To search for calmodulin (CaM) targets, we performed affinity chromatography purification of a rat brain extract using CaM fused with GST as the affinity ligand. Proteomic analysis was then carried out to identify CaM-binding proteins. In addition to identifying 36 known CaM-binding proteins, including CaM kinases, calcineurin, nNOS, the IP(3) receptor, and Ca(2+)-ATPase, we identified an ER transmembrane protein, wolframin [the product of the Wolfram syndrome gene (WFS1)] as interacting. A CaM overlay and an immunoprecipitation assay revealed that wolframin is capable of binding the Ca(2+)/CaM complex in vitro and in transfected cells. Surface plasmon resonance analysis and zero-length cross-linking showed that the N-terminal cytoplasmic domain (residues 2-285) of wolframin binds to an equimolar unit of CaM in a Ca(2+)-dependent manner with a K(D) for CaM of 0.15 muM. Various truncation and deletion mutants showed that the Ca(2+)/CaM binding region in wolframin is located from Glu90 to Trp186. Furthermore, we demonstrated that three mutations (Ala127Thr, Ala134Thr, and Arg178Pro) associated with Wolfram syndrome completely abolished CaM binding of wolframin. This observation may indicate that CaM binding is important for wolframin function and that impairment of this interaction by mutation contributes to the pathology seen in Wolfram syndrome.

A rapid coarse residue-based computational method for x-ray solution scattering characterization of protein folds and multiple conformational states of large protein complexes

We present a coarse residue-based computational method to rapidly compute the solution scattering profile from a protein with dynamical fluctuations. The method is built upon a coarse-grained (CG) representation of the protein. This CG representation takes advantage of the intrinsic low-resolution and CG nature of solution scattering data. It allows rapid scattering determination from a large number of conformations that can be extracted from CG simulations to obtain scattering characterization of protein conformations. The method includes several important elements, effective residue structure factors derived from the Protein Data Bank, explicit treatment of water molecules in the hydration layer at the surface of the protein, and an ensemble average of scattering from a variety of appropriate conformations to account for macromolecular flexibility. This simplified method is calibrated and illustrated to accurately reproduce the experimental scattering curve of Hen egg white lysozyme. We then illustrated the applications of this CG method by computing the solution scattering patterns of several representative protein folds and multiple conformational states. The results suggest that solution scattering data, when combined with the reliable computational method that we developed, show great potential for a better structural description of multidomain complexes in different functional states, and for recognizing structural folds when sequence similarity to a protein of known structure is low.

Calculation of proteins' total side-chain torsional entropy and its influence on protein-ligand interactions

Despite the high density within a typical protein fold, the ensemble of sterically permissible side-chain repackings is vast. Here, we examine the extent of this variability that survives energetic biases due to van der Waals interactions, hydrogen bonding, salt bridges, and solvation. Monte Carlo simulations of an atomistic model exhibit thermal fluctuations among a diverse set of side-chain arrangements, even with the peptide backbone fixed in its crystallographic conformation. We have quantified the torsional entropy of this native-state ensemble, relative to that of a noninteracting reference system, for 12 small proteins. The reduction in entropy per rotatable bond due to each kind of interaction is remarkably consistent across this set of molecules. To assess the biophysical importance of these fluctuations, we have estimated side-chain entropy contributions to the binding affinity of several peptide ligands with calmodulin. Calculations for our fixed-backbone model correlate very well with experimentally determined binding entropies over a range spanning more than 80 kJ/(mol x 308 K).

Determinants in CaV1 channels that regulate the Ca2+ sensitivity of bound calmodulin

Calmodulin binds to IQ motifs in the alpha(1) subunit of Ca(V)1.1 and Ca(V)1.2, but the affinities of calmodulin for the motif and for Ca(2+) are higher when bound to Ca(V)1.2 IQ. The Ca(V)1.1 IQ and Ca(V)1.2 IQ sequences differ by four amino acids. We determined the structure of calmodulin bound to Ca(V)1.1 IQ and compared it with that of calmodulin bound to Ca(V)1.2 IQ. Four methionines in Ca(2+)-calmodulin form a hydrophobic binding pocket for the peptide, but only one of the four nonconserved amino acids (His-1532 of Ca(V)1.1 and Tyr-1675 of Ca(V)1.2) contacts this calmodulin pocket. However, Tyr-1675 in Ca(V)1.2 contributes only modestly to the higher affinity of this peptide for calmodulin; the other three amino acids in Ca(V)1.2 contribute significantly to the difference in the Ca(2+) affinity of the bound calmodulin despite having no direct contact with calmodulin. Those residues appear to allow an interaction with calmodulin with one lobe Ca(2+)-bound and one lobe Ca(2+)-free. Our data also provide evidence for lobe-lobe interactions in calmodulin bound to Ca(V)1.2.

Domain swapping and different oligomeric States for the complex between calmodulin and the calmodulin-binding domain of calcineurin a

BACKGROUND

Calmodulin (CaM) is a ubiquitously expressed calcium sensor that engages in regulatory interactions with a large number of cellular proteins. Previously, a unique mode of CaM target recognition has been observed in the crystal structure of a complex between CaM and the CaM-binding domain of calcineurin A.

METHODOLOGY/PRINCIPAL FINDINGS

We have solved a high-resolution crystal structure of a complex between CaM and the CaM-binding domain of calcineurin A in a novel crystal form, which shows a dimeric assembly of calmodulin, as observed before in the crystal state. We note that the conformation of CaM in this complex is very similar to that of unliganded CaM, and a detailed analysis revels that the CaM-binding motif in calcineurin A is of a novel '1-11' type. However, using small-angle X-ray scattering (SAXS), we show that the complex is fully monomeric in solution, and a structure of a canonically collapsed CaM-peptide complex can easily be fitted into the SAXS data. This result is also supported by size exclusion chromatography, where the addition of the ligand peptide decreases the apparent size of CaM. In addition, we studied the energetics of binding by isothermal titration calorimetry and found them to closely resemble those observed previously for ligand peptides from CaM-dependent kinases.

CONCLUSIONS/SIGNIFICANCE

Our results implicate that CaM can also form a complex with the CaM-binding domain of calcineurin in a 1 ratio 1 stoichiometry, in addition to the previously observed 2 ratio 2 arrangement in the crystal state. At the structural level, going from 2 ratio 2 association to two 1 ratio 1 complexes will require domain swapping in CaM, accompanied by the characteristic bending of the central linker helix between the two lobes of CaM.

Calmodulin mediates the Ca2+-dependent regulation of Cx44 gap junctions

We have shown previously that the Ca2+-dependent inhibition of lens epithelial cell-to-cell communication is mediated in part by the direct association of calmodulin (CaM) with connexin43 (Cx43), the major connexin in these cells. We now show that elevation of [Ca2+](i) in HeLa cells transfected with the lens fiber cell gap junction protein sheep Cx44 also results in the inhibition of cell-to-cell dye transfer. A peptide comprising the putative CaM binding domain (aa 129-150) of the intracellular loop region of this connexin exhibited a high affinity, stoichiometric interaction with Ca2+-CaM. NMR studies indicate that the binding of Cx44 peptide to CaM reflects a classical embracing mode of interaction. The interaction is an exothermic event that is both enthalpically and entropically driven in which electrostatic interactions play an important role. The binding of the Cx44 peptide to CaM increases the CaM intradomain cooperativity and enhances the Ca2+-binding affinities of the C-domain of CaM more than twofold by slowing the rate of Ca2+ release from the complex. Our data suggest a common mechanism by which the Ca2+-dependent inhibition of the alpha-class of gap junction proteins is mediated by the direct association of an intracellular loop region of these proteins with Ca2+-CaM.

Separate elements within a single IQ-like motif in adenylyl cyclase type 8 impart ca2+/calmodulin binding and autoinhibition

The ubiquitous Ca²⁺-sensing protein calmodulin (CaM) fulfills its numerous signaling functions through a wide range of modular binding and activation mechanisms. By activating adenylyl cyclases (ACs) 1 and 8, Ca²⁺ acting via calmodulin impacts on the signaling of the other major cellular second messenger cAMP. In possessing two CaM-binding domains, a 1-5-8-14 motif at the N terminus and an IQ-like motif (IQlm) at the C terminus, AC8 offers particularly sophisticated regulatory possibilities. The IQlm has remained unexplored beyond the suggestion that it bound CaM, and the larger C2b region of which it is part was involved in the relief of autoinhibition of AC8. Here we attempt to distinguish the function of individual residues of the IQlm. From a complementary approach of in vitro and cell population AC activity assays, as well as CaM binding, we propose that the IQlm alone, and not the majority of the C2b, imparts CaM binding and autoinhibitory functions. Moreover, this duality of function is spatially separated and depends on amino acid side-chain character. Accordingly, residues critical for CaM binding are positively charged and clustered toward the C terminus, and those essential for the maintenance of autoinhibition are hydrophobic and more N-terminal. Secondary structure prediction of the IQlm supports this separation, with an ideally placed break in the α-helical nature of the sequence. We additionally find that the N and C termini of AC8 interact, which is an association specifically abrogated by fully Ca²⁺-bound, but not Ca²⁺-free, CaM. These data support a sophisticated activation mechanism of AC8 by CaM, in which the duality of the IQlm function is critical.

Crystal structures of the GCaMP calcium sensor reveal the mechanism of fluorescence signal change and aid rational design

The genetically encoded calcium indicator GCaMP2 shows promise for neural network activity imaging, but is currently limited by low signal-to-noise ratio. We describe x-ray crystal structures as well as solution biophysical and spectroscopic characterization of GCaMP2 in the calcium-free dark state, and in two calcium-bound bright states: a monomeric form that dominates at intracellular concentrations observed during imaging experiments and an unexpected domain-swapped dimer with decreased fluorescence. This series of structures provides insight into the mechanism of Ca2+-induced fluorescence change. Upon calcium binding, the calmodulin (CaM) domain wraps around the M13 peptide, creating a new domain interface between CaM and the circularly permuted enhanced green fluorescent protein domain. Residues from CaM alter the chemical environment of the circularly permuted enhanced green fluorescent protein chromophore and, together with flexible inter-domain linkers, block solvent access to the chromophore. Guided by the crystal structures, we engineered a series of GCaMP2 point mutants to probe the mechanism of GCaMP2 function and characterized one mutant with significantly improved signal-to-noise. The mutation is located at a domain interface and its effect on sensor function could not have been predicted in the absence of structural data.

Thermodynamic effects of noncoded and coded methionine substitutions in calmodulin

The methionine residues in the calcium (Ca2+) regulatory protein calmodulin (CaM) are structurally and functionally important. They are buried within the N- and C-domains of apo-CaM but become solvent-exposed in Ca2+-CaM, where they interact with numerous target proteins. Previous structural studies have shown that methionine substitutions to the noncoded amino acids selenomethionine, ethionine, or norleucine, or mutation to leucine do not impact the main chain structure of CaM. Here we used differential scanning calorimetry to show that these substitutions enhance the stability of both domains, with the largest increase in melting temperature (19-26 degrees C) achieved with leucine or norleucine in the apo-C-domain. Nuclear magnetic resonance spectroscopy experiments also revealed the loss of a slow conformational exchange process in the Leu-substituted apo-C-domain. In addition, isothermal titration calorimetry experiments revealed considerable changes in the enthalpy and entropy of target binding to apo-CaM and Ca2+-CaM, but the free energy of binding was largely unaffected due to enthalpy-entropy compensation. Collectively, these results demonstrate that noncoded and coded methionine substitutions can be accommodated in CaM because of the structural plasticity of the protein. However, adjustments in side-chain packing and dynamics lead to significant differences in protein stability and the thermodynamics of target binding.