Evaluating evolutionary constraint on the rapidly evolving gene matK using protein composition

The rapidly evolving chloroplast matK gene has nucleotide and amino acid substitution rates suggestive of progression toward a pseudogene state. However, molecular evidence has demonstrated that matK is expressed and functional. We explore in this paper the underlying factors behind the mode and tempo of matK evolution that allow this protein coding gene to accommodate such elevated rates of substitution and yet maintain functionality. Conservative amino acid replacement may reconcile the fast evolutionary rate in matK with conservation in protein function. Based on this premise, we have examined putative amino acid sequences for MATK from across green plants to determine constraint on this protein as indicated by variation in composition of amino acid side chain category. Amino acids in the MATK ORF were divided into six categories based on chemical properties of their side chains: nonpolar, uncharged (pH 7), basic, acidic, aromatic, and special (amino acids that specifically affect protein structure, i.e., proline, glycine, and cysteine). The amount of standard deviation (SD) in side chain composition was used as a measure of variation and constraint, where a low SD implied high evolutionary constraint and a high SD implied low constraint. Further, we used secondary structure prediction to evaluate if conservation observed in side chain composition was reflected in stable predicted structure. The results of this study demonstrate evolutionary constraint on MATK, identify three regions of functional importance, show highly conserved secondary structure, and support the putative function of MATK as a group II intron maturase.