Enthalpy screen of drug candidates

Screening approaches for the identification of Nrf2-Keap1 protein-protein interaction inhibitors targeting hot spot residues

Protein-protein interactions (PPIs) play a crucial role in most biological processes and are important targets in the development of therapeutic agents. However, small molecule drug discovery that targets PPIs remains very challenging. Targeting hot spot residues is considered the best option for inhibiting such interactions, but there are few examples of how knowledge of hot spots can be used in high throughput screening to find hit compounds. A substrate adaptor protein for a ubiquitin ligase complex, Kelch-like ECH-associated protein 1 (Keap1), negatively modulates the expression of genes involved in cellular protection against oxidative stress. Here, we focused on three arginine hot spot residues in the Keap1 substrate binding pocket (Arg380, Arg415, and Arg483), and screened the carboxylic acid library owned by Japan Tobacco Inc. for compounds that interact with the arginine residues in differential scanning fluorescence assays. Furthermore, we identified several small molecule compounds that specifically bind to the Keap1 Kelch domain hot spots by comparing binding to alanine mutant proteins (R380A, R415A, and R483A) with binding to the wild-type protein using surface plasmon resonance (SPR) screening. These compounds inhibited the protein-protein interaction between the Keap1 Kelch domain and the nuclear factor erythroid 2-related factor 2 (Nrf2) peptide, and the ubiquitination of Nrf2 catalyzed by the Cul3/RINGBox 1 E3 ligase. In addition, the binding mode of one compound (Compound 4) was determined by X-ray crystallography after validation of binding by isothermal titration calorimetry, native mass spectrometry, and nuclear magnetic resonance. Compound 4 had favorable thermodynamic properties, and noncovalently bound to Keap1 with a stoichiometry of 1:1. Our results suggest that Compound 4 could potentially be developed into effective therapeutic or preventive agents for a variety of diseases and conditions such as oxidative stress response, inflammation, and carcinogenesis. We believe that the use of a set of complementary biophysical techniques including the SPR assay with single alanine mutant of hot spots provides opportunities to identify hit compounds for developing inhibitors of PPIs.

Target-ligand binding affinity from single point enthalpy calculation and elemental composition

Reliable target-ligand binding thermodynamics data are essential for successful drug design and molecular engineering projects. Besides experimental methods, a number of theoretical approaches have been introduced for the generation of binding thermodynamics data. However, available approaches often neglect electronic effects or explicit water molecules influencing target-ligand interactions. To handle electronic effects within a reasonable time frame, we introduce a fast calculator QMH-L using a single target-ligand complex structure pre-optimized at the molecular mechanics level. QMH-L is composed of the semi-empirical quantum mechanics calculation of binding enthalpy with predicted explicit water molecules at the complex interface, and a simple descriptor based on the elemental composition of the ligand. QMH-L estimates the target-ligand binding free energy with a root mean square error (RMSE) of 0.94 kcal mol-1. The calculations also provide binding enthalpy values and they were compared with experimental binding thermodynamics data collected from the most reliable isothermal titration calorimetry studies of systems including various protein targets and challenging, large peptide ligands with a molecular weight of up to 2-3 thousand. The single point enthalpy calculations of QMH-L require modest computational resources and are based on short runs with open source and/or free software like Gromacs, Mopac, MobyWat, and Fragmenter. QMH-L can be applied for fast, automated scoring of drug candidates during a virtual screen, enthalpic engineering of new ligands or thermodynamic explanation of complex interactions.

Resveratrol Derivatives Inhibit Transthyretin Fibrillization: Structural Insights into the Interactions between Resveratrol Derivatives and Transthyretin

Hereditary ATTR amyloidosis is a disease caused by the deposition of amyloid fibrils formed by mutated transthyretin (TTR), a protein that binds to thyroid hormone in the serum, in the organs. The development of a small molecule that binds to and stabilizes TTR is a promising strategy for the treatment of ATTR amyloidosis. In the present study, we demonstrated that the resveratrol derivatives including pterostilbene available as a dietary supplement inhibit the fibrillization of V30M-TTR to the same extent as the approved drug tafamidis. Furthermore, based on a thermodynamic and X-ray crystallographic analysis, the binding of the resveratrol derivative to TTR was shown to be enthalpy-driven, with the binding enthalpy being acquired by hydrogen bonding to S117. Moreover, direct observation of hydrogen atoms by neutron crystallography provided details of the hydrogen bond network by S117 and emphasized the importance of the CH···π interaction by L110 in the ligand binding.

Structure-based development and preclinical evaluation of the SARS-CoV-2 3C-like protease inhibitor simnotrelvir

The persistent pandemic of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its variants accentuates the great demand for developing effective therapeutic agents. Here, we report the development of an orally bioavailable SARS-CoV-2 3C-like protease (3CLpro) inhibitor, namely simnotrelvir, and its preclinical evaluation, which lay the foundation for clinical trials studies as well as the conditional approval of simnotrelvir in combination with ritonavir for the treatment of COVID-19. The structure-based optimization of boceprevir, an approved HCV protease inhibitor, leads to identification of simnotrelvir that covalently inhibits SARS-CoV-2 3CLpro with an enthalpy-driven thermodynamic binding signature. Multiple enzymatic assays reveal that simnotrelvir is a potent pan-CoV 3CLpro inhibitor but has high selectivity. It effectively blocks replications of SARS-CoV-2 variants in cell-based assays and exhibits good pharmacokinetic and safety profiles in male and female rats and monkeys, leading to robust oral efficacy in a male mouse model of SARS-CoV-2 Delta infection in which it not only significantly reduces lung viral loads but also eliminates the virus from brains. The discovery of simnotrelvir thereby highlights the utility of structure-based development of marked protease inhibitors for providing a small molecule therapeutic effectively combatting human coronaviruses.

Strong Binding of C-Glycosylic1,2-Thiodisaccharides to Galectin-3─Enthalpy-Driven Affinity Enhancement by Water-Mediated Hydrogen Bonds

Galectin-3 is involved in multiple pathways of many diseases, including cancer, fibrosis, and diabetes, and it is a validated pharmaceutical target for the development of novel therapeutic agents to address unmet medical needs. Novel 1,2-thiodisaccharides with a C-glycosylic functionality were synthesized by the photoinitiated thiol-ene click reaction of O-peracylated 1-C-substituted glycals and 1-thio-glycopyranoses. Subsequent global deprotection yielded test compounds, which were studied for their binding to human galectin-3 by fluorescence polarization and isothermal titration calorimetry to show low micromolar Kd values. The best inhibitor displayed a Kd value of 8.0 μM. An analysis of the thermodynamic binding parameters revealed that the binding Gibbs free energy (ΔG) of the new inhibitors was dominated by enthalpy (ΔH). The binding mode of the four most efficient 1,2-thiodisaccharides was also studied by X-ray crystallography that uncovered the unique role of water-mediated hydrogen bonds in conferring enthalpy-driven affinity enhancement for the new inhibitors. This 1,2-thiodisaccharide-type scaffold represents a new lead for galectin-3 inhibitor discovery and offers several possibilities for further development.

Forces Driving a Magic Bullet to Its Target: Revisiting the Role of Thermodynamics in Drug Design, Development, and Optimization

Drug discovery strategies have advanced significantly towards prioritizing target selectivity to achieve the longstanding goal of identifying “magic bullets” amongst thousands of chemical molecules screened for therapeutic efficacy. A myriad of emerging and existing health threats, including the SARS-CoV-2 pandemic, alarming increase in bacterial resistance, and potentially fatal chronic ailments, such as cancer, cardiovascular disease, and neurodegeneration, have incentivized the discovery of novel therapeutics in treatment regimens. The design, development, and optimization of lead compounds represent an arduous and time-consuming process that necessitates the assessment of specific criteria and metrics derived via multidisciplinary approaches incorporating functional, structural, and energetic properties. The present review focuses on specific methodologies and technologies aimed at advancing drug development with particular emphasis on the role of thermodynamics in elucidating the underlying forces governing ligand–target interaction selectivity and specificity. In the pursuit of novel therapeutics, isothermal titration calorimetry (ITC) has been utilized extensively over the past two decades to bolster drug discovery efforts, yielding information-rich thermodynamic binding signatures. A wealth of studies recognizes the need for mining thermodynamic databases to critically examine and evaluate prospective drug candidates on the basis of available metrics. The ultimate power and utility of thermodynamics within drug discovery strategies reside in the characterization and comparison of intrinsic binding signatures that facilitate the elucidation of structural–energetic correlations which assist in lead compound identification and optimization to improve overall therapeutic efficacy.

Structural and thermodynamic analyses of interactions between death-associated protein kinase 1 and anthraquinones

Death-associated protein kinase 1 (DAPK1) is a serine/threonine protein kinase that regulates apoptosis and autophagy. DAPK1 is considered to be a therapeutic target for amyloid-β deposition, endometrial adenocarcinomas and acute ischemic stroke. Here, the potent inhibitory activity of the natural anthraquinone purpurin against DAPK1 phosphorylation is shown. Thermodynamic analysis revealed that while the binding affinity of purpurin is similar to that of CPR005231, which is a DAPK1 inhibitor with an imidazopyridazine moiety, the binding of purpurin was more enthalpically favorable. In addition, the inhibition potencies were correlated with the enthalpic changes but not with the binding affinities. Crystallographic analysis of the DAPK1-purpurin complex revealed that the formation of a hydrogen-bond network is likely to contribute to the favorable enthalpic changes and that stabilization of the glycine-rich loop may cause less favorable entropic changes. The present findings indicate that purpurin may be a good lead compound for the discovery of inhibitors of DAPK1, and the observation of enthalpic changes could provide important clues for drug development.

Determination of the binding mechanism of histone deacetylase inhibitors

This article places its focus on methods and tools enabling the elucidation of the mechanism by which ligands, small-molecule inhibitors, or substrates interact with zinc-containing bacterial or human members of the histone deacetylase family (HDACs). These methods include biochemical and biophysical approaches and can be subdivided into equilibrium and kinetic methods. More information about the exact mode of action can be obtained by combining these methods with specific mutant variants of the enzymes and/or series of structural similar ligands. All available equilibrium and kinetic data including additional information from 3D structures of HDAC-ligand complexes can be beneficially combined in a data analysis procedure called Integrated Global-Fit analysis eventually providing the most likely binding mechanism.

Role of Hydrophobicity and Charge of Amyloid-Beta Oligomer Eliminating d-Peptides in the Interaction with Amyloid-Beta Monomers

Inhibition of the self-assembly process of amyloid-beta and even more the removal of already existing toxic amyloid-beta assemblies represent promising therapeutic strategies against Alzheimer's disease. To approach this aim, we selected a d-enantiomeric peptide by phage-display based on the interaction with amyloid-beta monomers. This lead compound was successfully optimized by peptide microarrays with respect to its affinity and specificity to the target resulting in d-peptides with both increased hydrophobicity and net charge. Here, we present a detailed biophysical characterization of the interactions between these optimized d-peptides and amyloid-beta monomers in comparison to the original lead compound in order to obtain a more thorough understanding of the physicochemical determinants of the interactions. Kinetics and apparent stoichiometry of complex formation were studied using surface plasmon resonance. Potential modes of binding to amyloid-beta were identified, and the influences of ionic strength on complex stability, as well as on the inhibitory effect on amyloid-beta aggregation were investigated. The results reveal a very different mode of interaction of the optimized d-peptides based on a combination of electrostatic and hydrophobic interactions as compared to the mostly electrostatically driven interaction of the lead compound. These conclusions were supported by the thermodynamic profiles of the interaction between optimized d-peptides and Aβ monomers, which indicate an increase in binding entropy with respect to the lead compound.

Design of Potent pan-IAP and Lys-Covalent XIAP Selective Inhibitors Using a Thermodynamics Driven Approach

Recently we reported that rapid determination of enthalpy of binding can be achieved for a large number of congeneric agents or in combinatorial libraries fairly efficiently. We show that using a thermodynamic Craig plot can be very useful in dissecting the enthalpy and entropy contribution of different substituents on a common scaffold, in order to design potent, selective, or pan-active compounds. In our implementation, the approach identified a critical Lys residue in the BIR3 domain of XIAP. We report for the first time that it is possible to target such residue covalently to attain potent and selective agents. Preliminary cellular studies in various models of leukemia, multiple myeloma, and pancreatic cancers suggest that the derived agents possess a potentially intriguing pattern of activity, especially for cell lines that are resistant to the pan-IAP antagonist and clinical candidate LCL161.

Enthalpy-Based Screening of Focused Combinatorial Libraries for the Identification of Potent and Selective Ligands

In modern drug discovery, the ability of biophysical methods, including nuclear magnetic resonance spectroscopy or surface plasmon resonance, to detect and characterize ligand-protein interactions accurately and unambiguously makes these approaches preferred versus conventional biochemical high-throughput screening of large collections of compounds. Nonetheless, ligand screening strategies that address simultaneously potency and selectivity have not yet been fully developed. In this work, we propose a novel method for screening large collections of combinatorial libraries using enthalpy measurements as a primary screening technique. We demonstrate that selecting binders that are driven by enthalpy (ΔH) results in agents that are not only potent but also more selective for a given target. This general and novel approach, we termed ΔH screening of fPOS (enthalpy screening of focused positional scanning library), combines the principles of focused combinatorial chemistry with rapid calorimetry measurements to efficiently identify potent and selective inhibitors.