Population density analysis for determining the protonation state of the catalytic dyad in BACE1-tertiary carbinamine-based inhibitor complex

Computational modelling of potent β-secretase (BACE1) inhibitors towards Alzheimer's disease treatment

Researchers have identified the β-amyloid precursor protein cleaving enzyme 1 (BACE1) in the multifactorial pathway of Alzheimer's disease (AD) as a drug target. The design and development of molecules to inhibit BACE1 as a potential cure for AD thus remained significant. Herein, we simulated two potent BACE1 inhibitors (AM-6494 and CNP-520) to understand their binding affinity at the atomistic level. AM-6494 is a newly reported potent BACE1 inhibitor with an IC₅₀ value of 0.4 nM in vivo and now picked for preclinical considerations. Umibecestat (CNP-520), which was discontinued at human trials lately, was considered to enable a reasonable evaluation of our results. Using density functional theory (DFT) and Our Own N-layered Integrated molecular Orbital and Molecular Mechanics (ONIOM), we achieved the aim of this investigation. These computational approaches enabled the prediction of the electronic properties of AM-6494 and CNP-520 plus their binding energies when complexed with BACE1. For AM-6494 and CNP-520 interaction with protonated BACE1, the ONIOM calculation gave binding free energy of -62.849 and -33.463 kcal/mol, respectively. In the unprotonated model, we observed binding free energy of -59.758 kcal/mol in AM-6494. Taken together thermochemistry of the process and molecular interaction plot, AM-6494 is more favourable than CNP-520 towards the inhibition of BACE1. The protonated model gave slightly better binding energy than the unprotonated form. However, both models could sufficiently describe ligand binding to BACE1 at the atomistic level. Understanding the detailed molecular interaction of these inhibitors could serve as a basis for pharmacophore exploration towards improved inhibitor design.

G-score: A function to solve the puzzle of modeling the protonation states of β-secretase binding pocket

The population density concept has emerged as a proposal for the analysis of molecular dynamics results, the key characteristic of population density is the evaluation of the simultaneous occurrence of a set of relevant parameters for a system. However, despite its statistical strength, selection of the tolerance level for the comparison of different models may appear as arbitrary. This work introduces the G-score, a function which summarizes and categorizes the results of population density analysis. Additionally, it incorporates parameters based on rmsd and dihedral angles, besides the protein-protein and protein-ligand interatomic distances conventionally used, which complement each other to provide a better description of the behavior of the system. These newly-proposed tools were applied to determine the most probable protonation state of the aspartic dyad of BACE1, Asp93 and Asp289, in the presence of three types of transition state inhibitors namely: reduced amides, tertiary carbinamines and hydroxyethylamines. The results show a full agreement between G-score values and population density charts, with the advantage of allowing a quick and direct comparison among all the considered models. We anticipate that the simplicity of calculating the parameters employed in this study will permit the extensive use of population density and the G-score for other molecular systems.