2-Amino-3,4-dihydroquinazolines as Inhibitors of BACE-1 (β-Site APP Cleaving Enzyme): Use of Structure Based Design to Convert a Micromolar Hit into a Nanomolar Lead

Rational design of novel diaryl ether-linked benzimidazole derivatives as potent and selective BACE1 inhibitors

Due to the large size and high flexibility of the catalytic active site of BACE1 enzyme, the development of nonpeptide inhibitors with optimal pharmacological properties is still highly demanding. In this work, we have discovered 2-aminobenzimidazole-containg ether scaffolds having potent and selective inhibitory potentials against BACE1 enzyme. We have synthesized novel 29 compounds and optimization of aryl linker region resulted in highly potent BACE1 inhibitory activities with EC50 values of 0.05-2.71 μM. The aryloxy-phenyl analogs 20j showed the EC50 value as low as 0.07 μM in the enzyme assay, whereas, the benzyloxyphenyl dervative 24b was comparatively less effective in the enzyme assay. But interestingly the latter was more effective in the cell assay (EC50 value 1.2 μM). While comparing synthesized derivatives in the cell assay using PC12-APPSW cell, compound 27f appeared as the most potent BACE1 inhibitor having EC50 value 0.7 μM. This scaffold also showed high selectivity over BACE2 enzyme and cathepsin D. Furthermore, the research findings were bolstered through the incorporation of molecular docking, molecular dynamics, and DFT studies. We firmly believe that these discoveries will pave the way for the development of a novel class of small-molecule selective BACE1 inhibitors.

To Design Scalable Free Energy Perturbation Networks, Optimal Is Not Enough

Drug discovery is accelerated with computational methods such as alchemical simulations to estimate ligand affinities. In particular, relative binding free energy (RBFE) simulations are beneficial for lead optimization. To use RBFE simulations to compare prospective ligands in silico, researchers first plan the simulation experiment, using graphs where nodes represent ligands and graph edges represent alchemical transformations between ligands. Recent work demonstrated that optimizing the statistical architecture of these perturbation graphs improves the accuracy of the predicted changes in the free energy of ligand binding. Therefore, to improve the success rate of computational drug discovery, we present the open-source software package High Information Mapper (HiMap)─a new take on its predecessor, Lead Optimization Mapper (LOMAP). HiMap removes heuristics decisions from design selection and instead finds statistically optimal graphs over ligands clustered with machine learning. Beyond optimal design generation, we present theoretical insights for designing alchemical perturbation maps. Some of these results include that for n number of nodes, the precision of perturbation maps is stable at n·ln(n) edges. This result indicates that even an "optimal" graph can result in unexpectedly high errors if a plan includes too few alchemical transformations for the given number of ligands and edges. And, as a study compares more ligands, the performance of even optimal graphs will deteriorate with linear scaling of the edge count. In this sense, ensuring an A- or D-optimal topology is not enough to produce robust errors. We additionally find that optimal designs will converge more rapidly than radial and LOMAP designs. Moreover, we derive bounds for how clustering reduces cost for designs with a constant expected relative error per cluster, invariant of the size of the design. These results inform how to best design perturbation maps for computational drug discovery and have broader implications for experimental design.

Synthesis, radiolabeling, and evaluation of a potent β-site APP cleaving enzyme (BACE1) inhibitor for PET imaging of BACE1 in vivo

The β-site APP-cleaving enzyme 1 (BACE1) plays important roles in the proteolytic processing of amyloid precursor protein, and can be regarded as an important target for the diagnosis and treatment of AD. This study aimed to report the synthesis and evaluation of an ¹⁸F-labeled 2-amino-3,4-dihydroquinazoline analog as a potential BACE1 radioligand. A fluoropropyl side chain was introduced to the phenyl of this 3,4-dihydroquinazoline scaffold to generate the radioligand. Our preliminary data indicated that although the 2-amino-3,4-dihydroquinazoline scaffold possessed favorable in-vitro properties as a PET ligand, its poor brain uptake hindered the in-vivo imaging of BACE1. Further investigation would be required to optimize the scaffold for the development of a blood-brain-barrier-permeable BACE1-targeted PET ligand.

Therapeutic potential of quinazoline derivatives for Alzheimer's disease: A comprehensive review

Quinazolines are considered as a promising class of bioactive heterocyclic compounds with broad properties. Particularly, the quinazoline scaffold has an impressive role in the design and synthesis of new CNS-active drugs. The drug-like properties and pharmacological characteristics of quinazoline could lead to different drugs with various targets. Among CNS disorders, Alzheimer's disease (AD) is a progressive neurodegenerative disorder with memory loss, cognitive decline and language dysfunction. AD is a complex and multifactorial disease therefore, the need for finding multi-target drugs against this devastative disease is urgent. A literature survey revealed that quinazoline derivatives have diverse therapeutic potential for AD as modulators/inhibitors of β-amyloid, tau protein, cholinesterases, monoamine oxidases, and phosphodiesterases as well as other protective effects. Thus, we describe here the most relevant and recent studies about anti-AD agents with quinazoline structure which can further aid the development and discovery of new anti-AD agents.

Synthesis and biological evaluation of benzimidazoles/1,3,5-triazine-2,4-diamine hybrid compounds: a new class of multifunctional alzheimer targeting agents

Twelve novel benzimidazole/1,3,5-triazine-2,4-diamine hybrids were synthesized and biologically studied as multifunctional Alzheimer-controlling agents.

Friedel-Crafts synthesis of bis(trifluoromethylated)-4-aryl-3,4-dihydroquinazolines, bis(trifluoromethylated)-3,4-dihydroquinazoline-4-ols and trifluoromethyl arylketoimines using N-aryltrifluoroacetimidoyl chlorides and benzene derivatives

A simple and straightforward approach to access biologically important N-heterocycles, namely bis(trifluoromethylated)-4-aryl-3,4-dihydroquinazolines, bis(trifluoromethylated)-3,4-dihydroquinazoline-4-ols and trifluoromethyl arylketoimines, as a group of important frameworks or initial substance from N-aryltrifluoroacetimidoyl chlorides with benzene derivatives via Friedel-Crafts reaction has been described. This novel strategy provides synthesis of trifluoromethylated dihydroquinazolines and trifluoromethyl arylketoimines in good to excellent yields.

Dynamic Docking Using Multicanonical Molecular Dynamics: Simulating Complex Formation at the Atomistic Level

Multicanonical molecular dynamics (McMD)-based dynamic docking has been applied to predict the native binding configurations for several protein receptors and their ligands. Due to the enhanced sampling capabilities of McMD, it can exhaustively sample bound and unbound ligand configurations, as well as receptor conformations, and thus enables efficient sampling of the conformational and configurational space, not possible using canonical MD simulations. As McMD samples a wide configurational space, extensive analysis is required to study the diverse ensemble consisting of bound and unbound structures. By projecting the reweighted ensemble onto the first two principal axes obtained via principal component analysis of the multicanonical ensemble, the free energy landscape (FEL) can be obtained. Further analysis produces representative structures positioned at the local minima of the FEL, where these structures are then ranked by their free energy. In this chapter, we describe our dynamic docking methodology, which has successfully reproduced the native binding configuration for small compounds, medium-sized compounds, and peptide molecules.

One Pot Synthesis of Micromolar BACE-1 Inhibitors Based on the Dihydropyrimidinone Scaffold and Their Thia and Imino Analogues

A library of dihydropyrimidinones was synthesized via a “one-pot” three component Biginelli reaction using different aldehydes in combination with β-dicarbonyl compounds and urea. Selected 2-thiooxo and 2-imino analogs were also obtained with the Biginelli reaction from thiourea and guanidine hydrochloride, respectively. The products were screened in vitro for their β-secretase inhibitory activity. The majority of the compounds resulted to be active, with IC₅₀ in the range 100 nM–50 μM.

4-Substituted 2-amino-3,4-dihydroquinazolines with a 3-hairpin turn side chain as novel inhibitors of BACE-1

Herein, we report the identification, design, and synthesis of a series of 4-substituted 2-amino-3,4-dihydroquinazolines with hairpin turn side chains as novel inhibitors of BACE-1. The dihydroquinazoline derivatives were rationally designed by modifying the amide group and relocating the α -hydrophobic substituent on the hairpin turn side chain of lead compound 2 to the C4-position on the 3,4-dihydroquinazoline scaffold to facilitate interactions with the S₁, S₂ and S₁' subsites of BACE-1. Among these derivatives, two compounds exhibited potent BACE-1 inhibitory activity: 4-methyl-substituted (22a, BACE-1 CFA IC₅₀ = 0.38 μM; BACE-1 WCA IC₅₀ = 0.14 μM) and 4-cyclohexylmethyl-substituted (22b, BACE-1 CFA IC₅₀ = 0.49 μM; BACE-1 WCA IC₅₀ = 0.14 μM) 2-amino-3,4-dihydroquinazoline, each bearing a side chain of N-cyclohexyl-N-((1-methyl-1H-pyrazol-4-yl)methyl amide. The results suggest that the structural modifications maintain the hairpin turn topology similar to that of compound 2 and provide an additional interaction with the S₂ subsite.

Triflic anhydride mediated synthesis of 3,4-dihydroquinazolines: a three-component one-pot tandem procedure

A one-pot three-component tandem reaction involving a key Pictet-Spengler-like annulation step has been developed, providing an efficient method for the synthesis of 3,4-dihydroquinazolines in moderate to good yields from amides, aldehydes, and amines. The multicomponent triflic anhydride mediated reaction tolerates the installation of numerous functional groups, affording extensive diversity about the heterocyclic scaffold.

Performance evaluation of molecular docking and free energy calculations protocols using the D3R Grand Challenge 4 dataset

Using the D3R Grand Challenge 4 dataset containing Beta-secretase 1 (BACE) and Cathepsin S (CatS) inhibitors, we have evaluated the performance of our in-house docking workflow that involves in the first step the selection of the most suitable docking software for the system of interest based on structural and functional information available in public databases, followed by the docking of the dataset to predict the binding modes and ranking of ligands. The macrocyclic nature of the BACE ligands brought additional challenges, which were dealt with by a careful preparation of the three-dimensional input structures for ligands. This provided top-performing predictions for BACE, in contrast with CatS, where the predictions in the absence of guiding constraints provided poor results. These results highlight the importance of previous structural knowledge that is needed for correct predictions on some challenging targets. After the end of the challenge, we also carried out free energy calculations (i.e. in a non-blinded manner) for CatS using the pmx software and several force fields (AMBER, Charmm). Using knowledge-based starting pose construction allowed reaching remarkable accuracy for the CatS free energy estimates. Interestingly, we show that the use of a consensus result, by averaging the results from different force fields, increases the prediction accuracy.

Evaluation of a Series of β-Secretase 1 Inhibitors Containing Novel Heteroaryl-Fused-Piperazine Amidine Warheads

Despite several years of research, only a handful of β-secretase (BACE) 1 inhibitors have entered clinical trials as potential therapeutics against Alzheimer's disease. The intrinsic basic nature of low molecular weight, amidine-containing BACE 1 inhibitors makes them far from optimal as central nervous system drugs. Herein we present a set of novel heteroaryl-fused piperazine amidine inhibitors designed to lower the basicity of the key, enzyme binding, amidine functionality. This study resulted in the identification of highly potent (IC₅₀ ≤ 10 nM), permeable lead compounds with a reduced propensity to suffer from P-glycoprotein-mediated efflux.

Protein aggregation linked to Alzheimer's disease revealed by saturation transfer MRI

The goal of this study was to develop a molecular biomarker for the detection of protein aggregation involved in Alzheimer's disease (AD) by exploiting the features of the water saturation transfer spectrum (Z-spectrum), the CEST signal of which is sensitive to the molecular configuration of proteins. A radial-sampling steady-state sequence based ultrashort echo time (UTE) readout was implemented to image the Z-spectrum in the mouse brain, especially the contributions from mobile proteins at the frequency offsets for the composite protein amide proton (+3.6 ppm) and aliphatic proton (-3.6 ppm) signals. Using a relatively weak radiofrequency (RF) saturation amplitude, contributions due to strong magnetization transfer contrast (MTC) from solid-like macromolecules and direct water saturation (DS) were minimized. For practical measure of the changes in the mobile protein configuration, we defined a saturation transfer difference (ΔST) by subtracting the Z-spectral signals at ±3.6 ppm from a control signal at 8 ppm. Phantom studies of glutamate solution, protein (egg white) and hair conditioner show the capability of the proposed scheme to minimize the contributions from amine protons, DS, and MTC, respectively. The ST signal at ±3.6 ppm of the cross-linked bovine serum albumin (BSA) solutions demonstrated that the ΔST signal can be used to monitor the aggregation process of the mobile proteins. High-resolution ΔST images of AD mouse brains at ±3.6 ppm of mouse brains showed significantly reduced ΔST (-3.6) signal compared to the age-matched wild-type (WT) mice. Thus, this signal has potential to serve as a molecular biomarker for monitoring protein aggregation in AD.

Bennett acceptance ratio method to calculate the binding free energy of BACE1 inhibitors: Theoretical model and design of new ligands of the enzyme

In recent years, the design, development, and evaluation of several inhibitors of the BACE1 enzyme, as part of Alzheimer's treatment, have gathered the scientific community's interest. Here, a linear regression model was built using binding free energy calculations through the Bennett acceptance ratio method for 20 known inhibitors of the BACE1 enzyme, with a Pearson coefficient of R = 0.88 and R²  = 0.78. The validation of this model was verified employing eight additional random inhibitors, which also gave a linear correlation with R = 0.97 and R²  = 0.93. Furthermore, this linear regression model was also used for proposing the structure of four potential BACE1 inhibitors, and the most active of them gave a theoretical K_d  = 10 nM. However, these molecules have not been synthesized yet. Our team used a total time of more than 800 ns for the Molecular Dynamics to carry out this study, and all the software used were freely available.

Synthesis of dihydroquinazolines from 2-aminobenzylamine: N 3 -aryl derivatives with electron-withdrawing groups

The sequential N-functionalization of 2-aminobenzylamine (2-ABA) followed by cyclodehydration allowed for a straightforward and efficient synthesis of 3,4-dihydroquinazolines with N-aryl substituents bearing electron-withdrawing groups. The sequence involves an initial S_NAr displacement, N-acylation and MW-assisted ring closure. Remarkably, the uncatalyzed N-arylation of 2-ABA led to the monosubstitution product using equimolar amounts of both reagents. The individual steps were optimized achieving good to excellent overall yields of the desired heterocycles, avoiding additional protection and deprotection steps. A mechanistic interpretation for the cyclodehydration reaction promoted by trimethylsilyl polyphosphate (PPSE) is also proposed on the basis of literature data and our experimental observations.

Reduced Aβ secretion by human neurons under conditions of strongly increased BACE activity

The initial step in the amyloidogenic cascade of amyloid precursor protein (APP) processing is catalyzed by beta-site APP-cleaving enzyme (BACE), and this protease has increased activities in affected areas of Alzheimer's disease brains. We hypothesized that altered APP processing, because of augmented BACE activity, would affect the actions of direct and indirect BACE inhibitors. We therefore compared post-mitotic human neurons (LUHMES) with their BACE-overexpressing counterparts (BLUHMES). Although β-cleavage of APP was strongly increased in BLUHMES, they produced less full-length and truncated amyloid beta (Aβ) than LUHMES. Moreover, low concentrations of BACE inhibitors decreased cellular BACE activity as expected, but increased Aβ_1-40 levels. Several other approaches to modulate BACE activity led to a similar, apparently paradoxical, behavior. For instance, reduction in intracellular acidification by bepridil increased Aβ production in parallel with decreased BACE activity. In contrast to BLUHMES, the respective control cells (LUHMES or BLUHMES with catalytically inactive BACE) showed conventional pharmacological responses. Other non-canonical neurochemical responses (so-called 'rebound effects') are well-documented for the Aβ pathway, especially for γ-secretase: a partial block of its activity leads to an increased Aβ secretion by some cell types. We therefore compared LUHMES and BLUHMES regarding rebound effects of γ-secretase inhibitors and found an Aβ rise in LUHMES but not in BLUHMES. Thus, different cellular factors are responsible for the γ-secretase- versus BACE-related Aβ rebound. We conclude that increased BACE activity, possibly accompanied by an altered cellular localization pattern, can dramatically influence Aβ generation in human neurons and affect pharmacological responses to secretase inhibitors. OPEN PRACTICES: Open Science: This manuscript was awarded with the Open Materials Badge. For more information see: https://cos.io/our-services/open-science-badges/.

High throughput evaluation of macrocyclization strategies for conformer stabilization

While macrocyclization of a linear compound to stabilize a known bioactive conformation can be a useful strategy to increase binding potency, the difficulty of macrocycle synthesis can limit the throughput of such strategies. Thus computational techniques may offer the higher throughput required to screen large numbers of compounds. Here we introduce a method for evaluating the propensity of a macrocyclic compound to adopt a conformation similar that of a known active linear compound in the binding site. This method can be used as a fast screening tool for prioritizing macrocycles by leveraging the assumption that the propensity for the known bioactive substructural conformation relates to the affinity. While this method cannot to identify new interactions not present in the known linear compound, it could quickly differentiate compounds where the three dimensional geometries imposed by the macrocyclization prevent adoption of conformations with the same contacts as the linear compound in their conserved region. Here we report the implementation of this method using an RMSD-based structural descriptor and a Boltzmann-weighted propensity calculation and apply it retrospectively to three macrocycle linker optimization design projects. We found the method performs well in terms of prioritizing more potent compounds.

Recent advances in transition-metal-catalyzed reactions of alkynes with isoxazoles

Isoxazoles, as masked 1,3-dicarbonyl equivalents, have proven to be versatile building blocks and pivotal intermediates for the construction of a variety of useful azacycles with molecular complexity. As a result, a range of new reactions have been discovered based on isoxazoles in the past decade. However, the relevant reactions of isoxazoles with alkynes have seldom been explored. In this review, we will focus on the recent progress in the transition-metal-catalyzed formal annulations for the efficient synthesis of N-heterocycles between alkynes and isoxazoles by highlighting their specificity and applicability, and the mechanistic rationale is presented where possible.

Palladium-catalyzed selective synthesis of 3,4-dihydroquinazolines from electron-rich arylamines, electron-poor arylamines and glyoxalates

Palladium-catalyzed selective synthesis of various 3,4-dihydroquinazolines from electron-rich arylamines, electron-poor arylamines and glyoxalates was developed.

Accurate and Reliable Prediction of the Binding Affinities of Macrocycles to Their Protein Targets

Macrocycles have been emerging as a very important drug class in the past few decades largely due to their expanded chemical diversity benefiting from advances in synthetic methods. Macrocyclization has been recognized as an effective way to restrict the conformational space of acyclic small molecule inhibitors with the hope of improving potency, selectivity, and metabolic stability. Because of their relatively larger size as compared to typical small molecule drugs and the complexity of the structures, efficient sampling of the accessible macrocycle conformational space and accurate prediction of their binding affinities to their target protein receptors poses a great challenge of central importance in computational macrocycle drug design. In this article, we present a novel method for relative binding free energy calculations between macrocycles with different ring sizes and between the macrocycles and their corresponding acyclic counterparts. We have applied the method to seven pharmaceutically interesting data sets taken from recent drug discovery projects including 33 macrocyclic ligands covering a diverse chemical space. The predicted binding free energies are in good agreement with experimental data with an overall root-mean-square error (RMSE) of 0.94 kcal/mol. This is to our knowledge the first time where the free energy of the macrocyclization of linear molecules has been directly calculated with rigorous physics-based free energy calculation methods, and we anticipate the outstanding accuracy demonstrated here across a broad range of target classes may have significant implications for macrocycle drug discovery.

Mitophagy Failure in Fibroblasts and iPSC-Derived Neurons of Alzheimer's Disease-Associated Presenilin 1 Mutation

Familial Alzheimer’s disease (FAD) is clearly related with the accumulation of amyloid-beta (Aβ) and its deleterious effect on mitochondrial function is well established. Anomalies in autophagy have also been described in these patients. In the present work, functional analyses have been performed to study mitochondrial recycling process in patient-derived fibroblasts and neurons from induced pluripotent stem cells harboring the presenilin 1 mutation A246E. Mitophagy impairment was observed due to a diminished autophagy degradation phase associated with lysosomal anomalies, thus causing the accumulation of dysfunctional mitochondria labeled by Parkin RBR E3 ubiquitin protein ligase (PARK2). The failure of mitochondrial recycling by autophagy was enhanced in the patient-derived neuronal model. Our previous studies have demonstrated similar mitophagy impairment in sporadic Alzheimer’s disease (AD); therefore, our data indicate that mitophagy deficiency should be considered a common nexus between familial and sporadic cases of the disease.

Syntheses of 3,4- and 1,4-dihydroquinazolines from 2-aminobenzylamine

A straightforward strategy for the synthesis of dihydroquinazolines is presented, which allows for the preparation of 3,4- and 1,4-dihydroquinazolines with different substitution patterns from 2-aminobenzylamine (2-ABA) as common precursor. The required functionalization of both amino groups present in 2-ABA was achieved by different routes involving selective N-acylation and cesium carbonate-mediated N-alkylation reactions, avoiding protection/deprotection steps. The heterocycles were efficiently synthesized in short reaction times by microwave-assisted ring closure of the corresponding aminoamides promoted by ethyl polyphosphate (PPE).

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

Discovery of imidazopyridines containing isoindoline-1,3-dione framework as a new class of BACE1 inhibitors: Design, synthesis and SAR analysis

Alzheimer's disease is characterized by chronic neurodegeneration leading to dementia. The main cause of neurodegeneration is considered to be the accumulation of amyloid-β. Inhibiting BACE1 is a well-studied approach to lower the burden of amyloid-β aggregates. We designed a series of imidazopyridines-based compounds bearing phthalimide moieties as inhibitors of BACE1. The compounds 8a-o were synthesized by the Groebke-Blackburn-Bienaymé three-component reaction of heteroaromatic amidines, aldehydes and isocyanides. Evaluating the BACE1 inhibitory effects of the synthesized compounds revealed that introducing an aminocyclohexyl moiety in the imidazopyridine core resulted in a significant improvement in its BACE1 inhibitory potential. In this regard, compound 8e was the most potent against BACE1 with an IC₅₀ value of 2.84 (±0.95) μM. Molecular docking revealed that the nitrogen atom of imidazopyridines and the oxygen atom of the phenoxypropyl linker were involved in hydrogen bound interactions with Asp228 and Asp32 of BACE1 active site, respectively. The phthalimide moiety oriented toward the flap pocket and interacted with phe108, lle110, Trp115, Ile118 through van der Waal's and hydrophobic interactions. These findings demonstrate that imidazopyridines-based compounds bearing phthalimide moiety have the potential to decrease amyloid-β levels and ameliorate the symptoms of Alzheimer's disease.

Effect of sampling on BACE-1 ligands binding free energy predictions via MM-PBSA calculations

The BACE-1 enzyme is a prime target to find a cure to Alzheimer's disease. In this article, we used the MM-PBSA approach to compute the binding free energies of 46 reported ligands to this enzyme. After showing that the most probable protonation state of the catalytic dyad is mono-protonated (on ASP32), we performed a thorough analysis of the parameters influencing the sampling of the conformational space (in total, more than 35 μs of simulations were performed). We show that ten simulations of 2 ns gives better results than one of 50 ns. We also investigated the influence of the protein force field, the water model, the periodic boundary conditions artifacts (box size), as well as the ionic strength. Amber03 with TIP3P, a minimal distance of 1.0 nm between the protein and the box edges and a ionic strength of I = 0.2 M provides the optimal correlation with experiments. Overall, when using these parameters, a Pearson correlation coefficient of R = 0.84 (R²  = 0.71) is obtained for the 46 ligands, spanning eight orders of magnitude of K_d (from 0.017 nm to 2000 μM, i.e., from -14.7 to -3.7 kcal/mol), with a ligand size from 22 to 136 atoms (from 138 to 937 g/mol). After a two-parameter fit of the binding affinities for 12 of the ligands, an error of RMSD = 1.7 kcal/mol was obtained for the remaining ligands. © 2017 Wiley Periodicals, Inc.

New directions for protease inhibitors directed drug discovery

Proteases play crucial roles in various biological processes, and their activities are essential for all living organisms-from viruses to humans. Since their functions are closely associated with many pathogenic mechanisms, their inhibitors or activators are important molecular targets for developing treatments for various diseases. Here, we describe drugs/drug candidates that target proteases, such as malarial plasmepsins, β-secretase, virus proteases, and dipeptidyl peptidase-4. Previously, we reported inhibitors of aspartic proteases, such as renin, human immunodeficiency virus type 1 protease, human T-lymphotropic virus type I protease, plasmepsins, and β-secretase, as drug candidates for hypertension, adult T-cell leukaemia, human T-lymphotropic virus type I-associated myelopathy, malaria, and Alzheimer's disease. Our inhibitors are also described in this review article as examples of drugs that target proteases. © 2015 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 106: 563-579, 2016.

Active Site Mapping of an Aspartic Protease by Multiple Fragment Crystal Structures: Versatile Warheads To Address a Catalytic Dyad

Crystallography is frequently used as follow-up method to validate hits identified by biophysical screening cascades. The capacity of crystallography to directly screen fragment libraries is often underestimated, due to its supposed low-throughput and need for high-quality crystals. We applied crystallographic fragment screening to map the protein-binding site of the aspartic protease endothiapepsin by individual soaking experiments. Here, we report on 41 fragments binding to the catalytic dyad and adjacent specificity pockets. The analysis identifies already known warheads but also reveals hydrazide, pyrazole, or carboxylic acid fragments as novel functional groups binding to the dyad. A remarkable swapping of the S1 and S1' pocket between structurally related fragments is explained by either steric demand, required displacement of a well-bound water molecule, or changes of trigonal-planar to tetrahedral geometry of an oxygen functional group in a side chain. Some warheads simultaneously occupying both S1 and S1' are promising starting points for fragment-growing strategies.

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.

Fragment-Based Discovery of 2-Aminoquinazolin-4(3H)-ones As Novel Class Nonpeptidomimetic Inhibitors of the Plasmepsins I, II, and IV

2-Aminoquinazolin-4(3H)-ones were identified as a novel class of malaria digestive vacuole plasmepsin inhibitors by using NMR-based fragment screening against Plm II. Initial fragment hit optimization led to a submicromolar inhibitor, which was cocrystallized with Plm II to produce an X-ray structure of the complex. The structure showed that 2-aminoquinazolin-4(3H)-ones bind to the open flap conformation of the enzyme and provided clues to target the flap pocket. Further improvement in potency was achieved via introduction of hydrophobic substituents occupying the flap pocket. Most of the 2-aminoquinazolin-4(3H)-one based inhibitors show a similar activity against digestive Plms I, II, and IV and >10-fold selectivity versus CatD, although varying the flap pocket substituent led to one Plm IV selective inhibitor. In cell-based assays, the compounds show growth inhibition of Plasmodium falciparum 3D7 with IC50 ∼ 1 μM. Together, these results suggest 2-aminoquinazolin-4(3H)-ones as perspective leads for future development of an antimalarial agent.

BACE1 (β-secretase) inhibitors for the treatment of Alzheimer's disease

BACE1 (β-secretase, memapsin 2, Asp2) has emerged as a promising target for the treatment of Alzheimer's disease. BACE1 is an aspartic protease which functions in the first step of the pathway leading to the production and deposition of amyloid-β peptide (Aβ). Its gene deletion showed only mild phenotypes. BACE1 inhibition has direct implications in the Alzheimer's disease pathology without largely affecting viability. However, inhibiting BACE1 selectively in vivo has presented many challenges to medicinal chemists. Since its identification in 2000, inhibitors covering many different structural classes have been designed and developed. These inhibitors can be largely classified as either peptidomimetic or non-peptidic inhibitors. Progress in these fields resulted in inhibitors that contain many targeted drug-like characteristics. In this review, we describe structure-based design strategies and evolution of a wide range of BACE1 inhibitors including compounds that have been shown to reduce brain Aβ, rescue the cognitive decline in transgenic AD mice and inhibitor drug candidates that are currently in clinical trials.