COVID-BIOCHIP: A Web Tool to Analyse COVID-19 Antigen Microarrays.

During the COVID-19 pandemic, there was a growing need to characterise the disease. A very important aspect is the ability to measure the immunisation extent, which can be achieved using antigen microarrays that quantitively measure the presence of COVID-related antibodies. A significant limitation for these tests was the complexity of manually analysing the results, and the limited availability of software for its analysis. In this paper, we describe the development of COVID-BIOCHIP, an ad-hoc web-based solution for the automatic analysis and visualisation of COVID-19 antigen microarray data results.

An Integrated Systems Biology Approach Identifies the Proteasome as A Critical Host Machinery for ZIKV and DENV Replication.

The Zika virus (ZIKV) and dengue virus (DENV) flaviviruses exhibit similar replicative processes but have distinct clinical outcomes. A systematic understanding of virus-host protein-protein interaction networks can reveal cellular pathways critical to viral replication and disease pathogenesis. Here we employed three independent systems biology approaches toward this goal. First, protein array analysis of direct interactions between individual ZIKV/DENV viral proteins and 20,240 human proteins revealed multiple conserved cellular pathways and protein complexes, including proteasome complexes. Second, an RNAi screen of 10,415 druggable genes identified the host proteins required for ZIKV infection and uncovered that proteasome proteins were crucial in this process. Third, high-throughput screening of 6016 bioactive compounds for ZIKV inhibition yielded 134 effective compounds, including six proteasome inhibitors that suppress both ZIKV and DENV replication. Integrative analyses of these orthogonal datasets pinpoint proteasomes as critical host machinery for ZIKV/DENV replication. Our study provides multi-omics datasets for further studies of flavivirus-host interactions, disease pathogenesis, and new drug targets.

Global Identification of Small Ubiquitin-related Modifier (SUMO) Substrates Reveals Crosstalk between SUMOylation and Phosphorylation Promotes Cell Migration.

Proteomics studies have revealed that SUMOylation is a widely used post-translational modification (PTM) in eukaryotes. However, how SUMO E1/2/3 complexes use different SUMO isoforms and recognize substrates remains largely unknown. Using a human proteome microarray-based activity screen, we identified over 2500 proteins that undergo SUMO E3-dependent SUMOylation. We next constructed a SUMO isoform- and E3 ligase-dependent enzyme-substrate relationship network. Protein kinases were significantly enriched among SUMOylation substrates, suggesting crosstalk between phosphorylation and SUMOylation. Cell-based analyses of tyrosine kinase, PYK2, revealed that SUMOylation at four lysine residues promoted PYK2 autophosphorylation at tyrosine 402, which in turn enhanced its interaction with SRC and full activation of the SRC-PYK2 complex. SUMOylation on WT but not the 4KR mutant of PYK2 further elevated phosphorylation of the downstream components in the focal adhesion pathway, such as paxillin and Erk1/2, leading to significantly enhanced cell migration during wound healing. These studies illustrate how our SUMO E3 ligase-substrate network can be used to explore crosstalk between SUMOylation and other PTMs in many biological processes.

Multiplexed Biomarker Panels Discriminate Zika and Dengue Virus Infection in Humans.

Zika virus (ZIKV) and dengue virus (DENV) are closely related flaviviruses that cause widespread, acute febrile illnesses, notably microcephaly for fetuses of infected pregnant women. Detecting the viral cause of these illnesses is paramount to determine risks to patients, counsel pregnant women, and help fight outbreaks. A combined diagnostic algorithm for ZIKV and DENV requires Reverse transcription polymerase chain reaction (RT-PCR) and IgM antibody detection. RT-PCR differentiates between DENV and ZIKV infections during the acute phases of infection, but differentiation based on IgM antibodies is currently nearly impossible in endemic areas. We have developed a ZIKV/DENV protein array and tested it with serum samples collected from ZIKV- and DENV-infected patients and healthy subjects in Puerto Rico. Our analyses reveal a biomarker panel that are capable of discriminating ZIKV and DENV infections with high accuracy, including Capsid protein from African ZIKV strain MR766, and other 5 pair of proteins (NS1, NS2A, NS3, NS4B and NS5) from ZIKV and DENV respectively. Both sensitivity and specificity of the test for ZIKV from DENV are around 90%. We propose that the ZIKV/DENV protein array will be used in future studies to discriminate patients infected with ZIKV from DENV.

Zika-Virus-Encoded NS2A Disrupts Mammalian Cortical Neurogenesis by Degrading Adherens Junction Proteins.

Zika virus (ZIKV) directly infects neural progenitors and impairs their proliferation. How ZIKV interacts with the host molecular machinery to impact neurogenesis in vivo is not well understood. Here, by systematically introducing individual proteins encoded by ZIKV into the embryonic mouse cortex, we show that expression of ZIKV-NS2A, but not Dengue virus (DENV)-NS2A, leads to reduced proliferation and premature differentiation of radial glial cells and aberrant positioning of newborn neurons. Mechanistically, in vitro mapping of protein-interactomes and biochemical analysis suggest interactions between ZIKA-NS2A and multiple adherens junction complex (AJ) components. Functionally, ZIKV-NS2A, but not DENV-NS2A, destabilizes the AJ complex, resulting in impaired AJ formation and aberrant radial glial fiber scaffolding in the embryonic mouse cortex. Similarly, ZIKA-NS2A, but not DENV-NS2A, reduces radial glial cell proliferation and causes AJ deficits in human forebrain organoids. Together, our results reveal pathogenic mechanisms underlying ZIKV infection in the developing mammalian brain.

A Human Proteome Array Approach to Identifying Key Host Proteins Targeted by Toxoplasma Kinase ROP18.

Toxoplasma kinase ROP18 is a key molecule responsible for the virulence of Toxoplasma gondii; however, the mechanisms by which ROP18 exerts parasite virulence via interaction with host proteins remain limited to a small number of identified substrates. To identify a broader array of ROP18 substrates, we successfully purified bioactive mature ROP18 and used it to probe a human proteome array. Sixty eight new putative host targets were identified. Functional annotation analysis suggested that these proteins have a variety of functions, including metabolic process, kinase activity and phosphorylation, cell growth, apoptosis and cell death, and immunity, indicating a pleiotropic role of ROP18 kinase. Among these proteins, four candidates, p53, p38, UBE2N, and Smad1, were further validated. We demonstrated that ROP18 targets p53, p38, UBE2N, and Smad1 for degradation. Importantly, we demonstrated that ROP18 phosphorylates Smad1 Ser-187 to trigger its proteasome-dependent degradation. Further functional characterization of the substrates of ROP18 may enhance understanding of the pathogenesis of Toxoplasma infection and provide new therapeutic targets. Similar strategies could be used to identify novel host targets for other microbial kinases functioning at the pathogen-host interface.

Global analysis of phosphorylation networks in humans.

Phosphorylation-mediated signaling plays a crucial role in nearly every aspect of cellular physiology. A recent study based on protein microarray experiments identified a large number of kinase-substrate relationships (KSRs), and built a comprehensive and reliable phosphorylation network in humans. Analysis of this network, in conjunction with additional resources, revealed several key features. First, comparison of the human and yeast phosphorylation networks uncovered an evolutionarily conserved signaling backbone dominated by kinase-to-kinase relationships. Second, although most of the KSRs themselves are not conserved, the functions enriched in the substrates for a given kinase are often conserved. Third, the prevalence of kinase-transcription factor regulatory modules suggests that phosphorylation and transcriptional regulatory networks are inherently wired together to form integrated regulatory circuits. Overall, the phosphorylation networks described in this work promise to offer new insights into the properties of kinase signaling pathways, at both the global and the protein levels. This article is part of a Special Issue entitled: Computational Proteomics, Systems Biology & Clinical Implications. Guest Editor: Yudong Cai.

PhosphoNetworks: a database for human phosphorylation networks.

SUMMARY: Phosphorylation plays an important role in cellular signal transduction. Current phosphorylation-related databases often focus on the phosphorylation sites, which are mainly determined by mass spectrometry. Here, we present PhosphoNetworks, a phosphorylation database built on a high-resolution map of phosphorylation networks. This high-resolution map of phosphorylation networks provides not only the kinase-substrate relationships (KSRs), but also the specific phosphorylation sites on which the kinases act on the substrates. The database contains the most comprehensive dataset for KSRs, including the relationships from a recent high-throughput project for identification of KSRs using protein microarrays, as well as known KSRs curated from the literature. In addition, the database also includes several analytical tools for dissecting phosphorylation networks. PhosphoNetworks is expected to play a prominent role in proteomics and phosphorylation-related disease research. AVAILABILITY AND IMPLEMENTATION: http://www.phosphonetworks.org

DNA methylation presents distinct binding sites for human transcription factors.

DNA methylation, especially CpG methylation at promoter regions, has been generally considered as a potent epigenetic modification that prohibits transcription factor (TF) recruitment, resulting in transcription suppression. Here, we used a protein microarray-based approach to systematically survey the entire human TF family and found numerous purified TFs with methylated CpG (mCpG)-dependent DNA-binding activities. Interestingly, some TFs exhibit specific binding activity to methylated and unmethylated DNA motifs of distinct sequences. To elucidate the underlying mechanism, we focused on Kruppel-like factor 4 (KLF4), and decoupled its mCpG- and CpG-binding activities via site-directed mutagenesis. Furthermore, KLF4 binds specific methylated or unmethylated motifs in human embryonic stem cells in vivo. Our study suggests that mCpG-dependent TF binding activity is a widespread phenomenon and provides a new framework to understand the role and mechanism of TFs in epigenetic regulation of gene transcription. DOI:http://dx.doi.org/10.7554/eLife.00726.001.

Global expression profiling of transcription factor genes provides new insights into pathogenicity and stress responses in the rice blast fungus.

Because most efforts to understand the molecular mechanisms underpinning fungal pathogenicity have focused on studying the function and role of individual genes, relatively little is known about how transcriptional machineries globally regulate and coordinate the expression of a large group of genes involved in pathogenesis. Using quantitative real-time PCR, we analyzed the expression patterns of 206 transcription factor (TF) genes in the rice blast fungus Magnaporthe oryzae under 32 conditions, including multiple infection-related developmental stages and various abiotic stresses. The resulting data, which are publicly available via an online platform, provided new insights into how these TFs are regulated and potentially work together to control cellular responses to a diverse array of stimuli. High degrees of differential TF expression were observed under the conditions tested. More than 50% of the 206 TF genes were up-regulated during conidiation and/or in conidia. Mutations in ten conidiation-specific TF genes caused defects in conidiation. Expression patterns in planta were similar to those under oxidative stress conditions. Mutants of in planta inducible genes not only exhibited sensitive to oxidative stress but also failed to infect rice. These experimental validations clearly demonstrated the value of TF expression patterns in predicting the function of individual TF genes. The regulatory network of TF genes revealed by this study provides a solid foundation for elucidating how M. oryzae regulates its pathogenesis, development, and stress responses.

Construction of human activity-based phosphorylation networks.

The landscape of human phosphorylation networks has not been systematically explored, representing vast, unchartered territories within cellular signaling networks. Although a large number of in vivo phosphorylated residues have been identified by mass spectrometry (MS)-based approaches, assigning the upstream kinases to these residues requires biochemical analysis of kinase-substrate relationships (KSRs). Here, we developed a new strategy, called CEASAR, based on functional protein microarrays and bioinformatics to experimentally identify substrates for 289 unique kinases, resulting in 3656 high-quality KSRs. We then generated consensus phosphorylation motifs for each of the kinases and integrated this information, along with information about in vivo phosphorylation sites determined by MS, to construct a high-resolution map of phosphorylation networks that connects 230 kinases to 2591 in vivo phosphorylation sites in 652 substrates. The value of this data set is demonstrated through the discovery of a new role for PKA downstream of Btk (Bruton's tyrosine kinase) during B-cell receptor signaling. Overall, these studies provide global insights into kinase-mediated signaling pathways and promise to advance our understanding of cellular signaling processes in humans.

Phosphorylation of the chromatin binding domain of KSHV LANA.

The Kaposi sarcoma associated herpesvirus (KSHV) latency associated nuclear antigen (LANA) is expressed in all KSHV associated malignancies and is essential for maintenance of KSHV genomes in infected cells. To identify kinases that are potentially capable of modifying LANA, in vitro phosphorylation assays were performed using an Epstein Barr virus plus LANA protein microarray and 268 human kinases purified in active form from yeast. Interestingly, of the Epstein-Barr virus proteins on the array, the EBNA1 protein had the most similar kinase profile to LANA. We focused on nuclear kinases and on the N-terminus of LANA (amino acids 1-329) that contains the LANA chromatin binding domain. Sixty-three nuclear kinases phosphorylated the LANA N-terminus. Twenty-four nuclear kinases phosphorylated a peptide covering the LANA chromatin binding domain (amino acids 3-21). Alanine mutations of serine 10 and threonine 14 abolish or severely diminish chromatin and histone binding by LANA. However, conversion of these residues to the phosphomimetic glutamic acid restored histone binding suggesting that phosphorylation of serine 10 and threonine 14 may modulate LANA function. Serine 10 and threonine 14 were validated as substrates of casein kinase 1, PIM1, GSK-3 and RSK3 kinases. Short-term treatment of transfected cells with inhibitors of these kinases found that only RSK inhibition reduced LANA interaction with endogenous histone H2B. Extended treatment of PEL cell cultures with RSK inhibitor caused a decrease in LANA protein levels associated with p21 induction and a loss of PEL cell viability. The data indicate that RSK phosphorylation affects both LANA accumulation and function.

Rapid identification of monospecific monoclonal antibodies using a human proteome microarray.

To broaden the range of tools available for proteomic research, we generated a library of 16,368 unique full-length human ORFs that are expressible as N-terminal GST-His(6) fusion proteins. Following expression in yeast, these proteins were then individually purified and used to construct a human proteome microarray. To demonstrate the usefulness of this reagent, we developed a streamlined strategy for the production of monospecific monoclonal antibodies that used immunization with live human cells and microarray-based analysis of antibody specificity as its central components. We showed that microarray-based analysis of antibody specificity can be performed efficiently using a two-dimensional pooling strategy. We also demonstrated that our immunization and selection strategies result in a large fraction of monospecific monoclonal antibodies that are both immunoblot and immunoprecipitation grade. Our data indicate that the pipeline provides a robust platform for the generation of monoclonal antibodies of exceptional specificity.

Regulation of CK2 by phosphorylation and O-GlcNAcylation revealed by semisynthesis.

Protein serine-threonine kinase casein kinase II (CK2) is involved in a myriad of cellular processes including cell growth and proliferation through its phosphorylation of hundreds of substrates, yet how CK2 function is regulated is poorly understood. Here we report that the CK2 catalytic subunit CK2α is modified by O-linked ß-N-acetyl-glucosamine (O-GlcNAc) on Ser347, proximal to a cyclin-dependent kinase phosphorylation site (Thr344). We use protein semisynthesis to show that phosphorylation of Thr344 increases the cellular stability of CK2α by strengthening its interaction with Pin1, whereas glycosylation of Ser347 seems to be antagonistic to Thr344 phosphorylation and permissive to proteasomal degradation. By performing kinase assays with site-specifically phospho- and glyco-modified CK2α in combination with CK2ß and Pin1 binding partners on human protein microarrays, we show that the kinase substrate selectivity of CK2 is modulated by these specific post-translational modifications. This study suggests how a promiscuous protein kinase can be regulated at multiple levels to achieve particular biological outputs.

A functional protein microarray approach to characterizing posttranslational modifications on lysine residues.

Functional protein microarrays offer a versatile platform to address diverse biological questions. Printing individually purified proteins in a spatially addressable format makes it straightforward to investigating binary interactions. To connect substrates to their upstream modifying enzymes, such as kinases, ubiqutin (Ub) ligases, SUMOylation E3 ligases, and acetyltransferases, is an especially daunting task using traditional methodologies. In recent years, regulation via various types of posttranslational modifications (PTMs) on lysine residues is emerging as an important mechanism(s) underlining diverse biological -processes. Our group has been developing and applying functional protein microarrays constructed for different model organisms to globally identify enzyme-substrate interactions with a focus on lysine PTMs. In particular, we have characterized the pleiotropic functions of a ubiquitin E3 ligase, Rsp5, via identification of its downstream substrates using a yeast proteome chip. Also, we have identified nonhistone substrates of the acetyltransferase NuA4 complex in yeast, and revealed that reversible acetylation on a metabolic enzyme affects a glucose metabolism and contributes to life span. In this chapter, we will provide detailed protocols for the investigation of ubiquitylation and acetylation. These protocols are generally applicable for different organisms.

Differential roles of the phospholipase C genes in fungal development and pathogenicity of Magnaporthe oryzae.

Calcium plays a critical role in a variety of cellular processes in cells. However, relatively little is known about the biological effects of Ca²+ signaling on morphogenesis and pathogenesis in the rice blast fungus Magnaporthe oryzae compared to other signaling pathways. We have previously demonstrated that MoPLC1-mediated calcium regulation is important for infection-related development and pathogenicity in M. oryzae. In the present study, four genes encoding phospholipase C (PLC) isozymes (MoPLC2 to MoPLC5), which differ from MoPLC1 in their domain organization, were additionally identified. The C2 domain involved in Ca²+-dependent membrane binding is found only in MoPLC2 and MoPLC3. Detailed functional analysis using deletion mutants for MoPLC2 and MoPLC3 indicated that MoPLC2 and MoPLC3 play essential roles in development. The two deletion mutants for MoPLC2 and MoPLC3 showed reduced conidiation and a defect in appressorium-mediated penetration. Reintroduction of the genes restored defects of ΔMoplc2 and ΔMoplc3. Notably, ΔMoplc2 and ΔMoplc3 mutants developed multiple appressoria on separate germ tubes of a conidium, indicating that MoPLC2- and MoPLC3-regulated signaling suppresses a feedback loop of a pathway for appressorial development. The similarity in phenotypic defects between the two mutants indicates that both MoPLC2 and MoPLC3 are important for regulation of appropriate levels of signaling molecules in a similar manner. Comparative analysis indicated that the two MoPLCs-mediated signaling pathways have interrelated, but distinct, roles in the development of M. oryzae.

Homeobox transcription factors are required for conidiation and appressorium development in the rice blast fungus Magnaporthe oryzae.

The appropriate development of conidia and appressoria is critical in the disease cycle of many fungal pathogens, including Magnaporthe oryzae. A total of eight genes (MoHOX1 to MoHOX8) encoding putative homeobox transcription factors (TFs) were identified from the M. oryzae genome. Knockout mutants for each MoHOX gene were obtained via homology-dependent gene replacement. Two mutants, DeltaMohox3 and DeltaMohox5, exhibited no difference to wild-type in growth, conidiation, conidium size, conidial germination, appressorium formation, and pathogenicity. However, the DeltaMohox1 showed a dramatic reduction in hyphal growth and increase in melanin pigmentation, compared to those in wild-type. DeltaMohox4 and DeltaMohox6 showed significantly reduced conidium size and hyphal growth, respectively. DeltaMohox8 formed normal appressoria, but failed in pathogenicity, probably due to defects in the development of penetration peg and invasive growth. It is most notable that asexual reproduction was completely abolished in DeltaMohox2, in which no conidia formed. DeltaMohox2 was still pathogenic through hypha-driven appressoria in a manner similar to that of the wild-type. However, DeltaMohox7 was unable to form appressoria either on conidial germ tubes, or at hyphal tips, being non-pathogenic. These factors indicate that M. oryzae is able to cause foliar disease via hyphal appressorium-mediated penetration, and MoHOX7 is mutually required to drive appressorium formation from hyphae and germ tubes. Transcriptional analyses suggest that the functioning of M. oryzae homeobox TFs is mediated through the regulation of gene expression and is affected by cAMP and Ca(2+) signaling and/or MAPK pathways. The divergent roles of this gene set may help reveal how the genome and regulatory pathways evolved within the rice blast pathogen and close relatives.

Profiling the human protein-DNA interactome reveals ERK2 as a transcriptional repressor of interferon signaling.

Protein-DNA interactions (PDIs) mediate a broad range of functions essential for cellular differentiation, function, and survival. However, it is still a daunting task to comprehensively identify and profile sequence-specific PDIs in complex genomes. Here, we have used a combined bioinformatics and protein microarray-based strategy to systematically characterize the human protein-DNA interactome. We identified 17,718 PDIs between 460 DNA motifs predicted to regulate transcription and 4,191 human proteins of various functional classes. Among them, we recovered many known PDIs for transcription factors (TFs). We identified a large number of unanticipated PDIs for known TFs, as well as for previously uncharacterized TFs. We also found that over three hundred unconventional DNA-binding proteins (uDBPs)--which include RNA-binding proteins, mitochondrial proteins, and protein kinases--showed sequence-specific PDIs. One such uDBP, ERK2, acts as a transcriptional repressor for interferon gamma-induced genes, suggesting important biological roles for such proteins.

Phospholipase C-mediated calcium signalling is required for fungal development and pathogenicity in Magnaporthe oryzae.

Calcium signalling has profound implications in the fungal infection of plants and animals, during which a series of physiological and morphological transitions are required. In this article, using a model fungal pathogen, Magnaporthe oryzae, we demonstrate that the regulation of the intracellular calcium concentration ([Ca(2+)](int)) is essential for fungal development and pathogenesis. Imaging of [Ca(2+)](int) showed that infection-specific morphogenesis is highly correlated with the spatiotemporal regulation of calcium flux. Deletion of the fungal phospholipase C gene (M. oryzae phospholipase C 1, MoPLC1) suppressed calcium flux, resulting in a fungus defective in developmental steps, including appressorium formation and pathogenicity. Surprisingly, the PLC-delta1 gene of mouse was able to functionally substitute for MoPLC1 by restoring the calcium flux, suggesting the evolutionary conservation of the phospholipase C-mediated regulation of calcium flux. Our results reveal that MoPLC1 is a conserved modulator of calcium flux that is essential for the regulation of key steps in fungal development and pathogenesis.