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

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.

Molecular Multi-Target Approach for Human Acetylcholinesterase, Butyrylcholinesterase and β-Secretase 1: Next Generation for Alzheimer's Disease Treatment

Alzheimer's Disease (AD) is a neurodegenerative condition characterized by progressive memory loss and other affected cognitive functions. Pharmacological therapy of AD relies on inhibitors of the enzymes acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), offering only a palliative effect and being incapable of stopping or reversing the neurodegenerative process. However, recent studies have shown that inhibiting the enzyme β-secretase 1 (BACE-1) may be able to stop neurodegeneration, making it a promising target. Considering these three enzymatic targets, it becomes feasible to apply computational techniques to guide the identification and planning of molecules capable of binding to all of them. After virtually screening 2119 molecules from a library, 13 hybrids were built and further screened by triple pharmacophoric model, molecular docking, and molecular dynamics (t = 200 ns). The selected hybrid G meets all stereo-electronic requirements to bind to AChE, BChE, and BACE-1 and offers a promising structure for future synthesis, enzymatic testing, and validation.

Catalytic Mechanism of Human T-Cell Leukemia Virus Type 1 Protease Investigated by Combined QM/MM Molecular Dynamics Simulations

Combined quantum mechanical and molecular mechanical (QM/MM) molecular dynamics simulations were performed to investigate the catalytic mechanism of human T-cell leukemia virus type 1 (HTLV-1) protease, a retroviral aspartic protease that is a potential therapeutic target for curing HTLV-1-associated diseases. To elucidate the proteolytic cleavage mechanism, we determined the two-dimensional free energy surfaces of the HTLV-1 protease-catalyzed reactions through various possible pathways. The free energy simulations suggest that the catalytic reactions of the HTLV-1 protease occur in the following sequential steps: (1) a proton is transferred from the lytic water to Asp32', followed by the nucleophilic addition of the resulting hydroxyl to the carbonyl carbon of the scissile bond, forming a tetrahedral oxyanion intermediate, and (2) a proton is transferred from Asp32 to the peptide nitrogen of the scissile bond, leading to the spontaneous breakage of the scissile bond. The rate-limiting step of this catalytic process is the proton transfer from Asp32 to the peptide nitrogen of the scissile bond, with a free energy of activation of 21.1 kcal/mol. This free energy barrier is close to the experimentally determined free energy of activation (16.3 kcal/mol) calculated from the measured catalytic rate constant (k_cat). This mechanistic study provides detailed dynamic and structural information that will facilitate the design of mechanism-based inhibitors for the treatment of HTLV-1-associated diseases.

Probing the intermolecular interactions, binding affinity, charge density distribution and dynamics of silibinin in dual targets AChE and BACE1: QTAIM and molecular dynamics perspective

Alzheimer's disease (AD) is the grievous neurodegenerative disorder. Reportedly, many enzymes are responsible for this disease, in which notably, acetylcholinesterase (AChE) and β-secretase (BACE1) are largely involved for AD. An experimental study reports that silibinin molecule inhibits both AChE and BACE1 enzymes. Present study aims to understand the dual binding mechanism of silibinin in the active site of AChE and BACE1 from the intermolecular interactions, conformational flexibility, charge density distribution, binding energy and the stability of molecule. To obtain the above information, the molecular docking, molecular dynamics (MD) and QTAIM (quantum theory of atoms in molecules) calculations have been performed. The molecular docking reveals that silibinin molecule is forming strong and weak intermolecular interactions with the catalytic site of both enzymes. The QTAIM analysis for the binding pockets of both complexes shows the charge density distribution of intermolecular interactions. The electrostatic potential map displays the electronegative/positive regions at the interaction zone of silibinin with AChE and BACE1 complexes. The MD simulation confirms that the silibinin molecule is stable in the active site of AChE and BACE1 enzymes. The binding free energies of silibinin with both enzymes are more favorable to have the interactions.Communicated by Ramaswamy H. Sarma.

New insights into the catalytic mechanism of the SARS-CoV-2 main protease: an ONIOM QM/MM approach

SARS-CoV-2 M^pro, also known as the main protease or 3C-like protease, is a key enzyme involved in the replication process of the virus that is causing the COVID-19 pandemic. It is also the most promising antiviral drug target targeting SARS-CoV-2 virus. In this work, the catalytic mechanism of M^pro was studied using the full model of the enzyme and a computational QM/MM methodology with a 69/72-atoms QM region treated at DLPNO-CCSD(T)/CBS//B3LYP/6-31G(d,p):AMBER level and including the catalytic important oxyanion-hole residues. The transition state of each step was fully characterized and described together with the related reactants and products. The rate-limiting step of the catalytic process is the hydrolysis of the thioester-enzyme adduct, and the calculated barrier closely agrees with the available kinetic data. The calculated Gibbs free energy profile, together with the full atomistic detail of the structures involved in catalysis, can now serve as valuable models for the rational drug design of transition state analogs as new inhibitors targeting the SARS-CoV-2 virus.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

Computational modelling of potent β-secretase (BACE1) inhibitors towards Alzheimer's disease treatment

Researchers have identified the β-amyloid precursor protein cleaving enzyme 1 (BACE1) in the multifactorial pathway of Alzheimer's disease (AD) as a drug target. The design and development of molecules to inhibit BACE1 as a potential cure for AD thus remained significant. Herein, we simulated two potent BACE1 inhibitors (AM-6494 and CNP-520) to understand their binding affinity at the atomistic level. AM-6494 is a newly reported potent BACE1 inhibitor with an IC₅₀ value of 0.4 nM in vivo and now picked for preclinical considerations. Umibecestat (CNP-520), which was discontinued at human trials lately, was considered to enable a reasonable evaluation of our results. Using density functional theory (DFT) and Our Own N-layered Integrated molecular Orbital and Molecular Mechanics (ONIOM), we achieved the aim of this investigation. These computational approaches enabled the prediction of the electronic properties of AM-6494 and CNP-520 plus their binding energies when complexed with BACE1. For AM-6494 and CNP-520 interaction with protonated BACE1, the ONIOM calculation gave binding free energy of -62.849 and -33.463 kcal/mol, respectively. In the unprotonated model, we observed binding free energy of -59.758 kcal/mol in AM-6494. Taken together thermochemistry of the process and molecular interaction plot, AM-6494 is more favourable than CNP-520 towards the inhibition of BACE1. The protonated model gave slightly better binding energy than the unprotonated form. However, both models could sufficiently describe ligand binding to BACE1 at the atomistic level. Understanding the detailed molecular interaction of these inhibitors could serve as a basis for pharmacophore exploration towards improved inhibitor design.

Synthesis and biological evaluation of 2'-Aminochalcone: A multi-target approach to find drug candidates to treat Alzheimer's disease

Alzheimer's disease (AD) is a neurodegenerative process that compromises cognitive functions. The physiopathology of AD is multifactorial and is mainly supported by the cholinergic and amyloid hypotheses, which allows the identification the fundamental role of some markers, such as the enzymes acetylcholinesterase (AChE) and β-secretase (BACE-1), and the β-amyloid peptide (Aβ). In this work, we prepared a series of chalcones and 2'-aminochalcones, which were tested against AChE and BACE-1 enzymes and on the aggregation of Aβ. All compounds inhibited AChE activity with different potencies. We have found that the majority of chalcones having the amino group are able to inhibit BACE-1, which was not observed for chalcones without this group. The most active compound is the one derived from 2,3-dichlorobenzaldeyde, having an IC₅₀ value of 2.71 μM. A molecular docking study supported this result, showing a good interaction of the amino group with aspartic acid residues of the catalytic diade of BACE-1. Thioflavin-T fluorescence emission is reduced in 30 - 40%, when Aβ₄₂ is incubated in the presence of some chalcones under aggregation conditions. In vitro cytotoxicity and in silico prediction of pharmacokinetic properties were also conducted in this study.

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.

Structure Related Inhibition of Enzyme Systems in Cholinesterases and BACE1 In Vitro by Naturally Occurring Naphthopyrone and Its Glycosides Isolated from Cassia obtusifolia

Cassia obtusifolia Linn. have been used to improve vision, inflammatory diseases, and as hepatoprotective agents and to promote urination from ancient times. In the present study, we investigated the influence of glycosylation of components of C. obtusifolia and structure-activity relationships (SARs) with respect to the inhibition of acetylcholinesterase (AChE), butyrylcholinesterase (BChE), and β-site amyloid precursor protein (APP)-cleaving enzyme 1 (BACE1), which are related to Alzheimer’s disease (AD). All six C. obtusifolia-derived compounds, rubrofusarin (1), rubrofusarin 6-O-β-d-glucopyranoside (2), rubrofusarin 6-O-β-d-gentiobioside (3), nor-rubrofusarin 6-O-β-d-glucoside (4), isorubrofusarin 10-O-β-d-gentiobioside (5), and rubrofusarin 6-O-β-d-triglucoside (6) showed promising inhibitory activity against AChE/BACE1. Compounds 3 and 4 showed most significant inhibition against AChE and BACE1, respectively. The SARs results emphasized the importance of gentiobiosyl moiety in the rubrofusarin for AChE inhibition, whereas the presence of hydroxyl group at C-8 and the glucosyl moiety at the C-6 position in the nor-rubrofusarin appeared to largely determine BACE1 inhibition. Kinetics and docking studies showed the lowest binding energy and highest affinity for mixed-type inhibitors, 3 and 4. Hydrophobic bonds interactions and the number of hydrogen bonds determined the strength of the protein-inhibitor interaction. These results suggest that C. obtusifolia and its constituents have therapeutic potential, and that the SARs of its active components are further explored with a view towards developing a treatment for AD.

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.

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.

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.