PAIRpred: Partner‐specific prediction of interacting residues from sequence and structure

We present a novel partner‐specific protein–protein interaction site prediction method called PAIRpred. Unlike most existing machine learning binding site prediction methods, PAIRpred uses information from both proteins in a protein complex to predict pairs of interacting residues from the two proteins. PAIRpred captures sequence and structure information about residue pairs through pairwise kernels that are used for training a support vector machine classifier. As a result, PAIRpred presents a more detailed model of protein binding, and offers state of the art accuracy in predicting binding sites at the protein level as well as inter‐protein residue contacts at the complex level. We demonstrate PAIRpred's performance on Docking Benchmark 4.0 and recent CAPRI targets. We present a detailed performance analysis outlining the contribution of different sequence and structure features, together with a comparison to a variety of existing interface prediction techniques. We have also studied the impact of binding‐associated conformational change on prediction accuracy and found PAIRpred to be more robust to such structural changes than existing schemes. As an illustration of the potential applications of PAIRpred, we provide a case study in which PAIRpred is used to analyze the nature and specificity of the interface in the interaction of human ISG15 protein with NS1 protein from influenza A virus. Python code for PAIRpred is available at http://combi.cs.colostate.edu/supplements/pairpred/ . Proteins 2014; 82:1142–1155. © 2013 Wiley Periodicals, Inc.     

The complex structure of calmodulin bound to a calcineurin peptide

The activity of the protein phosphatase calcineurin (CN) is regulated by an autoinhibition mechanism wherein several domains from its catalytic A subunit, including the calmodulin binding domain (CaMBD), block access to its active site. Upon binding of Ca2+ and calmodulin (Ca2+/CaM) to CaMBD, the autoinhibitory domains dissociate from the catalytic groove, thus activating the enzyme. To date, the structure of the CN/CaM/Ca2+ complex has not been determined in its entirety. Previously, we determined the structure of a fusion protein consisting of CaM and a 25-residue peptide taken from the CaMBD, joined by a 5-glycine linker. This structure revealed a novel CaM binding motif. However, the presence of the extraneous glycine linker cast doubt on the authenticity of this structure as an accurate representation of CN/CaM binding in vivo. Thus, here, we have determined the crystal structure of CaM complexed with the 25-residue CaMBD peptide without the glycine linker at a resolution of 2.1 A. The structure is essentially identical to the fusion construction which displays CaM bound to the CaMBD peptide as a dimer with an open, elongated conformation. The N-lobe from one molecule and C-lobe from another encompass and bind the CaMBD peptide. Thus, it validates the existence of this novel CaM binding motif. Our experiments suggest that the dimeric CaM/CaMBD complex exists in solution, which is unambiguously validated using a carefully-designed CaM-sepharose pull-down experiment. We discuss structural features that produce this novel binding motif, including the role of the CaMBD peptide residues Arg-408, Val-409, and Phe-410, which work to provide rigidity to the otherwise flexible central CaM helix joining the N- and C-lobes, ultimately keeping these lobes apart and forcing "head-to-tail" dimerization to attain the requisite N- and C-lobe pairing for CaMBD binding.     

CaMELS: In silico prediction of calmodulin binding proteins and their binding sites

Due to Ca²⁺‐dependent binding and the sequence diversity of Calmodulin (CaM) binding proteins, identifying CaM interactions and binding sites in the wet‐lab is tedious and costly. Therefore, computational methods for this purpose are crucial to the design of such wet‐lab experiments. We present an algorithm suite called CaMELS (CalModulin intEraction Learning System) for predicting proteins that interact with CaM as well as their binding sites using sequence information alone. CaMELS offers state of the art accuracy for both CaM interaction and binding site prediction and can aid biologists in studying CaM binding proteins. For CaM interaction prediction, CaMELS uses protein sequence features coupled with a large‐margin classifier. CaMELS models the binding site prediction problem using multiple instance machine learning with a custom optimization algorithm which allows more effective learning over imprecisely annotated CaM‐binding sites during training. CaMELS has been extensively benchmarked using a variety of data sets, mutagenic studies, proteome‐wide Gene Ontology enrichment analyses and protein structures. Our experiments indicate that CaMELS outperforms simple motif‐based search and other existing methods for interaction and binding site prediction. We have also found that the whole sequence of a protein, rather than just its binding site, is important for predicting its interaction with CaM. Using the machine learning model in CaMELS, we have identified important features of protein sequences for CaM interaction prediction as well as characteristic amino acid sub‐sequences and their relative position for identifying CaM binding sites. Python code for training and evaluating CaMELS together with a webserver implementation is available at the URL: http://faculty.pieas.edu.pk/fayyaz/software.html#camels .     

Identification of the critical structural determinants of the EF-hand domain arrangements in calcium binding proteins

EF-hand calcium binding proteins (CaBPs) share strong sequence homology, but exhibit great diversity in structure and function. Thus although calmodulin (CaM) and calcineurin B (CNB) both consist of four EF hands, their domain arrangements are quite distinct. CaM and the CaM-like proteins are characterized by an extended architecture, whereas CNB and the CNB-like proteins have a more compact form. In this study, we performed structural alignments and molecular dynamics (MD) simulations on 3 CaM-like proteins and 6 CNB-like proteins, and quantified their distinct structural and dynamical features in an effort to establish how their sequences specify their structures and dynamics. Alignments of the EF2-EF3 region of these proteins revealed that several residues (not restricted to the linker between the EF2 and EF3 motifs) differed between the two groups of proteins. A customized inverse folding approach followed by structural assessments and MD simulations established the critical role of these residues in determining the structure of the proteins. Identification of the critical determinants of the two different EF-hand domain arrangements and the distinct dynamical features relevant to their respective functions provides insight into the relationships between sequence, structure, dynamics and function among these EF-hand CaBPs.     

Protein flexibility is key to cisplatin crosslinking in calmodulin

Chemical crosslinking in combination with Fourier transform ion cyclotron resonance mass spectrometry (FTICR MS) has significant potential for studying protein structures and protein-protein interactions. Previously, cisplatin has been shown to be a crosslinker and crosslinks multiple methionine (Met) residues in apo-calmodulin (apo-CaM). However, the inter-residue distances obtained from nuclear magnetic resonance structures are inconsistent with the measured distance constraints by crosslinking. Met residues lie too far apart to be crosslinked by cisplatin. Here, by combining FTICR MS with a novel computational flexibility analysis, the flexible nature of the CaM structure is found to be key to cisplatin crosslinking in CaM. It is found that the side chains of Met residues can be brought together by flexible motions in both apo-CaM and calcium-bound CaM (Ca(4) -CaM). The possibility of cisplatin crosslinking Ca(4) -CaM is then confirmed by MS data. Therefore, flexibility analysis as a fast and low-cost computational method can be a useful tool for predicting crosslinking pairs in protein crosslinking analysis and facilitating MS data analysis. Finally, flexibility analysis also indicates that the crosslinking of platinum to pairs of Met residues will effectively close the nonpolar groove and thus will likely interfere with the binding of CaM to its protein targets, as was proved by comparing assays for cisplatin-modified/unmodified CaM binding to melittin. Collectively, these results suggest that cisplatin crosslinking of apo-CaM or Ca(4) -CaM can inhibit the ability of CaM to recognize its target proteins, which may have important implications for understanding the mechanism of tumor resistance to platinum anticancer drugs.     

Fast methionine-based solution structure determination of calcium-calmodulin complexes

Here we present a novel NMR method for the structure determination of calcium-calmodulin (Ca²⁺-CaM)-peptide complexes from a limited set of experimental restraints. A comparison of solved CaM-peptide structures reveals invariability in CaM’s backbone conformation and a structural plasticity in CaM’s domain orientation enabled by a flexible linker. Knowing this, the collection and analysis of an extensive set of NOESY spectra is redundant. Although RDCs can define CaM domain orientation in the complex, they lack the translational information required to position the domains on the bound peptide and highlight the necessity of intermolecular NOEs. Here we employ a specific isotope labeling strategy in which the role of methionine in CaM-peptide interactions is exploited to collect these critical NOEs. By ¹H, ¹³C-labeling the methyl groups of deuterated methionine against a ²H, ¹²C background, we can acquire a ¹³C-edited NOESY characterized by simplified, easily analyzable spectra. Together with measured CaM backbone H_N-N RDCs and intrapeptide NOE-based distances, these intermolecular NOEs provide restraints for a low temperature torsion-angle dynamics and simulated annealing protocol used to calculate the complex structure. We have applied our method to a CaM complex previously solved through X-ray crystallography: Ca²⁺-CaM bound to the CaM kinase I peptide (PDB code: 1MXE). The resulting structure has a backbone RMSD of 1.6 Å to that previously published. We have also used this test complex to investigate the importance of homologous model selection on the calculated outcome. In addition to having application for fast complex structure determination, this method can be used to determine the structures of difficult complexes characterized by chemical shift overlap and broad signals for which the traditional method based on the use of fully ¹³C, ¹⁵N-labeled CaM fails.     

Local structure prediction with local structure-based sequence profiles

MOTIVATION: A large body of experimental and theoretical evidence suggests that local structural determinants are frequently encoded in short segments of protein sequence. Although the local structural information, once recognized, is particularly useful in protein structural and functional analyses, it remains a difficult problem to identify embedded local structural codes based solely on sequence information. RESULTS: In this paper, we describe a local structure prediction method aiming at predicting the backbone structures of nine-residue sequence segments. Two elements are the keys for this local structure prediction procedure. The first key element is the LSBSP1 database, which contains a large number of non-redundant local structure-based sequence profiles for nine-residue structure segments. The second key element is the consensus approach, which identifies a consensus structure from a set of hit structures. The local structure prediction procedure starts by matching a query sequence segment of nine consecutive amino acid residues to all the sequence profiles in the local structure-based sequence profile database (LSBSP1). The consensus structure, which is at the center of the largest structural cluster of the hit structures, is predicted to be the native state structure adopted by the query sequence segment. This local structure prediction method is assessed with a large set of random test protein structures that have not been used in constructing the LSBSP1 database. The benchmark results indicate that the prediction capacities of the novel local structure prediction procedure exceed the prediction capacities of the local backbone structure prediction methods based on the I-sites library by a significant margin. AVAILABILITY: All the computational and assessment procedures have been implemented in the integrated computational system PrISM.1 (Protein Informatics System for Modeling). The system and associated databases for LINUX systems can be downloaded from the website: http://www.columbia.edu/~ay1/.     

Solution X-ray scattering reveals a novel structure of calmodulin complexed with a binding domain peptide from the HIV-1 matrix protein p17

The solution structures of complexes between calcium-saturated calmodulin (Ca (2+)/CaM) and a CaM-binding domain of the HIV-1 matrix protein p17 have been determined by small-angle X-ray scattering with use of synchrotron radiation as an intense and stable X-ray source. We used three synthetic peptides of residues 11-28, 26-47, and 11-47 of p17 to demonstrate the diversity of CaM-binding conformation. Ca (2+)/CaM complexed with residues 11-28 of p17 adopts a dumbbell-like structure at a molar ratio of 1:2, suggesting that the two peptides bind each lobe of CaM, respectively. Ca (2+)/CaM complexed with residues 26-47 of p17 at a molar ratio of 1:1 adopts a globular structure similar to the NMR structure of Ca (2+)/CaM bound to M13, which adopted a compact globular structure. In contrast to these complexes, Ca (2+)/CaM binds directly with both CaM-binding sites of residues 11-47 of p17 at a molar ratio of 1:1, which induces a novel structure different from known structures previously reported between Ca (2+)/CaM and peptide. A tertiary structural model of the novel structure was constructed using the biopolymer module of Insight II 2000 on the basis of the scattering data. The two domains of CaM remain essentially unchanged upon complexation. The hinge motions, however, occur in a highly flexible linker of CaM, in which the electrostatic residues 74Arg, 78Asp, and 82Glu interact with N-terminal electrostatic residues of the peptide (residues 12Glu, 15Arg, and 18Lys). The acidic residues in the N-terminal domain of CaM interact with basic residues in a central part of the peptide, thereby enabling the central part to change the conformations, while an acidic residue in the C-terminal domain interacts with two basic residues in the two helical sites of the peptide. The overall structure of the complex adopts an extended structure with the radius of gyration of 20.5 A and the interdomain distance of 34.2 A. Thus, the complex is principally stabilized by electrostatic interactions. The hydrophobic patches of Ca (2+)/CaM are not responsible for the binding with the hydrophobic residues in the peptide, suggesting that CaM plays a role to sequester the myristic acid moiety of p17.     

Seventy-five percent accuracy in protein secondary structure prediction

In this study we present an accurate secondary structure prediction procedure by using a query and related sequences. The most novel aspect of our approach is its reliance on local pairwise alignment of the sequence to be predicted with each related sequence rather than utilization of a multiple alignment. The residue-by-residue accuracy of the method is 75% in three structural states after jack-knife tests. The gain in prediction accuracy compared with the existing techniques, which are at best 72%, is achieved by secondary structure propensities based on both local and long-range effects, utilization of similar sequence information in the form of carefully selected pairwise alignment fragments, and reliance on a large collection of known protein primary structures. The method is especially appropriate for large-scale sequence analysis efforts such as genome characterization, where precise and significant multiple sequence alignments are not available or achievable. Proteins 27:329–335, 1997. © 1997 Wiley-Liss, Inc.     

A Novel Trans Conformation of Ligand-Free Calmodulin

Calmodulin (CaM) is a highly conserved eukaryotic protein that binds specifically to more than 100 target proteins in response to calcium (Ca²⁺) signal. CaM adopts a considerable degree of structural plasticity to accomplish this physiological role; however, the nature and extent of this plasticity remain to be fully understood. Here, we report the crystal structure of a novel trans conformation of ligand-free CaM where the relative disposition of two lobes of CaM is different, a conformation to-date not reported. While no major structural changes were observed in the independent N- and C-lobes as compared with previously reported structures of Ca²⁺/CaM, the central helix was tilted by ∼90° at Arg75. This is the first crystal structure of CaM to show a drastic conformational change in the central helix, and reveals one of several possible conformations of CaM to engage with its binding partner.     

Position-specific residue preference features around the ends of helices and strands and a novel strategy for the prediction of secondary structures

It has been many years since position-specific residue preference around the ends of a helix was revealed. However, all the existing secondary structure prediction methods did not exploit this preference feature, resulting in low accuracy in predicting the ends of secondary structures. In this study, we collected a relatively large data set consisting of 1860 high-resolution, non-homology proteins from the PDB, and further analyzed the residue distributions around the ends of regular secondary structures. It was found that there exist position-specific residue preferences (PSRP) around the ends of not only helices but also strands. Based on the unique features, we proposed a novel strategy and developed a tool named E-SSpred that treats the secondary structure as a whole and builds models to predict entire secondary structure segments directly by integrating relevant features. In E-SSpred, the support vector machine (SVM) method is adopted to model and predict the ends of helices and strands according to the unique residue distributions around them. A simple linear discriminate analysis method is applied to model and predict entire secondary structure segments by integrating end-prediction results, tri-peptide composition, and length distribution features of secondary structures, as well as the prediction results of the most famous program PSIPRED. The results of fivefold cross-validation on a widely used data set demonstrate that the accuracy of E-SSpred in predicting ends of secondary structures is about 10% higher than PSIPRED, and the overall prediction accuracy (Q(3) value) of E-SSpred (82.2%) is also better than PSIPRED (80.3%). The E-SSpred web server is available at http://bioinfo.hust.edu.cn/bio/tools/E-SSpred/index.html.     

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

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

Munc13-Like skMLCK Variants Cannot Mimic the Unique Calmodulin Binding Mode of Munc13 as Evidenced by Chemical Cross-Linking and Mass Spectrometry

Among the neuronal binding partners of calmodulin (CaM) are Munc13 proteins as essential presynaptic regulators that play a key role in synaptic vesicle priming and are crucial for presynaptic short-term plasticity. Recent NMR structural investigations of a CaM/Munc13-1 peptide complex have revealed an extended structure, which contrasts the compact structures of most classical CaM/target complexes. This unusual binding mode is thought to be related to the presence of an additional hydrophobic anchor residue at position 26 of the CaM binding motif of Munc13-1, resulting in a novel 1-5-8-26 motif. Here, we addressed the question whether the 1-5-8-26 CaM binding motif is a Munc13-related feature or whether it can be induced in other CaM targets by altering the motif''s core residues. For this purpose, we chose skeletal muscle myosin light chain kinase (skMLCK) with a classical 1-5-8-14 CaM binding motif and constructed three skMLCK peptide variants mimicking Munc13-1, in which the hydrophobic anchor amino acid at position 14 was moved to position 26. Chemical cross-linking between CaM and skMLCK peptide variants combined with high-resolution mass spectrometry yielded insights into the peptides'' binding modes. This structural comparison together with complementary binding data from surface plasmon resonance experiments revealed that skMLCK variants with an artificial 1-5-8-26 motif cannot mimic CaM binding of Munc13-1. Apparently, additional features apart from the spacing of the hydrophobic anchor residues are required to define the functional 1-5-8-26 motif of Munc13-1. We conclude that Munc13 proteins display a unique CaM binding behavior to fulfill their role as efficient presynaptic calcium sensors over broad range of Ca²⁺ concentrations.     

Target-induced conformational adaptation of calmodulin revealed by the crystal structure of a complex with nematode Ca(2+)/calmodulin-dependent kinase kinase peptide.

Calmodulin (CaM) is a ubiquitous calcium (Ca(2+)) sensor which binds and regulates protein serine/threonine kinases along with many other proteins in a Ca(2+)-dependent manner. For this multi-functionality, conformational plasticity is essential; however, the nature and magnitude of CaM's plasticity still remains largely undetermined. Here, we present the 1.8 A resolution crystal structure of Ca(2+)/CaM, complexed with the 27-residue synthetic peptide corresponding to the CaM-binding domain of the nematode Caenorhabditis elegans Ca(2+)/CaM-dependent kinase kinase (CaMKK). The peptide bound in this crystal structure is a homologue of the previously NMR-derived complex with rat CaMKK, but benefits from improved structural resolution. Careful comparison of the present structure to previous crystal structures of CaM complexed with unrelated peptides derived from myosin light chain kinase and CaM kinase II, allow a quantitative analysis of the differences in the relative orientation of the N and C-terminal domains of CaM, defined as a screw axis rotation angle ranging from 156 degrees to 196 degrees. The principal differences in CaM interaction with various peptides are associated with the N-terminal domain of CaM. Unlike the C-terminal domain, which remains unchanged internally, the N-terminal domain of CaM displays significant differences in the EF-hand helix orientation between this and other CaM structures. Three hydrogen bonds between CaM and the peptide (E87-R336, E87-T339 and K75-T339) along with two salt bridges (E11-R349 and E114-K334) are the most probable determinants for the binding direction of the CaMKK peptide to CaM.     

Myristoylation-regulated direct interaction between calcium-bound calmodulin and N-terminal region of pp60v-src

pp60v-src tyrosine protein kinase was suggested to interact with Ca2+-bound calmodulin (Ca2+/CaM) through the N-terminal region based on its structural similarities to CAP-23/NAP-22, a myristoylated neuron-specific protein, whose myristoyl group is essential for interaction with Ca2+/CaM; (1) the N terminus of pp60v-src is myristoylated like CAP-23/NAP-22; (2) both lysine residues are required for the myristoylation-dependent interaction and serine residues that are thought to regulate the interaction through the phosphorylations located in the N-terminal region of pp60v-src. To verify this possibility, we investigated the direct interaction between pp60v-src and Ca2+/CaM using a myristoylated peptide corresponding to the N-terminal region of pp60v-src. The binding assay indicated that only the myristoylated peptide binds to Ca2+/CaM, and the non-myristoylated peptide is not able to bind to Ca2+/CaM. Analyses of the binding kinetics revealed two independent reactions with the dissociation constants (KD) of 2.07 x 10(-9)M (KD1) and 3.93 x 10(-6)M (KD2), respectively. Two serine residues near the myristoyl moiety of the peptide (Ser2, Ser11) were phosphorylated by protein kinase C in vitro, and the phosphorylation drastically reduced the interaction. NMR experiments indicated that two molecules of the myristoylated peptide were bound around the hydrophobic clefts of a Ca2+/CaM molecule. The small-angle X-ray scattering analyses showed that the size of the peptide-Ca2+/CaM complex is 2-3A smaller than that of the known Ca2+/CaM-target molecule complexes. These results demonstrate clearly the direct interaction between pp60v-src and Ca2+/CaM in a novel manner different from that of known Ca2+/CaM, the target molecules, interactions.     

The binding of myristoylated N-terminal nonapeptide from neuro-specific protein CAP-23/NAP-22 to calmodulin does not induce the globular structure observed for the calmodulin-nonmyristylated peptide complex

CAP-23/NAP-22, a neuron-specific protein kinase C substrate, is Nalpha-myristoylated and interacts with calmodulin (CaM) in the presence of Ca2+ ions. Takasaki et al. (1999, J Biol Chem 274:11848-11853) have recently found that the myristoylated N-terminal nonapeptide of CAP-23/NAP-22 (mC/N9) binds to Ca2+ -bound CaM (Ca2+/CaM). In the present study, small-angle X-ray scattering was used to investigate structural changes of Ca2+/CaM induced by its binding to mC/N9 in solution. The binding of one mC/N9 molecule induced an insignificant structural change in Ca2+/CaM. The 1:1 complex appeared to retain the extended conformation much like that of Ca2+/CaM in isolation. However, it could be seen that the binding of two mC/N9 molecules induced a drastic structural change in Ca2+/CaM, followed by a slight structural change by the binding of more than two but less than four mC/N9 molecules. Under the saturated condition (the molar ratio of 1:4), the radius of gyration (Rg) for the Ca2+/CaM-mC/N9 complex was 19.8 +/- 0.3 A. This value was significantly smaller than that of Ca2+/CaM (21.9 +/- 0.3 A), which adopted a dumbbell structure and was conversely 2-3 A larger than those of the complexes of Ca2+/CaM with the nonmyristoylated target peptides of myosin light chain kinase or CaM kinase II, which adopted a compact globular structure. The pair distance distribution function had no shoulder peak at around 40 A, which was mainly due to the dumbbell structure. These results suggest that Ca2+/CaM interacts with Nalpha-myristoylated CAP-23/NAP-22 differently than it does with other nonmyristoylated target proteins. The N-terminal amino acid sequence alignment of CAP-23/NAP-22 and other myristoylated proteins suggests that the protein myristoylation plays important roles not only in the binding of CAP-23/NAP-22 to Ca2+/CaM, but also in the protein-protein interactions related to other myristoylated proteins.     

A structural comparison of ‘real’ and ‘model’ calmodulin clarified allosteric interactions regulating domain motion

Calmodulin (CaM) is a multifunctional calcium-binding protein, which regulates various biochemical processes. CaM acts via structural changes and complex forming with its target enzymes. CaM has two globular domains (N-lobe and C-lobe) connected by a long linker region. Upon calcium binding, the N-lobe and C-lobe undergo local conformational changes, after that, entire CaM wraps the target enzyme through a large conformational change. However, the regulation mechanism, such as allosteric interactions regulating the conformational changes, is still unclear. In order to clarify the allosteric interactions, in this study, experimentally obtained ‘real’ structures are compared to ‘model’ structures lacking the allosteric interactions. As the allosteric interactions would be absent in calcium-free CaM (apo-CaM), allostery-eliminated calcium-bound CaM (holo-CaM) models were constructed by combining the apo-CaM’s linker and the holo-CaM’s N- and C-lobe. Before the comparison, the ‘real’ and ‘model’ structures were clustered and cluster–cluster relationship was determined by a principal component analysis. The structures were compared based on the relationship, then, a distance map and a contact probability analysis clarified that the inter-domain motion is regulated by several groups of inter-domain contacting residue pairs. The analyses suggested that these residues cause inter-domain translation and rotation, and as a consequence, the motion encourage structural diversity. The resultant diversity would contribute to the functional versatility of CaM.     

Different conformational switches underlie the calmodulin-dependent modulation of calcium pumps and channels

We have used fluorescence spectroscopy to investigate the structure of calmodulin (CaM) bound with CaM-binding sequences of either the plasma membrane Ca-ATPase or the skeletal muscle ryanodine receptor (RyR1) calcium release channel. Following derivatization with N-(1-pyrene)maleimide at engineered sites (T34C and T110C) within the N- and C-domains of CaM, contact interactions between these opposing domains of CaM resulted in excimer fluorescence that permits us to monitor conformational states of bound CaM. Complementary measurements take advantage of the unique conserved Trp within CaM-binding sequences that functions as a hydrophobic anchor in CaM binding and permits measurements of both a local and global peptide structure. We find that CaM binds with high affinity in a collapsed structure to the CaM-binding sequences of both the Ca-ATPase and RyR1, resulting in excimer formation that is indicative of contact interactions between the N- and the C-domains of CaM in complex with these CaM-binding peptides. There is a 4-fold larger amount of excimer formation for CaM bound to the CaM-binding sequence of the Ca-ATPase in comparison to RyR1, indicating a closer structural coupling between CaM domains in this complex. Prior to CaM association, the CaM-binding sequences of the Ca-ATPase and RyR1 are conformationally disordered. Upon CaM association, the CaM-binding sequence of the Ca-ATPase assumes a highly ordered structure. In comparison, the CaM-binding sequence of RyR1 remains conformationally disordered irrespective of CaM binding. These results suggest an important role for interdomain contact interactions between the opposing domains of CaM in stabilizing the structure of the peptide complex. The substantially different structural responses associated with CaM binding to Ca-ATPase and RyR1 indicates a plasticity in their respective binding mechanisms that accomplishes different physical mechanisms of allosteric regulation, involving either the dissociation of a C-terminal regulatory domain necessary for pump activation or the modulation of intersubunit interactions to diminish RyR1 channel activity.     

Bayesian segmentation of protein secondary structure.

We present a novel method for predicting the secondary structure of a protein from its amino acid sequence. Most existing methods predict each position in turn based on a local window of residues, sliding this window along the length of the sequence. In contrast, we develop a probabilistic model of protein sequence/structure relationships in terms of structural segments, and formulate secondary structure prediction as a general Bayesian inference problem. A distinctive feature of our approach is the ability to develop explicit probabilistic models for alpha-helices, beta-strands, and other classes of secondary structure, incorporating experimentally and empirically observed aspects of protein structure such as helical capping signals, side chain correlations, and segment length distributions. Our model is Markovian in the segments, permitting efficient exact calculation of the posterior probability distribution over all possible segmentations of the sequence using dynamic programming. The optimal segmentation is computed and compared to a predictor based on marginal posterior modes, and the latter is shown to provide significant improvement in predictive accuracy. The marginalization procedure provides exact secondary structure probabilities at each sequence position, which are shown to be reliable estimates of prediction uncertainty. We apply this model to a database of 452 nonhomologous structures, achieving accuracies as high as the best currently available methods. We conclude by discussing an extension of this framework to model nonlocal interactions in protein structures, providing a possible direction for future improvements in secondary structure prediction accuracy.