A Comprehensive In Silico Analysis of the Functional and Structural Impact of Nonsynonymous SNPs in the ABCA1 Transporter Gene

Abstract

Disease phenotypes and defects in function can be traced to nonsynonymous single nucleotide polymorphisms (nsSNPs), which are important indicators of action sites and effective potential therapeutic approaches. Identification of deleterious nsSNPs is crucial to characterize the genetic basis of diseases, assess individual susceptibility to disease, determinate molecular and therapeutic targets, and predict clinical phenotypes. In this study using PolyPhen2 and MutPred in silico algorithms, we analyzed the genetic variations that can alter the expression and function of the ABCA1 gene that causes the allelic disorders familial hypoalphalipoproteinemia and Tangier disease. Predictions were validated with published results from in vitro, in vivo, and human studies. Out of a total of 233 nsSNPs, 80 (34.33%) were found deleterious by both methods. Among these 80 deleterious nsSNPs found, 29 (12.44%) rare variants resulted highly deleterious with a probability >0.8. We have observed that mostly variants with verified functional effect in experimental studies are correctly predicted as damage variants by MutPred and PolyPhen2 tools. Still, the controversial results of experimental approaches correspond to nsSNPs predicted as neutral by both methods, or contradictory predictions are obtained for them. A total of seventeen nsSNPs were predicted as deleterious by PolyPhen2, which resulted neutral by MutPred. Otherwise, forty two nsSNPs were predicted as deleterious by MutPred, which resulted neutral by PolyPhen2.

Figure 1

Predicted structure of ABCA1 [21, 23, 105]. The protein consists of 2261 amino acids and comprises 2 halves of similar structure. Each half encodes a transmembrane domain containing 6 helices (1–6 and 7–12) and 1 nucleotide binding domain (NBD-1 and NBD-2), where ATP is bound and utilized as energy for substrate transport across the membrane. Each NBD contains 2 highly conserved peptide motifs known as Walker A and Walker B, which are present in many proteins that use ATP and a Walker C signature unique to ABC transporters. ABCA1 is predicted to have an N terminus oriented into the cytosol and two large extracellular loops that are highly glycosylated and linked by cysteine bonds (Y indicates the glycosylation sites, and S–S indicates predicted disulfide bonds).

Figure 2

Workflow for predicting the effect of nsSNPs on ABCA1 protein. Gene specific variations were retrieved from the main variation databases. Mutations were annotated and nsSNPs belonging to the canonical transcript were selected in order to run the algorithms for functional prediction.

Figure 3

ROC (receiver operating characteristics) curve comparing the performance of MutPred and PolyPhen2 methods in predicting the outcome of nsSNPs. The dashed curve corresponds to MutPred predictions. As the area under the curve (AUC) for the MutPred method predictions is larger than the areas over the curve (AOC) which corresponds to the PolyPhen2 predictions curve, we can confirm that MutPred outperforms PolyPhen2 in predicting the outcome of nsSNPs. The dataset used for the performance comparison was obtained from Varibench as stated in Methods section.

Figure 4

Distribution of ABCA1 nonsynonymous SNPs (nsSNPs), synonymous SNPs (sSNPs), 3′-UTR and 5′-UTR, and intronic SNPs.

Figure 5

Scatter plot shows the correlation between the predictions made by PolyPhen2 (x-axis) and MutPred (y-axis) for 233 amino acid substitutions on ABCA1. The diagonal line represents a perfect correlation (=1) between both prediction methods for every mutation. The overall correlation of the predictions made by both methods is high (0.57). The majority of mutations classified as pathogenic by PolyPhen2 with the highest score (=1) are also classified as pathogenic by MutPred but within a score range between 0.51 and 1.