Signatures of natural selection on mutations of residues with multiple posttranslational modifications

Evolutionary Signals in Coronaviral Structural Proteins Suggest Possible Complex Mechanisms of Post-Translational Regulation in SARS-CoV-2 Virus

Post-translational regulation of proteins has emerged as a central topic of research in the field of functional proteomics. Post-translational modifications (PTMs) dynamically control the activities of proteins and are involved in a wide range of biological processes. Crosstalk between different types of PTMs represents a key mechanism of regulation and signaling. Due to the current pandemic of the novel and dangerous SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) virus, here we present an in silico analysis of different types of PTMs in structural proteins of coronaviruses. A dataset of PTM sites was studied at three levels: conservation analysis, mutational analysis and crosstalk analysis. We identified two sets of PTMs which could have important functional roles in the regulation of the structural proteins of coronaviruses. Additionally, we found seven interesting signals of potential crosstalk events. These results reveal a higher level of complexity in the mechanisms of post-translational regulation of coronaviral proteins and provide new insights into the adaptation process of the SARS-CoV-2 virus.

mUSP: a high-accuracy map of the in situ crosstalk of ubiquitylation and SUMOylation proteome predicted via the feature enhancement approach

Reversible post-translational modification (PTM) orchestrates various biological processes by changing the properties of proteins. Since many proteins are multiply modified by PTMs, identification of PTM crosstalk site has emerged to be an intriguing topic and attracted much attention. In this study, we systematically deciphered the in situ crosstalk of ubiquitylation and SUMOylation that co-occurs on the same lysine residue. We first collected 3363 ubiquitylation-SUMOylation (UBS) crosstalk site on 1302 proteins and then investigated the prime sequence motifs, the local evolutionary degree and the distribution of structural annotations at the residue and sequence levels between the UBS crosstalk and the single modification sites. Given the properties of UBS crosstalk sites, we thus developed the mUSP classifier to predict UBS crosstalk site by integrating different types of features with two-step feature optimization by recursive feature elimination approach. By using various cross-validations, the mUSP model achieved an average area under the curve (AUC) value of 0.8416, indicating its promising accuracy and robustness. By comparison, the mUSP has significantly better performance with the improvement of 38.41 and 51.48% AUC values compared to the cross-results by the previous single predictor. The mUSP was implemented as a web server available at http://bioinfo.ncu.edu.cn/mUSP/index.html to facilitate the query of our high-accuracy UBS crosstalk results for experimental design and validation.

New insights into the heat responses of grape leaves via combined phosphoproteomic and acetylproteomic analyses

Heat stress is a serious and widespread threat to the quality and yield of many crop species, including grape (Vitis vinifera L.), which is cultivated worldwide. Here, we conducted phosphoproteomic and acetylproteomic analyses of leaves of grape plants cultivated under four distinct temperature regimes. The phosphorylation or acetylation of a total of 1011 phosphoproteins with 1828 phosphosites and 96 acetyl proteins with 148 acetyl sites changed when plants were grown at 35 °C, 40 °C, and 45 °C in comparison with the proteome profiles of plants grown at 25 °C. The greatest number of changes was observed at the relatively high temperatures. Functional classification and enrichment analysis indicated that phosphorylation, rather than acetylation, of serine/arginine-rich splicing factors was involved in the response to high temperatures. This finding is congruent with previous observations by which alternative splicing events occurred more frequently in grapevine under high temperature. Changes in acetylation patterns were more common than changes in phosphorylation patterns in photosynthesis-related proteins at high temperatures, while heat-shock proteins were associated more with modifications involving phosphorylation than with those involving acetylation. Nineteen proteins were identified with changes associated with both phosphorylation and acetylation, which is consistent with crosstalk between these posttranslational modification types.

g2pDB: A Database Mapping Protein Post-Translational Modifications to Genomic Coordinates

Large scale proteomics have made it possible to broadly screen samples for the presence of many types of post-translational modifications, such as phosphorylation, acetylation, and ubiquitination. This type of data has allowed the localization of these modifications to either a specific site on a proteolytically generated peptide or to within a small domain on the peptide. The resulting modification acceptor sites can then be mapped onto the appropriate protein sequences and the information archived. This paper describes the usage of a very large archive of experimental observations of human post-translational modifications to create a map of the most reproducible modification observations onto the complete set of human protein sequences. This set of modification acceptor sites was then directly translated into the genomic coordinates for the codons for the residues at those sites. We constructed the database g2pDB using this protein-to-codon site mapping information. The information in g2pDB has been made available through a RESTful-style API, allowing researchers to determine which specific protein modifications would be perturbed by a set of observed nucleotide variants determined by high throughput DNA or RNA sequencing.

A comprehensive catalog of the lysine-acetylation targets in rice (Oryza sativa) based on proteomic analyses

Lysine acetylation is a dynamic and reversible post-translational modification that plays an important role in the gene transcription regulation. Here, we report high quality proteome-scale data for lysine-acetylation (Kac) sites and Kac proteins in rice (Oryza sativa). A total of 1337 Kac sites in 716 Kac proteins with diverse biological functions and subcellular localizations were identified in rice seedlings. About 42% of the sites were predicted to be localized in the chloroplast. Seven putative acetylation motifs were detected. Phenylalanine, located in both the upstream and downstream of the Kac sites, is the most conserved amino acid surrounding the regions. In addition, protein interaction network analysis revealed that a variety of signaling pathways are modulated by protein acetylation. KEGG pathway category enrichment analysis indicated that glyoxylate and dicarboxylate metabolism, carbon metabolism, and photosynthesis pathways are significantly enriched. Our results provide an in-depth understanding of the acetylome in rice seedlings, and the method described here will facilitate the systematic study of how Kac functions in growth, development, and abiotic and biotic stress responses in rice and other plants.

BIOLOGICAL SIGNIFICANCE

Rice is one of the most important crops consumption and is a model monocot for research. In this study, we combined a highly sensitive immune-affinity purification method (used pan anti-acetyl-lysine antibody conjugated agarose for immunoaffinity acetylated peptide enrichment) with high-resolution LC-MS/MS. In total, we identified 1337 Kac sites on 716 Kac proteins in rice cells. Bioinformatic analysis of the acetylome revealed that the acetylated proteins are involved in a variety of cellular functions and have diverse subcellular localizations. We also identified seven putative acetylation motifs in the acetylated proteins of rice. In addition, protein interaction network analysis revealed that a variety of signaling pathways were modulated by protein acetylation. KEGG pathway category enrichment analysis indicated that glyoxylate and dicarboxylate metabolism, carbon metabolism, and photosynthesis pathways were significantly enriched. To our knowledge, the number of Kac sites we identified was 23-times greater and the number of Kac proteins was 16-times greater than in a previous report. Our results provide an in-depth understanding of the acetylome in rice seedlings, and the method described here will facilitate the systematic study of how Kac functions in growth, development and responses to abiotic and biotic stresses in rice or other plants.

Suberoylanilide hydroxamic acid treatment reveals crosstalks among proteome, ubiquitylome and acetylome in non-small cell lung cancer A549 cell line

Suberoylanilide hydroxamic acid (SAHA) is a well-known histone deacetylase (HDAC) inhibitor and has been used as practical therapy for breast cancer and non-small cell lung cancer (NSCLC). It is previously demonstrated that SAHA treatment could extensively change the profile of acetylome and proteome in cancer cells. However, little is known about the impact of SAHA on other protein modifications and the crosstalks among different modifications and proteome, hindering the deep understanding of SAHA-mediated cancer therapy. In this work, by using SILAC technique, antibody-based affinity enrichment and high-resolution LC-MS/MS analysis, we investigated quantitative proteome, acetylome and ubiquitylome as well as crosstalks among the three datasets in A549 cells toward SAHA treatment. In total, 2968 proteins, 1099 acetylation sites and 1012 ubiquitination sites were quantified in response to SAHA treatment, respectively. With the aid of intensive bioinformatics, we revealed that the proteome and ubiquitylome were negatively related upon SAHA treatment. Moreover, the impact of SAHA on acetylome resulted in 258 up-regulated and 99 down-regulated acetylation sites at the threshold of 1.5 folds. Finally, we identified 55 common sites with both acetylation and ubiquitination, among which ubiquitination level in 43 sites (78.2%) was positive related to acetylation level.

Systematic analysis of the in situ crosstalk of tyrosine modifications reveals no additional natural selection on multiply modified residues

Recent studies have indicated that different post-translational modifications (PTMs) synergistically orchestrate specific biological processes by crosstalks. However, the preference of the crosstalk among different PTMs and the evolutionary constraint on the PTM crosstalk need further dissections. In this study, the in situ crosstalk at the same positions among three tyrosine PTMs including sulfation, nitration and phosphorylation were systematically analyzed. The experimentally identified sulfation, nitration and phosphorylation sites were collected and integrated with reliable predictions to perform large-scale analyses of in situ crosstalks. From the results, we observed that the in situ crosstalk between sulfation and nitration is significantly under-represented, whereas both sulfation and nitration prefer to co-occupy with phosphorylation at same tyrosines. Further analyses suggested that sulfation and nitration preferentially co-occur with phosphorylation at specific positions in proteins, and participate in distinct biological processes and functions. More interestingly, the long-term evolutionary analysis indicated that multi-PTM targeting tyrosines didn't show any higher conservation than singly modified ones. Also, the analysis of human genetic variations demonstrated that there is no additional functional constraint on inherited disease, cancer or rare mutations of multiply modified tyrosines. Taken together, our systematic analyses provided a better understanding of the in situ crosstalk among PTMs.