Structure of the Roc-COR domain tandem of C. tepidum, a prokaryotic homologue of the human LRRK2 Parkinson kinase

Ras of complex proteins (Roc) belongs to the superfamily of Ras-related small G-proteins that always occurs in tandem with the C-terminal of Roc (COR) domain. This Roc-COR tandem is found in the bacterial and eukaryotic world. Its most prominent member is the leucine-rich repeat kinase LRRK2, which is mutated and activated in Parkinson patients. Here, we investigated biochemically and structurally the Roco protein from Chlorobium tepidum. We show that Roc is highly homologous to Ras, whereas the COR domain is a dimerisation device. The juxtaposition of the G-domains and mutational analysis suggest that the Roc GTPase reaction is stimulated and/or regulated by dimerisation in a nucleotide-dependent manner. The region most conserved between bacteria and man is the interface between Roc and COR, where single-point Parkinson mutations of the Roc and COR domains are in close proximity. The analogous mutations in C. tepidum Roc-COR decrease the GTPase reaction rate, most likely due to a modification of the interaction between the Roc and COR domains.     

Conformational heterogeneity of the Roc domains in C. tepidum Roc–COR and implications for human LRRK2 Parkinson mutations

Ras of complex proteins (Roc) is a Ras-like GTP-binding domain that always occurs in tandem with the C-terminal of Roc (COR) domain and is found in bacteria, plants and animals. Recently, it has been shown that Roco proteins belong to the family of G-proteins activated by nucleotide (nt)-dependent dimerization (GADs). We investigated the RocCOR tandem from the bacteria Chlorobium tepidum with site-directed spin labelling and pulse EPR distance measurements to follow conformational changes during the Roco G-protein cycle. Our results confirm that the COR domains are a stable dimerization device serving as a scaffold for the Roc domains that, in contrast, are structurally heterogeneous and dynamic entities. Contrary to other GAD proteins, we observed only minor structural alterations upon binding and hydrolysis of GTP, indicating significant mechanistic variations within this protein class. Mutations in the most prominent member of the Roco family of proteins, leucine-rich repeat (LRR) kinase 2 (LRRK2), are the most frequent cause of late-onset Parkinson''s disease (PD). Using a stable recombinant LRRK2 Roc-COR-kinase fragment we obtained detailed kinetic data for the G-protein cycle. Our data confirmed that dimerization is essential for efficient GTP hydrolysis and PD mutations in the Roc domain result in decreased GTPase activity. Previous data have shown that these LRRK2 PD-mutations are located in the interface between Roc and COR. Importantly, analogous mutations in the conserved C. tepidum Roc/COR interface significantly influence the structure and nt-induced conformational changes of the Roc domains.     

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.     

Intramolecular activation mechanism of the Dictyostelium LRRK2 homolog Roco protein GbpC

GbpC is a large multidomain protein involved in cGMP-mediated chemotaxis in the cellular slime mold Dictyostelium discoideum. GbpC belongs to the Roco family of proteins that often share a central core region, consisting of leucine-rich repeats, a Ras domain (Roc), a Cor domain, and a MAPKKKinase domain. In addition to this core, GbpC contains a RasGEF domain and two cGMP-binding domains. Here, we report on an intramolecular signaling cascade of GbpC. In vitro, the RasGEF domain of GbpC specifically accelerates the GDP/GTP exchange of the Roc domain. Moreover, cGMP binding to GbpC strongly stimulates the binding of GbpC to GTP-agarose, suggesting cGMP-stimulated GDP/GTP exchange at the Roc domain. The function of the protein in vivo was investigated by rescue analysis of the chemotactic defect of gbpC null cells. Mutants that lack a functional guanine exchange factor (GEF), Roc, or kinase domain are inactive in vivo. Together, the results suggest a four-step intramolecular activation mechanism of the Roco protein GbpC: cGMP binding to the cyclic nucleotide-binding domains, activation of the GEF domain, GDP/GTP exchange of Roc, and activation of the MAPKKK domain.     

Leucine-rich repeat kinase 2: relevance to Parkinson's disease

Human leucine-rich repeat kinase 2 (LRRK2) is a novel kinase belonging to the ROCO protein superfamily (Ras of complex proteins (Roc) with a C-terminal of Roc domain). This large complex protein of 280kDa contains several functional domains including leucine-rich repeats, Ras-related GTPase, mitogen-activated protein kinase kinase kinase (MAPKKK), and WD40 repeats. While definitive functions of LRRK2 have yet to be described, the domain structure of LRRK2 suggests that it plays an important role in the regulation of signal transduction cascades through its dual enzymatic activities of GTPase and MAPKKK. Moreover, mutations in LRRK2 have been found to be thus far the most frequent cause of late-onset familial and idiopathic Parkinson's disease. Further investigations should allow for the elucidation of how pathogenic mutations trigger changes in the structure and function of LRRK2 that lead to aberrant signal transduction and neurodegeneration in Parkinson's disease.     

Sequence conservation between porcine and human LRRK2

Leucine-rich repeat kinase 2 (LRRK2) is a member of the ROCO protein superfamily (Ras of complex proteins (Roc) with a C-terminal Roc domain). Mutations in the LRRK2 gene lead to autosomal dominant Parkinsonism. We have cloned the porcine LRRK2 cDNA in an attempt to characterize conserved and therefore likely functional domains. The LRRK2 cDNA contains an open reading frame of 7,578 bp. The predicted LRRK2 protein consists of 2,526 amino acids of 86–93% identity with its mammalian couterparts. The deduced amino acid sequence of encoded porcine LRRK2 protein displays extensive homology with its human counterpart, with greatest similarities in those regions that contain the kinase domain, the Roc domain and the COR motif. Expression of porcine LRRK2 mRNA in various organs and tissues is similar to its human counterpart and not limited to the brain. The obtained data show that the LRRK2 sequence and expression patterns are conserved across species. The porcine LRRK2 gene was mapped to chromosome 5q25. The results obtained suggest that the LRRK2 gene might be of particular interest in our attempt to generate a transgenic porcine model for Parkinson’s disease. The sequence of the porcine LRRK2 cDNA, encoding the LRRK2/dardarin protein, and the genomic sequence of LRRK2 have been submitted to DDBJ/EMBL/GenBank under the Accessions Numbers EU019992, and EU019994, respectively.     

Structural and functional in silico analysis of LRRK2 missense substitutions

The LRRK2 gene (Leucine-Rich Repeat Kinase 2, PARK8) is mutated in a significant number of cases of autosomal dominant Parkinson’s disease (PD) and in some sporadic cases of late-onset PD. LRRK2 is a large, complex protein that comprises several interaction domains: armadillo, ankyrin, leucine-rich repeats and WD40 domains; two catalytic domains: ROC-GTPase and serine/threonine kinase; and a COR domain (unknown function). Pathogenic mutations are scattered all over the domains of LRRK2, although the prevalence of mutations in some domains is higher (ROC-GTPase, COR and kinase). In this work, we model the structure of each domain to predict and explore the effects of described missense mutations and polymorphisms. The results allow us to postulate the possible effects of pathogenic mutations in the function of the protein, and hypothesize the importance of some polymorphisms that have not been linked directly to PD, but act as risk factors for the disease. In our analysis, we also study the effects of PD-related mutations in the kinase domain structure and in the phosphorylation of the activation loop to determine effects on kinase activity.     

Identifying determinants of cullin binding specificity among the three functionally different Drosophila melanogaster Roc proteins via domain swapping

Background

Cullin-dependent E3 ubiquitin ligases (CDL) are key regulators of protein destruction that participate in a wide range of cell biological processes. The Roc subunit of CDL contains an evolutionarily conserved RING domain that binds ubiquitin charged E2 and is essential for ubiquitylation. Drosophila melanogaster contains three highly related Roc proteins: Roc1a and Roc2, which are conserved in vertebrates, and Roc1b, which is specific to Drosophila. Our previous genetic data analyzing Roc1a and Roc1b mutants suggested that Roc proteins are functionally distinct, but the molecular basis for this distinction is not known.

Methodology/Principal Findings

Using co-immunoprecipitation studies we show that Drosophila Roc proteins bind specific Cullins: Roc1a binds Cul1-4, Roc1b binds Cul3, and Roc2 binds Cul5. Through domain swapping experiments, we demonstrate that Cullin binding specificity is strongly influenced by the Roc NH₂-terminal domain, which forms an inter-molecular β sheet with the Cullin. Substitution of the Roc1a RING domain with that of Roc1b results in a protein with similar Cullin binding properties to Roc1a that is active as an E3 ligase but cannot complement Roc1a mutant lethality, indicating that the identity of the RING domain can be an important determinant of CDL function. In contrast, the converse chimeric protein with a substitution of the Roc1b RING domain with that of Roc1a can rescue the male sterility of Roc1b mutants, but only when expressed from the endogenous Roc1b promoter. We also identified mutations of Roc2 and Cul5 and show that they cause no overt developmental phenotype, consistent with our finding that Roc2 and Cul5 proteins are exclusive binding partners, which others have observed in human cells as well.

Conclusions

The Drosophila Roc proteins are highly similar, but have diverged during evolution to bind a distinct set of Cullins and to utilize RING domains that have overlapping, but not identical, function in vivo.     

LRRK2 pathobiology in Parkinson's disease

Mutations in the catalytic Roc‐COR and kinase domains of leucine‐rich repeat kinase 2 (LRRK2) are a common cause of familial Parkinson's disease (PD). LRRK2 mutations cause PD with age‐related penetrance and clinical features identical to late‐onset sporadic PD. Biochemical studies support an increase in LRRK2 kinase activity and a decrease in GTPase activity for kinase domain and Roc‐COR mutations, respectively. Strong evidence exists that LRRK2 toxicity is kinase dependent leading to extensive efforts to identify selective and brain‐permeable LRRK2 kinase inhibitors for clinical development. Cell and animal models of PD indicate that LRRK2 mutations affect vesicular trafficking, autophagy, protein synthesis, and cytoskeletal function. Although some of these biological functions are affected consistently by most disease‐linked mutations, others are not and it remains currently unclear how mutations that produce variable effects on LRRK2 biochemistry and function all commonly result in the degeneration and death of dopamine neurons. LRRK2 is typically present in Lewy bodies and its toxicity in mammalian models appears to be dependent on the presence of α‐synuclein, which is elevated in human iPS‐derived dopamine neurons from patients harboring LRRK2 mutations. Here, we summarize biochemical and functional studies of LRRK2 and its mutations and focus on aberrant vesicular trafficking and protein synthesis as two leading mechanisms underlying LRRK2‐linked disease.

     

Identification and analysis of novel tandem repeats in the cell surface proteins of archaeal and bacterial genomes using computational tools

We have identified four novel repeats and two domains in cell surface proteins encoded by the Methanosarcina acetivorans genome and in some archaeal and bacterial genomes. The repeats correspond to a certain number of amino acid residues present in tandem in a protein sequence and each repeat is characterized by conserved sequence motifs. These correspond to: (a) a 42 amino acid (aa) residue RIVW repeat; (b) a 45 aa residue LGxL repeat; (c) a 42 aa residue LVIVD repeat; and (d) a 54 aa residue LGFP repeat. The domains correspond to a certain number of aa residues in a protein sequence that do not comprise internal repeats. These correspond to: (a) a 200 aa residue DNRLRE domain; and (b) a 70 aa residue PEGA domain. We discuss the occurrence of these repeats and domains in the different proteins and genomes analysed in this work.     

FIST: a sensory domain for diverse signal transduction pathways in prokaryotes and ubiquitin signaling in eukaryotes

MOTIVATION: Sensory domains that are conserved among Bacteria, Archaea and Eucarya are important detectors of common signals detected by living cells. Due to their high sequence divergence, sensory domains are difficult to identify. We systematically look for novel sensory domains using sensitive profile-based searches initiated with regions of signal transduction proteins where no known domains can be identified by current domain models. RESULTS: Using profile searches followed by multiple sequence alignment, structure prediction and domain architecture analysis, we have identified a novel sensory domain termed FIST, which is present in signal transduction proteins from Bacteria, Archaea and Eucarya. Chromosomal proximity of FIST-encoding genes to those coding for proteins involved in amino acid metabolism and transport suggest that FIST domains bind small ligands, such as amino acids.     

Posttranslational processing and differential glycosylation of Tractin, an Ig-superfamily member involved in regulation of axonal outgrowth

Tractin is a novel member of the Ig-superfamily which has a highly unusual structure. It contains six Ig domains, four FNIII-like domains, an acidic domain, 12 repeats of a novel proline- and glycine-rich motif with sequence similarity to collagen, a transmembrane domain, and an intracellular tail with an ankyrin and a PDZ domain binding motif. By generating domain-specific antibodies, we show that Tractin is proteolytically processed at two cleavage sites, one located in the third FNIII domain, and a second located just proximal to the transmembrane domain resulting in the formation of four fragments. The most NH(2)-terminal fragment which is glycosylated with the Lan3-2, Lan4-2, and Laz2-369 glycoepitopes is secreted, and we present evidence which supports a model in which the remaining fragments combine to form a secreted homodimer as well as a transmembrane heterodimer. The extracellular domain of the dimers is mostly made up of the collagen-like PG/YG-repeat domain but also contains 11/2 FNIII domain and the acidic domain. The collagen-like PG/YG-repeat domain could be selectively digested by collagenase and we show by yeast two-hybrid analysis that the intracellular domain of Tractin can interact with ankyrin. Thus, the transmembrane heterodimer of Tractin constitutes a novel protein domain configuration where sequence that has properties similar to that of extracellular matrix molecules is directly linked to the cytoskeleton through interactions with ankyrin.     

Cloning of a novel member of the low-density lipoprotein receptor family

A gene encoding a novel transmembrane protein was identified by DNA sequence analysis within the insulin-dependent diabetes mellitus (IDDM) locus IDDM4 on chromosome 11q13. Based on its chromosomal position, this gene is a candidate for conferring susceptibility to diabetes. The gene, termed low-density lipoprotein receptor related protein 5 (LRP5), encodes a protein of 1615 amino acids that contains conserved modules which are characteristic of the low-density lipoprotein (LDL) receptor family. These modules include a putative signal peptide for protein export, four epidermal growth factor (EGF) repeats with associated spacer domains, three LDL-receptor (LDLR) repeats, a single transmembrane spanning domain, and a cytoplasmic domain. The encoded protein has a unique organization of EGF and LDLR repeats; therefore, LRP5 likely represents a new category of the LDLR family. Both human and mouse LRP5 cDNAs have been isolated and the encoded mature proteins are 95% identical, indicating a high degree of evolutionary conservation.     

Conservation of inter-residue interactions and prediction of folding rates of domain repeats

Domains are the main structural and functional units of larger proteins. They tend to be contiguous in primary structure and can fold and function independently. It has been observed that 10–20% of all encoded proteins contain duplicated domains and the average pairwise sequence identity between them is usually low. In the present study, we have analyzed the structural similarity between domain repeats of proteins with known structures available in the Protein Data Bank using structure-based inter-residue interaction measures such as the number of long-range contacts, surrounding hydrophobicity, and pairwise interaction energy. We used RADAR program for detecting the repeats in a protein sequence which were further validated using Pfam domain assignments. The sequence identity between the repeats in domains ranges from 20 to 40% and their secondary structural elements are well conserved. The number of long-range contacts, surrounding hydrophobicity calculations and pairwise interaction energy of the domain repeats clearly reveal the conservation of 3-D structure environment in the repeats of domains. The proportions of mainchain–mainchain hydrogen bonds and hydrophobic interactions are also highly conserved between the repeats. The present study has suggested that the computation of these structure-based parameters will give better clues about the tertiary environment of the repeats in domains. The folding rates of individual domains in the repeats predicted using the long-range order parameter indicate that the predicted folding rates correlate well with most of the experimentally observed folding rates for the analyzed independently folded domains.     

Identification and analysis of novel amino acid sequence repeats and domains in Pyrobaculum aerophilum using computational tools

We have identified four repeats and five domains that are novel in proteins encoded by the Pyrobaculum aerophilum str. IM2 proteome using automated in silico methods. A "repeat" corresponds to a region comprising less than 55 amino acid residues that occurs more than once in the protein sequence and sometimes present in tandem. A "domain" corresponds to a conserved region comprising greater than 55 amino acid residues and may be present as single or multiple copies in the protein sequence. These correspond to (1) 85 amino acid residues AAG domain, (2) 72 amino acid residues GFGN domain, (3) 43 amino acid residues KGG repeat, (4) 25 amino acid residues RWE repeat, (5) 25 amino acid residues RID repeat, (6) 108 amino acid residues NDFA domain, (7) 140 amino acid residues VxY domain, (8) 35 amino acid residues LLPN repeat and (9) 98 amino acid residues GxY domain. A repeat or domain is characterized by specific conserved sequence motifs. We discuss the presence of these repeats and domains in proteins from other genomes and their probable secondary structure.     

The glycosaminoglycan attachment regions of human aggrecan

Aggrecan possesses both chondroitin sulfate (CS) and keratan sulfate (KS) chains attached to its core protein, which reside mainly in the central region of the molecule termed the glycosaminoglycan-attachment region. This region is further subdivided into the KS-rich domain and two adjacent CS-rich domains (CS1 and CS2). The CS1 domain of the human is unique in exhibiting length polymorphism due to a variable number of tandem amino acid repeats. The focus of this work was to determine how length polymorphism affects the structure of the CS1 domain and whether CS and KS chains can coexist in the different glycosaminoglycan-attachment domains. The CS1 domain possesses several amino acid repeat sequences that divide it into three subdomains. Variation in repeat number may occur in any of these domains, with the consequence that CS1 domains of the same length may possess different amino acid sequences. There was no evidence to support the presence of KS in either the CS1 or the CS2 domains nor the presence of CS in the KS-rich domain. The structure of the CS chains was shown to vary between the CS1 and CS2 domains, particularly in the adult, with variation occurring in chain length and the sulfation of the non-reducing terminal N-acetyl galactosamine residue. CS chains in the adult CS2 domain were shorter than those in the CS1 domain and possessed disulfated terminal residues in addition to monosulfated residues. There was, however, no change in the sulfation pattern of the disaccharide repeats in the CS chains from the two domains.     

LRRK2 phosphorylates moesin at threonine-558: characterization of how Parkinson's disease mutants affect kinase activity

Mutations in the LRRK2 (leucine-rich repeat kinase-2) gene cause late-onset PD (Parkinson's disease). LRRK2 contains leucine-rich repeats, a GTPase domain, a COR [C-terminal of Roc (Ras of complex)] domain, a kinase and a WD40 (Trp-Asp 40) motif. Little is known about how LRRK2 is regulated, what its physiological substrates are or how mutations affect LRRK2 function. Thus far LRRK2 activity has only been assessed by autophosphorylation and phosphorylation of MBP (myelin basic protein), which is catalysed rather slowly. We undertook a KESTREL (kinase substrate tracking and elucidation) screen in rat brain extracts to identify proteins that were phosphorylated by an activated PD mutant of LRRK2 (G2019S). This led to the discovery that moesin, a protein which anchors the actin cytoskeleton to the plasma membrane, is efficiently phosphorylated by LRRK2, at Thr558, a previously identified in-vivo-phosphorylation site that regulates the ability of moesin to bind actin. LRRK2 also phosphorylated ezrin and radixin, which are related to moesin, at the residue equivalent to Thr558, as well as a peptide (LRRKtide: RLGRDKYKTLRQIRQ) encompassing Thr558. We exploited these findings to determine how nine previously reported PD mutations of LRRK2 affected kinase activity. Only one of the mutations analysed, namely G2019S, stimulated kinase activity. Four mutations inhibited LRRK2 kinase activity (R1941H, I2012T, I2020T and G2385R), whereas the remainder (R1441C, R1441G, Y1699C and T2356I) did not influence activity. Therefore the manner in which LRRK2 mutations induce PD is more complex than previously imagined and is not only caused by an increase in LRRK2 kinase activity. Finally, we show that the minimum catalytically active fragment of LRRK2 requires an intact GTPase, COR and kinase domain, as well as a WD40 motif and a C-terminal tail. The results of the present study suggest that moesin, ezrin and radixin may be LRRK2 substrates, findings that have been exploited to develop the first robust quantitative assay to measure LRRK2 kinase activity.