Regeneration Mechanisms of Arabidopsis thaliana Methionine Sulfoxide Reductases B by Glutaredoxins and Thioredoxins

Methionine oxidation leads to the formation of S- and R-diastereomers of methionine sulfoxide (MetSO), which are reduced back to methionine by methionine sulfoxide reductases (MSRs) A and B, respectively. MSRBs are classified in two groups depending on the conservation of one or two redox-active Cys; 2-Cys MSRBs possess a catalytic Cys-reducing MetSO and a resolving Cys, allowing regeneration by thioredoxins. The second type, 1-Cys MSRBs, possess only the catalytic Cys. The biochemical mechanisms involved in activity regeneration of 1-Cys MSRBs remain largely elusive. In the present work we used recombinant plastidial Arabidopsis thaliana MSRB1 and MSRB2 as models for 1-Cys and 2-Cys MSRBs, respectively, to delineate the Trx- and glutaredoxin-dependent reduction mechanisms. Activity assays carried out using a series of cysteine mutants and various reductants combined with measurements of free thiols under distinct oxidation conditions and mass spectrometry experiments show that the 2-Cys MSRB2 is reduced by Trx through a dithiol-disulfide exchange involving both redox-active Cys of the two partners. Regarding 1-Cys MSRB1, oxidation of the enzyme after substrate reduction leads to the formation of a stable sulfenic acid on the catalytic Cys, which is subsequently glutathionylated. The deglutathionylation of MSRB1 is achieved by both mono- and dithiol glutaredoxins and involves only their N-terminal conserved catalytic Cys. This study proposes a detailed mechanism of the regeneration of 1-Cys MSRB activity by glutaredoxins, which likely constitute physiological reductants for this type of MSR.Proteins are prone to oxidative modifications due to the action of reactive oxygen species. Methionine (Met), one of the most susceptible amino acids to oxidation (1), is converted into methionine sulfoxide (MetSO),³ resulting in altered conformation and activity for many proteins (1). Methionine sulfoxide reductases (MSRs), which catalyze the reduction of MetSO back to Met, are present in most living organisms. MSRA, the first MSR isolated (2), is specific of the MetSO S-diastereomer and participates in protection against oxidative stress (3). A second MSR type, MSRB, which catalytically reduces the MetSO R-diastereomer, was identified later (4). MSRA and MSRB are monomeric enzymes that display no sequence or structural homologies but share a similar three-step catalytic mechanism, (i) reduction of MetSO by MSR and formation of a sulfenic acid intermediate on the “catalytic” cysteine (Cys), (ii) formation of a disulfide bond between catalytic and “resolving” Cys and release of H₂O, and (iii) reduction of the disulfide bond by a reductant (5, 6). Thioredoxins (Trxs) have been proposed to be the biological reductant for MSRs (2, 7). Trxs are small and ubiquitous disulfide reductases with a WC(G/P)PC active site. They function as electron donors and play essential roles in many processes through control of protein conformation and activity by supplying the reducing power needed to reduce disulfide bonds in target proteins.Most MSRBs, named 2-Cys MSRBs, possess two conserved Cys and are actually reduced by Trxs (7, 8). However, in a second class of MSRBs, termed 1-Cys MSRBs and representing ∼40% of known MSRBs, the resolving Cys residue corresponding to Cys-63 in Escherichia coli is replaced by Thr or Ser (8, 9). Although some of these MSRBs possess another potential resolving Cys (9), most 1-Cys MSRBs do not have any additional Cys, indicating that an alternative mechanism, which does not involve the formation of an intramolecular disulfide reduction, is needed for their regeneration (7). Contrasting data concerning the role of Trxs in providing electrons to these MSRBs have been reported. Several studies showed that cytosolic Trx is not an efficient reductant for human 1-Cys MSRBs (10–12), whereas mitochondrial Trxs were recently reported to efficiently regenerate mitochondrial 1-Cys MSRBs (13). It has been proposed that regeneration of mammalian and plant 1-Cys MSRBs could involve direct reduction of the cysteine sulfenic acid form generated during catalysis (10, 13–15).Arabidopsis thaliana possesses two plastidial MSRBs referred to as MSRB1 and MSRB2 and related to 1-Cys and 2-Cys MSRB types, respectively. MSRB2 possesses two CXXC motifs potentially implicated in the coordination of a zinc atom, a Cys in position 187 corresponding to the catalytic Cys-117 of E. coli MSRB, a potential resolving Cys in position 134, and an additional Cys in position 115. MSRB1 also contains the four Cys residues potentially coordinating zinc, the potential catalytic Cys-186, and a Cys in position 116, whereas the potential resolving Cys is replaced by a Thr in position 132. Previously, we showed that various types of canonical Trxs are efficient electron suppliers to MSRB2, whereas MSRB1 can only be reduced by the peculiar Trx CDSP32 (chloroplastic drought-induced stress protein of 32 kDa) and by Grxs (15–17). Grxs are oxidoreductases of the Trx superfamily possessing either a monothiol CXXS or a dithiol CXXC active site and are generally reduced by glutathione (18). Grxs are able to reduce protein disulfides, but also glutathione-mixed disulfides, a reaction termed deglutathionylation, for which Trxs are not efficient catalysts (19, 20). Classical dithiol Grxs can reduce disulfide bonds using both active site Cys residues, as shown for E. coli ribonucleotide reductase, but can also reduce glutathione-mixed disulfides through a monothiol mechanism that requires only the N-terminal active site Cys (21). CXXS-type Grxs catalyze deglutathionylation either through a monothiol mechanism, as recently shown for chloroplastic GrxS12 (CSYS active site) (22), or through a dithiol mechanism as suggested for Grxs with a CGFS active site (20, 23).We reported recently the involvement of Grxs in the regeneration of MSRB activity (15). Nevertheless, the precise biochemical mechanism underlying regeneration by Grxs remains unknown. In this study we performed a detailed analysis of the roles of redox-active Cys in reductants (Trxs and Grxs) and in acceptors (plastidial Arabidopsis MSRBs). We provide evidence that reduction of MSRB2 by Trxs is achieved through a classical dithiol-disulfide exchange. The data on MSRB1 reveal that 1-Cys MSRBs are regenerated by Grxs through a glutathionylation step of the catalytic Cys.     

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]     

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]     

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.     

Analysis of Proteins That Rapidly Change Upon Mechanistic/Mammalian Target of Rapamycin Complex 1 (mTORC1) Repression Identifies Parkinson Protein 7 (PARK7) as a Novel Protein Aberrantly Expressed in Tuberous Sclerosis Complex (TSC)

Many biological processes involve the mechanistic/mammalian target of rapamycin complex 1 (mTORC1). Thus, the challenge of deciphering mTORC1-mediated functions during normal and pathological states in the central nervous system is challenging. Because mTORC1 is at the core of translation, we have investigated mTORC1 function in global and regional protein expression. Activation of mTORC1 has been generally regarded to promote translation. Few but recent works have shown that suppression of mTORC1 can also promote local protein synthesis. Moreover, excessive mTORC1 activation during diseased states represses basal and activity-induced protein synthesis. To determine the role of mTORC1 activation in protein expression, we have used an unbiased, large-scale proteomic approach. We provide evidence that a brief repression of mTORC1 activity in vivo by rapamycin has little effect globally, yet leads to a significant remodeling of synaptic proteins, in particular those proteins that reside in the postsynaptic density. We have also found that curtailing the activity of mTORC1 bidirectionally alters the expression of proteins associated with epilepsy, Alzheimer''s disease, and autism spectrum disorder—neurological disorders that exhibit elevated mTORC1 activity. Through a protein–protein interaction network analysis, we have identified common proteins shared among these mTORC1-related diseases. One such protein is Parkinson protein 7, which has been implicated in Parkinson''s disease, yet not associated with epilepsy, Alzheimers disease, or autism spectrum disorder. To verify our finding, we provide evidence that the protein expression of Parkinson protein 7, including new protein synthesis, is sensitive to mTORC1 inhibition. Using a mouse model of tuberous sclerosis complex, a disease that displays both epilepsy and autism spectrum disorder phenotypes and has overactive mTORC1 signaling, we show that Parkinson protein 7 protein is elevated in the dendrites and colocalizes with the postsynaptic marker postsynaptic density-95. Our work offers a comprehensive view of mTORC1 and its role in regulating regional protein expression in normal and diseased states.The mechanistic/mammalian target of rapamycin complex 1 (mTORC1)¹ is a serine/threonine protein kinase that is highly expressed in many cell types (1). In the brain, mTORC1 tightly coordinates different synaptic plasticities — long-term potentiation (LTP) and long-term depression (LTD) — the molecular correlates of learning and memory (2–5). Because mTORC1 is at the core of many synaptic signaling pathways downstream of glutamate and neurotrophin receptors, many hypothesize that dysregulated mTORC1 signaling underlies cognitive deficits observed in several neurodegenerative diseases (3, 6–17). For example, mTORC1 and its downstream targets are hyperactive in human brains diagnosed with Alzheimer''s disease (AD) (18–20). Additionally in animal models of autism spectrum disorder (ASD), altered mTORC1 signaling contributes to the observed synaptic dysfunction and aberrant network connectivity (13, 15, 21–27). Furthermore, epilepsy, which is common in AD and ASD, has enhanced mTORC1 activity (28–32).Phosphorylation of mTORC1, considered the active form, is generally regarded to promote protein synthesis (33). Thus, many theorize that diseases with overactive mTORC1 arise from excessive protein synthesis (14). Emerging data, however, show that suppressing mTORC1 activation can trigger local translation in neurons (34, 35). Pharmacological antagonism of N-methyl-d-aspartate (NMDA) receptors, a subtype of glutamate receptors that lies upstream of mTOR activation, promotes the synthesis of the voltage-gated potassium channel, Kv1.1, in dendrites (34, 35). Consistent with these results, in models of temporal lobe epilepsy there is a reduction in the expression of voltage-gated ion channels including Kv1.1 (30, 31, 36). Interestingly in a model of focal neocortical epilepsy, overexpression of Kv1.1 blocked seizure activity (37). Because both active and inactive mTORC1 permit protein synthesis, we sought to determine the proteins whose expression is altered when mTORC1 phosphorylation is reduced in vivo.Rapamycin is an FDA-approved, immunosuppressive drug that inhibits mTORC1 activity (38). We capitalized on the ability of rapamycin to reduce mTORC1 activity in vivo and the unbiased approach of mass spectrometry to identify changes in protein expression. Herein, we provide evidence that mTORC1 activation bidirectionally regulates protein expression, especially in the PSD where roughly an equal distribution of proteins dynamically appear and disappear. Remarkably, using protein–protein interaction networks facilitated the novel discovery that PARK7, a protein thus far only implicated in Parkinson''s disease, (1) is up-regulated by increased mTORC1 activity, (2) resides in the PSD only when mTORC1 is active, and (3) is aberrantly expressed in a rodent model of TSC, an mTORC1-related disease that has symptoms of epilepsy and autism. Collectively, these data provide the first comprehensive list of proteins whose abundance or subcellular distributions are altered with acute changes in mTORC1 activity in vivo.     

A Pipeline for Determining Protein–Protein Interactions and Proximities in the Cellular Milieu

It remains extraordinarily challenging to elucidate endogenous protein-protein interactions and proximities within the cellular milieu. The dynamic nature and the large range of affinities of these interactions augment the difficulty of this undertaking. Among the most useful tools for extracting such information are those based on affinity capture of target bait proteins in combination with mass spectrometric readout of the co-isolated species. Although highly enabling, the utility of affinity-based methods is generally limited by difficulties in distinguishing specific from nonspecific interactors, preserving and isolating all unique interactions including those that are weak, transient, or rapidly exchanging, and differentiating proximal interactions from those that are more distal. Here, we have devised and optimized a set of methods to address these challenges. The resulting pipeline involves flash-freezing cells in liquid nitrogen to preserve the cellular environment at the moment of freezing; cryomilling to fracture the frozen cells into intact micron chunks to allow for rapid access of a chemical reagent and to stabilize the intact endogenous subcellular assemblies and interactors upon thawing; and utilizing the high reactivity of glutaraldehyde to achieve sufficiently rapid stabilization at low temperatures to preserve native cellular interactions. In the course of this work, we determined that relatively low molar ratios of glutaraldehyde to reactive amines within the cellular milieu were sufficient to preserve even labile and transient interactions. This mild treatment enables efficient and rapid affinity capture of the protein assemblies of interest under nondenaturing conditions, followed by bottom-up MS to identify and quantify the protein constituents. For convenience, we have termed this approach Stabilized Affinity Capture Mass Spectrometry. Here, we demonstrate that Stabilized Affinity Capture Mass Spectrometry allows us to stabilize and elucidate local, distant, and transient protein interactions within complex cellular milieux, many of which are not observed in the absence of chemical stabilization.Insights into many cellular processes require detailed information about interactions between the participating proteins. However, the analysis of such interactions can be challenging because of the often-diverse physicochemical properties and the abundances of the constituent proteins, as well as the sometimes wide range of affinities and complex dynamics of the interactions. One of the key challenges has been acquiring information concerning transient, low affinity interactions in highly complex cellular milieux (3, 4).Methods that allow elucidation of such information include co-localization microscopy (5), fluorescence protein Förster resonance energy transfer (4), immunoelectron microscopy (5), yeast two-hybrid (6), and affinity capture (7, 8). Among these, affinity capture (AC)¹ has the unique potential to detect all specific in vivo interactions simultaneously, including those that interact both directly and indirectly. In recent times, the efficacy of such affinity isolation experiments has been greatly enhanced through the use of sensitive modern mass spectrometric protein identification techniques (9). Nevertheless, AC suffers from several shortcomings. These include the problem of 1) distinguishing specific from nonspecific interactors (10, 11); 2) preserving and isolating all unique interactions including those that are weak and/or transient, as well as those that exchange rapidly (10, 12, 13); and 3) differentiating proximal from more distant interactions (14).We describe here an approach to address these issues, which makes use of chemical stabilization of protein assemblies in the complex cellular milieu prior to AC. Chemical stabilization is an emerging technique for stabilizing and elucidating protein associations both in vitro (15–20) and in vivo (3, 12, 14, 21–29), with mass spectrometric (MS) readout of the AC proteins and their connectivities. Such chemical stabilization methods are indeed well-established and are often used in electron microscopy for preserving complexes and subcellular structures both in the cellular milieu (3) and in purified complexes (30, 31), wherein the most reliable, stable, and established stabilization reagents is glutaraldehyde. Recently, glutaraldehyde has been applied in the “GraFix” protocol in which purified protein complexes are subjected to centrifugation through a density gradient that also contains a gradient of glutaraldehyde (30, 31), allowing for optimal stabilization of authentic complexes and minimization of nonspecific associations and aggregation. GraFix has also been combined with mass spectrometry on purified complexes bound to EM grids to obtain a compositional analysis of the complexes (32), thereby raising the possibility that glutaraldehyde can be successfully utilized in conjunction with AC in complex cellular milieux directly.In this work, we present a robust pipeline for determining specific protein-protein interactions and proximities from cellular milieux. The first steps of the pipeline involve the well-established techniques of flash freezing the cells of interest in liquid nitrogen and cryomilling, which have been known for over a decade (33, 34) to preserve the cellular environment, as well as having shown outstanding performance when used in analysis of macromolecular interactions in yeast (35–39), bacterial (40, 41), trypanosome (42), mouse (43), and human (44–47) systems. The resulting frozen powder, composed of intact micron chunks of cells that have great surface area and outstanding solvent accessibility, is well suited for rapid low temperature chemical stabilization using glutaraldehyde. We selected glutaraldehyde for our procedure based on the fact that it is a very reactive stabilizing reagent, even at lower temperatures, and because it has already been shown to stabilize enzymes in their functional state (48–50). We employed highly efficient, rapid, single stage affinity capture (36, 51) for isolation and bottom-up MS for analysis of the macromolecular assemblies of interest (52–54). For convenience, we have termed this approach Stabilized Affinity-Capture Mass Spectrometry (SAC-MS).     

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]     

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.     

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.     

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]