Rare Variant Analysis and Molecular Dynamics Simulation in Alzheimer’s Disease Identifies Exonic Variants in FLG

Background: Although an increasing number of common variants contributing to Alzheimer’s disease (AD) are uncovered by genome-wide association studies, they can only explain less than half of the heritability of AD. Rare variant association studies (RVAS) has become an increasingly important area to explain the risk or trait variability of AD. Method: To investigate the potential rare variants that cause AD, we screened 70,209 rare variants from two cohorts of a 175 AD cohort and a 214 cognitively normal cohort from the Alzheimer’s Disease Neuroimaging Initiative database. MIRARE, a novel RVAS method, was performed on 232 non-synonymous variants selected by ANNOVAR annotation. Molecular docking and molecular dynamics (MD) simulation were adopted to verify the interaction between the chosen functional variants and BACE1. Results: MIRAGE analysis revealed significant associations between AD and six potential pathogenic genes, including PREX2, FLG, DHX16, NID2, ZnF585B and ZnF875. Only interactions between FLG (including wild type and rs3120654(SER742TYR)) and BACE1 were verified by molecular docking and MD simulation. The interaction of FLG(SER742TYR) with BACE1 was greater than that of wildtype FLG with BACE1. Conclusions: According to the literature search, bio-informatics analysis, and molecular docking and MD simulation, we find non-synonymous rare variants in six genes, especially FLG(rs3120654), that may play key roles in AD.