Fred A. Hamprecht

DEVELOPMENT AND ASSESSMENT OF NEW EXCHANGE-CORRELATION FUNCTIONALS

We recently presented a new method for developing generalized gradient approximation (GGA) exchange-correlation energy functionals, using a least-squares procedure involving numerical exchange-correlation potentials and experimental energetics and nuclear gradients. In this paper we use the same method to develop a new GGA functional, denoted HCTH, based on an expansion recently suggested by Becke [J. Chem. Phys. 107, 8554 (1997)]. For our extensive training set, the new functional yields improved energetics compared to both the BLYP and B3LYP functionals [Phys. Rev. A 38, 3098 (1988); Phys. Rev. B 37, 785 (1988); J. Chem. Phys. 98, 5648 (1993); J. Phys. Chem. 98, 11623 (1994)]. The geometries of these systems, together with those of a set of transition metal compounds, are shown to be an improvement over the BLYP functional, while the reaction barriers for six hydrogen abstraction reactions are comparable to those of B3LYP. These improvements are achieved without introducing any fraction of exact orbital...

Ilastik: Interactive learning and segmentation toolkit

Segmentation is the process of partitioning digital images into meaningful regions. The analysis of biological high content images often requires segmentation as a first step. We propose ilastik as an easy-to-use tool which allows the user without expertise in image processing to perform segmentation and classification in a unified way. ilastik learns from labels provided by the user through a convenient mouse interface. Based on these labels, ilastik infers a problem specific segmentation. A random forest classifier is used in the learning step, in which each pixels neighborhood is characterized by a set of generic (nonlinear) features. ilastik supports up to three spatial plus one spectral dimension and makes use of all dimensions in the feature calculation. ilastik provides realtime feedback that enables the user to interactively refine the segmentation result and hence further fine-tune the classifier. An uncertainty measure guides the user to ambiguous regions in the images. Real time performance is achieved by multi-threading which fully exploits the capabilities of modern multi-core machines. Once a classifier has been trained on a set of representative images, it can be exported and used to automatically process a very large number of images (e.g. using the CellProfiler pipeline). ilastik is an open source project and released under the BSD license at www.ilastik.org.

A Comparative Study of Modern Inference Techniques for Discrete Energy Minimization Problems

Even years ago, Szeliski et al. published an influential study on energy minimization methods for Markov random fields (MRF). This study provided valuable insights in choosing the best optimization technique for certain classes of problems. While these insights remain generally useful today, the phenominal success of random field models means that the kinds of inference problems we solve have changed significantly. Specifically, the models today often include higher order interactions, flexible connectivity structures, large label-spaces of different cardinalities, or learned energy tables. To reflect these changes, we provide a modernized and enlarged study. We present an empirical comparison of 24 state-of-art techniques on a corpus of 2,300 energy minimization instances from 20 diverse computer vision applications. To ensure reproducibility, we evaluate all methods in the OpenGM2 framework and report extensive results regarding runtime and solution quality. Key insights from our study agree with the results of Szeliski et al. for the types of models they studied. However, on new and challenging types of models our findings disagree and suggest that polyhedral methods and integer programming solvers are competitive in terms of runtime and solution quality over a large range of model types.

Visualizing a Homogeneous Blend in Bulk Heterojunction Polymer Solar Cells by Analytical Electron Microscopy

To increase efficiency of bulk heterojunctions for photovoltaic devices, the functional morphology of active layers has to be understood, requiring visualization and discrimination of materials with very similar characteristics. Here we combine high-resolution spectroscopic imaging using an analytical transmission electron microscope with nonlinear multivariate statistical analysis for classification of multispectral image data. We obtain a visual representation showing homogeneous phases of donor and acceptor, connected by a third composite phase, depending in its extent on the way the heterojunction is fabricated. For the first time we can correlate variations in nanoscale morphology determined by material contrast with measured solar cell efficiency. In particular we visualize a homogeneously blended phase, previously discussed to diminish charge separation in solar cell devices.

Robust Prediction of the MASCOT Score for an Improved Quality Assessment in Mass Spectrometric Proteomics

Protein identification by tandem mass spectrometry is based on the reliable processing of the acquired data. Unfortunately, the generation of a large number of poor quality spectra is commonly observed in LC-MS/MS, and the processing of these mostly noninformative spectra with its associated costs should be avoided. We present a continuous quality score that can be computed very quickly and that can be considered an approximation of the MASCOT score in case of a correct identification. This score can be used to reject low quality spectra prior to database identification, or to draw attention to those spectra that exhibit a (supposedly) high information content, but could not be identified. The proposed quality score can be calibrated automatically on site without the need for a manually generated training set. When this score is turned into a classifier and when features are used that are independent of the instrument, the proposed approach performs equally to previously published classifiers and feature sets and also gives insights into the behavior of the MASCOT score.

Automated Detection and Segmentation of Synaptic Contacts in Nearly Isotropic Serial Electron Microscopy Images

We describe a protocol for fully automated detection and segmentation of asymmetric, presumed excitatory, synapses in serial electron microscopy images of the adult mammalian cerebral cortex, taken with the focused ion beam, scanning electron microscope (FIB/SEM). The procedure is based on interactive machine learning and only requires a few labeled synapses for training. The statistical learning is performed on geometrical features of 3D neighborhoods of each voxel and can fully exploit the high z-resolution of the data. On a quantitative validation dataset of 111 synapses in 409 images of 1948×1342 pixels with manual annotations by three independent experts the error rate of the algorithm was found to be comparable to that of the experts (0.92 recall at 0.89 precision). Our software offers a convenient interface for labeling the training data and the possibility to visualize and proofread the results in 3D. The source code, the test dataset and the ground truth annotation are freely available on the website http://www.ilastik.org/synapse-detection.

A Comparative Study of Modern Inference Techniques for Structured Discrete Energy Minimization Problems

Szeliski et al. published an influential study in 2006 on energy minimization methods for Markov random fields. This study provided valuable insights in choosing the best optimization technique for certain classes of problems. While these insights remain generally useful today, the phenomenal success of random field models means that the kinds of inference problems that have to be solved changed significantly. Specifically, the models today often include higher order interactions, flexible connectivity structures, large label-spaces of different cardinalities, or learned energy tables. To reflect these changes, we provide a modernized and enlarged study. We present an empirical comparison of more than 27 state-of-the-art optimization techniques on a corpus of 2453 energy minimization instances from diverse applications in computer vision. To ensure reproducibility, we evaluate all methods in the OpenGM 2 framework and report extensive results regarding runtime and solution quality. Key insights from our study agree with the results of Szeliski et al. for the types of models they studied. However, on new and challenging types of models our findings disagree and suggest that polyhedral methods and integer programming solvers are competitive in terms of runtime and solution quality over a large range of model types.

Segmentation of SBFSEM Volume Data of Neural Tissue by Hierarchical Classification

Three-dimensional electron-microscopic image stacks with almost isotropic resolution allow, for the first time, to determine the complete connection matrix of parts of the brain. In spite of major advances in staining, correct segmentation of these stacks remains challenging, because very few local mistakes can lead to severe global errors. We propose a hierarchical segmentation procedure based on statistical learning and topology-preserving grouping. Edge probability maps are computed by a random forest classifier (trained on hand-labeled data) and partitioned into supervoxels by the watershed transform. Over-segmentation is then resolved by another random forest. Careful validation shows that the results of our algorithm are close to human labelings.

Concise Representation of Mass Spectrometry Images by Probabilistic Latent Semantic Analysis

Imaging mass spectrometry (IMS) is a promising technology which allows for detailed analysis of spatial distributions of (bio)molecules in organic samples. In many current applications, IMS relies heavily on (semi)automated exploratory data analysis procedures to decompose the data into characteristic component spectra and corresponding abundance maps, visualizing spectral and spatial structure. The most commonly used techniques are principal component analysis (PCA) and independent component analysis (ICA). Both methods operate in an unsupervised manner. However, their decomposition estimates usually feature negative counts and are not amenable to direct physical interpretation. We propose probabilistic latent semantic analysis (pLSA) for non-negative decomposition and the elucidation of interpretable component spectra and abundance maps. We compare this algorithm to PCA, ICA, and non-negative PARAFAC (parallel factors analysis) and show on simulated and real-world data that pLSA and non-negative PARAFAC are superior to PCA or ICA in terms of complementarity of the resulting components and reconstruction accuracy. We further combine pLSA decomposition with a statistical complexity estimation scheme based on the Akaike information criterion (AIC) to automatically estimate the number of components present in a tissue sample data set and show that this results in sensible complexity estimates.

Different phosphorylation states of the anaphase promoting complex in response to antimitotic drugs: A quantitative proteomic analysis

The anaphase promoting complex (APC) controls the degradation of proteins during exit from mitosis and entry into S-phase. The activity of the APC is regulated by phosphorylation during mitosis. Because the phosphorylation pattern provides insights into the complexity of regulation of the APC, we studied in detail the phosphorylation patterns at a single mitotic state of arrest generated by various antimitotic drugs. We examined the phosphorylation patterns of the APC in HeLa S3 cells after they were arrested in prometaphase with taxol, nocodazole, vincristine, or monastrol. There were 71 phosphorylation sites on nine of the APC subunits. Despite the common state of arrest, the various antimitotic drug treatments resulted in differences in the phosphorylation patterns and phosphorylation stoichiometries. The relative phosphorylation stoichiometries were determined by using a method adapted from the isotope-free quantitation of the extent of modification (iQEM). We could show that during drug arrest the phosphorylation state of the APC changes, indicating that the mitotic arrest is not a static condition. We discuss these findings in terms of the variable efficacy of antimitotic drugs in cancer chemotherapy.