Maximum likelihood estimation of subsequence conservation

A statistical method is presented for comparing protein sequences by partitioning the polymers and estimating each subsegment's degree of conservation. Conservation is measured as a function of the number of transitions occurring in the underlying time homogeneous Markov process assumed to govern amino acid mutations. The Markovian assumption also permits estimation of the ancestral sequence. Partitioning and estimation are carried out via maximum likelihood. The method is contrasted with the commonly utilized percent homology measure. A moving likelihood ratio plot to aid in identifying regions of high conservation is suggested as an analogue to moving hydrophobicity plots. An application is presented which identifies highly conserved regions in thymidylate synthase from L. casei and E. coli.