Adding calorimetric data to decision making in lead discovery: a hot tip

A review of biophysical strategies to investigate protein-ligand binding: What have we employed?

The protein-ligand binding frequently occurs in living organisms and plays a crucial role in the execution of the functions of proteins and drugs. It is also an indispensable part of drug discovery and screening. While the methods for investigating protein-ligand binding are diverse, each has its own objectives, strengths, and limitations, which all influence the choice of method. Many studies concentrate on one or a few specific methods, suggesting that comprehensive summaries are lacking. Therefore in this review, these methods are comprehensively summarized and are discussed in detail: prediction and simulation methods, thermal and thermodynamic methods, spectroscopic methods, methods of determining three-dimensional structures of the complex, mass spectrometry-based methods and others. It is also important to integrate these methods based on the specific objectives of the research. With the aim of advancing pharmaceutical research, this review seeks to deepen the understanding of the protein-ligand binding process.

Ensembles of interconverting protein complexes with multiple interaction domains

Many critical biological processes depend on protein complexes that exist as ensembles of subcomplexes rather than a discrete complex. The subcomplexes dynamically interconvert with one another, and the ability to accurately resolve the composition of the diverse molecular species in the ensemble is crucial for understanding the contribution of each subcomplex to the overall function of the protein complex. Advances in computational programs have made it possible to predict the various molecular species in these ensembles, but experimental approaches to identify the pool of subcomplexes and associated stoichiometries are often challenging. This review highlights some experimental approaches that can be used to resolve the diverse molecular species in protein complexes that exist as ensembles of sub complexes.

Recent advances in computational and experimental protein-ligand affinity determination techniques

INTRODUCTION

Modern drug discovery revolves around designing ligands that target the chosen biomolecule, typically proteins. For this, the evaluation of affinities of putative ligands is crucial. This has given rise to a multitude of dedicated computational and experimental methods that are constantly being developed and improved.

AREAS COVERED

In this review, the authors reassess both the industry mainstays and the newest trends among the methods for protein - small-molecule affinity determination. They discuss both computational affinity predictions and experimental techniques, describing their basic principles, main limitations, and advantages. Together, this serves as initial guide to the currently most popular and cutting-edge ligand-binding assays employed in rational drug design.

EXPERT OPINION

The affinity determination methods continue to develop toward miniaturization, high-throughput, and in-cell application. Moreover, the availability of data analysis tools has been constantly increasing. Nevertheless, cross-verification of data using at least two different techniques and careful result interpretation remain of utmost importance.

Heliorhodopsin-mediated light-modulation of ABC transporter

Heliorhodopsins (HeRs) have been hypothesized to have widespread functions. Recently, the functions for few HeRs have been revealed; however, the hypothetical functions remain largely unknown. Herein, we investigate light-modulation of heterodimeric multidrug resistance ATP-binding cassette transporters (OmrDE) mediated by Omithinimicrobium cerasi HeR. In this study, we classifiy genes flanking the HeR-encoding genes and identify highly conservative residues for protein-protein interactions. Our results reveal that the interaction between OcHeR and OmrDE shows positive cooperatively sequential binding through thermodynamic parameters. Moreover, light-induced OcHeR upregulates OmrDE drug transportation. Hence, the binding may be crucial to drug resistance in O. cerasi as it survives in a drug-containing habitat. Overall, we unveil a function of HeR as regulatory rhodopsin for multidrug resistance. Our findings suggest potential applications in optogenetic technology.

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin

Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.

3-Amide-β-carbolines block the cell cycle by targeting CDK2 and DNA in tumor cells potentially as anti-mitotic agents

β-Carboline alkaloids are natural and synthetic products with outstanding antitumor activity. C3 substituted and dimerized β-carbolines exert excellent antitumor activity. In the present research, 37 β-carboline derivatives were synthesized and characterized. Their cytotoxicity, cell cycle, apoptosis, and CDK2- and DNA-binding affinity were evaluated. β-Carboline monomer M3 and dimer D4 showed selective activity and higher cytotoxicity in tumor cells than in normal cells. Structure-activity relationships (SAR) indicated that the amide group at C3 enhanced the antitumor activity. M3 blocked the A549 (IC50 = 1.44 ± 1.10 μM) cell cycle in the S phase and inhibited A549 cell migration, while D4 blocked the HepG2 (IC50 = 2.84 ± 0.73 μM) cell cycle in the G0/G1 phase, both of which ultimately induced apoptosis. Furthermore, associations of M3 and D4 with CDK2 and DNA were proven by network pharmacology analysis, molecular docking, and western blotting. The expression level of CDK2 was downregulated in M3-treated A549 cells and D4-treated HepG2 cells. Moreover, M3 and D4 interact with DNA and CDK2 at sub-micromolar concentrations in endothermic interactions caused by entropy-driven adsorption processes, which means that the favorable entropy change (ΔS > 0) overcomes the unfavorable enthalpy change (ΔH > 0) and drives the spontaneous reaction (ΔG < 0). Overall, these results clarified the antitumor mechanisms of M3 and D4 through disrupting the cell cycle by binding DNA and CDK2, which demonstrated the potential of M3 and D4 as novel antiproliferative drugs targeting mitosis.

Accurate calculation of affinity changes to the close state of influenza A M2 transmembrane domain in response to subtle structural changes of adamantyl amines using free energy perturbation methods in different lipid bilayers

Experimental binding free energies of 27 adamantyl amines against the influenza M2(22-46) WT tetramer, in its closed form at pH 8, were measured by ITC in DPC micelles. The measured Kd's range is ~44 while the antiviral potencies (IC50) range is ~750 with a good correlation between binding free energies computed with Kd and IC50 values (r = 0.76). We explored with MD simulations (ff19sb, CHARMM36m) the binding profile of complexes with strong, moderate and weak binders embedded in DMPC, DPPC, POPC or a viral mimetic membrane and using different experimental starting structures of M2. To predict accurately differences in binding free energy in response to subtle changes in the structure of the ligands, we performed 18 alchemical perturbative single topology FEP/MD NPT simulations (OPLS2005) using the BAR estimator (Desmond software) and 20 dual topology calculations TI/MD NVT simulations (ff19sb) using the MBAR estimator (Amber software) for adamantyl amines in complex with M2(22-46) WT in DMPC, DPPC, POPC. We observed that both methods with all lipids show a very good correlation between the experimental and calculated relative binding free energies (r = 0.77-0.87, mue = 0.36-0.92 kcal mol-1) with the highest performance achieved with TI/MBAR and lowest performance with FEP/BAR in DMPC bilayers. When antiviral potencies are used instead of the Kd values for computing the experimental binding free energies we obtained also good performance with both FEP/BAR (r = 0.83, mue = 0.75 kcal mol-1) and TI/MBAR (r = 0.69, mue = 0.77 kcal mol-1).

Functional study of a lytic polysaccharide monooxygenase MsLPMO3 from Morchella sextelata in the oxidative degradation of cellulose

Lytic polysaccharide monooxygenases (LPMOs) can improve the effectiveness with which agricultural waste is utilized. This study described the potent AA9 family protein MsLPMO3, derived from Morchella sextelata. It exhibited strong binding to phosphoric acid swollen cellulose (PASC), and had the considerable binding ability to Cu2+ with a Kd value of 2.70 μM by isothermal titration calorimetry (ITC). MsLPMO3 could also act on PASC at the C1 carbon via MALDI-TOF-MS results. Moreover, MsLPMO3 could boost the hydrolysis efficiency of corncob and wheat bran in combination with glycoside hydrolases. MsLPMO3 also exhibited strong oxidizing ability for 2,6-dimethoxyphenol (2,6-DMP), achieving the best Vmax value of 443.36 U·g-1 for pH 7.4 with a H2O2 concentration of 300 µM. The structure of MsLPMO3 was obtained using AlphaFold2, and the molecular docking results elucidated the specific interactions and key residues involved in the recognition process between MsLPMO3 and cellulose. Altogether, this study expands the knowledge of AA9 family proteins in cellulose degradation, providing valuable insights into the mechanisms of synergistic degradation of lignocellulose with cellulases.

A miniaturized mode-of-action profiling platform enables high throughput characterization of the molecular and cellular dynamics of EZH2 inhibition

The market approval of Tazemetostat (TAZVERIK) for the treatment of follicular lymphoma and epithelioid sarcoma has established "enhancer of zeste homolog 2" (EZH2) as therapeutic target in oncology. Despite their structural similarities and common mode of inhibition, Tazemetostat and other EZH2 inhibitors display differentiated pharmacological profiles based on their target residence time. Here we established high throughput screening methods based on time-resolved fluorescence energy transfer, scintillation proximity and high content analysis microscopy to quantify the biochemical and cellular binding of a chemically diverse collection of EZH2 inhibitors. These assays allowed to further characterize the interplay between EZH2 allosteric modulation by methylated histone tails (H3K27me3) and inhibitor binding, and to evaluate the impact of EZH2's clinically relevant mutant Y641N on drug target residence times. While all compounds in this study exhibited slower off-rates, those with clinical candidate status display significantly slower target residence times in wild type EZH2 and disease-related mutants. These inhibitors interact in a more entropy-driven fashion and show the most persistent effects in cellular washout and antiproliferative efficacy experiments. Our work provides mechanistic insights for the largest cohort of EZH2 inhibitors reported to date, demonstrating that-among several other binding parameters-target residence time is the best predictor of cellular efficacy.

Water network chemistry to exploit the nature of catalytic water molecules in Mtb DNA gyrase: a computational study to understand the binding mechanism of fluoroquinolones

The dynamics of DNA gyrase and mutants of DNA gyrA such as G88A, A90V, S91P, D94A, D94G, D94N, D94Y; and double-point mutant (S91P-D94G), are meticulously investigated using computational approaches. Molecular dynamics (MD) and hydration thermodynamics have shed light on the fundamental, mechanistic basis of mutations on the conformational stability of Quinolone Binding Pocket (QBP) of DNA gyrase. Analysis of MD results revealed the displacement of a single crystal water molecule (HOH201) from the catalytic site of wild-type (WT) and mutants of DNA gyrA. This prompted our research group to probe the five crystal water molecules present in the QBP of the enzyme using water thermodynamics. Hydration thermodynamics analysis revealed the displacement of HOH201 due to unstable thermodynamic signatures. Further, the analysis highlighted significant changes in thermodynamic signatures and locations of five crystal water hydration sites upon mutation. Integrated MD simulations and water thermodynamics provided promising insights into the conformational changes and inaccessibility of the catalytic water molecule that can influence the design of DNA gyrase inhibitors.Communicated by Ramaswamy H. Sarma.

Alcohol-Induced Conformation Changes and Thermodynamic Signatures in the Binding of Polyphenols to Proline-Rich Salivary Proteins

The first contact of polyphenols (tannins) with the human body occurs in the mouth, where they are known to interact with proline-rich proteins (PRPs). These interactions are important at a sensory level, especially for the development of astringency, but affect also various other biochemical processes. Employing thermodynamic measurements, fluorescence and CD spectroscopy, we investigated the binding process of the prototypical polyphenol ellagic acid (EA) to different IB-PRPs and BSA, also in the presence of ethanol, which is known to influence tannin-protein interactions. Binding of EA to BSA and the small peptide IB7-14 is weak, but very strong to IB9-37. The differences in binding strength and stoichiometry are due to differences in the binding motifs, which also lead to differences in the thermodynamic signatures of the binding process. EA binding to BSA is enthalpy-driven, whereas binding to both IB7-14 and IB9-37 is entropy-driven. The presence of 10 vol.% EtOH, as present in wines, increases the binding constant of EA with BSA and IB7-14 drastically, but not that with IB9-37; however, it changes the binding stoichiometry. These differences can be attributed to the effect of EtOH on the conformation dynamics of the proteins and to changes in hydration properties in alcoholic solution.

The Time and Place for Nature in Drug Discovery

The case for a renewed focus on Nature in drug discovery is reviewed; not in terms of natural product screening, but how and why biomimetic molecules, especially those produced by natural processes, should deliver in the age of artificial intelligence and screening of vast collections both in vitro and in silico. The declining natural product-likeness of licensed drugs and the consequent physicochemical implications of this trend in the context of current practices are noted. To arrest these trends, the logic of seeking new bioactive agents with enhanced natural mimicry is considered; notably that molecules constructed by proteins (enzymes) are more likely to interact with other proteins (e.g., targets and transporters), a notion validated by natural products. Nature's finite number of building blocks and their interactions necessarily reduce potential numbers of structures, yet these enable expansion of chemical space with their inherent diversity of physical characteristics, pertinent to property-based design. The feasible variations on natural motifs are considered and expanded to encompass pseudo-natural products, leading to the further logical step of harnessing bioprocessing routes to access them. Together, these offer opportunities for enhancing natural mimicry, thereby bringing innovation to drug synthesis exploiting the characteristics of natural recognition processes. The potential for computational guidance to help identifying binding commonalities in the route map is a logical opportunity to enable the design of tailored molecules, with a focus on "organic/biological" rather than purely "synthetic" structures. The design and synthesis of prototype structures should pay dividends in the disposition and efficacy of the molecules, while inherently enabling greener and more sustainable manufacturing techniques.

Heterocyclic Compounds as Dipeptidyl Peptidase-IV Inhibitors with Special Emphasis on Oxadiazoles as Potent Anti-Diabetic Agents

Dipeptidyl peptidase-IV (DPP-IV) inhibitors, often known as gliptins, have been used to treat type 2 diabetes mellitus (T2DM). They may be combined with other medications as an additional treatment or used alone as a monotherapy. In addition to insulin, sulfonylureas, thiazolidinediones, and metformin, these molecules appear as possible therapeutic options. Oxadiazole rings have been employed in numerous different ways during drug development efforts. It has been shown that including them in the pharmacophore increases the amount of ligand that may be bound. The exceptional hydrogen bond acceptor properties of oxadiazoles and the distinct hydrocarbon bonding potential of their regioisomers have been established. Beside their anti-diabetic effects, oxadiazoles display a wide range of pharmacological properties. In this study, we made the assumption that molecules containing oxadiazole rings may afford a different approach to the treatment of diabetes, not only for controlling glycemic levels but also for preventing atherosclerosis progression and other complications associated with diabetes. It was observed that oxadiazole fusion with benzothiazole, 5-(2,5,2-trifluoroethoxy) phenyl, β-homophenylalanine, 2-methyl-2-{5-(4-chlorophenyl), diamine-bridged bis-coumarinyl, 5-aryl-2-(6′-nitrobenzofuran-2′-yl), nitrobenzofuran, and/or oxindole leads to potential anti-diabetic activity.

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.

Insight into Thermodynamic and Kinetic Profiles in Small-Molecule Optimization

Structure-activity relationships (SARs) and structure-property relationships (SPRs) have been considered the most important factors during the drug optimization process. For medicinal chemists, improvements in the potencies and druglike properties of small molecules are regarded as their major goals. Among them, the binding affinity and selectivity of small molecules on their targets are the most important indicators. In recent years, there has been growing interest in using thermodynamic and kinetic profiles to analyze ligand-receptor interactions, which could provide not only binding affinities but also detailed binding parameters for small-molecule optimization. In this perspective, we are trying to provide an insight into thermodynamic and kinetic profiles in small-molecule optimization. Through a highlight of strategies on the small-molecule optimization with specific cases, we aim to put forward the importance of structure-thermodynamic relationships (STRs) and structure-kinetic relationships (SKRs), which could provide more guidance to find safe and effective small-molecule drugs.

Quantifying Energetic and Entropic Pathways in Molecular Systems

When examining dynamics occurring at nonzero temperatures, both energy and entropy must be taken into account to describe activated barrier crossing events. Furthermore, good reaction coordinates need to be constructed to describe different metastable states and the transition mechanisms between them. Here we use a physics-based machine learning method called state predictive information bottleneck (SPIB) to find nonlinear reaction coordinates for three systems of varying complexity. SPIB is able to correctly predict an entropic bottleneck for an analytical flat-energy double-well system and identify the entropy- and energy-dominated pathways for an analytical four-well system. Finally, for a simulation of benzoic acid permeation through a lipid bilayer, SPIB is able to discover the entropic and energetic barriers to the permeation process. Given these results, we thus establish that SPIB is a reasonable and robust method for finding the important entropy, energy, and enthalpy barriers in physical systems, which can then be used to enhance the understanding and sampling of different activated mechanisms.

Identification and characterization of TYK2 pseudokinase domain stabilizers that allosterically inhibit TYK2 signaling

Kinase inhibition continues to be a major focus of pharmaceutical research and discovery due to the central role of these proteins in the regulation of cellular processes. One family of kinases of pharmacological interest, due to its role in activation of immunostimulatory pathways, is the Janus kinase family. Small molecule inhibitors targeting the individual kinase proteins within this family have long been sought-after therapies. High sequence and structural similarity of the family members makes selective inhibitors difficult to identify but critical because of their inter-related multiple cellular regulatory pathways. Herein, we describe the identification of inhibitors of the important Janus kinase, TYK2, a regulator of type I interferon response. In addition, the biochemical and structural confirmation of the direct interaction of these small molecules with the TYK2 pseudokinase domain is described and a potential mechanism of allosteric regulation of TYK2 activity through stabilization of the pseudokinase domain is proposed.

Calorimetric Measurements of Biological Interactions and Their Relationships to Finite Time Thermodynamics Parameters

A description and examination of the potential for calorimetry for use in exploring the entropy flows in biological and or reacting systems is presented. A calorimeter operation background is provided, and two case studies are investigated using a transient numerical simulation. The first case describes a single cell calorimeter containing a single phase material excited by heat generation source function such as joule heating. The second case is a reacting system. The basic observation parameter, the temperature, cannot be used to separate the entropy property changes and the rate of entropy production in the second case. The calculated transient response can be further analyzed to determine the equilibrium constant once the reaction equation and stoichiometric constants are specified which allows entropy property changes and the rate of entropy production to be determined. In a biological community, the equivalent of the reaction equation and a definition of an equilibrium constant are not available for all systems. The results for the two cases illustrate that using calorimetry measurements to identify the entropy flows in biological community activities requires further work to establish a framework similar to that chemical reacting systems that are based on an equilibrium type parameter.

Intragenic suppressors unravel the role of the SCREAM ACT-like domain for bHLH partner selectivity in stomatal development

Multicellular organisms develop specialized cell types to achieve complex functions of tissues and organs. The basic helix-loop-helix (bHLH) proteins act as master regulatory transcription factors of such specialized cell types. Plant stomata are cellular valves in the aerial epidermis for efficient gas exchange and water control. Stomatal differentiation is governed by sequential actions of three lineage-specific bHLH proteins, SPEECHLESS (SPCH), MUTE, and FAMA, specifying initiation and proliferation, commitment, and terminal differentiation, respectively. A broadly expressed bHLH, SCREAM (SCRM), heterodimerizes with SPCH/MUTE/FAMA and drives stomatal differentiation via switching its partners. Yet nothing is known about its heterodimerization properties or partner preference. Here, we report the role of the SCRM C-terminal ACT-like (ACTL) domain for heterodimerization selectivity. Our intragenic suppressor screen of a dominant scrm-D mutant identified the ACTL domain as a mutation hotspot. Removal of this domain or loss of its structural integrity abolishes heterodimerization with MUTE, but not with SPCH or FAMA, and selectively abrogates the MUTE direct target gene expression. Consequently, the scrm-D ACTL mutants confer massive clusters of arrested stomatal precursor cells that cannot commit to differentiation when redundancy is removed. Structural and biophysical studies further show that SPCH, MUTE, and FAMA also possess the C-terminal ACTL domain, and that ACTL•ACTL heterodimerization is sufficient for partner selectivity. Our work elucidates a role for the SCRM ACTL domain in the MUTE-governed proliferation-differentiation switch and suggests mechanistic insight into the biological function of the ACTL domain, a module uniquely associated with plant bHLH proteins, as a heterodimeric partner selectivity interface.

Machine learning reveals key ion selectivity mechanisms in polymeric membranes with subnanometer pores

Designing single-species selective membranes for high-precision separations requires a fundamental understanding of the molecular interactions governing solute transport. Here, we comprehensively assess molecular-level features that influence the separation of 18 different anions by nanoporous cellulose acetate membranes. Our analysis identifies the limitations of bulk solvation characteristics to explain ion transport, highlighted by the poor correlation between hydration energy and the measured permselectivity (R² = 0.37). Entropy-enthalpy compensation, spanning 40 kilojoules per mole, leads to a free-energy barrier (∆G^‡) variation of only ~8 kilojoules per mole across all anions. We apply machine learning to elucidate descriptors for energetic barriers from a set of 126 collected features. Notably, electrostatic features account for 75% of the overall features used to describe ∆G^‡, despite the relatively uncharged state of cellulose acetate. Our work presents an approach for studying ion transport across nanoporous membranes that could enable the design of ion-selective membranes.

Characterizing Thermodynamics of Protein-Glycosaminoglycan Interactions Using Isothermal Titration Calorimetry

It has now become increasingly clear that a complete atomic description of how biomacromolecules recognize each other requires knowledge not only of the structures of the complexes but also of how kinetics and thermodynamics drive the binding process. In particular, such knowledge is lacking for protein-glycosaminoglycan (GAG) complexes. Isothermal titration calorimetry (ITC) is the only technique that can provide all of the thermodynamic parameters-enthalpy, entropy, free energy (binding constant), and stoichiometry-from a single experiment. Here we describe different factors that must be taken into consideration in carrying out ITC titrations to obtain meaningful thermodynamic data of protein-GAG interactions.

Dopamine D3 receptor ligand suppresses the expression of levodopa-induced dyskinesia in nonhuman primate model of parkinson's disease

Parkinson's disease (PD) is a complex multisystem, chronic and so far incurable disease with significant unmet medical needs. The incidence of PD increases with aging and the expected burden will continue to escalate with our aging population. Since its discovery in the 1961 levodopa has remained the gold standard pharmacotherapy for PD. However, the progressive nature of the neurodegenerative process in and beyond the nigrostriatal system causes a multitude of side effects, including levodopa-induced dyskinesia within 5 years of therapy. Attenuating dyskinesia has been a significant challenge in the clinical management of PD. We report on a small molecule that eliminates the expression of levodopa-induced dyskinesia and significantly improves PD-like symptoms. The lead compound PD13R we discovered is a dopamine D3 receptor partial agonist with high affinity and selectivity, orally active and with desirable drug-like properties. Future studies are aimed at developing this lead compound for treating PD patients with dyskinesia.

In silico to In vivo development of a polyherbal against Haemonchus contortus

Haemonchus contortus is a major constraint in the development of small ruminant subsector due to significant production losses incurred by it. The present study explores the antiparasitic potential of three anthelmintic plants (Butea monosperma, Vitex negundo and Catharanthus roseus (L.) G.Don) against H. contortus taking albendazole as the standard. In silico molecular docking and pharmacokinetic prediction studies were conducted with known bioactive molecules of these plants (palasonin, vinblastine, vincristine, betulinic acid and ursolic acid) against Glutamate Dehydrogenase (GDH) and tubulin molecules of the parasite. Methanolic extracts of these herbs were fractionated (hexane, ethyl acetate, chloroform and methanol) and used in in vitro larvicidal studies. Based on the in vitro data, two herbal prototypes were developed and clinically tested. All the 5 ligand molecules showed better binding affnity for GDH and tubulin protein as compared with albendazole and shared similar binding site in the core of the GDH hexamer with slight variations. Albendazole approximately stacked against GLY190A residue, showing hydrophobic interactions with PRO157A and a Pi-cation electrostatic interaction with ARG390 along with four hydrogen bonds. Vincristine formed 2 pi-anionic electrostatic bonds with ASP158 of B and C subunits alongwith hydrogen bonding and hydrophobic interaction and an additional pi-anion electrostatic interaction at ASP158A for vinblastine. Albendazole bound to α-tubulin next to colchicine site whereas vinblastine is bound at the nearby laulimalide/peloruside site of the dimer. Betulinic acid showed lateral interaction between the H2-H3 loop of one alpha subunit and H10 of the adjacent alpha subunit of two tubulin dimers. Ursolic acid and palasonin bound at the intradimer N site of microtubulin involving the H1-H7 and H1-H2 zone, respectively. The in vitro studies demonstrated good dose dependent anthelmintic potential. Both the prototypes were quite efficacious in clearing the infection, keeping it to a minimal for more than 5 months, probably, through direct anthelmintic effect through GDH, tubulin depolymerization and uncoupling as well as indirectly through immunomodulation along with antioxidant and anti-inflammatory properties.

Melatonin Binding to Human NQO2 by Isothermal Titration Calorimetry

To ensure the physical interaction between a protein and its ligand, many techniques can be applied. One of them, isothermal titration calorimetry (ITC), measures the heat exchange between a forming molecular complex and its milieu. From this heat exchange, it is possible to acquire the thermodynamic parameters, the binding stoichiometry and the affinity constant (K_a) between the two interacting binding partners, which can then be used to determine the dissociation constant (K_d). We made use of ITC to determine the true K_d of melatonin for its putative receptor MT3, also known as the enzyme quinone reductase 2 (NQO2). In this chapter, we describe the step-by-step procedure for performing this experiment and extend it to 2-iodomelatonin, a melatonin derivative that was used in the initial identification and characterization of MT3. The dissociation constants of melatonin and 2-iodomelatonin toward NQO2 derived from these experiments are in line with data reported previously, albeit using alternative techniques.

Current Trends in SPR Biosensing of SARS-CoV-2 Entry Inhibitors

The emerging risk of viral diseases has triggered the search for preventive and therapeutic agents. Since the beginning of the COVID-19 pandemic, greater efforts have been devoted to investigating virus entry mechanisms into host cells. The feasibility of plasmonic sensing technologies for screening interactions of small molecules in real time, while providing the pharmacokinetic drug profiling of potential antiviral compounds, offers an advantageous approach over other biophysical methods. This review summarizes recent advancements in the drug discovery process of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) inhibitors using Surface Plasmon Resonance (SPR) biosensors. A variety of SPR assay formats are discussed according to the binding kinetics and drug efficacies of both natural products and repurposed drugs. Special attention has been given to the targeting of antiviral agents that block the receptor binding domain of the spike protein (RBD-S) and the main protease (3CLpro) of SARS-CoV-2. The functionality of plasmonic biosensors for high-throughput screening of entry virus inhibitors was also reviewed taking into account experimental parameters (binding affinities, selectivity, stability), potential limitations and future applications.

Unraveling the Thermodynamics of Ultra-tight Binding of Intrinsically Disordered Proteins

Protein interactions mediated by the intrinsically disordered proteins (IDPs) are generally associated with lower affinities compared to those between globular proteins. Here, we characterize the association between the intrinsically disordered HigA2 antitoxin and its globular target HigB2 toxin from Vibrio cholerae using competition ITC experiments. We demonstrate that this interaction reaches one of the highest affinities reported for IDP-target systems (K_D = 3 pM) and can be entirely attributed to a short, 20-residue-long interaction motif that folds into α-helix upon binding. We perform an experimentally based decomposition of the IDP-target association parameters into folding and binding contributions, which allows a direct comparison of the binding contribution with those from globular ultra-high affinity binders. We find that the HigA2-HigB2 interface is energy optimized to a similar extent as the interfaces of globular ultra-high affinity complexes, such as barnase-barstar. Evaluation of other ultra-high affinity IDP-target systems shows that a strategy based on entropy optimization can also achieve comparably high, picomolar affinities. Taken together, these examples show how IDP-target interactions achieve picomolar affinities either through enthalpy optimization (HigA2-HigB2), resembling the ultra-high affinity binding of globular proteins, or via bound-state fuzziness and entropy optimization (CcdA-CcdB, histone H1-prothymosin α).

Protein-Small Molecule Interactions in Native Mass Spectrometry

Small molecule drug discovery has been propelled by the continual development of novel scientific methodologies to occasion therapeutic advances. Although established biophysical methods can be used to obtain information regarding the molecular mechanisms underlying drug action, these approaches are often inefficient, low throughput, and ineffective in the analysis of heterogeneous systems including dynamic oligomeric assemblies and proteins that have undergone extensive post-translational modification. Native mass spectrometry can be used to probe protein-small molecule interactions with unprecedented speed and sensitivity, providing unique insights into polydisperse biomolecular systems that are commonly encountered during the drug discovery process. In this review, we describe potential and proven applications of native MS in the study of interactions between small, drug-like molecules and proteins, including large multiprotein complexes and membrane proteins. Approaches to quantify the thermodynamic and kinetic properties of ligand binding are discussed, alongside a summary of gas-phase ion activation techniques that have been used to interrogate the structure of protein-small molecule complexes. We additionally highlight some of the key areas in modern drug design for which native mass spectrometry has elicited significant advances. Future developments and applications of native mass spectrometry in drug discovery workflows are identified, including potential pathways toward studying protein-small molecule interactions on a whole-proteome scale.

Nonspecific Binding-Fundamental Concepts and Consequences for Biosensing Applications

Nature achieves differentiation of specific and nonspecific binding in molecular interactions through precise control of biomolecules in space and time. Artificial systems such as biosensors that rely on distinguishing specific molecular binding events in a sea of nonspecific interactions have struggled to overcome this issue. Despite the numerous technological advancements in biosensor technologies, nonspecific binding has remained a critical bottleneck due to the lack of a fundamental understanding of the phenomenon. To date, the identity, cause, and influence of nonspecific binding remain topics of debate within the scientific community. In this review, we discuss the evolution of the concept of nonspecific binding over the past five decades based upon the thermodynamic, intermolecular, and structural perspectives to provide classification frameworks for biomolecular interactions. Further, we introduce various theoretical models that predict the expected behavior of biosensors in physiologically relevant environments to calculate the theoretical detection limit and to optimize sensor performance. We conclude by discussing existing practical approaches to tackle the nonspecific binding challenge in vitro for biosensing platforms and how we can both address and harness nonspecific interactions for in vivo systems.

Applications of isothermal titration calorimetry in pure and applied research from 2016 to 2020

The last 5 years have seen a series of advances in the application of isothermal titration microcalorimetry (ITC) and interpretation of ITC data. ITC has played an invaluable role in understanding multiprotein complex formation including proteolysis-targeting chimeras (PROTACS), and mitochondrial autophagy receptor Nix interaction with LC3 and GABARAP. It has also helped elucidate complex allosteric communication in protein complexes like trp RNA-binding attenuation protein (TRAP) complex. Advances in kinetics analysis have enabled the calculation of kinetic rate constants from pre-existing ITC data sets. Diverse strategies have also been developed to study enzyme kinetics and enzyme-inhibitor interactions. ITC has also been applied to study small molecule solvent and solute interactions involved in extraction, separation, and purification applications including liquid-liquid separation and extractive distillation. Diverse applications of ITC have been developed from the analysis of protein instability at different temperatures, determination of enzyme kinetics in suspensions of living cells to the adsorption of uremic toxins from aqueous streams.

Uncertainty in protein-ligand binding constants: asymmetric confidence intervals versus standard errors

Equilibrium binding constants (K_b) between chemical compounds and target proteins or between interacting proteins provide a quantitative understanding of biological interaction mechanisms. Reported uncertainties of measured experimental parameters are critical for decision-making in many scientific areas, e.g., in lead compound discovery processes and in comparing computational predictions with experimental results. Uncertainties in measured K_b values are commonly represented by a symmetric normal distribution, often quoted in terms of the experimental value plus-minus the standard deviation. However, in general, the distributions of measured K_b (and equivalent K_d) values and the corresponding free energy change ΔG_b are all asymmetric to varying degree. Here, using a simulation approach, we illustrate the effect of asymmetric K_b distributions within the realm of isothermal titration calorimetry (ITC) experiments. Further we illustrate the known, but perhaps not widely appreciated, fact that when distributions of any of K_b, K_d and ΔG_b are transformed into each other, their degree of asymmetry is changed. Consequently, we recommend that a more accurate way of expressing the uncertainties of K_b, K_d, and ΔG_b values is to consistently report 95% confidence intervals, in line with other authors' suggestions. The ways to obtain such error ranges are discussed in detail and exemplified for a binding reaction obtained by ITC.

A Micromachined Picocalorimeter Sensor for Liquid Samples with Application to Chemical Reactions and Biochemistry

Calorimetry has long been used to probe the physical state of a system by measuring the heat exchanged with the environment as a result of chemical reactions or phase transitions. Application of calorimetry to microscale biological samples, however, is hampered by insufficient sensitivity and the difficulty of handling liquid samples at this scale. Here, a micromachined calorimeter sensor that is capable of resolving picowatt levels of power is described. The sensor consists of low-noise thermopiles on a thin silicon nitride membrane that allow direct differential temperature measurements between a sample and four coplanar references, which significantly reduces thermal drift. The partial pressure of water in the ambient around the sample is maintained at saturation level using a small hydrogel-lined enclosure. The materials used in the sensor and its geometry are optimized to minimize the noise equivalent power generated by the sensor in response to the temperature field that develops around a typical sample. The experimental response of the sensor is characterized as a function of thermopile dimensions and sample volume, and its capability is demonstrated by measuring the heat dissipated during an enzymatically catalyzed biochemical reaction in a microliter-sized liquid droplet. The sensor offers particular promise for quantitative measurements on biological systems.

Investigations of Lipid Binding to Acyl-CoA-Binding Proteins (ACBP) Using Isothermal Titration Calorimetry (ITC)

Isothermal titration calorimetry (ITC) is a quantitative, biophysical method to investigate intermolecular binding between biomolecules by directly measuring the heat exchange in the binding reaction. The assay is carried out in solution when the molecules interact in vitro. This allows to determine values for binding affinity (K_d), binding stoichiometry (n), as well as changes in Gibbs free energy (ΔG), entropy (ΔS), and enthalpy (ΔH). This method also addresses the kinetics of enzymatic reactions for a substrate during conversion to a product. ITC has been used to study the interactions between proteins and ligands such as those of acyl-CoA-binding proteins (ACBPs) and acyl-CoA thioesters or ACBPs with protein partners. ITC has also been used in investigating interactions between antiserum and antigen, as well as those involving RNA and DNA and other macromolecules. We describe the methods used to isolate and purify a recombinant rice ACBP (OsACBP) for ITC. To study OsACBP binding to long-chain acyl-CoA thioesters, a microcalorimeter was used at 30 °C, and the ligand (acyl-CoA thioesters or a protein partner in the first cell), was mixed with the ACBP protein solution in a second cell, for more than 40 min comprising 20 injections. Subsequently, the binding parameters including the heat-release data were analyzed and various thermodynamic parameters were calculated.

A Familiar Protein-Ligand Interaction Revisited with Multiple Methods

The interaction of hen egg white lysozyme with the trisaccharide tri-N-acetyl glucosamine has been well-characterized by biophysical methods and structural biology. In this chapter, we present a series of experiments designed to detect and quantify that interaction using several commonly available biophysical methods: thermal shift assay, fluorescence intensity, microscale thermophoresis, isothermal titration calorimetry, and surface plasmon resonance.These experiments have been used for teaching and troubleshooting in a core facility. By taking a set of representative data from several years of practical courses, we are able to demonstrate the robustness of the protocols, calculate confidence intervals for the dissociation constant from each method, and illustrate the degree of consistency between those methods when applied to a simple system in a single location by different experimenters.

Systematic Screening of Viral Entry Inhibitors Using Surface Plasmon Resonance

Surface plasmon resonance (SPR) analytical method was initially used as biosensor for analyzing diverse biomolecular interactions and recently gained important place in the drug discovery. Here, I describe the procedures for screening of inhibitors against the viral proteins using the SPR. Using the described procedures, in the past, we were able to identify several antiviral products that interfere viral-host receptor proteins interactions.

Exploring the chemistry behind protein-glycosaminoglycan conjugate: A steady-state and kinetic spectroscopy based approach

The impact of glycosaminoglycan (chondroitin sulphate, CS) on bone morphogenetic protein - 2 (BMP - 2) structure, stability (thermal and chemical), association kinetics and conformation was monitored by multiple spectroscopic techniques (UV-Visible, fluorescence and circular dichroism). The absorbance in peptide region and fluorescence intensity of BMP - 2 was quenched in presence of CS; thus, confirming the formation of a ground-state complex. As there was an increase in Stern-Volmer constant observed as a function of temperature, idea of dynamic quenching was established. However, the negligible changes in lifetime indicated static quenching; thus, making the process a combination of static-dynamic quenching. Basically, the protein - glycan interaction was driven by entropy of the system and mediated by hydrophobic interactions. Secondary structure (CD spectroscopy) of native protein was significantly affected (intensity became more negative) in presence of CS, thus, introducing more compactness in the protein. CS infused thermal and chemical stability into BMP - 2 via alteration in its conformation. The rate of association was inversely proportional to concentration of quencher (CS), which confirmed the correlation between large size (~ 5 times the size of protein) and structural complexity of CS with fewer binding sites present in BMP - 2. The rate of association in presence of urea, suggested a decrease in association rate as a function of urea concentration for 15 μM CS. Experimental evidences suggested an interaction between protein and glycan mediated by hydrophobic interactions, which deciphers structural, thermal and chemical stability into protein.

Genetic design of co-expressed Mesorhizobium loti carbonic anhydrase and chaperone GroELS to enhancing carbon dioxide sequestration

Mesorhizobium loti carbonic anhydrase (MlCA), an intrinsically high catalytic enzyme, has been employed for carbon dioxide capture and sequestration. However, recombinant expression of MlCA in Escherichia coli often forms inclusion bodies. Hence, protein partners such as fusion-tags and molecular chaperones are involved in regarding reduce the harshness of protein folding. TrxA-tag and GroELS have been chosen to co-express with MlCA in E. coli under an inducible T7 promoter or a constitutive J23100 promoter to compare productivity and activity. The results possessed that coupling protein partners effectively increased soluble MlCA up to 2.9-folds under T7 promoter, thus enhancing the CA activity by 120% and achieving a 5.2-folds turnover rate. Besides, it has also shifted the optimum temperature from 40 °C to 50 °C, promoted stability in the broad pH range (4.5 to 9.5) and the presence of various metal ions. Based on the in vitro assay and isothermal titration calorimetry (ITC) analysis, GroELS enhancing CA activity was due to change the intrinsic thermodynamic properties of the enzyme from endothermic to exothermic reaction (i.e., ∆H = 89.8 to -121.8 kJ/mol). Therefore, the collaboration of TrxA-MlCA with GroELS successfully augmented CO₂ biomineralization.

Handling the Hurdles on the Way to Anti-tuberculosis Drug Development

The global epidemic of tuberculosis (TB) imposes a sustained epidemiologic vigilance and investments in research by governments. Mycobacterium tuberculosis, the main causative agent of TB in human beings, is a very successful pathogen, being the main cause of death in the population among infectious agents. In 2018, ~10 million individuals were contaminated with this bacillus and became ill with TB, and about 1.2 million succumbed to the disease. Most of the success of the M. tuberculosis to linger in the population comes from its ability to persist in an asymptomatic latent state into the host and, in fact, the majority of the individuals are unaware of being contaminated. Even though TB is a treatable disease and is curable in most cases, the treatment is lengthy and laborious. In addition, the rise of resistance to first-line anti-TB drugs elicits a response from TB research groups to discover new chemical entities, preferably with novel mechanisms of action. The pathway to find a new TB drug, however, is arduous and has many barriers that are difficult to overcome. Fortunately, several approaches are available today to be pursued by scientists interested in anti-TB drug development, which goes from massively testing chemical compounds against mycobacteria, to discovering new molecular targets by genetic manipulation. This review presents some difficulties found along the TB drug development process and illustrates different approaches that might be used to try to identify new molecules or targets that are able to impair M. tuberculosis survival.

Simultaneous determination of thermodynamic and kinetic data by isothermal titration calorimetry

BACKGROUND

Thermodynamic and binding kinetic data increasingly support and guide the drug optimization process.

METHODS

Because ITC thermograms contain binding thermodynamic and kinetic information, an efficient protocol for the simultaneous extraction of thermodynamic and kinetic data for 1:1 protein ligand reactions from AFFINImeter kinITC in one single experiment are presented.

RESULTS

The effort to apply this protocol requires the same time as for the standard protocol but increases the precision of both thermodynamic and kinetic data.

CONCLUSIONS

The protocol enables reliable extraction of both thermodynamic and kinetic data for 1:1 protein-ligand binding reactions with improved precision compared to the 'standard protocol'.

GENERAL SIGNIFICANCE

Thermodynamic and kinetic data are recorded under exactly the same conditions in solution without any labeling or immobilization from a protein sample that is not 100% active and would otherwise render the extraction of kinetic parameters impossible.

Binding Thermodynamics to Intrinsically Disordered Protein Domains

Many proteins are intrinsically disordered or contain one or more disordered domains. These domains can participate in binding interactions with other proteins or small ligands. Binding to intrinsically disordered protein domains requires the folding or structuring of those regions such that they can establish well-defined stoichiometric interactions. Since, in such a situation binding is coupled to folding, the energetics of those two events is reflected in the measured binding thermodynamics. In this protocol, we illustrate the thermodynamic differences between binding coupled to folding and binding independent of folding for the same protein. As an example, we use the HIV-1 envelope glycoprotein gp120 that contains structured as well as disordered domains. In the experiments presented, the binding of gp120 to molecules that bind to disordered regions and trigger structuring (CD4 or MAb 17b) and to molecules that bind to structured regions and do not induce conformational structuring (MAb b12) is discussed.

Lead Identification of 8-(Methylamino)-2-oxo-1,2-dihydroquinoline Derivatives as DNA Gyrase Inhibitors: Hit-to-Lead Generation Involving Thermodynamic Evaluation

DNA gyrase and topoisomerase IV are well-validated pharmacological targets, and quinolone antibacterial drugs are marketed as their representative inhibitors. However, in recent years, resistance to these existing drugs has become a problem, and new chemical classes of antibiotics that can combat resistant strains of bacteria are strongly needed. In this study, we applied our hit-to-lead (H2L) chemistry for the identification of a new chemical class of GyrB/ParE inhibitors by efficient use of thermodynamic parameters. Investigation of the core fragments obtained by fragmentation of high-throughput screening hit compounds and subsequent expansion of the hit fragment was performed using isothermal titration calorimetry (ITC). The 8-(methylamino)-2-oxo-1,2-dihydroquinoline derivative 13e showed potent activity against Escherichia coli DNA gyrase with an IC₅₀ value of 0.0017 μM. In this study, we demonstrated the use of ITC for primary fragment screening, followed by structural optimization to obtain lead compounds, which advanced into further optimization for creating novel antibacterial agents.

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.

Dynamic control of plant water use using designed ABA receptor agonists

Drought causes crop losses worldwide, and its impact is expected to increase as the world warms. This has motivated the development of small-molecule tools for mitigating the effects of drought on agriculture. We show here that current leads are limited by poor bioactivity in wheat, a widely grown staple crop, and in tomato. To address this limitation, we combined virtual screening, x-ray crystallography, and structure-guided design to develop opabactin (OP), an abscisic acid (ABA) mimic with up to an approximately sevenfold increase in receptor affinity relative to ABA and up to 10-fold greater activity in vivo. Studies in Arabidopsis thaliana reveal a role of the type III receptor PYRABACTIN RESISTANCE-LIKE 2 for the antitranspirant efficacy of OP. Thus, virtual screening and structure-guided optimization yielded newly discovered agonists for manipulating crop abiotic stress tolerance and water use.

Theory and applications of differential scanning fluorimetry in early-stage drug discovery

Differential scanning fluorimetry (DSF) is an accessible, rapid, and economical biophysical technique that has seen many applications over the years, ranging from protein folding state detection to the identification of ligands that bind to the target protein. In this review, we discuss the theory, applications, and limitations of DSF, including the latest applications of DSF by ourselves and other researchers. We show that DSF is a powerful high-throughput tool in early drug discovery efforts. We place DSF in the context of other biophysical methods frequently used in drug discovery and highlight their benefits and downsides. We illustrate the uses of DSF in protein buffer optimization for stability, refolding, and crystallization purposes and provide several examples of each. We also show the use of DSF in a more downstream application, where it is used as an in vivo validation tool of ligand-target interaction in cell assays. Although DSF is a potent tool in buffer optimization and large chemical library screens when it comes to ligand-binding validation and optimization, orthogonal techniques are recommended as DSF is prone to false positives and negatives.

Nanocalorimeters for biomolecular analysis and cell metabolism monitoring

Nanocalorimeters, or microfabricated calorimeters, provide a promising way to characterize the thermal process of biological processes, such as biomolecule interactions and cellular metabolic activities. They enabled miniaturized heat measurement onto a chip device with potential benefits including low sample consumption, low cost, portability, and high throughput. Over the past few decades, researchers have tried to improve nanocalorimeters' performance, in terms of sensitivity, accuracy, and detection resolution, by exploring different sensing methods, thermal insulation techniques, and liquid handling methods. The enhanced devices resulted in new applications in recent years, and here we have summarized the performance parameters and applications based on categories. Finally, we have listed the current technical difficulties in nanocalorimeter research and hope for future solutions to overcome them.