Detailed computational study of p53 and p16: using evolutionary sequence analysis and disease-associated mutations to predict the functional consequences of allelic variants

Deciding whether a missense allelic variant affects protein function is important in many contexts. We previously demonstrated that a detailed analysis of p53 intragenic conservation correlates with somatic mutation hotspots. Here we refine these evolutionary studies and expand them to the p16/Ink4a gene. We calculated that in order for ‘absolute conservation’ of a codon across multiple species to achieve P<0.05, the evolutionary substitution database must contain at least 3(M) variants, where M equals the number of codons in the gene. Codons in p53 were divided into high (73% of codons), intermediate (29% of codons), and low (0 codons) likelihood of being mutation hotspots. From a database of 263 somatic missense p16 mutations, we identified only four codons that are mutational hotspots at P<0.05 (8 mutations). However, data on function, structure, and disease association support the conclusion that 11 other codons with ≥5 somatic mutations also likely indicate functionally critical residues, even though P0.05. We calculated p16 evolution using amino acid substitution matrices and nucleotide substitution distances. We looked for evolutionary parameters at each codon that would predict whether missense mutations were disease associated or disrupted function. The current p16 evolutionary substitution database is too small to determine whether observations of ‘absolute conservation’ are statistically significant. Increasing the number of sequences from three to seven significantly improved the predictive value of evolutionary computations. The sensitivity and specificity for conservation scores in predicting disease association of p16 codons is 70–80%. Despite the small p16 sequence database, our calculations of high conservation correctly predicted loss of cell cycle arrest function in 75% of tested codons, and low conservation correctly predicted wild-type function in 80–90% of codons. These data validate our hypothesis that detailed evolutionary analyses help predict the consequences of missense amino-acid variants.