Sampling From the Proteome to the Human Leukocyte Antigen-DR (HLA-DR) Ligandome Proceeds Via High Specificity

Comprehensive analysis of the complex nature of the Human Leukocyte Antigen (HLA) class II ligandome is of utmost importance to understand the basis for CD4⁺ T cell mediated immunity and tolerance. Here, we implemented important improvements in the analysis of the repertoire of HLA-DR-presented peptides, using hybrid mass spectrometry-based peptide fragmentation techniques on a ligandome sample isolated from matured human monocyte-derived dendritic cells (DC). The reported data set constitutes nearly 14 thousand unique high-confident peptides, i.e. the largest single inventory of human DC derived HLA-DR ligands to date. From a technical viewpoint the most prominent finding is that no single peptide fragmentation technique could elucidate the majority of HLA-DR ligands, because of the wide range of physical chemical properties displayed by the HLA-DR ligandome. Our in-depth profiling allowed us to reveal a strikingly poor correlation between the source proteins identified in the HLA class II ligandome and the DC cellular proteome. Important selective sieving from the sampled proteome to the ligandome was evidenced by specificity in the sequences of the core regions both at their N- and C- termini, hence not only reflecting binding motifs but also dominant protease activity associated to the endolysosomal compartments. Moreover, we demonstrate that the HLA-DR ligandome reflects a surface representation of cell-compartments specific for biological events linked to the maturation of monocytes into antigen presenting cells. Our results present new perspectives into the complex nature of the HLA class II system and will aid future immunological studies in characterizing the full breadth of potential CD4⁺ T cell epitopes relevant in health and disease.Human Leukocyte Antigen (HLA)¹ class II molecules on professional antigen presenting cells such as dendritic cells (DC) expose peptide fragments derived from exogenous and endogenous proteins to be screened by CD4⁺ T cells (1, 2). The activation and recruitment of CD4⁺ T cells recognizing disease-related peptide antigens is critical for the development of efficient antipathogen or antitumor immunity. Furthermore, the presentation of self-peptides and their interaction with CD4⁺ T cells is essential to maintain immunological tolerance and homeostasis (3). Knowledge of the nature of HLA class II-presented peptides on DC is of great importance to understand the rules of antigen processing and peptide binding motifs (4), whereas the identity of disease-related antigens may provide new knowledge on immunogenicity and leads for the development of vaccines and immunotherapy (5, 6).Mass spectrometry (MS) has proven effective for the analysis HLA class II-presented peptides (4, 7, 8). MS-based ligandome studies have demonstrated that HLA class II molecules predominantly present peptides derived from exogenous proteins that entered the cells by endocytosis and endogenous proteins that are associated with the endo-lysosomal compartments (4). Yet proteins residing in the cytosol, nucleus or mitochondria can also be presented by HLA class II molecules, primarily through autophagy (9–11). Multiple studies have mapped the HLA class II ligandome of antigen presenting cells in the context of infectious pathogens (12), autoimmune diseases (13–17) or cancer (14, 18, 19), or those that are essential for self-tolerance in the human thymus (3, 20). Notwithstanding these efforts, and certainly not in line with the extensive knowledge on the HLA class I ligandome (21), the nature of the HLA class II-presented peptide repertoire and particular its relationship to the cellular source proteome remains poorly understood.To advance our knowledge on the HLA-DR ligandome on activated DC without having to deal with limitations in cell yield from peripheral human blood (12, 21, 22) or tissue isolates (3), we explored the use of MUTZ-3 cells. This cell line has been used as a model of human monocyte-derived DCs. MUTZ-3 cells can be matured to act as antigen presenting cells and express then high levels of HLA class II molecules, and can be propagated in vitro to large cell densities (23–25). We also evaluated the performance of complementary and hybrid MS fragmentation techniques electron-transfer dissociation (ETD), electron-transfer/higher-energy collision dissociation (EThcD) (26), and higher-energy collision dissociation (HCD) to sequence and identify the HLA class II ligandome. Together this workflow allowed for the identification of an unprecedented large set of about 14 thousand unique peptide sequences presented by DC derived HLA-DR molecules, providing an in-depth view of the complexity of the HLA class II ligandome, revealing underlying features of antigen processing and surface-presentation to CD4⁺ T cells.     

Epstein-Barr Virus Lacking Glycoprotein gp42 Can Bind to B Cells but Is Not Able To Infect

The Epstein-Barr virus gH-gL complex includes a third glycoprotein, gp42, which is the product of the BZLF2 open reading frame (ORF). gp42 has been implicated as critical to infection of the B lymphocyte by virtue of its interaction with HLA class II on the B-cell surface. A neutralizing antibody that reacts with gp42 inhibits virus-cell fusion and blocks binding of gp42 to HLA class II; antibody to HLA class II can inhibit infection, and B cells that lack HLA class II can only be infected if HLA class II expression is restored. To confirm whether gp42 is an essential component of the virion, we derived a recombinant virus with a selectable marker inserted into the BZLF2 ORF to interrupt expression of the protein. A complex of gH and gL was expressed by the recombinant virus in the absence of gp42. Recombinant virus egressed from the cell normally and could bind to receptor-positive cells. It had, however, lost the ability to infect or transform B lymphocytes. Treatment with polyethylene glycol restored the infectivity of recombinant virus, confirming that gp42 is essential for penetration of the B-cell membrane.Entry of enveloped viruses into mammalian cells requires that the virion envelope fuse with the cell membrane after attachment to the cell surface. Herpesviruses require the functions of multiple protein species to mediate this event, and in keeping with the common origin and diverse habitats of these viruses, some of the proteins involved in penetration appear to be conserved throughout the family and some appear to be restricted to individual members or more closely related members with similar tropism. The two glycoproteins gH and gL fall into the first category of conserved proteins. Glycoprotein gH has been implicated as a major player in virus-cell fusion in many herpesviruses (8, 10, 11, 22, 28, 32, 34), and gL is an essential partner which is required for folding and transport of gH out of the endoplasmic reticulum (6, 19, 21, 27, 28, 35, 38, 45). The gH and gL homologs of Epstein-Barr virus (EBV) are gp85, the product of the BXLF2 open reading frame (ORF) (13, 31), and gp25, the product of the BKRF2 ORF (45), and these homologs appear to behave much as their counterparts in other herpesviruses do (45). However, a third glycoprotein, gp42, associates with the EBV gH-gL complex and falls into the second category of proteins, those with a more restricted distribution.Glycoprotein gp42 is the product of the BZLF2 ORF (26), and although there may be a functionally similar protein in cytomegalovirus (18, 24), it is not predicted to have a homolog in other human herpesviruses. It does, however, have a homolog in ORF51 of equine herpes virus 2 (43). Both EBV and equine herpes virus 2 infect B lymphocytes (1), and several lines of evidence suggest that, at least in the case of EBV, gp42 is critical to the infection of this cell type. A monoclonal antibody (MAb) called F-2-1 that reacts with gp42 has no affect on EBV attachment to its receptor, complement receptor type 2 (CR2) (CD21), but inhibits fusion of the virus with the B-cell membrane and neutralizes infection (29). Glycoprotein gp42 interacts with the β₁ domain of the HLA class II protein HLA-DR (39), and MAb F-2-1 interferes with this interaction (25). Like F-2-1, a MAb to HLA-DR or a soluble form of gp42 can block B-cell transformation, and B-cell lines which lack expression of HLA class II are not susceptible to superinfection with EBV unless expression of HLA class II is restored (25). Collectively these observations suggest that gp42, probably by virtue of its interaction with HLA class II, is essential to infection of the B lymphocyte. To answer directly the question of whether gp42 is an indispensable glycoprotein, we derived a virus that could be definitively shown to lack expression of the molecule and examined its ability to infect normal resting B lymphocytes. We report here that virus with expression of gp42 blocked can exit cells normally and can bind to receptor-positive target cells. However, it is unable to penetrate into cells and initiate infection.     

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.     

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]     

The Effect of Proteasome Inhibition on the Generation of the Human Leukocyte Antigen (HLA) Peptidome

The Major histocompatibility complex (MHC) class I peptidome is thought to be generated mostly through proteasomal degradation of cellular proteins, a notion that is based on the alterations in presentation of selected peptides following proteasome inhibition. We evaluated the effects of proteasome inhibitors, epoxomicin and bortezomib, on human cultured cancer cells. Because the inhibitors did not reduce the level of presentation of the cell surface human leukocyte antigen (HLA) molecules, we followed their effects on the rates of synthesis of both HLA peptidome and proteome of the cells, using dynamic stable isotope labeling in tissue culture (dynamic-SILAC). The inhibitors reduced the rates of synthesis of most cellular proteins and HLA peptides, yet the synthesis rates of some of the proteins and HLA peptides was not decreased by the inhibitors and of some even increased. Therefore, we concluded that the inhibitors affected the production of the HLA peptidome in a complex manner, including modulation of the synthesis rates of the source proteins of the HLA peptides, in addition to their effect on their degradation. The collected data may suggest that the current reliance on proteasome inhibition may overestimate the centrality of the proteasome in the generation of the MHC peptidome. It is therefore suggested that the relative contribution of the proteasomal and nonproteasomal pathways to the production of the MHC peptidome should be revaluated in accordance with the inhibitors effects on the synthesis rates of the source proteins of the MHC peptides.The repertoires and levels of peptides, presented by the major histocompatibility complex (MHC)¹ class I molecules at the cells'' surface, are modulated by multiple factors. These include the rates of synthesis and degradation of their source proteins, the transport efficacy of the peptides through the transporter associated with antigen processing (TAP) into the endoplasmic reticulum (ER), their subsequent processing and loading onto the MHC molecules within the ER, and the rates of transport of the MHC molecules with their peptide cargo to the cell surface. The off-rates of the presented peptides, the residence time of the MHC complexes at the cell surface, and their retrograde transport back into the cytoplasm, influence, as well, the presented peptidomes (reviewed in (1)). Even though significant portions of the MHC class I peptidomes are thought to be derived from newly synthesized proteins, including misfolded proteins, defective ribosome products (DRiPs), and short lived proteins (SLiPs), most of the MHC peptidome is assumed to originate from long-lived proteins, which completed their functional cellular roles or became defective (retirees), (reviewed in (2)).The main protease, supplying the MHC peptidome production pipeline, is thought to be the proteasome (3). It is also responsible for generation of the final carboxyl termini of the MHC peptides (4), (reviewed in (5)). The final trimming of the n-termini of the peptides, within the endoplasmic reticulum (ER), is thought to be performed by amino peptidases, such as ERAP1/ERAAP, which fit the peptides into their binding groove on the MHC molecules (6) (reviewed in (7)). Nonproteasomal proteolytic pathways were also suggested as possible contributors to the MHC peptidome, including proteolysis by the ER resident Signal peptide peptidase (8, 9), the cytoplasmic proteases Insulin degrading enzyme (10), Tripeptidyl peptidase (11–13), and a number of proteases within the endolysosome pathway (reviewed recently in (14–17)). In contrast to the mostly cytoplasmic and ER production of the MHC class I peptidome, the class II peptidome is produced in a special compartment, associated with the endolysosome pathway (18–20). This pathway is also thought to participate in the cross presentation of class I peptides, derived from proteins up-taken by professional antigen presenting cells (21), (reviewed in (15–17, 22)).The centrality of the proteasomes in the generation of the MHC peptidome was deduced mostly from the observed change in presentation levels of small numbers of selected peptides, following proteasome inhibition (3, 23). Even the location of some of the genes encoding the catalytic subunits of the immunoproteasome (LMP2 and LMP7) (24) within the MHC class II genomic locus, was suggested to support the involvement of the proteasome in the generation of the MHC class I peptidome (25). Similar conclusions were deduced from alterations in peptide presentation, following expression of the catalytic subunits of the immunoproteasome (26), (reviewed in (5)). Yet, although most of the reports indicated reductions in presentation of selected peptides by proteasome inhibition (3, 27–29), others have observed only limited, and sometimes even opposite effects (23, 30–32).The matter is further complicated by the indirect effects of proteasome inhibition used for such studies on the arrest of protein synthesis by the cells (33–35), on the transport rates of the MHC molecules to the cell surface, and on their retrograde transport back to the vesicular system (36) (reviewed in (37)). Proteasome inhibition likely causes shortage of free ubiquitin, reduced supply of free amino acids, and induces an ER unfolded protein response (UPR), which signals the cells to block most (but not all) cellular protein synthesis (reviewed in (38)). Because a significant portion of the MHC peptidome originates from degradation of DRiPs and SLiPs (reviewed in (2)), arrest of new protein synthesis should influence the presentation of their derived MHC peptides. Taken together, these arguments may suggest that merely following the changes in the presentation levels of the MHC molecules, or even of specific MHC peptides, after proteasome inhibition, does not provide the full picture for deducing the relative contribution of the proteasomal pathway to the production of the MHC peptidome (reviewed in (7)).MHC peptidome analysis can be performed relatively easily by modern capillary chromatography combined with mass spectrometry (reviewed in (39)). The peptides are recovered from immunoaffinity purified MHC molecules after detergent solubilization of the cells (40, 41), from soluble MHC molecules secreted to the cells'' growth medium (42, 43) or from patients'' plasma (44). The purified peptides pools are resolved by capillary chromatography and the individual peptides are identified and quantified by tandem mass spectrometry (40), (reviewed in (45–47)). In cultured cells, quantitative analysis can also be followed by metabolic incorporation of stable isotope labeled amino acids (SILAC) (48). Furthermore, the rates of de novo synthesis of both MHC peptides and their proteins of origin can be followed using the dynamic-SILAC proteomics approach (49) with its further adaptation to HLA peptidomics (50–52).This study attempts to define the relative contribution of the proteasomes to the production of HLA class I peptidome by simultaneously following the effects of proteasome inhibitors, epoxomicin and bortezomib (Velcade), on the rates of de novo synthesis of both the HLA class I peptidome and the cellular proteome of the same MCF-7 human breast cancer cultured cells. The proteasome inhibitors did not reduce the levels of HLA presentations, yet affected the rates of production of both the HLA peptidome and cellular proteome, mostly decreasing, but also increasing some of the synthesis rates of the HLA peptides and cellular proteins. Thus, we suggest that the degree of contribution of the proteasomal pathway to the production of the HLA-I peptidome should be re-evaluated in accordance with their effects on the entire HLA class-I peptidome, while taking into consideration the inhibitors'' effects on the synthesis (and degradation) rates of the source proteins of each of the studied HLA peptides.     

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.     

Proinsulin and the Genetics of Diabetes Mellitus

Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.     

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]     

Tyrosine Phosphorylation of Growth Factor Receptor-bound Protein-7 by Focal Adhesion Kinase in the Regulation of Cell Migration, Proliferation, and Tumorigenesis

We have previously reported that growth factor receptor-bound protein-7 (Grb7), an Src-homology 2 (SH2)-containing adaptor protein, enables interaction with focal adhesion kinase (FAK) to regulate cell migration in response to integrin activation. To further elucidate the signaling events mediated by FAK·Grb7 complexes in promoting cell migration and other cellular functions, we firstly examined the phos pho ryl a ted tyrosine site(s) of Grb7 by FAK using an in vivo mutagenesis. We found that FAK was capable of phos pho rylating at least 2 of 12 tyrosine residues within Grb7, Tyr-188 and Tyr-338. Moreover, mutations converting the identified Tyr to Phe inhibited integrin-dependent cell migration as well as impaired cell proliferation but not survival compared with the wild-type control. Interestingly, the above inhibitory effects caused by the tyrosine phos pho ryl a tion-deficient mutants are probably attributed to their down-regulation of phospho-Tyr-397 of FAK, thereby implying a mechanism by competing with wild-type Grb7 for binding to FAK. Consequently, these tyrosine phos pho ryl a tion-deficient mutants evidently altered the phospho-Tyr-118 of paxillin and phos pho ryl a tion of ERK1/2 but less on phospho-Ser-473 of AKT, implying their involvement in the FAK·Grb7-mediated cellular functions. Additionally, we also illustrated that the formation of FAK·Grb7 complexes and Grb7 phos pho ryl a tion by FAK in an integrin-dependent manner were essential for cell migration, proliferation and anchorage-independent growth in A431 epidermal carcinoma cells, indicating the importance of FAK·Grb7 complexes in tumorigenesis. Our data provide a better understanding on the signal transduction event for FAK·Grb7-mediated cellular functions as well as to shed light on a potential therapeutic in cancers.Growth factor receptor bound protein-7 (Grb7)² is initially identified as a SH2 domain-containing adaptor protein bound to the activated EGF receptor (1). Grb7 is composed of an N-terminal proline-rich region, following a putative RA (Ras-associating) domain and a central PH (pleckstrin homology) domain and a BPS motif (between PH and SH2 domains), and a C-terminal SH2 domain (2–6). Despite the lack of enzymatic activity, the presence of multiple protein-protein interaction domains allows Grb7 family adaptor proteins to participate in versatile signal transduction pathways and, therefore, to regulate many cellular functions (4–6). A number of signaling molecules has been reported to interact with these featured domains, although most of the identified Grb7 binding partners are mediated through its SH2 domain. For example, the SH2 domain of Grb7 has been demonstrated to be capable of binding to the phospho-tyrosine sites of EGF receptor (1), ErbB2 (7), ErbB3 and ErbB4 (8), Ret (9), platelet-derived growth factor receptor (10), insulin receptor (11), SHPTP2 (12), Tek/Tie2 (13), caveolin (14), c-Kit (15), EphB1 (16), G6f immunoreceptor protein (17), Rnd1 (18), Shc (7), FAK (19), and so on. The proceeding α-helix of the PH domain of Grb7 is the calmodulin-binding domain responsible for recruiting Grb7 to plasma membrane in a Ca²⁺-dependent manner (20), and the association between the PH domain of Grb7 and phosphoinositides is required for the phosphorylation by FAK (21). Two additional proteins, NIK (nuclear factor κB-inducing kinase) and FHL2 (four and half lim domains isoform 2), in association with the GM region (Grb and Mig homology region) of Grb7 are also reported, although the physiological functions for these interactions remain unknown (22, 23). Recently, other novel roles in translational controls and stress responses through the N terminus of Grb7 are implicated for the findings of Grb7 interacting with the 5′-untranslated region of capped targeted KOR (kappa opioid receptor) mRNA and the Hu antigen R of stress granules in an FAK-mediated phosphorylation manner (24, 25).Unlike its member proteins Grb10 and Grb14, the role of Grb7 in cell migration is unambiguous and well documented. This is supported by a series of studies. Firstly, Grb7 family members share a significantly conserved molecular architecture with the Caenorhabditis elegans Mig-10 protein, which is involved in neuronal cell migration during embryonic development (4, 5, 26), suggesting that Grb7 may play a role in cell migration. Moreover, Grb7 is often co-amplified with Her2/ErbB2 in certain human cancers and tumor cell lines (7, 27, 28), and its overexpression resulted in invasive and metastatic consequences of various cancers and tumor cells (23, 29–33). On the contrary, knocking down Grb7 by RNA interference conferred to an inhibitory outcome of the breast cancer motility (34). Furthermore, interaction of Grb7 with autophosphorylated FAK at Tyr-397 could promote integrin-mediated cell migration in NIH 3T3 and CHO cells, whereas overexpression of its SH2 domain, an dominant negative mutant of Grb7, inhibited cell migration (19, 35). Recruitment and phosphorylation of Grb7 by EphB1 receptors enhanced cell migration in an ephrin-dependent manner (16). Recently, G7–18NATE, a selective Grb7-SH2 domain affinity cyclic peptide, was demonstrated to efficiently block cell migration of tumor cells (32, 36). In addition to cell migration, Grb7 has been shown to play a role in a variety of physiological and pathological events, for instance, kidney development (37), tumorigenesis (7, 14, 38–41), angiogenic activity (20), proliferation (34, 42, 43), anti-apoptosis (44), gene expression regulation (24), Silver-Russell syndrome (45), rheumatoid arthritis (46), atopic dermatitis (47), and T-cell activation (17, 48). Nevertheless, it remains largely unknown regarding the downstream signaling events of Grb7-mediated various functions. In particular, given the role of Grb7 as an adaptor molecule and its SH2 domain mainly interacting with upstream regulators, it will be interesting to identify potential downstream effectors through interacting with the functional GM region or N-terminal proline-rich region.In this report, we identified two tyrosine phosphorylated sites of Grb7 by FAK and deciphered the signaling targets downstream through these phosphorylated tyrosine sites to regulate various cellular functions such as cell migration, proliferation, and survival. In addition, our study sheds light on tyrosine phosphorylation of Grb7 by FAK involved in tumorigenesis.     

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]     

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]