The Parkinson disease gene LRRK2: evolutionary and structural insights

Mutations in the human leucine-rich repeat kinase 2 (LRRK2) gene are associated with both familial and sporadic Parkinson disease (PD). LRRK2 belongs to a gene family known as Roco. Roco genes encode for large proteins with several protein domains. Particularly, all Roco proteins have a characteristic GTPase domain, named Roc, plus a domain of unknown function called COR. In addition, LRRK2 and several other Roco proteins also contain a protein kinase domain. In this study, I use a combination of phylogenetic and structural analyses of the COR, Roc, and kinase domains present in Roco proteins to describe the origin and evolutionary history of LRRK2. Phylogenetic analyses using these domains demonstrate that LRRK2 emerged from a duplication that occurred after the protostome-deuterostome split. The duplication was followed by the acquisition by LRRK2 proteins of a specific type of N-terminal repeat, described here for the first time. This repeat is absent in the proteins encoded by the paralogs of LRRK2, called LRRK1 or in protostome LRRK proteins. These results suggest that Drosophila or Caenorhabditis LRRK genes may not be good models to understand human LRRK2 function. Genes in the slime mold Dictyostelium discoideum with structures very similar to those found in animal LRRK genes, including the protein kinase domain, have been described. However, phylogenetic analyses suggest that this structural similarity is due to independent acquisitions of distantly related protein kinase domains. Finally, I confirm in an extensive sequence analysis that the Roc GTPase domain is related but still substantially different from small GTPases, such as Rab, Ras, or Rho. Modeling based on known kinase structures suggests that mutations in LRRK2 that cause familiar PD may alter the local 3-dimensional folding of the LRRK2 protein without affecting its overall structure.     

Gene structure and function of tyrosine kinases in the marine sponge Geodia cydonium: autapomorphic characters in Metazoa

Porifera (sponges) represent the most ancient, extant metazoan phylum. They existed already prior to the 'Cambrian Explosion'. Based on the analysis of aa sequences of informative proteins, it is highly likely that all metazoan phyla evolved from only one common ancestor (monophyletic origin). As 'autapomorphic' proteins which are restricted to Metazoa only, integrin receptors, receptors with scavenger receptor cysteine-rich repeats, neuronal-like receptors and protein-tyrosine kinases (PTKs) have been identified in Porifera. From the marine sponge Geodia cydonium, a receptor tyrosine kinase (RTK) has been cloned that comprises the characteristic structural topology known from other metazoan RTKs; an extracellular domain, the transmembrane region, the juxtamembrane region and the TK domain. Only two introns, within the coding region of the RTK gene, could be found, which separate the two highly polymorphic immunoglobulin-like domains, found in the extracellular region of the enzyme. The functional role of this sponge RTK could be demonstrated both in situ (grafting experiments) and in vitro (increase of intracellular Ca2+ level). Upstream of this RTK gene, two further genes coding for tyrosine kinases (TK) have been identified. Both are intron-free. The deduced aa sequence of the first gene shows no transmembrane segment; from the second gene--so far--only half of its catalytic domain is known. A phylogenetic analysis with the TK domains from these sequences and a fourth, from a novel scavenger RTK (all domains comprise the signature for the TK class II receptors), showed that they are distantly related to the insulin and insulin-like receptors. The presented findings support the 'introns-late' hypothesis for such genes that encode 'metazoan' proteins. It is proposed that the TKs evolved from protein-serine/threonine kinases through modularization and subsequent exon shuffling. After formation of the ancestral TKs, the modules lost the framing introns to protect the evolutionary novelty. Since cell culture systems of sponges are now available, it can be expected that soon also those mechanisms that control the developmental programs will be unravelled.     

Structural characterization of the LEM motif common to three human inner nuclear membrane proteins.

BACKGROUND: Integral membrane proteins of the inner nuclear membrane are involved in chromatin organization and postmitotic reassembly of the nucleus. The discovery that mutations in the gene encoding emerin causes X-linked Emery-Dreifuss muscular dystrophy has enhanced interest in such proteins. A common structural domain of 50 residues, called the LEM domain, has been identified in emerin MAN1, and lamina-associated polypeptide (LAP) 2. In particular, all LAP2 isoforms share an N-terminal segment composed of such a LEM domain that is connected to a highly divergent LEM-like domain by a linker that is probably unstructured. RESULTS: We have determined the three-dimensional structures of the LEM and LEM-like domains of LAP2 using nuclear magnetic resonance and molecular modeling. Both domains adopt the same fold, mainly composed of two large parallel alpha helices. CONCLUSIONS: The structural LEM motif is found in human inner nuclear membrane proteins and in protein-protein interaction domains from bacterial multienzyme complexes. This suggests that LEM and LEM-like domains are protein-protein interaction domains. A region conserved in all LEM domains, at the surface of helix 2, could mediate interaction between LEM domains and a common protein partner.     

Diversity and evolution of the thyroglobulin type-1 domain superfamily

Multidomain proteins are gaining increasing consideration for their puzzling, flexible utilization in nature. The presence of the characteristic thyroglobulin type-1 (Tg1) domain as a protein module in a variety of multicellular organisms suggests pivotal roles for this building block. To gain insight into the evolution of Tg1 domains, we performed searches of protein, expressed sequence tag, and genome databases. Tg1 domains were found to be Metazoa specific, and we retrieved a total of 170 Tg1 domain-containing protein sequences. Their architectures revealed a wide taxonomic distribution of proteins containing Tg1 domains followed or preceded by secreted protein, acidic, rich in cysteines (SPARC)-type extracellular calcium-binding domains. Other proteins contained lineage-specific domain combinations of peptidase inhibitory modules or domains with different biological functions. Phylogenetic analysis showed that Tg1 domains are highly conserved within protein structures, whereas insertion into novel proteins is followed by rapid diversification. Seven different basic types of protein architecture containing the Tg1 domain were identified in vertebrates. We examined the evolution of these protein groups by combining Tg1 domain phylogeny with additional analyses based on other characteristic domains. Testicans and secreted modular calcium binding protein (SMOCs) evolved from invertebrate homologs by introduction of vertebrate-specific domains, nidogen evolved by insertion of a Tg1 domain into a preexisting architecture, and the remaining four have unique architectures. Thyroglobulin, Trops, and the major histocompatibility complex class II-associated invariant chain are vertebrate specific, while an insulin-like growth factor-binding protein and nidogen were also identified in urochordates. Among vertebrates, we observed differences in protein repertoires, which result from gene duplication and domain duplication. Members of five groups have been characterized at the molecular level. All exhibit subtle differences in their specificities and function either as peptidase inhibitors (thyropins), substrates, or both. As far as the sequence is concerned, only a few conserved residues were identified. In combination with structural data, our analysis shows that the Tg1 domain fold is highly adaptive and comprises a relatively well-conserved core surrounded by highly variable loops that account for its multipurpose function in the animal kingdom.     

Solution structure of the RWD domain of the mouse GCN2 protein

GCN2 is the alpha-subunit of the only translation initiation factor (eIF2alpha) kinase that appears in all eukaryotes. Its function requires an interaction with GCN1 via the domain at its N-terminus, which is termed the RWD domain after three major RWD-containing proteins: RING finger-containing proteins, WD-repeat-containing proteins, and yeast DEAD (DEXD)-like helicases. In this study, we determined the solution structure of the mouse GCN2 RWD domain using NMR spectroscopy. The structure forms an alpha + beta sandwich fold consisting of two layers: a four-stranded antiparallel beta-sheet, and three side-by-side alpha-helices, with an alphabetabetabetabetaalphaalpha topology. A characteristic YPXXXP motif, which always occurs in RWD domains, forms a stable loop including three consecutive beta-turns that overlap with each other by two residues (triple beta-turn). As putative binding sites with GCN1, a structure-based alignment allowed the identification of several surface residues in alpha-helix 3 that are characteristic of the GCN2 RWD domains. Despite the apparent absence of sequence similarity, the RWD structure significantly resembles that of ubiquitin-conjugating enzymes (E2s), with most of the structural differences in the region connecting beta-strand 4 and alpha-helix 3. The structural architecture, including the triple beta-turn, is fundamentally common among various RWD domains and E2s, but most of the surface residues on the structure vary. Thus, it appears that the RWD domain is a novel structural domain for protein-binding that plays specific roles in individual RWD-containing proteins.     

The generation of new protein functions by the combination of domains

During evolution, many new proteins have been formed by the process of gene duplication and combination. The genes involved in this process usually code for whole domains. Small proteins contain one domain; medium and large proteins contain two or more domains. We have compared homologous domains that occur in both one-domain proteins and multidomain proteins. We have determined (1) how the functions of the individual domains in the multidomain proteins combine to produce their overall functions and (2) the extent to which these functions are similar to those in the one-domain homologs. We describe how domain combinations increase the specificity of enzymes; act as links between domains that have functional roles; regulate activity; combine within one chain functions that can act either independently, in concert or in new contexts; and provide the structural framework for the evolution of entirely new functions.     

Intra-chain 3D segment swapping spawns the evolution of new multidomain protein architectures

Multidomain proteins form in evolution through the concatenation of domains, but structural domains may comprise multiple segments of the chain. In this work, we demonstrate that new multidomain architectures can evolve by an apparent three-dimensional swap of segments between structurally similar domains within a single-chain monomer. By a comprehensive structural search of the current Protein Data Bank (PDB), we identified 32 well-defined segment-swapped proteins (SSPs) belonging to 18 structural families. Nearly 13% of all multidomain proteins in the PDB may have a segment-swapped evolutionary precursor as estimated by more permissive searching criteria. The formation of SSPs can be explained by two principal evolutionary mechanisms: (i) domain swapping and fusion (DSF) and (ii) circular permutation (CP). By large-scale comparative analyses using structural alignment and hidden Markov model methods, it was found that the majority of SSPs have evolved via the DSF mechanism, and a much smaller fraction, via CP. Functional analyses further revealed that segment swapping, which results in two linkers connecting the domains, may impart directed flexibility to multidomain proteins and contributes to the development of new functions. Thus, inter-domain segment swapping represents a novel general mechanism by which new protein folds and multidomain architectures arise in evolution, and SSPs have structural and functional properties that make them worth defining as a separate group.     

Tau binds and organizes Escherichia coli replication proteins through distinct domains. Domain III, shared by gamma and tau, binds delta delta ' and chi psi

The DnaX complex of the DNA polymerase holoenzyme assembles the beta(2) processivity factor onto the primed template enabling highly processive replication. The key ATPases within this complex are tau and gamma, alternative frameshift products of the dnaX gene. Of the five domains of tau, I-III are shared with gamma In vivo, gamma binds the auxiliary subunits deltadelta' and chipsi (Glover, B. P., and McHenry, C. S. (2000) J. Biol. Chem. 275, 3017-3020). To localize deltadelta' and chipsi binding domains within gamma domains I-III, we measured the binding of purified biotin-tagged DnaX proteins lacking specific domains to deltadelta' and chipsi by surface plasmon resonance. Fusion proteins containing either DnaX domains I-III or domains III-V bound deltadelta' and chipsi subunits. A DnaX protein only containing domains I and II did not bind deltadelta' or chipsi. The binding affinity of chipsi for DnaX domains I-III and domains III-V was the same as that of chipsi for full-length tau, indicating that domain III contained all structural elements required for chipsi binding. Domain III of tau also contained deltadelta' binding sites, although the interaction between deltadelta' and domains III-V of tau was 10-fold weaker than the interaction between deltadelta' and full length tau. The presence of both delta and chipsi strengthened the delta'-C(0)tau interaction by at least 15-fold. Domain III was the only domain common to all of tau fusion proteins whose interaction with delta' was enhanced in the presence of delta and chipsi. Thus, domain III of the DnaX proteins not only contains the deltadelta' and chipsi binding sites but also contains the elements required for the positive cooperative assembly of the DnaX complex.     

The structure and dynamics of tandem WW domains in a negative regulator of notch signaling, Suppressor of deltex

WW domains mediate protein recognition, usually though binding to proline-rich sequences. In many proteins, WW domains occur in tandem arrays. Whether or how individual domains within such arrays cooperate to recognize biological partners is, as yet, poorly characterized. An important question is whether functional diversity of different WW domain proteins is reflected in the structural organization and ligand interaction mechanisms of their multiple domains. We have determined the solution structure and dynamics of a pair of WW domains (WW3-4) from a Drosophila Nedd4 family protein called Suppressor of deltex (Su(dx)), a regulator of Notch receptor signaling. We find that the binding of a type 1 PPPY ligand to WW3 stabilizes the structure with effects propagating to the WW4 domain, a domain that is not active for ligand binding. Both WW domains adopt the characteristic triple-stranded beta-sheet structure, and significantly, this is the first example of a WW domain structure to include a domain (WW4) lacking the second conserved Trp (replaced by Phe). The domains are connected by a flexible linker, which allows a hinge-like motion of domains that may be important for the recognition of functionally relevant targets. Our results contrast markedly with those of the only previously determined three-dimensional structure of tandem WW domains, that of the rigidly oriented WW domain pair from the RNA-splicing factor Prp40. Our data illustrate that arrays of WW domains can exhibit a variety of higher order structures and ligand interaction mechanisms.     

Functional fingerprints of folds: evidence for correlated structure-function evolution

Using structural similarity clustering of protein domains: protein domain universe graph (PDUG), and a hierarchical functional annotation: gene ontology (GO) as two evolutionary lenses, we find that each structural cluster (domain fold) exhibits a distribution of functions that is unique to it. These functional distributions are functional fingerprints that are specific to characteristic structural clusters and vary from cluster to cluster. Furthermore, as structural similarity threshold for domain clustering in the PDUG is relaxed we observe an influx of earlier-diverged domains into clusters. These domains join clusters without destroying the functional fingerprint. These results can be understood in light of a divergent evolution scenario that posits correlated divergence of structural and functional traits in protein domains from one or few progenitors.     

Structure and function of human tissue-type plasminogen activator (t-PA)

Full-length tissue-type plasminogen activator (t-PA) cDNA served to construct deletion mutants within the N-terminal "heavy" (H)-chain of the t-PA molecule. The H-chain cDNA consists of an array of structural domains homologous to domains present on other plasma proteins ("finger," "epidermal growth factor," "kringles"). These structural domains have been located on an exon or a set of exons. The endpoints of the deletions nearly coincide with exon-intron junctions of the chromosomal t-PA gene. Recombinant t-PA deletion mutant proteins were obtained after transient expression in mouse Ltk- cells, transfected with SV40-pBR322-derived t-PA cDNA plasmids. It is demonstrated that the serine protease moiety of t-PA and its substrate specificity for plasminogen is entirely contained within the C-terminal "light" (L)-chain of the protein. The presence of cDNA, encoding the t-PA signal peptide preceding the remaining portion of t-PA, suffices to achieve secretion of (mutant) t-PA into the medium. The stimulatory effect of fibrin on the plasminogen activator activity of t-PA was shown to be mediated by the kringle K2 domain and, to a lesser extent, by the finger domain. The other domains on the H-chain, kringle K1, and the epidermal growth-factor-like domain, do not contribute to this property of t-PA. These findings correlate well with the fibrin-binding properties of the rt-PA deletion-mutant proteins, indicating that stimulation of the activity is based on aligning of the substrate plasminogen and its enzyme t-PA on the fibrin matrix. The primary target for endothelial plasminogen activator inhibitor (PAI) is located within the L-chain of t-PA. Deleting specific segments of t-PA H-chain cDNA and subsequent transient expression in mouse Ltk- cells of t-PA deletion-mutant proteins did not affect the formation of a stable complex between mutant t-PA and PAI.     

Cross-species analyses identify the BNIP-2 and Cdc42GAP homology (BCH) domain as a distinct functional subclass of the CRAL_TRIO/Sec14 superfamily

The CRAL_TRIO protein domain, which is unique to the Sec14 protein superfamily, binds to a diverse set of small lipophilic ligands. Similar domains are found in a range of different proteins including neurofibromatosis type-1, a Ras GTPase-activating Protein (RasGAP) and Rho guanine nucleotide exchange factors (RhoGEFs). Proteins containing this structural protein domain exhibit a low sequence similarity and ligand specificity while maintaining an overall characteristic three-dimensional structure. We have previously demonstrated that the BNIP-2 and Cdc42GAP Homology (BCH) protein domain, which shares a low sequence homology with the CRAL_TRIO domain, can serve as a regulatory scaffold that binds to Rho, RhoGEFs and RhoGAPs to control various cell signalling processes. In this work, we investigate 175 BCH domain-containing proteins from a wide range of different organisms. A phylogenetic analysis with ~100 CRAL_TRIO and similar domains from eight representative species indicates a clear distinction of BCH-containing proteins as a novel subclass within the CRAL_TRIO/Sec14 superfamily. BCH-containing proteins contain a hallmark sequence motif R(R/K)h(R/K)(R/K)NL(R/K)xhhhhHPs ('h' is large and hydrophobic residue and 's' is small and weekly polar residue) and can be further subdivided into three unique subtypes associated with BNIP-2-N, macro- and RhoGAP-type protein domains. A previously unknown group of genes encoding 'BCH-only' domains is also identified in plants and arthropod species. Based on an analysis of their gene-structure and their protein domain context we hypothesize that BCH domain-containing genes evolved through gene duplication, intron insertions and domain swapping events. Furthermore, we explore the point of divergence between BCH and CRAL-TRIO proteins in relation to their ability to bind small GTPases, GAPs and GEFs and lipid ligands. Our study suggests a need for a more extensive analysis of previously uncharacterized BCH, 'BCH-like' and CRAL_TRIO-containing proteins and their significance in regulating signaling events involving small GTPases.     

Light-modulated abundance of an mRNA encoding a calmodulin-regulated, chromatin-associated NTPase in pea

A cDNA encoding a 47 kDa nucleoside triphosphatase (NTPase) that is associated with the chromatin of pea nuclei has been cloned and sequenced. The translated sequence of the cDNA includes several domains predicted by known biochemical properties of the enzyme, including five motifs characteristic of the ATP-binding domain of many proteins, several potential casein kinase II phosphorylation sites, a helix-turn-helix region characteristic of DNA-binding proteins, and a potential calmodulin-binding domain. The deduced primary structure also includes an N-terminal sequence that is a predicted signal peptide and an internal sequence that could serve as a bipartite-type nuclear localization signal. Both in situ immunocytochemistry of pea plumules and immunoblots of purified cell fractions indicate that most of the immunodetectable NTPase is within the nucleus, a compartment proteins typically reach through nuclear pores rather than through the endoplasmic reticulum pathway. The translated sequence has some similarity to that of human lamin C, but not high enough to account for the earlier observation that IgG against human lamin C binds to the NTPase in immunoblots. Northern blot analysis shows that the NTPase mRNA is strongly expressed in etiolated plumules, but only poorly or not at all in the leaf and stem tissues of light-grown plants. Accumulation of NTPase mRNA in etiolated seedlings is stimulated by brief treatments with both red and far-red light, as is characteristic of very low-fluence phytochrome responses. Southern blotting with pea genomic DNA indicates the NTPase is likely to be encoded by a single gene.     

Schistosoma mansoni: the dicer gene and its expression

RNA interference (RNAi) is a gene silencing mechanism that plays an important role in regulating gene expression in many eukaryotes and has become a valuable molecular tool for analyzing gene function. Multi-domain nucleases called Dicer proteins play pivotal roles in RNAi. In this paper, we characterize the structure and expression of the Dicer gene from the platyhelminth parasite Schistosoma mansoni. The gene (SmDicer) is over 54kb long and comprises 30 exons that potentially encode a 2641 amino acid protein. This is the largest Dicer protein yet described. SmDicer contains all domains that are characteristic of metazoan dicers including an amino terminal helicase domain, DUF283, a PAZ domain, two RNAse III domains and an RNA binding domain. An examination of the available S. mansoni genome sequence suggests that the Dicer gene described here is the only Dicer gene in the parasite genome. SmDicer is expressed throughout schistosome development suggesting that RNAi technologies might be employed in deciphering gene function in all life stages of this parasite.     

Identifying protein domains by global analysis of soluble fragment data

The production and analysis of individual structural domains is a common strategy for studying large or complex proteins, which may be experimentally intractable in their full-length form. However, identifying domain boundaries is challenging if there is little structural information concerning the protein target. One experimental procedure for mapping domains is to screen a library of random protein fragments for solubility, since truncation of a domain will typically expose hydrophobic groups, leading to poor fragment solubility. We have coupled fragment solubility screening with global data analysis to develop an effective method for identifying structural domains within a protein. A gene fragment library is generated using mechanical shearing, or by uracil doping of the gene and a uracil-specific enzymatic digest. A split green fluorescent protein (GFP) assay is used to screen the corresponding protein fragments for solubility when expressed in Escherichia coli. The soluble fragment data are then analyzed using two complementary approaches. Fragmentation “hotspots” indicate possible interdomain regions. Clustering algorithms are used to group related fragments, and concomitantly predict domain location. The effectiveness of this Domain Seeking procedure is demonstrated by application to the well-characterized human protein p85α.