RESEARCH ARTICLE


Interpreting the Mechanism of APOE (p.Leu167del) Mutation in the Incidence of Familial Hypercholesterolemia; An In-silico Approach



Omran Mohammed Rashidi1, #, Fatima Amanullah H.Nazar2, #, Mohamed Nabil Alama3, Zuhier Ahmed Awan1, *
1 Department of Clinical Biochemistry. Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
2 Department of Biology, Genomic and Biotechnology Section. Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia
3 Adult interventional cardiology, Cardiology unit, King Abdulaziz University Hospital (KAUH), Jeddah, Saudi Arabia


Article Metrics

CrossRef Citations:
0
Total Statistics:

Full-Text HTML Views: 1483
Abstract HTML Views: 492
PDF Downloads: 706
ePub Downloads: 537
Total Views/Downloads: 3218
Unique Statistics:

Full-Text HTML Views: 735
Abstract HTML Views: 274
PDF Downloads: 279
ePub Downloads: 165
Total Views/Downloads: 1453



© 2017 Rashidi et al.

open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC-BY 4.0), a copy of which is available at: https://creativecommons.org/licenses/by/4.0/legalcode. This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

* Address correspondence to this author at the Division of Clinical Biochemistry, King Abdulaziz University, Abdullah Sulayman, Jeddah 22254, Saudi Arabia, Tel: +966-12-640-1000# 73717 / 22123; Fax: 640-8451; E-mails: zawan@kau.edu.sa; zuhier.awan@gmail.com

# Equal contribution and Shared first authorship


Abstract

Background:

Apolipoprotein E (APOE) gene is a ligand protein in humans which mediates the metabolism of cholesterol by binding to the low-density lipoprotein receptor (LDLR). P.Leu167del mutation in APOE gene was recently connected with Familial Hypercholesterolemia, a condition associated with premature cardiovascular disease. The consequences of this mutation on the protein structure and its receptor binding capacity remain largely unknown.

Objective:

The current study aims to further decipher the underlying mechanism of this mutation using advanced software-based algorithms. The consequences of disrupting the leucine zipper by this mutation was studied at the structural and functional level of the APOE protein.

Methods:

3D protein modeling for both APOE and LDLR (wild types), along with APOE (p.Leu167del) mutant type were generated using homology modeling template-based alignment. Structural deviation analysis was performed to evaluate the spatial orientation and the stability of the mutant APOE structure. Molecular docking analysis simulating APOE-LDLR protein interaction was carried out, in order to evaluate the impact of the mutation on the binding affinity.

Result:

Structural deviation analysis for APOE mutated model showed low degree of deviance scoring root-mean-square deviation, (RMSD) = 0.322 Å. Whereas Docking simulation revealed an enhanced molecular interaction towards the LDLR with an estimation of +171.03 kJ/mol difference in binding free energy.

Conclusion:

This in-silico study suggests that p.Leu167del is causing the protein APOE to associate strongly with its receptor, LDLR. This gain-of-function is likely hindering the ability of LDLR to be effectively recycled back to the surface of the hepatocytes to clear cholesterol from the circulation therefore leading to FH.



Keywords: Apolipoprotein E (APOE), P. Leu167del, Familial Hypercholesterolemia (FH), Homology modeling, Structural deviation, Molecular Docking Simulation.



1. INTRODUCTION

Familial Hypercholesterolemia (FH) is a genetic lipoprotein disorder characterized by elevated serum low-density lipoprotein cholesterol (LDL-C). It is inherited either as an autosomal dominant (AD) or as an autosomal recessive (AR) hypercholesterolemia. The worldwide prevalence of the heterozygous form is 1 in 500, while the homozygous is 1 in 1,000,000 [1]. In special populations there is an evidence that it is even more common due to inbreeding (French-Canadian) and consanguinity (Lebanese) [1]. Excessive levels of LDL-C can increase the risk of premature coronary heart disease through acceleration of atherosclerosis [2]. Genes associated with AD-FH phenotype involve: the low density lipoprotein receptor (LDLR) [2, 3], apolipoprotein B (APOB) [4] and proprotein convertase subtilisin/kexin type 9 (PCSK9) [5]. The AR form of FH was found to be caused by the very rare loss-of-function mutations in the low-density lipoprotein receptor adapter protein 1 (LDLRAP1)gene [6].

Genome-wide association studies (GWAS) have reported that the Apolipoprotein E (APOE) gene as strongly associated with LDL-C levels [7]. APOE is a multifunctional plasma protein that is regularly synthesized in several organs such as the liver, brain, kidney and spleen [8]. It is an essential component that mediates binding of all lipoproteins [9, 10]. The APOE-LDLR protein complex contributes significantly to cholesterol levels in the regular cycle of lipid metabolism [11]. The APOE gene exists in 3 major allelic variants; E2, E3 and E4, where E3 is considered to be the normal form, while other alleles have been associated with various diseases.

The p.Leu167del mutation was first reported by Marduel et al. to be strongly associated with a large family of positive autosomal dominant hypercholesterolemia (ADH) occurrence [12]. Later on, Awan et al. presented how p.Leu167del structurally altered the LDLR binding domain within the APOE protein. This mutation disrupted the leucine zipper near its 4th α-helical structure, which is structurally kown to be highly critical region in the protein. The previous observation was registered via computational biology using established template-based protein modeling, followed by 3D structure alignment of both the reference and Leu167del proteins [13].

The purpose of the current study is to decipher the underlying mechanism of p.Leu167del mutation using advanced software-based algorithms. We are aiming to further elucidate the consequences of disrupting the leucine zipper at the structural and functional levels of the APOE protein.

2. MATERIALS AND METHODS

2.1. Data Set

The reference genomic data for human APOE and LDLR genes were obtained from the National Center for Biological Information (NCBI). Corresponding amino acid sequences of APOE protein (CCDS Transcript ID: ENST00000252486.8) and LDLR protein (CCDS Transcript ID: ENST00000558518.5) were retrieved from Ensemble genome browser (http://www.ensembl.org/index.html).

2.2. 3D Protein Modeling

Crystalized protein structures of both candidate genes were absent from all protein and allied genomic databases. Therefore, computational modeling based on homology prediction was utilized to construct the reference protein structures. I-Tasser web server (http://zhanglab.ccmb.med.umich.edu/I-TASSER/) was employed in the modeling procedure. Protein models were generated based on multiple sequence alignment and fragment assembly simulations. I-Tasser generated 5 protein data bank (PDB) formatted models for each sequence. Only the top model with the highest confidence score (C-score), accompanied with the most reliable measurement of template modeling score (TM-scores) were selected for the follow-up analysis.

Predicted reference models of APOE and LDLR were subsequently subjected to energy minimization using DeepView-Swiss PDB Viewer (http://spdbv.vital-it.ch/). This software depends on GROningen MOlecular Simulation (GROMOS) 96 force field techniques to release local constraints at their local residues [14]. It is also noteworthy that the p.Leu167del mutation was manually curated in the referenced protein sequence of APOE before it was submitted for automodeling using Swiss model: automated homology modeling (https://swissmodel.expasy.org/) [15].

2.3. Protein Structural Divergence Analysis

The level of structural deviance of the mutant APOE protein was homologically established using protein structural alignment in accordance with the wild-type model. Using Pymol molecular viewer(https://www.pymol.org/), sequence alignment was performed followed by a structural superposition of the mutant and wild-type PDB files. The bonds of alpha carbon (Cα) atoms of the backbone peptide are more often considered for their structural alignment. They usually experience minimal variant confirmations. Pymol proceeded in multiple cycles of refinement and structural screening, ending with the calculation of Root Mean Square Deviation (RMSD). This score can illustrate the structural drift to the polypeptide chain level when both mutant and the wild type models are structurally aligned. Structural deformities induced by genetic mutations as p.Leu167del are usually affecting RMSD scores, and therefore model structures [16].

2.4. Molecular Docking Simulation

The Hex dock server (http://hex.loria.fr/) was used as an interactive molecular program. This server is capable of calculating docking models of complex protein pairs from their PDB format [17]. Predicted APOE and LDLR protein models were both submitted to the server into 2 separate protein complex sets; wild (LDLR + wild-type APOE) and mutant (LDLR + mutated APOE). Initially, polar hydrogen atoms were turned off, while the intermolecular axis was activated in order to identify the centroids of the protein models. As Hex reads the PDB structure, it uses its all-atom centre of mass as its centroid. These centroids are used as the local coordinate origin for docking [17].

Ellipsoid axes were assigned to recognise the difference in structural parameters, especially when protein complexes are in binding mode. Hex performed a 6-dimensional search over the full rotational ranges around each protein complex. Energy calculations were set at 180°, allowing for total rotational increments for both receptor (LDLR) and ligand (APOE) around their own centroids. The default grid spacing was set at a high-resolution scoring of 0.5 Å. The steric scan (n=16) phase of the docking was calculated to proceed at (1+40)/0.75=53 intermolecular separations, in +/- steps of 0.75 Å. The highest scoring scan orientations were applied to the final search phase (n=25) and the final energies were calculated in Kilojoules units. In the light of the above, any abnormalities regarding the binding affinity of mutant APOE protein model will be distinguished.

3. RESULTS

3.1. Protein Models

I-Tasser server predicted the models of complete peptide chains of APOE wild and LDLR via homological protein structures. The criteria for selecting the best models for further analysis depends on the confidence C-scores that typically ranges from -5 to 2. In addition to the highest TM-score, which usually ranges from 0 to 1. The score >0.5 resembles a model featuring similar structural characteristics to its template-based model.

C-score predicted by I-Tasser were equal to -0.50 and 0.04 for wild APOE and LDLR models respectively. However, the estimated TM-score for wild APOE resulted in 0.65±0.13, while 0.72±0.11 for the LDLR model. Energy minimization helped to repair distorted geometries and eliminate bad contacts in targeted proteins. Such analysis is usually recommended prior to performing computational docking simulations. Energy minimized models of APOE wild and LDLR are shown in (Fig. 1).

The mutated file of APOE model was generated instantly and by automation using wild type PDB file of wild APOE.

3.2. Structural Deviation Analysis

Aligned APOE models (wild/mutant) have yielded RMSD scores. Normally, RMSD is estimated to be <2.0 Å in the deviation analysis of protein structure. Computed RMSD values represent the average distance between the Cα backbone atoms (Fig. 2). A considerable degree of structural deviance with RMSD equalling to 0.322 Å was obtained from this analysis.

Fig. (1). (A). Protein 3D model of wild type APOE with emphasis on the Leu167aa position (B). Protein 3D model of the APOE p.Leu167del with emphasis on Leu167aa deletion from the model (C). Protein 3D model of the LDLR.

Fig. (2). APOE wild and mutated form are in a superimposed position.

3.3. Docking Analysis

Docking the wild type APOE model revealed an enhanced binding activity toward its receptor; LDLR, with an estimated binding energy of -643.33 kJ/mol. However, docking simulation of the mutated APOE model to the LDLR revealed an estimation binding energy of -472.30 kJ/mol. Interestingly, the docking simulation has revealed that the mutated APOE model is showing an increased molecular interaction towards the LDLR receptor with an estimation of +171.03 kJ/mol difference in binding free energy. Hex binding calculation results for both protein complexes are shown in (Table 1).

Table 1. Calculated binding energy of both wild-type and mutated protein complexes as predicted by the hex doc server.
Interacting Protein Complex Total Docking Energy (kJ/mol) Differences in Docking Energy Level of Binding Affinity
Receptor Protein Ligand Protein
LDLR APOE wild type -643.33 - -
APOE p.Leu167del -472.3 171.03 ↑↑

4. DISCUSSION

The field of structural analysis has significantly evolved in the quest to propose a systematic and cost-effective answer to some ambiguous issues related to biological systems. The consequence of some genetic mutations on these systems is one of the main issues that have been addressed [18]. Nowadays, the computational pipeline is widely supported with software-based algorithms. Protein-protein interaction simulation programs can confidently predict functional abrogation.

Regarding FH, computational biology has proven to be effective to illustrate a plausible variation in the genome which might indicate some pathogenicity to cause the disease phenotype [19, 20] [21, 22]. APOE protein, in particular, has been the subject of structural analysis in multiple studies [23, 24]. Awhile after the GWAS report of a locus in APOE gene to be highly associated with ADH, several separate studies investigated cases and families with the similar condition [12, 13, 25]. However, p.Leu167del mutation did not undergo extensive computational interpretation. The impact of this mutation on the structure and function of the protein can go through further investigation using advanced in-silico software specifically trained to test the effects of genetic mutations at the protein level. The current study provides additional insights on p.Leu167del mutation which has been reported to cause ADH and FH. Following our previous work [13], we adopted herein an approach, alternative to the traditional in-vivo and in-vitro studies, to elucidate the consequences of p.Leu167del mutation in APOE protein attributes at the molecular and cellular levels.

Although in-frame genetic mutations (small indels or insertions and deletions) do not alter the sequence of amino acid in the protein chain, Human Gene Mutation Database (HGMD Professional release 2015.4) reported >4000 disease-causing in-frame indels, corresponding to 2.2% of all mutations. The 1000 Genomes Project Consortium reported that 1.4% of detected exonic variants were indels [26]. Notably, functional and population annotations for these in-frame indels are becoming increasingly available [27]. Our current research also revealed that in-frame p.Leu167del has a significant, possibly pathogenic impact on its protein.

Fig. (3). An illustration of the wild-type protein complex of APOE-LDLR docking simulation.

In a former study [13], we have initially assessed the impact of this deletion on the 3D structure of APOE using online protein function prediction tools (SIFT and PolyPhen). The study hypothesized that the changed leucine zippers motif might weaken the lipoprotein particle, and prevent it from binding properly with its receptor, LDLR [13]. Leucine zippers are common helical motifs with periodic segments containing leucine residues in protein compounds [28]. Mainly they are featured in eukaryotic transcription factors as they contain highly conservative regions for DNA-protein binding activity [29]. However, they are also being annotated as one of the critical protein modules for protein-protein interactions [30]. Many tools such as SIFT and PolyPhen can predict the likely pathogenicity of certain mutations. However, the low specificity of these programmes undermines the level of their accuracy prediction [31].

Structural deviation analysis is often utilised to properly evaluate the spatial orientation of the mutant proteins, where certain degrees of residual deviance can compromise the stability and the dynamical state of protein structures, which in turn, can affect their biological function [32]. Herein, the structural deformity induced by this mutation, although disrupting the leucine zipper at the 4th α-helix, had a very minimal effect on the protein. RMSD is only a quantitative measurement of the degree of similarity of protein structures [33]. The threshold of structural deviance is usually considered as <2.0 Å [34]. Therefore, p.Leu167del does not seem to impact the APOE at a critical level.

The presence of deformed abnormalities can also affect the binding affinity of targeted protein when interacting with partner molecules [35]. Therefore, computational molecular docking has been a highly desirable method. It is specifically aimed to predict a 3-dimensional simulation of complex protein-protein interaction [36]. As shown above, this approach can simulate the actual ligand-receptor interaction to predict calculated energies that define the strength of their association [37].

The increased difference, around +171 kJ/mol in the binding energy between the mutated form of the APOE protein and LDLR was remarkable Fig. (3). While most deletion mutations altering proteins binding affinity may cause significant disturbance to the function of the protein, p.Leu167del on the other hand was interestingly found to influence APOE oppositely. The detected protein is strongly attached to its receptor protein; LDLR. Note that the LDLR protein is the only member of the LDL receptor superfamily to demonstrate high affinity for binding with the wild type APOE protein [38]. This makes our mutation of interest to be classified as a gain-of-function (GOF) mutation.

The interaction activity of the APOE-LDLR complex has been pointed as an essential cellular activity for mediation of cholesterol metabolism, as the functional interaction between both proteins had already been associated with the risk of CVD [39].

The results of the present study are consistent with a recently published report [40]. Cenarro et al. proposed that the association of p.Leu167del with ADH appears to be due to the down-regulation of the LDL receptor [40]. This can indicate that p.Leu167del resembles the characteristics of a GOF mutation in APOE protein. Similar to APOE (E4) isoform, that also negatively modulate the recycling capacity of the LDLR preventing cholesterol clearance from the circulation [41].

Eventually, the impact of p.Leu167del mutation on the binding activity of APOE might initiate a quest for investigating a possible therapeutic target. Further research is needed to better understand APOE gene normal and variant-related biological aspects leading to FH.

CONCLUSION

The mechanism of p.Leu167del in regards to hypercholesterolemia was demonstrated here by this advance computational analysis. The p.Leu167del mutation was shown to exert an increased force of attraction, which triggers the mutated APOE protein to strongly associate with the LDLR protein. This can hinder the ability of the LDLR to be effectively recycled back to the surface of hepatic cells. Subsequently, a decrease in LDL endocytosis results in an increase in LDL particles circulating in the plasma. These high cholesterol levels will increase the risk of atherosclerosis.

ETHICS APPROVAL AND CONSENT TO PARTICIPATE

Not applicable.

HUMAN AND ANIMAL RIGHTS

No Animals/Humans were used for studies that are base of this research.

CONSENT FOR PUBLICATION

Not applicable.

CONFLICT OF INTEREST

The author (editor) declares no conflict of interest, financial or otherwise.

ACKNOWLEDGEMENTS

The project was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, under grant No (G-613-114-37). ZA was also supported by Seraj Kaki Chair for the study of genetic polymorphisms in cardiovascular disease and diabetes, King Abdulaziz University, Jeddah, Saudi Arabia [grant numbers MBK/03/434]. Authors would like to thank both Dr Babajan Banaganapali and Dr Noor Shaik from Princess Al-Jawhara Center of Excellence in Research of Hereditary Disorders, Department of Genetic Medicine, King Abdul-Aziz University, Jeddah, Kingdom of Saudi Arabia for their inspirational guidance. Authors would also like to thank Mr.Almuatazbellah Awan for his contribution in proofreading the final composed version of the article.

REFRENCES

[1] Bamimore MA, Zaid A, Banerjee Y, Al-Sarraf A, Abifadel M, Seidah NG, et al. Familial hypercholesterolemia mutations in the Middle Eastern and North African region: A need for a national registry. J Clin Lipidol 2015; 9(2): 187-94.
[2] Rader DJ, Cohen J, Hobbs HH. Monogenic hypercholesterolemia: new insights in pathogenesis and treatment. J Clin Invest 2003; 111(12): 1795-803.
[3] Talmud PJ, Shah S, Whittall R, Futema M, Howard P, Cooper JA, et al. Use of low-density lipoprotein cholesterol gene score to distinguish patients with polygenic and monogenic familial hypercholesterolaemia: a case-control study. Lancet 2013; 381(9874): 1293-301.
[4] Vega GL, Grundy SM. In vivo evidence for reduced binding of low density lipoproteins to receptors as a cause of primary moderate hypercholesterolemia. J Clin Invest 1986; 78(5): 1410-4.
[5] Abifadel M, Varret M, Rabès J-P, Allard D, Ouguerram K, Devillers M, et al. Mutations in PCSK9 cause autosomal dominant hypercholesterolemia. Nat Genet 2003; 34(2): 154-6.
[6] Garcia CK, Wilund K, Arca M, Zuliani G, Fellin R, Maioli M, et al. Autosomal Recessive Hypercholesterolemia Caused by Mutations in a Putative LDL Receptor Adaptor Protein. Science (80) 2001; 292(5520): 1394-8.
[7] Asselbergs FW, Guo Y, van Iperen EP, Sivapalaratnam S, Tragante V, Lanktree MB, et al. Large-Scale Gene-Centric Meta-analysis across 32 Studies Identifies Multiple Lipid Loci. Am J Hum Genet 2012; 91(5): 823-38. Nov 2 [cited 2017 Jun 19]
[8] Mahley RW. cholesterol transport protein with expanding role in cell biology. Science 1988; 240(4852): 622-30. Apr 29 [cited 2017 Jun 19]
[9] Hauser PS, Narayanaswami V, Ryan RO. Apolipoprotein E: From lipid transport to neurobiology. Prog Lipid Res 2011; 50(1): 62-74. Jan [cited 2017 Jun 19]
[10] Jeon H, Blacklow SC. Structure and physiologic function of the low-Density lipoprotein receptor. Annu Rev Biochem 2005; 74(1): 535-62. Jun [cited 2017 Jun 19]
[11] Go G-W, Mani A. Low-density lipoprotein receptor (LDLR) family orchestrates cholesterol homeostasis. Yale J Biol Med 2012; 85(1): 19-28. Mar [cited 2017 Jun 19]
[12] Marduel M, Ouguerram K, Serre V, Bonnefont-Rousselot D, Marques-Pinheiro A, Erik Berge K, et al. Description of a Large Family with Autosomal Dominant Hypercholesterolemia Associated with the APOE p.Leu167del Mutation. Hum Mutat 2013; 34(1): 83-7. Jan [cited 2017 Jun 19]
[13] Awan Z, Choi HY, Stitziel N, Ruel I, Bamimore MA, Husa R, et al. APOE p APOE p.Leu167del mutation in familial hypercholesterolemia. Atherosclerosis 2013; 231(22): 218-2. Dec [cited 2017 Jun 19]
[14] Lindahl E, Azuara C, Koehl P, Delarue M. NOMAD-Ref: visualization, deformation and refinement of macromolecular structures based on all-atom normal mode analysis. Nucleic Acids Res 2006; 34(22): 218-2. 34(Web Server):W52–6. Jul 1 [cited 2017 Jun 19]
[15] Schwede T, Kopp J, Guex N, Peitsch MC. SWISS-MODEL: An automated protein homology-modeling server. Nucleic Acids Res 2006; 31(13): 3381-5. Jul 1 cited 2017 Jun 19
[16] Shaik NA, Kaleemuddin M, Banaganapalli B, Khan F, Shaik NS, Ajabnoor G, et al. Structural and functional characterization of pathogenic non- synonymous genetic mutations of human insulin-degrading enzyme by in silico methods. CNS Neurol Disord Drug Targets 2014; 13(3): 517-32. Apr cited 2017 Jun 19
[17] Macindoe G, Mavridis L, Venkatraman V, Devignes M-D, Ritchie DW. HexServer: an FFT-based protein docking server powered by graphics processors. Nucleic Acids Res 2010; 38(3): 517-32. Jul cited 2017 Jun 19
[18] Zhang Y, Skolnick J. The protein structure prediction problem could be solved using the current PDB library. Proc Natl Acad Sci U S A 2005; 102(4): 1029-34. Jan 25 [cited 2017 Jun 19
[19] Cuesta-López S, Falo F, Sancho J. Computational diagnosis of protein conformational diseases: Short molecular dynamics simulations reveal a fast unfolding of r-LDL mutants that cause familial hypercholesterolemia. Proteins Struct Funct Bioinforma 2006; 66(1): 87-95. Oct 16 [cited 2017 Jun 20
[20] Mollaki V, Progias P, Drogari E. Novel LDLR Variants in Patients with Familial Hypercholesterolemia: In Silico Analysis as a Tool to Predict Pathogenic Variants in Children and Their Families. Ann Hum Genet 2013; 77(5): 426-34. Sep [cited 2017 Jun 20
[21] Al-Rasadi K, Al-Waili K, Al-Zidi WA, Al-Abri AR, Al-Hinai AT, Al-Sabti HA, et al. Low-Density Lipoprotein Receptor Gene Mutation Analysis and Structure?Function Correlation in an Omani Arab Family With Familial Hypercholesterolemia. Angiology 2014; 65(10): 911-8. Nov [cited 2017 Jun 20]
[22] Angarica VE, Orozco M, Sancho J. Exploring the complete mutational space of the LDL receptor LA5 domain using molecular dynamics: linking SNPs with disease phenotypes in familial hypercholesterolemia. Hum Mol Genet 2016; 25(6): 1233-46.
[23] Frieden C, Garai K. Concerning the structure of apoE. Protein Sci 2013; 22(12): 1820-5. Dec [cited 2017 Jun 20]
[24] Hsieh Y-H, Chou C-Y. Structural and functional characterization of human apolipoprotein E 72-166 peptides in both aqueous and lipid environments. J Biomed Sci 2011; 18(1): 4.
[25] Wintjens R, Bozon D, Belabbas K. Global molecular analysis and APOE mutations in a cohort of autosomal dominant hypercholesterolemia patients in France. J Lipid Res 2011; 57(3): 482-91.
[26] Auton A, Abecasis GR, Altshuler DM, Durbin RM, Abecasis GR, Bentley DR, et al. A global reference for human genetic variation. Nature 2015; 526(7571): 68-74. Sep 30 [cited 2017 Jun 20]
[27] Douville C, Masica DL, Stenson PD, Cooper DN, Gygax DM, Kim R, et al. Assessing the Pathogenicity of Insertion and Deletion Variants with the Variant Effect Scoring Tool (VEST-Indel). Hum Mutat 2016; 37(1): 28-35. Jan [cited 2017 Jun 20]
[28] Landschulz WH, Johnson PF, McKnight SL. The leucine zipper: a hypothetical structure common to a new class of DNA binding proteins. Science 1988; 240(4860): 1759-64. Jun 24 [cited 2017 Jun 20]
[29] E ZG. Mini Review Roles of the bZIP gene family in rice. Genet Mol Res 2014; 13(2): 3025-6. cited 2017 Jun 20]
[30] Hakoshima T. Leucine Zippers. In: Encyclopedia of Life Sciences. Chichester: John Wiley & Sons, Ltd 2005.
[31] Flanagan SE, Patch A-M, Ellard S. Using SIFT and PolyPhen to Predict Loss-of-Function and Gain-of-Function Mutations. Genet Test Mol Biomarkers 2010; 14(4): 533-7. 2010 Aug [cited 2017 Jun 20]
[32] Godzik A. The structural alignment between two proteins: Is there a unique answer? Protein Sci 1996; 5(7): 1325-38. Jul [cited 2017 Jun 20]
[33] Carugo O, Pongor S. A normalized root-mean-spuare distance for comparing protein three-dimensional structures. Protein Sci 2008; 10(7): 1470-3. Dec 31 [cited 2017 Jun 20]
[34] Banaganapalli B, Mohammed K, Khan IA, Al-Aama JY, Elango R, Shaik NA. A Computational Protein Phenotype Prediction Approach to Analyze the Deleterious Mutations of Human MED12 Gene. J Cell Biochem 2016; 117(9): 2023-35. Sep [cited 2017 Jun 20]
[35] Li M, Petukh M, Alexov E, Panchenko AR. Predicting the Impact of Missense Mutations on Protein?Protein Binding Affinity. J Chem Theory Comput 2014; 10(4): 1770-80. Apr 8 [cited 2017 Jun 20]
[36] Lengauer T, Rarey M. Computational methods for biomolecular docking. Curr Opin Struct Biol 1996; 6(3): 402-6. Jun [cited 2017 Jun 20]
[37] Feig M, Onufriev A, Lee MS, Im W, Case DA, Brooks CL. Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures. J Comput Chem 2004; 25(2): 265-84.
[38] Johnson LA, Olsen RH, Merkens LS, DeBarber A, Steiner RD, Sullivan PM, et al. Apolipoprotein E?low density lipoprotein receptor interaction affects spatial memory retention and brain ApoE levels in an isoform-dependent manner. Neurobiol Dis 2014; 64: 150-62. Apr [cited 2017 Jun 20]
[39] Holtzman DM, Herz J, Bu G. Apolipoprotein E and Apolipoprotein E Receptors: Normal Biology and Roles in Alzheimer Disease. Cold Spring Harb Perspect Med 2012; 2: a006312-2. Mar 1 [cited 2017 Jun 20]
[40] Cenarro A, Etxebarria A. The p.Leu167del Mutation in APOE Gene Causes Autosomal Dominant Hypercholesterolemia by Down-regulation of LDL Receptor Expression in Hepatocytes. J Clin Endocrinol Metab 2016; 101(21): 150-62. May [cited 2017 Jun 20]
[41] Mahley RW, Weisgraber KH, Huang Y. Apolipoprotein E: structure determines function, from atherosclerosis to Alzheimer’s disease to AIDS. J Lipid Res 2008; 64(50(Supplement)): S183-8. Dec 19 [cited 2017 Jun 20]