The Evolution of Guanylyl Cyclases as Multidomain Proteins: Conserved Features of Kinase-Cyclase Domain Fusions

Guanylyl cyclases (GCs) are enzymes that generate cyclic GMP and regulate different physiologic and developmental processes in a number of organisms. GCs possess sequence similarity to class III adenylyl cyclases (ACs) and are present as either membrane-bound receptor GCs or cytosolic soluble GCs. We sought to determine the evolution of GCs using a large-scale bioinformatic analysis and found multiple lineage-specific expansions of GC genes in the genomes of many eukaryotes. Moreover, a few GC-like proteins were identified in prokaryotes, which come fused to a number of different domains, suggesting allosteric regulation of nucleotide cyclase activity Eukaryotic receptor GCs are associated with a kinase homology domain (KHD), and phylogenetic analysis of these proteins suggest coevolution of the KHD and the associated cyclase domain as well as a conservation of the sequence and the size of the linker region between the KHD and the associated cyclase domain. Finally, we also report the existence of mimiviral proteins that contain putative active kinase domains associated with a cyclase domain, which could suggest early evolution of the fusion of these two important domains involved in signa transduction.

The Linker Region in Receptor Guanylyl Cyclases Is a Key Regulatory Module MUTATIONAL ANALYSIS OF GUANYLYL CYCLASE C

Receptor guanylyl cyclases are multidomain proteins, and ligand binding to the extracellular domain increases the levels of intracellular cGMP. The intracellular domain of these receptors is composed of a kinase homology domain (KHD), a linker of similar to 70 amino acids, followed by the C-terminal guanylyl cyclase domain. Mechanisms by which these receptors are allosterically regulated by ligand binding to the extracellular domain and ATP binding to the KHD are not completely understood. Here we examine the role of the linker region in receptor guanylyl cyclases by a series of point mutations in receptor guanylyl cyclase C. The linker region is predicted to adopt a coiled coil structure and aid in dimerization, but we find that the effects of mutations neither follow a pattern predicted for a coiled coil peptide nor abrogate dimerization. Importantly, this region is critical for repressing the guanylyl cyclase activity of the receptor in the absence of ligand and permitting ligand-mediated activation of the cyclase domain. Mutant receptors with high basal guanylyl cyclase activity show no further activation in the presence of non-ionic detergents, suggesting that hydrophobic interactions in the basal and inactive conformation of the guanylyl cyclase domain are disrupted by mutation. Equivalent mutations in the linker region of guanylyl cyclase A also elevated the basal activity and abolished ligand-and detergent-mediated activation. We, therefore, have defined a key regulatory role for the linker region of receptor guanylyl cyclases which serves as a transducer of information from the extracellular domain via the KHD to the catalytic domain.

Src Homology 3-interacting Domain of Rv1917c of Mycobacterium tuberculosis Induces Selective Maturation of Human Dendritic Cells by Regulating PI3K-MAPK-NF-kappa B Signaling and Drives Th2 Immune Responses

Mycobacterium tuberculosis, an etiological agent of pulmonary tuberculosis, causes significant morbidity and mortality worldwide. Pathogenic mycobacteria survive in the host by subverting host innate immunity. Dendritic cells (DCs) are professional antigen-presenting cells that are vital for eliciting immune responses to infectious agents, including pathogenic mycobacteria. DCs orchestrate distinct Th responses based on the signals they receive. In this perspective, deciphering the interactions of the proline-glutamic acid/proline-proline-glutamic acid (PE/PPE) family of proteins of M. tuberculosis with DCs assumes significant pathophysiological attributes. In this study, we demonstrate that Rv1917c (PPE34), a representative member of the proline-proline-glutamic-major polymorphic tandem repeat family, interacts with TLR2 and triggers functional maturation of human DCs. Signaling perturbations implicated a critical role for integrated cross-talk among PI3K-MAPK and NF-kappa B signaling cascades in Rv1917c-induced maturation of DCs. However, this maturation of DCs was associated with a secretion of high amounts of anti-inflammatory cytokine IL-10, whereas Th1-polarizing cytokine IL-12 was not induced. Consistent with these results, Rv1917c-matured DCs favored secretion of IL-4, IL-5, and IL-10 from CD4(+) T cells and contributed to Th2-skewed cytokine balance ex vivo in healthy individuals and in patients with pulmonary tuberculosis. Interestingly, the Rv1917c-skewed Th2 immune response involved induced expression of cyclooxygenase-2 (COX-2) in DCs. Taken together, these results indicate that Rv1917c facilitates a shift in the ensuing immunity toward the Th2 phenotype and could aid in immune evasion by mycobacteria.

Genetic studies of the PRP17 gene of Saccharomyces cerevisiae: a domain essential for function maps to a nonconserved region of the protein

The PRP17 gene product is required for the second step of pre-mRNA splicing reactions. The C-terminal half of this protein bears four repeat units with homology to the beta transducin repeat. Missense mutations in three temperature-sensitive prp17 mutants map to a region in the N-terminal half of the protein. We have generated, in vitro, 11 missense alleles at the beta transducin repeat units and find that only one affects function in vivo. A phenotypically silent missense allele at the fourth repeat unit enhances the slow-growing phenotype conferred by an allele at the third repeat, suggesting an interaction between these domains. Although many missense mutations in highly conserved amino acids lack phenotypic effects, deletion analysis suggests an essential role for these units. Only mutations in the N-terminal nonconserved domain of PRP17 are synthetically lethal in combination with mutations in PRP16 and PRP18, two other gene products required for the second splicing reaction. A mutually allele-specific interaction between Prp17 and snr7, with mutations in U5 snRNA, was observed. We therefore suggest that the functional region of Prp17p that interacts with Prp18p, Prp16p, and U5 snRNA is the N terminal region of the protein.

Biochemical characterization of the intracellular domain of the human guanylyl cyclase C receptor provides evidence for a catalytically active homotrimer

Guanylyl cyclase C (GCC) is the receptor for the family of guanylin peptides and bacterial heat-stable enterotoxins (ST). The receptor is composed of an extracellular, ligand-binding domain and an intracellular domain with a region of homology to protein kinases and a guanylyl cyclase catalytic domain. We have expressed the entire intracellular domain of GCC in insect cells and purified the recombinant protein, GCC-IDbac, to study its catalytic activity and regulation. Kinetic properties of the purified protein were similar to that of full-length GCC, and high activity was observed when MnGTP was used as the substrate. Nonionic detergents, which stimulate the guanylyl cyclase activity of membrane-associated GCC, did not appreciably increase the activity of GCC-IDbac, indicating that activation of the receptor by Lubrol involved conformational changes that required the transmembrane and/or the extracellular domain. The guanylyl cyclase activity of GCC-IDbac was inhibited by Zn2+, at concentrations shown to inhibit adenylyl cyclase, suggesting a structural homology between the two enzymes. Covalent crosslinking of GCC-IDbac indicated that the protein could associate as a dimer, but a large fraction was present as a trimer. Gel filtration analysis also showed that the major fraction of the protein eluted at a molecular size of a trimer, suggesting that the dimer detected by cross-linking represented subtle differences in the juxtaposition of the individual polypeptide chains. We therefore provide evidence that the trimeric state of GCC is catalytically active, and sequences required to generate the trimer are present in the intracellular domain of GCC.

The repertoire of protein kinases encoded in the draft version of the human genome: atypical variations and uncommon domain combinations

Background: Phosphorylation by protein kinases is central to cellular signal transduction. Abnormal functioning of kinases has been implicated in developmental disorders and malignancies. Their activity is regulated by second messengers and by the binding of associated domains, which are also influential in translocating the catalytic component to their substrate sites, in mediating interaction with other proteins and carrying out their biological roles. Results: Using sensitive profile-search methods and manual analysis, the human genome has been surveyed for protein kinases. A set of 448 sequences, which show significant similarity to protein kinases and contain the critical residues essential for kinase function, have been selected for an analysis of domain combinations after classifying the kinase domains into subfamilies. The unusual domain combinations in particular kinases suggest their involvement in ubiquitination pathways and alternative modes of regulation for mitogen-activated protein kinase kinases (MAPKKs) and cyclin-dependent kinase (CDK)-like kinases. Previously unexplored kinases have been implicated in osteoblast differentiation and embryonic development on the basis of homology with kinases of known functions from other organisms. Kinases potentially unique to vertebrates are involved in highly evolved processes such as apoptosis, protein translation and tyrosine kinase signaling. In addition to coevolution with the kinase domain, duplication and recruitment of non-catalytic domains is apparent in signaling domains such as the PH, DAG-PE, SH2 and SH3 domains. Conclusions: Expansion of the functional repertoire and possible existence of alternative modes of regulation of certain kinases is suggested by their uncommon domain combinations. Experimental verification of the predicted implications of these kinases could enhance our understanding of their biological roles.

Filling-in Void and Sparse Regions in Protein Sequence Space by Protein-Like Artificial Sequences Enables Remarkable Enhancement in Remote Homology Detection Capability

Protein functional annotation relies on the identification of accurate relationships, sequence divergence being a key factor. This is especially evident when distant protein relationships are demonstrated only with three-dimensional structures. To address this challenge, we describe a computational approach to purposefully bridge gaps between related protein families through directed design of protein-like ``linker'' sequences. For this, we represented SCOP domain families, integrated with sequence homologues, as multiple profiles and performed HMM-HMM alignments between related domain families. Where convincing alignments were achieved, we applied a roulette wheel-based method to design 3,611,010 protein-like sequences corresponding to 374 SCOP folds. To analyze their ability to link proteins in homology searches, we used 3024 queries to search two databases, one containing only natural sequences and another one additionally containing designed sequences. Our results showed that augmented database searches showed up to 30% improvement in fold coverage for over 74% of the folds, with 52 folds achieving all theoretically possible connections. Although sequences could not be designed between some families, the availability of designed sequences between other families within the fold established the sequence continuum to demonstrate 373 difficult relationships. Ultimately, as a practical and realistic extension, we demonstrate that such protein-like sequences can be ``plugged-into'' routine and generic sequence database searches to empower not only remote homology detection but also fold recognition. Our richly statistically supported findings show that complementary searches in both databases will increase the effectiveness of sequence-based searches in recognizing all homologues sharing a common fold. (C) 2013 Elsevier Ltd. All rights reserved.

Enriching the annotation of Mycobacterium tuberculosis H37Rv proteome using remote homology detection approaches: Insights into structure and function

The availability of the genome sequence of Mycobacterium tuberculosis H37Rv has encouraged determination of large numbers of protein structures and detailed definition of the biological information encoded therein; yet, the functions of many proteins in M. tuberculosis remain unknown. The emergence of multidrug resistant strains makes it a priority to exploit recent advances in homology recognition and structure prediction to re-analyse its gene products. Here we report the structural and functional characterization of gene products encoded in the M. tuberculosis genome, with the help of sensitive profile-based remote homology search and fold recognition algorithms resulting in an enhanced annotation of the proteome where 95% of the M. tuberculosis proteins were identified wholly or partly with information on structure or function. New information includes association of 244 proteins with 205 domain families and a separate set of new association of folds to 64 proteins. Extending structural information across uncharacterized protein families represented in the M. tuberculosis proteome, by determining superfamily relationships between families of known and unknown structures, has contributed to an enhancement in the knowledge of structural content. In retrospect, such superfamily relationships have facilitated recognition of probable structure and/or function for several uncharacterized protein families, eventually aiding recognition of probable functions for homologous proteins corresponding to such families. Gene products unique to mycobacteria for which no functions could be identified are 183. Of these 18 were determined to be M. tuberculosis specific. Such pathogen-specific proteins are speculated to harbour virulence factors required for pathogenesis. A re-annotated proteome of M. tuberculosis, with greater completeness of annotated proteins and domain assigned regions, provides a valuable basis for experimental endeavours designed to obtain a better understanding of pathogenesis and to accelerate the process of drug target discovery. (C) 2014 Elsevier Ltd. All rights reserved.

NrichD database: sequence databases enriched with computationally designed protein-like sequences aid in remote homology detection

NrichD ( ext-link-type=''uri'' xlink:href=''http://proline.biochem.iisc.ernet.in/NRICHD/'' xlink:type=''simple''>http://proline.biochem.iisc.ernet.in/NRICHD/)< /named-content> is a database of computationally designed protein-like sequences, augmented into natural sequence databases that can perform hops in protein sequence space to assist in the detection of remote relationships. Establishing protein relationships in the absence of structural evidence or natural `intermediately related sequences' is a challenging task. Recently, we have demonstrated that the computational design of artificial intermediary sequences/linkers is an effective approach to fill naturally occurring voids in protein sequence space. Through a large-scale assessment we have demonstrated that such sequences can be plugged into commonly employed search databases to improve the performance of routinely used sequence search methods in detecting remote relationships. Since it is anticipated that such data sets will be employed to establish protein relationships, two databases that have already captured these relationships at the structural and functional domain level, namely, the SCOP database and the Pfam database, have been `enriched' with these artificial intermediary sequences. NrichD database currently contains 3 611 010 artificial sequences that have been generated between 27 882 pairs of families from 374 SCOP folds. The data sets are freely available for download. Additional features include the design of artificial sequences between any two protein families of interest to the user.

De-DUFing the DUFs: Deciphering distant evolutionary relationships of Domains of Unknown Function using sensitive homology detection methods

Background: In the post-genomic era where sequences are being determined at a rapid rate, we are highly reliant on computational methods for their tentative biochemical characterization. The Pfam database currently contains 3,786 families corresponding to ``Domains of Unknown Function'' (DUF) or ``Uncharacterized Protein Family'' (UPF), of which 3,087 families have no reported three-dimensional structure, constituting almost one-fourth of the known protein families in search for both structure and function. Results: We applied a `computational structural genomics' approach using five state-of-the-art remote similarity detection methods to detect the relationship between uncharacterized DUFs and domain families of known structures. The association with a structural domain family could serve as a start point in elucidating the function of a DUF. Amongst these five methods, searches in SCOP-NrichD database have been applied for the first time. Predictions were classified into high, medium and low-confidence based on the consensus of results from various approaches and also annotated with enzyme and Gene ontology terms. 614 uncharacterized DUFs could be associated with a known structural domain, of which high confidence predictions, involving at least four methods, were made for 54 families. These structure-function relationships for the 614 DUF families can be accessed on-line at http://proline.biochem.iisc.ernet.in/RHD_DUFS/. For potential enzymes in this set, we assessed their compatibility with the associated fold and performed detailed structural and functional annotation by examining alignments and extent of conservation of functional residues. Detailed discussion is provided for interesting assignments for DUF3050, DUF1636, DUF1572, DUF2092 and DUF659. Conclusions: This study provides insights into the structure and potential function for nearly 20 % of the DUFs. Use of different computational approaches enables us to reliably recognize distant relationships, especially when they converge to a common assignment because the methods are often complementary. We observe that while pointers to the structural domain can offer the right clues to the function of a protein, recognition of its precise functional role is still `non-trivial' with many DUF domains conserving only some of the critical residues. It is not clear whether these are functional vestiges or instances involving alternate substrates and interacting partners. Reviewers: This article was reviewed by Drs Eugene Koonin, Frank Eisenhaber and Srikrishna Subramanian.

Homology between jacalin and artocarpin from jackfruit (Artocarpus integrifolia) seeds. Partial sequence and preliminary crystallographic studies of artocarpin

Jacalin and artocarpin, the two lectins from jackfruit (Artocarpus integrifolia) seeds, have different physicochemical properties and carbohydrate-binding specificities. However, comparison of the partial amino-acid sequence of artocarpin with the known sequence of jacalin indicates close to 50% sequence identity. Artocarpin crystallizes in two forms, both monoclinic P2(1), with one and two tetramic molecules, respectively, in the asymmetric units of form I (a = 69.9, b = 73.7, c = 60.6 Angstrom and beta = 95.1 degrees) and form II (a = 87.6, b = 72.2, c = 92.6 Angstrom and beta = 101.1 degrees). Both the crystal structures have been solved by the molecular replacement method using the known structure of jacalin as the search model and ope of them partially refined, confirming that the two lectins are indeed homologous.

Analysis on conservation of disulphide bonds and their structural features in homologous protein domain families

Background: Disulphide bridges are well known to play key roles in stability, folding and functions of proteins. Introduction or deletion of disulphides by site-directed mutagenesis have produced varying effects on stability and folding depending upon the protein and location of disulphide in the 3-D structure. Given the lack of complete understanding it is worthwhile to learn from an analysis of extent of conservation of disulphides in homologous proteins. We have also addressed the question of what structural interactions replaces a disulphide in a homologue in another homologue. Results: Using a dataset involving 34,752 pairwise comparisons of homologous protein domains corresponding to 300 protein domain families of known 3-D structures, we provide a comprehensive analysis of extent of conservation of disulphide bridges and their structural features. We report that only 54% of all the disulphide bonds compared between the homologous pairs are conserved, even if, a small fraction of the non-conserved disulphides do include cytoplasmic proteins. Also, only about one fourth of the distinct disulphides are conserved in all the members in protein families. We note that while conservation of disulphide is common in many families, disulphide bond mutations are quite prevalent. Interestingly, we note that there is no clear relationship between sequence identity between two homologous proteins and disulphide bond conservation. Our analysis on structural features at the sites where cysteines forming disulphide in one homologue are replaced by non-Cys residues show that the elimination of a disulphide in a homologue need not always result in stabilizing interactions between equivalent residues. Conclusion: We observe that in the homologous proteins, disulphide bonds are conserved only to a modest extent. Very interestingly, we note that extent of conservation of disulphide in homologous proteins is unrelated to the overall sequence identity between homologues. The non-conserved disulphides are often associated with variable structural features that were recruited to be associated with differentiation or specialisation of protein function.

Topological mimicry and epitope duplication in the guanylyl cyclase C receptor

Guanylyl cyclase C (GCC) is the receptor for the gastrointestinal hormones, guanylin, and uroguanylin, in addition to the bacterial heat-stable enterotoxins, which are one of the major causes of watery diarrhea the world over. GCC is expressed in intestinal cells, colorectal tumor tissue and tumors originating from metastasis of the colorectal carcinoma. We have earlier generated a monoclonal antibody to human GCC, GCC:B10, which was useful for the immunohistochemical localization of the receptor in the rat intestine (Nandi A et al., 1997, J Cell Biochem 66:500-511), and identified its epitope to a 63-amino acid stretch in the intracellular domain of GCC. In view of the potential that this antibody has for the identification of colorectal tumors, we have characterized the epitope for GCC:B10 in this study. Overlapping peptide synthesis indicated that the epitope was contained in the sequence HIPPENIFPLE. This sequence was unique to GCC, and despite a short stretch of homology with serum amyloid protein and pertussis toxin, no cross reactivity was detected. The core epitope was delineated using a random hexameric phage display library, and two categories of sequences were identified, containing either a single, or two adjacent proline residues. No sequence identified by phage display was identical to the epitope present in GCC, indicating that phage sequences represented mimotopes of the native epitope. Alignment of these sequences with HIPPENIFPLE suggested duplication of the recognition motif, which was confirmed by peptide synthesis. These studies allowed us not only to define the requirements of epitope recognition by GCC:B10 monoclonal antibody, but also to describe a novel means of epitope recognition involving topological mimicry and probable duplication of the cognate epitope in the native guanylyl cyclase C receptor sequence.

Expression of the extracellular domain of the human heat-stable enterotoxin receptor in Escherichia coli and generation of neutralizing antibodies

The entire extracellular domain of the human heat-stable enterotoxin (ST) receptor as well as a truncated N-terminal domain were cloned as glutathione S-transferase fusion proteins and expressed in Escherichia coli. The recombinant fusion proteins were purified from both the cytosol and the inclusion body fractions by selective detergent extraction followed by glutathione-agarose affinity chromatography. The purified protein, corresponding to the entire extracellular domain, bound the stable toxin peptide with an affinity comparable to that of the native receptor characterized from the human colonic T84 cell line. No binding was observed with the N-terminal truncated fragment of the receptor under similar conditions, Polyclonal antibodies were raised to the entire extracellular domain fusion protein as well as the truncated extracellular domain fusion protein, and the antibodies were purified by affinity chromatography. Addition of the purified antibodies to T84 cells inhibited ST binding and abolished ST-mediated cGMP production, indicating that critical epitopes involved in ligand interaction are present in the N-terminal fragment of the receptor, Purified antibodies recognized a single protein of M(r) 160,000 Da on Western blotting with T84 membranes, corresponding to a size of the native glycosylated receptor in T84 cells. These studies are the first report of the expression, purification, and characterization of any member of the guanylyl cyclase family of receptors in E. coli and show that binding of the toxin to the extracellular domain of the receptor is possible in the absence of any posttranslational modifications such as glycosylation. The recombinant fusion proteins as well as the antibodies that we have generated could serve as useful tools in the identification of critical residues of the extracellular domain involved in ligand interaction.

Comparative kinomics of human and chimpanzee reveal unique kinship and functional diversity generated by new domain combinations

Background: Phosphorylation by protein kinases is a common event in many cellular processes. Further, many kinases perform specialized roles and are regulated by non-kinase domains tethered to kinase domain. Perturbation in the regulation of kinases leads to malignancy. We have identified and analysed putative protein kinases encoded in the genome of chimpanzee which is a close evolutionary relative of human. Result: The shared core biology between chimpanzee and human is characterized by many orthologous protein kinases which are involved in conserved pathways. Domain architectures specific to chimp/human kinases have been observed. Chimp kinases with unique domain architectures are characterized by deletion of one or more non-kinase domains in the human kinases. Interestingly, counterparts of some of the multi-domain human kinases in chimp are characterized by identical domain architectures but with kinase-like non-kinase domain. Remarkably, out of 587 chimpanzee kinases no human orthologue with greater than 95% sequence identity could be identified for 160 kinases. Variations in chimpanzee kinases compared to human kinases are brought about also by differences in functions of domains tethered to the catalytic kinase domain. For example, the heterodimer forming PB1 domain related to the fold of ubiquitin/Ras-binding domain is seen uniquely tethered to PKC-like chimpanzee kinase. Conclusion: Though the chimpanzee and human are evolutionary very close, there are chimpanzee kinases with no close counterpart in the human suggesting differences in their functions. This analysis provides a direction for experimental analysis of human and chimpanzee protein kinases in order to enhance our understanding on their specific biological roles.