Prediction and Analysis of Post-Translational Pyruvoyl Residue Modification Sites from Internal Serines in Proteins

Most of pyruvoyl-dependent proteins observed in prokaryotes and eukaryotes are critical regulatory enzymes, which are primary targets of inhibitors for anti-cancer and anti-parasitic therapy. These proteins undergo an autocatalytic, intramolecular self-cleavage reaction in which a covalently bound pyruvoyl group is generated on a conserved serine residue. Traditional detections of the modified serine sites are performed by experimental approaches, which are often labor-intensive and time-consuming. In this study, we initiated in an attempt for the computational predictions of such serine sites with Feature Selection based on a Random Forest. Since only a small number of experimentally verified pyruvoyl-modified proteins are collected in the protein database at its current version, we only used a small dataset in this study. After removing proteins with sequence identities >60%, a non-redundant dataset was generated and was used, which contained only 46 proteins, with one pyruvoyl serine site for each protein. Several types of features were considered in our method including PSSM conservation scores, disorders, secondary structures, solvent accessibilities, amino acid factors and amino acid occurrence frequencies. As a result, a pretty good performance was achieved in our dataset. The best 100.00% accuracy and 1.0000 MCC value were obtained from the training dataset, and 93.75% accuracy and 0.8441 MCC value from the testing dataset. The optimal feature set contained 9 features. Analysis of the optimal feature set indicated the important roles of some specific features in determining the pyruvoyl-group-serine sites, which were consistent with several results of earlier experimental studies. These selected features may shed some light on the in-depth understanding of the mechanism of the post-translational self-maturation process, providing guidelines for experimental validation. Future work should be made as more pyruvoyl-modified proteins are found and the method should be evaluated on larger datasets. At last, the predicting software can be downloaded from http://www.nkbiox.com/sub/pyrupred/index.html.     

蛋白酶体相关翻译后修饰的研究进展

蛋白酶体是具有多种蛋白水解酶活性的蛋白质降解系统,由于细胞内许多关键信号调控分子都是蛋白酶体的降解底物,因此蛋白酶体在细胞周期调控、基因表达、炎症反应等各种关键的生物学活动中都发挥着极其重要的调节作用。蛋白酶体的降解活性也同时受多种机制的调控,其中的翻译后修饰是蛋白酶体降解途径中一个不可忽视的方面。着重阐述蛋白酶体自身及其底物的几种重要的翻译后修饰,探讨最新的进展及其生物学意义。     

Computational Analysis of the Predicted Evolutionary Conservation of Human Phosphorylation Sites

Protein kinase-mediated phosphorylation is among the most important post-translational modifications. However, few phosphorylation sites have been experimentally identified for most species, making it difficult to determine the degree to which phosphorylation sites are conserved. The goal of this study was to use computational methods to characterize the conservation of human phosphorylation sites in a wide variety of eukaryotes. Using experimentally-determined human sites as input, homologous phosphorylation sites were predicted in all 432 eukaryotes for which complete proteomes were available. For each pair of species, we calculated phosphorylation site conservation as the number of phosphorylation sites found in both species divided by the number found in at least one of the two species. A clustering of the species based on this conservation measure was concordant with phylogenies based on traditional genomic measures. For a subset of the 432 species, phosphorylation site conservation was compared to conservation of both protein kinases and proteins in general. Protein kinases exhibited the highest degree of conservation, while general proteins were less conserved and phosphorylation sites were least conserved. Although preliminary, these data tentatively suggest that variation in phosphorylation sites may play a larger role in explaining phenotypic differences among organisms than differences in the complements of protein kinases or general proteins.     

Evolutionary Dynamics of Variant Genomes of Human Papillomavirus Types 18, 45, and 97

Human papillomavirus type 18 (HPV18) and HPV45 account for approximately 20% of all cervix cancers. We show that HPV18, HPV45, and the recently discovered HPV97 comprise a clade sharing a most recent common ancestor within HPV α7 species. Variant lineages of these HPV types were classified by sequence analysis of the upstream regulatory region/E6 region among cervical samples from a population-based study in Costa Rica, and 27 representative genomes from each major variant lineage were sequenced. Nucleotide variation within HPV18 and HPV45 was 3.82% and 2.39%, respectively, and amino acid variation was 4.73% and 2.87%, respectively. Only 18 nucleotide variations, of which 10 were nonsynonymous, were identified among three HPV97 genomes. Full-genome comparisons revealed maximal diversity between HPV18 African and non-African variants (2.6% dissimilarity), whereas HPV18 Asian-American [E1 (AA)] and European (E2) variants were closely related (less than 0.5% dissimilarity); HPV45 genomes had a maximal difference of 1.6% nucleotides. Using a Bayesian Markov chain Monte Carlo (MCMC) method, the divergence times of HPV18, -45, and -97 from their most recent common ancestors indicated that HPV18 diverged approximately 7.7 million years (Myr) ago, whereas HPV45 and HPV97 split off around 5.7 Myr ago, in a period encompassing the divergence of the great ape species. Variants within the HPV18/45/97 lineages were estimated to have diverged from their common ancestors in the genus Homo within the last 1 Myr (<0.7 Myr). To investigate the molecular basis of HPV18, HPV45, and HPV97 evolution, regression models of codon substitution were used to identify lineages and amino acid sites under selective pressure. The E5 open reading frame (ORF) of HPV18 and the E4 ORFs of HPV18, HPV45, and HPV18/45/97 had nonsynonymous/synonymous substitution rate ratios (d_N/d_S) over 1 indicative of positive Darwinian selection. The L1 ORF of HPV18 genomes had an increased proportion of nonsynonymous substitutions (4.93%; average d_N/d_S ratio [M3] = 0.3356) compared to HPV45 (1.86%; M3 = 0.1268) and HPV16 (2.26%; M3 = 0.1330) L1 ORFs. In contrast, HPV18 and HPV16 genomes had similar amino acid substitution rates within the E1 ORF (2.89% and 3.24%, respectively), while HPV45 E1 was highly conserved (amino acid substitution rate was 0.77%). These data provide an evolutionary history of this medically important clade of HPVs and identify an unexpected divergence of the L1 gene of HPV18 that may have clinical implications for the long-term use of an L1-virus-like particle-based prophylactic vaccine.Papillomaviruses (PVs) are a large family of related viruses with circular double-stranded DNA genomes 8 kb in size. Some PV types cause epithelial hyperplasias ranging from benign exophytic warts to premalignant lesions that can progress to invasive cancer. Among the 61 currently designated alpha human PVs (HPVs), the majority have been isolated from the mucosal surface of the genital or oral region (8, 14). Of these, a select group have oncogenic potential and are associated with cervical cancer (11). Specifically, HPV type 16 (HPV16) and HPV18 have been identified in approximately two-thirds of cervical cancers, this tumor is the second most common cancer in women, and it is the principal cancer of women in developing countries (5, 24, 25, 30, 37).To date, studies of HPV18 variants have identified three lineages corresponding to the continental locations where the viral samples were obtained: European (E), Asian-American (AA), and African (Af) (29). The phylogeny of HPV18 variants is reflective of the migration patterns of Homo sapiens and suggests that HPV18 variant lineages might have diverged through genetic isolation at approximately the same time as Homo sapiens began establishing residence in different continental regions. Previous HPV18 intratypic phylogenetic analyses were limited to partial regions of the genome (3, 7, 29). Nevertheless, studies also suggest that HPV18 variants are associated with different levels of oncogenic potential and persistence and histological tumor types (1, 6, 35, 36, 46).HPV45 and HPV97 are the viral types most closely related to HPV18 and taken together form a clade and share a most recent common ancestor (MRCA). HPV97 is a recently described rare type (8, 17). HPV18 and HPV45 account for approximately 20% of all cervix cancers (25). Although HPV45 is a common type found in cervical cancer, its evolutionary history and sequence variability have not been extensively studied.In this report, 27 complete genomes representing the major variant lineages of HPV18, HPV45, and HPV97 were cloned and/or sequenced from clinical samples. Based on full genomes, the intratype/intertype evolutionary trees of HPV18, HPV45, and HPV97 were constructed. By examining the rate ratio of nonsynonymous (d_N) to synonymous (d_S) substitutions per site, diversifying selection acting on each of the eight protein-encoding regions of HPV18, HPV45, and HPV97 was evaluated. In addition, the times of divergence of HPV18/45/97 variants from their MRCA were investigated. These data provide an evolutionary history of this medically important clade of HPVs.     

硫巯基化：一种新的蛋白质翻译后修饰

随着对硫化氢（hydrogen sulfide,H2S）生理效应的研究,蛋白质硫巯基化（S-sulfhydration）修饰已进入人们的视野。已知依赖于H2S的蛋白质硫巯基化是继磷酸化（phosphorylation）、泛素化（ubiquitylation）、乙酰化（acetylation）和S-亚硝基化（S-nitrosylation）等之后的一种新的蛋白质翻译后修饰方式。对动物的研究表明,蛋白质硫巯基化修饰通过影响蛋白质活性和功能,从而在细胞内信号通路中发挥重要的调控作用。最近的研究结果提示,硫巯基化修饰还参与调节植物新陈代谢和形态建成。本文阐述了依赖于H2S的蛋白质硫巯基化的作用机制、检测方法和生理功能,并提出硫巯基化修饰也可能参与植物细胞信号转导的观点。     

Conservation and Periodicity of DNA Bend Sites in Eukaryotic Genomes

DNA bend sites appear every 680 bp on average in the human -and ß-globin gene regions. Although most of theirmolecular nature has not been unraveled, a potential bend coresequence A₂N₈A₂N₈A₂ (A/A/A) and its complementary T₂N₈T₂N₈T₂(T/T/T) appeared preferentially either in or very close to mostof the bend sites, whereas other combinations of A₂ and T₂ dinucleotides,A/T/T + A/A/T, T/T/A + T/A/A and A/T/A + T/A/T, did not. Thedistances between any two of the core sequences in the entireß-globin locus showed a strong bias to a length of701–800 bp and multiples thereof, suggesting that thereis periodicity throughout the locus. This bias was not foundfor other combinations of A₂ and T₂. Again, this periodicitywas identified in many eukaryotic genes, whereas the tendencywas absent in mRNAs and prokaryotic as well as viral genomes.     

Post-Translational Tubulin Modifications in Human Astrocyte Cultures

The cytoskeletal protein tubulin plays an integral role in the functional specialization of many cell types. In the central nervous system, post-translational modifications and the expression of specific tubulin isotypes in neurons have been analyzed in greater detail than in their astrocytic counterparts. In this study, we characterized post-translational specifications of tubulin in human astrocytes using the normal human astrocyte (NHA; Lonza) commercial cell line of fetal origin. Immunocytochemical techniques were implemented in conjunction with confocal microscopy to image class III β-tubulin (βIII-tubulin), acetylated tubulin, and polyglutamylated tubulin using fluorescent antibody probes. Fluorescent probe intensity differences and colocalization were quantitatively assessed with the ‘EBImage’ package for the statistical programming language R. Colocalization analysis revealed that, although both acetylated tubulin and polyglutamylated tubulin showed a high degree of correlation with βIII-tubulin, the correlation with acetylated tubulin was stronger. Quantification and statistical analysis of fluorescence intensity demonstrated that the fluorescence probe intensity ratio for acetylated tubulin/βIII-tubulin was greater than the ratio for polyglutamylated tubulin/βIII-tubulin. The open source GEODATA set GSE819950, comprising RNA sequencing data for the NHA cell line, was mined for the expression of enzymes responsible for tubulin modifications. Our analysis uncovered greater expression at the mRNA level for enzymes reported to function in acetylation and deacetylation as compared to enzymes implicated in glutamylation and deglutamylation. Taken together, the results represent a step toward unraveling the tubulin isotypic expression profile and post-translational modification patterns in astrocytes during human brain development.     

Comparisons of Two Large Phaeoviral Genomes and Evolutionary Implications

The evolution of viral genomes has recently attracted considerable attention. We compare the sequences of two large viral genomes, EsV-1 and FirrV-1, belonging to the family of phaeoviruses which infect different species of marine brown algae. Although their genomes differ substantially in size, these viruses share similar morphologies and similar latent infection cycles. In fact, sequence comparisons show that the viruses have more than 60% of their genes in common. However, the order of genes is completely different in the two genomes, suggesting that extensive recombinational events in addition to several large deletions had occurred during the separate evolutionary routes from a common ancestor. We investigated genes encoding components of signal transduction pathways and genes encoding replicative functions in more detail. We found that the two genomes possess different, although overlapping, sets of genes in both classes, suggesting that different genes from each class were lost, perhaps randomly, after the separate evolution from an ancestral genome. Random loss would also account for the fact that more than one-third of the genes in one viral genome has no counterparts in the other genome. We speculate that the ancestral genome belonged to a cellular organism that had once invaded a primordial brown algal host.     

Open Modification Searching of SARS-CoV-2–Human Protein Interaction Data Reveals Novel Viral Modification Sites

The outbreak of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the coronavirus 2019 disease, has led to an ongoing global pandemic since 2019. Mass spectrometry can be used to understand the molecular mechanisms of viral infection by SARS-CoV-2, for example, by determining virus–host protein–protein interactions through which SARS-CoV-2 hijacks its human hosts during infection, and to study the role of post-translational modifications. We have reanalyzed public affinity purification–mass spectrometry data using open modification searching to investigate the presence of post-translational modifications in the context of the SARS-CoV-2 virus–host protein–protein interaction network. Based on an over twofold increase in identified spectra, our detected protein interactions show a high overlap with independent mass spectrometry-based SARS-CoV-2 studies and virus–host interactions for alternative viruses, as well as previously unknown protein interactions. In addition, we identified several novel modification sites on SARS-CoV-2 proteins that we investigated in relation to their interactions with host proteins. A detailed analysis of relevant modifications, including phosphorylation, ubiquitination, and S-nitrosylation, provides important hypotheses about the functional role of these modifications during viral infection by SARS-CoV-2.     

Post-Translational Modification of Xanthine Dehydrogenase in a Natural Population of DROSOPHILA MELANOGASTER

Second chromosomes of D. melanogaster were isolated from a single natural population, and 40 were analyzed by gel-sieving electrophoresis for the presence of polymorphic loci on chromosome 2 that act to modify xanthine dehydrogenase and/or aldehyde oxidase, whose structural genes map to chromosome 3. Clear evidence of polymorphism for one or more xanthine dehydrogenase modifier loci was obtained.     

Pathway Length and Evolutionary Constraint in Amino Acid Biosynthesis

The evolutionary properties of a metabolic network may be determined by the topology of the network. One attribute of pathways that make up the network is the number of enzymatic steps between initial substrates and final products. To determine the effect of pathway length on evolutionary lability of pathway structure, we examined amino acid biosynthetic pathways across 48 sequenced organisms. We demonstrate that longer pathways exhibit lower rates of change in pathway structure than shorter pathways. This finding suggests that increasing complexity may increase constraint on evolutionary change. (Matthew T. Rutter and Rebecca A. Zufall) Both authors contributed equally to this work.