Evaluation of fluorescence-based thermal shift assays for hit identification in drug discovery

Thermodynamic studies for drug design and screening

INTRODUCTION

A key part of drug design and development is the optimization of molecular interactions between an engineered drug candidate and its binding target. Thermodynamic characterization provides information about the balance of energetic forces driving binding interactions and is essential for understanding and optimizing molecular interactions.

AREAS COVERED

This review discusses the information that can be obtained from thermodynamic measurements and how this can be applied to the drug development process. Current approaches for the measurement and optimization of thermodynamic parameters are presented, specifically higher throughput and calorimetric methods. Relevant literature for this review was identified in part by bibliographic searches for the period 2004 - 2011 using the Science Citation Index and PUBMED and the keywords listed below.

EXPERT OPINION

The most effective drug design and development platform comes from an integrated process utilizing all available information from structural, thermodynamic and biological studies. Continuing evolution in our understanding of the energetic basis of molecular interactions and advances in thermodynamic methods for widespread application are essential to realize the goal of thermodynamically driven drug design. Comprehensive thermodynamic evaluation is vital early in the drug development process to speed drug development toward an optimal energetic interaction profile while retaining good pharmacological properties. Practical thermodynamic approaches, such as enthalpic optimization, thermodynamic optimization plots and the enthalpic efficiency index, have now matured to provide proven utility in the design process. Improved throughput in calorimetric methods remains essential for even greater integration of thermodynamics into drug design.

Changes in channel trafficking and protein stability caused by LQT2 mutations in the PAS domain of the HERG channel

Inherited human long-QT2 syndrome (LQTS) results from mutations in the gene encoding the HERG channel. Several LQT2-associated mutations have been mapped to the amino terminal cytoplasmic Per-Arnt-Sim (PAS) domain of the HERG1a channel subunit. Here we have characterized the trafficking properties of some LQT2-associated PAS domain mutants and analyzed rescue of the trafficking mutants by low temperature (27°C) or by the pore blocker drug E4031. We show that the LQT2-associated mutations in the PAS domain of the HERG channel display molecular properties that are distinct from the properties of LQT2-associated mutations in the trans-membrane region. Unlike the latter, many of the tested PAS domain LQT2-associated mutations do not result in trafficking deficiency of the channel. Moreover, the majority of the PAS domain mutations that cause trafficking deficiencies are not rescued by a pore blocking drug. We have also explored the in vitro folding stability properties of isolated mutant PAS domain proteins using a thermal unfolding fluorescence assay and a chemical unfolding assay.

Botulinum neurotoxin is shielded by NTNHA in an interlocked complex

Botulinum neurotoxins (BoNTs) are highly poisonous substances that are also effective medicines. Accidental BoNT poisoning often occurs through ingestion of Clostridium botulinum-contaminated food. Here, we present the crystal structure of a BoNT in complex with a clostridial nontoxic nonhemagglutinin (NTNHA) protein at 2.7 angstroms. Biochemical and functional studies show that NTNHA provides large and multivalent binding interfaces to protect BoNT from gastrointestinal degradation. Moreover, the structure highlights key residues in BoNT that regulate complex assembly in a pH-dependent manner. Collectively, our findings define the molecular mechanisms by which NTNHA shields BoNT in the hostile gastrointestinal environment and releases it upon entry into the circulation. These results will assist in the design of small molecules for inhibiting oral BoNT intoxication and of delivery vehicles for oral administration of biologics.

Structural and functional characterization of the product of disease-related factor H gene conversion

Numerous complement factor H (FH) mutations predispose patients to atypical hemolytic uremic syndrome (aHUS) and other disorders arising from inadequately regulated complement activation. No unifying structural or mechanistic consequences have been ascribed to these mutants beyond impaired self-cell protection. The S1191L and V1197A mutations toward the C-terminus of FH, which occur in patients singly or together, arose from gene conversion between CFH encoding FH and CFHR1 encoding FH-related 1. We show that neither single nor double mutations structurally perturbed recombinant proteins consisting of the FH C-terminal modules, 19 and 20 (FH19-20), although all three FH19-20 mutants were poor, compared to wild-type FH19-20, at promoting hemolysis of C3b-coated erythrocytes through competition with full-length FH. Indeed, our new crystal structure of the S1191L mutant of FH19-20 complexed with an activation-specific complement fragment, C3d, was nearly identical to that of the wild-type FH19-20:C3d complex, consistent with mutants binding to C3b with wild-type-like affinity. The S1191L mutation enhanced thermal stability of module 20, whereas the V1197A mutation dramatically decreased it. Thus, although mutant proteins were folded at 37 °C, they differ in conformational rigidity. Neither single substitutions nor double substitutions increased measurably the extent of FH19-20 self-association, nor did these mutations significantly affect the affinity of FH19-20 for three glycosaminoglycans, despite critical roles of module 20 in recognizing polyanionic self-surface markers. Unexpectedly, FH19-20 mutants containing Leu1191 self-associated on a heparin-coated surface to a higher degree than on surfaces coated with dermatan or chondroitin sulfates. Thus, potentially disease-related functional distinctions between mutants, and between FH and FH-related 1, may manifest in the presence of specific glycosaminoglycans.

Structure of the effector-binding domain of the arabinose repressor AraR from Bacillus subtilis

In Bacillus subtilis, the arabinose repressor AraR negatively controls the expression of genes in the metabolic pathway of arabinose-containing polysaccharides. The protein is composed of two domains of different phylogenetic origin and function: an N-terminal DNA-binding domain belonging to the GntR family and a C-terminal effector-binding domain that shows similarity to members of the GalR/LacI family. The crystal structure of the C-terminal effector-binding domain of AraR in complex with the effector L-arabinose has been determined at 2.2 Å resolution. The L-arabinose binding affinity was characterized by isothermal titration calorimetry and differential scanning fluorimetry; the K(d) value was 8.4 ± 0.4 µM. The effect of L-arabinose on the protein oligomeric state was investigated in solution and detailed analysis of the crystal identified a dimer organization which is distinctive from that of other members of the GalR/LacI family.

Advancements in high throughput biophysical technologies: applications for characterization and screening during early formulation development of monoclonal antibodies

The formulation development of monoclonal antibodies is extremely challenging, due to the diversity and complexity contained within this class of molecules. The physical and chemical properties of a monoclonal antibody dictate the behavior of the protein drug during manufacturing, storage and clinical administration. In the past few years, the use of high throughput technologies has been widely adapted to delineate unique properties of individual immunoglobulin G's (IgG's) important for their development. Numerous screening techniques have been designed to reveal physical and chemical characteristics of a protein relevant to stability under production, formulation and delivery conditions. In addition, protein stability under accelerated stresses has been utilized to predict long-term storage behavior for monoclonal antibodies in the formulation. In this review, we summarize the recent advancements in the field of biophysical technology, with a specific focus on the techniques that can be directly applied to the formulation development of monoclonal antibodies. Several case studies are also presented here to provide examples of combining existing biophysical methods with high throughput screening technology in the formulation development of monoclonal antibody drugs.

Structure of S. aureus HPPK and the discovery of a new substrate site inhibitor

The first structural and biophysical data on the folate biosynthesis pathway enzyme and drug target, 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (SaHPPK), from the pathogen Staphylococcus aureus is presented. HPPK is the second essential enzyme in the pathway catalysing the pyrophosphoryl transfer from cofactor (ATP) to the substrate (6-hydroxymethyl-7,8-dihydropterin, HMDP). In-silico screening identified 8-mercaptoguanine which was shown to bind with an equilibrium dissociation constant, K_d, of ∼13 µM as measured by isothermal titration calorimetry (ITC) and surface plasmon resonance (SPR). An IC₅₀ of ∼41 µM was determined by means of a luminescent kinase assay. In contrast to the biological substrate, the inhibitor has no requirement for magnesium or the ATP cofactor for competitive binding to the substrate site. The 1.65 Å resolution crystal structure of the inhibited complex showed that it binds in the pterin site and shares many of the key intermolecular interactions of the substrate. Chemical shift and ¹⁵N heteronuclear NMR measurements reveal that the fast motion of the pterin-binding loop (L2) is partially dampened in the SaHPPK/HMDP/α,β-methylene adenosine 5′-triphosphate (AMPCPP) ternary complex, but the ATP loop (L3) remains mobile on the µs-ms timescale. In contrast, for the SaHPPK/8-mercaptoguanine/AMPCPP ternary complex, the loop L2 becomes rigid on the fast timescale and the L3 loop also becomes more ordered – an observation that correlates with the large entropic penalty associated with inhibitor binding as revealed by ITC. NMR data, including ¹⁵N-¹H residual dipolar coupling measurements, indicate that the sulfur atom in the inhibitor is important for stabilizing and restricting important motions of the L2 and L3 catalytic loops in the inhibited ternary complex. This work describes a comprehensive analysis of a new HPPK inhibitor, and may provide a foundation for the development of novel antimicrobials targeting the folate biosynthetic pathway.

Probing the druggability of protein-protein interactions: targeting the Notch1 receptor ankyrin domain using a fragment-based approach

In order to achieve greater selectivity in drug discovery, researchers in both academia and industry are targeting cell regulatory systems. This often involves targeting the protein-protein interactions of regulatory multiprotein assemblies. Protein-protein interfaces are widely recognized to be challenging targets as they tend to be large and relatively flat, and therefore usually do not have the concave binding sites that characterize the so-called 'druggable genome'. One such prototypic multiprotein target is the Notch transcription complex, where an extensive network of protein-protein interactions stabilize the ternary complex comprising the ankyrin domain, CSL (CBF1/suppressor of Hairless/Lag-1) and MAML (Mastermind-like). Enhanced Notch activity is implicated in the development of T-ALL (T-cell acute lymphoblastic leukaemia) and selective inhibitors of Notch would be useful cancer medicines. In the present paper, we describe a fragment-based approach to explore the druggability of the ankyrin domain. Using biophysical methods and X-ray crystal structure analyses, we demonstrate that molecules can bind to the surface of the ankyrin domain at the interface region with CSL and MAML. We show that they probably represent starting points for designing larger compounds that can inhibit important protein-protein interactions that stabilize the Notch complex. Given the relatively featureless topography of the ankyrin domain, this unexpected development should encourage others to explore the druggability of such challenging multiprotein systems using fragment-based approaches.

Double-strand break-induced targeted mutagenesis in plants

Double-strand breaks are very potent inducers of DNA recombination. There is no recombination between DNA molecules unless one or two DNA strands are broken. It has become feasible to introduce double-strand breaks at specific chromosomal loci by using dedicated, redesigned endonucleases with altered DNA-binding specificities. Such breaks are mainly repaired by error-prone nonhomologous recombination pathways in somatic cells, thus frequently producing mutations at the preselected chromosomal sites. Although the art and science of reengineering protein properties have been advancing quickly, an empirical validation of new endonucleases in a particular experimental environment is essential for successful targeted mutagenesis experiments. This chapter presents methods that were developed for a comprehensive evaluation of the DNA-binding and DNA-cutting activities of homing endonucleases in maize cells; however, they can be adopted for similar evaluation studies of other endonucleases and other plant species that are amenable for Agrobacterium-mediated transformation.

Phenylalanine hydroxylase misfolding and pharmacological chaperones

Phenylketonuria (PKU) is a loss-of-function inborn error of metabolism. As many other inherited diseases the main pathologic mechanism in PKU is an enhanced tendency of the mutant phenylalanine hydroxylase (PAH) to misfold and undergo ubiquitin-dependent degradation. Recent alternative approaches with therapeutic potential for PKU aim at correcting the PAH misfolding, and in this respect pharmacological chaperones are the focus of increasing interest. These compounds, which often resemble the natural ligands and show mild competitive inhibition, can rescue the misfolded proteins by stimulating their renaturation in vivo. For PKU, a few studies have proven the stabilization of PKU-mutants in vitro, in cells, and in mice by pharmacological chaperones, which have been found either by using the tetrahydrobiopterin (BH(4)) cofactor as query structure for shape-focused virtual screening or by high-throughput screening of small compound libraries. Both approaches have revealed a number of compounds, most of which bind at the iron-binding site, competitively with respect to BH(4). Furthermore, PAH shares a number of ligands, such as BH(4), amino acid substrates and inhibitors, with the other aromatic amino acid hydroxylases: the neuronal/neuroendocrine enzymes tyrosine hydroxylase (TH) and the tryptophan hydroxylases (TPHs). Recent results indicate that the PAH-targeted pharmacological chaperones should also be tested on TH and the TPHs, and eventually be derivatized to avoid unwanted interactions with these other enzymes. After derivatization and validation in animal models, the PAH-chaperoning compounds represent novel possibilities in the treatment of PKU.

Targeting MgrA-mediated virulence regulation in Staphylococcus aureus

Increasing antibiotic resistance in human pathogens necessitates the development of new approaches against infections. Targeting virulence regulation at the transcriptional level represents a promising strategy yet to be explored. A global transcriptional regulator, MgrA in Staphylococcus aureus, was identified previously as a key virulence determinant. We have performed a fluorescence anisotropy (FA)-based high-throughput screen that identified 5, 5-methylenedisalicylic acid (MDSA), which blocks the DNA binding of MgrA. MDSA represses the expression of α-toxin that is up-regulated by MgrA and activates the transcription of protein A, a gene down-regulated by MgrA. MDSA alters bacterial antibiotic susceptibilities via an MgrA-dependent pathway. A mouse model of infection indicated that MDSA could attenuate S. aureus virulence. This work is a rare demonstration of utilizing small molecules to block protein-DNA interaction, thus tuning important biological regulation at the transcriptional level.

DMAN: a Java tool for analysis of multi-well differential scanning fluorimetry experiments

SUMMARY

Differential scanning fluorimetry (DSF) is a rapid technique that can be used in structural biology to study protein-ligand interactions. We have developed DMAN, a novel tool to analyse multi-well plate data obtained in DSF experiments. DMAN is easy to install and provides a user-friendly interface. Multi-well plate layouts can be designed by the user and experimental data can be annotated and analysed by DMAN according to the specified plate layout. Statistical tests for significance are performed automatically, and graphical tools are also provided to assist in data analysis. The modular concept of this software will allow easy development of other multi-well plate analysis applications in the future.

AVAILABILITY AND IMPLEMENTATION

DMAN is implemented in Java to provide a cross-platform compatibility. It is freely available to academic users at http://www.structuralchemistry.org/pcsb/. To download DMAN, users will be asked for their name, institution and email address. A manual can also be downloaded from this site.

CONTACT

conan.wang@griffith.edu.au; a.hofmann@griffith.edu.au.

Screening, identification, and characterization of mechanistically diverse inhibitors of the Mycobacterium tuberculosis enzyme, pantothenate kinase (CoaA)

The authors describe the discovery of anti-mycobacterial compounds through identifying mechanistically diverse inhibitors of the essential Mycobacterium tuberculosis (Mtb) enzyme, pantothenate kinase (CoaA). Target-driven drug discovery technologies often work with purified enzymes, and inhibitors thus discovered may not optimally inhibit the form of the target enzyme predominant in the bacterial cell or may not be available at the desired concentration. Therefore, in addition to addressing entry or efflux issues, inhibitors with diverse mechanisms of inhibition (MoI) could be prioritized before hit-to-lead optimization. The authors describe a high-throughput assay based on protein thermal melting to screen large numbers of compounds for hits with diverse MoI. Following high-throughput screening for Mtb CoaA enzyme inhibitors, a concentration-dependent increase in protein thermal stability was used to identify true binders, and the degree of enhancement or reduction in thermal stability in the presence of substrate was used to classify inhibitors as competitive or non/uncompetitive. The thermal shift-based MoI assay could be adapted to screen hundreds of compounds in a single experiment as compared to traditional biochemical approaches for MoI determination. This MoI was confirmed through mechanistic studies that estimated K(ie) and K(ies) for representative compounds and through nuclear magnetic resonance-based ligand displacement assays.

Thermal denaturation assays in chemical biology

Thermal denaturation-based methods are becoming increasingly used to characterize protein stability and interactions. Recent technical advances have made these methods more suitable for high throughput screening. Reasonable throughput and the ability to perform these screens using commonly used instruments, such as RT-PCR machines or simple plate readers equipped with heating devices, facilitate these experiments in almost any laboratory. Introducing an aggregation-based monitoring approach as well as alternative fluorophores has allowed the screening of a wider range of proteins, including membrane proteins, against large chemical libraries. Thermal denaturation-based methods are independent of protein function, which is especially useful for the identification of orphan protein function. Here, we review applications of thermal denaturation-based methods in characterizing protein stability and ligand binding, and also provide information on protocol modifications that may further increase throughput.

GPCR stabilization using the bicelle-like architecture of mixed sterol-detergent micelles

The biophysical characterization of purified membrane proteins typically requires detergent mediated extraction from native lipid membrane environments. In the case of human G protein-coupled receptors (GPCRs), this process has been complicated by their conformational heterogeneity and the general lack of understanding the composition and interactions within the diverse human cellular membrane environment. Several successful GPCR structure determination efforts have shown that the addition of cholesterol analogs is often critical for maintaining protein stability. We have identified sterols that substantially increase the stability of the NOP receptor (ORL-1), a member of the opioid GPCR family, in a mixed micelle environment. Using dynamic light scattering and small-angle X-ray scattering, we have determined that the most thermal stabilizing sterol, cholesteryl hemisuccinate, induces the formation of a bicelle-like micelle architecture when mixed with dodecyl maltoside detergent. Together with mutagenesis studies and recent GPCR structures, our results provide indications that stabilization is attained through a combination of specific sterol binding to GPCRs and modulation of micelle morphology.

Protein-excipient interactions: mechanisms and biophysical characterization applied to protein formulation development

The purpose of this review is to demonstrate the critical importance of understanding protein-excipient interactions as a key step in the rational design of formulations to stabilize and deliver protein-based therapeutic drugs and vaccines. Biophysical methods used to examine various molecular interactions between solutes and protein molecules are discussed with an emphasis on applications to pharmaceutical excipients in terms of their effects on protein stability. Key mechanisms of protein-excipient interactions such as electrostatic and cation-pi interactions, preferential hydration, dispersive forces, and hydrogen bonding are presented in the context of different physical states of the formulation such as frozen liquids, solutions, gels, freeze-dried solids and interfacial phenomenon. An overview of the different classes of pharmaceutical excipients used to formulate and stabilize protein therapeutic drugs is also presented along with the rationale for use in different dosage forms including practical pharmaceutical considerations. The utility of high throughput analytical methodologies to examine protein-excipient interactions is presented in terms of expanding formulation design space and accelerating experimental timelines.

High-throughput measurement, correlation analysis, and machine-learning predictions for pH and thermal stabilities of Pfizer-generated antibodies

Generating stable antibodies is an important goal in the development of antibody-based drugs. Often, thermal stability is assumed predictive of overall stability. To test this, we used different internally created antibodies and first studied changes in antibody structure as a function of pH, using the dye ANS. Comparison of the pH(50) values, the midpoint of the transition from the high-pH to the low-pH conformation, allowed us for the first time to rank antibodies based on their pH stability. Next, thermal stability was probed by heating the protein in the presence of the dye Sypro Orange. A new data analysis method allowed extraction of all three antibody unfolding transitions and showed close correspondence to values obtained by differential scanning calorimetry. T(1%) , the temperature at which 1% of the protein is unfolded, was also determined. Importantly, no correlations could be found between thermal stability and pH(50) , suggesting that to accurately quantify antibody stability, different measures of protein stability are necessary. The experimental data were further analyzed using a machine-learning approach with a trained model that allowed the prediction of biophysical stability using primary sequence alone. The pH stability predictions proved most successful and were accurate to within pH ±0.2.

Structural and biochemical analysis of nuclease domain of clustered regularly interspaced short palindromic repeat (CRISPR)-associated protein 3 (Cas3)

RNA transcribed from clustered regularly interspaced short palindromic repeats (CRISPRs) protects many prokaryotes from invasion by foreign DNA such as viruses, conjugative plasmids, and transposable elements. Cas3 (CRISPR-associated protein 3) is essential for this CRISPR protection and is thought to mediate cleavage of the foreign DNA through its N-terminal histidine-aspartate (HD) domain. We report here the 1.8 Å crystal structure of the HD domain of Cas3 from Thermus thermophilus HB8. Structural and biochemical studies predict that this enzyme binds two metal ions at its active site. We also demonstrate that the single-stranded DNA endonuclease activity of this T. thermophilus domain is activated not by magnesium but by transition metal ions such as manganese and nickel. Structure-guided mutagenesis confirms the importance of the metal-binding residues for the nuclease activity and identifies other active site residues. Overall, these results provide a framework for understanding the role of Cas3 in the CRISPR system.

The combined use of the Thermofluor assay and ThermoQ analytical software for the determination of protein stability and buffer optimization as an aid in protein crystallization

The Thermofluor assay, also referred to as a thermal shift assay or differential scanning fluorescence (DSF), is a fast and simple method that is based upon high-throughput measurements of protein stability. The Thermofluor method can be performed on nearly all qPCR machines, meaning no special instrumentation is necessary, and can be used to validate the quality of protein preparations, screen for ligands or cofactors, and discover buffers and additives that maximize protein stability. This unit describes how to set up a Thermofluor method on several common models of qPCR instruments, how to prepare the samples for the assay, and how to run and analyze the resulting data. The unit also describes a 96-well screen to determine optimal buffer conditions for protein stability, which may assist in protein crystallization, and details the use of custom software, developed specifically for the analysis of data from Thermofluor screens.

Ligands of thermophilic ABC transporters encoded in a newly sequenced genomic region of Thermotoga maritima MSB8 screened by differential scanning fluorimetry

The chromosome of Thermotoga maritima strain MSB8 was found to have an 8,870-bp region that is not present in its published sequence. The isolate that was sequenced by The Institute for Genomic Research (TIGR) in 1999 is apparently a laboratory variant of the isolate deposited at the Deutsche Sammlung von Mikroorganismen und Zellkulturen (DSM 3109) in 1986. This newly sequenced region from the DSMZ culture was located between TM1848 (cbp, cellobiose phosphorylase) and TM1847 (the 3' end of a truncated ROK regulator). The new region contained seven genes: a beta glucosidase gene (bglA), three trehalose ABC transporter genes (treEFG), three xylose ABC transporter genes (xylE2F2K2), and the 5' end of a gene encoding the ROK regulator TM1847. We present a new differential scanning fluorimetry method using a low pH that was necessary to screen potential ligands of these exceptionally thermostable periplasmic substrate-binding proteins. This method showed that trehalose, sucrose, and glucose stabilized TreE, and their binding was confirmed by measuring changes in intrinsic fluorescence upon ligand binding. Binding constants of 0.024 μM, 0.300 μM, and 56.78 μM at 60°C, respectively, were measured. XylE2 ligands were similarly determined and xylose, glucose, and fucose bound with K(d) (dissociation constant) values of 0.042 μM, 0.059 μM, and 1.436 μM, respectively. Since there is no discernible phenotypic difference between the TIGR isolate and the DSMZ isolate despite the variance in their genomes, we propose that they be called genomovars: T. maritima MSB8 genomovar TIGR and T. maritima MSB8 genomovar DSM 3109, respectively.

Calcium modulates force sensing by the von Willebrand factor A2 domain

von Willebrand factor (VWF) multimers mediate primary adhesion and aggregation of platelets. VWF potency critically depends on multimer size, which is regulated by a feedback mechanism involving shear-induced unfolding of the VWF-A2 domain and cleavage by the metalloprotease ADAMTS-13. Here we report crystallographic and single-molecule optical tweezers data on VWF-A2 providing mechanistic insight into calcium-mediated stabilization of the native conformation that protects A2 from cleavage by ADAMTS-13. Unfolding of A2 requires higher forces when calcium is present and primarily proceeds through a mechanically stable intermediate with non-native calcium coordination. Calcium further accelerates refolding markedly, in particular, under applied load. We propose that calcium improves force sensing by allowing reversible force switching under physiologically relevant hydrodynamic conditions. Our data show for the first time the relevance of metal coordination for mechanical properties of a protein involved in mechanosensing.

Quantitation of protein-protein interactions by thermal stability shift analysis

Thermal stability shift analysis is a powerful method for examining binding interactions in proteins. We demonstrate that under certain circumstances, protein-protein interactions can be quantitated by monitoring shifts in thermal stability using thermodynamic models and data analysis methods presented in this work. This method relies on the determination of protein stabilities from thermal unfolding experiments using fluorescent dyes such as SYPRO Orange that report on protein denaturation. Data collection is rapid and straightforward using readily available real-time polymerase chain reaction instrumentation. We present an approach for the analysis of the unfolding transitions corresponding to each partner to extract the affinity of the interaction between the proteins. This method does not require the construction of a titration series that brackets the dissociation constant. In thermal shift experiments, protein stability data are obtained at different temperatures according to the affinity- and concentration-dependent shifts in unfolding transition midpoints. Treatment of the temperature dependence of affinity is, therefore, intrinsic to this method and is developed in this study. We used the interaction between maltose-binding protein (MBP) and a thermostable synthetic ankyrin repeat protein (Off7) as an experimental test case because their unfolding transitions overlap minimally. We found that MBP is significantly stabilized by Off7. High experimental throughput is enabled by sample parallelization, and the ability to extract quantitative binding information at a single partner concentration. In a single experiment, we were able to quantify the affinities of a series of alanine mutants, covering a wide range of affinities (∼ 100 nM to ∼ 100 μM).

The double-length tyrosyl-tRNA synthetase from the eukaryote Leishmania major forms an intrinsically asymmetric pseudo-dimer

The single tyrosyl-tRNA synthetase (TyrRS) gene in trypanosomatid genomes codes for a protein that is twice the length of TyrRS from virtually all other organisms. Each half of the double-length TyrRS contains a catalytic domain and an anticodon-binding domain; however, the two halves retain only 17% sequence identity to each other. The structural and functional consequences of this duplication and divergence are unclear. TyrRS normally forms a homodimer in which the active site of one monomer pairs with the anticodon-binding domain from the other. However, crystal structures of Leishmania major TyrRS show that, instead, the two halves of a single molecule form a pseudo-dimer resembling the canonical TyrRS dimer. Curiously, the C-terminal copy of the catalytic domain has lost the catalytically important HIGH and KMSKS motifs characteristic of class I aminoacyl-tRNA synthetases. Thus, the pseudo-dimer contains only one functional active site (contributed by the N-terminal half) and only one functional anticodon recognition site (contributed by the C-terminal half). Despite biochemical evidence for negative cooperativity between the two active sites of the usual TyrRS homodimer, previous structures have captured a crystallographically-imposed symmetric state. As the L. major TyrRS pseudo-dimer is inherently asymmetric, conformational variations observed near the active site may be relevant to understanding how the state of a single active site is communicated across the dimer interface. Furthermore, substantial differences between trypanosomal TyrRS and human homologs are promising for the design of inhibitors that selectively target the parasite enzyme.

Role of disulfide bridges in archaeal family-B DNA polymerases

The family-B DNA polymerases obtained from the order Thermococcales, for example, Pyrococcus furiosus (Pfu-Pol) are commonly used in the polymerase chain reaction (PCR) because of their high thermostability and low error rates. Most of these polymerases contain four cysteines, arranged as two disulfide bridges. With Pfu-Pol C429-C443 forms one of the disulfides (DB1) and C507-C510 (DB2) the other. Although the disulfides are well conserved in the enzymes from the hyperthermophilic Thermococcales, they are less prevalent in euryarchaeal polymerases from other orders, and tend to be only found in other hyperthermophiles. Here, we report on the effects of deleting the disulfide bridges by mutating the relevant cysteines to serines. A variety of techniques, including differential scanning calorimetry and differential scanning fluorimetry, have shown that both disulfides make a contribution to thermostability, with DB1 being more important than DB2. However, even when both disulfides are removed, sufficient thermostability remains for normal (identical to the wild type) performance in PCR and quantitative (real-time) PCR. Therefore, polymerases totally lacking cysteine are fully compatible with most PCR-based applications. This observation opens the way to further engineering of polymerases by introduction of a single cysteine followed by appropriate chemical modification.

Screening the MayBridge Rule of 3 Fragment Library for Compounds That Interact with the Trypanosoma brucei myo-Inositol-3-Phosphate Synthase and/or Show Trypanocidal Activity

Inositol-3-phosphate synthase (INO1) has previously been genetically validated as a drug target against Trypanosoma brucei, the causative agent of African sleeping sickness. Chemical intervention of this essential enzyme could lead to new therapeutic agents. Unfortunately, no potent inhibitors of INO1 from any organism have been reported, so a screen for potential novel inhibitors of T. brucei INO1was undertaken. Detection of inhibition of T. brucei INO1 is problematic due to the nature of the reaction. Direct detection requires differentiation between glucose-6-phosphate and inositol-3-phosphate. Coupled enzyme assays could give false positives as potentially they could inhibit the coupling enzyme. Thus, an alternative approach of differential scanning fluorimetry to identify compounds that interact with T. brucei INO1 was employed to screen ~670 compounds from the MayBridge Rule of 3 Fragment Library. This approach identified 38 compounds, which significantly altered the T_m of TbINO1. Four compounds showed trypanocidal activity with ED50s in the tens of micromolar range, with 2 having a selectivity index in excess of 250. The trypanocidal and general cytotoxicity activities of all of the compounds in the library are also reported, with the best having ED50S of ~20 μM against T. brucei.

Identification of inhibitors that block vaccinia virus infection by targeting the DNA synthesis processivity factor D4

Smallpox was globally eradicated 30 years ago by vaccination. The recent threat of bioterrorism demands the development of improved vaccines and novel therapeutics to effectively preclude a reemergence of smallpox. One new therapeutic target is the vaccinia poxvirus processivity complex, comprising D4 and A20 proteins that enable the viral E9 DNA polymerase to synthesize extended strands. Five compounds identified from an AlphaScreen assay designed to disrupt A20:D4 binding were shown to be effective in: (i) blocking vaccinia processive DNA synthesis in vitro, (ii) preventing cellular infection with minimal cytotoxicity, and (iii) binding to D4, as evidenced by ThermoFluor. The EC(50) values for inhibition of viral infectivity ranged from 9.6 to 23 μM with corresponding selectivity indices (cytotoxicity CC(50)/viral infectivity EC(50)) of 3.9 to 17.8. The five compounds are thus potential therapeutics capable of halting smallpox DNA synthesis and infectivity through disruptive action against a component of the vaccinia processivity complex.

Leveraging the contribution of thermodynamics in drug discovery with the help of fluorescence-based thermal shift assays

The development of new drugs with better pharmacological and safety properties mandates the optimization of several parameters. Today, potency is often used as the sole biochemical parameter to identify and select new molecules. Surprisingly, thermodynamics, which is at the core of any interaction, is rarely used in drug discovery, even though it has been suggested that the selection of scaffolds according to thermodynamic criteria may be a valuable strategy. This poor integration of thermodynamics in drug discovery might be due to difficulties in implementing calorimetry experiments despite recent technological progress in this area. In this report, the authors show that fluorescence-based thermal shift assays could be used as prescreening methods to identify compounds with different thermodynamic profiles. This approach allows a reduction in the number of compounds to be tested in calorimetry experiments, thus favoring greater integration of thermodynamics in drug discovery.

A novel mode of protein kinase inhibition exploiting hydrophobic motifs of autoinhibited kinases: discovery of ATP-independent inhibitors of fibroblast growth factor receptor

Protein kinase inhibitors with enhanced selectivity can be designed by optimizing binding interactions with less conserved inactive conformations because such inhibitors will be less likely to compete with ATP for binding and therefore may be less impacted by high intracellular concentrations of ATP. Analysis of the ATP-binding cleft in a number of inactive protein kinases, particularly in the autoinhibited conformation, led to the identification of a previously undisclosed non-polar region in this cleft. This ATP-incompatible hydrophobic region is distinct from the previously characterized hydrophobic allosteric back pocket, as well as the main pocket. Generalized hypothetical models of inactive kinases were constructed and, for the work described here, we selected the fibroblast growth factor receptor (FGFR) tyrosine kinase family as a case study. Initial optimization of a FGFR2 inhibitor identified from a library of commercial compounds was guided using structural information from the model. We describe the inhibitory characteristics of this compound in biophysical, biochemical, and cell-based assays, and have characterized the binding mode using x-ray crystallographic studies. The results demonstrate, as expected, that these inhibitors prevent activation of the autoinhibited conformation, retain full inhibitory potency in the presence of physiological concentrations of ATP, and have favorable inhibitory activity in cancer cells. Given the widespread regulation of kinases by autoinhibitory mechanisms, the approach described herein provides a new paradigm for the discovery of inhibitors by targeting inactive conformations of protein kinases.

An RNAi-based chemical genetic screen identifies three small-molecule inhibitors of the Wnt/wingless signaling pathway

Misregulated β-catenin responsive transcription (CRT) has been implicated in the genesis of various malignancies, including colorectal carcinomas, and it is a key therapeutic target in combating various cancers. Despite significant effort, successful clinical implementation of CRT inhibitory therapeutics remains a challenging goal. This is, in part, because of the challenge of identifying inhibitory compounds that specifically modulate the nuclear transcriptional activity of β-catenin while not affecting its cytoskeletal function in stabilizing adherens junctions at the cell membrane. Here, we report an RNAi-based modifier screening strategy for the identification of CRT inhibitors. Our data provide support for the specificity of these inhibitory compounds in antagonizing the transcriptional function of nuclear β-catenin. We show that these inhibitors efficiently block Wnt/β-catenin-induced target genes and phenotypes in various mammalian and cancer cell lines. Importantly, these Wnt inhibitors are specifically cytotoxic to human colon tumor biopsy cultures as well as colon cancer cell lines that exhibit deregulated Wnt signaling.

Sortase-catalyzed transformations that improve the properties of cytokines

Recombinant protein therapeutics often suffer from short circulating half-life and poor stability, necessitating multiple injections and resulting in limited shelf-life. Conjugation to polyethylene glycol chains (PEG) extends the circulatory half-life of many proteins, but the methods for attachment often lack specificity, resulting in loss of biological activity. Using four-helix bundle cytokines as an example, we present a general platform that uses sortase-mediated transpeptidation to facilitate site-specific attachment of PEG to extend cytokine half-life with full retention of biological activity. Covalently joining the N and C termini of proteins to obtain circular polypeptides, again executed using sortase, increases thermal stability. We combined both PEGylation and circularization by exploiting two distinct sortase enzymes and the use of a molecular suture that allows both site-specific PEGylation and covalent closure. The method developed is general, uses a set of easily accessible reagents, and should be applicable to a wide variety of proteins, provided that their termini are not involved in receptor binding or function.

Small-molecule inhibitors of FABP4/5 ameliorate dyslipidemia but not insulin resistance in mice with diet-induced obesity

Fatty acid binding protein-4 (FABP4) and FABP5 are two closely related FA binding proteins expressed primarily in adipose tissue and/or macrophages. The small-molecule FABP4 inhibitor BMS309403 was previously reported to improve insulin sensitivity in leptin-deficient Lep(ob)/Lep(ob) (ob/ob) mice. However, this compound was not extensively characterized in the more physiologically relevant animal model of mice with diet-induced obesity (DIO). Here, we report the discovery and characterization of a novel series of FABP4/5 dual inhibitors represented by Compounds 1-3. Compared with BMS309403, the compounds had significant in vitro potency toward both FABP4 and FABP5. In cell-based assays, Compounds 2 and 3 were more potent than BMS309403 to inhibit lipolysis in 3T3-L1 adipocytes and in primary human adipocytes. They also inhibited MCP-1 release from THP-1 macrophages as well as from primary human macrophages. When chronically administered to DIO mice, BMS309403 and Compound 3 reduced plasma triglyceride and free FA levels. Compound 3 reduced plasma free FAs at a lower dose level than BMS309403. However, no significant change was observed in insulin, glucose, or glucose tolerance. Our results indicate that the FABP4/5 inhibitors ameliorate dyslipidemia but not insulin resistance in DIO mice.

The polypeptide binding conformation of calreticulin facilitates its cell-surface expression under conditions of endoplasmic reticulum stress

We define two classes of calreticulin mutants that retain glycan binding activity; those that display enhanced or reduced polypeptide-specific chaperone activity, due to conformational effects. Under normal conditions, neither set of mutants significantly impacts the ability of calreticulin to mediate assembly and trafficking of major histocompatibility complex class I molecules, which are calreticulin substrates. However, in cells treated with thapsigargin, which depletes endoplasmic reticulum calcium, major histocompatibility complex class I trafficking rates are accelerated coincident with calreticulin secretion, and detection of cell-surface calreticulin is dependent on its polypeptide binding conformations. Together, these findings identify a site on calreticulin that is an important determinant of the induction of its polypeptide binding conformation and demonstrate the relevance of the polypeptide binding conformations of calreticulin to endoplasmic reticulum stress-induced interactions.

Crystal structures of three protozoan homologs of tryptophanyl-tRNA synthetase

Tryptophanyl-tRNA synthetase (TrpRS) is an essential enzyme that is recognizably conserved across all forms of life. It is responsible for activating and attaching tryptophan to a cognate tRNA(Trp) molecule for use in protein synthesis. In some eukaryotes this original core function has been supplemented or modified through the addition of extra domains or the expression of variant TrpRS isoforms. The three TrpRS structures from pathogenic protozoa described here represent three illustrations of this malleability in eukaryotes. The Cryptosporidium parvum genome contains a single TrpRS gene, which codes for an N-terminal domain of uncertain function in addition to the conserved core TrpRS domains. Sequence analysis indicates that this extra domain, conserved among several apicomplexans, is related to the editing domain of some AlaRS and ThrRS. The C. parvum enzyme remains fully active in charging tRNA(Trp) after truncation of this extra domain. The crystal structure of the active, truncated enzyme is presented here at 2.4Å resolution. The Trypanosoma brucei genome contains separate cytosolic and mitochondrial isoforms of TrpRS that have diverged in their respective tRNA recognition domains. The crystal structure of the T. brucei cytosolic isoform is presented here at 2.8Å resolution. The Entamoeba histolytica genome contains three sequences that appear to be TrpRS homologs. However one of these, whose structure is presented here at 3.0Å resolution, has lost the active site motifs characteristic of the Class I aminoacyl-tRNA synthetase catalytic domain while retaining the conserved features of a fully formed tRNA(Trp) recognition domain. The biological function of this variant E. histolytica TrpRS remains unknown, but, on the basis of a completely conserved tRNA recognition region and evidence for ATP but not tryptophan binding, it is tempting to speculate that it may perform an editing function. Together with a previously reported structure of an unusual TrpRS from Giardia, these protozoan structures broaden our perspective on the extent of structural variation found in eukaryotic TrpRS homologs.

The TB Structural Genomics Consortium: a decade of progress

The TB Structural Genomics Consortium is a worldwide organization of collaborators whose mission is the comprehensive structural determination and analyses of Mycobacterium tuberculosis proteins to ultimately aid in tuberculosis diagnosis and treatment. Congruent to the overall vision, Consortium members have additionally established an integrated facilities core to streamline M. tuberculosis structural biology and developed bioinformatics resources for data mining. This review aims to share the latest Consortium developments with the TB community, including recent structures of proteins that play significant roles within M. tuberculosis. Atomic resolution details may unravel mechanistic insights and reveal unique and novel protein features, as well as important protein-protein and protein-ligand interactions, which ultimately lead to a better understanding of M. tuberculosis biology and may be exploited for rational, structure-based therapeutics design.

Oligomerization of EpsE coordinates residues from multiple subunits to facilitate ATPase activity

EpsE is an ATPase that powers transport of cholera toxin and hydrolytic enzymes through the Type II secretion (T2S) apparatus in the gram-negative bacterium, Vibrio cholerae. On the basis of structures of homologous Type II/IV secretion ATPases and our biochemical data, we believe that EpsE is active as an oligomer, likely a hexamer, and the binding, hydrolysis, and release of nucleotide cause EpsE to undergo dynamic structural changes, thus converting chemical energy to mechanical work, ultimately resulting in extracellular secretion. The conformational changes that occur as a consequence of nucleotide binding would realign conserved arginines (Arg(210), Arg(225), Arg(320), Arg(324), Arg(336), and Arg(369)) from adjoining domains and subunits to complete the active site around the bound nucleotide. Our data suggest that these arginines are essential for ATP hydrolysis, although their roles in shaping the active site of EpsE are varied. Specifically, we have shown that replacements of these arginine residues abrogate the T2S process due to a reduction of ATPase activity yet do not have any measurable effect on nucleotide binding or oligomerization of EpsE. We have further demonstrated that point mutations in the EpsE intersubunit interface also reduce ATPase activity without disrupting oligomerization, strengthening the idea that residues from multiple subunits must precisely interact in order for EpsE to be sufficiently active to support T2S. Our findings suggest that the action of EpsE is similar to that of other Type II/IV secretion ATPase family members, and thus these results may be widely applicable to the family as a whole.

Combining biophysical screening and X-ray crystallography for fragment-based drug discovery

Over the past decade, fragment-based drug discovery (FBDD) has gained importance for the generation of novel ideas to inspire synthetic chemistry. In order to identify small molecules that bind to a target protein, multiple approaches have been utilized by various groups in the pharmaceutical industry and by academic groups. The combination of fragment screening by biophysical methods and in particular with surface plasmon resonance technologies (SPR) together with the visualization of the binding properties by X-ray crystallography offers a number of benefits. Screening by SPR identifies ligands for a target protein as well as provides an assessment of the binding properties with respect to affinity, stoichiometry, and specificity of the interaction. Despite the huge technology advances of the past years, X-ray crystallography is still a resource-intensive technology, and SPR binding data provides excellent measures to prioritize X-ray experiments and consequently enable a better success rate in obtaining structural information. Information on the chemical structures of fragments binding to a protein can be used to perform similarity searches in compound libraries in order to establish structure-activity relationships as well as to explore particular scaffolds. At Roche we have applied this workflow for a number of targets and the experiences will be outlined in this review.

Protein thermal shifts to identify low molecular weight fragments

Measuring the strength of binding of low molecular weight ligands to a target protein is a significant challenge to fragment-based drug discovery that must be solved. Thermal shift assays are uniquely suited for this purpose, due to the thermodynamic effects of a ligand on protein thermal stability. We show here how to implement a thermal shift assay, describing the basic features and analysis of the protein unfolding data. We then describe in detail the effects of a ligand on the observed stability of the protein to produce a shift in stability. The anatomy of ligand-induced thermal shift data is discussed in detail. We describe the unique aspects of concentration-response curves, the effect of protein unfolding energetics, and the stoichiometry of the interaction. We outline a typical assay development strategy for optimizing dye type and concentration, protein concentration, and buffer conditions. Guidelines are presented to demonstrate the limits of detection for weak-binding ligands, as applied to sulfonamide-based inhibitors of carbonic anhydrase II and applied to nucleotide binding to the death-associated protein kinase 1 catalytic domain.

The use of flow cytometry for the detection of subvisible particles in therapeutic protein formulations

The amount, identity, and size distribution of particles in parenteral therapeutic protein formulations are of immense interest due to potential safety and efficacy-related implications. In this communication, we describe the use of a flow cytometer equipped with forward- and side-scattering as well as fluorescence detectors, to determine the number of subvisible particles in monoclonal antibody formulations. The method appears to detect particles of size 1 μ and larger, requiring relatively small sample volumes to estimate subvisible particle counts. Additionally, it facilitates differentiation of proteinaceous particles after staining with a fluorescent hydrophobic dye. The method is expected to be particularly well suited for pharmaceutical development, because it provides increased throughput due to the use of a 96-well autosampler.

Thermodynamic analysis of ligand-induced changes in protein thermal unfolding applied to high-throughput determination of ligand affinities with extrinsic fluorescent dyes

The quantification of protein-ligand interactions is essential for systems biology, drug discovery, and bioengineering. Ligand-induced changes in protein thermal stability provide a general, quantifiable signature of binding and may be monitored with dyes such as Sypro Orange (SO), which increase their fluorescence emission intensities upon interaction with the unfolded protein. This method is an experimentally straightforward, economical, and high-throughput approach for observing thermal melts using commonly available real-time polymerase chain reaction instrumentation. However, quantitative analysis requires careful consideration of the dye-mediated reporting mechanism and the underlying thermodynamic model. We determine affinity constants by analysis of ligand-mediated shifts in melting-temperature midpoint values. Ligand affinity is determined in a ligand titration series from shifts in free energies of stability at a common reference temperature. Thermodynamic parameters are obtained by fitting the inverse first derivative of the experimental signal reporting on thermal denaturation with equations that incorporate linear or nonlinear baseline models. We apply these methods to fit protein melts monitored with SO that exhibit prominent nonlinear post-transition baselines. SO can perturb the equilibria on which it is reporting. We analyze cases in which the ligand binds to both the native and denatured state or to the native state only and cases in which protein:ligand stoichiometry needs to treated explicitly.

A novel p53 phosphorylation site within the MDM2 ubiquitination signal: II. a model in which phosphorylation at SER269 induces a mutant conformation to p53

The p53 DNA-binding domain harbors a conformationally flexible multiprotein binding site that regulates p53 ubiquitination. A novel phosphorylation site exists within this region at Ser(269), whose phosphomimetic mutation inactivates p53. The phosphomimetic p53 (S269D) exhibits characteristics of mutant p53: stable binding to Hsp70 in vivo, elevated ubiquitination in vivo, inactivity in DNA binding and transcription, increased thermoinstability using thermal shift assays, and λ(max) of intrinsic tryptophan fluorescence at 403 nm rather than 346 nm, characteristic of wild type p53. These data indicate that p53 conformational stability is regulated by a phosphoacceptor site within an exposed flexible surface loop and that this can be destabilized by phosphorylation. To test whether other motifs within p53 have similarly evolved, we analyzed the effect of Ser(215) mutation on p53 function because Ser(215) is another inactivating phosphorylation site in the conformationally flexible PAb240 epitope. The p53(S215D) protein is inactive like p53(S269D), whereas p53(S215A) is as active as p53(S269A). However, the double mutant p53(S215A/S269A) was transcriptionally inactive and more thermally unstable than either individual Ser-Ala loop mutant. Molecular dynamics simulations suggest that (i) solvation of phospho-Ser(215) and phospho-Ser(269) by positive charged residues or solvent water leads to local unfolding, which is accompanied by local destabilization of the N-terminal loop and global destabilization of p53, and (ii) the double alanine 215/269 mutation disrupts hydrogen bonding normally stabilized by both Ser(215) and Ser(269). These data indicate that p53 has evolved two serine phosphoacceptor residues within conformationally flexible epitopes that normally stabilize the p53 DNA-binding domain but whose phosphorylation induces a mutant conformation to wild type p53.

High-affinity inhibitors of human NAD-dependent 15-hydroxyprostaglandin dehydrogenase: mechanisms of inhibition and structure-activity relationships

Background

15-hydroxyprostaglandin dehydrogenase (15-PGDH, EC 1.1.1.141) is the key enzyme for the inactivation of prostaglandins, regulating processes such as inflammation or proliferation. The anabolic pathways of prostaglandins, especially with respect to regulation of the cyclooxygenase (COX) enzymes have been studied in detail; however, little is known about downstream events including functional interaction of prostaglandin-processing and -metabolizing enzymes. High-affinity probes for 15-PGDH will, therefore, represent important tools for further studies.

Principal Findings

To identify novel high-affinity inhibitors of 15-PGDH we performed a quantitative high-throughput screen (qHTS) by testing >160 thousand compounds in a concentration-response format and identified compounds that act as noncompetitive inhibitors as well as a competitive inhibitor, with nanomolar affinity. Both types of inhibitors caused strong thermal stabilization of the enzyme, with cofactor dependencies correlating with their mechanism of action. We solved the structure of human 15-PGDH and explored the binding modes of the inhibitors to the enzyme in silico. We found binding modes that are consistent with the observed mechanisms of action.

Conclusions

Low cross-reactivity in screens of over 320 targets, including three other human dehydrogenases/reductases, suggest selectivity of the present inhibitors for 15-PGDH. The high potencies and different mechanisms of action of these chemotypes make them a useful set of complementary chemical probes for functional studies of prostaglandin-signaling pathways.

Enhanced version

This article can also be viewed as an enhanced version in which the text of the article is integrated with interactive 3D representations and animated transitions. Please note that a web plugin is required to access this enhanced functionality. Instructions for the installation and use of the web plugin are available in Text S2.