Quantum mechanics study of the hydroxyethylamines-BACE-1 active site interaction energies

The identification of BACE-1, a key enzyme in the production of Amyloid-β (Aβ) peptides, generated by the proteolytic processing of amyloid precursor protein, was a major advance in the field of Alzheimer's disease as this pathology is characterized by the presence of extracellular senile plaques, mainly comprised of Aβ peptides. Hydroxyethylamines have demonstrated a remarkable potential, like candidate drugs, for this disease using BACE-1 as target. Density Functional Theory calculations were employed to estimate interaction energies for the complexes formed between the hydroxyethylamine derivated inhibitors and 24 residues in the BACE-1 active site. The collected data offered not only a general but a particular quantitative description that gives a deep insight of the interactions in the active site, showing at the same time how ligand structural variations affect them. Polar interactions are the major energetic contributors for complex stabilization and those ones with charged aspartate residues are highlighted, as they contribute over 90% of the total attractive interaction energy. Ligand-ARG296 residue interaction reports the most repulsive value and decreasing of the magnitude of this repulsion can be a key feature for the design of novel and more potent BACE-1 inhibitors. Also it was explained why sultam derivated BACE-1 inhibitors are better ones than lactam based. Hydrophobic interactions concentrated at S1 zone and other relevant repulsions and attractions were also evaluated. The comparison of two different theory levels (X3LYP and M062X) allowed to confirm the relevance of the detected interactions as each theory level has its own strength to depict the forces involved, as is the case of M062X which is better describing the hydrophobic interactions. Those facts were also evaluated and confirmed by comparing the quantitative trend, of selected ligand-residue interactions, with MP2 theory level as reference standard method for electrostatic plus dispersion energies.