Discrimination of Mammalian GPI-Anchored Proteins by Hydropathy and Amino Acid Propensities

A 1.7-kb DNA fragment cloned from Zymomonas mobilis genomic DNA complemented the inability to grow on sucrose of a Sue ^? mutant of Z. mobilis that was deficient in the production of both extracellular levansucrase and invertase. Analysis of the nucleotide sequence of the fragment found two open reading frames (ORFs), both of which did not correspond to the structural gene for the levansucrase or the invertase. By subcloning each ORF into two different Suc ^? mutants of Z. mobilis, it has been found that the first ORF (gene zliE) activates the production of the extracellular levansucrase and invertase, and the second ORF (gene zliS) stimulates the secretion of the two enzymes. Gene zliS might contribute to the secretion of proteins having no signal peptide. The expression of zliE and zliS seemed to be under the control of the same promoter.     

A novel method for protein secondary structure prediction using dual-layer SVM and profiles

A high-performance method was developed for protein secondary structure prediction based on the dual-layer support vector machine (SVM) and position-specific scoring matrices (PSSMs). SVM is a new machine learning technology that has been successfully applied in solving problems in the field of bioinformatics. The SVM's performance is usually better than that of traditional machine learning approaches. The performance was further improved by combining PSSM profiles with the SVM analysis. The PSSMs were generated from PSI-BLAST profiles, which contain important evolution information. The final prediction results were generated from the second SVM layer output. On the CB513 data set, the three-state overall per-residue accuracy, Q3, reached 75.2%, while segment overlap (SOV) accuracy increased to 80.0%. On the CB396 data set, the Q3 of our method reached 74.0% and the SOV reached 78.1%. A web server utilizing the method has been constructed and is available at http://www.bioinfo.tsinghua.edu.cn/pmsvm.     

中国在翻译后修饰的生物信息学研究领域的进展与前瞻

翻译后修饰在调控蛋白质构象变化、活性以及功能方面具有重要作用,并参与了几乎所有细胞通路和过程。蛋白质翻译后修饰的鉴定是阐明细胞内分子机理的基础。相对于劳动密集的、耗费时间的实验工作,利用各种生物信息学方法开展翻译后修饰预测,能够提供准确、简便和快速的研究方案,并产生有价值的信息为进一步实验研究提供参考。文章主要综述了中国生物信息学者在翻译后修饰生物信息学领域所取得的研究进展,包括修饰底物与位点预测的计算方法学设计与完善、在线或本地化工具的设计与维护、修饰相关数据库及数据资源的构建及基于修饰蛋白质组学数据的生物信息学分析。通过比较国内外的同类研究,发现优势和不足,并对未来的研究作出前瞻。     

基于间隔二肽组分和递归特征消除法的DNA结合蛋白的鉴定

DNA结合蛋白(DNA-binding proteins,DBPs)的鉴定在原核和真核生物的基因和蛋白质功能注释研究中具有十分重要的意义.本研究首次运用间隔二肽组分(gapped-dipeptide composition,Gap DPC)结合递归特征消除法(recursive feature elimination,RFE)鉴定DBPs.首先获得待测蛋白质氨基酸序列的位置特异性得分矩阵(position specific scoring matrix,PSSM),在此基础上提取蛋白质的Gap DPC特征,通过RFE法选择最优特征,然后利用支持向量机(support vector machine,SVM)作为分类器,在蛋白质序列数据集PDB396和LB1068中进行夹克刀交叉验证(jackknife cross validation test).研究结果显示,基于PDB396和LB1068数据集,DBPs预测的准确率、Matthews相关系数、敏感性和特异性分别达到93.43%、0.86、89.04%和96.00%,以及86.33%、0.73、86.49%和86.18%,明显优于文献报道中的相关方法,为DBPs的鉴定提供了新的模型.     

Post-translational modifications of apolipoprotein A-I and Po proteins in the avian peripheral nerve

Apolipoprotein A-I (apo A-I), a soluble lipid transporter, and Po, the major glycoprotein of myelin, are actively synthesized during myelination. To explore the status of post-translational modifications of these proteins in the avian PNS during rapid myelination, endoneurial slices from one day old chick sciatic nerves were incubated with various radioactive precursors that could serve as indicators of such processes. The proteins were isolated from the incubation medium (secreted fraction), the 1% Triton-X-100-soluble intracellular-endoneurial (intracellular) fraction, and myelin-related and purified compact myelin fractions by immunoprecipitation with monospecific anti-apo A-I or anti-Po antisera. Our results demonstrated that secreted apo A-I is fatty acylated, but not phosphorylated or sulfated. Avian Po protein was phosphorylated by a phorbol ester sensitive protein kinase. Sulfation, as well as fatty acylation, of avian Po protein was observed in organ culture using highly sensitive methods of detection. These results indicate that fatty acylation of secreted apo A-I and phosphorylation, sulfation and fatty acylation of Po have been conserved during evolution, and that these post-translational modifications may play a common function in various species.     

Mapping of Post-translational Modifications of Transition Proteins,TP1 and TP2, and Identification of Protein Arginine Methyltransferase 4 and Lysine Methyltransferase 7 as Methyltransferase for TP2

In a unique global chromatin remodeling process during mammalian spermiogenesis, 90% of the nucleosomal histones are replaced by testis-specific transition proteins, TP1, TP2, and TP4. These proteins are further substituted by sperm-specific protamines, P1 and P2, to form a highly condensed sperm chromatin. In spermatozoa, a small proportion of chromatin, which ranges from 1 to 10% in mammals, retains the nucleosomal architecture and is implicated to play a role in transgenerational inheritance. However, there is still no mechanistic understanding of the interaction of chromatin machinery with histones and transition proteins, which facilitate this selective histone replacement from chromatin. Here, we report the identification of 16 and 19 novel post-translational modifications on rat endogenous transition proteins, TP1 and TP2, respectively, by mass spectrometry. By in vitro assays and mutational analysis, we demonstrate that protein arginine methyltransferase PRMT4 (CARM1) methylates TP2 at Arg⁷¹, Arg⁷⁵, and Arg⁹² residues, and lysine methyltransferase KMT7 (Set9) methylates TP2 at Lys⁸⁸ and Lys⁹¹ residues. Further studies with modification-specific antibodies that recognize TP2K88me1 and TP2R92me1 modifications showed that they appear in elongating to condensing spermatids and predominantly associated with the chromatin-bound TP2. This work establishes the repertoire of post-translational modifications that occur on TP1 and TP2, which may play a significant role in various chromatin-templated events during spermiogenesis and in the establishment of the sperm epigenome.     

Strategies for the effective identification of remotely related sequences in multiple PSSM search approach

Searches using position specific scoring matrices (PSSMs) have been commonly used in remote homology detection procedures such as PSI-BLAST and RPS-BLAST. A PSSM is generated typically using one of the sequences of a family as the reference sequence. In the case of PSI-BLAST searches the reference sequence is same as the query. Recently we have shown that searches against the database of multiple family-profiles, with each one of the members of the family used as a reference sequence, are more effective than searches against the classical database of single family-profiles. Despite relatively a better overall performance when compared with common sequence-profile matching procedures, searches against the multiple family-profiles database result in a few false positives and false negatives. Here we show that profile length and divergence of sequences used in the construction of a PSSM have major influence on the performance of multiple profile based search approach. We also identify that a simple parameter defined by the number of PSSMs corresponding to a family that is hit, for a query, divided by the total number of PSSMs in the family can distinguish effectively the true positives from the false positives in the multiple profiles search approach.     

PBOND： Web Server for the Prediction of Proline and Non-Proline cis / trans Isomerization

PBOND is a web server that predicts the conformation of the peptide bond between any two amino acids. PBOND classifies the peptide bonds into one out of four classes, namely cis imide (cis-Pro), cis amide (cis-nonPro), trans imide (trans-Pro) and trans amide (trans-nonPro). Moreover, for every prediction a reliability index is computed. The underlying structure of the server consists of three stages: (1) feature extraction, (2) feature selection and (3) peptide bond clas- sification. PBOND can handle both s...     

Functional annotation by sequence-weighted structure alignments: statistical analysis and case studies from the Protein 3000 structural genomics project in Japan

A method to functionally annotate structural genomics targets, based on a novel structural alignment scoring function, is proposed. In the proposed score, position-specific scoring matrices are used to weight structurally aligned residue pairs to highlight evolutionarily conserved motifs. The functional form of the score is first optimized for discriminating domains belonging to the same Pfam family from domains belonging to different families but the same CATH or SCOP superfamily. In the optimization stage, we consider four standard weighting functions as well as our own, the maximum substitution probability, and combinations of these functions. The optimized score achieves an area of 0.87 under the receiver-operating characteristic curve with respect to identifying Pfam families within a sequence-unique benchmark set of domain pairs. Confidence measures are then derived from the benchmark distribution of true-positive scores. The alignment method is next applied to the task of functionally annotating 230 query proteins released to the public as part of the Protein 3000 structural genomics project in Japan. Of these queries, 78 were found to align to templates with the same Pfam family as the query or had sequence identities > or = 30%. Another 49 queries were found to match more distantly related templates. Within this group, the template predicted by our method to be the closest functional relative was often not the most structurally similar. Several nontrivial cases are discussed in detail. Finally, 103 queries matched templates at the fold level, but not the family or superfamily level, and remain functionally uncharacterized.     

Post-translational excision of the carboxyl-terminal segment of CaM kinase phosphatase N and its cytosolic occurrence in the brain

Ca2+/Calmodulin-dependent protein kinase (CaM kinase) phosphatase, occurring in the cytoplasm of all tissues, dephosphorylates and thereby deactivates multifunctional CaM kinases, such as CaM kinases I, II and IV. In contrast, CaM kinase phosphatase N has been reported to occur almost exclusively in the brain and to be localized in the nucleus in the transfected COS-7 cells, as examined immunocytochemically with antibodies against the carboxyl-terminal segment of the enzyme, indicating its involvement in the deactivation of CaM kinase IV. Here, we show that the majority of the naturally occurring CaM kinase phosphatase N in the brain exists not in the intact form of the enzyme (83.4 kDa) but in a form (61.1 kDa) in which the carboxyl-terminal segment containing nuclear localization signals is deleted, and that it is present mostly in the cytoplasm but a little in the nucleus throughout the central nervous system, although occurring mostly in the nucleus in some large neurons. Strong immunostaining of the enzyme was also observed at postsynaptic density. These findings suggest that CaM kinase phosphatase N is involved in the regulation of not only CaM kinase IV but also CaM kinases II and I.     

Post-translational modification of barley 14-3-3A is isoform-specific and involves removal of the hypervariable C-terminus

The 14-3-3 protein family is a family of regulatory proteins involved in diverse cellular processes. In a previous study of regulation of individual 14-3-3 isoforms in the germinating barley embryo, we found that a post-translationally modified, 28 kDa form of 14-3-3A was present in specific cell fractions of the germinated embryo. In the present study, we identify the nature of the modification of 14-3-3A, and show that the 28 kDa doublet is the result of cleavage of the C-terminus. The 28 kDa forms of 14-3-3A lack ten or twelve amino acid residues at the non-conserved C-terminus of the protein, respectively. Barley 14-3-3B and 14-3-3C are not modified in a similar way. Like the 30 kDa form, in vitro produced 28 kDa 14-3-3A is still capable of binding AHA2 H⁺-ATPase in an overlay assay. Our results show a novel isoform-specific post-translational modification of 14-3-3 proteins that is regulated in a tissue-specific and developmental way.     

YASSPP: better kernels and coding schemes lead to improvements in protein secondary structure prediction

The accurate prediction of a protein's secondary structure plays an increasingly critical role in predicting its function and tertiary structure, as it is utilized by many of the current state-of-the-art methods for remote homology, fold recognition, and ab initio structure prediction. We developed a new secondary structure prediction algorithm called YASSPP, which uses a pair of cascaded models constructed from two sets of binary SVM-based models. YASSPP uses an input coding scheme that combines both position-specific and nonposition-specific information, utilizes a kernel function designed to capture the sequence conservation signals around the local window of each residue, and constructs a second-level model by incorporating both the three-state predictions produced by the first-level model and information about the original sequence. Experiments on three standard datasets (RS126, CB513, and EVA common subset 4) show that YASSPP is capable of producing the highest Q3 and SOV scores than that achieved by existing widely used schemes such as PSIPRED, SSPro 4.0, SAM-T99sec, as well as previously developed SVM-based schemes. On the EVA dataset it achieves a Q3 and SOV score of 79.34 and 78.65%, which are considerably higher than the best reported scores of 77.64 and 76.05%, respectively.     

PSLDoc: Protein subcellular localization prediction based on gapped-dipeptides and probabilistic latent semantic analysis

Prediction of protein subcellular localization (PSL) is important for genome annotation, protein function prediction, and drug discovery. Many computational approaches for PSL prediction based on protein sequences have been proposed in recent years for Gram-negative bacteria. We present PSLDoc, a method based on gapped-dipeptides and probabilistic latent semantic analysis (PLSA) to solve this problem. A protein is considered as a term string composed by gapped-dipeptides, which are defined as any two residues separated by one or more positions. The weighting scheme of gapped-dipeptides is calculated according to a position specific score matrix, which includes sequence evolutionary information. Then, PLSA is applied for feature reduction, and reduced vectors are input to five one-versus-rest support vector machine classifiers. The localization site with the highest probability is assigned as the final prediction. It has been reported that there is a strong correlation between sequence homology and subcellular localization (Nair and Rost, Protein Sci 2002;11:2836-2847; Yu et al., Proteins 2006;64:643-651). To properly evaluate the performance of PSLDoc, a target protein can be classified into low- or high-homology data sets. PSLDoc's overall accuracy of low- and high-homology data sets reaches 86.84% and 98.21%, respectively, and it compares favorably with that of CELLO II (Yu et al., Proteins 2006;64:643-651). In addition, we set a confidence threshold to achieve a high precision at specified levels of recall rates. When the confidence threshold is set at 0.7, PSLDoc achieves 97.89% in precision which is considerably better than that of PSORTb v.2.0 (Gardy et al., Bioinformatics 2005;21:617-623). Our approach demonstrates that the specific feature representation for proteins can be successfully applied to the prediction of protein subcellular localization and improves prediction accuracy. Besides, because of the generality of the representation, our method can be extended to eukaryotic proteomes in the future. The web server of PSLDoc is publicly available at http://bio-cluster.iis.sinica.edu.tw/~ bioapp/PSLDoc/.     

FTMove: A Web Server for Detection and Analysis of Cryptic and Allosteric Binding Sites by Mapping Multiple Protein Structures

Protein mapping distributes many copies of different molecular probes on the surface of a target protein in order to determine binding hot spots, regions that are highly preferable for ligand binding. While mapping of X-ray structures by the FTMap server is inherently static, this limitation can be overcome by the simultaneous analysis of multiple structures of the protein. FTMove is an automated web server that implements this approach. From the input of a target protein, by PDB code, the server identifies all structures of the protein available in the PDB, runs mapping on them, and combines the results to form binding hot spots and binding sites. The user may also upload their own protein structures, bypassing the PDB search for similar structures. Output of the server consists of the consensus binding sites and the individual mapping results for each structure - including the number of probes located in each binding site, for each structure. This level of detail allows the users to investigate how the strength of a binding site relates to the protein conformation, other binding sites, and the presence of ligands or mutations. In addition, the structures are clustered on the basis of their binding properties. The use of FTMove is demonstrated by application to 22 proteins with known allosteric binding sites; the orthosteric and allosteric binding sites were identified in all but one case, and the sites were typically ranked among the top five. The FTMove server is publicly available at https://ftmove.bu.edu.     

Chromophore: apoprotein interactions in phytochrome A*

Phytochromes are molecular light switches by virtue of their photochromic red/far-red reversibility. The His-324 residue next to the chromophore-linked Cys-323 plays a critical role in conferring photochromism to the tetrapyrrole chromophore in native phytochrome A. The chromophore appears to be enclosed between the amphiphilic α-helical chains in a hydrophobic pocket. The absorbance maxima of both the Pr and the Pfr forms of pea phytochrome A are blue-shifted by 10 and 20 nm, respectively, upon C-terminal truncation. We speculate that the quaternary structure of the phytochrome A molecule involves some interactions of the C-terminal half with the chromophore domain. The Pfr conformation of phytochrome includes an amphiphilic α-helix of the amino terminal chain, which occurs in 113 ms after picosecond photoisomerization of the Pr form. Compared to α-helical folding, unfolding of the α-helix occurs faster in about 310 μs upon phototransformation of the Pfr form of phytochrome A. The photochromic transformation of phytochrome A modulates protein kinase-catalysed phosphorylation sites in vivo and in vitro, but only a subtle local change in conformation is detectable in the phosphorylated phytochromes. This suggests that the post-translational modification serves as a surface label, rather than a transducer-activating trigger, for the recognition of a putative phytochrome receptor.     

Expanding the Paradigm: Intrinsically Disordered Proteins and Allosteric Regulation

Allosteric regulatory processes are implicated at all levels of biological function. Recent advances in our understanding of the diverse and functionally significant class of intrinsically disordered proteins have identified a multitude of ways in which disordered proteins function within the confines of the allosteric paradigm. Allostery within or mediated by intrinsically disordered proteins ensures robust and efficient signal integration through mechanisms that would be extremely unfavorable or even impossible for globular protein interaction partners. Here, we highlight recent examples that indicate the breadth of biological outcomes that can be achieved through allosteric regulation by intrinsically disordered proteins. Ongoing and future work in this rapidly evolving area of research will expand our appreciation of the central role of intrinsically disordered proteins in ensuring the fidelity and efficiency of cellular regulation.