Solution structures of the N-terminal domain of yeast calmodulin: Ca2+-dependent conformational change and its functional implication

Novel kinetoplastid-specific cAMP binding proteins identified by RNAi screening for cAMP resistance in Trypanosoma brucei

Cyclic AMP signalling in trypanosomes differs from most eukaryotes due to absence of known cAMP effectors and cAMP independence of PKA. We have previously identified four genes from a genome-wide RNAi screen for resistance to the cAMP phosphodiesterase (PDE) inhibitor NPD-001. The genes were named cAMP Response Protein (CARP) 1 through 4. Here, we report an additional six CARP candidate genes from the original sample, after deep sequencing of the RNA interference target pool retrieved after NPD-001 selection (RIT-seq). The resistance phenotypes were confirmed by individual RNAi knockdown. Highest level of resistance to NPD-001, approximately 17-fold, was seen for knockdown of CARP7 (Tb927.7.4510). CARP1 and CARP11 contain predicted cyclic AMP binding domains and bind cAMP as evidenced by capture and competition on immobilised cAMP. CARP orthologues are strongly enriched in kinetoplastid species, and CARP3 and CARP11 are unique to Trypanosoma. Localization data and/or domain architecture of all CARPs predict association with the T. brucei flagellum. This suggests a crucial role of cAMP in flagellar function, in line with the cell division phenotype caused by high cAMP and the known role of the flagellum for cytokinesis. The CARP collection is a resource for discovery of unusual cAMP pathways and flagellar biology.

Enhanced conformational sampling technique provides an energy landscape view of large-scale protein conformational transitions

A novel in silico approach (NMA–ITS) is introduced to rapidly and effectively sample the configuration space and give quantitative data for exploring the conformational changes of proteins.

Linear Response Path Following: A Molecular Dynamics Method To Simulate Global Conformational Changes of Protein upon Ligand Binding

Molecular functions of proteins are often fulfilled by global conformational changes that couple with local events such as the binding of ligand molecules. High molecular complexity of proteins has, however, been an obstacle to obtain an atomistic view of the global conformational transitions, imposing a limitation on the mechanistic understanding of the functional processes. In this study, we developed a new method of molecular dynamics (MD) simulation called the linear response path following (LRPF) to simulate a protein's global conformational changes upon ligand binding. The method introduces a biasing force based on a linear response theory, which determines a local reaction coordinate in the configuration space that represents linear coupling between local events of ligand binding and global conformational changes and thus provides one with fully atomistic models undergoing large conformational changes without knowledge of a target structure. The overall transition process involving nonlinear conformational changes is simulated through iterative cycles consisting of a biased MD simulation with an updated linear response force and a following unbiased MD simulation for relaxation. We applied the method to the simulation of global conformational changes of the yeast calmodulin N-terminal domain and successfully searched out the end conformation. The atomistically detailed trajectories revealed a sequence of molecular events that properly lead to the global conformational changes and identified key steps of local-global coupling that induce the conformational transitions. The LRPF method provides one with a powerful means to model conformational changes of proteins such as motors and transporters where local-global coupling plays a pivotal role in their functional processes.

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

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

Surface dynamics in allosteric regulation of protein-protein interactions: modulation of calmodulin functions by Ca2+

Knowledge of the structural basis of protein-protein interactions (PPI) is of fundamental importance for understanding the organization and functioning of biological networks and advancing the design of therapeutics which target PPI. Allosteric modulators play an important role in regulating such interactions by binding at site(s) orthogonal to the complex interface and altering the protein's propensity for complex formation. In this work, we apply an approach recently developed by us for analyzing protein surfaces based on steered molecular dynamics simulation (SMD) to the study of the dynamic properties of functionally distinct conformations of a model protein, calmodulin (CaM), whose ability to interact with target proteins is regulated by the presence of the allosteric modulator Ca(2+). Calmodulin is a regulatory protein that acts as an intracellular Ca(2+) sensor to control a wide variety of cellular processes. We demonstrate that SMD analysis is capable of pinpointing CaM surfaces implicated in the recognition of both the allosteric modulator Ca(2+) and target proteins. Our analysis of changes in the dynamic properties of the CaM backbone elicited by Ca(2+) binding yielded new insights into the molecular mechanism of allosteric regulation of CaM-target interactions.

Specific conformation and Ca(2+)-binding mode of yeast calmodulin: insight into evolutionary development

The vertebrate calmodulin is configured with two structurally independent globular lobes in N- and C-terminus, and a flexible central linker. Distinctly, two lobes of calmodulin from Saccharomyces cerevisiae (yCaM) interact and influence the Ca(2+)-binding profile of each other. We explored this further using the mutant proteins with eliminated Ca(2+)-binding ability in one of the lobes and found that the Ca(2+)-bound N-lobe associates with the Ca(2+)-free C-lobe to gain the Ca(2+) affinity of a wild-type level. Next, analysing series of C-terminal residue truncation mutant, we found that the truncation of C-terminal three residues induce the hyper Ca(2+) affinity. These residues are also important for the general structural behaviour of calmodulin, such as Ca(2+)-induced slow mobility shift in polyacrylamide gel electrophoresis and for the ability to activate Cmk1p (yeast calmodulin kinase). These suggest: (i) when Ca(2+) occupies only N-lobe, two lobes interact and form the stable intermediate leading to a proper level of Ca(2+) affinity; (ii) the C-terminal three residues are required to prohibit abnormal stabilization of the intermediate promoting abnormally high Ca(2+) affinity and for recognition of target enzymes. A model for Ca(2+) and target bindings of yCaM is proposed. Evolutional aspect concerning the biological significance of this model was discussed.

Activation of the edema factor of Bacillus anthracis by calmodulin: evidence of an interplay between the EF-calmodulin interaction and calcium binding

Calmodulin (CaM) is a remarkably flexible protein which can bind multiple targets in response to changes in intracellular calcium concentration. It contains four calcium-binding sites, arranged in two globular domains. The calcium affinity of CaM N-terminal domain (N-CaM) is dramatically reduced when the complex with the edema factor (EF) of Bacillus anthracis is formed. Here, an atomic explanation for this reduced affinity is proposed through molecular dynamics simulations and free energy perturbation calculations of the EF-CaM complex starting from different crystallographic models. The simulations show that electrostatic interactions between CaM and EF disfavor the opening of N-CaM domains usually induced by calcium binding. Relative calcium affinities of the N-CaM binding sites are probed by free energy perturbation, and dissociation probabilities are evaluated with locally enhanced sampling simulations. We show that EF impairs calcium binding on N-CaM through a direct conformational restraint on Site 1, by an indirect destabilization of Site 2, and by reducing the cooperativity between the two sites.

Approaches to the assignment of (19)F resonances from 3-fluorophenylalanine labeled calmodulin using solution state NMR

Traditional single site replacement mutations (in this case, phenylalanine to tyrosine) were compared with methods which exclusively employ (15)N and (19)F-edited two- and three-dimensional NMR experiments for purposes of assigning (19)F NMR resonances from calmodulin (CaM), biosynthetically labeled with 3-fluorophenylalanine (3-FPhe). The global substitution of 3-FPhe for native phenylalanine was tolerated in CaM as evidenced by a comparison of (1)H-(15)N HSQC spectra and calcium binding assays in the presence and absence of 3-FPhe. The (19)F NMR spectrum reveals six resolved resonances, one of which integrates to three 3-FPhe species, making for a total of eight fluorophenylalanines. Single phenylalanine to tyrosine mutants of five phenylalanine positions resulted in (19)F NMR spectra with significant chemical shift perturbations of the remaining resonances, and provided only a single definitive assignment. Although (1)H-(19)F heteronucleclear NOEs proved weak, (19)F-edited (1)H-(1)H NOESY connectivities were relatively easy to establish by making use of the (3)J(FH) coupling between the fluorine nucleus and the adjacent fluorophenylalanine delta proton. (19)F-edited NOESY connectivities between the delta protons and alpha and beta nuclei in addition to (15)N-edited (1)H, (1)H NOESY crosspeaks proved sufficient to assign 4 of 8 (19)F resonances. Controlled cleavage of the protein into two fragments using trypsin, and a repetition of the above 2D and 3D techniques resulted in unambiguous assignments of all 8 (19)F NMR resonances. Our studies suggest that (19)F-edited NOESY NMR spectra are generally adequate for complete assignment without the need to resort to mutational analysis.

Specific coiled-coil interactions contribute to a global model of the structure of the spindle pole body

As the microtubule-organizing center of yeast, the spindle pole body (SPB) is essential for cell viability. Structural studies of the SPB are limited by its low copy number in the cell, its large size and heterogeneous composition, and its association with the nuclear membrane. However, low-resolution or indirect structural information about the SPB may be deciphered through a variety of techniques. Interestingly, a large proportion of SPB proteins are predicted to contain one or more coiled coils, a common protein interaction motif. The high frequency of coiled coils suggests that this structure is important for establishing the overall architecture of the complex. Support for this hypothesis was reported previously for coiled coils from some SPB proteins. Here, we extend this approach of isolating and characterizing additional SPB coiled coils to improve our understanding of SPB structure and organization. Self-associating coiled coils from Bbp1, Mps2, and Nbp1 were observed to form stable parallel homodimers in solution. Coiled-coil peptides from Bbp1 and Mps2 were also observed to hetero-associate. Experimental coiled-coil interaction data from this work and previous studies, as well as predicted and experimental structures for other SPB protein fragments and domains, were integrated to generate a model of the SPB structure.

Dynamics and energetics: a consensus analysis of the impact of calcium on EF-CaM protein complex

We have studied the relationship between dynamical correlations and energetic contributions in an attempt to model the transmission of information inside protein-protein complexes. The complex formed between the edema factor (EF) of Bacillus anthracis and calmodulin (CaM) was taken as an example, as the formation and stability of the complex depend on the calcium complexation level. The effect of calcium through EF-CaM residue network has been investigated with various approaches: 1), the elastic network model; 2), the local feature analysis; 3), the generalized correlations; and 4), the energetic dependency maps (EDMs), on 15-ns molecular dynamics simulations of the complex loaded with 0, 2, or 4 Ca2+ ions. The elastic network model correctly describes the basic architecture of the complex but is poorly sensitive to the level of calcium compared to the other methods. The local feature analysis allows us to characterize the local dynamics of the complex and the propagation of the calcium signal through CaM. The analyses of global dynamics and energetics--through generalized correlations and EDMs--provide a comprehensive picture of EF-CaM architecture and can be unified by using the concept of residue network connectedness. A medium connectedness, defined as the ability of each residue to communicate with all remaining parts of the complex, is observed for the 2Ca2+ level, which was experimentally identified as the most stable form of EF-CaM. The hierarchy of relative stabilities given by the EDMs sheds a new light on the EF-CaM interaction mechanism described experimentally and supports an organization of the complex architecture centered around nucleation points.

Analysis of coiled-coil interactions between core proteins of the spindle pole body

The spindle pole body (SPB) is a multiprotein complex that organizes microtubules in yeast. Due to its large size and association with the nuclear membrane, little is known about its detailed structure. In particular, although many SPB components and some of the interactions between them have been identified, the molecular details of how most of these interactions occur are not known. The prevalence of predicted coiled-coil regions in SPB proteins suggests that some interactions may occur via coiled coils. Here this hypothesis is supported by biochemical characterization of isolated coiled-coil peptides derived from SPB proteins. Formation of four strongly self-associating coiled-coil complexes from Spc29, Spc42, and Spc72 was demonstrated by circular dichroism (CD) spectroscopy and a fluorescence resonance energy transfer (FRET) assay. Many weaker self- and heteroassociations were also detected by CD, FRET, and/or cross-linking. The thermal stabilities of nine candidate homooligomers were assessed; six unfolded cooperatively with melting temperatures ranging from <11 to >50 degrees C. Solution studies established that coiled-coil peptides derived from Spc42 and Spc72 form parallel dimers, and this was confirmed for Spc42 by a high-resolution crystal structure. These data contribute to a growing body of knowledge that will ultimately provide a detailed model of the SPB structure.

A 1.3-A structure of zinc-bound N-terminal domain of calmodulin elucidates potential early ion-binding step

Calmodulin (CaM) is a 16.8-kDa calcium-binding protein involved in calcium-signal transduction. It is the canonical member of the EF-hand family of proteins, which are characterized by a helix-loop-helix calcium-binding motif. CaM is composed of N- and C-terminal globular domains (N-CaM and C-CaM), and within each domain there are two EF-hand motifs. Upon binding calcium, CaM undergoes a significant, global conformational change involving reorientation of the four helix bundles in each of its two domains. This conformational change upon ion binding is a key component of the signal transduction and regulatory roles of CaM, yet the precise nature of this transition is still unclear. Here, we present a 1.3-A structure of zinc-bound N-terminal calmodulin (N-CaM) solved by single-wavelength anomalous diffraction phasing of a selenomethionyl N-CaM. Our zinc-bound N-CaM structure differs from previously reported CaM structures and resembles calcium-free apo-calmodulin (apo-CaM), despite the zinc binding to both EF-hand motifs. Structural comparison with calcium-free apo-CaM, calcium-loaded CaM, and a cross-linked calcium-loaded CaM suggests that our zinc-bound N-CaM reveals an intermediate step in the initiation of metal ion binding at the first EF-hand motif. Our data also suggest that metal ion coordination by two key residues in the first metal-binding site represents an initial step in the conformational transition induced by metal binding. This is followed by reordering of the N-terminal region of the helix exiting from this first binding loop. This conformational switch should be incorporated into models of either stepwise conformational transition or flexible, dynamic energetic state sampling-based transition.

An Approximate Method in Using Molecular Mechanics Simulations To Study Slow Protein Conformational Changes

The broad range of characteristic motions in proteins has limited the applicability of molecular dynamics simulations in studying large-scale conformational transitions. We present an approximate method, making use of standard MD simulations and using a much larger integration time step, to obtain the structural changes for slow systematic motions of large complex systems. We show the applicability of this method by simulating the open to closed Calmodulin calcium binding domain conformational changes. Starting with the Ca2+-bound X-ray structure, and after the removal of the Ca2+ ions, our calculation yielded intermediate conformations during the rearrangement of helices in each Ca2+ binding pocket, leading to a structure with a lowest rmsd of 1.56 A compared to the NMR apo-calmodulin structure.

Calmodulin's flexibility allows for promiscuity in its interactions with target proteins and peptides

The small bilobal calcium regulatory protein calmodulin (CaM) activates numerous target enzymes in response to transient changes in intracellular calcium concentrations. Binding of calcium to the two helix-loop-helix calcium-binding motifs in each of the globular domains induces conformational changes that expose a methionine-rich hydrophobic patch on the surface of each domain of the protein, which it uses to bind to peptide sequences in its target enzymes. Although these CaM-binding domains typically have little sequence identity, the positions of several bulky hydrophobic residues are often conserved, allowing for classification of CaM-binding domains into recognition motifs, such as the 1-14 and 1-10 motifs. For calcium-independent binding of CaM, a third motif known as the IQ motif is also common. Many CaM-peptide complexes have globular conformations, where CaM's central linker connecting the two domains unwinds, allowing the protein to wrap around a single predominantly alpha-helical target peptide sequence. However, novel structures have recently been reported where the conformation of CaM is highly dissimilar to these globular complexes, in some instances with less than a full compliment of bound calcium ions, as well as novel stoichiometries. Furthermore, many divergent CaM isoforms from yeast and plant species have been discovered with unique calcium-binding and enzymatic activation characteristics compared to the single CaM isoform found in mammals.

Basic interdomain boundary residues in calmodulin decrease calcium affinity of sites I and II by stabilizing helix-helix interactions

Calmodulin is an EF-hand calcium-binding protein (148 a.a.) essential in intracellular signal transduction. Its homologous N- and C-terminal domains are separated by a linker that appears disordered in NMR studies. In a study of an N-domain fragment of Paramecium CaM (PCaM1-75), the addition of linker residues 76 to 80 (MKEQD) raised the Tm by 9 degrees C and lowered calcium binding by 0.54 kcal/mol (Sorensen et al., [Biochemistry 2002;41:15-20]), showing that these tether residues affect energetics as well as being a barrier to diffusion. To determine the individual contributions of residues 74 through 80 (RKMKEQD) to stability and calcium affinity, we compared a nested series of 7 fragments (PCaM1-74 to PCaM1-80). For the first 4, PCaM1-74 through PCaM1-77, single amino acid additions at the C-terminus corresponded to stepwise increases in thermostability and decreases in calcium affinity with a net change of 13.5 degrees C in Tm and 0.55 kcal/mol in free energy. The thermodynamic properties of fragments PCaM1-77 through PCaM1-80 were nearly identical. We concluded that the 3 basic residues in the sequence from 74 to 77 (RKMK) are critical to the increased stability and decreased calcium affinity of the longer N-domain fragments. Comparisons of NMR (HSQC) spectra of 15N-PCaM1-74 and 15N-PCaM1-80 and analysis of high-resolution structural models suggest these residues are latched to amino acids in helix A of CaM. The addition of residues E78, Q79, and D80 had a minimal effect on sites I and II, but they may contribute to the mechanism of energetic communication between the domains.

Identification of calmodulin isoform-specific binding peptides from a phage-displayed random 22-mer peptide library

Plants express numerous calmodulin (CaM) isoforms that exhibit differential activation or inhibition of CaM-dependent enzymes in vitro; however, their specificities toward target enzyme/protein binding are uncertain. A random peptide library displaying a 22-mer peptide on a bacteriophage surface was constructed to screen peptides that specifically bind to plant CaM isoforms (soybean calmodulin (ScaM)-1 and SCaM-4 were used in this study) in a Ca2+-dependent manner. The deduced amino acid sequence analyses of the respective 80 phage clones that were independently isolated via affinity panning revealed that SCaM isoforms require distinct amino acid sequences for optimal binding. SCaM-1-binding peptides conform to a 1-5-10 ((FILVW)XXX(FILV) XXXX(FILVW)) motif (where X denotes any amino acid), whereas SCaM-4-binding peptide sequences conform to a 1-8-14 ((FILVW)XXXXXX(FAILVW)XXXXX(FILVW)) motif. These motifs are classified based on the positions of conserved hydrophobic residues. To examine their binding properties further, two representative peptides from each of the SCaM isoform-binding sequences were synthesized and analyzed via gel mobility shift assays, Trp fluorescent spectra analyses, and phosphodiesterase competitive inhibition experiments. The results of these studies suggest that SCaM isoforms possess different binding sequences for optimal target interaction, which therefore may provide a molecular basis for CaM isoform-specific function in plants. Furthermore, the isolated peptide sequences may serve not only as useful CaM-binding sequence references but also as potential reagents for studying CaM isoform-specific function in vivo.