Bayesian Protein Sequence and Structure Alignment (1404.1556v3)

Abstract: The structure of a protein is crucial in determining its functionality, and is much more conserved than sequence during evolution. A key task in structural biology is to compare protein structures in order to determine evolutionary relationships, estimate the function of newly-discovered structures, and predict unknown structures. We propose a Bayesian method for protein structure alignment, with the prior on alignments based on functions which penalise ``gaps'' in the aligned sequences. We show how a broad class of penalty functions fits into this framework, and how the resulting posterior distribution can be efficiently sampled. A commonly-used gap penalty function is shown to be a special case, and we propose a new penalty function which alleviates an undesirable feature of the commonly-used penalty. We illustrate our method on benchmark data sets, and find it competes well with popular tools from computational biology. Our method has the benefit of being able to potentially explore multiple competing alignments and quantify their merits probabilistically. The framework naturally allows for further information such as amino acid sequence to be included, and could be adapted to other situations such as flexible proteins or domain swaps.