Energetics of the native energy landscape of a two-domain calcium sensor protein: distinct folding features of the two domains

The protein folding energy landscape allows a thorough understanding of the protein folding problem which in turn helps in understanding various aspects of biological functions. Characterizing the cooperative unfolding units and the intermediates along the folding funnel of a protein is a challenging task. In this paper, we investigated the native energy landscape of EhCaBP, a calcium sensor, belonging to the same EF-hand superfamily as calmodulin. EhCaBP is a two-domain EF-hand protein consisting of two EF-hands in each domain and binding to four Ca2+ cations. Native-state hydrogen exchange (HX) was used to assess the folding features of the landscape and also to throw light on the structure-folding function paradigm of calcium sensor proteins. HX measurements under the EX2 regime provided the thermodynamic information about the protein folding events under native conditions. HX studies revealed that the unfolding of EhCaBP is not a two-state process. Instead, it proceeds through cooperative units. The C-terminal domain exhibits less denaturant dependence than the N-terminal domain, suggesting that the former is dominated by local fluctuations. It is interesting to note that the N- and C-terminal domains of EhCaBP have distinct folding features. In fact, these observed differences can regulate the domain-dependent target recognition of two-domain Ca2+ sensor proteins.     

Dynamics alignment: comparison of protein dynamics in the SCOP database

A novel methodology for comparison of protein dynamics is presented. Protein dynamics is calculated using the Gaussian network model and the modes of motion are globally aligned using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. The alignment is fast and can be used to analyze large sets of proteins. The methodology is applied to the four major classes of the SCOP database: "all alpha proteins," "all beta proteins," "alpha and beta proteins," and "alpha/beta proteins". We show that different domains may have similar global dynamics. In addition, we report that the dynamics of "all alpha proteins" domains are less specific to structural variations within a given fold or superfamily compared with the other classes. We report that domain pairs with the most similar and the least similar global dynamics tend to be of similar length. The significance of the methodology is that it suggests a new and efficient way of mapping between the global structural features of protein families/subfamilies and their encoded dynamics.     

Interaction interfaces of protein domains are not topologically equivalent across families within superfamilies: Implications for metabolic and signaling pathways

Using a data set of aligned protein domain superfamilies of known three-dimensional structure, we compared the location of interdomain interfaces on the tertiary folds between members of distantly related protein domain superfamilies. The data set analyzed is comprised of interdomain interfaces, with domains occurring within a polypeptide chain and those between two polypeptide chains. We observe that, in general, the interfaces between protein domains are formed entirely in different locations on the tertiary folds in such pairs. This variation in the location of interface happens in protein domains involved in a wide range of functions, such as enzymes, adapters, and domains that bind protein ligands, or cofactors. While basic biochemical functionality is preserved at the domain superfamily level, the effect of biochemical function on protein assemblies is different in these protein domains related by superfamily. The divergence between proteins, in most cases, is coupled with domain recruitment, with different modes of interaction with the recruited domain. This is in complete contrast to the observation that in closely related homologous protein domains, almost always the interaction interfaces are topologically equivalent. In a small subset of interacting domains within proteins related by remote homology, we observe that the relative positioning of domains with respect to one another is preserved. Based on the analysis of multidomain proteins of known or unknown structure, we suggest that variation in protein-protein interactions in members within a superfamily could serve as diverging points in otherwise parallel metabolic or signaling pathways. We discuss a few representative cases of diverging pathways involving domains in a superfamily.     

Diversity of conformational states and changes within the EF-hand protein superfamily

The EF-hand motif, which assumes a helix-loop-helix structure normally responsible for Ca2+ binding, is found in a large number of functionally diverse Ca2+ binding proteins collectively known as the EF-hand protein superfamily. In many superfamily members, Ca2+ binding induces a conformational change in the EF-hand motif, leading to the activation or inactivation of target proteins. In calmodulin and troponin C, this is described as a change from the closed conformational state in the absence of Ca2+ to the open conformational state in its presence. It is now clear from structures of other EF-hand proteins that this "closed-to-open" conformational transition is not the sole model for EF-hand protein structural response to Ca2+. More complex modes of conformational change are observed in EF-hand proteins that interact with a covalently attached acyl group (e.g., recoverin) and in those that dimerize (e.g., S100B, calpain). In fact, EF-hand proteins display a multitude of unique conformational states, together constituting a conformational continuum. Using a quantitative 3D approach termed vector geometry mapping (VGM), we discuss this tertiary structural diversity of EF-hand proteins and its correlation with target recognition.     

Probing the S100 protein family through genomic and functional analysis

The EF-hand superfamily of calcium binding proteins includes the S100, calcium binding protein, and troponin subfamilies. This study represents a genome, structure, and expression analysis of the S100 protein family, in mouse, human, and rat. We confirm the high level of conservation between mammalian sequences but show that four members, including S100A12, are present only in the human genome. We describe three new members of the S100 family in the three species and their locations within the S100 genomic clusters and propose a revised nomenclature and phylogenetic relationship between members of the EF-hand superfamily. Two of the three new genes were induced in bone-marrow-derived macrophages activated with bacterial lipopolysaccharide, suggesting a role in inflammation. Normal human and murine tissue distribution profiles indicate that some members of the family are expressed in a specific manner, whereas others are more ubiquitous. Structure-function analysis of the chemotactic properties of murine S100A8 and human S100A12, particularly within the active hinge domain, suggests that the human protein is the functional homolog of the murine protein. Strong similarities between the promoter regions of human S100A12 and murine S100A8 support this possibility. This study provides insights into the possible processes of evolution of the EF-hand protein superfamily. Evolution of the S100 proteins appears to have occurred in a modular fashion, also seen in other protein families such as the C2H2-type zinc-finger family.     

Characterization of Neurospora crassa ��-Actinin

α-Actinin, an actin-binding protein of the spectrin superfamily, is present in most eukaryotes except plants. It is composed of three domains: N-terminal CH-domains, C-terminal calcium-binding domain (with EF-hand motifs), and a central rod domain. We have cloned and expressed Neurospora crassa α-actinin as GST and GFP fusion proteins for biochemical characterization and in vivo localization, respectively. The intracellular localization pattern of α-actinin suggests that this protein is intimately associated with actin filaments and plays an important role in the processes of germination, hyphal elongation, septum formation, and conidiation. These functions were confirmed by the experiments on the effect of α-actinin gene deletion in N. crassa.     

Transconformations of the SERCA1 Ca-ATPase: a normal mode study

The transport of Ca(2+) by Ca-ATPase across the sarcoplasmic reticulum membrane is accompanied by several transconformations of the protein. Relying on the already established functional importance of low-frequency modes in dynamics of proteins, we report here a normal mode analysis of the Ca(2+)-ATPase based on the crystallographic structures of the E1Ca(2) and E2TG forms. The lowest-frequency modes reveal that the N and A(+Nter) domains undergo the largest amplitude movements. The dynamical domain analysis performed with the DomainFinder program suggests that they behave as rigid bodies, unlike the highly flexible P domain. We highlight two types of movements of the transmembrane helices: i), a concerted movement around an axis perpendicular to the membrane which "twists open" the lumenal side of the protein and ii), an individual translational and rotational mobility which is of lower amplitude for the helices hosting the calcium binding sites. Among all modes calculated for E1Ca, only three are enough to describe the transition to E2TG; the associated movements involve almost exclusively the A and N domains, reflecting the closure of the cytoplasmic headpiece and high displacement of the L7-8 lumenal loop. Subsequently, we discuss the potential contribution of the remaining low-frequency normal modes to the transconformations occurring within the overall calcium transport cycle.     

Amino acid sequence of the Ca2(+)-triggered luciferin binding protein of Renilla reniformis

The complete amino acid sequence of the Ca2(+)-triggered luciferin binding protein (LBP) of Renilla reniformis has been determined. The apoprotein has an unblocked amino terminus and contains 184 residues with a calculated Mr of 20,541. LBP is a member of the EF-hand superfamily of Ca2(+)-binding proteins and bears three predicted EF-hand domains. The sequence and organization of EF-hand domains are similar to those of the Ca2(+)-dependent photoprotein, aequorin.     

Domain stability and metal-induced folding of calcium- and integrin-binding protein 1

It is widely accepted that a pair of EF-hands is the functional unit of typical four EF-hand proteins such as calmodulin or troponin C. In this work we investigate the structure and stability of the four EF-hand domains in the related protein calcium- and integrin-binding protein 1 (CIB1) in the presence and absence of Mg2+ or Ca2+, to determine if similar EF-hand interactions occur. The backbone structure and flexibility of CIB1 were first studied by NMR spectroscopy, and these studies were complimented with steady-state fluorescence spectroscopy and chemical denaturation experiments using mutant CIB1 proteins having single Trp reporter groups in each of the four EF-hand domains EF-I (F34W), EF-II (F91W), EF-III (L128W), and EF-IV (F173W). We find that Mg2+-CIB1 adopts a well-folded structure similar to Ca2+-CIB1, except for some conformational heterogeneity in the C-terminal EF-IV domain. The structure of apo-CIB1 is significantly more dynamic, especially within EF-II, EF-III, and a partially unfolded EF-IV region, but the N-terminal EF-I region of apo-CIB1 has a well-ordered and more stable structure. The data reveal significant communication between the N- and C-lobes of CIB1, and show that transient intermediate conformations are formed along the unfolding pathway for each form of the protein. Collectively the data demonstrate that the communication between the paired EF-hand domains as well as between the N- and C-lobes of CIB1 is distinct from the ancestral proteins calmodulin and troponin C, which might be important for the unique function of CIB1 in numerous biological processes.     

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.     

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.     

Microbial βγ-crystallins

βγ-Crystallins have emerged as a superfamily of structurally homologous proteins with representatives across the domains of life. A major portion of this superfamily is constituted by members from microorganisms. This superfamily has also been recognized as a novel group of Ca²⁺-binding proteins with huge diversity. The βγ domain shows variable properties in Ca²⁺ binding, stability and association with other domains. The various members present a series of evolutionary adaptations culminating in great diversity in properties and functions. Most of the predicted βγ-crystallins are yet to be characterized experimentally. In this review, we outline the distinctive features of microbial βγ-crystallins and their position in the βγ-crystallin superfamily.     

How important are transmembrane helices of bitopic membrane proteins?

The topology of a bitopic membrane protein consists of a single transmembrane helix connecting two extra-membranous domains. As opposed to helices from polytopic proteins, the transmembrane helices of bitopic proteins were initially considered as merely hydrophobic anchors, while more recent studies have begun to shed light on their role in the protein's function. Herein the overall importance of transmembrane helices from bitopic membrane proteins was analyzed using a relative conservation analysis. Interestingly, the transmembrane domains of bitopic proteins are on average, significantly more conserved than the remainder of the protein, even when taking into account their smaller amino acid repertoire. Analysis of highly conserved transmembrane domains did not reveal any unifying consensus, pointing to a great diversity in their conservation patterns. However, Fourier power spectrum analysis was able to show that regardless of the conservation motif, in most sequences a significant conservation moment was observed, in that one side of the helix was conserved while the other was not. Taken together, it may be possible to conclude that a significant proportion of transmembrane helices from bitopic membrane proteins participate in specific interactions, in a variety of modes in the plane of the lipid bilayer.     

How important are transmembrane helices of bitopic membrane proteins?

The topology of a bitopic membrane protein consists of a single transmembrane helix connecting two extra-membranous domains. As opposed to helices from polytopic proteins, the transmembrane helices of bitopic proteins were initially considered as merely hydrophobic anchors, while more recent studies have begun to shed light on their role in the protein's function. Herein the overall importance of transmembrane helices from bitopic membrane proteins was analyzed using a relative conservation analysis. Interestingly, the transmembrane domains of bitopic proteins are on average, significantly more conserved than the remainder of the protein, even when taking into account their smaller amino acid repertoire. Analysis of highly conserved transmembrane domains did not reveal any unifying consensus, pointing to a great diversity in their conservation patterns. However, Fourier power spectrum analysis was able to show that regardless of the conservation motif, in most sequences a significant conservation moment was observed, in that one side of the helix was conserved while the other was not. Taken together, it may be possible to conclude that a significant proportion of transmembrane helices from bitopic membrane proteins participate in specific interactions, in a variety of modes in the plane of the lipid bilayer.     

Novel fold revealed by the structure of a FAS1 domain pair from the insect cell adhesion molecule fasciclin I

Fasciclin I is an insect neural cell adhesion molecule consisting of four FAS1 domains, homologs of which are present in many bacterial, plant, and animal proteins. The crystal structure of FAS1 domains 3 and 4 of Drosophila fasciclin I reveals a novel domain fold, consisting of a seven-stranded beta wedge and a number of alpha helices. The two domains are arranged in a linear fashion and interact through a substantial polar interface. Missense mutations in the FAS1 domains of the human protein betaig-h3 cause corneal dystrophies. Many mutations alter highly conserved core residues, but the two most common mutations, affecting Arg-124 and Arg-555, map to exposed alpha-helical regions, suggesting reduced protein solubility as the disease mechanism.     

Solution structure of human Mts1 (S100A4) as determined by NMR spectroscopy

Mts1 is a member of the S100 family of Ca2+-binding proteins and is implicated in promoting tumor progression and metastasis. To better understand the structure-function relationships of this protein and to begin characterizing its Ca2+-dependent interaction with protein binding targets, the three-dimensional structure of mts1 was determined in the apo state by NMR spectroscopy. As with other S100 protein family members, mts1 is a symmetric homodimer held together by noncovalent interactions between two helices from each subunit (helices 1, 4, 1', and 4') to form an X-type four-helix bundle. Each subunit of mts1 has two EF-hand Ca2+-binding domains: a pseudo-EF-hand (or S100-hand) and a typical EF-hand that are brought into proximity by a small two-stranded antiparallel beta-sheet. The S100-hand is formed by helices 1 and 2, and is similar in conformation to other members of the S100 family. In the typical EF-hand, the position of helix 3 is similar to that of another member of the S100 protein family, calcyclin (S100A6), and less like that of other S100 family members for which three-dimensional structures are available in the calcium-free state (e.g., S100B and S100A1). The differences in the position of helix 3 in the apo state of these four S100 proteins are likely due to variations in the amino acid sequence in the C-terminus of helix 4 and in loop 2 (the hinge region) and could potentially be used to subclassify the S100 protein family.     

The DOMON domains are involved in heme and sugar recognition

We expand the functionally uncharacterized DOMON domain superfamily to identify several novel families, including the first prokaryotic representatives. Using several computational tools we show that it is involved in ligand binding--either as heme- or sugar-binding domains. We present evidence that the DOMON domain along with the DM13 domain comprises a novel electron-transfer system potentially involved in oxidative modification of animal cell-surface proteins. Other novel versions might function as sugar sensors of histidine kinases of bacterial two component systems. Supplementary information: Supplementary data are available at Bioinformatics online and also at ftp://ftp.ncbi.nih.gov/pub/aravind/domon/.     

Structural trees for protein superfamilies

Structural trees for large protein superfamilies, such as β proteins with the aligned β sheet packing, β proteins with the orthogonal packing of α helices, two-layer and three-layer α/β proteins, have been constructed. The structural motifs having unique overall folds and a unique handedness are taken as root structures of the trees. The larger protein structures of each superfamily are obtained by a stepwise addition of α helices and/or β strands to the corresponding root motif, taking into account a restricted set of rules inferred from known principles of the protein structure. Among these rules, prohibition of crossing connections, attention to handedness and compactness, and a requirement for α helices to be packed in α-helical layers and β strands in β layers are the most important. Proteins and domains whose structures can be obtained by stepwise addition of α helices and/or β strands to the same root motif can be grouped into one structural class or a superfamily. Proteins and domains found within branches of a structural tree can be grouped into subclasses or subfamilies. Levels of structural similarity between different proteins can easily be observed by visual inspection. Within one branch, protein structures having a higher position in the tree include the structures located lower. Proteins and domains of different branches have the structure located in the branching point as the common fold. Proteins 28:241–260, 1997. © 1997 Wiley-Liss Inc.