Generative AI impact on protein stability prediction in breast cancer genes

Abstract

AbstractThe functional classification of a missense variant in cancer predisposition genes is often challenging due to how rare the variant is observed in the population. When available, clinicians utilize a combination of family history,in vitrofunctional assays andin silicomethods to infer protein function.In silicomethods, such as missense predictors (predict changes in protein function) and protein stability predictors (predict changes in free energy) have been used to help classify a missense variant in accordance with the American College of Medical Genetics and Genomics (ACMG) guideline. To measure protein stability, manyin silicoalgorithms predict stability based on the change of free energy and most accurate protein stability predictors require a wild-type protein template. In this study, we examine the use of generative AI to predict high-resolution protein structures as templates analyzed with protein stability methods to evaluate loss of function (LOF) activity in cancer predisposition genesBRCA1, BRCA2, PALB2andRAD51Cupon the presence of missense variant. Utilizing multiplexed assay of variant effect measurements and variant classifications from ClinVar, we find that prediction of Gibbs free energy (ΔΔG) from AlphaFold2 (AF2) structures analyzed with FoldX predicts LOF better than experimental-derived wild type structures in the BRCT domain ofBRCA1and the DNA binding domain (DBD) ofBRCA2, but not inPALB2andRAD51C. We also find that AF2 structures in the BRCT domain of BRCA1 and DBD-Dss1 domain of BRCA2 analyzed with FoldX measure homologous DNA recombination (HDR) activity significantly better than Rosetta and DDGun3D. Our study also revealed that there are other factors that contribute to predicting loss of function activity other than protein stability, with AlphaMissense ranking the best overall predictor of LOF activity in these tumor suppressor breast cancer genes.Author SummaryThe stability of a protein, often expressed in terms of Gibbs free energy (ΔΔG), is a critical factor in predicting loss of function (LOF) activity when a missense variant is present. The effect is higher in haploinsufficient genes like the tumor suppressor genesBRCA1, BRCA2, PALB2andRAD51C. Protein stability predictors that utilizes a wild-type structure to make its predictions is often limited by the availability of experimentally-derived protein structures. Here, in our study we show that generative AI, like AlphaFold2 (AF2) can predict structures similar to experimentally-derived structures with high similarity. Furthermore, protein stability tools such as FoldX, Rosetta, and DDGun3D can be used in conjunction to measure changes in stability. From our study, we find that complex AF2 structures representing the BRCT domain ofBRCA1and DBD domain ofBRCA2analyzed by FoldX predicts function significantly better than the experimentally-derived structures. However, predicted |ΔΔG| does not predict function better than purpose-builtin silicomissense predictors for protein function. Overall, we find the AlphaMissense is the best predictor to predict function in these tumor suppressor breast cancer genes.