fuNTRp: Identifying protein positions for variation driven functional tuning

Evaluating the impact of non-synonymous genetic variants is essential for uncovering disease associations and mechanisms of evolution. Understanding corresponding sequence changes is also fundamental for synthetic protein design and stability assessments. However, the performance gain of variant effect predictors observed in recent years is not in line with the increased complexity of new methods. One likely reason for this might be that most approaches use similar sets of gene/protein features for modeling variant effect, often emphasizing sequence conservation. While high levels of conservation highlight residues essential for protein activity, much of the in vivo observable variation is arguably weaker in its impact and, thus, requires evaluation at a higher level of resolution. Here we describe function Neutral/Toggle/Rheostat predictor (funtrp), a novel computational method that categorizes protein positions based on the position-specific expected range of mutational impacts: Neutral (weak/no effects), Rheostat (function-tuning positions), or Toggle (on/off switches). We show that position types do not correlate strongly with familiar protein features such as conservation or protein disorder. We also find that position type distribution varies across different protein functions. Finally, we demonstrate that position types reflect experimentally determined functional effects and can thus improve performance of existing variant effect predictors and suggest a way forward for the development of new ones.