Olav Zimmermann

Comprehensive large-scale assessment of intrinsic protein disorder

MOTIVATION Intrinsically disordered regions are key for the function of numerous proteins. Due to the difficulties in experimental disorder characterization, many computational predictors have been developed with various disorder flavors. Their performance is generally measured on small sets mainly from experimentally solved structures, e.g. Protein Data Bank (PDB) chains. MobiDB has only recently started to collect disorder annotations from multiple experimental structures. RESULTS MobiDB annotates disorder for UniProt sequences, allowing us to conduct the first large-scale assessment of fast disorder predictors on 25 833 different sequences with X-ray crystallographic structures. In addition to a comprehensive ranking of predictors, this analysis produced the following interesting observations. (i) The predictors cluster according to their disorder definition, with a consensus giving more confidence. (ii) Previous assessments appear over-reliant on data annotated at the PDB chain level and performance is lower on entire UniProt sequences. (iii) Long disordered regions are harder to predict. (iv) Depending on the structural and functional types of the proteins, differences in prediction performance of up to 10% are observed. AVAILABILITY The datasets are available from Web site at URL: http://mobidb.bio.unipd.it/lsd. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Support vector machines for prediction of dihedral angle regions

MOTIVATION Most secondary structure prediction programs target only alpha helix and beta sheet structures and summarize all other structures in the random coil pseudo class. However, such an assignment often ignores existing local ordering in so-called random coil regions. Signatures for such ordering are distinct dihedral angle pattern. For this reason, we propose as an alternative approach to predict directly dihedral regions for each residue as this leads to a higher amount of structural information. RESULTS We propose a multi-step support vector machine (SVM) procedure, dihedral prediction (DHPRED), to predict the dihedral angle state of residues from sequence. Trained on 20,000 residues our approach leads to dihedral region predictions, that in regions without alpha helices or beta sheets is higher than those from secondary structure prediction programs. AVAILABILITY DHPRED has been implemented as a web service, which academic researchers can access from our webpage http://www.fz-juelich.de/nic/cbb

Simulation of Top7-CFr: a transient helix extension guides folding.

Protein structures often feature β-sheets in which adjacent β-strands have large sequence separation. How the folding process orchestrates the formation and correct arrangement of these strands is not comprehensively understood. Particularly challenging are proteins in which β-strands at the N and C termini are neighbors in a β-sheet. The N-terminal β-strand is synthesized early on, but it can not bind to the C terminus before the chain is fully synthesized. During this time, there is a danger that the β-strand at the N terminus interacts with nearby molecules, leading to potentially harmful aggregates of incompletely folded proteins. Simulations of the C-terminal fragment of Top7 show that this risk of misfolding and aggregation can be avoided by a “caching” mechanism that relies on the “chameleon” behavior of certain segments.

LOCUSTRA: accurate prediction of local protein structure using a two-layer support vector machine approach.

Constraint generation for 3d structure prediction and structure-based database searches benefit from fine-grained prediction of local structure. In this work, we present LOCUSTRA, a novel scheme for the multiclass prediction of local structure that uses two layers of support vector machines (SVM). Using a 16-letter structural alphabet from de Brevern et al. (Proteins: Struct., Funct., Bioinf. 2000, 41, 271-287), we assess its prediction ability for an independent test set of 222 proteins and compare our method to three-class secondary structure prediction and direct prediction of dihedral angles. The prediction accuracy is Q16=61.0% for the 16 classes of the structural alphabet and Q3=79.2% for a simple mapping to the three secondary classes helix, sheet, and coil. We achieve a mean phi(psi) error of 24.74 degrees (38.35 degrees) and a median RMSDA (root-mean-square deviation of the (dihedral) angles) per protein chain of 52.1 degrees. These results compare favorably with related approaches. The LOCUSTRA web server is freely available to researchers at http://www.fz-juelich.de/nic/cbb/service/service.php.

A new optimization phase for scientific workflow management systems

Scientific workflows have emerged as an important tool for combining computational power with data analysis for all scientific domains in e-science. They help scientists to design and execute complex in silico experiments. However, with increasing complexity it becomes more and more infeasible to optimize scientific workflows by trial and error. To address this issue, this paper describes the design of a new optimization phase integrated in the established scientific workflow life cycle. We have also developed a flexible optimization application programming interface (API) and have integrated it into a scientific workflow management system. A sample plugin for parameter optimization based on genetic algorithms illustrates, how the API enables rapid implementation of concrete workflow optimization methods. Finally, a use case taken from the area of structural bioinformatics validates how the optimization approach facilitates setup, execution and monitoring of workflow parameter optimization in high performance computing e-science environments.

The effect of composition on Diffusion of macromolecules in a crowded environment

We study diffusion of macromolecules in a crowded cytoplasm-like environment, focusing on its dependence on composition and its crossover to the anomalous subdiffusion. The crossover and the diffusion itself depend on both the volume fraction and the relative concentration of macromolecules. In accordance with previous theoretical and experimental studies, diffusion slows down when the volume fraction increases. Contrary to expectations, however, the diffusion is also strongly dependent on the molecular composition. The crossover time decreases and diffusion slows down when the smaller macromolecules start to dominate. Interestingly, diffusion is faster in a cytoplasm-like (more polydisperse) system than it is in a two-component system, at comparable packing fractions, or even when the cytoplasm packing fraction is larger.

Folding of Top7 in unbiased all-atom Monte Carlo simulations

For computational studies of protein folding, proteins with both helical and β‐sheet secondary structure elements are very challenging, as they expose subtle biases of the physical models. Here, we present reproducible folding of a 92 residue α/β protein (residues 3–94 of Top7, PDB ID: 1QYS) in computer simulations starting from random initial conformations using a transferable physical model which has been previously shown to describe the folding and thermodynamic properties of about 20 other smaller proteins of different folds. Top7 is a de novo designed protein with two α‐helices and a five stranded β‐sheet. Experimentally, it is known to be unusually stable for its size, and its folding transition distinctly deviates from the two‐state behavior commonly seen in natural single domain proteins. In our all‐atom implicit solvent parallel tempering Monte Carlo simulations, Top7 shows a rapid transition to a group of states with high native‐like secondary structure, and a much slower subsequent transition to the native state with a root mean square deviation of about 3.5 Å from the experimentally determined structure. Consistent with experiments, we find Top7 to be thermally extremely stable, although the simulations also find a large number of very stable non‐native states with high native‐like secondary structure. Proteins 2013; 81:1446–1456.

A UNICORE Plugin for HPC-Enabled Scientific Workflows in Taverna 2.2

As scientific workflows are becoming more complex and apply compute-intensive methods to increasingly large data volumes, access to HPC resources is becoming mandatory. We describe the development of a novel plug in for the Tavern a workflow system, which provides transparent and secure access to HPC/grid resources via the UNICORE grid middleware, while maintaining the ease of use that has been the main reason for the success of scientific workflow systems. A use case from the bioinformatics domain demonstrates the potential of the UNICORE plug in for Tavern a by creating a scientific workflow that executes the central parts in parallel on a cluster resource.

On specifying and sharing scientific workflow optimization results using research objects

Reusing and repurposing scientific workflows for novel scientific experiments is nowadays facilitated by workflow repositories. Such repositories allow scientists to find existing workflows and re-execute them. However, workflow input parameters often need to be adjusted to the research problem at hand. Adapting these parameters may become a daunting task due to the infinite combinations of their values in a wide range of applications. Thus, a scientist may preferably use an automated optimization mechanism to adjust the workflow set-up and improve the result. Currently, automated optimizations must be started from scratch as optimization meta-data are not stored together with workflow provenance data. This important meta-data is lost and can neither be reused nor assessed by other researchers. In this paper we present a novel approach to capture optimization meta-data by extending the Research Object model and reusing the W3C standards. We validate our proposal through a real-world use case taken from the biodivertsity domain, and discuss the exploitation of our solution in the context of existing e-Science infrastructures.

Discrete-continuous reaction-diffusion model with mobile point-like sources and sinks

Abstract.In many applications in soft and biological physics, there are multiple time and length scales involved but often with a distinct separation between them. For instance, in enzyme kinetics, enzymes are relatively large, move slowly and their copy numbers are typically small, while the metabolites (being transformed by these enzymes) are often present in abundance, are small in size and diffuse fast. It seems thus natural to apply different techniques to different time and length levels and couple them. Here we explore this possibility by constructing a stochastic-deterministic discrete-continuous reaction-diffusion model with mobile sources and sinks. Such an approach allows in particular to separate different sources of stochasticity. We demonstrate its application by modelling enzyme-catalysed reactions with freely diffusing enzymes and a heterogeneous source of metabolites. Our calculations suggest that using a higher amount of less active enzymes, as compared to fewer more active enzymes, reduces the metabolite pool size and correspondingly the lag time, giving rise to a faster response to external stimuli. The methodology presented can be extended to more complex systems and offers exciting possibilities for studying problems where spatial heterogeneities, stochasticity or discreteness play a role.Graphical abstract