Isothermal titration calorimetry in drug discovery

Isothermal titration calorimetry analysis of the binding between the maltodextrin binding protein malE of Staphylococcus aureus with maltodextrins of various lengths

Staphylococcus aureus (S. aureus), a Gram-positive bacterium, causes a wide range of infections, and diagnosis at an early stage is challenging. Targeting the maltodextrin transporter has emerged as a promising strategy for imaging bacteria and has been able to image a wide range of bacteria including S. aureus. However, little is known about the maltodextrin transporter in S. aureus, and this prevents new S. aureus specific ligands for the maltodextrin transporter from being developed. In Gram-positive bacteria, including S. aureus, the first step of maltodextrin transport is the binding of the maltodextrin-binding protein malE to maltodextrins. Thus, understanding the binding affinity and characteristics of malE from S. aureus is important to developing efficient maltodextrin-based imaging probes. We evaluated the affinity of malE of S. aureus to maltodextrins of various lengths. MalE of S. aureus (SAmalE) was expressed in E. coli BL21(DE3) and purified by Ni-NTA resin. The affinities of SAmalE to maltodextrins were evaluated with isothermal titration calorimetry. SAmalE has low affinity to maltose but binds to maltotriose and longer maltodextrins up to maltoheptaose with affinities up to Ka = 9.02 ± 0.49 × 105 M-1. SAmalE binding to maltotriose-maltoheptaose was exothermic and fit a single-binding site model. The van't Hoff enthalpy in the binding reaction of SAmalE with maltotriose was 9.9 ± 1.3 kcal/mol, and the highest affinity of SAmalE was observed with maltotetraose with Ka = 9.02 ± 0.49 × 105 M-1. In the plot of ΔH-T*ΔS, the of Enthalpy-Entropy Compensation effect was observed in binding reaction of SAmalE to maltodextrins. Acarbose and maltotetraiol bind with SAmalE indicating that SAmalE is tolerant of modifications on both the reducing and non-reducing ends of maltodextrins. Our results show that unlike ECmalE and similar to the maltodextrin binding protein of Streptococci, SAmalE primarily binds to maltodextrins via hydrogen bonds. This is distinct from the maltodextrin binding protein of Streptococci, SAmalE that binds to maltotetraiol with high affinity. Understanding the binding characteristics and tolerance to maltodextrins modifications by maltodextrin binding proteins will hopefully provide the basis for developing bacterial species-specific maltodextrin-based imaging probes.

Solid-phase extraction with the functionalization of calcium-sensing receptors onto magnetic microspheres as an affinity probe can capture ligands selectively from herbal extract

Magnetic solid phase extraction with the functionalization of protein onto micro- or nano-particles as a probe is favorable for the discovery of new drugs from complicated natural products. Herein, we aimed to develop a rapid method by immobilizing halogenated alkane dehalogenase (Halo)-tagged calcium-sensing receptor (CaSR) directly out of crude cell lysates onto the surface of magnetic microspheres (MM) with no need to purify protein. Thereby we achieved CaSR-functionalized MM for revealing adsorption characteristics of agonist neomycin and screening ligands from herbal medicine Radix Astragali (RA). About 43.87 mg CaSR could be immobilized per 1 g MM within 30 min, and the acquired CaSR-functionalized MM showed good stability and activity for 4 weeks. The maximum adsorption capacity of neomycin on CaSR-functionalized MM was determined as 4.70 × 10-4 ~ 3.96 × 10-4 mol/g within 277 ~ 310 K, and its adsorption isotherm characteristics described best by the Temkin model were further validated using isothermal titration calorimetry. It was inferred that CaSR's affinity for neomycin was driven by electrostatic forces in a spontaneous process when the system reached an equilibrium state. Moreover, the ligands from the RA extract were screened, three of which were assigned as astragaloside IV, ononin, and calycosin based on HPLC-MS. Our findings demonstrated that the functionalization of a receptor onto magnetic materials designed as an affinity probe has the capability to recognize its agonist and capture the ligands selectively from complex matrices like herbs.

Histidine N1-position-specific methyltransferase CARNMT1 targets C3H zinc finger proteins and modulates RNA metabolism

Histidine (His) residues are methylated in various proteins, but their roles and regulation mechanisms remain unknown. Here, we show that carnosine N-methyltransferase 1 (CARNMT1), a known His methyltransferase of dipeptide carnosine (βAla-His), is a major His N1-position-specific methyltransferase. We found that 52 His sites in 20 proteins underwent CARNMT1-mediated methylation. The consensus methylation site for CARNMT1 was identified as Cx(F/Y)xH, a C3H zinc finger (C3H ZF) motif. CARNMT1-deficient and catalytically inactive mutant mice showed embryonic lethality. Among the CARNMT1 target C3H ZF proteins, RNA degradation mediated by Roquin and tristetraprolin (TTP) was affected by CARNMT1 and its enzymatic activity. Furthermore, the recognition of the 3' splice site of the CARNMT1 target C3H ZF protein U2AF1 was perturbed, and pre-mRNA alternative splicing (AS) was affected by CARNMT1 deficiency. These findings indicate that CARNMT1-mediated protein His methylation, which is essential for embryogenesis, plays roles in diverse aspects of RNA metabolism by targeting C3H ZF-type RNA-binding proteins and modulating their functions, including pre-mRNA AS and mRNA degradation regulation.

A Phenylboronic Acid-Based Transition State Analogue Yields Nanomolar Inhibition of Mandelate Racemase

Mandelate racemase (MR) catalyzes the Mg²⁺-dependent interconversion of (R)- and (S)-mandelate by stabilizing the altered substrate in the transition state (TS) by ∼26 kcal/mol. The enzyme has been employed as a model to explore the limits to which the free energy of TS stabilization may be captured by TS analogues to effect strong binding. Herein, we determined the thermodynamic parameters accompanying binding of a series of bromo-, chloro-, and fluoro-substituted phenylboronic acids (PBAs) by MR and found that binding was predominately driven by favorable entropy changes. 3,4-Dichloro-PBA was discovered to be the most potent inhibitor yet identified for MR, binding with a K_d^app value of 11 ± 2 nM and exceeding the binding of the substrate by ∼72,000-fold. The ΔC_p value accompanying binding (-488 ± 18 cal·mol^-1 K^-1) suggested that dispersion forces contribute significantly to the binding. The pH-dependence of the inhibition revealed that MR preferentially binds the anionic, tetrahedral form of 3,4-dichloro-PBA with a pH-independent K_i value of 5.7 ± 0.5 nM, which was consistent with the observed upfield shift of the ¹¹B NMR signal. The linear free energy relationship between log(k_cat/K_m) and log(1/K_i) for wild-type and 11 MR variants binding 3,4-dichloro-PBA had a slope of 0.8 ± 0.2, indicating that MR recognizes the inhibitor as an analogue of the TS. Hence, halogen substitution may be utilized to capture additional free energy of TS stabilization arising from dispersion forces to enhance the binding of boronic acid inhibitors by MR.

Binding of common organic UV-filters to the thyroid hormone transport protein transthyretin using in vitro and in silico studies: Potential implications in health

Several anthropogenic contaminants have been identified as competing with the thyroid hormone thyroxine (T₄) for binding to transport proteins as transthyretin (TTR). This binding can potentially create toxicity mechanisms posing a threat to human health. Many organic UV filters (UVFs) and paraben preservatives (PBs), widely used in personal care products, are chemicals of emerging concern due to their adverse effects as potential thyroid-disrupting compounds. Recently, organic UVFs have been found in paired maternal and fetal samples and PBs have been detected in placenta, which opens the possibility of the involvement of TTR in the transfer of these chemicals across physiological barriers. We aimed to investigate a discrete set of organic UVFs and PBs to identify novel TTR binders. The binding affinities of target UVFs towards TTR were evaluated using in vitro T₄ competitive binding assays. The ligand-TTR affinities were determined by isothermal titration calorimetry (ITC) and compared with known TTR ligands. In parallel, computational studies were used to predict the 3-D structures of the binding modes of these chemicals to TTR. Some organic UVFs, compounds 2,2',4,4'-tetrahydroxybenzophenone (BP2, Kd = 0.43 μM); 2,4-dihydroxybenzophenone (BP1, Kd = 0.60 μM); 4,4'-dihydroxybenzophenone (4DHB, Kd = 0.83 μM), and 4-hydroxybenzophenone (4HB, Kd = 0.93 μM), were found to display a high affinity to TTR, being BP2 the strongest TTR binder (ΔH = -14.93 Kcal/mol). Finally, a correlation was found between the experimental ITC data and the TTR-ligand docking scores obtained by computational studies. The approach integrating in vitro assays and in silico methods constituted a useful tool to find TTR binders among common organic UVFs. Further studies on the involvement of the transporter protein TTR in assisting the transplacental transfer of these chemicals across physiological barriers and the long-term consequences of prenatal exposure to them should be pursued.

The effect of Mg2+ on Ca2+ binding to cardiac troponin C in hypertrophic cardiomyopathy associated TNNC1 variants

Cardiac troponin C (cTnC) is the critical Ca2+ -sensing component of the troponin complex. Binding of Ca2+ to cTnC triggers a cascade of conformational changes within the myofilament that culminate in force production. Hypertrophic cardiomyopathy (HCM)-associated TNNC1 variants generally induce a greater degree and duration of Ca2+ binding, which may underly the hypertrophic phenotype. Regulation of contraction has long been thought to occur exclusively through Ca2+ binding to site II of cTnC. However, work by several groups including ours suggest that Mg2+ , which is several orders of magnitude more abundant in the cell than Ca2+ , may compete for binding to the same cTnC regulatory site. We previously used isothermal titration calorimetry (ITC) to demonstrate that physiological concentrations of Mg2+ may decrease site II Ca2+ -binding in both N-terminal and full-length cTnC. Here, we explore the binding of Ca2+ and Mg2+ to cTnC harbouring a series of TNNC1 variants thought to be causal in HCM. ITC and thermodynamic integration (TI) simulations show that A8V, L29Q and A31S elevate the affinity for both Ca2+ and Mg2+ . Further, L48Q, Q50R and C84Y that are adjacent to the EF hand binding motif of site II have a more significant effect on affinity and the thermodynamics of the binding interaction. To the best of our knowledge, this work is the first to explore the role of Mg2+ in modifying the Ca2+ affinity of cTnC mutations linked to HCM. Our results indicate a physiologically significant role for cellular Mg2+ both at baseline and when elevated on modifying the Ca2+ binding properties of cTnC and the subsequent conformational changes which precede cardiac contraction.

Exiting the tunnel of uncertainty: crystal soak to validated hit

Crystallographic fragment screens provide an efficient and effective way to identify small-molecule ligands of a crystallized protein. Due to their low molecular weight, such hits tend to have low, often unquantifiable, affinity for their target, complicating the twin challenges of validating the hits as authentic solution-phase ligands of the target and identifying the `best' hit(s) for further elaboration. In this article, approaches that address these challenges are assessed. Using retrospective analysis of a recent ATAD2 hit-identification campaign, alongside other examples of successful fragment-screening campaigns, it is suggested that hit validation and prioritization are best achieved by a `triangulation' approach in which the results of multiple available biochemical and biophysical techniques are correlated to develop qualitative structure-activity relationships (SARs). Such qualitative SARs may indeed be the only means by which to navigate a project through the tunnel of uncertainty that prevails before on-scale biophysical, biochemical and/or biological measurements become possible.

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.

Stability of the Retinoid X Receptor-α Homodimer in the Presence and Absence of Rexinoid and Coactivator Peptide

Differential scanning calorimetry and differential scanning fluorimetry were used to measure the thermal stability of human retinoid X receptor-α ligand binding domain (RXRα LBD) homodimer in the absence or presence of rexinoid and coactivator peptide, GRIP-1. The apo-RXRα LBD homodimer displayed a single thermal unfolding transition with a T_m of 58.7 °C and an unfolding enthalpy (ΔH) of 673 kJ/mol (12.5 J/g), much lower than average value (35 J/g) of small globular proteins. Using a heat capacity change (ΔC_p) of 15 kJ/(mol K) determined by measurements at different pH values, the free energy of unfolding (ΔG) of the native state was 33 kJ/mol at 37 °C. Rexinoid binding to the apo-homodimer increased T_m by 5 to 9 °C and increased the ΔG of the native homodimer by 12 to 20 kJ/mol at 37 °C, consistent with the nanomolar dissociation constant (K_d) of the rexinoids. GRIP-1 binding to holo-homodimers containing rexinoid resulted in additional increases in ΔG of 14 kJ/mol, a value that was the same for all three rexinoids. Binding of rexinoid and GRIP-1 resulted in a combined 50% increase in unfolding enthalpy, consistent with reduced structural fluidity and more compact folding observed in other published structural studies. The complexes of UAB110 and UAB111 are each more stable than the UAB30 complex by 8 kJ/mol due to enhanced hydrophobic interactions in the binding pocket because of their larger end groups. This increase in thermodynamic stability positively correlates with their improved RXR activation potency. Thermodynamic measurements are thus valuable in predicting agonist potency.

Structural stability of SARS-CoV-2 3CLpro and identification of quercetin as an inhibitor by experimental screening

The global health emergency generated by coronavirus disease 2019 (COVID-19) has prompted the search for preventive and therapeutic treatments for its pathogen, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). There are many potential targets for drug discovery and development to tackle this disease. One of these targets is the main protease, Mpro or 3CLpro, which is highly conserved among coronaviruses. 3CLpro is an essential player in the viral replication cycle, processing the large viral polyproteins and rendering the individual proteins functional. We report a biophysical characterization of the structural stability and the catalytic activity of 3CLpro from SARS-CoV-2, from which a suitable experimental in vitro molecular screening procedure has been designed. By screening of a small chemical library consisting of about 150 compounds, the natural product quercetin was identified as reasonably potent inhibitor of SARS-CoV-2 3CLpro (Ki ~ 7 μM). Quercetin could be shown to interact with 3CLpro using biophysical techniques and bind to the active site in molecular simulations. Quercetin, with well-known pharmacokinetic and ADMET properties, can be considered as a good candidate for further optimization and development, or repositioned for COVID-19 therapeutic treatment.

Water Bridges Play a Key Role in Affinity and Selectivity for Malarial Protease Falcipain-2

Falcipain-2 (FP-2) is hemoglobinase considered an attractive drug target of Plasmodium falciparum. Recently, it has been shown that peptidomimetic nitriles containing a 3-pyridyl (3Pyr) moiety at P2 display high affinity and selectivity for FP-2 with respect to human cysteine cathepsins (hCats), outperforming other P2-Pyr isomers and analogs. Further characterization demonstrated that certain P3 variants of these compounds possess micromolar inhibition of parasite growth in vitro and no cytotoxicity against human cell lines. However, the structural determinants underlying the selectivity of the 3Pyr-containing nitriles for FP-2 remain unknown. In this work, we conduct a thorough computational study combining MD simulations and free energy calculations to decipher the bases of the selectivity of the aforementioned nitriles. Our results reveal that water bridges involving the nitrogen and one carboxyl oxygen of I85 and D234 of FP-2, respectively, and the nitrogen of the neutral 3Pyr moiety, which are either less prevalent or nonexistent in the other complexes, explain the experimental activity profiles. The presence of crystallographic waters close to the bridging water positions in the studied proteases strongly supports the occurrence of such interactions. Overall, our findings suggest that selective FP-2 inhibitors can be designed by promoting water bridge formation at the bottom of the S2 subsite and/or by introducing complementary groups that displace the bridging water.

The peptide PROTAC modality: a novel strategy for targeted protein ubiquitination

Despite dramatic advances in drug discovery over the decades, effective therapeutic strategies for cancers treatment are still in urgent demands. PROteolysis TArgeting Chimera (PROTAC), a novel therapeutic modality, has been vigorously promoted in preclinical and clinical applications. Unlike small molecule PROTAC, peptide PROTAC (p-PROTAC) with advantages of high specificity and low toxicity, while avoiding the limitations of shallow binding pockets through large interacting surfaces, provides promising substitutions for E3 ubiquitin ligase complex-mediated ubiquitination of "undruggable proteins". It is worth noting that successful applications of p-PROTAC still have some obstacles, including low stability and poor membrane permeability. Hence, we highlight that p-PROTAC combined with cell-penetrating peptides, constrained conformation technique, and targeted delivery systems could be the future efforts for potential translational research.

Quantitative Determination of Ca2+-binding to Ca2+-sensor Proteins by Isothermal Titration Calorimetry

Diverse and complex molecular recognitions are central elements of signal transduction cascades. The strength and nature of these interaction modes can be determined by different experimental approaches. Among those, Isothermal titration calorimetry (ITC) offers certain advantages by providing binding constants and thermodynamic parameters from titration series without a need to label or immobilize one or more interaction partners. Furthermore, second messenger homeostasis involving Ca²⁺-ions requires in particular knowledge about stoichiometries and affinities of Ca²⁺-binding to Ca²⁺-sensor proteins or Ca²⁺-dependent regulators, which can be obtained by employing ITC. We used ITC to measure these parameters for a set of neuronal Ca²⁺-sensor proteins operating in photoreceptor cells. Here, we present a step wise protocol to (a) measure Ca²⁺ interaction with the Ca²⁺-sensor guanylate cyclase-activating protein 1, (b) to design an ITC experiment and prepare samples, (c) to remove Ca²⁺ nearly completely from Ca²⁺ binding proteins without using a chelating agent like EGTA.

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.

Concepts and Core Principles of Fragment-Based Drug Design

In this review, a general introduction to fragment-based drug design and the underlying concepts is given. General considerations and methodologies ranging from library selection/construction over biophysical screening and evaluation methods to in-depth hit qualification and subsequent optimization strategies are discussed. These principles can be generally applied to most classes of drug targets. The examples given for fragment growing, merging, and linking strategies at the end of the review are set in the fields of enzyme-inhibitor design and macromolecule–macromolecule interaction inhibition. Building upon the foundation of fragment-based drug discovery (FBDD) and its methodologies, we also highlight a few new trends in FBDD.

High-resolution structure of RGS17 suggests a role for Ca2+ in promoting the GTPase-activating protein activity by RZ subfamily members

Regulator of G protein signaling (RGS) proteins are negative regulators of G protein-coupled receptor (GPCR) signaling through their ability to act as GTPase-activating proteins (GAPs) for activated Gα subunits. Members of the RZ subfamily of RGS proteins bind to activated Gα_o, Gα_z, and Gα_i1-3 proteins in the nervous system and thereby inhibit downstream pathways, including those involved in Ca²⁺-dependent signaling. In contrast to other RGS proteins, little is known about RZ subfamily structure and regulation. Herein, we present the 1.5-Å crystal structure of RGS17, the most complete and highest-resolution structure of an RZ subfamily member to date. RGS17 cocrystallized with Ca²⁺ bound to conserved positions on the predicted Gα-binding surface of the protein. Using NMR chemical shift perturbations, we confirmed that Ca²⁺ binds in solution to the same site. Furthermore, RGS17 had greater than 55-fold higher affinity for Ca²⁺ than for Mg²⁺ Finally, we found that Ca²⁺ promotes interactions between RGS17 and activated Gα and decreases the K_m for GTP hydrolysis, potentially by altering the binding mechanism between these proteins. Taken together, these findings suggest that Ca²⁺ positively regulates RGS17, which may represent a general mechanism by which increased Ca²⁺ concentration promotes the GAP activity of the RZ subfamily, leading to RZ-mediated inhibition of Ca²⁺ signaling.

Analysis of Isothermal Titration Calorimetry Data for Complex Interactions Using I2CITC

I2CITC allows the analysis of isothermal titration calorimetry (ITC) data for complex coupled equilibria. Here we describe how, using I2CITC, ITC data for systems involving a self-aggregating ligand and a host offering one or two binding sites can be analyzed, how interaction models can be tested, and how confidence intervals for the optimized parameters can be determined.

First Report of Plant-Derived β-Sitosterol with Antithrombotic, in Vivo Anticoagulant, and Thrombus-Preventing Activities in a Mouse Model

Inhibitors of thrombin, a key enzyme in the blood coagulation cascade, are of great interest because of their selective specificity and effectiveness in anticoagulation therapy against cardiovascular disorders. The natural soybean phytosterol, β-sitosterol (BSS) demonstrated anticoagulant activity by dose-dependent inhibition of thrombin in an uncompetitive manner with a K_i value of 0.267 μM as well as by partial inhibition of thrombin-catalyzed platelet aggregation with a half-maximal inhibitory concentration (IC₅₀) value of 10.45 ± 2.88 μM against platelet-rich plasma and 9.2 ± 1.2 μM against washed platelets. An in silico study indicated binding of BSS to thrombin, which was experimentally verified by spectrofluorometric and isothermal calorimetric analyses. Under in vitro conditions, BSS demonstrated thrombolytic activity by activating plasminogen, albeit it is devoid of protease (fibrinogenolytic) activity. BSS was noncytotoxic to mammalian cells, nonhemolytic, demonstrated its in vivo anticoagulant activity when administered orally, and inhibited k-carrageen-induced thrombus formation in the tails of mice. Our results suggest that dietary supplementation of BSS may help to prevent thrombosis-associated cardiovascular disorders.

Protein-ligand binding enthalpies from near-millisecond simulations: Analysis of a preorganization paradox

Calorimetric studies of protein-ligand binding sometimes yield thermodynamic data that are difficult to understand. Today, molecular simulations can be used to seek insight into such calorimetric puzzles, and, when simulations and experiments diverge, the results can usefully motivate further improvements in computational methods. Here, we apply near-millisecond duration simulations to estimate the relative binding enthalpies of four peptidic ligands with the Grb2 SH2 domain. The ligands fall into matched pairs, where one member of each pair has an added bond that preorganizes the ligand for binding and thus may be expected to favor binding entropically, due to a smaller loss in configurational entropy. Calorimetric studies have shown that the constrained ligands do in fact bind the SH2 domain more tightly than the flexible ones, but, paradoxically, the improvement in affinity for the constrained ligands is enthalpic, rather than entropic. The present enthalpy calculations yield the opposite trend, as they suggest that the flexible ligands bind more exothermically. Additionally, the small relative binding enthalpies are found to be balances of large differences in the energies of structural components such as ligand and the binding site residues. As a consequence, the deviations from experiment in the relative binding enthalpies represent small differences between these large numbers and hence may be particularly susceptible to error, due, for example, to approximations in the force field. We also computed first-order estimates of changes in configurational entropy on binding. These too are, arguably, paradoxical, as they tend to favor binding of the flexible ligands. The paradox is explained in part by the fact that the more rigid constrained ligands reduce the entropy of binding site residues more than their flexible analogs do, at least in the simulations. This result offers a rather general counterargument to the expectation that preorganized ligands should be associated with more favorable binding entropies, other things being equal.

Overview of the detection methods for equilibrium dissociation constant KD of drug-receptor interaction

Drug-receptor interaction plays an important role in a series of biological effects, such as cell proliferation, immune response, tumor metastasis, and drug delivery. Therefore, the research on drug-receptor interaction is growing rapidly. The equilibrium dissociation constant (KD) is the basic parameter to evaluate the binding property of the drug-receptor. Thus, a variety of analytical methods have been established to determine the KD values, including radioligand binding assay, surface plasmon resonance method, fluorescence energy resonance transfer method, affinity chromatography, and isothermal titration calorimetry. With the invention and innovation of new technology and analysis method, there is a deep exploration and comprehension about drug-receptor interaction. This review discusses the different methods of determining the KD values, and analyzes the applicability and the characteristic of each analytical method. Conclusively, the aim is to provide the guidance for researchers to utilize the most appropriate analytical tool to determine the KD values.

Anticoagulant mechanism, pharmacological activity, and assessment of preclinical safety of a novel fibrin(ogen)olytic serine protease from leaves of Leucas indica

The harnessing of medicinal plants containing a plethora of bioactive molecules may lead to the discovery of novel, potent and safe therapeutic agents to treat thrombosis-associated cardiovascular diseases. A 35 kDa (m/z 34747.5230) serine protease (lunathrombase) showing fibrin(ogen)olytic activity and devoid of N- and O- linked oligosaccharides was purified from an extract of aqueous leaves from L. indica. The LC-MS/MS analysis, de novo sequencing, secondary structure, and amino acid composition determination suggested the enzyme’s novel characteristic. Lunathrombase is an αβ-fibrinogenase, demonstrating anticoagulant activity with its dual inhibition of thrombin and FXa by a non-enzymatic mechanism. Spectrofluorometric and isothermal calorimetric analyses revealed the binding of lunathrombase to fibrinogen, thrombin, and/or FXa with the generation of endothermic heat. It inhibited collagen/ADP/arachidonic acid-induced mammalian platelet aggregation, and demonstrated antiplatelet activity via COX-1 inhibition and the upregulation of the cAMP level. Lunathrombase showed in vitro thrombolytic activity and was not inhibited by endogenous protease inhibitors α₂ macroglobulin and antiplasmin. Lunathrombase was non-cytotoxic to mammalian cells, non-hemolytic, and demonstrated dose-dependent (0.125–0.5 mg/kg) in vivo anticoagulant and plasma defibrinogenation activities in a rodent model. Lunathrombase (10 mg/kg) did not show toxicity or adverse pharmacological effects in treated animals.

Acute induction of anomalous and amyloidogenic blood clotting by molecular amplification of highly substoichiometric levels of bacterial lipopolysaccharide

It is well known that a variety of inflammatory diseases are accompanied by hypercoagulability, and a number of more-or-less longer-term signalling pathways have been shown to be involved. In recent work, we have suggested a direct and primary role for bacterial lipopolysaccharide (LPS) in this hypercoagulability, but it seems never to have been tested directly. Here, we show that the addition of tiny concentrations (0.2 ng l⁻¹) of bacterial LPS to both whole blood and platelet-poor plasma of normal, healthy donors leads to marked changes in the nature of the fibrin fibres so formed, as observed by ultrastructural and fluorescence microscopy (the latter implying that the fibrin is actually in an amyloid β-sheet-rich form that on stoichiometric grounds must occur autocatalytically). They resemble those seen in a number of inflammatory (and also amyloid) diseases, consistent with an involvement of LPS in their aetiology. These changes are mirrored by changes in their viscoelastic properties as measured by thromboelastography. As the terminal stages of coagulation involve the polymerization of fibrinogen into fibrin fibres, we tested whether LPS would bind to fibrinogen directly. We demonstrated this using isothermal calorimetry. Finally, we show that these changes in fibre structure are mirrored when the experiment is done simply with purified fibrinogen and thrombin (±0.2 ng l⁻¹ LPS). This ratio of concentrations of LPS : fibrinogen in vivo represents a molecular amplification by the LPS of more than 10⁸-fold, a number that is probably unparalleled in biology. The observation of a direct effect of such highly substoichiometric amounts of LPS on both fibrinogen and coagulation can account for the role of very small numbers of dormant bacteria in disease progression in a great many inflammatory conditions, and opens up this process to further mechanistic analysis and possible treatment.

Design and testing of bicyclic inhibitors of Grb7--are two cycles better than one?

Grb7 is an adapter protein involved in the propagation of signals in cancer cell migration and proliferation, and is thus a target for the development of novel anti-cancer agents. An 11-residue thioether-cyclized peptide known as G7-18NATE has previously been developed, that inhibits Grb7 via specific interactions with its SH2 domain with micromolar affinity. Here we explore whether the peptide binding is enhanced by the addition of a second linkage designed to restrain the peptide in its bound conformation and thus reduce the entropic loss upon binding. The use of an O-ally ser covalent linkage between residue positions 1 and 8 successfully enhanced the affinity, and ITC showed that the entropic loss was reduced. A peptide with thioether-cyclization exchanged for an amide linkage showed reduce affinity, though the formation of a disulfide bond between positions 1 and 8 in this peptide enhanced its binding. This study paves the way for improving the G7-18NATE scaffold for second generation inhibitors of Grb7.

Choosing a suitable method for the identification of replication origins in microbial genomes

As the replication of genomic DNA is arguably the most important task performed by a cell and given that it is controlled at the initiation stage, the events that occur at the replication origin play a central role in the cell cycle. Making sense of DNA replication origins is important for improving our capacity to study cellular processes and functions in the regulation of gene expression, genome integrity in much finer detail. Thus, clearly comprehending the positions and sequences of replication origins which are fundamental to chromosome organization and duplication is the first priority of all. In view of such important roles of replication origins, tremendous work has been aimed at identifying and testing the specificity of replication origins. A number of computational tools based on various skew types have been developed to predict replication origins. Using various in silico approaches such as Ori-Finder, and databases such as DoriC, researchers have predicted the locations of replication origins sites for thousands of bacterial chromosomes and archaeal genomes. Based on the predicted results, we should choose an effective method for identifying and confirming the interactions at origins of replication. Here we describe the main existing experimental methods that aimed to determine the replication origin regions and list some of the many the practical applications of these methods.

A combinatorial biophysical approach; FTSA and SPR for identifying small molecule ligands and PAINs

Biophysical methods have emerged as attractive screening techniques in drug discovery both as primary hit finding methodologies, as in the case of weakly active compounds such as fragments, and as orthogonal methods for hit validation for compounds discovered through conventional biochemical or cellular assays. Here we describe a dual method employing fluorescent thermal shift assay (FTSA), also known as differential scanning fluorimetry (DSF) and surface plasmon resonance (SPR), to interrogate ligands of the kinase p38α as well as several known pan-assay interference compounds (PAINs) such as aggregators, redox cyclers, and fluorescence quenchers. This combinatorial approach allows for independent verification of several biophysical parameters such as KD, kon, koff, ΔG, ΔS, and ΔH, which may further guide chemical development of a ligand series. Affinity values obtained from FTSA curves allow for insight into compound binding compared with reporting shifts in melting temperature. Ligand-p38 interaction data were in good agreement with previous literature. Aggregators and fluorescence quenchers appeared to reduce fluorescence signal in the FTSAs, causing artificially high shifts in Tm values, whereas redox compounds caused either shifts in affinity that did not agree between FTSA and SPR or a depression of FTSA signal.

A New-Class Antibacterial-Almost. Lessons in Drug Discovery and Development: A Critical Analysis of More than 50 Years of Effort toward ATPase Inhibitors of DNA Gyrase and Topoisomerase IV

The introduction into clinical practice of an ATPase inhibitor of bacterial DNA gyrase and topoisomerase IV (topo IV) would represent a new-class agent for the treatment of resistant bacterial infections. Novobiocin, the only historical member of this class, established the clinical proof of concept for this novel mechanism during the late 1950s, but its use declined rapidly and it was eventually withdrawn from the market. Despite significant and prolonged effort across the biopharmaceutical industry to develop other agents of this class, novobiocin remains the only ATPase inhibitor of gyrase and topo IV ever to progress beyond Phase I. In this review, we analyze the historical attempts to discover and develop agents within this class and highlight factors that might have hindered those efforts. Within the last 15 years, however, our technical understanding of the molecular details of the inhibition of the gyrase and topo IV ATPases, the factors governing resistance development to such inhibitors, and our knowledge of the physical properties required for robust clinical drug candidates have all matured to the point wherein the industry may now address this mechanism of action with greater confidence. The antibacterial spectrum within this class has recently been extended to begin to include serious Gram negative pathogens such as Pseudomonas aeruginosa, Acinetobacter baumannii, and Klebsiella pneumoniae. In spite of this recent technical progress, adverse economics associated with antibacterial R&D over the last 20 years has diminished industry's ability to commit the resources and perseverance needed to bring new-class agents to launch. Consequently, a number of recent efforts in the ATPase class have been derailed by organizational rather than scientific factors. Nevertheless, within this context we discuss the unique opportunity for the development of ATPase inhibitors of gyrase and topo IV as new-class antibacterial agents with broad spectrum potential.

Characterization of protein-protein interactions by isothermal titration calorimetry

The analysis of protein-protein interactions has attracted the attention of many researchers from both a fundamental point of view and a practical point of view. From a fundamental point of view, the development of an understanding of the signaling events triggered by the interaction of two or more proteins provides key information to elucidate the functioning of many cell processes. From a practical point of view, understanding protein-protein interactions at a quantitative level provides the foundation for the development of antagonists or agonists of those interactions. Isothermal Titration Calorimetry (ITC) is the only technique with the capability of measuring not only binding affinity but the enthalpic and entropic components that define affinity. Over the years, isothermal titration calorimeters have evolved in sensitivity and accuracy. Today, TA Instruments and MicroCal market instruments with the performance required to evaluate protein-protein interactions. In this methods paper, we describe general procedures to analyze heterodimeric (porcine pancreatic trypsin binding to soybean trypsin inhibitor) and homodimeric (bovine pancreatic α-chymotrypsin) protein associations by ITC.

Thermodynamic computational approach to capture molecular recognition in the binding of different inhibitors to the DNA gyrase B subunit from Escherichia coli

DNA gyrase subunit B, that catalyzes the hydrolysis of ATP, is an attractive target for the development of antibacterial drugs. This work is intended to rationalize molecular recognition at DNA gyrase B enzyme - inhibitor binding interface through the evaluation of different scoring functions in finding the correct pose and scoring properly 50 Escherichia coli DNA Gyrase B inhibitors belonging to five different classes. Improving the binding free energy calculation accuracy is further attempted by using rescoring schemes after short molecular dynamic simulations of the obtained docked complexes. These data are then compared with the corresponding experimental enzyme activity data. The results are analyzed from a structural point of view emphasizing the strengths and limitations of the techniques applied in the study.

Drug-binding energetics of human α-1-acid glycoprotein assessed by isothermal titration calorimetry and molecular docking simulations

This study utilizes sensitive, modern isothermal titration calorimetric methods to characterize the microscopic thermodynamic parameters that drive the binding of basic drugs to α-1-acid glycoprotein (AGP) and thereby rationalize the thermodynamic data in relation to docking models and crystallographic structures of the drug-AGP complexes. The binding of basic compounds from the tricyclic antidepressant series, together with miaserine, chlorpromazine, disopyramide and cimetidine, all displayed an exothermically driven binding interaction with AGP. The impact of protonation/deprotonation events, ionic strength, temperature and the individual selectivity of the A and F1*S AGP variants on drug-binding thermodynamics was characterized. A correlation plot of the thermodynamic parameters for all of the test compounds revealed that an enthalpy-entropy compensation is in effect. The exothermic binding energetics of the test compounds were driven by a combination of favorable (negative) enthalpic (∆Hº) and favorable (positive) entropic (∆Sº) contributions to the Gibbs free energy (∆Gº). Collectively, the data imply that the free energies that drive drug binding to AGP and its relationship to drug serum residency evolve from the complex interplay of enthalpic and entropic forces from interactions with explicit combinations of hydrophobic and polar side-chain sub-domains within the multi-lobed AGP ligand binding cavity.

Isothermal titration calorimetry for studying protein-ligand interactions

Isothermal titration calorimetry (ITC) is a biophysical technique that allows a thermodynamic characterization of an interactive system. It is a free in solution technique that requires no labeling, using heat as signal. ITC allows simultaneous determination of affinity K a, stoichiometry n, enthalpy change ΔH and calculation of free energy change ΔG and entropy change ΔS in one single experiment. It is the only technique that allows direct enthalpy change measurement. By accessing the enthalpy change, we get a step closer in estimating the driving forces that characterize the interaction of a protein with a ligand, information much needed in the drug discovery process.

Evaluation of the association of mercury(II) with some dicysteinyl tripeptides

The present study was undertaken to gain insight into the associations of mercury(II) with dicysteinyl tripeptides in buffered media at pH 7.4. We investigated the effects of increasing the distance between cysteinyl residues on mercury(II) associations and complex formations. The peptide-mercury(II) formation constants and their associated thermodynamic parameters in 3-(N-morpholino)propanesulfonic acid (MOPS) buffered solutions were evaluated by isothermal titration calorimetry. Complexes formed in different relative ratios of mercury(II) to cysteinyl peptides in ammonium formate buffered solutions were characterized by LTQ Orbitrap mass spectrometry. The results from these studies show that n-alkyl dicysteinyl peptides (CP 1-4), and an aryl dicysteinyl peptide (CP 5) can serve as effective "double anchors" to accommodate the coordination sites of mercury(II) to form predominantly one-to-one Hg(peptide) complexes. The aryl dicysteinyl peptide (CP 5) also forms the two-to-two Hg(2)(peptide)(2) complex. In the presence of excess peptide, Hg(peptide)(2) complexes are also detected. Notably, increasing the distance between the ligating groups or "anchor points" in CP 1-5 does not significantly affect their affinity for mercury(II). However, the enthalpy change (ΔH) values (ΔH(1) ~ -91 kJ mol(-1) and ΔH(2) ~ -66 kJ mol(-1)) for complex formation between CP 4 and 5 with mercury(II) are about one and a half times larger than the related values for CP 1, 2 and 3 (ΔH(1) ~ -66 kJ mol(-1) and ΔH(2) ~ 46 kJ mol(-1)). The corresponding entropy change (ΔS) values (ΔS(1) ~ -129 J K(-1) mol(-1) and ΔS(2) ~ -116 J K(-1) mol(-1)) of the structurally larger dicysteinyl peptides CP 4 and 5 are less entropically favorable than for CP 1, 2 and 3 (ΔS(1) ~ -48 J K(-1) mol(-1) and ΔS(2) ~ -44 J K(-1) mol(-1)). Generally, these associations result in a decrease in entropy, indicating that these peptide-mercury complexes potentially form highly ordered structures. The results from this study show that dicysteinyl tripeptides are effective in binding mercury(II) and they are promising motifs for the design of multi-cysteinyl peptides for binding more than one mercury(II) ion per peptide.

Interaction of the non-phosphorylated peptide G7-18NATE with Grb7-SH2 domain requires phosphate for enhanced affinity and specificity

Src-homology (SH2) domains are an attractive target for the inhibition of specific signalling pathways but pose the challenge of developing a truly specific inhibitor. The G7-18NATE cyclic peptide is reported to specifically inhibit the growth factor receptor bound protein 7 (Grb7) adapter protein, implicated in the progression of several cancer types, via interactions with its SH2 domain. G7-18NATE effectively inhibits the interaction of Grb7 with ErbB3 and focal adhesion kinase in cell lysates and, with the addition of a cell permeability sequence, inhibits the growth and migration of a number of breast cancer cell lines. It is thus a promising lead in the development of therapeutics targeted to Grb7. Here we investigate the degree to which G7-18NATE is specific for the Grb7-SH2 domain compared with closely related SH2 domains including those of Grb10, Grb14, and Grb2 using surface plasmon resonance. We demonstrate that G7-18NATE binds with micromolar binding affinity to Grb7-SH2 domain (K(D)  = 4-6 μm) compared with 50-200 times lower affinity for Grb10, Grb14, and Grb2 but that this specificity depends critically on the presence of phosphate in millimolar concentrations. Other differences in buffer composition, including use of Tris or 2-(N-Morpholino)ethanesulfonic acid or varying the pH, do not impact on the interaction. This suggests that under cellular conditions, G7-18NATE binds with highest affinity to Grb7. In addition, our findings demonstrate that the basis of specificity of G7-18NATE binding to the Grb7-SH2 domain is via other than intrinsic structural features of the protein, representing an unexpected mode of molecular recognition.

Combination of isothermal titration calorimetry and time-resolved luminescence for high affinity antibody-ligand interaction thermodynamics and kinetics

For experiments using synthetic ligands as probes for biological experiments, it is useful to determine the specificity and affinity of the ligands for their receptors. As ligands with higher affinities are developed (K(A)>10(8)M(-1); K(D)<10(-8)M), a new challenge arises: to measure these values accurately. Isothermal titration calorimetry measures heat produced or consumed during ligand binding, and also provides the equilibrium binding constant. However, as normally practiced, its range is limited. Displacement titration, where a competing weaker ligand is used to lower the apparent affinity of the stronger ligand, can be used to determine the binding affinity as well as the complete thermodynamic data for ligand-antibody complexes with very high affinity. These equilibrium data have been combined with kinetic measurements to yield the rate constants as well. We describe this methodology, using as an example antibody 2D12.5, which captures yttrium S-2-(4-aminobenzyl)-1, 4, 7, 10-tetraazacyclododecanetetraacetate.

In-depth biophysical analysis of interactions between therapeutic antibodies and the extracellular domain of the epidermal growth factor receptor

Targeting of the epidermal growth factor receptor (EGFR) with monoclonal antibodies has become an established antitumor strategy in clinical use or in late stages of drug development. The mAbs effector mechanisms have been widely analyzed based on in vivo or cell studies. Hereby we intend to complement these functional studies by investigating the mAb-EGFR interactions on a molecular level. Surface plasmon resonance, isothermal titration calorimetry, and static light scattering were employed to characterize the interactions of matuzumab, cetuximab, and panitumumab with the extracellular soluble form ecEGFR. The kinetic and thermodynamic determinants dissected the differences in mAbs binding mechanism toward ecEGFR. The quantitative stoichiometric data clearly demonstrated the bivalent binding of the mAbs to two ecEGFR molecules. Our results complement earlier studies on simultaneous binding of cetuximab and matuzumab. The antibodies retain their bivalent binding mode achieving a 1:2:1 complex formation. Interestingly the binding parameters remain nearly constant for the individual antibodies in this ternary assembly. In contrast the binding of panitumumab is almost exclusive either by directly blocking the accessibility for the second antibody or by negative allosteric modulation. Overall we provide a comprehensive biophysical dataset on binding parameters, the complex assembly, and relative epitope accessibility for therapeutic anti-EGFR antibodies.

Pathway and mechanism of drug binding to G-protein-coupled receptors

How drugs bind to their receptors--from initial association, through drug entry into the binding pocket, to adoption of the final bound conformation, or "pose"--has remained unknown, even for G-protein-coupled receptor modulators, which constitute one-third of all marketed drugs. We captured this pharmaceutically critical process in atomic detail using the first unbiased molecular dynamics simulations in which drug molecules spontaneously associate with G-protein-coupled receptors to achieve final poses matching those determined crystallographically. We found that several beta blockers and a beta agonist all traverse the same well-defined, dominant pathway as they bind to the β(1)- and β(2)-adrenergic receptors, initially making contact with a vestibule on each receptor's extracellular surface. Surprisingly, association with this vestibule, at a distance of 15 Å from the binding pocket, often presents the largest energetic barrier to binding, despite the fact that subsequent entry into the binding pocket requires the receptor to deform and the drug to squeeze through a narrow passage. The early barrier appears to reflect the substantial dehydration that takes place as the drug associates with the vestibule. Our atomic-level description of the binding process suggests opportunities for allosteric modulation and provides a structural foundation for future optimization of drug-receptor binding and unbinding rates.

Virtual screening for potential inhibitors of bacterial MurC and MurD ligases

Mur ligases are bacterial enzymes involved in the cytoplasmic steps of peptidoglycan biosynthesis and are viable targets for antibacterial drug discovery. We have performed virtual screening for potential ATP-competitive inhibitors targeting MurC and MurD ligases, using a protocol of consecutive hierarchical filters. Selected compounds were evaluated for inhibition of MurC and MurD ligases, and weak inhibitors possessing dual inhibitory activity have been identified. These compounds represent new scaffolds for further optimisation towards multiple Mur ligase inhibitors with improved inhibitory potency.

T-cell epitope prediction and immune complex simulation using molecular dynamics: state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.