Exploring the binding of BACE-1 inhibitors using comparative binding energy analysis (COMBINE)

Sulfonamido, amido heterocyclic adducts of tetrazole derivatives as BACE1 inhibitors: in silico exploration

Alzheimer's disease is a neurodegenerative disorder accounting for 60-80% of dementia cases and is accompanied by a high mortality rate in patients above 70 years of age. The formation of senile plaques composed of amyloid-β protein is a hallmark of Alzheimer's disease. Beta-site APP cleaving enzyme 1 (BACE1) is a proteolytic enzyme involved in the degradation of amyloid precursor protein, which further degrades to form toxic amyloid-β fragments. Hence, inhibition of BACE1 was stated to be an effective strategy for Alzheimer's therapeutics. Keeping in mind the structures of different BACE1 inhibitors that had reached the clinical trials, we designed a library of compounds (total 164) based on a substituted 5-amino tetrazole scaffold which was an isosteric replacement of the cyclic amidine moiety, a common component of the BACE1 inhibitors which reached the clinical trials. The scaffold was linked to different structural moieties with the aid of an amide or sulfonamide bond to design some novel molecules. Molecular docking was initially performed and the top 5 molecules were selected based on docking scores and protein-ligand interactions. Furthermore, molecular dynamic simulations were performed for these molecules (3g, 7k, 8n, 9d, 9g) for 100 ns and MM-GBSA calculations were performed for each of these complexes. After critical evaluation of the obtained results, three potential molecules (9d, 8n, and 7k) were forwarded for prolonged stability studies by performing molecular dynamic simulations for 250 ns and simultaneous MM-GBSA calculations. It was observed that the compounds (9d, 8n, and 7k) were forming good interactions with the amino acid residues of the catalytic site of the enzyme with multiple non-covalent interactions. In MD simulations, the compounds have shown better stability and better binding energy throughout the runtime.

Identification of BACE-1 Inhibitors against Alzheimer’s Disease through E-Pharmacophore-Based Virtual Screening and Molecular Dynamics Simulation Studies: An Insilco Approach

Alzheimer is a severe memory and cognitive impairment neurodegenerative disease that is the most common cause of dementia worldwide and characterized by the pathological accumulation of tau protein and amyloid-beta peptides. In this study, we have developed E-pharmacophore modeling to screen the eMolecules database with the help of a reported co-crystal structure bound with Beta-Site Amyloid Precursor Protein Cleaving Enzyme 1 (BACE-1). Flumemetamol, florbetaben, and florbetapir are currently approved drugs for use in the clinical diagnosis of Alzheimer’s disease. Despite the benefits of commercially approved drugs, there is still a need for novel diagnostic agents with enhanced physicochemical and pharmacokinetic properties compared to those currently used in clinical practice and research. In the E-pharmacophore modeling results, it is revealed that two aromatic rings (R19, R20), one donor (D12), and one acceptor (A8) are obtained, and also that similar pharmacophoric features of compounds are identified from pharmacophore-based virtual screening. The identified screened hits were filtered for further analyses using structure-based virtual screening and MM/GBSA. From the analyses, top hits such as ZINC39592220 and en1003sfl.46293 are selected based on their top docking scores (−8.182 and −7.184 Kcal/mol, respectively) and binding free energy (−58.803 and −56.951 Kcal/mol, respectively). Furthermore, a molecular dynamics simulation and MMPBSA study were performed, which revealed admirable stability and good binding free energy throughout the simulation period. Moreover, Qikprop results revealed that the selected, screened hits have good drug-likeness and pharmacokinetic properties. The screened hits ZINC39592220 and en1003sfl.46293 could be used to develop drug molecules against Alzheimer’s disease.

Design, Synthesis, and Pharmacological Evaluation of Embelin-Aryl/alkyl Amine Hybrids as Orally Bioavailable Blood-Brain Barrier Permeable Multitargeted Agents with Therapeutic Potential in Alzheimer's Disease: Discovery of SB-1448

The complex and multifaceted nature of Alzheimer's disease has brought about a pressing demand to develop ligands targeting multiple pathways to combat its outrageous prevalence. Embelin is a major secondary metabolite of Embelia ribes Burm f., one of the oldest herbs in Indian traditional medicine. It is a micromolar inhibitor of cholinesterases (ChEs) and β-site amyloid precursor protein cleaving enzyme 1 (BACE-1) with poor absorption, distribution, metabolism, and excretion (ADME) properties. Herein, we synthesize a series of embelin-aryl/alkyl amine hybrids to improve its physicochemical properties and therapeutic potency against targeted enzymes. The most active derivative, 9j (SB-1448), inhibits human acetylcholinesterase (hAChE), human butyrylcholinesterase (hBChE), and human BACE-1 (hBACE-1) with IC₅₀ values of 0.15, 1.6, and 0.6 μM, respectively. It inhibits both ChEs noncompetitively with k_i values of 0.21 and 1.3 μM, respectively. It is orally bioavailable, crosses blood-brain barrier (BBB), inhibits Aβ self-aggregation, possesses good ADME properties, and protects neuronal cells from scopolamine-induced cell death. The oral administration of 9j at 30 mg/kg attenuates the scopolamine-induced cognitive impairments in C57BL/6J mice.

BACE-1 Inhibitors Targeting Alzheimer's Disease

The accumulation of amyloid-β (Aβ) is the main event related to Alzheimer's disease (AD) progression. Over the years, several disease-modulating approaches have been reported, but without clinical success. The amyloid cascade hypothesis evolved and proposed essential targets such as tau protein aggregation and modulation of β-secretase (β-site amyloid precursor protein cleaving enzyme 1 - BACE-1) and γ-secretase proteases. BACE-1 cuts the amyloid precursor protein (APP) to release the C99 fragment, giving rise to several Aβ peptide species during the subsequent γ-secretase cleavage. In this way, BACE-1 has emerged as a clinically validated and attractive target in medicinal chemistry, as it plays a crucial role in the rate of Aβ generation. In this review, we report the main results of candidates in clinical trials such as E2609, MK8931, and AZD-3293, in addition to highlighting the pharmacokinetic and pharmacodynamic-related effects of the inhibitors already reported. The current status of developing new peptidomimetic, non-peptidomimetic, naturally occurring, and other class inhibitors are demonstrated, considering their main limitations and lessons learned. The goal is to provide a broad and complete approach to the subject, exploring new chemical classes and perspectives.

Prediction of the Drug-Target Binding Kinetics for Flexible Proteins by Comparative Binding Energy Analysis

There is growing consensus that the optimization of the kinetic parameters for drug-protein binding leads to improved drug efficacy. Therefore, computational methods have been developed to predict kinetic rates and to derive quantitative structure-kinetic relationships (QSKRs). Many of these methods are based on crystal structures of ligand-protein complexes. However, a drawback is that each ligand-protein complex is usually treated as having a single structure. Here, we present a modification of COMparative BINding Energy (COMBINE) analysis, which uses the structures of ligand-protein complexes to predict binding parameters. We introduce the option of using multiple structures to describe each ligand-protein complex in COMBINE analysis and apply this to study the effects of protein flexibility on the derivation of dissociation rate constants (k_off) for inhibitors of p38 mitogen-activated protein (MAP) kinase, which has a flexible binding site. Multiple structures were obtained for each ligand-protein complex by performing docking to an ensemble of protein configurations obtained from molecular dynamics simulations. Coefficients to scale ligand-protein interaction energies determined from energy-minimized structures of ligand-protein complexes were obtained by partial least squares regression, and they allowed for the computation of k_off values. The QSKR model obtained using single, energy-minimized crystal structures for each ligand-protein complex had higher predictive power than the QSKR model obtained with multiple structures from ensemble docking. However, incorporation of ligand-protein flexibility helped to highlight additional ligand-protein interactions that lead to longer residence times, such as interactions with residues Arg67 and Asp168, which are close to the ligand in many crystal structures. These results show that COMBINE analysis is a promising method to guide the design of compounds that bind to flexible proteins with improved binding kinetics.

A Protocol to Use Comparative Binding Energy Analysis to Estimate Drug-Target Residence Time

Comparative Binding Energy (COMBINE) analysis is an approach for deriving a target-specific scoring function to compute binding free energy, drug-binding kinetics, or a related property by exploiting the information contained in the three-dimensional structures of receptor-ligand complexes. Here, we describe the process of setting up and running COMBINE analysis to derive a Quantitative Structure-Kinetics Relationship (QSKR) for the dissociation rate constants (k_off) of inhibitors of a drug target. The derived QSKR model can be used to estimate residence times (τ, τ=1/k_off) for similar inhibitors binding to the same target, and it can also help to identify key receptor-ligand interactions that distinguish inhibitors with short and long residence times. Herein, we demonstrate the protocol for the application of COMBINE analysis on a dataset of 70 inhibitors of heat shock protein 90 (HSP90) belonging to 11 different chemical classes. The procedure is generally applicable to any drug target with known structural information on its complexes with inhibitors.

Molecular insights into the inhibitory mechanism of bi-functional bis-tryptoline triazole against β-secretase (BACE1) enzyme

The β-site amyloid precursor protein-cleaving enzyme 1 (β-secretase, BACE1) is involved in the formation of amyloid-β (Aβ) peptide that aggregates into soluble oligomers, amyloid fibrils, and plaques responsible for the neurodegeneration in Alzheimer disease (AD). BACE1 is one of the prime therapeutic targets for the design of inhibitors against AD as BACE1 participate in the rate-limiting step in Aβ production. Jiaranaikulwanitch et al. reported bis-tryptoline triazole (BTT) compound as a potent inhibitor against BACE1, Aβ aggregation as well as possessing metal chelation and antioxidant activity. However, the molecular mechanism of BACE1 inhibition by BTT remains unclear. Thus, molecular docking and molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of BTT against BACE1. MD simulations highlight that BTT interact with catalytic aspartic dyad residues (Asp32 and Asp228) and active pocket residues of BACE1. The hydrogen-bond interactions, hydrophobic contacts, and π-π stacking interactions of BTT with flap residues (Val67-Asp77) of BACE1 confine the movement of the flap and help to achieve closed (non-active) conformation. The PCA analysis highlights lower conformational fluctuations for BACE1-BTT complex, which suggests enhanced conformational stability in comparison to apo-BACE1. The results of the present study provide key insights into the underlying inhibitory mechanism of BTT against BACE1 and will be helpful for the rational design of novel inhibitors with enhanced potency against BACE1.

Multi-Target β-Protease Inhibitors from Andrographis paniculata: In Silico and In Vitro Studies

Natural products derived from plants play a vital role in the discovery of new drug candidates, and these are used for novel therapeutic drug development. Andrographis paniculata and Spilanthes paniculata are used extensively as medicinal herbs for the treatment of various ailments, and are reported to have neuroprotective properties. β-amyloid is a microscopic brain protein whose significant aggregation is detected in mild cognitive impairment and Alzheimer’s disease (AD) brains. The accumulation of β-amyloid disrupts cell communication and triggers inflammation by activating immune cells, leading to neuronal cell death and cognitive disabilities. The proteases acetylcholinesterase (AChE), butyrylcholinesterase (BChE), and beta secretase-1 (BACE-1) have been reported to be correlated with the synthesis and growth of β-amyloid plaques in the brains of AD patients. In the present study, the phenolic compounds from A. paniculata and S. paniculata that have been reported in the literature were selected for the current investigation. Furthermore, we employed molecular docking and molecular dynamics studies of the phenolic compounds with the proteins AChE, BChE, and BACE-1 in order to evaluate the binding characteristics and identify potent anti-amyloid agents against the neurodegenerative diseases such as AD. In this investigation, we predicted three compounds from A. paniculata with maximum binding affinities with cholinesterases and BACE-1. The computational investigations predicted that these compounds follow the rule of five. We further evaluated these molecules for in vitro inhibition activity against all the enzymes. In the in vitro investigations, 3,4-di-o-caffeoylquinic acid (5281780), apigenin (5280443), and 7-o-methylwogonin (188316) were found to be strong inhibitors of AChE, BChE, and BACE-1. These findings suggest that these compounds can be potent multi-target inhibitors of the proteases that might cumulatively work and inhibit the initiation and formation of β-amyloid plaques, which is a prime cause of neurotoxicity and dementia. According to our knowledge, these findings are the first report on natural compounds isolated from A. paniculata as multi-target potent inhibitors and anti-amyloid agents.

Hybrid approach to sieve out natural compounds against dual targets in Alzheimer's Disease

Excess Aβ production by the key protease BACE1, results in Aβ aggregation, forming amyloid plaques, all of which contribute to the pathogenesis of Alzheimer’s disease. Besides the multi-factorial nature of the disease, the diversity in the size and shape of known ligands that bind to the active site of BACE1, that is the flexibility of the enzyme, pose a serious challenge for the identification of drug candidates. To address the issue of receptor flexibility we have carried out ensemble docking with multiple receptor conformations. Therein, two representative structures each from closed and semi-open BACE1 conformations were selected for virtual screening to identify compounds that bind to the active site of both the conformations. These outperformed compounds were ranked using pharmacophore models generated by a ligand-based approach, for the identification of BACE1 inhibitors. The inhibitors were further predicted for anti-amyloidogenic activity using a QSAR model already established by our group thus enlisting compounds with dual potency. BACE1 inhibitory and anti-amyloidogenic activity for the commercially available compounds were validated using in vitro studies. Thus, incorporation of receptor flexibility in BACE1 through ensemble docking in conjunction with structure and ligand-based approach for screening might act as an effective protocol for obtaining promising scaffolds against AD.

Multidisciplinary approaches for targeting the secretase protein family as a therapeutic route for Alzheimer's disease

The continual increase of the aging population worldwide renders Alzheimer's disease (AD) a global prime concern. Several attempts have been focused on understanding the intricate complexity of the disease's development along with the on- andgoing search for novel therapeutic strategies. Incapability of existing AD drugs to effectively modulate the pathogenesis or to delay the progression of the disease leads to a shift in the paradigm of AD drug discovery. Efforts aimed at identifying AD drugs have mostly focused on the development of disease-modifying agents in which effects are believed to be long lasting. Of particular note, the secretase enzymes, a group of proteases responsible for the metabolism of the β-amyloid precursor protein (βAPP) and β-amyloid (Aβ) peptides production, have been underlined for their promising therapeutic potential. This review article attempts to comprehensively cover aspects related to the identification and use of drugs targeting the secretase enzymes. Particularly, the roles of secretases in the pathogenesis of AD and their therapeutic modulation are provided herein. Moreover, an overview of the drug development process and the contribution of computational (in silico) approaches for facilitating successful drug discovery are also highlighted along with examples of relevant computational works. Promising chemical scaffolds, inhibitors, and modulators against each class of secretases are also summarized herein. Additionally, multitarget secretase modulators are also taken into consideration in light of the current growing interest in the polypharmacology of complex diseases. Finally, challenging issues and future outlook relevant to the discovery of drugs targeting secretases are also discussed.

Alzheimer: A Decade of Drug Design. Why Molecular Topology can be an Extra Edge?

BACKGROUND

The last decade was characterized by a growing awareness about the severity of dementia in the field of age-related and no age-related diseases and about the importance to invest resources in the research of new, effective treatments. Among the dementias, Alzheimer's plays a substantial role because of its extremely high incidence and fatality. Several pharmacological strategies have been tried but still now, Alzheimer keeps being an untreatable disease. In literature, the number of QSAR related drug design attempts about new treatments for Alzheimer is huge, but only few results can be considered noteworthy. Providing a detailed analysis of the actual situation and reporting the most notable results in the field of drug design and discovery, the current review focuses on the potential of molecular topology as a reliable tool in finding new anti-Alzheimer lead compounds.

METHODS

Published works on QSAR applied to the search of anti-Alzheimer's drugs during the last 10 years has been tracked. 2D and 3D-QSAR, HQSAR, topological indexes, etc. have been analyzed, as well as different mechanisms of action, such as MAO, AchE, etc. An example of topological indexes' application to the search of potential anti-Alzheimer drugs is reported.

RESULTS

Results show that QSAR methods during the last decade represented an excellent approach to the search of new effective drugs against Alzheimer's. In particular, QSAR based on molecular topology allows the establishment of a direct structure-property link that results in the identification of new hits and leads.

CONCLUSION

Molecular topology is a powerful tool for the discovery of new anti-Alzheimer drugs covering simultaneously different mechanisms of action, what may help to find a definitive cure for the disease.

In Silico Studies in Drug Research Against Neurodegenerative Diseases

Background

Neurodegenerative diseases such as Alzheimer's disease (AD), amyotrophic lateral sclerosis, Parkinson's disease (PD), spinal cerebellar ataxias, and spinal and bulbar muscular atrophy are described by slow and selective degeneration of neurons and axons in the central nervous system (CNS) and constitute one of the major challenges of modern medicine. Computer-aided or in silico drug design methods have matured into powerful tools for reducing the number of ligands that should be screened in experimental assays.

Methods

In the present review, the authors provide a basic background about neurodegenerative diseases and in silico techniques in the drug research. Furthermore, they review the various in silico studies reported against various targets in neurodegenerative diseases, including homology modeling, molecular docking, virtual high-throughput screening, quantitative structure activity relationship (QSAR), hologram quantitative structure activity relationship (HQSAR), 3D pharmacophore mapping, proteochemometrics modeling (PCM), fingerprints, fragment-based drug discovery, Monte Carlo simulation, molecular dynamic (MD) simulation, quantum-mechanical methods for drug design, support vector machines, and machine learning approaches.

Results

Detailed analysis of the recently reported case studies revealed that the majority of them use a sequential combination of ligand and structure-based virtual screening techniques, with particular focus on pharmacophore models and the docking approach.

Conclusion

Neurodegenerative diseases have a multifactorial pathoetiological origin, so scientists have become persuaded that a multi-target therapeutic strategy aimed at the simultaneous targeting of multiple proteins (and therefore etiologies) involved in the development of a disease is recommended in future.

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.

Investigation of naphthofuran moiety as potential dual inhibitor against BACE-1 and GSK-3β: molecular dynamics simulations, binding energy, and network analysis to identify first-in-class dual inhibitors against Alzheimer's disease

BACE-1 and GSK-3β are potential therapeutic drug targets for Alzheimer's disease. Recently, both the targets received attention for designing dual inhibitors for Alzheimer's disease. Until now, only two-scaffold triazinone and curcumin have been reported as BACE-1 and GSK-3β dual inhibitors. Docking, molecular dynamics, clustering, binding energy, and network analysis of triazinone derivatives with BACE-1 and GSK-3β was performed to get molecular insight into the first reported dual inhibitor. Further, we designed and evaluated a naphthofuran series for its ability to inhibit BACE-1 and GSK-3β with the computational approaches. Docking study of naphthofuran series showed a good binding affinity towards both the targets. Molecular dynamics, binding energy, and network analysis were performed to compare their binding with the targets and amino acids responsible for binding. Naphthofuran series derivatives showed good interaction within the active site residues of both of the targets. Hydrogen bond occupancy and binding energy suggested strong binding with the targets. Dual-inhibitor binding was mostly governed by the hydrophobic interactions for both of the targets. Per residue energy decomposition and network analysis identified the key residues involved in the binding and inhibiting BACE-1 and GSK-3β. The results indicated that naphthofuran series derivative 11 may be a promising first-in-class dual inhibitor against BACE-1 and GSK-3β. This naphthofuran series may be further explored to design better dual inhibitors. Graphical abstract Naphthofuran derivative as a dual inhibitor for BACE-1 and GSK-3β.

Identification of new BACE1 inhibitors using Pharmacophore and Molecular dynamics simulations approach

Inhibition of β-Secretase (BACE1) is crucial for the treatment of Alzheimer's disease (AD). Availability of BACE1 crystal structures in both apo and complexed forms enables to find structure-based BACE1 inhibitors for controlling AD. There are two catalytic aspartates (ASP32 and ASP228) presents in the active domain of BACE1. In order to understand the binding mechanism and structure-activity relationship of amidine-containing BACE1 inhibitors, molecular docking, and pharmacophore and 3D-QSAR studies have been carried out with 34 amidine derivatives to develop a pharmacophore model. Pharmacophore-based virtual screening (PBVS) has been performed against BACE1 (PDB ID: 2FDP), using three chemical databases (CoCoCo, Enamine, Zinc), which yielded 6000 hit compounds. These compounds were further analyzed using structure-based docking in hierarchical filtering approaches of Glide such as HTVS, SP, and XP precision modes. The docking results show that binding orientations of the inhibitors at Asp dyad active site amino acid residues of β-Secretase. Results from glide XP docking and induced fit docking showed that four leads (Lead1, Lead3, Lead4 and Lead5) have good interactions with the target protein in comparison with cocrystal (amino-ethylene BACE1 inhibitor). Further, molecular dynamics (MD) simulation for these leads bound with BACE1 shows conformational stability and difference in dynamical flap behaviors of the active site with cocrystal inhibitor. Binding free energetic using MM-GB/SA approaches suggest lead 1 and lead 3 has comparably favorable binding to cocrystal inhibitor. Thus, the present study emphasizes these leads for an effective drug to treat Alzheimer disease.

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

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.

Structural basis for the development of SARS 3CL protease inhibitors from a peptide mimic to an aza-decaline scaffold

Design of inhibitors against severe acute respiratory syndrome (SARS) chymotrypsin-like protease (3CL(pro) ) is a potentially important approach to fight against SARS. We have developed several synthetic inhibitors by structure-based drug design. In this report, we reveal two crystal structures of SARS 3CL(pro) complexed with two new inhibitors based on our previous work. These structures combined with six crystal structures complexed with a series of related ligands reported by us are collectively analyzed. To these eight complexes, the structural basis for inhibitor binding was analyzed by the COMBINE method, which is a chemometrical analysis optimized for the protein-ligand complex. The analysis revealed that the first two latent variables gave a cumulative contribution ratio of r(2)  = 0.971. Interestingly, scores using the second latent variables for each complex were strongly correlated with root mean square deviations (RMSDs) of side-chain heavy atoms of Met(49) from those of the intact crystal structure of SARS-3CL(pro) (r = 0.77) enlarging the S2 pocket. The substantial contribution of this side chain (∼10%) for the explanation of pIC50 s was dependent on stereochemistry and the chemical structure of the ligand adapted to the S2 pocket of the protease. Thus, starting from a substrate mimic inhibitor, a design for a central scaffold for a low molecular weight inhibitor was evaluated to develop a further potent inhibitor. © 2015 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 106: 391-403, 2016.

Study on the active mechanism of β-secretase inhibitors by molecular simulations

The proteolytic enzyme β-secretase (BACE-1) is one of potential drug targets for treating Alzheimers's disease. First, the reliable and accurate models of three-dimensional quantitative structure-activity relationship for the BACE-1 inhibitors were established, and the several important structural factors that mainly influence the inhibitory activity were obtained. Second, the results of molecular docking presented the binding mode between BACE-1 and its inhibitors, and molecular dynamic simulations provided the details of the receptor-ligand interactions. Furthermore, several new derivatives were designed and validated based on these theoretical analyses. Our studies revealed the binding mechanism between BACE-1 and its inhibitors, and provide some insights into the further structural modification and the design of new inhibitors with higher activity.