Phosphorylation is universally used for controlling protein function, but knowledge of the phosphoproteome in vertebrate embryos has been limited. However, recent technical advances make it possible to define an organism's phosphoproteome at a more comprehensive level. Xenopus laevis offers established advantages for analyzing the regulation of protein function by phosphorylation. Functionally unbiased, comprehensive information about the Xenopus phosphoproteome would provide a powerful guide for future studies of phosphorylation in a developmental context. To this end, we performed a phosphoproteomic analysis of Xenopus oocytes, eggs, and embryos using recently developed mass spectrometry methods. We identified 1,441 phosphorylation sites present on 654 different Xenopus proteins, including hundreds of previously unknown phosphorylation sites. This approach identified several phosphorylation sites described in the literature and/or evolutionarily conserved in other organisms, validating the data's quality. These data will serve as a powerful resource for the exploration of phosphorylation and protein function within a developmental context. Developmental Dynamics 238:1433-1443, 2009. (c) 2009 Wiley-Liss, Inc.

Protein post-translational modifications (PTMs) are a fast and versatility mechanism used by the cell to regulate the function of proteins in response to changing conditions. PTMs can alter the activity of proteins by allosteric regulation or by controlling protein interactions, localization and abundance. Recent advances in proteomics have revealed the extent of regulation by PTMs and the different mechanisms used in nature to exert control over protein function via PTMs. These developments can serve as the foundation for the rational design of protein regulation. Here we review the advances in methods to determine the function of PTMs, protein allosteric control and examples of rational design of PTM regulation. These advances create an opportunity to move synthetic biology forward by making use of a level of regulation that is of yet unexplored.

The African clawed frog Xenopus laevis is an important model organism for studies in developmental and cell biology, including cell-signaling. However, our knowledge of X. laevis protein post-translational modifications remains scarce. Here, we used a mass spectrometry-based approach to survey the phosphoproteome of this species, compiling a list of 2636 phosphosites. We used structural information and phosphoproteomic data for 13 other species in order to predict functionally important phospho-regulatory events. We found that the degree of conservation of phosphosites across species is predictive of sites with known molecular function. In addition, we predicted kinase-protein interactions for a set of cell-cycle kinases across all species. The degree of conservation of kinase-protein interactions was found to be predictive of functionally relevant regulatory interactions. Finally, using comparative protein structure models, we find that phosphosites within structured domains tend to be located at positions with high conformational flexibility. Our analysis suggests that a small class of phosphosites occurs in positions that have the potential to regulate protein conformation.

Various physiological stimuli, such as cold environment, diet, and hormones, trigger brown adipose tissue (BAT) to produce heat through sympathetic nervous system (SNS)- and β-adrenergic receptors (βARs). The βAR stimulation increases intracellular cAMP levels through heterotrimeric G proteins and adenylate cyclases, but the processes by which cAMP modulates brown adipocyte function are not fully understood. Here we described that specific ablation of cAMP production in brown adipocytes led to reduced lipolysis, mitochondrial biogenesis, uncoupling protein 1 (Ucp1) expression, and consequently defective adaptive thermogenesis. Elevated cAMP signaling by sympathetic activation inhibited Salt-inducible kinase 2 (Sik2) through protein kinase A (PKA)-mediated phosphorylation in brown adipose tissue. Inhibition of SIKs enhanced Ucp1 expression in differentiated brown adipocytes and Sik2 knockout mice exhibited enhanced adaptive thermogenesis at thermoneutrality in an Ucp1-dependent manner. Taken together, our data indicate that suppressing Sik2 by PKA-mediated phosphorylation is a requisite for SNS-induced Ucp1 expression and adaptive thermogenesis in BAT, and targeting Sik2 may present a novel therapeutic strategy to ramp up BAT thermogenic activity in humans.

Enteroviruses (EVs) comprise a large genus of positive-sense, single-stranded RNA viruses whose members cause a number of important and widespread human diseases, including poliomyelitis, myocarditis, acute flaccid myelitis and the common cold. How EVs co-opt cellular functions to promote replication and spread is incompletely understood. Here, using genome-scale CRISPR screens, we identify the actin histidine methyltransferase SET domain containing 3 (SETD3) as critically important for viral infection by a broad panel of EVs, including rhinoviruses and non-polio EVs increasingly linked to severe neurological disease such as acute flaccid myelitis (EV-D68) and viral encephalitis (EV-A71). We show that cytosolic SETD3, independent of its methylation activity, is required for the RNA replication step in the viral life cycle. Using quantitative affinity purification-mass spectrometry, we show that SETD3 specifically interacts with the viral 2A protease of multiple enteroviral species, and we map the residues in 2A that mediate this interaction. 2A mutants that retain protease activity but are unable to interact with SETD3 are severely compromised in RNA replication. These data suggest a role of the viral 2A protein in RNA replication beyond facilitating proteolytic cleavage. Finally, we show that SETD3 is essential for in vivo replication and pathogenesis in multiple mouse models for EV infection, including CV-A10, EV-A71 and EV-D68. Our results reveal a crucial role of a host protein in viral pathogenesis, and suggest targeting SETD3 as a potential mechanism for controlling viral infections.

The genus Escherichia is composed of several species and cryptic clades, including E. coli, which behaves as a vertebrate gut commensal, but also as an opportunistic pathogen involved in both diarrheic and extra-intestinal diseases. To characterize the genetic determinants of extra-intestinal virulence within the genus, we carried out an unbiased genome-wide association study (GWAS) on 370 commensal, pathogenic and environmental strains representative of the Escherichia genus phylogenetic diversity and including E. albertii (n = 7), E. fergusonii (n = 5), Escherichia clades (n = 32) and E. coli (n = 326), tested in a mouse model of sepsis. We found that the presence of the high-pathogenicity island (HPI), a ~35 kbp gene island encoding the yersiniabactin siderophore, is highly associated with death in mice, surpassing other associated genetic factors also related to iron uptake, such as the aerobactin and the sitABCD operons. We confirmed the association in vivo by deleting key genes of the HPI in E. coli strains in two phylogenetic backgrounds. We then searched for correlations between virulence, iron capture systems and in vitro growth in a subset of E. coli strains (N = 186) previously phenotyped across growth conditions, including antibiotics and other chemical and physical stressors. We found that virulence and iron capture systems are positively correlated with growth in the presence of numerous antibiotics, probably due to co-selection of virulence and resistance. We also found negative correlations between virulence, iron uptake systems and growth in the presence of specific antibiotics (i.e. cefsulodin and tobramycin), which hints at potential "collateral sensitivities" associated with intrinsic virulence. This study points to the major role of iron capture systems in the extra-intestinal virulence of the genus Escherichia.

Complex networks of regulatory relationships between protein kinases comprise a major component of intracellular signaling. Although many kinase-kinase regulatory relationships have been described in detail, these tend to be limited to well-studied kinases whereas the majority of possible relationships remains unexplored. Here, we implement a data-driven, supervised machine learning method to predict human kinase-kinase regulatory relationships and whether they have activating or inhibiting effects. We incorporate high-throughput data, kinase specificity profiles, and structural information to produce our predictions. The results successfully recapitulate previously annotated regulatory relationships and can reconstruct known signaling pathways from the ground up. The full network of predictions is relatively sparse, with the vast majority of relationships assigned low probabilities. However, it nevertheless suggests denser modes of inter-kinase regulation than normally considered in intracellular signaling research. A record of this paper's transparent peer review process is included in the Supplemental Information.

Amino acids fulfil a diverse range of roles in proteins, each utilising its chemical properties in different ways in different contexts to create required functions. For example, cysteines form disulphide or hydrogen bonds in different circumstances and charged amino acids do not always make use of their charge. The repertoire of amino acid functions and the frequency at which they occur in proteins remains understudied. Measuring large numbers of mutational consequences, which can elucidate the role an amino acid plays, was prohibitively time-consuming until recent developments in deep mutational scanning. In this study, we gathered data from 28 deep mutational scanning studies, covering 6,291 positions in 30 proteins, and used the consequences of mutation at each position to define a mutational landscape. We demonstrated rich relationships between this landscape and biophysical or evolutionary properties. Finally, we identified 100 functional amino acid subtypes with a data-driven clustering analysis and studied their features, including their frequencies and chemical properties such as tolerating polarity, hydrophobicity or being intolerant of charge or specific amino acids. The mutational landscape and amino acid subtypes provide a foundational catalogue of amino acid functional diversity, which will be refined as the number of studied protein positions increases.

Cells inherit two centrioles, the older of which is uniquely capable of generating a cilium. Using proteomics and superresolved imaging, we identify a module that we term DISCO (distal centriole complex). The DISCO components CEP90, MNR, and OFD1 underlie human ciliopathies. This complex localizes to both distal centrioles and centriolar satellites, proteinaceous granules surrounding centrioles. Cells and mice lacking CEP90 or MNR do not generate cilia, fail to assemble distal appendages, and do not transduce Hedgehog signals. Disrupting the satellite pools does not affect distal appendage assembly, indicating that it is the centriolar populations of MNR and CEP90 that are critical for ciliogenesis. CEP90 recruits the most proximal known distal appendage component, CEP83, to root distal appendage formation, an early step in ciliogenesis. In addition, MNR, but not CEP90, restricts centriolar length by recruiting OFD1. We conclude that DISCO acts at the distal centriole to support ciliogenesis by restraining centriole length and assembling distal appendages, defects in which cause human ciliopathies.

Genome-wide association studies have discovered numerous genomic loci associated with Alzheimer's disease (AD); yet the causal genes and variants are incompletely identified. We performed an updated genome-wide AD meta-analysis, which identified 37 risk loci, including new associations near CCDC6, TSPAN14, NCK2 and SPRED2. Using three SNP-level fine-mapping methods, we identified 21 SNPs with >50% probability each of being causally involved in AD risk and others strongly suggested by functional annotation. We followed this with colocalization analyses across 109 gene expression quantitative trait loci datasets and prioritization of genes by using protein interaction networks and tissue-specific expression. Combining this information into a quantitative score, we found that evidence converged on likely causal genes, including the above four genes, and those at previously discovered AD loci, including BIN1, APH1B, PTK2B, PILRA and CASS4.

Drug repurposing provides a rapid approach to meet the urgent need for therapeutics to address COVID-19. To identify therapeutic targets relevant to COVID-19, we conducted Mendelian randomization analyses, deriving genetic instruments based on transcriptomic and proteomic data for 1,263 actionable proteins that are targeted by approved drugs or in clinical phase of drug development. Using summary statistics from the Host Genetics Initiative and the Million Veteran Program, we studied 7,554 patients hospitalized with COVID-19 and >1 million controls. We found significant Mendelian randomization results for three proteins (ACE2, P = 1.6 × 10-6; IFNAR2, P = 9.8 × 10-11 and IL-10RB, P = 2.3 × 10-14) using cis-expression quantitative trait loci genetic instruments that also had strong evidence for colocalization with COVID-19 hospitalization. To disentangle the shared expression quantitative trait loci signal for IL10RB and IFNAR2, we conducted phenome-wide association scans and pathway enrichment analysis, which suggested that IFNAR2 is more likely to play a role in COVID-19 hospitalization. Our findings prioritize trials of drugs targeting IFNAR2 and ACE2 for early management of COVID-19.

The Hedgehog (Hh) pathway is essential for organ development, homeostasis, and regeneration. Dysfunction of this cascade drives several cancers. To control expression of pathway target genes, the G protein-coupled receptor (GPCR) Smoothened (SMO) activates glioma-associated (GLI) transcription factors via an unknown mechanism. Here, we show that, rather than conforming to traditional GPCR signaling paradigms, SMO activates GLI by binding and sequestering protein kinase A (PKA) catalytic subunits at the membrane. This sequestration, triggered by GPCR kinase (GRK)-mediated phosphorylation of SMO intracellular domains, prevents PKA from phosphorylating soluble substrates, releasing GLI from PKA-mediated inhibition. Our work provides a mechanism directly linking Hh signal transduction at the membrane to GLI transcription in the nucleus. This process is more fundamentally similar between species than prevailing hypotheses suggest. The mechanism described here may apply broadly to other GPCR- and PKA-containing cascades in diverse areas of biology.

The epidermal growth factor receptor (EGFR) inhibitor cetuximab is the only FDA-approved oncogene-targeting therapy for head and neck squamous cell carcinoma (HNSCC). Despite variable treatment response, no biomarkers exist to stratify patients for cetuximab therapy in HNSCC. Here, we applied unbiased hierarchical clustering to reverse-phase protein array molecular profiles from patient-derived xenograft (PDX) tumors and revealed 2 PDX clusters defined by protein networks associated with EGFR inhibitor resistance. In vivo validation revealed unbiased clustering to classify PDX tumors according to cetuximab response with 88% accuracy. Next, a support vector machine classifier algorithm identified a minimalist biomarker signature consisting of 8 proteins - caveolin-1, Sox-2, AXL, STING, Brd4, claudin-7, connexin-43, and fibronectin - with expression that strongly predicted cetuximab response in PDXs using either protein or mRNA. A combination of caveolin-1 and Sox-2 protein levels was sufficient to maintain high predictive accuracy, which we validated in tumor samples from patients with HNSCC with known clinical response to cetuximab. These results support further investigation into the combined use of caveolin-1 and Sox-2 as predictive biomarkers for cetuximab response in the clinic.

Cold-induced thermogenesis in endotherms demands adaptive thermogenesis fueled by mitochondrial respiration and Ucp1-mediated uncoupling in multilocular brown adipocytes (BAs). However, dietary regulation of thermogenesis in BAs isn't fully understood. Here, we describe that the deficiency of Leucine-rich pentatricopeptide repeat containing-protein (Lrpprc) in BAs reduces mtDNA-encoded ETC gene expression, causes ETC proteome imbalance, and abolishes the mitochondria-fueled thermogenesis. BA-specific Lrpprc knockout mice are cold resistant in a 4°C cold-tolerance test in the presence of food, which is accompanied by the activation of transcription factor 4 (ATF4) and proteome turnover in BAs. ATF4 activation genetically by BA-specific ATF4 overexpression or physiologically by a low-protein diet feeding can improve cold tolerance in wild-type and Ucp1 knockout mice. Furthermore, ATF4 activation in BAs improves systemic metabolism in obesogenic environment regardless of Ucp1's action. Therefore, our study reveals a diet-dependent but Ucp1-independent thermogenic mechanism in BAs that is relevant to systemic thermoregulation and energy homeostasis.

Reversible protein phosphorylation is an important mechanism for regulating (dis)assembly of biomolecular condensates. However, condensate-specific phosphosites remain largely unknown, thereby limiting our understanding of the underlying mechanisms. Here, we combine solubility proteome profiling with phosphoproteomics to quantitatively map several hundred phosphosites enriched in either soluble or condensate-bound protein subpopulations, including a subset of phosphosites modulating protein-RNA interactions. We show that multi-phosphorylation of the C-terminal disordered segment of heteronuclear ribonucleoprotein A1 (HNRNPA1), a key RNA-splicing factor, reduces its ability to locate to nuclear clusters. For nucleophosmin 1 (NPM1), an essential nucleolar protein, we show that phosphorylation of S254 and S260 is crucial for lowering its partitioning to the nucleolus and additional phosphorylation of distal sites enhances its retention in the nucleoplasm. These phosphorylation events decrease RNA and protein interactions of NPM1 to regulate its condensation. Our dataset is a rich resource for systematically uncovering the phosphoregulation of biomolecular condensates.

Cellular functions are governed by molecular machines that assemble through protein-protein interactions. Their atomic details are critical to studying their molecular mechanisms. However, fewer than 5% of hundreds of thousands of human protein interactions have been structurally characterized. Here we test the potential and limitations of recent progress in deep-learning methods using AlphaFold2 to predict structures for 65,484 human protein interactions. We show that experiments can orthogonally confirm higher-confidence models. We identify 3,137 high-confidence models, of which 1,371 have no homology to a known structure. We identify interface residues harboring disease mutations, suggesting potential mechanisms for pathogenic variants. Groups of interface phosphorylation sites show patterns of co-regulation across conditions, suggestive of coordinated tuning of multiple protein interactions as signaling responses. Finally, we provide examples of how the predicted binary complexes can be used to build larger assemblies helping to expand our understanding of human cell biology.

Genetic alterations that activate protein kinase A (PKA) are found in many tumor types. Yet, their downstream oncogenic signaling mechanisms are poorly understood. We used global phosphoproteomics and kinase activity profiling to map conserved signaling outputs driven by a range of genetic changes that activate PKA in human cancer. Two signaling networks were identified downstream of PKA: RAS/MAPK components and an Aurora Kinase A (AURKA)/glycogen synthase kinase (GSK3) sub-network with activity toward MYC oncoproteins. Findings were validated in two PKA-dependent cancer models: a novel, patient-derived fibrolamellar carcinoma (FLC) line that expresses a DNAJ-PKAc fusion and a PKA-addicted melanoma model with a mutant type I PKA regulatory subunit. We identify PKA signals that can influence both de novo translation and stability of the proto-oncogene c-MYC. However, the primary mechanism of PKA effects on MYC in our cell models was translation and could be blocked with the eIF4A inhibitor zotatifin. This compound dramatically reduced c-MYC expression and inhibited FLC cell line growth in vitro. Thus, targeting PKA effects on translation is a potential treatment strategy for FLC and other PKA-driven cancers.

Apolipoprotein (apo) E4 is the major genetic risk factor for Alzheimer's disease. While neurons generally produce a minority of the apoE in the central nervous system, neuronal expression of apoE increases dramatically in response to stress and is sufficient to drive pathology. Currently, the molecular mechanisms of how apoE4 expression may regulate pathology are not fully understood. Here, we expand upon our previous studies measuring the impact of apoE4 on protein abundance to include the analysis of protein phosphorylation and ubiquitylation signaling in isogenic Neuro-2a cells expressing apoE3 or apoE4. ApoE4 expression resulted in a dramatic increase in vasodilator-stimulated phosphoprotein (VASP) S235 phosphorylation in a protein kinase A (PKA)-dependent manner. This phosphorylation disrupted VASP interactions with numerous actin cytoskeletal and microtubular proteins. Reduction of VASP S235 phosphorylation via PKA inhibition resulted in a significant increase in filopodia formation and neurite outgrowth in apoE4-expressing cells, exceeding levels observed in apoE3-expressing cells. Our results highlight the pronounced and diverse impact of apoE4 on multiple modes of protein regulation and identify protein targets to restore apoE4-related cytoskeletal defects.

A widespread strategy employed by pathogens to establish infection is to inhibit host-cell protein synthesis. Legionella pneumophila, an intracellular bacterial pathogen and the causative organism of Legionnaires' disease, secretes a subset of protein effectors into host cells that inhibit translation elongation. Mechanistic insights into how the bacterium targets translation elongation remain poorly defined. We report here that the Legionella effector SidI functions in an unprecedented way as a transfer-RNA mimic that directly binds to and glycosylates the ribosome. The 3.1 Å cryo-electron microscopy structure of SidI reveals an N-terminal domain with an 'inverted L' shape and surface-charge distribution characteristic of tRNA mimicry, and a C-terminal domain that adopts a glycosyl transferase fold that licenses SidI to utilize GDP-mannose as a sugar precursor. This coupling of tRNA mimicry and enzymatic action endows SidI with the ability to block protein synthesis with a potency comparable to ricin, one of the most powerful toxins known. In Legionella-infected cells, the translational pausing activated by SidI elicits a stress response signature mimicking the ribotoxic stress response, which is activated by elongation inhibitors that induce ribosome collisions. SidI-mediated effects on the ribosome activate the stress kinases ZAKα and p38, which in turn drive an accumulation of the protein activating transcription factor 3 (ATF3). Intriguingly, ATF3 escapes the translation block imposed by SidI, translocates to the nucleus and orchestrates the transcription of stress-inducible genes that promote cell death, revealing a major role for ATF3 in the response to collided ribosome stress. Together, our findings elucidate a novel mechanism by which a pathogenic bacterium employs tRNA mimicry to hijack a ribosome-to-nuclear signalling pathway that regulates cell fate.

The Cullin-RING E3 ligase (CRL) family is commonly hijacked by pathogens to redirect the host ubiquitin proteasome machinery to specific targets. During HIV infection, CRL5 is hijacked by HIV Vif to target viral restriction factors of the APOBEC3 family for ubiquitination and degradation. Here, using a quantitative proteomics approach, we identify the E3 ligase ARIH2 as a regulator of CRL5-mediated APOBEC3 degradation. The CUL5Vif/CBFß complex recruits ARIH2 where it acts to transfer ubiquitin directly to the APOBEC3 targets. ARIH2 is essential for CRL5-dependent HIV infectivity in primary CD4+ T cells. Furthermore, we show that ARIH2 cooperates with CRL5 to prime other cellular substrates for polyubiquitination, suggesting this may represent a general mechanism beyond HIV infection and APOBEC3 degradation. Taken together, these data identify ARIH2 as a co-factor in the Vif-hijacked CRL5 complex that contributes to HIV infectivity and demonstrate the operation of the E1-E2-E3/E3-substrate ubiquitination mechanism in a viral infection context.