Determining Structural Changes for Ligand Recognition between Human and Rat Phosphorylated BACE1 in Silico and Its Phosphorylation by GSK3β at Thr252 by in Vitro Studies

Alzheimer's disease (AD) is a neurodegenerative disease affecting older adults. AD pathogenesis involves the production of the highly neurotoxic amyloid-β peptide 1-42 (Aβ1-42) from β-site amyloid precursor protein cleaving enzyme 1 (BACE1). The phosphorylation of BACE1 at Thr252 increases its enzymatic activity. This study examined the phosphorylation of BACE1 from human and rat BACE1 in silico through phosphorylation predictors. Besides, we explored how phosphorylation at various sites affected the BACE1 structure and its affinity with amyloid precursor protein (APP) and six BACE1 inhibitors. Additionally, we evaluated the phosphorylation of Thr252-BACE1 by glycogen synthase kinase 3 β (GSK3β) in vitro. The phosphorylation predictors showed that Thr252, Ser59, Tyr76, Ser71, and Ser83 could be phosphorylated. Also, Ser127 in rat BACE1 can be phosphorylated, but human BACE1 has a Gly at this position. Molecular dynamics simulations showed that Ser127 plays an important role in the open and closed BACE1 conformational structures. Docking studies and the molecular mechanics generalized Born surface area (MMGBSA) approach showed that human BACE1 phosphorylated at Thr252 and rat BACE1 phosphorylated at Ser71 have the best binding and free energy with APP, forming hydrogen bonds with Asp672. Importantly, inhibitors have a higher affinity for the phosphorylated rat BACE1 than for its human counterpart, which could explain their failure during clinical trials. Finally, in vitro experiments showed that GSK3β could phosphorylate BACE1. In conclusion, BACE1 phosphorylation influences the BACE1 conformation and its recognition of ligands and substrates. Thus, these features should be carefully considered in the design of BACE1 inhibitors.

Molecular dynamics simulations depict structural motions of the whole human aryl hydrocarbon receptor influencing its binding of ligands and HSP90

The aryl hydrocarbon receptor (AhR) has broad biological functions when its ligands activate it; the non-binding interactions with AhR have not been fully elucidated due to the absence of a complete tridimensional (3D) structure. Therefore, utilization of the whole 3D structure from Homo sapiens AhR by in silico studies will allow us to better study and analyze the binding mode of its full and partial agonists, and antagonists, as well as its interaction with the HSP90 chaperone. The 3D AhR structure was obtained from I-TASSER and subjected to molecular dynamics (MD) simulations to obtain different structural conformations and determine the most populated AhR conformer by clustering analyses. The AhR-3D structures selected from MD simulations and those from clustering analyses were used to achieve docking studies with some of its ligands and protein-protein docking with HSP90. Once the AhR-3D structure was built, its Ramachandran maps and energy showed a well-qualified 3D model. MD simulations showed that the per-Arnt-Sim homology (PAS) PAS A, PAS B, and Q domains underwent conformational changes, identifying the conformation when agonists were binding also, and HSP90 was binding near the PAS A, PAS B, and Q domains. However, when antagonists are binding, HSP90 does not bind near the PAS A, PAS B, and Q domains. These studies show that the complex agonist-AhR-HSP90 can be formed, but this complex is not formed when an antagonist is binding. Knowing the conformations when the ligands bind to AHR and the behavior of HSP90 allows for an understanding of its activity.Communicated by Ramaswamy H. Sarma.