Geodesic slice sampling on the sphere

Probability measures on the sphere form an important class of statistical models and are used, for example, in modeling directional data or shapes. Due to their widespread use, but also as an algorithmic building block, efficient sampling of distributions on the sphere is highly desirable. We propose a shrinkage based and an idealized geodesic slice sampling Markov chain, designed to generate approximate samples from distributions on the sphere. In particular, the shrinkage based algorithm works in any dimension, is straight-forward to implement and has no tuning parameters. We verify reversibility and show that under weak regularity conditions geodesic slice sampling is uniformly ergodic. Numerical experiments show that the proposed slice samplers achieve excellent mixing on challenging targets including the Bingham distribution and mixtures of von Mises-Fisher distributions. In these settings our approach outperforms standard samplers such as random-walk Metropolis Hastings and Hamiltonian Monte Carlo.Comment: 53 pages, 10 figures in the main text, 3 figures and 1 table in the appendi

Graphical analysis of NMR structural quality and interactive contact map of NOE assignments in ARIA

BACKGROUND: The Ambiguous Restraints for Iterative Assignment (ARIA) approach is widely used for NMR structure determination. It is based on simultaneously calculating structures and assigning NOE through an iterative protocol. The final solution consists of a set of conformers and a list of most probable assignments for the input NOE peak list. RESULTS: ARIA was extended with a series of graphical tools to facilitate a detailed analysis of the intermediate and final results of the ARIA protocol. These additional features provide (i) an interactive contact map, serving as a tool for the analysis of assignments, and (ii) graphical representations of structure quality scores and restraint statistics. The interactive contact map between residues can be clicked to obtain information about the restraints and their contributions. Profiles of quality scores are plotted along the protein sequence, and contact maps provide information of the agreement with the data on a residue pair level. CONCLUSIONS: The g raphical tools and outputs described here significantly extend the validation and analysis possibilities of NOE assignments given by ARIA as well as the analysis of the quality of the final structure ensemble. These tools are included in the latest version of ARIA, which is available at http://aria.pasteur.fr. The Web site also contains an installation guide, a user manual and example calculations

Bayesian Evidence and Model Selection

In this paper we review the concepts of Bayesian evidence and Bayes factors, also known as log odds ratios, and their application to model selection. The theory is presented along with a discussion of analytic, approximate and numerical techniques. Specific attention is paid to the Laplace approximation, variational Bayes, importance sampling, thermodynamic integration, and nested sampling and its recent variants. Analogies to statistical physics, from which many of these techniques originate, are discussed in order to provide readers with deeper insights that may lead to new techniques. The utility of Bayesian model testing in the domain sciences is demonstrated by presenting four specific practical examples considered within the context of signal processing in the areas of signal detection, sensor characterization, scientific model selection and molecular force characterization.Comment: Arxiv version consists of 58 pages and 9 figures. Features theory, numerical methods and four application

A unifying probabilistic framework for analyzing residual dipolar couplings

Residual dipolar couplings provide complementary information to the nuclear Overhauser effect measurements that are traditionally used in biomolecular structure determination by NMR. In a de novo structure determination, however, lack of knowledge about the degree and orientation of molecular alignment complicates the analysis of dipolar coupling data. We present a probabilistic framework for analyzing residual dipolar couplings and demonstrate that it is possible to estimate the atomic coordinates, the complete molecular alignment tensor, and the error of the couplings simultaneously. As a by-product, we also obtain estimates of the uncertainty in the coordinates and the alignment tensor. We show that our approach encompasses existing methods for determining the alignment tensor as special cases, including least squares estimation, histogram fitting, and elimination of an explicit alignment tensor in the restraint energy

Cooperative binding: a multiple personality

Cooperative binding has been described in many publications and has been related to or defined by several different properties of the binding behavior of the ligand to the target molecule. In addition to the commonly used Hill coefficient, other characteristics such as a sigmoidal shape of the overall titration curve in a linear plot, a change of ligand affinity of the other binding sites when a site of the target molecule becomes occupied, or complex roots of the binding polynomial have been used to define or to quantify cooperative binding. In this work, we analyze how the different properties are related in the most general model for binding curves based on the grand canonical partition function and present several examples which highlight differences between the cooperativity characterizing properties which are discussed. Our results mainly show that among the presented definitions there are not two which fully coincide. Moreover, this work poses the question whether it can make sense to distinguish between positive and negative cooperativity based on the macroscopic binding isotherm only. This article shall emphasize that scientists who investigate cooperative effects in biological systems could help avoiding misunderstandings by stating clearly which kind of cooperativity they discuss.Facultad de Ciencias ExactasCentro Regional de Estudios Genómico

Does the act of surgery provoke activation of "latent" metastases in early breast cancer?

This paper is written in support of the challenging article by Retsky and colleagues in this issue of Breast Cancer Research, and develops on the idea that the act of surgery can provoke the outgrowth of dormant micrometastases, which often leads to the failure of screening to deliver its promise. The therapeutic consequence of this idea involves the use of antiangiogenic drugs before surgery