BACE1 inhibitory activities of peptidic BACE1 inhibitors.
Abstract
Alzheimer’s disease is the most common cause of dementia. According to the amyloid hypothesis, β-secretase (BACE1) is a promising molecular target for the development of anti-Alzheimer’s disease drugs. BACE1 triggers the formation of the amyloid-β (Aβ) peptides that are the main component of senile plaques in the brain of patients with Alzheimer’s disease. As BACE1 cleaves the amyloid precursor protein at the N-terminus of the Aβ domain, BACE1 inhibitors reduce the Aβ level in the brain. Previously, we designed a series of peptidic inhibitors that possessed a substrate transition-state analogue, and the structure-activity relationship of our inhibitors was evaluated, based on docking and scoring, using the docking simulation software Molecular Operating Environment (MOE). However, there was no association between the scoring values and the inhibitory activities at the P2 position. Hence, we hypothesized that the interaction of the P2 position of the inhibitor with the S2 site of BACE1 was critical for the mechanism of inhibition, and we proposed the novel concept of ‘electron donor bioisostere’ for drug discovery. Using this concept, we designed potent small molecule non-peptidic BACE1 inhibitors.
Keywords
- Alzheimer’s disease
- BACE1 inhibitor
- docking simulation
- electron donor bioisostere
- in-silico conformational structure-based design
1. Introduction
Alzheimer’s disease (AD), which is the most common cause of dementia, is characterized by progressive intellectual deterioration. In 1901, Alois Alzheimer, a psychiatrist and neuropathologist, observed a 51-year-old female patient at Frankfurt Asylum. The patient showed strange behavioural symptoms and the loss of short-term memory, which was later called ‘AD’. Unfortunately, the cause of AD was unclear until recently, and there have been no treatment approaches since that first report by Dr. Alzheimer over 100 years ago. Recently, the development of many drug candidates based on the amyloid hypothesis has been reported. β-secretase (BACE1; β-site amyloid precursor protein (APP) cleaving enzyme 1) is a promising molecular target for the development of anti-Alzheimer’s drugs. BACE1 triggers the formation of the amyloid-β (Aβ) peptide that is the main component of the senile plaques found in the brain of AD patients. Previously, we had designed a series of peptidic inhibitors that possessed a substrate transition-state analogue, and evaluated the structure-activity relationship of our inhibitors, based on docking and scoring, using the docking simulation software Molecular Operating Environment ((MOE), Chemical Computing Group Inc., Canada).
1.1. Pathology of Alzheimer’s disease
Although the cause of AD was unclear until recently, a breakthrough was obtained from the genetic study of some patients with familial AD. Certain mutations of the amyloid precursor protein (APP) or presenilin gene increased amyloid-β peptides (Aβs) in the brain, which indicated their involvement in the pathogenesis of AD [1–4]. Aβ is produced from APP by two processing enzymes, β-secretase and γ-secretase, which are potential molecular targets for the development of anti-AD drugs [5–7]. The cleavage sites of APP are shown in Figure 1A. BACE1, one of the processing enzymes of APP, triggers Aβ formation in the rate-limiting first step by the cleavage of APP at the Aβ domain N-terminus (β-site). BACE1 is a type-I transmembrane aspartic protease with 501 amino acids. BACE1 and APP are located in the same intracellular granules, the endoplasmic reticulum, Golgi, and trans-Golgi networks, which are acidic environments, which suggested that Aβs are produced in these locations [8]. Next, another aspartic protease, the γ-secretase complex, cleaved at the C-terminus of the Aβ domain and released Aβ peptides. The γ-secretase complex that contains a protein

Figure 1.
(A) Amyloid precursor protein (APP) and its cleavage site. (B) Early peptidic BACE1 inhibitor by Ghosh et al.
1.2. Early peptidic BACE1 inhibitors
An early inhibitor of BACE1, an aspartic protease, was designed on the basis of a substrate transition-state concept, as well as that of other aspartic proteases, such as renin and HIV protease, which have a substrate transition-state analogue at the P1 position [10–16]. It is well-known that the Swedish mutant APP (K670N and M671L double mutation, Figure 1A) is cleaved faster than wild-type APP by BACE1, which results in increased Aβ1-42 and Aβ1-40 levels. Early BACE1 inhibitors were designed based on the Swedish-mutant APP amino acid sequence. In 1999, Sinha et al. from Elan Pharmaceuticals purified the BACE1 enzyme from the human brain using a transition-state analogue based on the Swedish-mutant sequence, and succeeded in cloning the BACE1 enzyme [17]. Ghosh et al. reported the potent inhibitors
We have reported a series of peptidic BACE1 inhibitors that possessed a norstatine-type transition-state analogue [22–30]. Our early inhibitors are shown in Table 1. Octapeptide
Compound | P4 | P3 | P2 | P1 | P1’ | P2’ | P3’ | P4’ | BACE1 inhibition (%)1 |
---|---|---|---|---|---|---|---|---|---|
E | V | N2 | Nst | D | A | E | F | <20 | |
E | V | N2 | Pns | D | A | E | F | 24 | |
E | V | D | Pns | D | A | E | F | 25 | |
E | V | M | Pns | D | A | E | F | 42 | |
E | V | K3 | Pns | D | A | E | F | <20 | |
E | V | L | Pns | D | A | E | F | >90 | |
V | L | Pns | D | A | E | F | <20 | ||
L | Pns | D | A | E | F | <20 | |||
Pns | D | A | E | F | <20 | ||||
E | V | L | Pns | D | A | E | 60 | ||
E | V | L | Pns | D | A | 46 | |||
E | V | L | Pns | D | 61 | ||||
E | V | L | Pns | 34 |
Table 1.
1BACE1 inhibition activities at 2 μM.
2P2 residue corresponding to the Swedish-mutant APP sequence.
3P2 residue corresponding to the wild-type APP sequence.
2. Docking simulation and design of pentapeptidic BACE1 inhibitors
Early BACE1 inhibitors were designed using the coordinates of the first reported X-ray crystal structure (1FKN) of a complex between BACE1 and OM99-2. OM99-2 has an Asn residue, which corresponded to the P2 residue of Swedish mutant sequence. As OM99-2 in 1FKN appeared to interact with the BACE1-Arg235 side chain
Compound | Xaa | U_ele1 | U_vdw2 | U_str3 | U_dock4 | BACE1 inhibition (%)5 |
---|---|---|---|---|---|---|
−284.7 | 28.0 | 14.6 | −242.1 | — | ||
N | −184.4 | 4.2 | 10.2 | −170.0 | 28.4 | |
M | −196.1 | 14.8 | 11.4 | −169.9 | 63.9 | |
Y | −195.2 | 18.5 | 13.6 | −163.1 | 36.9 | |
D | −174.1 | 2.0 | 11.2 | −160.9 | 33.2 | |
I | −196.3 | 25.7 | 11.1 | −159.5 | 65.8 | |
F | −176.8 | 9.0 | 11.6 | −156.2 | 47.3 | |
E | −179.8 | 11.9 | 22.3 | −145.6 | 36.3 | |
L | −155.1 | 3.3 | 10.2 | −141.6 | 83.7 | |
W | −187.1 | 19.6 | 26.6 | −140.9 | 71.3 | |
Q | −148.9 | −0.3 | 9.8 | −139.4 | 14.1 | |
Cha6 | −147.1 | −1.2 | 9.7 | −138.6 | 84.1 |
Table 2.
Docking simulation of pentapeptidic BACE1 inhibitors and their score using the 1FKN X-ray crystal structure.
1Electrostatic energy between BACE1 and inhibitor (kcal/mol).
2van der Waals energy between BACE1 and inhibitor (kcal/mol).
3Strain energy of inhibitor (kcal/mol).
4Docking score (kcal/mol); U_dock = U_ele + U_vdw + U_str.
5BACE1 inhibition % at 2 μM.
6Cyclohexylalanine (Cha).

Figure 2.
(A) Calculated active sites using the Alpha Site Finder application in MOE software. (B) OM99-2 docked in BACE1. Aqua and magenta colours indicate the X-ray crystal structure 1FKN and the energy-minimized structure under the MMFF94x force field, respectively.
As the docking model using the coordinate set of 1FKN appeared to give an unfavourable score for the BACE1 inhibitor, we compared the publicly available X-ray crystal structures of BACE1-inhibitor complexes. Surprisingly, the guanidine group of BACE1-Arg235 in most crystal structures, except 1FKN, showed similar figures flopping over the P2 region of the inhibitors, and the nearest distances between the guanidino-plane of Arg235 side chain and the P2 region of the inhibitor showed similar values of approximately 3Å [31]. The P2 moieties in most of the crystal structures found to interact with the BACE1-Arg235 side chain were a methyl group, carbonyl oxygen atom, or aromatic ring, which were bound to the guanidine-plane of Arg235 side chain by CH-π, O-π, or π-π stacking interactions. This suggested that the π-orbital on the guanidino-plane interacted with the P2 region by a weak quantum force such as stacking or σ-π interaction. The only exception was the interaction in the first reported X-ray crystal structure, 1FKN. The P2 moiety of OM99-2 in the crystal structure 1FKN appeared to interact with the BACE1-Arg235 side chain

Figure 3.
Interaction of BACE1-Arg235 with BACE1 inhibitors in X-ray crystal structures. (A) PDB ID: 1FKN, (B) PDB ID: 1M4H, (C) PDB ID: 2B8L and (D) PDB ID: 2P83.

Figure 4.
Location of BACE1-Arg235 in the respective crystal structures. The blue, magenta, green and red stick models represent the X-ray crystal structures of the BACE1-inhibitor complexes, 2B8L, 1M4H, 1W51 and 2IQG, respectively.
For the reasons mentioned above, we performed a docking calculation using the X-ray crystal structure 1M4H, in which the P2 moiety of the inhibitor (OM00-3) had a similar size to that of our inhibitor (Table 3). OM00-3 appears to show a high docking score value owing to its large molecular size: OM00-3 has many more amide bonds than our peptapeptidic inhibitors, and can closely interact with the BACE1 active site
Compound | Xaa | U_ele1 | U_vdw2 | U_str3 | U_dock4 | BACE1 inhibition (%)5 |
---|---|---|---|---|---|---|
−233.4 | 0.7 | 15.7 | −217.0 | — | ||
L | −195.7 | 11.3 | 10.3 | −174.1 | 83.7 | |
W | −189.3 | 13.0 | 10.5 | −165.8 | 71.3 | |
I | −194.6 | 18.5 | 11.1 | −165.0 | 65.8 | |
F | −195.6 | 23.0 | 11.7 | −160.9 | 47.3 | |
M | −195.6 | 24.6 | 10.1 | −160.9 | 63.9 | |
Cha6 | −194.4 | 26.6 | 11.6 | −156.2 | 84.1 | |
Y | −196.7 | 26.5 | 14.1 | −156.1 | 36.9 | |
E | −190.7 | 24.2 | 10.6 | −155.9 | 36.3 | |
N | −190.7 | 24.2 | 11.3 | −155.2 | 28.4 | |
Q | −196.9 | 27.6 | 15.7 | −153.6 | 14.1 | |
D | −185.5 | 28.5 | 9.7 | −147.3 | 33.2 |
Table 3.
Docking simulation of pentapeptidic BACE1 inhibitors and their scoring using the 1M4H X-ray crystal structure.
1Electrostatic energy between BACE1 and inhibitor (kcal/mol).
2van der Waals energy between BACE1 and inhibitor (kcal/mol).
3Strain energy of inhibitor (kcal/mol).
4Docking score (kcal/mol); U_dock = U_ele + U_vdw + U_str.
5BACE1 inhibition % at 2 μM.
6Cyclohexylalanine (Cha).

Figure 5.
The correlation between BACE1 inhibitory activities and the docking score values.
Furthermore, we designed a series of BACE1 inhibitors that possessed one or more bioisosteres of carboxylic acid from pentapeptide

Figure 6.
Design of practical BACE1 inhibitors using a carboxylic acid bioisostere.
3. Design of small-sized non-peptidic BACE1 inhibitors
At present, many non-peptidic BACE1 inhibitors have been discovered. The research of Elan pharmaceuticals, Merck (MSD), Pfizer, and Schering-Plough, and Ghosh et al. reported the BACE1 inhibitors

Figure 7.
BACE1 inhibitors with an isophthalic scaffold.

Figure 8.
Design of small-sized non-peptidic BACE1 inhibitors with a 2,6-pyridinedicarboxylic scaffold.
Although
Based on the electron-donor bioisostere concept, we speculated that an electron-rich halogen atom could interact with the electron-poor guanidine π-orbital by Coulomb interaction. Using the
4. Conclusion
Although we calculated the docking scoring values by a docking calculation between BACE1 and its inhibitors using the first reported X-ray crystal structure 1FKN, we found no association between the scoring values and BACE1 inhibition. We found that a specific interaction, a quantum chemical interaction between the Arg235 side chain and the P2 region of the inhibitor, played a critical role in the inhibitory mechanism of BACE1. Whereas most BACE1 inhibitors, except OM99-2, interacted with BACE1-Arg235 by a quantum chemical interaction, such as stacking and σ-π interaction, many early BACE1 inhibitors were designed using the 1FKN coordinate set. As the crystal structure 1FKN showed a hydrogen bonding between the BACE1-Arg235 side chain and OM99-2, the early studies on BACE1 inhibitor design might have misdirected, as a docking simulation using 1FKN appears to be meaningless. In fact, unlike 1FKN, there is no hydrogen bonding interaction present in most of the X-ray crystal structures. We selected the peptide sequence that showed potent inhibitory activity by a docking simulation using the X-ray crystal structure 1M4H, and designed potent peptidic BACE1 inhibitors with one or more carboxylic acid bioisosteres. Moreover, we focused on a quantum chemical interaction, and designed the potent non-peptidic BACE1 inhibitor
Acknowledgments
This study was supported in part by the Grants-in-Aid for Scientific Research from MEXT (Ministry of Education, Culture, Sports, Science and Technology), Japan (KAKENHI No. 23590137 and No. 26460163), and a donation from Professor Emeritus Tetsuro Fujita of Kyoto University. At the time of writing, we received word that Prof. Fujita had passed away on January 1, 2017. Prof. Fujita was the teacher of one of the authors, Y. Hamada, and was known as the inventor of a treatment agent for multiple sclerosis. We dedicate this article to Prof. Fujita.
References
- 1.
Selkoe DJ. Toward a comprehensive theory for Alzheimer’s disease. Hypothesis: Alzheimer’s disease is caused by the cerebral accumulation and cytotoxicity of amyloid β-protein. Annals of the New York Academy of Sciences. 2000; 924 :17-25. DOI: 10.1111/j.1749-6632.2000.tb05554.x - 2.
Selkoe DJ. The deposition of amyloid proteins in the aging mammalian brain: Implications for Alzheimer’s disease. Annals of Medicine. 1989; 21 :73-76. DOI: 10.3109/07853898909149187 - 3.
Selkoe DJ. Translating cell biology into therapeutic advances in Alzheimer’s disease. Nature. 1999; 399 :A23-A31. DOI: 10.1038/399a023 - 4.
Sinha S, Lieberburg I. Cellular mechanisms of β-amyloid production and secretion. Proceedings of the National Academy of Sciences of the United States of America. 1999; 96 :11049-11053. DOI: 10.1073/pnas.96.20.11049 - 5.
Vassar R, Bennett BD, Babu-Khan S, et al. β-secretase cleavage of Alzheimer’s amyloid precursor protein by the transmembrane aspartic protease BACE. Science. 1999: 286 ;735-741. DOI: 10.1523/JNEUROSCI.3657-09.2009 - 6.
Hussain I, Powell D, Howlett DR, et al. Identification of a novel aspartic protease (Asp 2) as β-secretase. Neuroscience. 1999: 14 ;419-427. DOI: 10.1523/JNEUROSCI.3657-09.2009 - 7.
Yan R, Bienkowski MJ, Shuck ME, et al. Membrane-anchored aspartyl protease with Alzheimer’s disease β-secretase activity. Nature. 1999: 402 ;533-537. DOI: 10.1038/990107 - 8.
Turner PR, O’Connor K, Tate WP, Abraham WC. Roles of amyloid precursor protein and its fragment in regulating neural activity, plasticity and memory. Progress in Neurobiology. 2003: 70 ;1-32. DOI: 10.1016/S0301-0082(03)00089-3 - 9.
Roberds SL, Anderson J, Basi G, et al. BACE knockout mice are healthy despite lacking the primary β-secretase activity in brain: Implications for Alzheimer's disease therapeutics. Human Molecular Genetics. 2001: 10 ;1317-1324. DOI: 10.1093/hmg/10.12.1317 - 10.
Hamada Y, Kiso Y. New directions for protease inhibitors directed drug discovery. Biopolymers. 2016: 106 ;563-579. DOI: 10.1002/bip.22780 - 11.
Hamada Y, Kiso Y. Aspartic protease inhibitors as drug candidates for treating various difficult-to-treat diseases. In: Ryadnov M, Farkas E, editors. Amino Acids, Peptides and Proteins. Vol. 39. London: Royal Society of Chemistry; 2015. pp. 114-147. DOI: 10.1039/9781849739962-00114 - 12.
Hamada Y. Drug discovery of b-secretase inhibitors based on quantum chemical interactions for the treatment of Alzheimer’s disease. SOJ Pharmacy & Pharmaceutical Sciences. 2014: 1 ;1-8. DOI: 10.15226/2374-6866/1/3/00118 - 13.
Hamada Y, Kiso Y. Recent progress in the drug discovery of non-peptidic BACE1 inhibitors. Expert Opinion on Drug Discovery. 2009: 4 ;391-416. DOI: 10.1517/17460440902806377 - 14.
Hamada Y, Kiso Y. The application of bioisosteres in drug design for novel drug discovery: Focusing on acid protease inhibitors. Expert Opinion on Drug Discovery. 2012: 7 ;903-922. DOI: 10.1517/17460441.2012.712513 - 15.
Hamada Y, Kiso Y. Advances in the identification of β-secretase inhibitors. Expert Opinion on Drug Discovery. 2013: 8 ;709-731. DOI: 10.1517/17460441.2013.784267 - 16.
Nguyen J-T, Hamada Y, Kimura T, Kiso Y. Design of potent aspartic protease inhibitors to treat various diseases. Archiv der Pharmazie—Chemistry in Life Sciences. 2008: 341 ;523-535. DOI: 10.1002/ardp.200700267 - 17.
Sinha S, Anderson JP, Barbour R, et al. Purification and cloning of amyloid precursor protein β-secretase from human brain. Nature. 1999: 402 ;537-540. DOI: 10.1038/990114 - 18.
Ghosh AK, Shin D, Downs D, et al. Design of potent inhibitors for human brain memapsin 2 (β-secretase). Journal of the American Chemical Society. 2000: 122 ;3522-3523. DOI: 10.1021/ja000300g - 19.
Hong L, Koelsch G, Lin X, et al. Structure of the protease domain of memapsin 2 (β-secretase) complexed with inhibitor. Science. 2000: 290 ;150-153. DOI: 10.1126/science.290.5489.150 - 20.
Ghosh AK, Bilcer G, Harwood, C et al. Structure-based design: Potent inhibitors of human brain memapsin 2 (β-secretase). Journal of Medicinal Chemistry. 2001: 44 ;2865-2868. DOI: 10.1021/jm0101803 - 21.
Ghosh AK, Devasamudram T, Hong L, et al. Structure-based design of cycloamide-urethane-derived novel inhibitors of human brain memapsin 2 (β-secretase). Bioorganic & Medicinal Chemistry Letters. 2005: 15 ;15-20. DOI: 10.1016/j.bmcl.2004.10.084 - 22.
Shuto D, Kasai K, Hidaka K, et al. KMI-008, a novel β-secretase inhibitor containing a hydroxymethylcarbonyl isostere as a transition-state mimic: Design and synthesis of substrate-based octapeptides. Bioorganic & Medicinal Chemistry Letters. 2003: 13 ;4273-4276. DOI: 10.1016/j.bmcl.2003.09.053 - 23.
Kimura T, Shuto D, Kasai K, et al. KMI-358 and KMI-370, highly potent and small-sized BACE1 inhibitors containing phenylnorstatine. Bioorganic & Medicinal Chemistry Letters. 2004: 14 ;1527-1531. DOI: 10.1016/j.bmcl.2003.12.088 - 24.
Kimura T, Shuto D, Hamada Y, et al. Design and synthesis of highly active Alzheimer’s β-secretase (BACE1) inhibitors, KMI-420 and KMI-429, with enhanced chemical stability. Bioorganic & Medicinal Chemistry Letters. 2005: 15 ;211-215. DOI: 10.1016/j.bmcl.2004.09.090 - 25.
Asai M, Hattori C, Iwata N, et al. The novel β-secretase inhibitor KMI-429 reduces amyloid beta peptide production in amyloid precursor protein transgenic and wild-type mice. Journal of Neurochemistry. 2006: 96 ;533-540. DOI: 10.1111/j.1471-4159.2005.03576.x - 26.
Kimura T, Hamada Y, Stochaj M, et al. Design and synthesis of potent β-secretase (BACE1) inhibitors with P1’ carboxylic acid bioisostere. Bioorganic & Medicinal Chemistry Letters. 2006: 16 ;2380-2386. DOI: 10.1016/j.bmcl.2006.01.108 - 27.
Hamada Y, Igawa N, Ikari H, et al. β-Secretase inhibitors: Modification at the P4 position and improvement of inhibitory activity in cultured cells. Bioorganic & Medicinal Chemistry Letters. 2006: 16 ;4354-4359. DOI: 10.1016/j.bmcl.2006.05.046 - 28.
Hamada Y, Abdel-Rahman H, Yamani A, et al. BACE1 inhibitors: Optimization by replacing the P1’ residue with non-acidic moiety. Bioorganic & Medicinal Chemistry Letters. 2008: 18 ;1649-1653. DOI: 10.1016/j.bmcl.2008.01.058 - 29.
Tagad HD, Hamada Y, Nguyen J-T, et al. Design of pentapeptidic BACE1 inhibitors with carboxylic acid bioisosteres at P1' and P4 positions. Bioorganic & Medicinal Chemistry. 2010: 18 ;3175-3186. DOI: 10.1016/j.bmc.2010.03.032 - 30.
Tagad HD, Hamada Y, Nguyen J-T, et al. Structure-guided design and synthesis of P1’ position 1-phenylcycloalkylamine-derived pentapeptidic BACE1 inhibitors. Bioorganic & Medicinal Chemistry. 2011: 19 ;5238-5246. DOI: 10.1016/j.bmc.2011.07.002 - 31.
Hamada Y, Ohta H, Miyamoto N, et al. Significance of interaction of BACE1-Arg235 with its ligands and design of BACE1 inhibitors with P2 pyridine scaffold. Bioorganic & Medicinal Chemistry Letters. 2009: 19 ;2435-2439. DOI: 10.1016/j.bmcl.2009.03.049 - 32.
Hamada Y, Ohta H, Miyamoto N, et al. Novel non-peptidic and small-sized BACE1 inhibitors. Bioorganic & Medicinal Chemistry Letters. 2008: 18 ;1654-1658. DOI: 10.1016/j.bmcl.2008.01.056 - 33.
Hamada Y, Nakanishi T, Suzuki K, et al. Novel BACE1 inhibitors possessing a 5-nitroisophthalic scaffold at the P2 position. Bioorganic & Medicinal Chemistry Letters. 2012: 22 ;4640-4644. DOI: 10.1016/j.bmcl.2012.05.089 - 34.
Suzuki K, Hamada Y, Nguyen J-T, et al. Novel BACE1 inhibitors with a non-acidic heterocycle at the P1’ position. Bioorganic & Medicinal Chemistry. 2013: 21 ;6665-6673. DOI: 10.1016/j.bmc.2013.08.016 - 35.
Hamada Y, Suzuki K, Nakanishi T, et al. Structure-activity relationship study of BACE1 inhibitors possessing a chelidonic or 2,6-pyridinedicarboxylic scaffold at the P2 position. Bioorganic & Medicinal Chemistry Letters. 2014: 24 ;618-623. DOI: 10.1016/j.bmcl.2013.12.007 - 36.
Imai YN, Inoue Y, Nakanishi I, Kitaura K. Cl–π interactions in protein–ligand complexes. Protein Science. 2008: 17 ;1129-1137. DOI: 10.1110/ps.033910.107