
==== Front
Protein Sci
Protein Sci
10.1002/(ISSN)1469-896X
PRO
Protein Science : A Publication of the Protein Society
0961-8368
1469-896X
John Wiley & Sons, Inc. Hoboken, USA

10.1002/pro.5147
PRO5147
Research Article
Research Article
AlzDiscovery: A computational tool to identify Alzheimer's disease‐causing missense mutations using protein structure information
Pan et al.
Pan Qisheng 1 2
Parra Georgina Becerra 1 2
Myung Yoochan 1 2
Portelli Stephanie 1 2
Nguyen Thanh Binh 1 2
Ascher David B. https://orcid.org/0000-0003-2948-2413
1 2 d.ascher@uq.edu.au

1 The Australian Centre for Ecogenomics, School of Chemistry and Molecular Bioscience University of Queensland Brisbane Australia
2 Computational Biology and Clinical Informatics Baker Heart and Diabetes Institute Melbourne Australia
* Correspondence
David B. Ascher, The Australian Centre for Ecogenomics, School of Chemistry and Molecular Bioscience, University of Queensland, Brisbane, Queensland 4072, Australia.
Email: d.ascher@uq.edu.au

14 9 2024
10 2024
14 9 2024
33 10 10.1002/pro.v33.10 e514714 7 2024
25 3 2024
31 7 2024
© 2024 The Author(s). Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

Alzheimer's disease (AD) is one of the most common forms of dementia and neurodegenerative diseases, characterized by the formation of neuritic plaques and neurofibrillary tangles. Many different proteins participate in this complicated pathogenic mechanism, and missense mutations can alter the folding and functions of these proteins, significantly increasing the risk of AD. However, many methods to identify AD‐causing variants did not consider the effect of mutations from the perspective of a protein three‐dimensional environment. Here, we present a machine learning‐based analysis to classify the AD‐causing mutations from their benign counterparts in 21 AD‐related proteins leveraging both sequence‐ and structure‐based features. Using computational tools to estimate the effect of mutations on protein stability, we first observed a bias of the pathogenic mutations with significant destabilizing effects on family AD‐related proteins. Combining this insight, we built a generic predictive model, and improved the performance by tuning the sample weights in the training process. Our final model achieved the performance on area under the receiver operating characteristic curve up to 0.95 in the blind test and 0.70 in an independent clinical validation, outperforming all the state‐of‐the‐art methods. Feature interpretation indicated that the hydrophobic environment and polar interaction contacts were crucial to the decision on pathogenic phenotypes of missense mutations. Finally, we presented a user‐friendly web server, AlzDiscovery, for researchers to browse the predicted phenotypes of all possible missense mutations on these 21 AD‐related proteins. Our study will be a valuable resource for AD screening and the development of personalized treatment.

AlphaFold2
Alzheimer's disease
machine learning
missense mutation
protein structure
web server
National Health and Medical Research Council 10.13039/501100000925 GNT1174405 source-schema-version-number2.0
cover-dateOctober 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:14.09.2024
Pan Q , Parra GB , Myung Y , Portelli S , Nguyen TB , Ascher DB . AlzDiscovery: A computational tool to identify Alzheimer's disease‐causing missense mutations using protein structure information. Protein Science. 2024;33 (10 ):e5147. 10.1002/pro.5147

Review Editor: Nir Ben‐Tal
==== Body
pmc1 INTRODUCTION

Alzheimer's disease (AD), one of the most common forms of dementia and neurodegenerative diseases, has affected over 50 million elderly people in the last century (Knopman et al., 2021). This number is expected to double every 5 years and is projected to exceed 150 million by 2050 (Knopman et al., 2021; Scheltens et al., 2021), bringing a substantial burden on not only the AD‐affected individuals but also their families and society. Two main subcategories were identified based on an age cutoff of 65 years: early‐onset AD (EOAD), accounting for 5%–10% of patients, and late‐onset AD, accounting for 90%–95% of patients (Panegyres & Chen, 2013). AD patients usually first suffer from memory loss, and the symptoms can progressively deteriorate to cognitive disorder, problems with language, and change of personality (Ballard et al., 2011; Dubois et al., 2016; Wattmo et al., 2016). However, there is no cure for this complex disease, despite the expensive treatments to improve the symptoms.

Missense variants are one of the main sources of the dysfunctions of AD‐related proteins. Many rare variants of the amyloid precursor protein (APP), presenilin 1 (PSEN1), and presenilin 2 (PSEN2) increased the ratio of Aβ42/Aβ40 producing higher cellular toxicity (Kelleher 3rd & Shen, 2017). The loss‐of‐function mutations in triggering receptor expressed on myeloid cells‐2 (TREM2) in microglia could also increase the risk of AD up to 4‐fold (Hall‐Roberts et al., 2020; Jin et al., 2014), similar to the risk increment of heterozygous ε4‐apolipoprotein E (APOE) (Gratuze et al., 2018). However, considering the multifactorial pathogenic mechanism of AD, limited resources are available to understand the effects of the mutations of these genes and their phenotypes toward AD. A more comprehensive review was presented in the Supporting Information (Table S1).

Substantial computational methods have been established to estimate the risk of AD on variants since the last decade. The conventional, conservation‐based methods (Adzhubei et al., 2013; Choi & Chan, 2015; Hecht et al., 2015; Ng & Henikoff, 2003) have been used for annotating clinically reported mutations (Zhang et al., 2020), but these methods have poor specificity of diseases. The large deep mutational scanning data (Fowler et al., 2023; Notin et al., 2023) and the deep learning architecture have empowered more sophisticated methods (Brandes et al., 2023; Cheng et al., 2023; Gray et al., 2018; Jagota et al., 2023; Munro & Singh, 2021; Wu et al., 2021) in the recent 5 years, but many of them have limited focus on the impact of variants on protein structural environment. The latest AD‐specific tool, Alz‐Disc (Kulandaisamy et al., 2023), showed promising results in identifying disease‐causing mutations in different AD‐related proteins. This method, however, had no consideration of structural information, and their web server did not provide sufficient interpretation of the annotation of the missense variants.

We have previously shown that using computational tools could help improve the characterization of missense mutations and the prediction of disease phenotypes (Pan et al., 2023). In the current work, we implemented a similar machine learning‐based analysis (Portelli et al., 2021, 2023) on the missense mutations of 21 AD‐related proteins. Qualitative comparisons showed that disease‐causing mutations located in the buried regions tend to have a destabilizing effect on the protein structures. Considering these insights, we initially established a generalized model with comparable performance to current state‐of‐the‐art methods and further improved the performance at the protein‐specific level by tuning the weights in the training process. Our final models also achieved high performance on a clinical validation set, suggesting potential clinical utility. Finally, we developed a user‐friendly web tool to present all the possible missense mutations in these 21 AD‐related proteins. Our work proposed pathogenic molecular drivers toward AD, offering an invaluable resource for a better understanding of the pathology of AD.

2 RESULTS

The mutation analysis pipeline to identify AD‐causing variants was first adopted from previous works (Portelli et al., 2021, 2023) and was further scaled up for 21 AD‐related proteins, whose workflow is depicted in Figure 1. To simplify the name of the protein involved, we used the Gene symbol to refer to the proteins, which could be looked up in the Supporting Information (Table S1). The primary dataset for model training and validation was adapted from previous work with correction of labels of mutations and an independent clinical validation set was curated through ClinVar (Landrum et al., 2014) and literature reviews. We modeled the structures of the target proteins using AlphaFold2 (Jumper et al., 2021) and computed various features to describe the effects of mutations on both sequences and structures using computational biophysical measurements, followed by a qualitative comparison. We then performed a machine learning‐based analysis to accurately characterize these mutations.

FIGURE 1 Mutation analysis pipeline to identify Alzheimer's disease (AD)‐causing variants. Missense mutation data were first collected from multiple databases with their labels. Structural information was provided by AlphaFold2 models due to the unavailability of experiment‐determined structures. Then, a wide range of features were computed to describe the effect of mutations. To link these effects with the phenotypes, both qualitative comparisons and the development of machine learning models were implemented, followed by model assessment. Finally, AlzDiscovery, a user‐friendly, interactive web platform, was built for browsing the predicted phenotypes of all the possible missense mutations of each of these proteins.

2.1 Distribution of missense variants of mutation datasets

In the primary dataset, we curated a balanced dataset (n = 680) on 21 AD‐related proteins with a similar distribution of disease‐causing mutations (n = 320) and the neutral (non‐AD‐causing) counterparts (n = 360) (Table S2, Figure 2). All 21 targets, which are involved in different pathological pathways to AD (Table S1), have different sequences and structures (Figure 2c,d), suggesting a different environment for mutations, except PSEN1 and PSEN2, which share over 60% sequence identity. However, there is a large imbalanced distribution at the protein‐specific levels. Nine out of 21 proteins did not have both phenotypes, eight of which only contained the mutations labeled as neutral. Meanwhile, over 78% of the pathogenic mutations were found in PSEN1 (Figure 2). The other two familial AD (fAD)‐related genes, APP and PSEN2, showed balanced distributions across the different phenotypes.

FIGURE 2 Distribution of the neutral (a) and disease‐causing (b) mutations from 21 Alzheimer's disease (AD)‐related proteins in the primary dataset. Sequence identity (c) and structure similarity measured by all‐atom root mean squared deviation (RMSD) (d) between 21 AD‐related proteins in the dataset are also presented.

As for the clinical validation, a dataset with 125 mutations was curated on 11 proteins from the primary dataset, with most mutations labeled as neutral (88%) (Table S3). In addition, this validation set only contained two mutations from PSEN1, whose distribution is completely different from the one in the primary dataset, offering a more challenging clinical validation set for the final model.

2.2 Exploring molecular drivers leading to Alzheimer's disease

When comparing the biophysical measurements between disease‐causing and neutral mutations, a qualitative comparison was first performed on the protein structural environment, to investigate the potential impact of mutations. We noticed that the pathogenic mutations are predominantly located in the protein core, while neutral mutations are more solvent exposed (RSA: p‐value = 3.77 × 10−17, residue depth: p‐value = 9.33 × 10−16, Figure 3a,b). In addition, there are differences in secondary structure at the mutation site (Figure 3), particularly on the helices (p‐value = 1.16 × 10−18), as around 50% of target proteins are transmembrane proteins. Disease‐causing mutations are enriched in the structured regions, while the neutral ones are in the disordered regions (IUPred short score: p‐value = 2.26 × 10−43, IUPred long score: p‐value = 3.83 × 10−39, Figure 3f,g). Combining these differences in the structural environment at the mutation site, we suspect the pathogenic mutations tend to be located in the buried, helical regions of the protein structure, and they may alter the hydrophobic environment, which could cause a drastic effect on protein folding, such as the change of protein stability.

FIGURE 3 Qualitative comparisons on the protein structural environment to identify molecular drivers leading to AD. The two‐tailed Wilcoxon Sign‐rank tests were performed to compare the structural properties at the mutation site, including relative solvent accessibility (RSA) (a), residue depth (b), secondary structure type (c–e), protein disorderness (f–g), residue confidence score from AlphaFold2 (h), and torsion angle of peptide bonds (i, j), at 5% significance level.

We therefore compared and noticed that there is a significant difference in the change of protein stability between neutral and disease‐causing mutations (p‐value <0.002, Figure 4a–d), consistent across four predictors, namely mCSM‐Stability (Pires et al., 2014), SAAFEC‐SEQ (Li et al., 2021), DDMut (Zhou et al., 2023), and DynaMut2 (Rodrigues et al., 2021). The AD‐causing group displays a stronger destabilizing effect, containing more mutations with negative ΔΔG values. Further to that, we structurally analyzed the mutations in those three fAD‐related proteins (APP, PSEN1, and PSEN2) as a case study, and confirmed that a majority of mutations cause destabilizing effects (APP: 49.12%, PSEN1: 73.12%, and PSEN2: 52.78%), suggesting the change of stability upon mutations as a pathogenic molecular driver (Figure 4e–j), consistent with previous experimental findings (Grant et al., 2007; Rezaei‐Ghaleh et al., 2023; Roher et al., 1993). Our results provide a hypothesis that the missense mutations in APP monomers carrying destabilizing effects may introduce more protein aggregation of Aβ peptides to maintain protein stability.

FIGURE 4 Pathogenic molecular drivers leading to AD. Qualitative analyses were performed to compare the effect of mutations on the change of protein stability computed by mCSM‐Stability (a), SAAFEC‐SEQ (b), DDMut (c), and DynaMut2 (d), measured by ΔΔG (Kcal/mol) with zero as a cutoff (ΔΔG > 0: stabilizing; ΔΔG < 0: destabilizing). Qualitative analysis was implemented using the two‐tailed Wilcoxon Sign‐rank test at a 5% significance level. Detailed structural analysis was implemented to visualize the mutations with stabilizing/destabilizing effects on three familial AD‐related genes, APP (e and f), PSEN1 (g and h), and PSEN2 (i and j).

To further explore the intramolecular changes upon mutations, six representative mutations from the PSEN1 protein, which is well identified to be AD‐causing in the previous study (Do et al., 2023), were selected for further investigation (Figure 5). In general, these mutant structures present considerable differences in the structural environment. The change of biochemical properties caused by mutations introduced the loss of π–π interactions, especially in P117L, L166P, and L435F. The P177L and L286V were located in the flexible loop regions, which cause less impact on protein stability based on the prediction of DDMut (ΔΔGP117L = 0.05 Kcal/mol and ΔΔGL286V = 0.01 Kcal/mol, respectively), while the other four mutations (ΔΔGI143T = −0.07 Kcal/mol, ΔΔGL166P = −0.08 Kcal/mol, ΔΔGG384A = −0.18 Kcal/mol, and ΔΔGL435F = −2.11 Kcal/mol) could potentially disrupt protein folding, negatively affect the conformational changes to active state of γ‐secretase, and further reduce the activity of cleavage of APP, consistent with previous findings (Do et al., 2023; Somavarapu & Kepp, 2016).

FIGURE 5 Analysis of the intramolecular residue interaction of six representative mutations of PSEN1 protein, namely P117L, I143T, L166P, L286V, G384A, and L435F. Wild‐type residues were colored as cyan, while the mutant residues were colored as magenta. The hydrogen bonds were represented by red dashed lines. The polar interactions were represented by orange dashed lines. The π–π interactions were represented by green dashed lines.

2.3 Predictive performance to identify disease‐causing mutations using machine learning

Given the bias on the destabilizing effect of disease‐causing mutations, we tried to build a machine‐learning model to better characterize the pathogenic variants leading to AD. The primary dataset was first divided into training (60%), validation (20%), and blind test (20%) sets, ensuring non‐redundancy between datasets according to unique protein position. Seven machine learning algorithms were first employed to elucidate the association between different molecular features and phenotypes of mutations, and the algorithm with the best predictive performance was selected and further optimized using greedy feature selection (details in Section 4). Our initial generic model based on the Gradient Boosting algorithm achieved high performance on both 10‐fold cross‐validation (10‐CV) (area under the curve (AUC) = 0.94, Matthew's correlation coefficient (MCC) = 0.79) and the non‐redundant blind test (AUC = 0.94, MCC = 0.75) (Table S4). After configuring the sample weights in the training process to improve the attention for target proteins, our enhanced model obtained a higher overall performance on the blind test (AUC = 0.95, MCC = 0.79, Table 1). The consistent performance of the blind test on sensitivity (0.88), specificity (0.92), and precision (0.90) suggests our model was not biased toward a specific phenotype.

TABLE 1 Performance of AlzDiscovery of different validation.

Test type	BACC	F1‐score	MCC	Sensitivity a (recall)	Precision a	Specificity a	AUC	
10‐CV	0.89	0.88	0.79	0.86	0.91	0.92	0.95	
Blind test	0.90	0.89	0.79	0.88	0.90	0.92	0.95	
Clinical validation	0.68	0.41	0.32	0.47	0.37	0.89	0.70	
a In this work, the label “Disease‐causing” was considered as positive, while “Neutral” was considered as negative.

In addition to a small improvement in the generalized performance, the accuracy at the protein‐specific level for the target protein was significantly improved (Figure S1), especially on those proteins with small sample size of mutations in the primary dataset (alternative PRNP and NOS3). The performance on APOE has a small drop after weight tuning, while the performance on the other proteins remains the same. Our model fails to predict the pathogenic mutations in ADAM10, which is likely due to the extremely small number of samples (n = 2) in the dataset.

To assess the generalizability of our model at the protein‐specific level, a leave‐one‐protein‐out validation was conducted (Table S5). We noticed that there is a large overall performance deterioration, with AUC down to 0.65. This decrease in performance can be attributed to the imbalanced distribution of data within our primary dataset, as over 78% of all mutations were present in one protein, PSEN1. In excluding this protein during leave‐one‐protein‐out validation, our previous accuracy of 0.97 dropped to 0.21, showing the significant contribution of these mutations to the overall performance. Apart from this, in considering 21 different genes, different mechanisms contributing to AD have been assessed in our initial model. During leave‐one‐protein‐out, removing specific proteins could underrepresent minor mechanisms of AD, which would also be reflected in lower accuracy.

We then applied the enhanced model to the clinical validation set to assess model generalizability and the clinical utility, excluding APOE, for which we applied the generic model due to the performance drop obtained through the enhanced model. Despite an overall performance deterioration, our model still achieved a relatively reasonable performance (AUC = 0.70 and MCC = 0.32) (Table 3), indicating a certain degree of generalizability of the model. The model still achieved a relatively balanced capability on classifying two different phenotypes (sensitivity = 0.47 and precision = 0.37), but the high specificity (0.92) suggests a certain bias in the predictions toward Neutral phenotypes. At the protein‐specific level, the performance of some proteins, such as CSF1R, PRNP, and VCP, was improved by around 4%–22% (Figure S1D). However, there was a performance drop when predicting the mutations from PSEN2.

2.4 Performance comparison with the state‐of‐the‐art methods

When comparing our method with the current advanced variant effect predictors, we found that our tool achieved comparable performance on the blind test with the AD‐specific approach, Alz‐Disc, across different metrics (Table 2, Figure 6). Though the comparison of performance of the blind test between our method and Alz‐Disc may not be fair due to the redundancy caused by different divisions of training and test sets, our model achieved a much higher AUC (AUC = 0.95) than that of Alz‐Disc (AUC = 0.55), suggesting a higher confidence of the predicted phenotypes from our model. At the protein‐specific level, while the Alz‐Disc achieved a better accuracy on three proteins, APP, MAPT, and PRNP, our method is capable of accurately identifying the mutations from TREM2, a newly proposed AD biomarker, as we corrected the mutation labels of TREM2 in the primary dataset (Figure 5B). In terms of the independent validation, our model presented a much better performance (AUC = 0.70 and MCC = 0.32) in the clinical benchmark, compared with Alz‐Disc (AUC = 0.54 and MCC = 0.18) (Table 3). Particularly, our model could better classify mutations in well‐known biomarkers, such as APP, APOE, MAPT, and TREM2 (Table S6). As we trained our model by integrating structure‐based features, this may suggest the importance of considering the effects of mutation in a three‐dimensional environment to enhance the robustness.

TABLE 2 Performance comparison on the blind test of the primary dataset.

Method	Type	BACC	F1‐score	MCC	Sensitivity (recall)	Specificity	Precision	AUC	
AlzDiscovery	AD‐specific	0.90	0.89	0.79	0.88	0.92	0.90	0.95	
Alz‐Disc	AD‐specific	0.92	0.91	0.84	0.88	0.96	0.95	0.55	
DeMaSk	Deep mutation scanning‐based	0.73	0.71	0.47	0.69	0.77	0.74	0.78	
Envision	Deep mutation scanning‐based	0.57	0.64	0.16	0.83	0.31	0.52	0.60	
PROVEAN	Conservation‐based	0.73	0.74	0.46	0.80	0.66	0.68	0.77	
SIFT	Conservation‐based	0.66	0.67	0.32	0.72	0.59	0.62	0.51	
SNAP2	Conservation‐based	0.67	0.71	0.37	0.88	0.46	0.60	0.77	
AlphaMissense	Deep learning‐based	0.82	0.80	0.65	0.75	0.89	0.86	0.88	
ESM1b	Deep learning‐based	0.80	0.79	0.59	0.82	0.77	0.77	0.85	
CPT‐1	Deep learning‐based	0.84	0.83	0.68	0.80	0.87	0.85	0.87	

FIGURE 6 Predictive performance comparisons and model interpretation. Predictive performance of AlzDiscovery and the other variant effect predictors on the blind test is shown in the receiver operating characteristic (ROC) curve (a), and a performance at the protein‐specific levels is compared using accuracy due to the imbalance of mutation phenotypes of each protein (b). Feature importance of the AlzDiscovery model was interpreted using SHapley Additive exPlanations (SHAP) method (c). Features were ordered from top to bottom based on their contributions to the models. To interpret the plot, the contribution of features toward different labels is represented by the direction of the x‐axis (negative: disease‐causing; positive: neutral), while the numeric value of the feature is represented by the color of each data point (red: High; blue: Low). Combining these two aspects, the contributions of each feature could be interpreted.

TABLE 3 Performance comparison on the clinical validation.

Method	Type	BACC	F1‐score	MCC	Sensitivity (recall)	Precision	Specificity	AUC	
AlzDiscovery	AD‐specific	0.68	0.41	0.32	0.47	0.37	0.89	0.70	
Alz‐Disc	AD‐specific	0.61	0.30	0.18	0.40	0.24	0.83	0.54	
DeMaSk	Deep mutation scanning‐based	0.42	0.16	−0.11	0.53	0.10	0.31	0.64	
Envision	Deep mutation scanning‐based	0.48	0.20	−0.05	0.87	0.12	0.09	0.59	
PROVEAN	Conservation‐based	0.54	0.23	0.07	0.93	0.13	0.15	0.57	
SIFT	Conservation‐based	0.41	0.10	−0.12	0.2	0.07	0.62	0.61	
SNAP2	Conservation‐based	0.52	0.21	0.03	0.60	0.13	0.45	0.54	
AlphaMissense	Deep learning‐based	0.39	0.14	−0.15	0.47	0.09	0.32	0.40	
ESM1b	Deep learning‐based	0.53	0.22	0.04	0.73	0.13	0.33	0.50	
CPT‐1	Deep learning‐based	0.61	0.27	0.14	0.73	0.16	0.48	0.58	

Our enhanced model outperformed the conservation‐based methods (SIFT (Ng & Henikoff, 2003), SNAP2 (Hecht et al., 2015), and PROVEAN (Choi & Chan, 2015)) on the blind test (Table 2, Figure 6). Meanwhile, our approach also presents a better performance on both the blind test and independent validation set than the ones of deep mutation scanning (DMS)‐based methods and deep learning‐based methods, including the latest approach AlphaMissense (AUCblind‐test = 0.88 and AUCvalidation = 0.40) (Cheng et al., 2023). These generalized variant effect predictors were established to evaluate the fitness of the residue substitution, and therefore may not have enough specificity on neurodegenerative diseases such as AD. However, it should be emphasized that those new deep‐learning‐based methods, including AlphaMissense, ESM1b (Brandes et al., 2023), and CPT‐1 (Jagota et al., 2023), present a higher generalizability to identify AD‐causing mutations, which could be due to the use of feature vectors cooperated with protein structural information.

2.5 Model interpretation

Since both our generic and protein‐specific models obtained good performance, we then tried to explore the contributions of each feature to gain a better understanding of the instrumental factors to identify AD‐causing mutations (Figure 6c). Our models leverage 12 features describing the mutation from different aspects. The feature with the highest contributions is the neighboring aliphatic amino acid frequency (aliphatic_freq), consistent with one of the features of Alz‐Disc (Kulandaisamy et al., 2023). Combining this with another three features, namely the structural polar contact at theWT mutation site (ITPolar), the difference of structural weak polar contact (ITd_WeakPolar), and the hydrophobic mutant (hydrophobic_mt), we suspect that the increment of hydrophobicity at the mutation site may not only impact protein stability but promote protein aggregation, consistent with previous studies (Dasuri et al., 2010; Kim & Hecht, 2006). There are five other features derived from AAindex (MEIH800101_delta_2, VASM830102_delta_3_norm, CHAM830107_mt, KLEP840101_mt, and FAUJ880112_mt), which account for different biochemical properties of 20 natural amino acids, namely the change of hydrophobic behavior in two neighboring amino acids (Meirovitch et al., 1980), the change of computed conformational state in neighboring region (Vasquez et al., 1983), and three properties of mutant amino acids (the intermolecular force and steric parameters (Charton & Charton, 1983), the net charge (Klein et al., 1984), and the negative charged (Klein et al., 1984), respectively). Two structural features from graph‐based signatures (Aro:Sul‐11.00 and Neg:Sul‐6.50) explain the structural environment of atoms with aromatic, sulfur, and negative pharmacophores, while the feature Neg represents the difference of negative pharmacophores between WT and mutant amino acids.

In brief, these features mainly inform the importance of the hydrophobic environment at the mutation site, the biochemical properties in the neighboring region, and the change of polar residue contacts. The change of hydrophobicity at the mutation environment could be particularly crucial to the aggregation of Aβ and the formation of amyloids, agreeing with our qualitative analysis and other previous works (Sgourakis et al., 2007; van Gils et al., 2020).

2.6 Saturation mutagenesis and web server

Finally, we applied our final models with the best performance at the protein‐specific level to generate the predictions of all possible missense mutations for these 21 AD‐related proteins. To be specific, we applied the enhanced protein‐specific models to all proteins for in silico saturation mutagenesis, excluding APOE and PSEN2 because of the performance deterioration of their weighted model, on which we applied the initial generic model. Here we present the mutation landscape of three fAD‐related proteins, APP, PSEN1, and PSEN2 (Figure 7), and three strong risk factors, APOE, MAPT, and TREM2 (Figure S2), while the rest are available on the AlzDiscovery web server https://biosig.lab.uq.edu.au/alzdiscovery/.

FIGURE 7 Mutation landscape leading to AD of three familial AD‐related proteins, APP, PSEN1, and PSEN2. Protein structures are colored based on two color schemes, namely the average probability leading to diseases of all possible missense mutations at each mutation site (a–c) and the change of protein stability upon mutations at the mutation site (ΔΔG, d–f), while the heatmaps present detailed results of different mutants at each site of three proteins (g–i).

We initially noticed that the pathogenic variants are enriched in the region where mutations tend to destabilize protein folding (Figures 7a–f and S2A–F), consistent with the qualitative analyses of computational biophysical measurements. Following detailed analysis of the intramolecular residue contacts in the fAD‐related genes (Figure S3) suggests that the destabilizing effect is due to the change in interaction profiles, such as the alteration of hydrophobic and polar contacts. Our structure‐based model, however, may not be capable of a thorough explanation of all mutations (Figure S3A) due to the complexity of AD. A similar analysis on the change of residue contacts is found through the web server.

The pathogenic area tends to form clusters at some specific mutation sites, where most residue substitutions could cause pathogenicity. This is mainly because a number of features in the models are capturing the residue environment at the mutation loci, rather than a more detailed change between WT and mutant amino acid. These selections of features are associated with the bias of mutation distribution in the dataset, as there are not enough mutations at the same mutation site of a protein with different phenotypes, resulting in the model trying to identify the disease‐causing regions instead of the disease‐causing mutations. In addition, there is a bias to neutral phenotypes of those mutations toward the negative charge amino acids, particularly obvious in APP, PSEN2, and TREM2 (Figures 7g,j and S2F), due to the use of feature FAUJ880112 (Figure 6c). The limited mutation types may also lead to a neutral bias on negative charge mutants. In terms of APOE and TREM2, disease‐causing mutations are mainly clustered at the helical regions (Figure S2), consistent with our results on qualitative analysis. A different pattern was noticed in the tau pathology. The mutation landscape on the MAPT shows only a relatively small number of pathogenic mutations clustered at mutation site P587 (Figure S2B). This indicates that missense mutations on tau protein may not be the main drivers of AD but other dementia, such as frontotemporal dementia with Parkinsonism linked to chromosome 17 (FTDP‐17), consistent with experimental studies (Alonso Adel et al., 2004; Bhaskar et al., 2005).

3 DISCUSSION

The complicated pathogenic mechanism of AD has been significantly studied over the last century, which manifests different biomarkers, including the accumulation of Aβ peptides, tau pathology, dysregulation of lipidology, and the activation of innate immune reaction in the human central nervous system. Missense mutations on these AD‐related proteins have been reported as one of the key pathogenic drivers (Bossaerts et al., 2022; Brownjohn et al., 2018; Sherrington et al., 1996; Wakutani et al., 2004). In the current work, we have characterized the effects of mutations on 21 AD‐related proteins with a structural insight using computational biophysical methods and machine learning. Our model achieved a relatively robust performance on different benchmarks, and the interpretation of the models revealed that missense mutations could alter the hydrophobic environment, protein stability, and polar residue interactions on AD‐related proteins, which could be important molecular drivers leading to disease. A user‐friendly and interactive web platform was established for all researchers to browse the predicted phenotypes and visualize the changes in the protein structures.

We propose two important mutational risk factors, namely the hydrophobic environment at the mutation site and the changes in protein stability caused by mutation. These two factors are both associated with the protein intact folding, which is crucial for protein functions. The authors of Alz‐Disc also proposed a similar idea (Kulandaisamy et al., 2023) by interpreting the neighboring amino acid frequency, but we directly used the biophysical measurements on the protein structures and further supported this hypothesis. In addition to the effect on a monomeric structure, the change in hydrophobicity and protein stability are also key factors in protein aggregation. We could have implemented more biophysical measurements, such as the protein–protein interaction of Aβ peptides, and the lipid binding of APOE (Frieden et al., 2017) and TREM2 (Wang et al., 2015). These in‐depth analyses, however, were hindered by the complexity and scale of this work.

By improving the implementation of the mutation analysis pipeline in the previous work (Portelli et al., 2021; Portelli et al., 2023), we managed to investigate the effect of missense mutations and predict the disease‐related phenotypes on multiple proteins, by tuning the sample weights of the mutation from a target protein. This method could further improve the performance at the protein‐specific level, but we also noticed the performance deterioration on two proteins, APOE and PSEN2, due to the bias of features, such as the biased representation of negative‐charged mutant amino acids. We therefore use the generic models to perform saturation mutagenesis for these two proteins. This is because while we believe the weight‐tuning approach could improve the predictive performance (Fang et al., 2022), we also aimed to give the best predictions on these mutations for clinical use.

Further to this, there were three limitations in this work that were faced due to the lack of variant data availability. Firstly, the imbalanced distribution of mutations across proteins studied in the primary dataset resulted in a bias of our model, which is exemplified by the performance deteriorations on the leave‐one‐protein‐out test (Table S5). Similarly, the insufficient coverage of different types of residue substitutions led to a relative deficiency in performance, especially the performance in predicting the mutations toward negative‐charged amino acids (Figure 7g–i). Meanwhile, there is no data to describe the interactions among the mutations and AD‐related genes, which means that our model did not capture the complex gene regulation network associated with the pathogenesis of AD. Therefore, the mutation analysis toward AD can be improved when more high‐quality data becomes available in the future.

Understanding the potential risk of AD could be useful for patient monitoring in the early stages, especially for EOAD and fAD. Clinical researchers could easily access our web tool and retrieve the results from computational saturation mutagenesis. Further to that, an interactive structure viewer is available to present the residue interactions, the regions with decreased protein stability, and other biochemical properties. Thus, we believe that our work could be an invaluable resource for not only further AD research but also clinical application. Personalizing treatment could be designed by leveraging the mutation information from our saturation mutagenesis, as a pioneering study has shown the feasibility of using small‐molecule drugs to stabilize Aβ (Lohr et al., 2022).

4 MATERIALS AND METHODOLOGY

4.1 Data curation and protein structural information

4.1.1 Primary dataset

We leveraged the dataset curated by previous work (Kulandaisamy et al., 2023) and implemented data cleaning. According to the previous work, mutations associated with AD from multiple resources, including Alzforum (Kinoshita & Clark, 2007), Humsavar from UniProt (Apweiler et al., 2004), ClinVar (Landrum et al., 2018), dbSNP (Smigielski et al., 2000), and HuVarBase (Ganesan et al., 2019), were extracted and labeled as disease‐causing and neutral respectively according to the annotations from the databases. We next performed a data cleaning on this original dataset, which involved the removal of synonymous mutations, and missense mutations that did not map to UniProt canonical sequences. All mutations located in the fAD‐related genes (APP, PSEN1, and PSEN2) and AD risk factors (APOE, MAPT, and TREM2) were further manually checked to ascertain their reliability. Notably, all TREM2 mutations were labeled as Neutral in the original dataset, so we carefully corrected the mutation phenotypes by checking the contradictory labels in these resources. In the final dataset, six TREM2 mutations (R47H, R62H, T66M, R136Q, H157Y, and T223I) were marked as Disease‐causing, as reported in the literature (Jin et al., 2014; Sayed et al., 2021; Sirkis et al., 2016). After data cleaning, the final dataset consisted of 680 mutations on 21 AD‐related proteins (available in the Supplementary Files), which was used for model development and optimization.

4.1.2 Clinical validation

To assess the clinical utility of our model, missense mutations of these 21 AD‐related proteins were collected from ClinVar via precise phenotyping based on clinical conditions from literature or the American College of Medical Genetics and Genomics (ACMG) guidelines (Richards et al., 2015). Redundant mutations between the primary dataset and clinical validation set were removed. The final clinical validation set contains 125 mutations on 11 proteins, which is available in the Supplementary Files.

4.1.3 Structure curation

Protein structural models were generated to improve the performance of identifying pathogenic mutations toward AD. Due to the unavailability of the experiment‐determined structure of most proteins of interest, protein structures were generated by the latest AlphaFold2 (Jumper et al., 2021) using the canonical protein sequences with template date 2022‐03‐02. Models with the highest confidence score (predicted Local Distance Difference Test, pLDDT (Mariani et al., 2013)) were selected for feature generations.

4.2 Feature generations

Numerous sequence‐ and structure‐based features were computed to annotate the effects of mutations in various aspects, including biochemical properties, protein biophysical properties, local residue environment, and functionality assessment.

4.2.1 Biochemical properties

To annotate the wild‐type (WT) and mutant amino acids, all 20 amino acids were assigned to one of the five groups based on their biochemical properties, including hydrophobic (A, F, I, L, M, V, W, and Y), polar (N, Q, S, and T), negatively charged (D and E), positively charged (H, K, and R), and special (C, G, and P). The properties from AAindex 1 (Kawashima et al., 1999) and three basic properties, isoelectric point (pI), molecular weight, and molecular volume (Annu Rev Biophys Bioeng, 1984; Zamyatnin, 1972), were also included to describe both the WT and mutant residues.

To annotate the mutations, the biochemical scores from multiple amino acid substitution matrices and the statistical interpretations of protein contact potentials from the biochemical databases, AAindex 2 and 3 (Kawashima et al., 1999), were extracted based on the amino acid substitutions. To consider the effect on position‐specific conservation, scores were extracted from the Position Specific Substitution Matrix generated by PSI‐BLAST (Altschul et al., 1997; Camacho et al., 2009) against nr databases (Sayers et al., 2022) with three iterations.

To further consider the neighboring environment at the mutation site, the difference of the properties from AAindex 1 between WT and mutant amino acid was calculated, whose idea originated from previous work (Kulandaisamy et al., 2023). A sliding window of one, three, five, and seven residues at the mutation site was considered. WT properties were the average properties of the amino acids inside the window, while the mutant properties were unprocessed and extracted from the resources. This could consider the changes in biochemical properties between the mutant and the surrounding WT. Both the raw difference and the one normalized by the length of the protein sequence were included as features.

4.2.2 Protein biophysical properties

The mutation effect on protein stability was computed by a sequence‐based method, SAAFEC‐SEQ (Li et al., 2021), and three structure‐based methods, mCSM‐Stability (Pires et al., 2014), DynaMut2 (Rodrigues et al., 2021), and DDMut (Zhou et al., 2023), using default settings. The reliability of these methods on AlphaFold2 models has been systematically studied in a previous work (Pan et al., 2022). All the calculations were measured using ΔΔG (Kcal/mol), with zero as a cutoff (ΔΔG >0: stabilizing; ΔΔG <0: destabilizing; and ΔΔG = 0: neutral).

4.2.3 Local residue environment

The local environment at the mutation site was described from three aspects, namely neighboring amino acid frequency, general structural properties, and the residue contacts and atomic distance pattern, to identify local differences between neutral and disease‐causing mutations.

Neighboring amino acid frequency

To model the sequence‐based residue environment, the neighboring amino acid frequency was computed by considering a 15‐residue window at the mutation site (Kulandaisamy et al., 2023). All 20 amino acids were assigned into another seven different groups based on previous works (Caldararu et al., 2020), namely hydrophobic (A, L, and M), aliphatic (I and V), aromatic (F, Y, and W), long polar (E, Q, K, and R), short polar (H, S, T, and C), short charged/polar (D and N), and structure‐breaking (G and P). The count of amino acids in each group of the seven leading and following residues was computed and normalized using the length of the protein sequence to generate the frequency. This feature could be used to represent the biochemical property of the mutation environment.

General structural properties

Relative solvent accessibility (RSA) and residue depth were measured via Biopython packages (Cock et al., 2009) to quantify solvent exposure in the mutation site. Secondary structure and torsion angles of peptide bonds (phi and psi) were computed using the DSSP program (Joosten et al., 2011; Kabsch & Sander, 1983). To quantify the disorderliness of the protein structure, both predictions of IUPred (Dosztanyi, 2018; Erdos et al., 2021) and the pLDDT scores from the AlphaFold2 model were introduced.

Residue contacts and atomic distance pattern

Local residue contact at the mutation site of both WT and mutant structures (generated by MODELLER Marti‐Renom et al., 2000; Sali & Blundell, 1993; Webb & Sali, 2016) was computed using Arpeggio packages (Jubb et al., 2017). The atomic distance pattern at the mutation site was captured by graph‐based signatures (Pires et al., 2014). In graph‐based signatures, atoms with eight different pharmacophores are modeled as nodes and their contacts within a certain distance cutoff are modeled as edges. The residue environmental context is described using a cumulative distribution of different types of contacts via different configurations of distance steps and cutoffs from the mutation site.

4.2.4 Functionality assessment

Two conventional variant effect predictors, Sorting Intolerant From Tolerant (SIFT) (Ng & Henikoff, 2003) and SNAP2 (Hecht et al., 2015), were used to estimate the impact of mutations on protein function based on conservation.

4.3 Machine learning analysis

To build the generalized model including all mutations in different proteins, the primary dataset was first divided into training (60%), validation (20%), and blind test (20%) sets, ensuring non‐redundancy between data sets according to unique protein position. Seven machine learning algorithms, including Decision Trees, Random Forest, Extra Trees, Adaptive Boosting, Gradient Boosting, Extreme Gradient Boosting, and Logistic Regression, were trained and tested on the training and validation sets first (Table S7), followed by the optimization using greedy feature selection (Figure S4A), whose algorithm has been well described in previous works (Myung et al., 2020; Rodrigues et al., 2021). After hyperparameter tuning using HalvingRandomSearch, three boosting‐based and tree‐based parameters were set, namely n_estimators = 80, learning_rate = 0.1, and max_depth = 7 (Figure S4B, Supplementary File halving_search.csv). Predictive performance was evaluated using accuracy, balanced accuracy (BACC), recall, precision, F1‐score, MCC, and area under the receiver operating characteristic curve (AUC).

To further enhance the capability of identifying AD‐causing mutations at the protein‐specific levels, the sample weights were tuned during the training process for 21 AD‐related proteins (Fang et al., 2022). Specifically, mutations from the target protein were assigned with weights = 11.0, while the mutations from non‐target proteins were assigned with weights = 1.0, by setting the “sample weights” of the training process in the Scikit‐Learn packages. This choice of weight was generated by optimizing the generalized performance. Protein‐specific performance was evaluated using accuracy, as not all the target proteins had both neutral and disease‐causing mutations.

Predictive performance of both our enhanced models was compared with different state‐of‐the‐art methods, including an AD‐specific method, Alz‐Disc (Kulandaisamy et al., 2023), two DMS‐based methods, namely DeMaSk (Munro & Singh, 2021) and Envision (Gray et al., 2018), three conventional predictors, including SIFT (Ng & Henikoff, 2003), SNAP2 (Hecht et al., 2015), and PROVEAN (Choi & Chan, 2015), and three deep learning‐based methods, which are AlphaMissense (Cheng et al., 2023), ESM1b (Brandes et al., 2023), and CPT‐1 (Jagota et al., 2023), respectively. A cutoff was used to convert the fitness score into disease‐causing and neutral phenotypes. The suggested cutoffs of 0.05, 0.0, and −2.5 were used for SIFT, SNAP2, and PROVEAN, respectively. In terms of the DMS‐based and deep learning‐based methods, we used the training set to determine a cutoff to maximize the performance of MCC, which has been described in previous work (Pan et al., 2023), and the final cutoffs used in this work were illustrated here (DeMaSk: −0.22, Envision: 0.97, AlphaMissense: 0.61, ESM1b: −8.18, CPT‐1: 0.39).

4.4 Web server

AlzDiscovery has been implemented as a user‐friendly interactive web server freely available to the research community at: https://biosig.lab.uq.edu.au/alzdiscovery/. Users can browse the effects of pathogenicity of all the possible missense mutations of 21 AD‐related proteins by choosing the mutation of interest. Detailed structure information at the mutation site, such as interatomic interactions, is displayed on the 3D viewers (Rose et al., 2018). Predicted phenotypes of saturation mutagenesis using our best performing models at the protein‐specific levels can be downloaded from the website for further analysis and patient monitoring.

AUTHOR CONTRIBUTIONS

Qisheng Pan: Methodology; formal analysis; writing – original draft; visualization; validation. Georgina Becerra Parra: Methodology; data curation; formal analysis; writing – review and editing. Yoochan Myung: Software; writing – review and editing. Stephanie Portelli: Methodology; formal analysis; writing – review and editing. Thanh Binh Nguyen: Methodology; formal analysis; writing – review and editing. David B. Ascher: Methodology; writing – review and editing; supervision; conceptualization; formal analysis.

CONFLICT OF INTEREST STATEMENT

The authors declared no conflict of interest.

Supporting information

Appendix S1. Introduction of 21 AD‐related proteins and their mechanism toward AD (Table S1). Distributions of mutation data in the primary dataset (Table S2) and validation set (Table S3). Performance of development of the initial generic models and enhanced models after weight tuning (Table S4, Table S7, Figure S1, Figure S4). Performance of the leave‐one‐protein‐out evaluation (Table S5). Comparison of performance on the protein‐specific level between AlzDiscovery and Alz‐Disc in the independent validation (Table S6). Mutation landscape leading to AD of three strong risk factors, APOE, MAPT, and TREM2 (Figure S2). Analysis of intramolecular residue contacts of three mutations in APP, PSEN1, and PSEN2, respectively (Figure S3).

Appendix S2. Primary dataset in Excel spreadsheet. Independent validation set in Excel spreadsheet; All AlphaFold2 models analyzed in this work. Summary of HalvingRandomSearch of hyperparameter tuning in Excel spreadsheet.

ACKNOWLEDGMENTS

The authors thank the funding of Research Scholarship from University of Queensland (to Q.P and G.B.P) supported by the Investigator Grant from the National Health and Medical Research Council (NHMRC) of Australia (GNT1174405 to D.B.A); Victorian Government's Operational Infrastructure Support Program. Open access publishing facilitated by The University of Queensland, as part of the Wiley ‐ The University of Queensland agreement via the Council of Australian University Librarians.
==== Refs
REFERENCES

Adzhubei I , Jordan DM , Sunyaev SR . Predicting functional effect of human missense mutations using PolyPhen‐2. Curr Protoc Hum Genet. 2013; Chapter 7: Unit7 20. PMID: 23315928 {Medline}.
Alonso Adel C , Mederlyova A , Novak M , Grundke‐Iqbal I , Iqbal K . Promotion of hyperphosphorylation by frontotemporal dementia tau mutations. J Biol Chem. 2004;279 :34873–34881. PMID: 15190058 {Medline}.15190058
Altschul SF , Madden TL , Schaffer AA , Zhang J , Zhang Z , Miller W , et al. Gapped BLAST and PSI‐BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25 :3389–3402. PMID: 9254694 {Medline}.9254694
Amino acid, peptide, and protein volume in solution. Annu Rev Biophys Bioeng. 1984;13 :145–165. PMID: 6378067 {Medline}, 191, 219.6378067
Apweiler R , Bairoch A , Wu CH , Barker WC , Boeckmann B , Ferro S , et al. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2004;32 :D115–D119. PMID: 14681372 {Medline}.14681372
Ballard C , Gauthier S , Corbett A , Brayne C , Aarsland D , Jones E . Alzheimer's disease. Lancet. 2011;377 :1019–1031. PMID: 21371747 {Medline}.21371747
Bhaskar K , Yen SH , Lee G . Disease‐related modifications in tau affect the interaction between Fyn and tau. J Biol Chem. 2005;280 :35119–35125. PMID: 16115884 {Medline}.16115884
Bossaerts L , Hendrickx Van de Craen E , Cacace R , Asselbergh B , van Broeckhoven C . Rare missense mutations in ABCA7 might increase Alzheimer's disease risk by plasma membrane exclusion. Acta Neuropathol Commun. 2022;10 :43 PMID: 35361255 {Medline}.35361255
Brandes N , Goldman G , Wang CH , Ye CJ , Ntranos V . Genome‐wide prediction of disease variant effects with a deep protein language model. Nat Genet. 2023;55 :1512–1522. PMID: 37563329 {Medline}.37563329
Brownjohn PW , Smith J , Solanki R , Lohmann E , Houlden H , Hardy J , et al. Functional studies of missense TREM2 mutations in human stem cell‐derived microglia. Stem Cell Rep. 2018;10 :1294–1307. PMID: 29606617 {Medline}.
Caldararu O , Mehra R , Blundell TL , Kepp KP . Systematic investigation of the data set dependency of protein stability predictors. J Chem Inf Model. 2020;60 :4772–4784. PMID: 32786698 {Medline}.32786698
Camacho C , Coulouris G , Avagyan V , Ma N , Papadopoulos J , Bealer K , et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10 :421 PMID: 20003500 {Medline}.20003500
Charton M , Charton BI . The dependence of the Chou‐Fasman parameters on amino acid side chain structure. J Theor Biol. 1983;102 :121–134. PMID: 6876837 {Medline}.6876837
Cheng J , Novati G , Pan J , Bycroft C , Zemgulyte A , Applebaum T , et al. Accurate proteome‐wide missense variant effect prediction with AlphaMissense. Science. 2023;381 :eadg7492 PMID: 37733863 {Medline}.37733863
Choi Y , Chan AP . PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics. 2015;31 :2745–2747. PMID: 25851949 {Medline}.25851949
Cock PJ , Antao T , Chang JT , Chapman BA , Cox CJ , Dalke A , et al. Biopython: freely available python tools for computational molecular biology and bioinformatics. Bioinformatics. 2009;25 :1422–1423. PMID: 19304878 {Medline}.19304878
Dasuri K , Ebenezer P , Zhang L , Fernandez‐Kim SO , Bruce‐Keller AJ , Markesbery WR , et al. Increased protein hydrophobicity in response to aging and Alzheimer disease. Free Radic Biol Med. 2010;48 :1330–1337. PMID: 20188163 {Medline}.20188163
Do HN , Devkota S , Bhattarai A , Wolfe MS , Miao Y . Effects of presenilin‐1 familial Alzheimer's disease mutations on gamma‐secretase activation for cleavage of amyloid precursor protein. Commun Biol. 2023;6 :174 PMID: 36788318 {Medline}.36788318
Dosztanyi Z . Prediction of protein disorder based on IUPred. Protein Sci. 2018;27 :331–340. PMID: 29076577 {Medline}.29076577
Dubois B , Hampel H , Feldman HH , Scheltens P , Aisen P , Andrieu S , et al. Preclinical Alzheimer's disease: Definition, natural history, and diagnostic criteria. Alzheimers Dement. 2016;12 :292–323. PMID: 27012484 {Medline}.27012484
Erdos G , Pajkos M , Dosztanyi Z . IUPred3: prediction of protein disorder enhanced with unambiguous experimental annotation and visualization of evolutionary conservation. Nucleic Acids Res. 2021;49 :W297–W303. PMID: 34048569 {Medline}.34048569
Fang M , Su Z , Abolhassani H , Itan Y , Jin X , Hammarstrom L . VIPPID: a gene‐specific single nucleotide variant pathogenicity prediction tool for primary immunodeficiency diseases. Brief Bioinform. 2022;23 : PMID: 35598327 {Medline}.
Fowler DM , Adams DJ , Gloyn AL , Hahn WC , Marks DS , Muffley LA , et al. An atlas of variant effects to understand the genome at nucleotide resolution. Genome Biol. 2023;24 :147 PMID: 37394429 {Medline}.37394429
Frieden C , Wang H , Ho CMW . A mechanism for lipid binding to apoE and the role of intrinsically disordered regions coupled to domain‐domain interactions. Proc Natl Acad Sci U S A. 2017;114 :6292–6297. PMID: 28559318 {Medline}.28559318
Ganesan K , Kulandaisamy A , Binny Priya S , Gromiha MM . HuVarBase: a human variant database with comprehensive information at gene and protein levels. PLoS One. 2019;14 :e0210475. PMID: 30703169 {Medline}.30703169
Grant MA , Lazo ND , Lomakin A , Condron MM , Arai H , Yamin G , et al. Familial Alzheimer's disease mutations alter the stability of the amyloid beta‐protein monomer folding nucleus. Proc Natl Acad Sci U S A. 2007;104 :16522–16527. PMID: 17940047 {Medline}.17940047
Gratuze M , Leyns CEG , Holtzman DM . New insights into the role of TREM2 in Alzheimer's disease. Mol Neurodegener. 2018;13 :66 PMID: 30572908 {Medline}.30572908
Gray VE , Hause RJ , Luebeck J , Shendure J , Fowler DM . Quantitative missense variant effect prediction using large‐scale mutagenesis data. Cell Syst. 2018;6 :116–124 e113. PMID: 29226803 {Medline}.29226803
Hall‐Roberts H , Agarwal D , Obst J , Smith TB , Monzon‐Sandoval J , di Daniel E , et al. TREM2 Alzheimer's variant R47H causes similar transcriptional dysregulation to knockout, yet only subtle functional phenotypes in human iPSC‐derived macrophages. Alzheimers Res Ther. 2020;12 :151 PMID: 33198789 {Medline}.33198789
Hecht M , Bromberg Y , Rost B . Better prediction of functional effects for sequence variants. BMC Genomics. 2015;16 (Suppl 8 ):S1 PMID: 26110438 {Medline}.
Jagota M , Ye C , Albors C , Rastogi R , Koehl A , Ioannidis N , et al. Cross‐protein transfer learning substantially improves disease variant prediction. Genome Biol. 2023;24 :182 PMID: 37550700 {Medline}.37550700
Jin SC , Benitez BA , Karch CM , Cooper B , Skorupa T , Carrell D , et al. Coding variants in TREM2 increase risk for Alzheimer's disease. Hum Mol Genet. 2014;23 :5838–5846. PMID: 24899047 {Medline}.24899047
Joosten RP , te Beek TA , Krieger E , Hekkelman ML , Hooft RW , Schneider R , et al. A series of PDB related databases for everyday needs. Nucleic Acids Res. 2011;39 :D411–D419. PMID: 21071423 {Medline}.21071423
Jubb HC , Higueruelo AP , Ochoa‐Montano B , Pitt WR , Ascher DB , Blundell TL . Arpeggio: a web server for calculating and visualising interatomic interactions in protein structures. J Mol Biol. 2017;429 :365–371. PMID: 27964945 {Medline}.27964945
Jumper J , Evans R , Pritzel A , Green T , Figurnov M , Ronneberger O , et al. Highly accurate protein structure prediction with AlphaFold. Nature. 2021;596 :583–589. PMID: 34265844 {Medline}.34265844
Kabsch W , Sander C . Dictionary of protein secondary structure: pattern recognition of hydrogen‐bonded and geometrical features. Biopolymers. 1983;22 :2577–2637. PMID: 6667333 {Medline}.6667333
Kawashima S , Ogata H , Kanehisa M . AAindex: amino acid index database. Nucleic Acids Res. 1999;27 :368–369. PMID: 9847231 {Medline}.9847231
Kelleher RJ 3rd , Shen J . Presenilin‐1 mutations and Alzheimer's disease. Proc Natl Acad Sci U S A. 2017;114 :629–631. PMID: 28082723 {Medline}.28082723
Kim W , Hecht MH . Generic hydrophobic residues are sufficient to promote aggregation of the Alzheimer's Abeta42 peptide. Proc Natl Acad Sci U S A. 2006;103 :15824–15829. PMID: 17038501 {Medline}, 19164, 19169.17038501
Kinoshita J , Clark T . Alzforum. Methods Mol Biol. 2007;401 :365–381. PMID: 18368375 {Medline}.18368375
Klein P , Kanehisa M , DeLisi C . Prediction of protein function from sequence properties. Discriminant analysis of a data base. Biochim Biophys Acta. 1984;787 :221–226. PMID: 6547351 {Medline}.6547351
Knopman DS , Amieva H , Petersen RC , Chetelat G , Holtzman DM , Hyman BT , et al. Alzheimer disease. Nat Rev Dis Primers. 2021;7 :33 PMID: 33986301 {Medline}.33986301
Kulandaisamy A , Parvathy Dharshini SA , Gromiha MM . Alz‐disc: a tool to discriminate disease‐causing and neutral mutations in Alzheimer's disease. Comb Chem High Throughput Screen. 2023;26 :769–777. PMID: 35619290 {Medline}.35619290
Landrum MJ , Lee JM , Benson M , Brown GR , Chao C , Chitipiralla S , et al. ClinVar: improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 2018;46 :D1062–D1067. PMID: 29165669 {Medline}.29165669
Landrum MJ , Lee JM , Riley GR , Jang W , Rubinstein WS , Church DM , et al. ClinVar: public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res. 2014;42 :D980–D985. PMID: 24234437 {Medline}.24234437
Li G , Panday SK , Alexov E . SAAFEC‐SEQ: a sequence‐based method for predicting the effect of single point mutations on protein thermodynamic stability. Int J Mol Sci. 2021;22 : PMID: 33435356 {Medline}.
Mariani V , Biasini M , Barbato A , Schwede T . lDDT: a local superposition‐free score for comparing protein structures and models using distance difference tests. Bioinformatics. 2013;29 :2722–2728. PMID: 23986568 {Medline}.23986568
Marti‐Renom MA , Stuart AC , Fiser A , Sanchez R , Melo F , Sali A . Comparative protein structure modeling of genes and genomes. Annu Rev Biophys Biomol Struct. 2000;29 :291–325. PMID: 10940251 {Medline}.10940251
Meirovitch H , Rackovsky S , Scheraga HA . Empirical studies of hydrophobicity. 1. Effect of protein size on the hydrophobic behavior of amino acids. Macromolecules. 1980;13 :1398–1405.
Munro D , Singh M . DeMaSk: a deep mutational scanning substitution matrix and its use for variant impact prediction. Bioinformatics. 2021;36 :5322–5329. PMID: 33325500 {Medline}.33325500
Myung Y , Rodrigues CHM , Ascher DB , Pires DEV . mCSM‐AB2: guiding rational antibody design using graph‐based signatures. Bioinformatics. 2020;36 :1453–1459. PMID: 31665262 {Medline}.31665262
Ng PC , Henikoff S . SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31 :3812–3814. PMID: 12824425 {Medline}.12824425
Notin P , Kollasch AW , Ritter D , van Niekerk L , Paul S , Spinner H , et al. ProteinGym: large‐scale benchmarks for protein design and fitness prediction. bioRxiv. 2023; PMID: 38106144 {Medline}.
Pan Q , Nguyen TB , Ascher DB , Pires DEV . Systematic evaluation of computational tools to predict the effects of mutations on protein stability in the absence of experimental structures. Brief Bioinform. 2022;23 : PMID: 35189634 {Medline}.
Pan Q , Portelli S , Nguyen TB , Ascher DB . Characterization on the oncogenic effect of the missense mutations of p53 via machine learning. Brief Bioinform. 2023;25 : PMID: 38018912 {Medline}.
Panegyres PK , Chen HY . Differences between early and late onset Alzheimer's disease. Am J Neurodegener Dis. 2013;2 :300–306. PMID: 24319647 {Medline}.24319647
Pires DE , Ascher DB , Blundell TL . mCSM: predicting the effects of mutations in proteins using graph‐based signatures. Bioinformatics. 2014;30 :335–342. PMID: 24281696 {Medline}.24281696
Portelli S , Albanaz A , Pires DEV , Ascher DB . Identifying the molecular drivers of ALS‐implicated missense mutations. J Med Genet. 2023;60 :484–490. PMID: 36180205 {Medline}.36180205
Portelli S , Barr L , de Sa AGC , Pires DEV , Ascher DB . Distinguishing between PTEN clinical phenotypes through mutation analysis. Comput Struct Biotechnol J. 2021;19 :3097–3109. PMID: 34141133 {Medline}.34141133
Rezaei‐Ghaleh N , Amininasab M , Giller K , Becker S . Familial Alzheimer's disease‐related mutations differentially alter stability of amyloid‐beta aggregates. J Phys Chem Lett. 2023;14 :1427–1435. PMID: 36734539 {Medline}.36734539
Richards S , Aziz N , Bale S , Bick D , Das S , Gastier‐Foster J , et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17 :405–424. PMID: 25741868 {Medline}.25741868
Rodrigues CHM , Pires DEV , Ascher DB . DynaMut2: assessing changes in stability and flexibility upon single and multiple point missense mutations. Protein Sci. 2021;30 :60–69. PMID: 32881105 {Medline}.32881105
Roher AE , Lowenson JD , Clarke S , Wolkow C , Wang R , Cotter RJ , et al. Structural alterations in the peptide backbone of beta‐amyloid core protein may account for its deposition and stability in Alzheimer's disease. J Biol Chem. 1993;268 :3072–3083. PMID: 8428986 {Medline}.8428986
Rose AS , Bradley AR , Valasatava Y , Duarte JM , Prlic A , Rose PW . NGL viewer: web‐based molecular graphics for large complexes. Bioinformatics. 2018;34 :3755–3758. PMID: 29850778 {Medline}.29850778
Sali A , Blundell TL . Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993;234 :779–815. PMID: 8254673 {Medline}.8254673
Sayed FA , Kodama L , Fan L , Carling GK , Udeochu JC , Le D , et al. AD‐linked R47H‐TREM2 mutation induces disease‐enhancing microglial states via AKT hyperactivation. Sci Transl Med. 2021;13 :eabe3947 PMID: 34851693 {Medline}.34851693
Sayers EW , Bolton EE , Brister JR , Canese K , Chan J , Comeau DC , et al. Database resources of the national center for biotechnology information. Nucleic Acids Res. 2022;50 :D20–D26. PMID: 34850941 {Medline}.34850941
Scheltens P , de Strooper B , Kivipelto M , Holstege H , Chetelat G , Teunissen CE , et al. Alzheimer's disease. Lancet. 2021;397 :1577–1590. PMID: 33667416 {Medline}.33667416
Sgourakis NG , Yan Y , McCallum SA , Wang C , Garcia AE . The Alzheimer's peptides Abeta40 and 42 adopt distinct conformations in water: a combined MD/NMR study. J Mol Biol. 2007;368 :1448–1457. PMID: 17397862 {Medline}.17397862
Sherrington R , Froelich S , Sorbi S , Campion D , Chi H , Rogaeva EA , et al. Alzheimer's disease associated with mutations in presenilin 2 is rare and variably penetrant. Hum Mol Genet. 1996;5 :985–988. PMID: 8817335 {Medline}.8817335
Sirkis DW , Bonham LW , Aparicio RE , Geier EG , Ramos EM , Wang Q , et al. Rare TREM2 variants associated with Alzheimer's disease display reduced cell surface expression. Acta Neuropathol Commun. 2016;4 :98 PMID: 27589997 {Medline}.27589997
Smigielski EM , Sirotkin K , Ward M , Sherry ST . dbSNP: a database of single nucleotide polymorphisms. Nucleic Acids Res. 2000;28 :352–355. PMID: 10592272 {Medline}.10592272
Somavarapu AK , Kepp KP . Loss of stability and hydrophobicity of presenilin 1 mutations causing Alzheimer's disease. J Neurochem. 2016;137 :101–111. PMID: 26756738 {Medline}.26756738
Lohr T , Kohlhoff K, Heller GT , Camilloni C , Vendruscolo M . A small molecule stabilises the disordered native state of the Alzheimer's Aβ peptide. ACS Chem Neurosci. 2022;13 :1738–1745. PMID: 35649268 {Medline}.35649268
van Gils JHM , van Dijk E , Peduzzo A , Hofmann A , Vettore N , Schutzmann MP , et al. The hydrophobic effect characterises the thermodynamic signature of amyloid fibril growth. PLoS Comput Biol. 2020;16 :e1007767 PMID: 32365068 {Medline}.32365068
Vasquez M , Nemethy G , Scheraga HA . Computed conformational states of the 20 naturally occurring amino acid residues and of the prototype residue α‐aminobutyric acid. Macromolecules. 1983;16 :1043–1049.
Wakutani Y , Watanabe K , Adachi Y , Wada‐Isoe K , Urakami K , Ninomiya H , et al. Novel amyloid precursor protein gene missense mutation (D678N) in probable familial Alzheimer's disease. J Neurol Neurosurg Psychiatry. 2004;75 :1039–1042. PMID: 15201367 {Medline}.15201367
Wang Y , Cella M , Mallinson K , Ulrich JD , Young KL , Robinette ML , et al. TREM2 lipid sensing sustains the microglial response in an Alzheimer's disease model. Cell. 2015;160 :1061–1071. PMID: 25728668 {Medline}.25728668
Wattmo C , Minthon L , Wallin AK . Mild versus moderate stages of Alzheimer's disease: three‐year outcomes in a routine clinical setting of cholinesterase inhibitor therapy. Alzheimers Res Ther. 2016;8 :7 PMID: 26883213 {Medline}.26883213
Webb B , Sali A . Comparative protein structure modeling using MODELLER. Curr Protoc Bioinformatics. 2016;54 :5–6. PMID: 27322406 {Medline}.
Wu Y , Li R , Sun S , Weile J , Roth FP . Improved pathogenicity prediction for rare human missense variants. Am J Hum Genet. 2021;108 :1891–1906. PMID: 34551312 {Medline}, 2389.34551312
Zamyatnin AA . Protein volume in solution. Prog Biophys Mol Biol. 1972;24 :107–123. PMID: 4566650 {Medline}.4566650
Zhang W , Jiao B , Xiao T , Liu X , Liao X , Xiao X , et al. Association of rare variants in neurodegenerative genes with familial Alzheimer's disease. Ann Clin Transl Neurol. 2020;7 :1985–1995. PMID: 32941707 {Medline}.32941707
Zhou Y , Pan Q , Pires DEV , Rodrigues CHM , Ascher DB . DDMut: predicting effects of mutations on protein stability using deep learning. Nucleic Acids Res PMID: 37283042 {Medline}. 2023;51 :W122–W128.37283042
