Structures of ubiquitin-like (Ubl) and Hsp90-like domains of sacsin provide insight into pathological mutations.

Autosomal recessive spastic ataxia of Charlevoix-Saguenay (ARSACS) is a neurodegenerative disease that is caused by mutations in the SACS gene. The product of this gene is a very large 520-kDa cytoplasmic protein, sacsin, with a ubiquitin-like (Ubl) domain at the N terminus followed by three large sacsin internal repeat (SIRPT) supradomains and C-terminal J and HEPN domains. The SIRPTs are predicted to contain Hsp90-like domains, suggesting a potential chaperone activity. In this work, we report the structures of the Hsp90-like Sr1 domain of SIRPT1 and the N-terminal Ubl domain determined at 1.55- and 2.1-Å resolutions, respectively. The Ubl domain crystallized as a swapped dimer that could be relevant in the context of full-length protein. The Sr1 domain displays the Bergerat protein fold with a characteristic nucleotide-binding pocket, although it binds nucleotides with very low affinity. The Sr1 structure reveals that ARSACS-causing missense mutations (R272H, R272C, and T201K) disrupt protein folding, most likely leading to sacsin degradation. This work lends structural support to the view of sacsin as a molecular chaperone and provides a framework for future studies of this protein.

Structural studies of the yeast DNA damage-inducible protein Ddi1 reveal domain architecture of this eukaryotic protein family.

The eukaryotic Ddi1 family is defined by a conserved retroviral aspartyl protease-like (RVP) domain found in association with a ubiquitin-like (UBL) domain. Ddi1 from Saccharomyces cerevisiae additionally contains a ubiquitin-associated (UBA) domain. The substrate specificity and role of the protease domain in the biological functions of the Ddi family remain unclear. Yeast Ddi1 has been implicated in the regulation of cell cycle progression, DNA-damage repair, and exocytosis. Here, we investigated the multi-domain structure of yeast Ddi1 using X-ray crystallography, nuclear magnetic resonance, and small-angle X-ray scattering. The crystal structure of the RVP domain sheds light on a putative substrate recognition site involving a conserved loop. Isothermal titration calorimetry confirms that both UBL and UBA domains bind ubiquitin, and that Ddi1 binds K48-linked diubiquitin with enhanced affinity. The solution NMR structure of a helical domain that precedes the protease displays tertiary structure similarity to DNA-binding domains from transcription regulators. Our structural studies suggest that the helical domain could serve as a landing platform for substrates in conjunction with attached ubiquitin chains binding to the UBL and UBA domains.

Efficient traversal of beta-sheet protein folding pathways using ensemble models.

Molecular dynamics (MD) simulations can now predict ms-timescale folding processes of small proteins; however, this presently requires hundreds of thousands of CPU hours and is primarily applicable to short peptides with few long-range interactions. Larger and slower-folding proteins, such as many with extended ß-sheet structure, would require orders of magnitude more time and computing resources. Furthermore, when the objective is to determine only which folding events are necessary and limiting, atomistic detail MD simulations can prove unnecessary. Here, we introduce the program tFolder as an efficient method for modelling the folding process of large ß-sheet proteins using sequence data alone. To do so, we extend existing ensemble ß-sheet prediction techniques, which permitted only a fixed anti-parallel ß-barrel shape, with a method that predicts arbitrary ß-strand/ß-strand orientations and strand-order permutations. By accounting for all partial and final structural states, we can then model the transition from random coil to native state as a Markov process, using a master equation to simulate population dynamics of folding over time. Thus, all putative folding pathways can be energetically scored, including which transitions present the greatest barriers. Since correct folding pathway prediction is likely determined by the accuracy of contact prediction, we demonstrate the accuracy of tFolder to be comparable with state-of-the-art methods designed specifically for the contact prediction problem alone. We validate our method for dynamics prediction by applying it to the folding pathway of the well-studied Protein G. With relatively very little computation time, tFolder is able to reveal critical features of the folding pathways which were only previously observed through time-consuming MD simulations and experimental studies. Such a result greatly expands the number of proteins whose folding pathways can be studied, while the algorithmic integration of ensemble prediction with Markovian dynamics can be applied to many other problems.

Self-association studies of the bifunctional N-acetylglucosamine-1-phosphate uridyltransferase from Escherichia coli.

The N-acetylglucosamine-1-phosphate uridyltransferase (GlmU) is a key bifunctional enzyme in the biosynthesis of UDP-GlcNAc, a precursor in the synthesis of cell wall peptidoglycan. Crystal structures of the enzyme from different bacterial strains showed that the polypeptide forms a trimer through a unique parallel left-handed beta helix domain. Here, we show that the GlmU enzyme from Escherichia coli forms a hexamer in solution. Sedimentation equilibrium analytical ultracentrifugation demonstrated that the enzyme is in a trimer/hexamer equilibrium. Small-angle X-ray scattering studies were performed to determine the structure of the hexameric assembly and showed that two trimers assemble through their N-terminal domains. The interaction is mediated by a loop that undergoes a large conformational change in the uridyl transferase reaction, a feature that may affect the enzymatic activity of GlmU.

The three-dimensional architecture of Hox cluster silencing.

Spatial chromatin organization is emerging as an important mechanism to regulate the expression of genes. However, very little is known about genome architecture at high-resolution in vivo. Here, we mapped the three-dimensional organization of the human Hox clusters with chromosome conformation capture (3C) technology. We show that computational modeling of 3C data sets can identify candidate regulatory proteins of chromatin architecture and gene expression. Hox genes encode evolutionarily conserved master regulators of development which strict control has fascinated biologists for over 25 years. Proper transcriptional silencing is key to Hox function since premature expression can lead to developmental defects or human disease. We now show that the HoxA cluster is organized into multiple chromatin loops that are dependent on transcription activity. Long-range contacts were found in all four silent clusters but looping patterns were specific to each cluster. In contrast to the Drosophila homeotic bithorax complex (BX-C), we found that Polycomb proteins are only modestly required for human cluster looping and silencing. However, computational three-dimensional Hox cluster modeling identified the insulator-binding protein CTCF as a likely candidate mediating DNA loops in all clusters. Our data suggest that Hox cluster looping may represent an evolutionarily conserved structural mechanism of transcription regulation.

Chromatin conformation signatures of cellular differentiation.

One of the major genomics challenges is to better understand how correct gene expression is orchestrated. Recent studies have shown how spatial chromatin organization is critical in the regulation of gene expression. Here, we developed a suite of computer programs to identify chromatin conformation signatures with 5C technology http://Dostielab.biochem.mcgill.ca. We identified dynamic HoxA cluster chromatin conformation signatures associated with cellular differentiation. Genome-wide chromatin conformation signature identification might uniquely identify disease-associated states and represent an entirely novel class of human disease biomarkers.