Collective estimation of multiple bivariate density functions with application to angular-sampling-based protein loop modeling

Mehdi Maadooliat, Lan Zhou, Seyed Morteza Najibi, Xin Gao, Jianhua Z. Huang

Research output: Contribution to journalArticlepeer-review

9 Scopus citations

Abstract

This paper develops a method for simultaneous estimation of density functions for a collection of populations of protein backbone angle pairs using a data-driven, shared basis that is constructed by bivariate spline functions defined on a triangulation of the bivariate domain. The circular nature of angular data is taken into account by imposing appropriate smoothness constraints across boundaries of the triangles. Maximum penalized likelihood is used to fit the model and an alternating blockwise Newton-type algorithm is developed for computation. A simulation study shows that the collective estimation approach is statistically more efficient than estimating the densities individually. The proposed method was used to estimate neighbor-dependent distributions of protein backbone dihedral angles (i.e., Ramachandran distributions). The estimated distributions were applied to protein loop modeling, one of the most challenging open problems in protein structure prediction, by feeding them into an angular-sampling-based loop structure prediction framework. Our estimated distributions compared favorably to the Ramachandran distributions estimated by fitting a hierarchical Dirichlet process model; and in particular, our distributions showed significant improvements on the hard cases where existing methods do not work well.
Original languageEnglish (US)
Pages (from-to)43-56
Number of pages14
JournalJournal of the American Statistical Association
Volume111
Issue number513
DOIs
StatePublished - May 5 2016

Fingerprint Dive into the research topics of 'Collective estimation of multiple bivariate density functions with application to angular-sampling-based protein loop modeling'. Together they form a unique fingerprint.

Cite this