Profile hidden Markov models for analyzing similarities and dissimilarities in the bacterial reaction center and photosystem II.

The bacterial photosynthetic reaction center is the evolutionary ancestor of the Photosystem II reaction center. These proteins share the same fold and perform the same biological function. Nevertheless, the details of their molecular reaction mechanism differ. It is of significant biological and biochemical interest to determine which functional characteristics are conserved at the level of the protein sequences. Since the level of sequence identity between the bacterial photosynthetic reaction center and Photosystem II is low, a progressive multiple-sequence alignment leads to errors in identifying the conserved residues. In such a situation, profile hidden Markov models (pHMM) can be used to obtain reliable multiple-sequence alignments. We therefore constructed the pHMM with the help of a sequence alignment based on a structural superposition of both proteins. To validate the multiple-sequence alignments obtained with the pHMM, the conservation of residues with known functional importance was examined. Having confirmed the correctness of the multiple-sequence alignments, we analyzed the conservation of residues involved in hydrogen bonding and redox potential tuning of the cofactors. Our analysis reveals similarities and dissimilarities between the bacterial photosynthetic reaction center and Photosystem II at the protein sequence level, hinting at different charge separation and charge transfer mechanisms. The conservation analysis that we perform in this paper can be considered as a model for analyzing the conservation in proteins with a low level of sequence identity.