André Leier

POSSUM: a bioinformatics toolkit for generating numerical sequence feature descriptors based on PSSM profiles

Summary: Evolutionary information in the form of a Position‐Specific Scoring Matrix (PSSM) is a widely used and highly informative representation of protein sequences. Accordingly, PSSM‐based feature descriptors have been successfully applied to improve the performance of various predictors of protein attributes. Even though a number of algorithms have been proposed in previous studies, there is currently no universal web server or toolkit available for generating this wide variety of descriptors. Here, we present POSSUM (Position‐Specific Scoring matrix‐based feature generator for machine learning), a versatile toolkit with an online web server that can generate 21 types of PSSM‐based feature descriptors, thereby addressing a crucial need for bioinformaticians and computational biologists. We envisage that this comprehensive toolkit will be widely used as a powerful tool to facilitate feature extraction, selection, and benchmarking of machine learning‐based models, thereby contributing to a more effective analysis and modeling pipeline for bioinformatics research. Availability and implementation: http://possum.erc.monash.edu/. Contact: trevor.lithgow@monash.edu or jiangning.song@monash.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Integrative physical oncology

Cancer is arguably the ultimate complex biological system. Solid tumors are microstructured soft matter that evolves as a consequence of spatio‐temporal events at the intracellular (e.g., signaling pathways, macromolecular trafficking), intercellular (e.g., cell–cell adhesion/communication), and tissue (e.g., cell–extracellular matrix interactions, mechanical forces) scales. To gain insight, tumor and developmental biologists have gathered a wealth of molecular, cellular, and genetic data, including immunohistochemical measurements of cell type‐specific division and death rates, lineage tracing, and gain‐of‐function/loss‐of‐function mutational analyses. These data are empirically extrapolated to a diagnosis/prognosis of tissue‐scale behavior, e.g., for clinical decision. Integrative physical oncology (IPO) is the science that develops physically consistent mathematical approaches to address the significant challenge of bridging the nano (nm)–micro (µm) to macro (mm, cm) scales with respect to tumor development and progression. In the current literature, such approaches are referred to as multiscale modeling. In the present article, we attempt to assess recent modeling approaches on each separate scale and critically evaluate the current ‘hybrid‐multiscale’ models used to investigate tumor growth in the context of brain and breast cancers. Finally, we provide our perspective on the further development and the impact of IPO. WIREs Syst Biol Med 2012, 4:1–14. doi: 10.1002/wsbm.158

Anomalous diffusion and multifractional Brownian motion: simulating molecular crowding and physical obstacles in systems biology

There have been many recent studies from both experimental and simulation perspectives in order to understand the effects of spatial crowding in molecular biology. These effects manifest themselves in protein organisation on the plasma membrane, on chemical signalling within the cell and in gene regulation. Simulations are usually done with lattice- or meshless-based random walks but insights can also be gained through the computation of the underlying probability density functions of these stochastic processes. Until recently much of the focus had been on continuous time random walks, but some very recent work has suggested that fractional Brownian motion may be a good descriptor of spatial crowding effects in some cases. The study compares both fractional Brownian motion and continuous time random walks and highlights how well they can represent different types of spatial crowding and physical obstacles. Simulated spatial data, mimicking experimental data, was first generated by using the package Smoldyn. We then attempted to characterise this data through continuous time anomalously diffusing random walks and multifractional Brownian motion (MFBM) by obtaining MFBM paths that match the statistical properties of our sample data. Although diffusion around immovable obstacles can be reasonably characterised by a single Hurst exponent, we find that diffusion in a crowded environment seems to exhibit multifractional properties in the form of a different short- and long-time behaviour.

PROSPERous: high-throughput prediction of substrate cleavage sites for 90 proteases with improved accuracy

SummarynProteases are enzymes that specifically cleave the peptide backbone of their target proteins. As an important type of irreversible post-translational modification, protein cleavage underlies many key physiological processes. When dysregulated, proteases actions are associated with numerous diseases. Many proteases are highly specific, cleaving only those target substrates that present certain particular amino acid sequence patterns. Therefore, tools that successfully identify potential target substrates for proteases may also identify previously unknown, physiologically relevant cleavage sites, thus providing insights into biological processes and guiding hypothesis-driven experiments aimed at verifying protease-substrate interaction. In this work, we present PROSPERous, a tool for rapid in silico prediction of protease-specific cleavage sites in substrate sequences. Our tool is based on logistic regression models and uses different scoring functions and their pairwise combinations to subsequently predict potential cleavage sites. PROSPERous represents a state-of-the-art tool that enables fast, accurate and high-throughput prediction of substrate cleavage sites for 90 proteases.nnnAvailability and implementationnhttp://prosperous.erc.monash.edu/.nnnContactnjiangning.song@monash.edu or geoff.webb@monash.edu or r.pike@latrobe.edu.au.nnnSupplementary informationnSupplementary data are available at Bioinformatics online.

Probability distributed time delays: integrating spatial effects into temporal models

BackgroundIn order to provide insights into the complex biochemical processes inside a cell, modelling approaches must find a balance between achieving an adequate representation of the physical phenomena and keeping the associated computational cost within reasonable limits. This issue is particularly stressed when spatial inhomogeneities have a significant effect on systems behaviour. In such cases, a spatially-resolved stochastic method can better portray the biological reality, but the corresponding computer simulations can in turn be prohibitively expensive.ResultsWe present a method that incorporates spatial information by means of tailored, probability distributed time-delays. These distributions can be directly obtained by single in silico or a suitable set of in vitro experiments and are subsequently fed into a delay stochastic simulation algorithm (DSSA), achieving a good compromise between computational costs and a much more accurate representation of spatial processes such as molecular diffusion and translocation between cell compartments. Additionally, we present a novel alternative approach based on delay differential equations (DDE) that can be used in scenarios of high molecular concentrations and low noise propagation.ConclusionsOur proposed methodologies accurately capture and incorporate certain spatial processes into temporal stochastic and deterministic simulations, increasing their accuracy at low computational costs. This is of particular importance given that time spans of cellular processes are generally larger (possibly by several orders of magnitude) than those achievable by current spatially-resolved stochastic simulators. Hence, our methodology allows users to explore cellular scenarios under the effects of diffusion and stochasticity in time spans that were, until now, simply unfeasible. Our methodologies are supported by theoretical considerations on the different modelling regimes, i.e. spatial vs. delay-temporal, as indicated by the corresponding Master Equations and presented elsewhere.

PhosphoPredict: A bioinformatics tool for prediction of human kinase-specific phosphorylation substrates and sites by integrating heterogeneous feature selection

Protein phosphorylation is a major form of post-translational modification (PTM) that regulates diverse cellular processes. In silico methods for phosphorylation site prediction can provide a useful and complementary strategy for complete phosphoproteome annotation. Here, we present a novel bioinformatics tool, PhosphoPredict, that combines protein sequence and functional features to predict kinase-specific substrates and their associated phosphorylation sites for 12 human kinases and kinase families, including ATM, CDKs, GSK-3, MAPKs, PKA, PKB, PKC, and SRC. To elucidate critical determinants, we identified feature subsets that were most informative and relevant for predicting substrate specificity for each individual kinase family. Extensive benchmarking experiments based on both five-fold cross-validation and independent tests indicated that the performance of PhosphoPredict is competitive with that of several other popular prediction tools, including KinasePhos, PPSP, GPS, and Musite. We found that combining protein functional and sequence features significantly improves phosphorylation site prediction performance across all kinases. Application of PhosphoPredict to the entire human proteome identified 150 to 800 potential phosphorylation substrates for each of the 12 kinases or kinase families. PhosphoPredict significantly extends the bioinformatics portfolio for kinase function analysis and will facilitate high-throughput identification of kinase-specific phosphorylation sites, thereby contributing to both basic and translational research programs.

Exact model reduction with delays: closed-form distributions and extensions to fully bi-directional monomolecular reactions

In order to systematically understand the qualitative and quantitative behaviour of chemical reaction networks, scientists must derive and analyse associated mathematical models. However, biochemical systems are often very large, with reactions occurring at multiple time scales, as evidenced by signalling pathways and gene expression kinetics. Owing to the associated computational costs, it is then many times impractical, if not impossible, to solve or simulate these systems with an appropriate level of detail. By consequence, there is a growing interest in developing techniques for the simplification or reduction of complex biochemical systems. Here, we extend our recently presented methodology on exact reduction of linear chains of reactions with delay distributions in two ways. First, we report that it is now possible to deal with fully bi-directional monomolecular systems, including degradations, synthesis and generalized bypass reactions. Second, we provide all derivations of associated delays in analytical, closed form. Both advances have a major impact on further reducing computational costs, while still retaining full accuracy. Thus, we expect our new methodology to respond to current simulation needs in pharmaceutical, chemical and biological research.

iFeature: a Python package and web server for features extraction and selection from protein and peptide sequences

Summary: Structural and physiochemical descriptors extracted from sequence data have been widely used to represent sequences and predict structural, functional, expression and interaction profiles of proteins and peptides as well as DNAs/RNAs. Here, we present iFeature, a versatile Python‐based toolkit for generating various numerical feature representation schemes for both protein and peptide sequences. iFeature is capable of calculating and extracting a comprehensive spectrum of 18 major sequence encoding schemes that encompass 53 different types of feature descriptors. It also allows users to extract specific amino acid properties from the AAindex database. Furthermore, iFeature integrates 12 different types of commonly used feature clustering, selection and dimensionality reduction algorithms, greatly facilitating training, analysis and benchmarking of machine‐learning models. The functionality of iFeature is made freely available via an online web server and a stand‐alone toolkit. Availability and implementation: http://iFeature.erc.monash.edu/; https://github.com/Superzchen/iFeature/. Supplementary information: Supplementary data are available at Bioinformatics online.

Bastion6: A bioinformatics approach for accurate prediction of type VI secreted effectors

Motivation: Many Gram‐negative bacteria use type VI secretion systems (T6SS) to export effector proteins into adjacent target cells. These secreted effectors (T6SEs) play vital roles in the competitive survival in bacterial populations, as well as pathogenesis of bacteria. Although various computational analyses have been previously applied to identify effectors secreted by certain bacterial species, there is no universal method available to accurately predict T6SS effector proteins from the growing tide of bacterial genome sequence data. Results: We extracted a wide range of features from T6SE protein sequences and comprehensively analyzed the prediction performance of these features through unsupervised and supervised learning. By integrating these features, we subsequently developed a two‐layer SVM‐based ensemble model with fine‐grain optimized parameters, to identify potential T6SEs. We further validated the predictive model using an independent dataset, which showed that the proposed model achieved an impressive performance in terms of ACC (0.943), F‐value (0.946), MCC (0.892) and AUC (0.976). To demonstrate applicability, we employed this method to correctly identify two very recently validated T6SE proteins, which represent challenging prediction targets because they significantly differed from previously known T6SEs in terms of their sequence similarity and cellular function. Furthermore, a genome‐wide prediction across 12 bacterial species, involving in total 54 212 protein sequences, was carried out to distinguish 94 putative T6SE candidates. We envisage both this information and our publicly accessible web server will facilitate future discoveries of novel T6SEs. Availability and implementation: http://bastion6.erc.monash.edu/ Supplementary information: Supplementary data are available at Bioinformatics online.

Knowledge-Transfer learning for prediction of matrix metalloprotease substrate-cleavage sites

Matrix Metalloproteases (MMPs) are an important family of proteases that play crucial roles in key cellular and disease processes. Therefore, MMPs constitute important targets for drug design, development and delivery. Advanced proteomic technologies have identified type-specific target substrates; however, the complete repertoire of MMP substrates remains uncharacterized. Indeed, computational prediction of substrate-cleavage sites associated with MMPs is a challenging problem. This holds especially true when considering MMPs with few experimentally verified cleavage sites, such as for MMP-2, -3, -7, and -8. To fill this gap, we propose a new knowledge-transfer computational framework which effectively utilizes the hidden shared knowledge from some MMP types to enhance predictions of other, distinct target substrate-cleavage sites. Our computational framework uses support vector machines combined with transfer machine learning and feature selection. To demonstrate the value of the model, we extracted a variety of substrate sequence-derived features and compared the performance of our method using both 5-fold cross-validation and independent tests. The results show that our transfer-learning-based method provides a robust performance, which is at least comparable to traditional feature-selection methods for prediction of MMP-2, -3, -7, -8, -9 and -12 substrate-cleavage sites on independent tests. The results also demonstrate that our proposed computational framework provides a useful alternative for the characterization of sequence-level determinants of MMP-substrate specificity.