Uncompetitive, adduct-forming SARM1 inhibitors are neuroprotective in preclinical models of nerve injury and disease.

Axon degeneration is an early pathological event in many neurological diseases. The identification of the nicotinamide adenine dinucleotide (NAD) hydrolase SARM1 as a central metabolic sensor and axon executioner presents an exciting opportunity to develop novel neuroprotective therapies that can prevent or halt the degenerative process, yet limited progress has been made on advancing efficacious inhibitors. We describe a class of NAD-dependent active-site SARM1 inhibitors that function by intercepting NAD hydrolysis and undergoing covalent conjugation with the reaction product adenosine diphosphate ribose (ADPR). The resulting small-molecule ADPR adducts are highly potent and confer compelling neuroprotection in preclinical models of neurological injury and disease, validating this mode of inhibition as a viable therapeutic strategy. Additionally, we show that the most potent inhibitor of CD38, a related NAD hydrolase, also functions by the same mechanism, further underscoring the broader applicability of this mechanism in developing therapies against this class of enzymes.

Group Additivity in Ligand Binding Affinity: An Alternative Approach to Ligand Efficiency.

Group additivity is a concept that has been successfully applied to a variety of thermochemical and kinetic properties. This includes drug discovery, where functional group additivity is often assumed in ligand binding. Ligand efficiency can be recast as a special case of group additivity where ΔG/HA is the group equivalent (HA is the number of non-hydrogen atoms in a ligand). Analysis of a large data set of protein-ligand binding affinities (Ki) for diverse targets shows that in general ligand binding is distinctly nonlinear. It is possible to create a group equivalent scheme for ligand binding, but only in the context of closely related proteins, at least with regard to size. This finding has broad implications for drug design from both experimental and computational points of view. It also offers a path forward for a more general scheme to assess the efficiency of ligand binding.

Protein-ligand cocrystal structures: we can do better.

There is a large body of evidence that many protein-ligand cocrystal structures contain poorly refined ligand geometries. These errors result in bound structures that have nonideal bond lengths and angles, are strained, contain improbable conformations, and have bad protein-ligand contacts. Many of these problems can be greatly reduced with better refinement models.

Validity of ligand efficiency metrics.

A recent viewpoint article (Improving the plausibility of success with inefficient metrics. ACS Med. Chem. Lett. 2014, 5, 2-5) argued that the standard definition of ligand efficiency (LE) is mathematically invalid. In this viewpoint, we address this criticism and show categorically that the definition of LE is mathematically valid. LE and other metrics such as lipophilic ligand efficiency (LLE) can be useful during the multiparameter optimization challenge faced by medicinal chemists.

Discovery of 2-[3,5-dichloro-4-(5-isopropyl-6-oxo-1,6-dihydropyridazin-3-yloxy)phenyl]-3,5-dioxo-2,3,4,5-tetrahydro[1,2,4]triazine-6-carbonitrile (MGL-3196), a Highly Selective Thyroid Hormone Receptor ß agonist in clinical trials for the treatment of dyslipidemia.

The beneficial effects of thyroid hormone (TH) on lipid levels are primarily due to its action at the thyroid hormone receptor ß (THR-ß) in the liver, while adverse effects, including cardiac effects, are mediated by thyroid hormone receptor α (THR-α). A pyridazinone series has been identified that is significantly more THR-ß selective than earlier analogues. Optimization of this series by the addition of a cyanoazauracil substituent improved both the potency and selectivity and led to MGL-3196 (53), which is 28-fold selective for THR-ß over THR-α in a functional assay. Compound 53 showed outstanding safety in a rat heart model and was efficacious in a preclinical model at doses that showed no impact on the central thyroid axis. In reported studies in healthy volunteers, 53 exhibited an excellent safety profile and decreased LDL cholesterol (LDL-C) and triglycerides (TG) at once daily oral doses of 50 mg or higher given for 2 weeks.

The role of ligand efficiency metrics in drug discovery.

The judicious application of ligand or binding efficiency metrics, which quantify the molecular properties required to obtain binding affinity for a drug target, is gaining traction in the selection and optimization of fragments, hits and leads. Retrospective analysis of recently marketed oral drugs shows that they frequently have highly optimized ligand efficiency values for their targets. Optimizing ligand efficiency metrics based on both molecular mass and lipophilicity, when set in the context of the specific target, has the potential to ameliorate the inflation of these properties that has been observed in current medicinal chemistry practice, and to increase the quality of drug candidates.

Impact of computational structure-based methods on drug discovery.

Structure-based drug design has become an indispensible tool in drug discovery. The emergence of structure-based design is due to gains in structural biology that have provided exponential growth in the number of protein crystal structures, new computational algorithms and approaches for modeling protein-ligand interactions, and the tremendous growth of raw computer power in the last 30 years. Computer modeling and simulation have made major contributions to the discovery of many groundbreaking drugs in recent years. Examples are presented that highlight the evolution of computational structure-based design methodology, and the impact of that methodology on drug discovery.

Thermodynamics of ligand binding and efficiency.

Analysis of the experimental binding thermodynamics for approximately 100 protein-ligand complexes provides important insights into the factors governing ligand affinity and efficiency. The commonly accepted correlation between enthalpy and -TΔS is clearly observed for this relatively diverse data set. It is also clear that affinity (i.e., ΔG) is not generally correlated to either enthalpy or -TΔS. This is a worrisome trend since the vast majority of computational structure-based design is carried out using interaction energies for one, or at most a few, ligand poses. As such, these energies are most closely comparable to enthalpies not free energies. Closer inspection of the data shows that in a few cases the enthalpy (or -TΔS) is correlated with free energy. It is tempting to speculate that this could be an important consideration as to why some targets are readily amenable to modeling and others are not. Additionally, analysis of the enthalpy and -TΔS efficiencies shows that the trends observed for ligand efficiencies with respect to molecular size are primarily a consequence of enthalpic, not entropic, effects.

Macrocyclic BACE inhibitors: Optimization of a micromolar hit to nanomolar leads.

We have identified macrocyclic inhibitors of the aspartic protease BACE, implicated in the etiology of Alzheimer's disease. An X-ray structure of screening hit 1 in the BACE active site revealed a hairpin conformation suggesting that constrained macrocyclic derivatives may also bind there. Several of the analogs we prepared were >100x more potent than 1, such as 7 (5 nM K(i)).

Quantum mechanical pairwise decomposition analysis of protein kinase B inhibitors: validating a new tool for guiding drug design.

Quantum mechanical semiempirical comparative binding energy analysis calculations have been carried out for a series of protein kinase B (PKB) inhibitors derived from fragment- and structure-based drug design. These protein-ligand complexes were selected because they represent a consistent set of experimental data that includes both crystal structures and affinities. Seven scoring functions were evaluated based on both the PM3 and the AM1 Hamiltonians. The optimal models obtained by partial least-squares analysis of the aligned poses are predictive as measured by a number of standard statistical criteria and by validation with an external data set. An algorithm has been developed that provides residue-based contributions to the overall binding affinity. These residue-based binding contributions can be plotted in heat maps so as to highlight the most important residues for ligand binding. In the case of these PKB inhibitors, the maps show that Met166, Thr97, Gly43, Glu114, Ala116, and Val50, among other residues, play an important role in determining binding affinity. The interaction energy map makes it easy to identify the residues that have the largest absolute effect on ligand binding. The structure-activity relationship (SAR) map highlights residues that are most critical to discriminating between more and less potent ligands. Taken together the interaction energy and the SAR maps provide useful insights into drug design that would be difficult to garner in any other way.

Hit triage using efficiency indices after screening of compound libraries in drug discovery.

Fragment-based drug discovery (FBDD) is an important new tool to understand the molecule basis of ligand-biological target interactions. By combining optimal fragments, it is often possible to construct larger molecular weight compounds that have greater potency in a shorter period of time than can been achieved by the initial screening of larger molecular weight compound libraries. Alternatively, if screening of more traditional larger libraries has occurred, then it may be possible to analyze the data during the process of hit triage in such as way as to essentially adopt a fragment-based approach in reverse. In this review, we highlight general principles associated with the efficiency indices such as Ligand Efficiency (LE) in which screening data is normalized for biophysical properties such as molecular size. We further focus on the concept of Fit Quality (FQ), which standardizes LE values across molecular weight for more realistic, direct comparison. Using these simple concepts, one can apply FBDD routinely in the stage of hit triage when evaluating the data obtained after screening of compound libraries in drug discovery.

Nonpeptide urotensin-II receptor antagonists: a new ligand class based on piperazino-phthalimide and piperazino-isoindolinone subunits.

We have discovered two related chemical series of nonpeptide urotensin-II (U-II) receptor antagonists based on piperazino-phthalimide (5 and 6) and piperazino-isoindolinone (7) scaffolds. These structure types are distinctive from those of U-II receptor antagonist series reported in the literature. Antagonist 7a exhibited single-digit nanomolar potency in rat and human cell-based functional assays, as well as strong binding to the human U-II receptor. In advanced pharmacological testing, 7a blocked the effects of U-II in vitro in a rat aortic ring assay and in vivo in a rat ear-flush model. A discussion of U-II receptor antagonist pharmacophores is presented, and a specifically defined model is suggested from tricycle 13, which has a high degree of conformational constraint.

Ligand efficiency and fragment-based drug discovery.

The use of fragment-based drug discovery (FBDD) has increased in recent years since it is more likely to produce a better optimized compound of lower molecular weight. Ligand efficiency (LE) has become important for assessing fragments, HTS hits, and resulting optimized ligands. LE is useful for comparing ligands of equal molecular weight, but is ineffective for comparisons of ligands of differing molecular weight. LE has a strong dependence on molecular size, which has led us to develop a size-independent efficiency score termed fit quality. Evaluating FBDD examples from the literature using LE and fit quality, we find that, in general, the LEs of starting fragments are greater than those of larger, more elaborated, structures. Fit quality scores, however, tend to improve upon optimization of the fragments.

Ligand binding efficiency: trends, physical basis, and implications.

Ligand efficiency (i.e., potency/size) has emerged as an important metric in drug discovery. In general, smaller, more efficient ligands are believed to have improved prospects for good drug properties (e.g., bioavailability). Our analysis of thousands of ligands across a variety of targets shows that ligand efficiency is dependent on ligand size with smaller ligands having greater efficiencies, on average, than larger ligands. We propose two primary causes for this size dependence: the inevitable reduction in the quality of fit between ligand and receptor as the ligand becomes larger and more complex and the reduction in accessible ligand surface area on a per atom basis as size increases. These results have far-ranging implications for analysis of high-throughput screening hits, fragment-based approaches to drug discovery, and even computational models of potency.

2-Amino-3,4-dihydroquinazolines as inhibitors of BACE-1 (beta-site APP cleaving enzyme): Use of structure based design to convert a micromolar hit into a nanomolar lead.

A new aspartic protease inhibitory chemotype bearing a 2-amino-3,4-dihydroquinazoline ring was identified by high-throughput screening for the inhibition of BACE-1. X-ray crystallography revealed that the exocyclic amino group participated in a hydrogen bonding array with the two catalytic aspartic acids of BACE-1 (Asp(32), Asp(228)). BACE-1 inhibitory potency was increased (0.9 microM to 11 nM K(i)) by substitution into the unoccupied S(1)' pocket.

The role of molecular size in ligand efficiency.

Ligand efficiency is a simple metric for assessing whether a ligand derives its potency from optimal fit with the protein target or simply by virtue of making many contacts. Comparison of protein-ligand binding affinities for over 8000 ligands with 28 protein targets shows conclusively that the average ligand binding affinities are not linear with molecular size. It is therefore important to scale ligand efficiencies by the size of the ligand, particularly where small ligands (e.g., fragments) are involved. We propose a simple 'fit quality' metric that removes this dependence.

Linear interaction energy models for beta-secretase (BACE) inhibitors: Role of van der Waals, electrostatic, and continuum-solvation terms.

Computing the binding affinity of a protein-ligand complex is one of the most fundamental and difficult tasks in computer-aided drug design. Many approaches for computing binding affinities can be classified as linear interaction energy (LIE) models as they rely on some type of linear fit of computed interaction energies between ligand and protein. We have examined the computed interaction energies of a series of beta-secretase (BACE) inhibitors in terms of van der Waals, coulombic, and continuum-solvation contributions to ligand binding. We have also systematically examined the effect of different protonation states of the protein and ligands. We find that the binding affinities are relatively insensitive to the protonation state of the protein when neutral ligands are considered. Inclusion of charged ligands leads to large deviations in the coulomb, solvation, and even van der Waals terms. The latter is due to increased repulsive van der Waals interactions in the complex due to the strong coulomb attraction found between oppositely charged functional groups in the protein and ligand. In general, we find that the best models are obtained when the protein is judiciously charged (e.g. Asp32-, Arg235+) and the potentially charged ligands are treated as neutral.

Assigning the protonation states of the key aspartates in ß-Secretase using QM/MM X-ray structure refinement.

ß-Secretase, a.k.a. ß-APP cleaving enzyme (BACE), is an aspartyl protease that has been implicated as a key target in the pathogenesis of Alzheimer's disease (AD). The identification of the protonation states of the key aspartates in ß-secretase is of great interest both in understanding the reaction mechanism and in guiding the design of drugs against AD. However, the resolutions of currently available crystal structures for BACE are not sufficient to determine the hydrogen atom locations. We have assigned the protonation states of the key aspartates using a novel method, QM/MM X-ray refinement. In our approach, an energy function is introduced to the refinement where the atoms in the active site are modeled by quantum mechanics (QM) and the other atoms are represented by molecular mechanics (MM). The gradients derived from the QM/MM energy function are combined with those from the X-ray target to refine the crystal structure of a complex containing BACE and an inhibitor. A total number of 8 protonation configurations of the aspartyl dyad were considered and QM/MM X-ray refinements were performed for all of them. The relative stability of the refined structures was scored by constructing the thermodynamic cycle using the energetics calculated by fully quantum mechanical self-consistent reaction field (QM/SCRF) calculations. While all 8 refined structures fit the observed electron density about equally well, we find the mono-protonated configurations to be strongly favored energetically, especially the configuration with the inner oxygen of Asp32 protonated and the hydroxyl of the inhibitor pointing towards Asp228. It was also found that these results depend on the constraints imposed by the X-ray data. We suggest that one of the strengths of this approach is that the resulting structures are a consensus of theoretical and experimental data and remark on the significance of our results in structure based drug design and mechanistic studies.