Comparative analysis of structural properties of the C-type-lectin-like domain (CTLD)

The superfamily of proteins containing the C-type-lectin-like domain (CTLD) is a group of abundant extracellular metazoan proteins characterized by evolutionary flexibility and functional versatility. Several CTLDs are also found in parasitic prokaryotes and viruses. The 37 distinct currently available CTLD structures demonstrate significant structural conservation despite low or undetectable sequence similarity. Our aim in this study was to perform an extensive comparative analysis of all available CTLD structures to establish the most conserved structural features of the fold, and to test and extend the early analysis of Drickamer. By implication, these features should be those critical for maintenance of integrity of the fold. By analyzing CTLD structures superimposed by several methods, we have established groups of conserved structural positions involved in fold maintenance but not in ligand binding; these are consistent with the fold's known functional flexibility. In addition to the well-recognized disulfide bridges, groups of conserved residues are involved in hydrophobic interactions stabilizing the core of the fold and the long loop region, and in an alpha2-beta1-beta5 polar interaction. Evaluation of the conclusions of the structure comparison study compared with alignments of all available human, mouse and Caenorhabditis elegans CTLD sequences showed that conservation patterns are preserved throughout the whole CTLD sequence space. Our observations provide an improved understanding of CTLD structure, and will help in identification of new CTLDs and the mechanisms that drive and constrain the coevolution of the structure and function of the fold.     

The C-type lectin-like receptors of Dectin-1 cluster in natural killer gene complex

Natural killer gene complex (NKC) encodes a group of proteins with a single C-type lectin-like domain, (CTLD) which can be subdivided several subfamilies according to their structures and expression patterns. The receptors containing the conserved calcium binding sites in the CTLD fold belong to group II of C-type lectin superfamily and are expressed on myeloid cells and non- myeloid cells. The receptors lacking conserved calcium binding sites in the CTLD fold have evolved to bind ligands other than carbohydrates independently on calcium and thereby are named as C-type lectin-like receptors. The C-type lectin-like receptors are previously thought to be exclusively expressed on natural killer (NK) cells and enable NK cells to discriminate self, missing self or altered self. However, some C-type lectin-like receptors are identified in myeloid cells and are intensely investigated, recently. These myeloid C-type lectin-like receptors, especially Dectin-1 cluster, have a wide variety of ligands, including those of exogenous origin, and play important roles in the physiological functions and pathological processes including immune homeostasis, immune defenses, and immune surveillance. In this review, we summarize each member of the Dectin-1 cluster, including their structural profiles, expression patterns, signaling properties as well as known physiological functions.     

Sweet Tooth, a novel receptor protein-tyrosine kinase with C-type lectin-like extracellular domains

A gene encoding a novel type of receptor protein-tyrosine kinase was identified in Hydra vulgaris. The extracellular portion of this receptor (which we have named Sweet Tooth) contains four C-type lectin-like domains (CTLDs). Comparison of the sequences of these domains with the sequences of the carbohydrate recognition domains of various vertebrate C-type lectins shows that Sweet Tooth CTLD1 and CTLD4 have amino acids in common with those shown to be involved in carbohydrate binding by the lectins. Comparison of sequences encoding CTLD1 from the Sweet Tooth genes from different species of Hydra shows variation in some of the conserved residues that participate in carbohydrate binding in C-type lectins. The Sweet Tooth gene is expressed widely in the Hydra polyp, and expression is particularly high in the endoderm of the tentacles. Treatment of polyps with peptides corresponding to sequences in the Sweet Tooth CTLDs results in the disintegration of the animal. These same peptides do not block adhesion or morphogenesis of Hydra cell aggregates.     

Charged residues at protein interaction interfaces: unexpected conservation and orchestrated divergence

Protein-protein interactions play an essential role in the functioning of cell. The importance of charged residues and their diverse role in protein-protein interactions have been well studied using experimental and computational methods. Often, charged residues located in protein interaction interfaces are conserved across the families of homologous proteins and protein complexes. However, on a large scale, it has been recently shown that charged residues are significantly less conserved than other residue types in protein interaction interfaces. The goal of this work is to understand the role of charged residues in the protein interaction interfaces through their conservation patterns. Here, we propose a simple approach where the structural conservation of the charged residue pairs is analyzed among the pairs of homologous binary complexes. Specifically, we determine a large set of homologous interactions using an interaction interface similarity measure and catalog the basic types of conservation patterns among the charged residue pairs. We find an unexpected conservation pattern, which we call the correlated reappearance, occurring among the pairs of homologous interfaces more frequently than the fully conserved pairs of charged residues. Furthermore, the analysis of the conservation patterns across different superkingdoms as well as structural classes of proteins has revealed that the correlated reappearance of charged residues is by far the most prevalent conservation pattern, often occurring more frequently than the unconserved charged residues. We discuss a possible role that the new conservation pattern may play in the long-range electrostatic steering effect.     

The C-type lectin-like domain superfamily

The superfamily of proteins containing C-type lectin-like domains (CTLDs) is a large group of extracellular Metazoan proteins with diverse functions. The CTLD structure has a characteristic double-loop ('loop-in-a-loop') stabilized by two highly conserved disulfide bridges located at the bases of the loops, as well as a set of conserved hydrophobic and polar interactions. The second loop, called the long loop region, is structurally and evolutionarily flexible, and is involved in Ca2+-dependent carbohydrate binding and interaction with other ligands. This loop is completely absent in a subset of CTLDs, which we refer to as compact CTLDs; these include the Link/PTR domain and bacterial CTLDs. CTLD-containing proteins (CTLDcps) were originally classified into seven groups based on their overall domain structure. Analyses of the superfamily representation in several completely sequenced genomes have added 10 new groups to the classification, and shown that it is applicable only to vertebrate CTLDcps; despite the abundance of CTLDcps in the invertebrate genomes studied, the domain architectures of these proteins do not match those of the vertebrate groups. Ca2+-dependent carbohydrate binding is the most common CTLD function in vertebrates, and apparently the ancestral one, as suggested by the many humoral defense CTLDcps characterized in insects and other invertebrates. However, many CTLDs have evolved to specifically recognize protein, lipid and inorganic ligands, including the vertebrate clade-specific snake venoms, and fish antifreeze and bird egg-shell proteins. Recent studies highlight the functional versatility of this protein superfamily and the CTLD scaffold, and suggest further interesting discoveries have yet to be made.     

Combining protein sequence,structure, and dynamics: A novel approach for functional evolution analysis of PAS domain superfamily

PAS domains are widespread in archaea, bacteria, and eukaryota, and play important roles in various functions. In this study, we aim to explore functional evolutionary relationship among proteins in the PAS domain superfamily in view of the sequence‐structure‐dynamics‐function relationship. We collected protein sequences and crystal structure data from RCSB Protein Data Bank of the PAS domain superfamily belonging to three biological functions (nucleotide binding, photoreceptor activity, and transferase activity). Protein sequences were aligned and then used to select sequence‐conserved residues and build phylogenetic tree. Three‐dimensional structure alignment was also applied to obtain structure‐conserved residues. The protein dynamics were analyzed using elastic network model (ENM) and validated by molecular dynamics (MD) simulation. The result showed that the proteins with same function could be grouped by sequence similarity, and proteins in different functional groups displayed statistically significant difference in their vibrational patterns. Interestingly, in all three functional groups, conserved amino acid residues identified by sequence and structure conservation analysis generally have a lower fluctuation than other residues. In addition, the fluctuation of conserved residues in each biological function group was strongly correlated with the corresponding biological function. This research suggested a direct connection in which the protein sequences were related to various functions through structural dynamics. This is a new attempt to delineate functional evolution of proteins using the integrated information of sequence, structure, and dynamics.     

ConSurf: identification of functional regions in proteins by surface-mapping of phylogenetic information

SUMMARY: We recently developed algorithmic tools for the identification of functionally important regions in proteins of known three dimensional structure by estimating the degree of conservation of the amino-acid sites among their close sequence homologues. Projecting the conservation grades onto the molecular surface of these proteins reveals patches of highly conserved (or occasionally highly variable) residues that are often of important biological function. We present a new web server, ConSurf, which automates these algorithmic tools. ConSurf may be used for high-throughput characterization of functional regions in proteins. AVAILABILITY: The ConSurf web server is available at:http://consurf.tau.ac.il. SUPPLEMENTARY INFORMATION: A set of examples is available at http://consurf.tau.ac.il under 'GALLERY'.     

Conserved core structure and active site residues in alkaline phosphatase superfamily enzymes.

Cofactor-independent phosphoglycerate mutase (iPGM) has been previously identified as a member of the alkaline phosphatase (AlkP) superfamily of enzymes, based on the conservation of the predicted metal-binding residues. Structural alignment of iPGM with AlkP and cerebroside sulfatase confirmed that all these enzymes have a common core structure and revealed similarly located conserved Ser (in iPGM and AlkP) or Cys (in sulfatases) residues in their active sites. In AlkP, this Ser residue is phosphorylated during catalysis, whereas in sulfatases the active site Cys residues are modified to formylglycine and sulfatated. Similarly located Thr residue forms a phosphoenzyme intermediate in one more enzyme of the AlkP superfamily, alkaline phosphodiesterase/nucleotide pyrophosphatase PC-1 (autotaxin). Using structure-based sequence alignment, we identified homologous Ser, Thr, or Cys residues in other enzymes of the AlkP superfamily, such as phosphopentomutase, phosphoglycerol transferase, phosphonoacetate hydrolase, and GPI-anchoring enzymes (glycosylphosphatidylinositol phosphoethanolamine transferases) MCD4, GPI7, and GPI13. We predict that catalytical cycles of all the enzymes of AlkP superfamily include phosphoenzyme (or sulfoenzyme) intermediates.     

AiCTL-6, a novel C-type lectin from bay scallop Argopecten irradians with a long C-type lectin-like domain

C-type lectins are a superfamily of carbohydrate-recognition proteins which play crucial roles in the innate immunity. In this study, a novel C-type lectin gene from scallop Argopecten irradians (designated as AiCTL-6) was cloned by rapid amplification of cDNA ends (RACE) approach based on expression sequence tag (EST) analysis. The full-length cDNA of AiCTL-6 was 1080 bp. The open reading frame encoded a polypeptide of 307 amino acids, including a signal sequence and a C-type lectin-like domain (CTLD) of 150 amino acid residues longer than any usual CTLD. It contained six conserved cysteine residues involved in the formation of three internal disulfide bridges and an EPD (Glu₂₆₉-Pro₂₇₀-Asp₂₇₁) motif at the Ca²⁺-binding site 2. The deduced amino acid sequence of AiCTL-6 showed high similarity to members of C-type lectin superfamily. By fluorescent quantitative real-time PCR, AiCTL-6 mRNA was found mainly in hepatopancreas and gill, and marginally expressed in other tissues. After the scallops were challenged by Listonella anguillarum for 6 h, the mRNA expression of AiCTL-6 was up-regulated significantly to 7.2-fold compared to the blank group. While at 9 h post Micrococcus luteus challenge, its expression level was 60.1 times higher than that of the blank group. The functional activity of AiCTL-6 was investigated by recombination and expression of the cDNA fragment encoding its mature peptide in Escherichia coli Rosetta gami (DE3). The recombinant AiCTL-6 could agglutinate Gram-negative bacteria E. coli TOP10F′, Gram-positive bacteria M. luteus and Staphylococcus aureus. These results collectively suggested that AiCTL-6, as a novel member of C-type lectin family, contributed to the host defense mechanisms against invading microorganism in A. irradians.     

Crystal structure of the C-type lectin-like domain from the human hematopoietic cell receptor CD69

CD69, one of the earliest specific antigens acquired during lymphoid activation, acts as a signal-transducing receptor involved in cellular activation events, including proliferation and induction of specific genes. CD69 belongs to a family of receptors that modulate the immune response and whose genes are clustered in the natural killer (NK) gene complex. The extracellular portion of these receptors represent a subfamily of C-type lectin-like domains (CTLDs), which are divergent from true C-type lectins and are referred to as NK-cell domains (NKDs). We have determined the three-dimensional structure of human CD69 NKD in two different crystal forms. CD69 NKD adopts the canonical CTLD fold but lacks the features involved in Ca(2+) and carbohydrate binding by C-type lectins. CD69 NKD dimerizes noncovalently, both in solution and in crystalline state. The dimer interface consists of a hydrophobic, loosely packed core, surrounded by polar interactions, including an interdomain beta sheet. The intersubunit core shows certain structural plasticity that may facilitate conformational rearrangements for binding to ligands. The surface equivalent to the binding site of other members of the CTLD superfamily reveals a hydrophobic patch surrounded by conserved charged residues that probably constitutes the CD69 ligand-binding site.     

Comparative sequence analysis of acid sensitive/resistance proteins in Escherichia coli and Shigella flexneri

The molecular basis for the survival of bacteria under extreme conditions in which growth is inhibited is a question of great current interest. A preliminary study was carried out to determine residue pattern conservation among the antiporters of enteric bacteria, responsible for extreme acid sensitivity especially in Escherichia coli and Shigella flexneri. Here we found the molecular evidence that proved the relationship between E. coli and S. flexneri. Multiple sequence alignment of the gadC coded acid sensitive antiporter showed many conserved residue patterns at regular intervals at the N-terminal region. It was observed that as the alignment approaches towards the C-terminal, the number of conserved residues decreases, indicating that the N-terminal region of this protein has much active role when compared to the carboxyl terminal. The motif, FHLVFFLLLGG, is well conserved within the entire gadC coded protein at the amino terminal. The motif is also partially conserved among other antiporters (which are not coded by gadC) but involved in acid sensitive/resistance mechanism. Phylogenetic cluster analysis proves the relationship of Escherichia coli and Shigella flexneri. The gadC coded proteins are converged as a clade and diverged from other antiporters belongs to the amino acid-polyamine-organocation (APC) superfamily.     

Ramachandran analysis of conserved glycyl residues in homologous proteins of known structure

High conservation of glycyl residues in homologous proteins is fairly frequent. It is commonly understood that glycine tends to be highly conserved either because of its unique Ramachandran angles or to avoid steric clash that would arise with a larger side chain. Using a database of aligned 3D structures of homologous proteins we identified conserved Gly in 288 alignment positions from 85 families. Ninety‐six of these alignment positions correspond to conserved Gly residue with (φ, ψ) values allowed for non‐glycyl residues. Reasons for this observation were investigated by in‐silico mutation of these glycyl residues to Ala. We found in 94% of the cases a short contact exists between the C^β atom of the introduced Ala with the atoms which are often distant in the primary structure. This suggests the lack of space even for a short side chain thereby explaining high conservation of glycyl residues even when they adopt (φ, ψ) values allowed for Ala. In 189 alignment positions, the conserved glycyl residues adopt (φ, ψ) values which are disallowed for Ala. In‐silico mutation of these Gly residues to Ala almost always results in steric hindrance involving C^β atom of Ala as one would expect by comparing Ramachandran maps for Ala and Gly. Rare occurrence of the disallowed glycyl conformations even in ultrahigh resolution protein structures are accompanied by short contacts in the crystal structures and such disallowed conformations are not conserved in the homologues. These observations raise the doubt on the accuracy of such glycyl conformations in proteins.     

Evolutionary conservation of protein vibrational dynamics

The aim of the present work is to study the evolutionary divergence of vibrational protein dynamics. To this end, we used the Gaussian Network Model to perform a systematic analysis of normal mode conservation on a large dataset of proteins classified into homologous sets of family pairs and superfamily pairs. We found that the lowest most collective normal modes are the most conserved ones. More precisely, there is, on average, a linear correlation between normal mode conservation and mode collectivity. These results imply that the previously observed conservation of backbone flexibility (B-factor) profiles is due to the conservation of the most collective modes, which contribute the most to such profiles. We discuss the possible roles of normal mode robustness and natural selection in the determination of the observed behavior. Finally, we draw some practical implications for dynamics-based protein alignment and classification and discuss possible caveats of the present approach.     

Sequence and comparative genomic analysis of actin-related proteins

Actin-related proteins (ARPs) are key players in cytoskeleton activities and nuclear functions. Two complexes, ARP2/3 and ARP1/11, also known as dynactin, are implicated in actin dynamics and in microtubule-based trafficking, respectively. ARP4 to ARP9 are components of many chromatin-modulating complexes. Conventional actins and ARPs codefine a large family of homologous proteins, the actin superfamily, with a tertiary structure known as the actin fold. Because ARPs and actin share high sequence conservation, clear family definition requires distinct features to easily and systematically identify each subfamily. In this study we performed an in depth sequence and comparative genomic analysis of ARP subfamilies. A high-quality multiple alignment of approximately 700 complete protein sequences homologous to actin, including 148 ARP sequences, allowed us to extend the ARP classification to new organisms. Sequence alignments revealed conserved residues, motifs, and inserted sequence signatures to define each ARP subfamily. These discriminative characteristics allowed us to develop ARPAnno (http://bips.u-strasbg.fr/ARPAnno), a new web server dedicated to the annotation of ARP sequences. Analyses of sequence conservation among actins and ARPs highlight part of the actin fold and suggest interactions between ARPs and actin-binding proteins. Finally, analysis of ARP distribution across eukaryotic phyla emphasizes the central importance of nuclear ARPs, particularly the multifunctional ARP4.     

A consensus prediction of the secondary structure for the 6-phospho-β-D-galactosidase superfamily

Two separate unrefined models for the secondary structure of two subfamilies of the 6-phospho-β-D -galactosidase superfamily were independently constructed by examining patterns of variation and conservation within homologous protein sequences, assigning surface, interior, parsing, and active site residues to positions in the alignment, and identifying periodicities in these. A consensus model for the secondary structure of the entire superfamily was then built. The prediction tests the limits of an unrefined prediction made using this approach in a large protein with substantial functional and sequence divergence within the family. The protein belongs to the (α–β class), with the core β strands aligned parallel. The supersecondary structural elements that are readily identified in this model is a parallel β sheet built by strands C, D, and E, with helices 2 and 3 connecting strands (C + D) and (D + E), respectively, and an analogous α–β unit (strand G and helix 7) toward the end of the sequence. The resemblance of the supersecondary model to the tertiary structure formed by 8-fold α–β barrel proteins is almost certainly not coincidental. © 1995 Wiley-Liss, Inc.     

Molecular architecture of mouse activating NKR-P1 receptors

Receptors belonging to NKR-P1 family and their specific Clr ligands form an alternative missing self recognition system critical in immunity against tumors and viruses, elimination of tumor cells subjected to genotoxic stress, activation of T cell dependent immune response, and hypertension. The three-dimensional structure of the extracellular domain of the mouse natural killer (NK) cell receptor mNKR-P1Aex has been determined by X-ray diffraction. The core of the C-type lectin domain (CTLD) is homologous to the other CTLD receptors whereas one quarter of the domain forms an extended loop interacting tightly with a neighboring loop in the crystal. This domain swapping mechanism results in a compact interaction interface. A second dimerization interface resembles the known arrangement of other CTLD NK receptors. A functional dimeric form of the receptor is suggested, with the loop, evolutionarily conserved within this family, proposed to participate in interactions with ligands.     

Molecular cloning and sequence analysis of aggretin, a collagen-like platelet aggregation inducer.

A cDNA library derived from the Malayan-pit-viper (Calloselasma rhodostoma) venom gland was constructed in the phagemid vector. Using the information of the N-terminal amino acid sequences of two subunits of aggretin, synthetic mixed-base oligonucleotides were employed as a screening probe for colony hybridization. Separate cDNA clones encoding for the alpha and beta chains of aggretin were isolated and sequenced. The results revealed that mature alpha and beta chains contain 136 and 123 amino acid residues, respectively. Aggretin subunits show high degrees of identity with respective subunits (50-60% for alpha, 49-58% for beta) of C-type lectin-like snake venoms. The identity to rattlesnake lectin is relatively lower (i.e., 39 and 30%). All cysteine residues in each chain of aggretin are well conserved and located at the positions corresponding to those of C-type lectins. Thus, three intracatenary disulfide bridges and an interchain disulfide bond between Cys83(alpha) and Cys75(beta) may be allocated. This is the first report regarding the entire sequence of venom GPIa/IIa agonist. According to the alignment of amino acid sequences, hypervariable regions among these C-type lectin-like proteins were revealed. These hypervariable regions are proposed to be the counterparts directly interacting with different receptors or different domains of a receptor on the surface of platelet.     

Structural Variations in Protein Superfamilies: Actin and Tubulin

Structures of homologous proteins are usually conserved during evolution, as are critical active site residues. This is the case for actin and tubulin, the two most important cytoskeleton proteins in eukaryotes. Actins and their related proteins (Arps) constitute a large superfamily whereas the tubulin family has fewer members. Unaligned sequences of these two protein families were analysed by searching for short groups of family-specific amino acid residues, that we call motifs, and by counting the number of residues from one motif to the next. For each sequence, the set of motif-to-motif residue counts forms a subfamily-specific pattern (landmark pattern) allowing actin and tubulin superfamily members to be identified and sorted into subfamilies. The differences between patterns of individual subfamilies are due to inserts and deletions (indels). Inserts appear to have arisen at an early stage in eukaryote evolution as suggested by the small but consistent kingdom-dependent differences found within many Arp subfamilies and in γ-tubulins. Inserts tend to be in surface loops where they can influence subfamily-specific function without disturbing the core structure of the protein. The relatively few indels found for tubulins have similar positions to established results, whereas we find many previously unreported indel positions and lengths for the metazoan Arps.     

The immunoglobulin fold family: sequence analysis and 3D structure comparisons.

Fifty-two 3D structures of Ig-like domains covering the immunoglobulin fold family (IgFF) were compared and classified according to the conservation of their secondary structures. Members of the IgFF are distantly related proteins or evolutionarily unrelated proteins with a similar fold, the Ig fold. In this paper, a multiple structural alignment of the conserved common core is described and the correlation between corresponding sequences is discussed. While the members of the IgFF exhibit wide heterogeneity in terms of tissue and species distribution or functional implications, the 3D structures of these domains are far more conserved than their sequences. We define topologically equivalent residues in the Ig-like domains, describe the hydrophobic common cores and discuss the presence of additional strands. The disulfide bridges, not necessary for the stability of the Ig fold, may have an effect on the compactness of the domains. Based upon sequence and structure analysis, we propose the introduction of two new subtypes (C3 and C4) to the previous classifications, in addition to a new global structural classification. The very low mean sequence identity between subgroups of the IgFF suggests the occurrence of both divergent and convergent evolutionary processes, explaining the wide diversity of the superfamily. Finally, this review suggest that hydrophobic residues constituting the common hydrophobic cores are important clues to explain how highly divergent sequences can adopt a similar fold.     

Development and validation of a consistency based multiple structure alignment algorithm

SUMMARY: We introduce an algorithm that uses the information gained from simultaneous consideration of an entire group of related proteins to create multiple structure alignments (MSTAs). Consistency-based alignment (CBA) first harnesses the information contained within regions that are consistently aligned among a set of pairwise superpositions in order to realign pairs of proteins through both global and local refinement methods. It then constructs a multiple alignment that is maximally consistent with the improved pairwise alignments. We validate CBA's alignments by assessing their accuracy in regions where at least two of the aligned structures contain the same conserved sequence motif. RESULTS: CBA correctly aligns well over 90% of motif residues in superpositions of proteins belonging to the same family or superfamily, and it outperforms a number of previously reported MSTA algorithms.