Investigation of intermolecular interactions and stability of verubecestat in the active site of BACE1: Development of first model from QM/MM-based charge density and MD analysis

Abstract

Alzheimer disease (AD) is a cruel neurodegenerative disorder caused by the deposition of amyloid β (Aβ) peptide inside the brain. The β-secretase (beta amyloid precursor protein (APP) cleaving enzyme 1, BACE1) is one of the enzymes involved in the cleavage of APP that leads to the Aβ formation and it is the primary target for the treatment of AD. Recent report outlines that verubecestat molecule strongly inhibits BACE1; however, its structure, binding mechanism and the stability in the active site of BACE1 are not yet known. The present study aims to determine the structure, binding affinity and the stability of verubecestat molecule in the active site of BACE1 from the molecular docking, quantum mechanics/molecular mechanics (QM/MM)-based charge density analysis and molecular dynamics simulation. Verubecestat molecule was docked at BACE1; it shows high binding affinity towards BACE1. Further, the conformational geometry and the intermolecular interactions of verubecestat in the active site of BACE1 were determined. The molecule forms strong interaction with the neighboring amino acids in the active site of BACE1. The onsite QM/MM-based charge density analysis reveals the nature of charge density distribution and the topological properties of intermolecular interactions of verubecestat molecule in the active site of BACE1. The calculated electrostatic potential (ESP) of verubecestat in the active site of BACE1 displays high negative and positive ESP regions of the molecule. This onsite QM/MM analysis is more relevant to the physiological situation. The molecular dynamics simulation has been performed, which confirms the high stability and compactness of verubecestat in the active site of BACE1. The MM-generalized Born surface area and MM-Poisson Boltzmann surface area free energy calculations of verubecestat–BACE1 also confirm the high binding affinity of verubecestat.

Communicated by Ramaswamy H. Sarma

Keywords:

• β-Secretase (BACE1)–Verubecestat
• intermolecular interaction
• QM/MM-based charge density and electrostatic potential
• molecular dynamics
• ADMET prediction

Acknowledgements

K.S. would like to thank Periyar University and University Research Fellowship (URF) to carry out this research work. The authors thank Schrödinger, Bangalore, India who allowed the use of software for this study and C-DAC, Bangalore for providing the GARUDA supercomputing facility.