Fredy Sussman

Effect of the protonation state of the titratable residues on the inhibitor affinity to BACE-1.

BACE-1 is one of the aspartic proteases involved in the cleavage of beta amyloid peptide, an initial step in the formation of amyloid plaques whose toxicity induces neuron death in Alzheimers disease patients. One of the central issues in the search of novel BACE-1 inhibitors is the optimum pH for the binding of inhibitors to the enzyme. It is known that the enzyme has optimal catalytic activity at acidic pH, while cell active inhibitors may bind optimally at higher pH. In this work we determine the effect of the pH on the affinities of a set of inhibitors, with a variety of chemical motifs, for the ectodomain region of BACE-1 by a surface plasmon resonance (SPR) biosensor based assay. In order to understand the molecular interactions that underlie the diverse optimum pH for the binding of the various inhibitors as observed experimentally, we have calculated the titration curves for a set of BACE-1 ligand complexes. The results indicate that the pK(a) values of the titratable residues of the protein depend on the nature of the ligand involved, in disagreement with previous work. The enzyme-inhibitor structures with the resulting protonation states at pH values 4.5 and 7.4 served as the starting point for the prediction of the pH-dependent binding ranking. Our calculations reproduced the entire affinity ranking observed upon pH increase and most of the binding trends among inhibitors, especially at low pH. Finally, our cell-based assays indicate a possible correlation between high inhibitor affinity at both acidic and neutral pH values, with optimal cell response, a result that may open new venues for the search of potent BACE-1 inhibitors that are active at the cellular level.

Computer-Aided Structure-Based Design of Multitarget Leads for Alzheimer’s Disease

Alzheimers disease is a neurodegenerative pathology with unmet clinical needs. A highly desirable approach to this syndrome would be to find a single lead that could bind to some or all of the selected biomolecules that participate in the amyloid cascade, the most accepted route for Alzheimer disease genesis. In order to circumvent the challenge posed by the sizable differences in the binding sites of the molecular targets, we propose a computer-assisted protocol based on a pharmacophore and a set of required interactions with the targets that allows for the automated screening of candidates. We used a combination of docking and molecular dynamics protocols in order to discard nonbinders, optimize the best candidates, and provide a rationale for their potential as inhibitors. To provide a proof of concept, we proceeded to screen the literature and databases, a task that allowed us to identify a set of carbazole-containing compounds that initially showed affinity only for the cholinergic targets in our experimental assays. Two cycles of design based on our protocol led to a new set of analogues that were synthesized and assayed. The assay results revealed that the designed inhibitors had improved affinities for BACE-1 by more than 3 orders of magnitude and also displayed amyloid aggregation inhibition and affinity for AChE and BuChE, a result that led us to a group of multitarget amyloid cascade inhibitors that also could have a positive effect at the cholinergic level.

On the active site protonation state in aspartic proteases: implications for drug design.

Aspartic proteases (AP) are a family of important hydrolytic enzymes in medicinal chemistry, since many of its members have become therapeutical targets for a wide variety of diseases from AIDS to Alzheimer. The enzymatic activity of these proteins is driven by the Asp dyad, a pair of active site Asp residues that participate in the hydrolysis of peptides. Hence, the protonation state of these and other acidic residues present in these enzymes determines the catalytic rate and the affinity for an inhibitor at a given pH. In the present work we have reviewed the effect of the protonation states of the titratable residues in APs both on catalysis and inhibition in this family of enzymes. The first section focuses on the details of the catalytic reaction mechanism picture brought about by a large number of kinetic, crystallographic and computational chemistry analyses. The results indicate that although the mechanism is similar in both retroviral and eukaryotic enzymes, there are some clear differences. For instance, while in the former family branch the binding of the substrate induces a mono-ionic charge state for the Asp dyad, this charge state seems to be already present in the unbound state of the eukaryotic enzymes. In this section we have explored as well the possible existence of low barrier hydrogen bonds (LBHBs) in the enzymatic path. Catalytic rate enhancement in APs could in part be explained by the lowering of the barrier for proton transfer in a hydrogen bond from donor to acceptor, which is a typical feature of LBHBs. Review of the published work indicates that the experimental support for this type of bonds is rather scarce and it may be more probable in the first stages of the hydrolytic mechanism in retroviral proteases. The second section deals with the effect of active site protonation state on inhibitor binding. The design of highly potent AP inhibitors, that could be the basis for drug leads require a deep knowledge of the protonation state of the active site residues induced by their presence. This vital issue has been tackled by experimental techniques like NMR, X-ray crystallography, calorimetric and binding kinetic techniques. Recently, we have developed a protocol that combines monitoring the pH effect on binding affinities by SPR methods and rationalization of the results by molecular mechanics based calculations. We have used this combined method on BACE-1 and HIV-1 PR, two important therapeutic targets. Our calculations are able to reproduce the inhibitor binding trends to either enzyme upon a pH increase. The results indicate that inhibitors that differ in the Asp dyad binding fragments will present different binding affinity trends upon a pH increase. Our calculations have enabled us to predict the protonation states at different pH values that underlie the above mentioned trends. We have found out that these results have many implications not only for in silico hit screening campaigns aimed at finding high affinity binders, but also (in the case of BACE-1) for the discovery of cell active compounds.

On a Possible Neutral Charge State for the Catalytic Dyad in β-Secretase When Bound to Hydroxyethylene Transition State Analogue Inhibitors

β-Secretase is one of the aspartic proteases involved in the formation of amyloid plaques in Alzheimers disease patients. Our previous results using a combination of surface plasmon resonance experiments with molecular modeling calculations suggested that the Asp dyad in β-secretase bound to hydroxylethylene containing inhibitors adopts a neutral charged state. In this work, we show that the Asp dyad diprotonated state reproduced the binding ranking of a set of these inhibitors better than alternative protonation states.

Synthesis and Characterization of a New RXR Agonist Based on the 6-tert-Butyl-1,1-dimethylindanyl Structure

A new ligand for RXR is described, which is based on a 6-tert-butyl-1,1-dimethylindanyl skeleton as bioisostere of the hydrophobic retinoid region. The Stille cross-coupling reaction allowed the attachment of the polyene side chain to the indanyl ring. Docking studies were carried out to explain the RXR binding profile of this analogue.

Effect of pH and ligand charge state on BACE-1 fragment docking performance

In this work we propose a protocol for estimating the effect of pH on the docking performance to BACE-1, which affords the charge state of the inhibitor as well as the protonation state of all ionisable residues in the protein at a given pH value. To the best of our knowledge, this is the first report of a protocol predicting the BACE-1 ligand docking poses not only at the neutral pH at which most crystallographic structures were obtained, but also at the optimal pH of the enzyme (in the acidic range), at which most of the BACE-1 binding affinity assays are performed. We have applied this protocol to a set of 23 fragment-like BACE-1 ligands that span four orders of magnitude in their binding affinities. The pKa values of the BACE-1 acidic residues deviate substantially from the estimates for model compounds in solution and display a ligand dependent variability, especially in the case of the catalytic Asp dyad residues. This outcome should have a strong bearing on the design of protocols for docking based BACE-1 screening campaigns. Finally, we were able to find an explanation for the poor docking success rate of some fragments based on the availability of anchoring points, a rationale that could help to improve hit rates in BACE-1 screening campaigns.

Synthesis and biological evaluation of 1α,25-dihydroxyvitamin D(3) analogues hydroxymethylated at C-26.

We designed by docking and synthesized two novel analogues of 1α,25-dihydroxyvitamin D(3) hydroxymethylated at C-26 (2 and 3). The syntheses were carried out by the convergent Wittig-Horner approach via epoxide 12a as a common key intermediate. The antiproliferative and transactivation potency of the compounds was evaluated in colon and breast cancer cell lines. The analogues showed a similar but reduced activity compared to 1,25(OH)(2)D(3). Analogue 3 was more potent than analogue 2, and in some assays it exhibited potency similar to that of the natural ligand.

Experimental and 'in silico' analysis of the effect of pH on HIV-1 protease inhibitor affinity: implications for the charge state of the protein ionogenic groups.

The pH dependence of the HIV-1 protease inhibitor affinity was studied by determining the interaction kinetics of a series of inhibitors at three pH values by surface plasmon resonance (SPR) biosensor analysis. The results were rationalized by molecular mechanics based protocols that have as a starting point the structures of the HIV-1 protease inhibitor complexes differing in the protonation states as predicted by our calculations. The SPR experiments indicate a variety of binding affinity pH dependencies which are rather well reproduced by our simulations. Moreover, our calculations are able to pinpoint the possible changes in the charged state of the protein binding site and of the inhibitor that underlie the observed effects of the pH on binding affinity. The combination of SPR and molecular mechanics calculations has afforded novel insights into the pH dependence of inhibitor interactions with their target. This work raises the possibility of designing inhibitors with different pH binding affinity profiles to the ones described here.

Prediction and analysis of binding affinities for chemically diverse HIV-1 PR inhibitors by the modified SAFE_p approach.

One of the biggest challenges in the “in silico” screening of enzyme ligands is to have a protocol that could predict the ligand binding free energies. In our group we have developed a very simple screening function (referred to as solvent accessibility free energy of binding predictor, SAFE_p) which we have applied previously to the study of peptidic HIV‐1 protease (HIV‐1 PR) inhibitors and later to cyclic urea type HIV‐1 PR inhibitors. In this work, we have extended the SAFE_p protocol to a chemically diverse set of HIV‐1 PR inhibitors with binding constants that differ by several orders of magnitude. The resulting function is able to reproduce the ranking and in many cases the value of the inhibitor binding affinities for the HIV‐1 PR, with accuracy comparable with that of costlier protocols. We also demonstrate that the binding pocket SAFE_p analysis can contribute to the understanding of the physical forces that participate in ligand binding. The analysis tools afforded by our protocol have allowed us to identify an induced fit phenomena mediated by the inhibitor and have demonstrated that larger fragments do not necessarily contribute the most to the binding free energy, an outcome partially brought about by the substantial role the desolvation penalty plays in the energetics of binding. Finally, we have revisited the effect of the Asp dyad protonation state on the predicted binding affinities.

Inhibitor docking screened by the modified SAFE_p scoring function: application to cyclic urea HIV-1 PR inhibitors.

Our laboratory has in the past developed a method for the prediction of ligand binding free energies to proteins, referred to as SAFE_p (Solvent free energy predictor). Previously, we have applied this protocol for the prediction of the binding free energy of peptidic and cyclic urea HIV‐1 PR inhibitors, whose X‐ray structures bound to enzyme are known. In this work, we present the first account of a docking simulation, where the ligand conformations were screened and inhibitor ranking was predicted on the basis of a modified SAFE_p approach, for a set of cyclic urea‐HIV‐1 PR complexes whose structures are not known. We show that the optimal dielectric constant for docking is rather high, in line with the values needed to reproduce some protein residue properties, like pKas. Our protocol is able to reproduce most of the observed binding ranking, even in the case that the components of the equation are not fitted to experimental data. Partition of the binding free energy into pocket and residue contributions sheds light into the importance of the inhibitors fragments and on the prediction of “hot spots” for resistance mutations.