Computational modeling of substrate specificity and catalysis of the β-secretase (BACE1) enzyme

Rational Design of Dual Inhibitors for Alzheimer's Disease: Insights from Computational Screening of BACE1 and GSK-3β

BACKGROUND

Alzheimer's disease (AD) is one of the most concerned neurodegenerative disorders across the world characterized by amyloid-beta (Aβ) plaques and neurofibrillary tangles (NFTs), leading to cognitive decline and memory loss. Targeting key pathways involved in AD like Aβ and NFT pathways, are crucial for the development of effective therapeutic strategies. In this study, we aimed to identify and establish promising dual inhibitors targeting BACE1 and GSK-3β, two proteins implicated in Aβ and NFT formation respectively.

METHODS

We have used molecular docking, ADME property analysis, and MMGBSA calculations for the identification of hit molecules and further evaluation of binding affinity, drug-like properties, and stability against BACE1 and GSK-3β.

RESULTS

Our results demonstrated strong binding affinities of ZINC000034853956 towards the active sites of both proteins, with favorable interactions involving key residues crucial for inhibitory activity. Additionally, ZINC000034853956 exhibited favorable drug-like properties. MD simulations revealed the stable binding of ZINC000034853956 to both BACE1 and GSK-3β over a 50 ns period, with consistent ligand-protein interactions, such as hydrogen bonding and hydrophobic contacts. These findings highlight the potential of ZINC000034853956 as a promising candidate for AD treatment, acting as a dual inhibitor targeting both BACE1 and GSK-3β. Overall, our study provides valuable insights into the potential of ZINC000034853956 as a dual inhibitor for AD. The strong binding affinity, favorable drug-like properties, and stability observed in MD simulations support its suitability for further optimization and preclinical studies.

CONCLUSION

Further investigations are warranted to elucidate the precise molecular mechanisms and therapeutic benefits of ZINC000034853956. Our findings offer hope for the development of novel therapeutic interventions targeting crucial pathways involved in AD neurodegeneration.

Differential effects of familial Alzheimer's disease-causing mutations on amyloid precursor protein (APP) trafficking, proteolytic conversion, and synaptogenic activity

The amyloid precursor protein (APP) is a key player in Alzheimer`s disease (AD) and the precursor of the Aβ peptide, which is generated by consecutive cleavages of β- and γ-secretases. Familial Alzheimer's disease (FAD) describes a hereditary subgroup of AD that represents a low percentage of AD cases with an early onset of the disease. Different APP FAD mutations are thought to have qualitatively different effects on its proteolytic conversion. However, few studies have explored the pathogenic and putative physiological differences in more detail. Here, we compared different FAD mutations, located at the β- (Swedish), α- (Flemish, Arctic, Iowa) or γ-secretase (Iberian) cleavage sites. We examined heterologous expression of APP WT and FAD mutants in non-neuronal cells and their impact on presynaptic differentiation in contacting axons of co-cultured neurons. To decipher the underlying molecular mechanism, we tested the subcellular localization, the endocytosis rate and the proteolytic processing in detail by immunoprecipitation-mass spectrometry. Interestingly, we found that only the Iberian mutation showed altered synaptogenic function. Furthermore, the APP Iowa mutant shows significantly decreased α-secretase processing which is in line with our results that APP carrying the Iowa mutation was significantly increased in early endosomes. However, most interestingly, immunoprecipitation-mass spectrometry analysis revealed that the amino acid substitutions of APP FAD mutants have a decisive impact on their processing reflected in altered Aβ profiles. Importantly, N-terminally truncated Aβ peptides starting at position 5 were detected preferentially for APP Flemish, Arctic, and Iowa mutants containing amino acid substitutions around the α-secretase cleavage site. The strongest change in the ratio of Aβ40/Aβ42 was observed for the Iberian mutation while APP Swedish showed a substantial increase in Aβ1-17 peptides. Together, our data indicate that familial AD mutations located at the α-, β-, and γ-secretase cleavage sites show considerable differences in the underlying pathogenic mechanisms.

Structural and molecular insights into tacrine-benzofuran hybrid induced inhibition of amyloid-β peptide aggregation and BACE1 activity

Amyloid-β (Aβ) aggregation and β-amyloid precursor protein cleaving enzyme 1 (BACE1) are the potential therapeutic drug targets for Alzheimer's disease (AD). A recent study highlighted that tacrine-benzofuran hybrid C1 displayed anti-aggregation activity against Aβ42 peptide and inhibit BACE1 activity. However, the inhibition mechanism of C1 against Aβ42 aggregation and BACE1 activity remains unclear. Thus, molecular dynamics (MD) simulations of Aβ42 monomer and BACE1 with and without C1 were performed to inspect the inhibitory mechanism of C1 against Aβ42 aggregation and BACE1 activity. In addition, a ligand-based virtual screening followed by MD simulations was employed to explore potent new small-molecule dual inhibitors of Aβ42 aggregation and BACE1 activity. MD simulations highlighted that C1 promotes the non aggregating helical conformation in Aβ42 and destabilizes D23-K28 salt bridge that plays a vital role in the self-aggregation of Aβ42. C1 displays a favourable binding free energy (-50.7 ± 7.3 kcal/mol) with Aβ42 monomer and preferentially binds to the central hydrophobic core (CHC) residues. MD simulations highlighted that C1 strongly interacted with the BACE1 active site (Asp32 and Asp228) and active pockets. The scrutiny of interatomic distances among key residues of BACE1 highlighted the close flap (non-active) position in BACE1 on the incorporation of C1. The MD simulations explain the observed high inhibitory activity of C1 against Aβ aggregation and BACE1 in the in vitro studies. The ligand-based virtual screening followed by MD simulations identified CHEMBL2019027 (C2) as a promising dual inhibitor of Aβ42 aggregation and BACE1 activity.Communicated by Ramaswamy H. Sarma.

Allostery Inhibition of BACE1 by Psychotic and Meroterpenoid Drugs in Alzheimer's Disease Therapy

In over a century since its discovery, Alzheimer’s disease (AD) has continued to be a global health concern due to its incurable nature and overwhelming increase among older people. In this paper, we give an overview of the efforts of researchers towards identifying potent BACE1 exosite-binding antibodies and allosteric inhibitors. Herein, we apply computer-aided drug design (CADD) methods to unravel the interactions of some proposed psychotic and meroterpenoid BACE1 allosteric site inhibitors. This study is aimed at validating the allosteric potentials of these selected compounds targeted at BACE1 inhibition. Molecular docking, molecular dynamic (MD) simulations, and post-MD analyses are carried out on these selected compounds, which have been experimentally proven to exhibit allosteric inhibition on BACE1. The SwissDock software enabled us to identify more than five druggable pockets on the BACE1 structural surface using docking. Besides the active site region, a melatonin derivative (compound 1) previously proposed as a BACE1 allostery inhibitor showed appreciable stability at eight different subsites on BACE1. Refinement with molecular dynamic (MD) simulations shows that the identified non-catalytic sites are potential allostery sites for compound 1. The allostery and binding mechanism of the selected potent inhibitors show that the smaller the molecule, the easier the attachment to several enzyme regions. This finding hereby establishes that most of these selected compounds failed to exhibit strong allosteric binding with BACE1 except for compound 1. We hereby suggest that further studies and additional identification/validation of other BACE1 allosteric compounds be done. Furthermore, this additional allosteric site investigation will help in reducing the associated challenges with designing BACE1 inhibitors while exploring the opportunities in the design of allosteric BACE1 inhibitors.

Screening Platform Based on Inductively Coupled Plasma Mass Spectrometry for β-Site Amyloid Protein Cleaving Enzyme 1 (BACE1) Inhibitors

β-Site amyloid protein cleaving enzyme 1 (BACE1) is a promising therapeutic target for developing inhibitors to alleviate Alzheimer's disease (AD). Herein, we established an inductively coupled plasma mass spectrometry (ICPMS)-based inhibitor screening platform. A biotin-labeled lanthanide-coded peptide probe (LCPP; biotin-PEG₂-EVNLDAEC-DOTA-Ln) was designed to determine the activity of BACE1 and evaluate the degree of inhibition of inhibitors. The platform was first validated with two commercially available inhibitors (BSI I and BSI IV) in terms of IC₅₀ values and then applied to two newly designed inhibitors (inhibitors II and III) based on the crystal structure of BACE1 interacting with inhibitor I, and each of them contained an acylguanidine core structure. We found that their inhibition effects were improved as evaluated by the sensitive and accurate LCPP-ICPMS platform, demonstrating its ability for new drug screening.

An automated protocol for modelling peptide substrates to proteases

BACKGROUND

Proteases are key drivers in many biological processes, in part due to their specificity towards their substrates. However, depending on the family and molecular function, they can also display substrate promiscuity which can also be essential. Databases compiling specificity matrices derived from experimental assays have provided valuable insights into protease substrate recognition. Despite this, there are still gaps in our knowledge of the structural determinants. Here, we compile a set of protease crystal structures with bound peptide-like ligands to create a protocol for modelling substrates bound to protease structures, and for studying observables associated to the binding recognition.

RESULTS

As an application, we modelled a subset of protease-peptide complexes for which experimental cleavage data are available to compare with informational entropies obtained from protease-specificity matrices. The modelled complexes were subjected to conformational sampling using the Backrub method in Rosetta, and multiple observables from the simulations were calculated and compared per peptide position. We found that some of the calculated structural observables, such as the relative accessible surface area and the interaction energy, can help characterize a protease's substrate recognition, giving insights for the potential prediction of novel substrates by combining additional approaches.

CONCLUSION

Overall, our approach provides a repository of protease structures with annotated data, and an open source computational protocol to reproduce the modelling and dynamic analysis of the protease-peptide complexes.

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.

Changing Paradigm from one Target one Ligand Towards Multi-target Directed Ligand Design for Key Drug Targets of Alzheimer Disease: An Important Role of In Silico Methods in Multi-target Directed Ligands Design

Alzheimer disease (AD) is now considered as a multifactorial neurodegenerative disorder and rapidly increasing to an alarming situation and causing higher death rate. One target one ligand hypothesis does not provide complete solution of AD due to multifactorial nature of the disease and one target one drug fails to provide better treatment against AD. Moreo-ver, currently available treatments are limited and most of the upcoming treatments under clinical trials are based on modulat-ing single target. So, the current AD drug discovery research is shifting towards a new approach for a better solution that simultaneously modulates more than one targets in the neurodegenerative cascade. This can be achieved by network pharma-cology, multi-modal therapies, multifaceted, and/or the more recently proposed term “multi-targeted designed drugs”. Drug discovery project is a tedious, costly and long-term project. Moreover, multi-target AD drug discovery added extra challeng-es such as the good binding affinity of ligands for multiple targets, optimal ADME/T properties, no/less off-target side effect and crossing of the blood-brain barrier. These hurdles may be addressed by insilico methods for an efficient solution in less time and cost as computational methods successfully applied to single target drug discovery project. Here, we are summariz-ing some of the most prominent and computationally explored single targets against AD and further, we discussed a success-ful example of dual or multiple inhibitors for same targets. Moreover, we focused on ligand and structure-based computa-tional approach to design MTDL against AD. However, it is not an easy task to balance dual activity in a single molecule but computational approach such as virtual screening docking, QSAR, simulation and free energy is useful in future MTDLs drug discovery alone or in combination with a fragment-based method. However, rational and logical implementations of computational drug designing methods are capable of assisting AD drug discovery and play an important role in optimizing multi-target drug discovery.

Computational insights into β-site amyloid precursor protein enzyme 1 (BACE1) inhibition by tanshinones and salvianolic acids from Salvia miltiorrhiza via molecular docking simulations

The rhizome of Salvia miltiorrhiza has emerged as a rich source of natural therapeutic agents, and its several compounds are supposed to exhibit favorable effects on Alzheimer's disease (AD). The present work investigate the anti-AD potentials of 12 tanshinones, three salvianolic acids and three caffeic acid derivatives from S. miltiorrhiza via the inhibition of β-site amyloid precursor protein cleaving enzyme 1 (BACE1). Among the tested compounds, deoxyneocryptotanshinone (1), salvianolic acid A (13) and salvianolic acid C (15) displayed good inhibitory effect on BACE1 with IC₅₀ values of 11.53 ± 1.13, 13.01 ± 0.32 and 9.18 ± 0.03 μM, respectively. Besides this, enzyme kinetic analysis on BACE1 revealed 13, a competitive type inhibitor while 1 and 15 showed mixed-type inhibition. Furthermore, molecular docking simulation displayed negative binding energies (AutoDock 4.2.6 = -10.0 to -7.1 kcal/mol) of 1, 13, and 15 for BACE1, indicating these compounds bound tightly to the active site of the enzyme with low energy and high affinity. The results of the present study clearly demonstrate that S. miltiorrhiza and its constituents have potential anti-AD activity and can be used as a therapeutic agent for the treatment of AD.

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β.

A PM7 dynamic residue-ligand interactions energy landscape of the BACE1 inhibitory pathway by hydroxyethylamine compounds. Part I: The flap closure process

BACE1 is an enzyme of scientific interest because it participates in the progression of Alzheimer's disease. Hydroxyethylamines (HEAs) are a family of compounds which exhibit inhibitory activity toward BACE1 at a nanomolar level, favorable pharmacokinetic properties and oral bioavailability. The first step in the inhibition of BACE1 by HEAs consists of their entrance into the protease active site and the resultant conformational change in the protein, from Apo to closed form. These two conformations differ in the position of an antiparallel loop (called the flap) which covers the entrance to the catalytic site. For BACE1, closure of this flap is vital to its catalytic activity and to inhibition of the enzyme due to the new interactions thereby formed with the ligand. In the present study a dynamic energy landscape of residue-ligand interaction energies (ReLIE) measured for 112 amino acids in the BACE1 active site and its immediate vicinity during the closure of the flap induced by 8 HEAs of different inhibitory power is presented. A total of 6.272 million ReLIE calculations, based on the PM7 semiempirical method, provided a deep and quantitative view of the first step in the inhibition of the aspartyl protease. The information suggests that residues Asp93, Asp289, Thr292, Thr293, Asn294 and Arg296 are anchor points for the ligand, accounting for approximately 45% of the total protein-ligand interaction. Additionally, flap closure improved the BACE1-HEA interaction by around 25%. Furthermore, the inhibitory activity of HEAs could be related to the capacity of these ligands to form said anchor point interactions and maintain them over time: the lack of some of these anchor interactions delayed flap closure or impeded it completely, or even caused the flap to reopen. The methodology employed here could be used as a tool to evaluate future structural modifications which lead to improvements in the favorability and stability of BACE1-HEA ReLIEs, aiding in the design of better inhibitors.

Fragment-based virtual screening approach and molecular dynamics simulation studies for identification of BACE1 inhibitor leads

Traditional structure-based virtual screening method to identify drug-like small molecules for BACE1 is so far unsuccessful. Location of BACE1, poor Blood Brain Barrier permeability and P-glycoprotein (Pgp) susceptibility of the inhibitors make it even more difficult. Fragment-based drug design method is suitable for efficient optimization of initial hit molecules for target like BACE1. We have developed a fragment-based virtual screening approach to identify/optimize the fragment molecules as a starting point. This method combines the shape, electrostatic, and pharmacophoric features of known fragment molecules, bound to protein conjugate crystal structure, and aims to identify both chemically and energetically feasible small fragment ligands that bind to BACE1 active site. The two top-ranked fragment hits were subjected for a 53 ns MD simulation. Principle component analysis and free energy landscape analysis reveal that the new ligands show the characteristic features of established BACE1 inhibitors. The potent method employed in this study may serve for the development of potential lead molecules for BACE1-directed Alzheimer's disease therapeutics.

Effects of Glycosylation on the Enzymatic Activity and Mechanisms of Proteases

Posttranslational modifications are an important feature of most proteases in higher organisms, such as the conversion of inactive zymogens into active proteases. To date, little information is available on the role of glycosylation and functional implications for secreted proteases. Besides a stabilizing effect and protection against proteolysis, several proteases show a significant influence of glycosylation on the catalytic activity. Glycans can alter the substrate recognition, the specificity and binding affinity, as well as the turnover rates. However, there is currently no known general pattern, since glycosylation can have both stimulating and inhibiting effects on activity. Thus, a comparative analysis of individual cases with sufficient enzyme kinetic and structural data is a first approach to describe mechanistic principles that govern the effects of glycosylation on the function of proteases. The understanding of glycan functions becomes highly significant in proteomic and glycomic studies, which demonstrated that cancer-associated proteases, such as kallikrein-related peptidase 3, exhibit strongly altered glycosylation patterns in pathological cases. Such findings can contribute to a variety of future biomedical applications.

A molecular dynamics study of the BACE1 conformational change from Apo to closed form induced by hydroxyethylamine derived compounds

BACE1 is an aspartyl protease which is a therapeutic target for Alzheimer's disease (AD) because of its participation in the rate-limiting step in the production of Aβ-peptide, the accumulation of which produces senile plaques and, in turn, the neurodegenerative effects associated with AD. The active site of this protease is composed in part by two aspartic residues (Asp93 and Asp289). Additionally, the catalytic site has been found to be covered by an antiparallel hairpin loop called the flap. The dynamics of this flap are fundamental to the catalytic function of the enzyme. When BACE1 is inactive (Apo), the flap adopts an open conformation, allowing a substrate or inhibitor to access the active site. Subsequent interaction with the ligand induces flap closure and the stabilization of the macromolecular complex. Further, the protonation state of the aspartic dyad is affected by the chemical nature of the species entering the active site, so that appropriate selection of protonation states for the ligand and the catalytic residues will permit the elucidation of the inhibitory pathway for BACE1. In the present study, comparative analysis of different combinations of protonation states for the BACE1-hydroxyethylamine (HEA) system is reported. HEAs are potent inhibitors of BACE1 with favorable pharmacological and kinetic properties, as well as oral bioavailability. The results of Molecular Dynamics (MD) simulations and population density calculations using 8 different parameters demonstrate that the LnAsp289 configuration (HEA with a neutral amine and the Asp289 residue protonated) is the only one which permits the expected conformational change in BACE1, from apo to closed form, after flap closure. Additionally, differences in their capacities to establish and maintain interactions with residues such as Asp93, Gly95, Thr133, Asp289, Gly291, and Asn294 during this step allow differentiation among the inhibitory activities of the HEAs. The results and methodology here reported will serve to elucidate the inhibitory pathway of other families of compounds that act as BACE1 inhibitors, as well as the design of better leader compounds for the treatment of AD.

Application of Free Energy Perturbation for the Design of BACE1 Inhibitors

Novel spiroaminodihydropyrroles probing for optimized interactions at the P3 pocket of β-secretase 1 (BACE1) were designed with the use of free energy perturbation (FEP) calculations. The resulting molecules showed pIC50 potencies in enzymatic BACE1 inhibition assays ranging from approximately 5 to 7. Good correlation was observed between the predicted activity from the FEP calculations and experimental activity. Simulations run with a default 5 ns approach delivered a mean unsigned error (MUE) between prediction and experiment of 0.58 and 0.91 kcal/mol for retrospective and prospective applications, respectively. With longer simulations of 10 and 20 ns, the MUE was in both cases 0.57 kcal/mol for the retrospective application, and 0.69 and 0.59 kcal/mol for the prospective application. Other considerations that impact the quality of the calculations are discussed. This work provides an example of the value of FEP as a computational tool for drug discovery.

Determination of the protonation state for the catalytic dyad in β-secretase when bound to hydroxyethylamine transition state analogue inhibitors: A molecular dynamics simulation study

BACE1 is an aspartyl protease of pharmacological interest for its direct participation in Alzheimer's disease (AD) through β-amyloid peptide production. Two aspartic acid residues are present in the BACE1 catalytic region which can adopt multiple protonation states depending on the chemical nature of its inhibitors, i.e., monoprotonated, diprotonated and di-deprotonated states. In the present study a series of protein-ligand molecular dynamics (MD) simulations was carried out to identify the most feasible protonation state adopted by the catalytic dyad in the presence of hydroxyethylamine transition state analogue inhibitors. The MD trajectories revealed that the di-deprotonated state is most prefered in the presence of hydroxyethilamine (HEA) family inhibitors. This appears as a result after evaluating, for all 9 protonation state configurations during the simulation time, the deviations of a set of distances and dihedral angles measured on the ligand, protein and protein-ligand complex with reference to an X-ray experimental BACE1/HEA crystallographic structure. These results will help to clarify the phenomena related to the HEAs inhibitory pathway, and improve HEAs databases' virtual screening and ligand design processes targeting β-secretase protein.

Specific Inhibition of β-Secretase Processing of the Alzheimer Disease Amyloid Precursor Protein

Development of disease-modifying therapeutics is urgently needed for treating Alzheimer disease (AD). AD is characterized by toxic β-amyloid (Aβ) peptides produced by β- and γ-secretase-mediated cleavage of the amyloid precursor protein (APP). β-secretase inhibitors reduce Aβ levels, but mechanism-based side effects arise because they also inhibit β-cleavage of non-amyloid substrates like Neuregulin. We report that β-secretase has a higher affinity for Neuregulin than it does for APP. Kinetic studies demonstrate that the affinities and catalytic efficiencies of β-secretase are higher toward non-amyloid substrates than toward APP. We show that non-amyloid substrates are processed by β-secretase in an endocytosis-independent manner. Exploiting this compartmentalization of substrates, we specifically target the endosomal β-secretase by an endosomally targeted β-secretase inhibitor, which blocked cleavage of APP but not non-amyloid substrates in many cell systems, including induced pluripotent stem cell (iPSC)-derived neurons. β-secretase inhibitors can be designed to specifically inhibit the Alzheimer process, enhancing their potential as AD therapeutics without undesired side effects.

Mechanisms of peptide hydrolysis by aspartyl and metalloproteases

Peptide hydrolysis has been involved in a wide range of biological, biotechnological, and industrial applications.

Computational Insights into Substrate and Site Specificities, Catalytic Mechanism, and Protonation States of the Catalytic Asp Dyad of β -Secretase

In this review, information regarding substrate and site specificities, catalytic mechanism, and protonation states of the catalytic Asp dyad of β-secretase (BACE1) derived from computational studies has been discussed. BACE1 catalyzes the rate-limiting step in the generation of Alzheimer amyloid beta peptide through the proteolytic cleavage of the amyloid precursor protein. Due to its biological functioning, this enzyme has been considered as one of the most important targets for finding the cure for Alzheimer's disease. Molecular dynamics (MD) simulations suggested that structural differences in the key regions (inserts A, D, and F and the 10s loop) of the enzyme are responsible for the observed difference in its activities towards the WT- and SW-substrates. The modifications in the flap, third strand, and insert F regions were found to be involved in the alteration in the site specificity of the glycosylphosphatidylinositol bound form of BACE1. Our QM and QM/MM calculations suggested that BACE1 hydrolyzed the SW-substrate more efficiently than the WT-substrate and that cleavage of the peptide bond occurred in the rate-determining step. The results from molecular docking studies showed that the information concerning a single protonation state of the Asp dyad is not enough to run an in silico screening campaign.

Cerebrospinal Aβ11-x and 17-x levels as indicators of mild cognitive impairment and patients' stratification in Alzheimer's disease

In the present work, the concentrations of Aβ11-x and Aβ17-x peptides (x=40 or 42), which result from the combined cleavages of β-amyloid precursor protein (AβPP) by β'/α or α/γ-secretases, respectively, were assessed in cerebrospinal fluid (CSF) samples from patients with Alzheimer's disease (AD) or mild cognitive impairment (MCI). Specific multiplexed assays were set up using new anti-40 and anti-42 monoclonal antibodies (mAbs) for the capture of these N-truncated Aβ peptides and anti-11 or anti-17 mAbs for their detection. The specificity, sensitivity and reproducibility of such assays were assessed using synthetic peptides and human cell models. Aβ11-x and Aβ17-x were then measured in CSF samples from patients with AD (n=23), MCI (n=23) and controls with normal cognition (n=21). Aβ11-x levels were significantly lower in patients with MCI than in controls. Compared with the combined quantification of Aβ1-42, total Tau (T-Tau) and phosphorylated Tau (P-Tau; AlzBio3, Innogenetics), the association of Aβ11-40, Aβ17-40 and T-Tau improved the discrimination between MCI and controls. Furthermore, when patients with MCI were classified into two subgroups (MCI ≤1.5 or ≥2 based on their CDR-SB (Cognitive Dementia Rating-Sum of Boxes) score), the CSF Aβ17-40/Aβ11-40 ratio was significantly higher in patients with CDR-SB ≤1.5 than in controls, whereas neither Aβ1-42, T-Tau nor P-Tau allowed the detection of this subpopulation. These results need to be confirmed in a larger clinical prospective cohort.

Molecular dynamics studies on both bound and unbound renin protease

The aspartic protease renin (REN) catalyses the rate-limiting step in the Renin-Angiotensin-Aldosterone System (RAAS), which regulates cardiovascular and renal homoeostasis in living organisms. Renin blockage is therefore an attractive therapeutic strategy for the treatment of hypertension. Herein, computational approaches were used to provide a structural characterization of the binding site, flap opening and dynamic rearrangements of REN in the key conserved residues and water molecules, with the binding of a dodecapeptide substrate or different inhibitors. All these structural insights during catalysis may assist future studies in developing novel strategies for REN inactivation. Our molecular dynamics simulations of several unbound-REN and bound-REN systems indicate similar flexible-segments plasticity with larger fluctuations in those belonging to the C-domain (exposed to the solvent). These segments are thought to assist the flap opening and closure to allow the binding of the substrate and catalytic water molecules. The unbound-REN simulation suggests that the flap can acquire three different conformations: closed, semi-open and open. Our results indicate that the semi-open conformation is already sufficient and appropriate for the binding of the angiotensinogen (Ang) tail, thus contributing to the high specificity of REN, and that both semi-open and open flap conformations are present in free and complexed enzymes. We additionally observed that the Tyr75-Trp39 H-bond has an important role in assisting flap movement, and we highlight several conserved water molecules and amino acids that are essential for the proper catalytic activity of REN.

Elucidating the catalytic mechanism of β-secretase (BACE1): a quantum mechanics/molecular mechanics (QM/MM) approach

In this quantum mechanics/molecular mechanics (QM/MM) study, the mechanisms of the hydrolytic cleavage of the Met2-Asp3 and Leu2-Asp3 peptide bonds of the amyloid precursor protein (WT-substrate) and its Swedish mutant (SW) respectively catalyzed by β-secretase (BACE1) have been investigated by explicitly including the electrostatic and steric effects of the protein environment in the calculations. BACE1 catalyzes the rate-determining step in the generation of Alzheimer amyloid beta peptides and is widely acknowledged as a promising therapeutic target. The general acid-base mechanism followed by the enzyme proceeds through the following two steps: (1) formation of the gem-diol intermediate and (2) cleavage of the peptide bond. The formation of the gem-diol intermediate occurs with the barriers of 19.6 and 16.1 kcal/mol for the WT- and SW-substrate respectively. The QM/MM energetics predict that with the barriers of 21.9 and 17.2 kcal/mol for the WT- and SW-substrate respectively the cleavage of the peptide bond occurs in the rate-determining step. The computed barriers are in excellent agreement with the measured barrier of ∼18.0 kcal/mol for the SW-substrate and in line with the experimental observation that the cleavage of this substrate is sixty times more efficient than the WT-substrate.

Computational Insights into Dynamics of Protein Aggregation and Enzyme-Substrate Interactions

In this Perspective, the roles of protein dynamics have been discussed in the aggregation of amyloid beta (Aβ) peptides and formation of enzyme-substrate complexes of beta-secretase (BACE1) and insulin-degrading enzyme (IDE). The studies regarding the influence of individual amino acid residues and specific regions on the structures and oligomerization of early Aβ aggregates and computations of their translational and rotational diffusion coefficients and order parameters exhibited that even the short-time-scale molecular dynamics simulations can reproduce certain experimental parameters with reasonable accuracy. The simulations elucidating the enzyme-substrate interactions of BACE1 and IDE successfully showed that the chemical nature and length of the substrates influence the dynamics and plasticity of both the enzyme and substrate. An atomic-level understanding of these processes will advance our efforts to develop therapeutic strategies for several deadly diseases through the design of small molecules with antiaggregation properties and substrate-specific "designer" forms of enzymes.

The catalytic mechanism of mouse renin studied with QM/MM calculations

Hypertension is a chronic condition that affects nearly 25% of adults worldwide. As the Renin-Angiotensin-Aldosterone System is implicated in the control of blood pressure and body fluid homeostasis, its combined blockage is an attractive therapeutic strategy currently in use for the treatment of several cardiovascular conditions. We have performed QM/MM calculations to study the mouse renin catalytic mechanism in atomistic detail, using the N-terminal His6-Asn14 segment of angiotensinogen as substrate. The enzymatic reaction (hydrolysis of the peptidic bond between residues in the 10th and 11th positions) occurs through a general acid/base mechanism and, surprisingly, it is characterized by three mechanistic steps: it begins with the creation of a first very stable tetrahedral gem-diol intermediate, followed by protonation of the peptidic bond nitrogen, giving rise to a second intermediate. In a final step the peptidic bond is completely cleaved and both gem-diol hydroxyl protons are transferred to the catalytic dyad (Asp32 and Asp215). The final reaction products are two separate peptides with carboxylic acid and amine extremities. The activation energy for the formation of the gem-diol intermediate was calculated as 23.68 kcal mol(-1), whereas for the other steps the values were 15.51 kcal mol(-1) and 14.40 kcal mol(-1), respectively. The rate limiting states were the reactants and the first transition state. The associated barrier (23.68 kcal mol(-1)) is close to the experimental values for the angiotensinogen substrate (19.6 kcal mol(-1)). We have also tested the influence of the density functional on the activation and reaction energies. All eight density functionals tested (B3LYP, B3LYP-D3, X3LYP, M06, B1B95, BMK, mPWB1K and B2PLYP) gave very similar results.

Protonation states of the catalytic dyad of β-secretase (BACE1) in the presence of chemically diverse inhibitors: a molecular docking study

In this molecular docking study, the protonation states of the catalytic Asp dyad of the beta-secretase (BACE1) enzyme in the presence of eight chemically diverse inhibitors have been predicted. BACE1 catalyzes the rate-determining step in the generation of Alzheimer amyloid beta peptides and is widely considered as a promising therapeutic target. All the inhibitors were redocked into their corresponding X-ray structures using a combination of eight different protonation states of the Asp dyad for each inhibitor. Five inhibitors were primarily found to favor two different monoprotonated states, and the remaining three favor a dideprotonated state. In addition, five of them exhibited secondary preference for a diprotonated state. These results show that the knowledge of a single protonation state of the Asp dyad is not sufficient to search for the novel inhibitors of BACE1 and the most plausible state for each inhibitor must be determined prior to conducting in-silico screening.

Loss of cleavage at β'-site contributes to apparent increase in β-amyloid peptide (Aβ) secretion by β-secretase (BACE1)-glycosylphosphatidylinositol (GPI) processing of amyloid precursor protein

Several lines of evidence implicate lipid raft microdomains in Alzheimer disease-associated β-amyloid peptide (Aβ) production. Notably, targeting β-secretase (β-site amyloid precursor protein (APP)-cleaving enzyme 1 (BACE1)) exclusively to lipid rafts by the addition of a glycosylphosphatidylinositol (GPI) anchor to its ectodomain has been reported to elevate Aβ secretion. Paradoxically, Aβ secretion is not reduced by the expression of non-raft resident S-palmitoylation-deficient BACE1 (BACE1-4C/A (C474A/C478A/C482A/C485A)). We addressed this apparent discrepancy in raft microdomain-associated BACE1 processing of APP in this study. As previously reported, we found that expression of BACE1-GPI elevated Aβ secretion as compared with wild-type BACE1 (WTBACE1) or BACE1-4C/A. However, this increase occurred without any difference in the levels of APP ectodomain released following BACE1 cleavage (soluble APPβ), arguing against an overall increase in BACE1 processing of APP per se. Further analysis revealed that WTBACE1 cleaves APP at β- and β'-sites, generating +1 and +11 β-C-terminal fragments and secreting intact as well as N-terminally truncated Aβ. In contrast, three different BACE1-GPI chimeras preferentially cleaved APP at the β-site, mainly generating +1 β-C-terminal fragment and secreting intact Aβ. As a consequence, cells expressing BACE1-GPI secreted relatively higher levels of intact Aβ without an increase in BACE1 processing of APP. Markedly reduced cleavage at β'-site exhibited by BACE1-GPI was cell type-independent and insensitive to subcellular localization of APP or the pathogenic KM/NL mutant. We conclude that the apparent elevation in Aβ secretion by BACE1-GPI is mainly attributed to preferential cleavage at the β-site and failure to detect +11 Aβ species secreted by cells expressing WTBACE1.