Unveiling a novel transient druggable pocket in BACE-1 through molecular simulations: Conformational analysis and binding mode of multisite inhibitors

PLoS One. 2017 May 15;12(5):e0177683. doi: 10.1371/journal.pone.0177683. eCollection 2017.

Abstract

The critical role of BACE-1 in the formation of neurotoxic ß-amyloid peptides in the brain makes it an attractive target for an efficacious treatment of Alzheimer's disease. However, the development of clinically useful BACE-1 inhibitors has proven to be extremely challenging. In this study we examine the binding mode of a novel potent inhibitor (compound 1, with IC50 80 nM) designed by synergistic combination of two fragments-huprine and rhein-that individually are endowed with very low activity against BACE-1. Examination of crystal structures reveals no appropriate binding site large enough to accommodate 1. Therefore we have examined the conformational flexibility of BACE-1 through extended molecular dynamics simulations, paying attention to the highly flexible region shaped by loops 8-14, 154-169 and 307-318. The analysis of the protein dynamics, together with studies of pocket druggability, has allowed us to detect the transient formation of a secondary binding site, which contains Arg307 as a key residue for the interaction with small molecules, at the edge of the catalytic cleft. The formation of this druggable "floppy" pocket would enable the binding of multisite inhibitors targeting both catalytic and secondary sites. Molecular dynamics simulations of BACE-1 bound to huprine-rhein hybrid compounds support the feasibility of this hypothesis. The results provide a basis to explain the high inhibitory potency of the two enantiomeric forms of 1, together with the large dependence on the length of the oligomethylenic linker. Furthermore, the multisite hypothesis has allowed us to rationalize the inhibitory potency of a series of tacrine-chromene hybrid compounds, specifically regarding the apparent lack of sensitivity of the inhibition constant to the chemical modifications introduced in the chromene unit. Overall, these findings pave the way for the exploration of novel functionalities in the design of optimized BACE-1 multisite inhibitors.

MeSH terms

  • Amyloid Precursor Protein Secretases / antagonists & inhibitors
  • Amyloid Precursor Protein Secretases / chemistry*
  • Amyloid Precursor Protein Secretases / metabolism
  • Aspartic Acid Endopeptidases / antagonists & inhibitors
  • Aspartic Acid Endopeptidases / chemistry*
  • Aspartic Acid Endopeptidases / metabolism
  • Binding Sites*
  • Drug Discovery
  • Enzyme Inhibitors / chemistry*
  • Enzyme Inhibitors / pharmacology
  • Molecular Conformation*
  • Molecular Docking Simulation
  • Molecular Dynamics Simulation*
  • Protein Binding
  • Structure-Activity Relationship

Substances

  • Enzyme Inhibitors
  • Amyloid Precursor Protein Secretases
  • Aspartic Acid Endopeptidases
  • BACE1 protein, human

Grants and funding

This work was supported by the Ministerio de Economía y Competitividad (SAF2014-57094-R), Generalitat de Catalunya (GC; 2014SGR1189, 2014SGR52), and ICREA Academia (FJL). The Consorci de Serveis Universitaris de Catalunya (CSUC; FJL) is acknowledged for providing computational resources. This work used the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk), via the HECBioSim consortium (EPSRC Grant EP/L000253/1). A fellowship from GC to ODP is gratefully acknowledged. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.