Peptide Handling by HLA-B27 Subtypes Influences Their Biological Behavior,Association with Ankylosing Spondylitis and Susceptibility to Endoplasmic Reticulum Aminopeptidase 1 (ERAP1)

HLA-B27 is strongly associated with ankylosing spondylitis (AS). We analyzed the relationship between structure, peptide specificity, folding, and stability of the seven major HLA-B27 subtypes to determine the role of their constitutive peptidomes in the pathogenicity of this molecule. Identification of large numbers of ligands allowed us to define the differences among subtype-bound peptidomes and to elucidate the peptide features associated with AS and molecular stability. The peptides identified only in AS-associated or high thermostability subtypes with identical A and B pockets were longer and had bulkier and more diverse C-terminal residues than those found only among non-AS-associated/lower-thermostability subtypes. Peptides sequenced from all AS-associated subtypes and not from non-AS-associated ones, thus strictly correlating with disease, were very rare. Residue 116 was critical in determining peptide binding, thermodynamic properties, and folding, thus emerging as a key feature that unified HLA-B27 biology. HLA-B27 ligands were better suited to TAP transport than their N-terminal precursors, and AS-associated subtype ligands were better than those from non-AS-associated subtypes, suggesting a particular capacity of AS-associated subtypes to bind epitopes directly produced in the cytosol. Peptides identified only from AS-associated/high-thermostability subtypes showed a higher frequency of ERAP1-resistant N-terminal residues than ligands found only in non-AS-associated/low-thermostability subtypes, reflecting a more pronounced effect of ERAP1 on the former group. Our results reveal the basis for the relationship between peptide specificity and other features of HLA-B27, provide a unified view of HLA-B27 biology and pathogenicity, and suggest a larger influence of ERAP1 polymorphism on AS-associated than non-AS-associated subtypes.The current ideas concerning the pathogenetic role of HLA-B27 in ankylosing spondylitis (AS) emphasize specific antigen presentation (1), misfolding (2), or immunomodulation mediated by heavy chain homodimers (3) expressed at the cell surface upon endosomal recycling (4). Recent research provided evidence that both misfolded HLA-B27 heavy chains and surface expressed B27 homodimers may activate the IL-23/IL-17 axis, a key inflammatory pathway in spondyloarthropathies, through distinct mechanisms, namely the unfolded protein response (5) and the stimulation of IL-17-producing T cells (6). In contrast, the fact that CD8+ T cells are not required for the HLA-B27-associated disease in transgenic rats (7, 8), and the failure to identify specific arthritogenic peptides, point out to a pathogenetic role of HLA-B27 based on its folding and/or non-canonical forms, rather than to an autoimmune mechanism based on molecular mimicry between foreign and self-derived peptides. Yet, on the basis of genetic and immunological studies (9, 10), an involvement of CD8+T cells in the human disease cannot be ruled out.Beyond the pathogenetic relevance of specific peptides, the constitutive HLA-B27-bound peptidome is related to the folding and stability of HLA-B27, because both features are peptide-dependent (11). This is strongly supported by the association of ERAP1, an aminopeptidase that trims peptides to their optimal size for MHC-I binding (12, 13), with ankylosing spondylitis (AS)¹ among HLA-B27-positive individuals (14), and by the demonstration that AS-associated ERAP1 polymorphism has a substantial effect on the HLA-B27 peptidome in live cells (15).Any pathogenetic mechanism must account for the differential association of HLA-B27 subtypes with AS. Whereas B*27:02, B*27:04 and B*27:05 are clearly associated with this disease, B*27:06 and B*27:09 are not (16, 17). B*27:07, a subtype present in multiple populations, is generally associated with AS, with one reported exception (18, 19). All these subtypes have the same structure in the A and B pockets of their peptide binding site, which accommodate the two N-terminal residues of their peptide ligands, but they differ in one or more positions in the F pocket, which binds the C-terminal peptide residue, as well as in other positions of the peptide binding site. In contrast, B*27:03, a subtype prevalent only in populations of Sub-Saharan African ancestry, differs from the B*27:05 prototype by a single Y59H change in the A pocket (20, 21), a difference that also sets it apart from all other subtypes (supplemental Table S1) and affects the binding preferences for N-terminal peptide residues (22–24). The nature of B*27:03 as a putative susceptibility factor for AS is unclear (19). In African populations in which this subtype is prevalent, neither this subtype nor B*27:05 are associated with this disease (25), presumably because of concurrent protective factor(s).In this study we carried out an extensive sequence analysis of HLA-B27 subtype-bound peptidomes to define their differential features as well as the extent and nature of peptide sharing among subtypes. The results revealed the basis for the intimate relationship between peptide specificity, folding, and stability of HLA-B27, provided a unified explanation on how subtype polymorphism alters the molecular biology of HLA-B27 and its association with AS, and demonstrated a differential influence of TAP and ERAP1 on AS-associated and non-AS-associated subtypes.     

Endoplasmic Reticulum Aminopeptidase 1 (ERAP1) Polymorphism Relevant to Inflammatory Disease Shapes the Peptidome of the Birdshot Chorioretinopathy-Associated HLA-A*29:02 Antigen

Birdshot chorioretinopathy is a rare ocular inflammation whose genetic association with HLA-A*29:02 is the highest between a disease and a major histocompatibility complex (MHC) molecule. It belongs to a group of MHC-I-associated inflammatory disorders, also including ankylosing spondylitis, psoriasis, and Behçet''s disease, for which endoplasmic reticulum aminopeptidases (ERAP) 1 and/or 2 have been identified as genetic risk factors. Since both enzymes are involved in the processing of MHC-I ligands, it seems reasonable that common peptide-mediated mechanisms may underlie the pathogenesis of these diseases. In this study, comparative immunopeptidomics was used to characterize >5000 A*29:02 ligands and quantify the effects of ERAP1 polymorphism and expression on the A*29:02 peptidome in human cells. The peptides predominant in an active ERAP1 context showed a higher frequency of nonamers and bulkier amino acid side chains at multiple positions, compared with the peptides predominant in a less active ERAP1 background. Thus, ERAP1 polymorphism has a large influence, shaping the A*29:02 peptidome through length-dependent and length-independent effects. These changes resulted in increased affinity and hydrophobicity of A*29:02 ligands in an active ERAP1 context. The results reveal the nature of the functional interaction between A*29:02 and ERAP1 and suggest that this enzyme may affect the susceptibility to birdshot chorioretinopathy by altering the A*29:02 peptidome. The complexity of these alterations is such that not only peptide presentation but also other potentially pathogenic features could be affected.Several major histocompatibility complex class I (MHC-I)¹ alleles are strongly associated with polygenic inflammatory diseases, including birdshot chorioretinopathy (BSCR: A*29:02), ankylosing spondylitis (AS: HLA-B*27), psoriasis (C*06:02), and Behçet''s disease (HLA-B*51). In the three latter disorders, ERAP1, an aminopeptidase of the endoplasmic reticulum performing the final trimming of MHC-I ligands (1, 2), is also a risk factor and is in epistasis with the predisposing MHC-I allele (3–5). These studies suggest common pathogenetic mechanisms involving the MHC-I bound peptidome. ERAP2, a related enzyme that acts in concert with ERAP1 (6, 7), influences the susceptibility to BSCR (8), AS (although not necessarily in epistasis with HLA-B*27) (9), Crohn′s disease (10), and preeclampsia (11–13).BSCR is a rare and severe form of bilateral posterior uveitis, showing a progressive inflammation of the choroid and retina, whose association with HLA-A*29 is the strongest for any disease and MHC. The frequency of this allele is about 7% in healthy individuals but >95% in BSCR patients (14, 15). This association specifically concerns A*29:02 and not the closely related allotype A*29:01 (8).Genetic studies on BSCR also showed a highly significant association within the LNPEP gene (rs7705093) in the 5q15 region, which includes the ERAP1 and ERAP2 genes. One single nucleotide polymorphism (SNP) in this region (rs10044354) correlated with ERAP2 expression. This was confirmed at the protein level, leading to the conclusion that ERAP2 expression predisposes to BSCR. Yet, an involvement of functional ERAP1 polymorphisms, not determining protein expression, was not excluded. These polymorphisms have a large influence on the HLA-B*27 peptidome (16, 17). In contrast, the effects of ERAP2 on MHC-I peptidomes are poorly understood and are probably dependent on the particular ERAP1 context since ERAP2 cooperates with ERAP1 in peptide processing. Thus, the present study was conducted to characterize A*29:02-bound peptidomes in various ERAP1 backgrounds and to determine the influence of ERAP1 polymorphism on the amounts and features of A*29:02 ligands in human cells.     

Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Structural Basis for T Cell Alloreactivity among Three HLA-B14 and HLA-B27 Antigens

The existence of cytotoxic T cells (CTL) cross-reacting with the human major histocompatibility antigens HLA-B14 and HLA-B27 suggests that their alloreactivity could be due to presentation of shared peptides in similar binding modes by these molecules. We therefore determined the crystal structures of the subtypes HLA-B*1402, HLA-B*2705, and HLA-B*2709 in complex with a proven self-ligand, pCatA (peptide with the sequence IRAAPPPLF derived from cathepsin A (residues 2–10)), and of HLA-B*1402 in complex with a viral peptide, pLMP2 (RRRWRRLTV, derived from latent membrane protein 2 (residues 236–244) of Epstein-Barr virus). Despite the exchange of 18 residues within the binding grooves of HLA-B*1402 and HLA-B*2705 or HLA-B*2709, the pCatA peptide is presented in nearly identical conformations. However, pLMP2 is displayed by HLA-B*1402 in a conformation distinct from those previously found in the two HLA-B27 subtypes. In addition, the complexes of HLA-B*1402 with the two peptides reveal a nonstandard, tetragonal mode of the peptide N terminus anchoring in the binding groove because of the exchange of the common Tyr-171 by His-171 of the HLA-B*1402 heavy chain. This exchange appears also responsible for reduced stability of HLA-B14-peptide complexes in vivo and slow assembly in vitro. The studies with the pCatA peptide uncover that CTL cross-reactive between HLA-B14 and HLA-B27 might primarily recognize the common structural features of the bound peptide, thus neglecting amino acid replacements within the rim of the binding grooves. In contrast, structural alterations between the three complexes with the pLMP2 peptide indicate how heavy chain polymorphisms can influence peptide display and prevent CTL cross-reactivity between HLA-B14 and HLA-B27 antigens.T cells possessing the ability to recognize major histocompatibility complex (MHC)² molecules from another individual of the same species, also termed alloreactive T cells, may constitute up to 10% of the T cell pool of an individual, and their precursor frequency can be 100–1,000-fold higher than that of self-restricted T cells directed against a foreign peptide (1, 2). The ability of alloreactive T cells to cross-react with nonself-MHC molecules is a major obstacle preventing successful organ transplantations (3–5). Two mechanisms, direct or indirect allorecognition, can be responsible for the rejection of a transplant by alloreactive T cells (6). In the first case, donor cells expressing MHC molecules are directly recognized by host T cells (7), whereas indirect allorecognition involves the presentation of peptides derived from donor proteins by MHC molecules of the host, followed by the detection of the complexes by the host T cells (8). However, although alloreactive T cells are very common and of great clinical importance, neither the primary basis for their existence nor the reasons underlying their cross-reactivity are sufficiently understood to draw general conclusions (9–11). Only very few studies have addressed the structural basis for the recognition of distinct MHC antigens by cross-reactive T cells (12–18). One of the most important questions regards the individual contribution of the bound peptide and binding groove residues of the heavy chain (HC) of MHC class I antigens to the interaction with T cell receptors (TCR).Here we analyze an HLA-B14 subtype, HLA-B*1402 (named B*1402), as well as two HLA-B27 subtypes, HLA-B*2705 and HLA-B*2709 (named B*2705 and B*2709), to shed light on the structural basis of peptide presentation and T cell alloreactivity among these HLA-B molecules. The amino acid sequences of B*1402 and B*2705 HC differ from each other at 18 positions, all of which are part of the peptide-binding groove (Fig. 1). These amino acid exchanges result in different repertoires of bound peptides; B*1402 and B*2705 share only about 4% of their peptides (19), whereas this value rises to 88% for the B*2705 and B*2709 subtypes (20), which are distinguished only by a single residue at the floor of the binding groove (B*2705, Asp-116; B*2709, His-116). The structural similarities between the two HLA-B27 subtypes (21–27) permit extensive cross-reactivity (up to 90%) of cytotoxic T cells (CTL) (28), whereas CTL alloreactivity between B*1402 and B*2705 is drastically reduced (to about 3%) (19), in line with the very limited overlap of their peptide repertoires.Open in a separate windowFIGURE 1.Amino acid sequence differences among B*1402 and B*2705 HC. The 18 residues distinguishing the two subtypes are all located in or in the immediate vicinity of the peptide-binding groove. B*2705 differs from B*2709 only by a D116H exchange (not shown). The residues are indicated by spheres with volumes roughly proportional to the volumes of the respective amino acid side chain in solution (77). The spheres are colored according to the biochemical properties of the respective amino acids, as indicated at the bottom of the image.The HLA-B14 and HLA-B27 subtypes are distinguished from most other HLA class I molecules in their requirement for an arginine at anchor position 2 of the bound peptide (p2) (20, 29, 30). This preference is nearly absolute in B*2705 and B*2709 (31), whereas B*1402 tolerates also glutamine, glutamate, and proline as p2 anchors (19, 29). Statistically significant differences between B*1402 and B*2705 are also found at several other peptide positions (19). Previous structural and cellular studies of the HLA-B27 subtypes have suggested that molecular mimicry between the viral peptide pLMP2 (RRRWRRLTV, derived from Epstein-Barr virus latent membrane protein 2, residues 236–244) and the self-peptide pVIPR (RRKWRRWHL, derived from vasoactive intestinal peptide type 1 receptor, residues 400–408), when bound to B*2705, serves as an example of how a cellular immune response could be triggered that might contribute to the onset of ankylosing spondylitis (AS) through an autoimmune mechanism (22, 24). CTL that recognize the B*2705 and the B*2709 subtypes in complex with the self-peptide pVIPR (22) exemplify alloreactivity in this system, although the D116H micropolymorphism is deeply buried and not directly accessible to a TCR.Alloreactive T cells are known to recognize a very diverse array of alloantigen-bound peptides (32, 33), so that virtually each T cell clone can be assumed to be specific for a distinct peptide. For this reason, the substantial correlation found in previous studies between peptide and the alloreactive T cell epitope sharing among HLA-B27 (reviewed in Ref. 34) or HLA-B14 subtypes (only 28.4% partial or full cross-reactivity, similar to peptide overlapping between the subtypes B*1402 and B*1403, see Ref. 19) supports a prominent role of peptides in determining alloreactive T cell cross-reaction, and it suggests that many shared ligands adopt antigenically similar conformations when bound to distinct HLA-B molecules. On the other hand, the results reported by Merino et al. (19) also demonstrate that the few CTL that cross-react with B*1402 and B*2705 did not exhibit cross-reactivity with B*1403, which is distinguished from B*1402 only by a single amino acid exchange in the α2-helix. Furthermore, they show that alloreactive CTL from various donors directed against B*2705 did not lyse cells expressing either B*1402 or B*1403, although the number of CTL tested might not have been high enough to detect a presumably low degree of cross-reactivity. Without structural data from HLA-B14 subtypes, however, these results are difficult to interpret.The pCatA peptide (IRAAPPPLF, derived from the signal sequence of cathepsin A, residues 2–10) is among the very few known common ligands of B*1402, B*2705 (19), and B*2709³ and can thus serve to study how a very different (B*1402) and two very similar subtypes (B*2705 and B*2709) handle a common ligand. On the other hand, the pLMP2 peptide is a proven natural ligand only of B*2705, whose possible presentation in vivo by B*2709 and HLA-B14 is not yet known, although this peptide can be complexed in vitro with B*2709 (24) and also with B*1402 (35). From previous crystallographic studies, it was known that pLMP2 is presented by the two HLA-B27 antigens in very different conformations (24). We expected that the pronounced sequence differences between B*1402 and the HLA-B27 alloantigens (Fig. 1) might even enhance the conformational dissimilarities that are observed when two very closely related subtypes such as B*2705 and B*2709 are compared. Discrepancies in peptide display could reasonably be expected to prevent CTL cross-reaction, so that pLMP2 might be considered as a representative of the vast majority of HLA-B14- and HLA-B27-presented ligands that must be responsible for the low degree of CTL cross-reactivity between these alloantigens. Despite these presumed differences between pCatA and pLMP2, both peptides may be seen as examples of ligands that could principally allow direct allorecognition.Here we report the crystal structures of B*1402·pCatA, B*2705·pCatA, B*2709·pCatA, and B*1402·pLMP2, and we compare them with each other and with the previously reported structures of B*2705·pLMP2 and B*2709·pLMP2 (24).     

Endogenous Processing and Presentation of T-cell Epitopes from Chlamydia trachomatis with Relevance in HLA-B27-associated Reactive Arthritis

Chlamydia trachomatis triggers reactive arthritis, a spondyloarthropathy linked to the human major histocompatibility complex molecule HLA-B27, through an unknown mechanism that might involve molecular mimicry between chlamydial and self-derived HLA-B27 ligands. Chlamydia-specific CD8⁺ T-cells are found in reactive arthritis patients, but the immunogenic epitopes are unknown. A previous screening of the chlamydial genome for putative HLA-B27 ligands predicted multiple peptides that were recognized in vitro by CD8⁺ T-lymphocytes from patients. Here stable transfectants expressing bacterial fusion proteins in human cells were generated to investigate the endogenous processing and presentation by HLA-B27 of two such epitopes through comparative immunoproteomics of HLA-B27-bound peptide repertoires. A predicted T-cell epitope, from the CT610 gene product, was presented by HLA-B27. This is, to our knowledge, the first endogenously processed epitope involved in HLA-B27-restricted responses against C. trachomatis in reactive arthritis. A second predicted epitope, from the CT634 gene product, was not detected. Instead a non-predicted nonamer from the same protein was identified. Both bacterial peptides showed very high homology with human sequences containing the HLA-B27 binding motif. Thus, expression and intracellular processing of chlamydial proteins into human cells allowed us to identify two bacterial HLA-B27 ligands, including the first endogenous T-cell epitope from C. trachomatis involved in spondyloarthropathy. That human proteins contain sequences mimicking chlamydial T-cell epitopes suggests a basis for an autoimmune component of Chlamydia-induced HLA-B27-associated disease.Chlamydia trachomatis is an obligate intracellular parasite that infects the urogenital epithelium. It is a very common pathogen and one frequently inducing reactive arthritis (ReA)¹ (1). Multiple strategies, including down-regulation of major histocompatibility complex (MHC) class I and class II expression (2–4) and persistence, have been developed by the bacteria to evade the immune system. Yet both CD4⁺ and CD8⁺ T-cell responses are activated upon infection (5). In particular, HLA-B27-restricted CD8⁺ T-lymphocytes are found in patients with Chlamydia-induced ReA (6, 7). The role of these cells in the pathogenesis and evolution of ReA to chronic disease is probably mediated by IFN-γ. Secretion of this cytokine is critical for the protective immunity against Chlamydia (8) because it inhibits the bacterial growth (9). However, this is often insufficient to promote complete clearance of C. trachomatis, and actually IFN-γ-induced depletion of the tryptophan pool induces the differentiation of the metabolically active reticular bodies to persistent forms (10), which sustain chronic infection and ReA. The high prevalence of HLA-B27 among patients with Chlamydia-induced ReA (11), especially in its chronic form, suggests a pathogenetic mechanism based on interactive effects of the bacteria and HLA-B27 that seems unrelated to the capacity of C. trachomatis to infect or replicate into HLA-B27-positive cells (12). One such mechanism could be T-cell-mediated autoimmunity elicited by molecular/antigenic mimicry between chlamydial and self-derived HLA-B27 ligands. Antigenic mimicry between chlamydial and homologous α-myosin-derived peptides is crucial to inducing autoimmune myocarditis in mice (13). Breakdown of cytotoxic T-lymphocyte (CTL) tolerance to HLA-B27 was observed in transgenic rats upon exposure to C. trachomatis (14). Cross-reactivity between HLA-B27-restricted Chlamydia-specific CTL with self-derived HLA-B27 epitopes has not been reported. However, a biochemical basis for it was provided by the finding of an endogenously processed and presented peptide from the DNA primase of C. trachomatis with high homology to a self-derived HLA-B27 ligand (15, 16).Because of the likely involvement of HLA-B27 in the pathogenesis of chronically evolving ReA, the role of CD8⁺ T-cell responses in the protective immunity against C. trachomatis and the presence of HLA-B27-restricted T-cells in patients with Chlamydia-induced ReA, the identification of relevant chlamydial epitopes becomes crucial to establish the pathogenetic mechanism of this disease. Unfortunately a direct analysis of chlamydial HLA-B27 ligands expressed on infected cells is exceedingly difficult because of their extremely low amounts, which challenge even the most sensitive techniques of MS. In the case of Chlamydia, the situation is further complicated by the down-regulation of MHC class I expression shortly after infection (3, 4). To our knowledge, only one MHC class I ligand was recently identified, in the mouse system, from Chlamydia muridarum-infected cells using state-of-the-art MS techniques (17). Due in part to this difficulty, alternative approaches, such as expression cloning and synthetic peptide epitope mapping (18, 19) or MHC class I tetramer arrays (20), have been used to identify MHC class I-restricted chlamydial T-cell epitopes in mice. In a previous study (6) predictive algorithms were used to screen the whole genome of C. trachomatis for nonamer peptide sequences containing the HLA-B*2705 binding motif and a high probability of being generated by proteasomal cleavage. This led to identifying multiple sequences that, when used as synthetic peptides in vitro, stimulated CD8⁺ T-cells from patients with Chlamydia-induced ReA. Such cells could also be detected in the synovial fluid of these patients using HLA-B27 tetramers complexed to some of these peptides (7).Although these strategies identify chlamydial sequences that are recognized by CD8⁺ T-cells they do not prove that these peptides are the endogenously processed epitopes that activated the natural T-cell responses to the bacteria in vivo. Because of the intrinsic cross-reactivity of T-cells (21, 22), it is conceivable that synthetic peptides recognized in vitro may be different from the natural epitopes generated by endogenous processing of the chlamydial proteins that elicit the HLA-B27-restricted T-cell responses in ReA patients. To investigate this issue we focused on two predicted epitopes (6). Stable transfectants expressing the corresponding chlamydial proteins fused to green fluorescent protein (GFP) were generated in a B*2705-positive cell line. The endogenous processing and presentation of the predicted epitopes or other peptides from the same bacterial protein were analyzed by comparative immunoproteomics analysis of the B*2705-bound peptide repertoires from transfected and untransfected cells and sequencing of peptides differentially presented on the bacterial protein transfectant.     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Phosphorylation Dependence of Hsp27 Multimeric Size and Molecular Chaperone Function

The molecular chaperone Hsp27 exists as a distribution of large oligomers that are disassembled by phosphorylation at Ser-15, -78, and -82. It is controversial whether the unphosphorylated Hsp27 or the widely used triple Ser-to-Asp phospho-mimic mutant is the more active molecular chaperone in vitro. This question was investigated here by correlating chaperone activity, as measured by the aggregation of reduced insulin or α-lactalbumin, with Hsp27 self-association as monitored by analytical ultracentrifugation. Furthermore, because the phospho-mimic is generally assumed to reproduce the phosphorylated molecule, the size and chaperone activity of phosphorylated Hsp27 were compared with that of the phospho-mimic. Hsp27 was triply phosphorylated by MAPKAP-2 kinase, and phosphorylation was tracked by urea-PAGE. An increasing degree of suppression of insulin or α-lactalbumin aggregation correlated with a decreasing Hsp27 self-association, which was the least for phosphorylated Hsp27 followed by the mimic followed by the unphosphorylated protein. It was also found that Hsp27 added to pre-aggregated insulin did not reverse aggregation but did inhibit these aggregates from assembling into even larger aggregates. This chaperone activity appears to be independent of Hsp27 phosphorylation. In conclusion, the most active chaperone of insulin and α-lactalbumin was the Hsp27 (elongated) dimer, the smallest Hsp27 subunit observed under physiological conditions. Next, the Hsp27 phospho-mimic is only a partial mimic of phosphorylated Hsp27, both in self-association and in chaperone function. Finally, the efficient inhibition of insulin aggregation by Hsp27 dimer led to the proposal of two models for this chaperone activity.Oligomeric heat shock protein 27 (Hsp27)² is a ubiquitous mammalian protein with a variety of functions in health and disease (1–8). These functions include ATP-independent chaperone activity in response to environmental stress, e.g. heat shock and oxidative stress, control of apoptosis, and regulation of actin cytoskeleton dynamics. Hsp27 is a member of the α-crystallin small heat shock protein family of which αB-crystallin is the archetype. These proteins are characterized by an α-crystallin domain of 80–90 residues consisting of roughly eight β-strands that form an intermolecular β-sheet interaction interface within a dimer, the basic building subunit of the oligomer (2, 4, 9–11).Hsp27 is in equilibrium between high molecular weight oligomers and much lower molecular weight multimers. It has been reported that unphosphorylated Hsp27 includes predominantly a distribution of high molecular species ranging in size from 12-mer to 35-mer (12–19). Phosphorylation of Hsp27 at serines 15, 78, and 82 by the p38-activated MAPKAP-2 kinase (20–22) or the use of the triple Ser-to-Asp phospho-mimic results in a major shift in the equilibrium toward much smaller multimers (23) and in an alteration of its function (1, 3, 6, 7, 24, 25). The size distribution of the smaller species has been reported to be between monomer and tetramer (12–16, 18, 19).Small heat shock proteins, including Hsp27, behave as ATP-independent molecular chaperones during cellular heat shock. They bind partially unfolded proteins and prevent their aggregation until the proteins can be refolded by larger ATP-dependent chaperones or are digested (7, 8, 26). This function includes the up-regulation and/or phosphorylation of Hsp27.It is not entirely clear what the role of Hsp27 size and phosphorylation state plays in its heat shock function because there are conflicting results in the literature. Some in vitro studies concluded that the unphosphorylated oligomeric Hsp27 (or the murine isoform Hsp25) protects proteins against aggregation better than does the phosphorylation mimic (13, 19, 27), whereas others found no difference (16, 28, 29), and still other studies found that the mimic protects better than does the unphosphorylated wild type (27, 30, 31). In-cell studies found that phosphorylation of Hsp27 was essential for thermo-protection of actin filaments (32), and the Hsp27 phosphorylation mimic decreased inclusion body formation better than did unphosphorylated Hsp27 (33). This study was undertaken to investigate the molecular chaperone function of Hsp27 by correlating chaperone activity with Hsp27 size and by comparing fully phosphorylated Hsp27 with its phospho-mimic.     

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

A variety of high-throughput methods have made it possible to generate detailed temporal expression data for a single gene or large numbers of genes. Common methods for analysis of these large data sets can be problematic. One challenge is the comparison of temporal expression data obtained from different growth conditions where the patterns of expression may be shifted in time. We propose the use of wavelet analysis to transform the data obtained under different growth conditions to permit comparison of expression patterns from experiments that have time shifts or delays. We demonstrate this approach using detailed temporal data for a single bacterial gene obtained under 72 different growth conditions. This general strategy can be applied in the analysis of data sets of thousands of genes under different conditions.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Multipattern Consensus Regions in Multiple Aligned Protein Sequences and Their Segmentation

Decomposing a biological sequence into its functional regions is an important prerequisite to understand the molecule. Using the multiple alignments of the sequences, we evaluate a segmentation based on the type of statistical variation pattern from each of the aligned sites. To describe such a more general pattern, we introduce multipattern consensus regions as segmented regions based on conserved as well as interdependent patterns. Thus the proposed consensus region considers patterns that are statistically significant and extends a local neighborhood. To show its relevance in protein sequence analysis, a cancer suppressor gene called p53 is examined. The results show significant associations between the detected regions and tendency of mutations, location on the 3D structure, and cancer hereditable factors that can be inferred from human twin studies.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]     

Printor, a Novel TorsinA-interacting Protein Implicated in Dystonia Pathogenesis

Early onset generalized dystonia (DYT1) is an autosomal dominant neurological disorder caused by deletion of a single glutamate residue (torsinA ΔE) in the C-terminal region of the AAA⁺ (ATPases associated with a variety of cellular activities) protein torsinA. The pathogenic mechanism by which torsinA ΔE mutation leads to dystonia remains unknown. Here we report the identification and characterization of a 628-amino acid novel protein, printor, that interacts with torsinA. Printor co-distributes with torsinA in multiple brain regions and co-localizes with torsinA in the endoplasmic reticulum. Interestingly, printor selectively binds to the ATP-free form but not to the ATP-bound form of torsinA, supporting a role for printor as a cofactor rather than a substrate of torsinA. The interaction of printor with torsinA is completely abolished by the dystonia-associated torsinA ΔE mutation. Our findings suggest that printor is a new component of the DYT1 pathogenic pathway and provide a potential molecular target for therapeutic intervention in dystonia.Early onset generalized torsion dystonia (DYT1) is the most common and severe form of hereditary dystonia, a movement disorder characterized by involuntary movements and sustained muscle spasms (1). This autosomal dominant disease has childhood onset and its dystonic symptoms are thought to result from neuronal dysfunction rather than neurodegeneration (2, 3). Most DYT1 cases are caused by deletion of a single glutamate residue at positions 302 or 303 (torsinA ΔE) of the 332-amino acid protein torsinA (4). In addition, a different torsinA mutation that deletes amino acids Phe³²³–Tyr³²⁸ (torsinA Δ323–328) was identified in a single family with dystonia (5), although the pathogenic significance of this torsinA mutation is unclear because these patients contain a concomitant mutation in another dystonia-related protein, ϵ-sarcoglycan (6). Recently, genetic association studies have implicated polymorphisms in the torsinA gene as a genetic risk factor in the development of adult-onset idiopathic dystonia (7, 8).TorsinA contains an N-terminal endoplasmic reticulum (ER)³ signal sequence and a 20-amino acid hydrophobic region followed by a conserved AAA⁺ (ATPases associated with a variety of cellular activities) domain (9, 10). Because members of the AAA⁺ family are known to facilitate conformational changes in target proteins (11, 12), it has been proposed that torsinA may function as a molecular chaperone (13, 14). TorsinA is widely expressed in brain and multiple other tissues (15) and is primarily associated with the ER and nuclear envelope (NE) compartments in cells (16–20). TorsinA is believed to mainly reside in the lumen of the ER and NE (17–19) and has been shown to bind lamina-associated polypeptide 1 (LAP1) (21), lumenal domain-like LAP1 (LULL1) (21), and nesprins (22). In addition, recent evidence indicates that a significant pool of torsinA exhibits a topology in which the AAA⁺ domain faces the cytoplasm (20). In support of this topology, torsinA is found in the cytoplasm, neuronal processes, and synaptic terminals (2, 3, 15, 23–26) and has been shown to bind cytosolic proteins snapin (27) and kinesin light chain 1 (20). TorsinA has been proposed to play a role in several cellular processes, including dopaminergic neurotransmission (28–31), NE organization and dynamics (17, 22, 32), and protein trafficking (27, 33). However, the precise biological function of torsinA and its regulation remain unknown.To gain insights into torsinA function, we performed yeast two-hybrid screens to search for torsinA-interacting proteins in the brain. We report here the isolation and characterization of a novel protein named printor (protein interactor of torsinA) that interacts selectively with wild-type (WT) torsinA but not the dystonia-associated torsinA ΔE mutant. Our data suggest that printor may serve as a cofactor of torsinA and provide a new molecular target for understanding and treating dystonia.     

Sip1, a Novel RS Domain-Containing Protein Essential for Pre-mRNA Splicing

Previous studies have shown that protein-protein interactions among splicing factors may play an important role in pre-mRNA splicing. We report here identification and functional characterization of a new splicing factor, Sip1 (SC35-interacting protein 1). Sip1 was initially identified by virtue of its interaction with SC35, a splicing factor of the SR family. Sip1 interacts with not only several SR proteins but also with U1-70K and U2AF65, proteins associated with 5′ and 3′ splice sites, respectively. The predicted Sip1 sequence contains an arginine-serine-rich (RS) domain but does not have any known RNA-binding motifs, indicating that it is not a member of the SR family. Sip1 also contains a region with weak sequence similarity to the Drosophila splicing regulator suppressor of white apricot (SWAP). An essential role for Sip1 in pre-mRNA splicing was suggested by the observation that anti-Sip1 antibodies depleted splicing activity from HeLa nuclear extract. Purified recombinant Sip1 protein, but not other RS domain-containing proteins such as SC35, ASF/SF2, and U2AF65, restored the splicing activity of the Sip1-immunodepleted extract. Addition of U2AF65 protein further enhanced the splicing reconstitution by the Sip1 protein. Deficiency in the formation of both A and B splicing complexes in the Sip1-depleted nuclear extract indicates an important role of Sip1 in spliceosome assembly. Together, these results demonstrate that Sip1 is a novel RS domain-containing protein required for pre-mRNA splicing and that the functional role of Sip1 in splicing is distinct from those of known RS domain-containing splicing factors.Pre-mRNA splicing takes place in spliceosomes, the large RNA-protein complexes containing pre-mRNA, U1, U2, U4/6, and U5 small nuclear ribonucleoprotein particles (snRNPs), and a large number of accessory protein factors (for reviews, see references 21, 22, 37, 44, and 48). It is increasingly clear that the protein factors are important for pre-mRNA splicing and that studies of these factors are essential for further understanding of molecular mechanisms of pre-mRNA splicing.Most mammalian splicing factors have been identified by biochemical fractionation and purification (3, 15, 19, 31–36, 45, 69–71, 73), by using antibodies recognizing splicing factors (8, 9, 16, 17, 61, 66, 67, 74), and by sequence homology (25, 52, 74).Splicing factors containing arginine-serine-rich (RS) domains have emerged as important players in pre-mRNA splicing. These include members of the SR family, both subunits of U2 auxiliary factor (U2AF), and the U1 snRNP protein U1-70K (for reviews, see references 18, 41, and 59). Drosophila alternative splicing regulators transformer (Tra), transformer 2 (Tra2), and suppressor of white apricot (SWAP) also contain RS domains (20, 40, 42). RS domains in these proteins play important roles in pre-mRNA splicing (7, 71, 75), in nuclear localization of these splicing proteins (23, 40), and in protein-RNA interactions (56, 60, 64). Previous studies by us and others have demonstrated that one mechanism whereby SR proteins function in splicing is to mediate specific protein-protein interactions among spliceosomal components and between general splicing factors and alternative splicing regulators (1, 1a, 6, 10, 27, 63, 74, 77). Such protein-protein interactions may play critical roles in splice site recognition and association (for reviews, see references 4, 18, 37, 41, 47 and 59). Specific interactions among the splicing factors also suggest that it is possible to identify new splicing factors by their interactions with known splicing factors.Here we report identification of a new splicing factor, Sip1, by its interaction with the essential splicing factor SC35. The predicted Sip1 protein sequence contains an RS domain and a region with sequence similarity to the Drosophila splicing regulator, SWAP. We have expressed and purified recombinant Sip1 protein and raised polyclonal antibodies against the recombinant Sip1 protein. The anti-Sip1 antibodies specifically recognize a protein migrating at a molecular mass of approximately 210 kDa in HeLa nuclear extract. The anti-Sip1 antibodies sufficiently deplete Sip1 protein from the nuclear extract, and the Sip1-depleted extract is inactive in pre-mRNA splicing. Addition of recombinant Sip1 protein can partially restore splicing activity to the Sip1-depleted nuclear extract, indicating an essential role of Sip1 in pre-mRNA splicing. Other RS domain-containing proteins, including SC35, ASF/SF2, and U2AF65, cannot substitute for Sip1 in reconstituting splicing activity of the Sip1-depleted nuclear extract. However, addition of U2AF65 further increases splicing activity of Sip1-reconstituted nuclear extract, suggesting that there may be a functional interaction between Sip1 and U2AF65 in nuclear extract.