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

The fluorescence-based thermal shift assay is a general method for identification of inhibitors of target proteins from compound libraries. Using an environmentally sensitive fluorescent dye to monitor protein thermal unfolding, the ligand-binding affinity can be assessed from the shift of the unfolding temperature (Delta Tm) obtained in the presence of ligands relative to that obtained in the absence of ligands. In this article, we report that the thermal shift assay can be conducted in an inexpensive, commercially available device for temperature control and fluorescence detection. The binding affinities obtained from thermal shift assays are compared with the binding affinities measured by isothermal titration calorimetry and with the IC(50) values from enzymatic assays. The potential pitfalls in the data analysis of thermal shift assays are also discussed.     

Ligand Binding Mechanics of Maltose Binding Protein

In the past decade, single-molecule force spectroscopy has provided new insights into the key interactions stabilizing folded proteins. A few recent studies probing the effects of ligand binding on mechanical protein stability have come to quite different conclusions. While some proteins seem to be stabilized considerably by a bound ligand, others appear to be unaffected. Since force acts as a vector in space, it is conceivable that mechanical stabilization by ligand binding is dependent on the direction of force application. In this study, we vary the direction of the force to investigate the effect of ligand binding on the stability of maltose binding protein (MBP). MBP consists of two lobes connected by a hinge region that move from an open to a closed conformation when the ligand maltose binds. Previous mechanical experiments, where load was applied to the N and C termini, have demonstrated that MBP is built up of four building blocks (unfoldons) that sequentially detach from the folded structure. In this study, we design the pulling direction so that force application moves the two MBP lobes apart along the hinge axis. Mechanical unfolding in this geometry proceeds via an intermediate state whose boundaries coincide with previously reported MBP unfoldons. We find that in contrast to N-C-terminal pulling experiments, the mechanical stability of MBP is increased by ligand binding when load is applied to the two lobes and force breaks the protein-ligand interactions directly. Contour length measurements indicate that MBP is forced into an open conformation before unfolding even if ligand is bound. Using mutagenesis experiments, we demonstrate that the mechanical stabilization effect is due to only a few key interactions of the protein with its ligand. This work illustrates how varying the direction of the applied force allows revealing important details about the ligand binding mechanics of a large protein.     

Integration of cell-free protein coexpression with an enzyme-linked immunosorbent assay enables rapid analysis of protein-protein interactions directly from DNA

Assays that integrate detection of binding with cell-free protein expression directly from DNA can dramatically increase the pace at which protein-protein interactions (PPIs) can be analyzed by mutagenesis. In this study, we present a method that combines in vitro protein production with an enzyme-linked immunosorbent assay (ELISA) to measure PPIs. This method uses readily available commodity instrumentation and generic antibody-affinity tag interactions. It is straightforward and rapid to execute, enabling many interactions to be assessed in parallel. In traditional ELISAs, reporter complexes are assembled stepwise with one layer at a time. In the method presented here, all the members of the reporter complex are present and assembled together. The signal strength is dependent on all the intercomponent interaction affinities and concentrations. Although this assay is straightforward to execute, establishing proper conditions and analysis of the results require a thorough understanding of the processes that determine the signal strength. The formation of the fully assembled reporter sandwich can be modeled as a competition between Langmuir adsorption isotherms for the immobilized components and binding equilibria of the solution components. We have shown that modeling this process provides semiquantitative understanding of the effects of affinity and concentration and can guide strategies for the development of experimental protocols. We tested the method experimentally using the interaction between a synthetic ankyrin repeat protein (Off7) and maltose-binding protein. Measurements obtained for a collection of alanine mutations in the interface between these two proteins demonstrate that a range of affinities can be analyzed.     

Effect of signal peptide on the stability and folding kinetics of maltose binding protein

While the role of the signal sequence in targeting proteins to specific subcellular compartments is well characterized, there are fewer studies that characterize its effects on the stability and folding kinetics of the protein. We report a detailed characterization of the folding kinetics and thermodynamic stabilities of maltose binding protein (MBP) and its precursor form, preMBP. Isothermal GdmCl and urea denaturation as a function of temperature and thermal denaturation studies have been carried out to compare stabilities of the two proteins. preMBP was found to be destabilized by about 2-6 kcal/mol (20-40%) with respect to MBP. Rapid cleavage of the signal peptide by various proteases shows that the signal peptide is accessible in the native form of preMBP. The observed rate constant of the major slow phase in folding was decreased 5-fold in preMBP relative to MBP. The rate constants of unfolding were similar at 25 degrees C, but preMBP also exhibited a large burst phase change in unfolding that was absent in MBP. At 10 degrees C, preMBP exhibited a higher unfolding rate than MBP as well as a large burst phase. The appreciable destabilization of MBP by signal peptide is functionally relevant, because it enhances the likelihood of finding the protein in an unfolded translocation-competent form and may influence the interactions of the protein with the translocation machinery. Destabilization is likely to result from favorable interactions between the hydrophobic signal peptide and other hydrophobic regions that are exposed in the unfolded state.     

Investigating protein unfolding kinetics by pulse proteolysis

Investigation of protein unfolding kinetics of proteins in crude samples may provide many exciting opportunities to study protein energetics under unconventional conditions. As an effort to develop a method with this capability, we employed “pulse proteolysis” to investigate protein unfolding kinetics. Pulse proteolysis has been shown to be an effective and facile method to determine global stability of proteins by exploiting the difference in proteolytic susceptibilities between folded and unfolded proteins. Electrophoretic separation after proteolysis allows monitoring protein unfolding without protein purification. We employed pulse proteolysis to determine unfolding kinetics of E. coli maltose binding protein (MBP) and E. coli ribonuclease H (RNase H). The unfolding kinetic constants determined by pulse proteolysis are in good agreement with those determined by circular dichroism. We then determined an unfolding kinetic constant of overexpressed MBP in a cell lysate. An accurate unfolding kinetic constant was successfully determined with the unpurified MBP. Also, we investigated the effect of ligand binding on unfolding kinetics of MBP using pulse proteolysis. On the basis of a kinetic model for unfolding of MBP•maltose complex, we have determined the dissociation equilibrium constant (K_d) of the complex from unfolding kinetic constants, which is also in good agreement with known K_d values of the complex. These results clearly demonstrate the feasibility and the accuracy of pulse proteolysis as a quantitative probe to investigate protein unfolding kinetics.     

How aggregation and conformational scrambling of unfolded states govern fluorescence emission spectra

In a case study on five homologous alpha-amylases we analyzed the properties of unfolded states as obtained from treatments with GndHCl and with elevated temperatures. In particular the wavelength of the tryptophan fluorescence emission peak (lambda(max)) is a valuable parameter to characterize properties of the unfolded state. In all cases with a typical red shift of the emission spectrum occurring during structural unfolding we observed a larger magnitude of this shift for GndHCl-induced unfolding as compared to thermal unfolding. Although a quantitative relation between aggregation and reduction of the unfolding induced red shifts cannot be given, our data indicate that protein aggregation contributes significantly to smaller magnitudes of red shifts as observed during thermal unfolding. In addition, other properties of the unfolded states, most probable structural compactness or simply differences in the conformational scrambling, also affect the magnitude of red shifts. For the irreversible unfolding alpha-amylases studied here, transition temperatures and magnitudes of red shifts are strongly depending on heating rates. Lower protein concentrations and smaller heating rates lead to larger red shifts upon thermal unfolding, indicating that under these conditions the protein aggregation is less pronounced.     

Native-specific stabilization of flavodoxin by the FMN cofactor: structural and thermodynamical explanation

Flavodoxins are useful models to investigate protein/cofactor interactions. The binding energy of the apoflavodoxin-FMN complex is high and therefore the holoflavodoxin is expected to be more stable than the apoprotein. This expectation has been challenged by reports on the stability of Desulfovibrio desulfuricans flavodoxin indicating that FMN binds to the unfolded polypeptide with similar affinity as to the native state, thus causing no net effect on protein stability. In previous work, we have analyzed in detail the stability of the apoflavodoxin from Anabaena PCC 7119 and the energetics of its functional complex with FMN. Here, we use the Anabaena holoprotein to directly investigate the contribution of the bound cofactor to protein stability through a detailed analysis of the chemical and thermal denaturation equilibria. Our data clearly shows that FMN binding largely stabilizes the protein towards both chemical and thermal denaturation, and that the stabilization observed at 25 degrees C in low ionic strength conditions is precisely the one expected if full release of the cofactor takes place upon flavodoxin unfolding. On the other hand, the binding of FMN to the native polypeptide is shown to simplify the thermal unfolding so that, while apoflavodoxin follows a three-state mechanism, the holoprotein unfolds in a two-state fashion. Comparison of the X-ray structure of native apoflavodoxin with the phi-structure of the thermal intermediate indicates that the increase in cooperativity driven by the cofactor originates in its preferential binding to the native state, which is a consequence of the disorganization in the intermediate of the FMN binding loops and of an adjacent longer loop.     

Computational protein design suggests that human PCNA‐partner interactions are not optimized for affinity

Increasing the affinity of binding proteins is invaluable for basic and applied biological research. Currently, directed protein evolution experiments are the main approach for generating such proteins through the construction and screening of large mutant libraries. Proliferating cell nuclear antigen (PCNA) is an essential hub protein that interacts with many different partners to tightly regulate DNA replication and repair in all eukaryotes. Here, we used computational design to generate human PCNA mutants with enhanced affinity for several different partners. We identified double mutations in PCNA, outside the main partner binding site, that were predicted to increase PCNA‐partner binding affinities compared to the wild‐type protein by forming additional hydrophobic interactions with conserved residues in the PCNA partners. Affinity increases were experimentally validated with four different PCNA partners, demonstrating that computational design can reveal unexpected regions where affinity enhancements in natural systems are possible. The designed PCNA mutants can be used as a valuable tool for further examination of the regulation of PCNA‐partner interactions during DNA replication and repair both in vitro and in vivo. More broadly, the ability to engineer affinity increases toward several PCNA partners suggests that interaction affinity is not an evolutionarily optimized trait of this system. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

A miniaturized technique for assessing protein thermodynamics and function using fast determination of quantitative cysteine reactivity

Protein thermodynamic stability is a fundamental physical characteristic that determines biological function. Furthermore, alteration of thermodynamic stability by macromolecular interactions or biochemical modifications is a powerful tool for assessing the relationship between protein structure, stability, and biological function. High-throughput approaches for quantifying protein stability are beginning to emerge that enable thermodynamic measurements on small amounts of material, in short periods of time, and using readily accessible instrumentation. Here we present such a method, fast quantitative cysteine reactivity, which exploits the linkage between protein stability, sidechain protection by protein structure, and structural dynamics to characterize the thermodynamic and kinetic properties of proteins. In this approach, the reaction of a protected cysteine and thiol-reactive fluorogenic indicator is monitored over a gradient of temperatures after a short incubation time. These labeling data can be used to determine the midpoint of thermal unfolding, measure the temperature dependence of protein stability, quantify ligand-binding affinity, and, under certain conditions, estimate folding rate constants. Here, we demonstrate the fQCR method by characterizing these thermodynamic and kinetic properties for variants of Staphylococcal nuclease and E. coli ribose-binding protein engineered to contain single, protected cysteines. These straightforward, information-rich experiments are likely to find applications in protein engineering and functional genomics.     

Structural,thermodynamic, and phosphatidylinositol 3‐phosphate binding properties of Phafin2

Phafin2 is a phosphatidylinositol 3‐phosphate (PtdIns(3)P) binding protein involved in the regulation of endosomal cargo trafficking and lysosomal induction of autophagy. Binding of Phafin2 to PtdIns(3)P is mediated by both its PH and FYVE domains. However, there are no studies on the structural basis, conformational stability, and lipid interactions of Phafin2 to better understand how this protein participates in signaling at the surface of endomembrane compartments. Here, we show that human Phafin2 is a moderately elongated monomer of ～28 kDa with an intensity‐average hydrodynamic diameter of ～7 nm. Circular dichroism (CD) analysis indicates that Phafin2 exhibits an α/β structure and predicts ～40% random coil content in the protein. Heteronuclear NMR data indicates that a unique conformation of Phafin2 is present in solution and dispersion of resonances suggests that the protein exhibits random coiled regions, in agreement with the CD data. Phafin2 is stable, displaying a melting temperature of 48.4°C. The folding‐unfolding curves, obtained using urea‐ and guanidine hydrochloride‐mediated denaturation, indicate that Phafin2 undergoes a two‐state native‐to‐denatured transition. Analysis of these transitions shows that the free energy change for urea‐ and guanidine hydrochloride‐induced Phafin2 denaturation in water is ～4 kcal mol^?1. PtdIns(3)P binding to Phafin2 occurs with high affinity, triggering minor conformational changes in the protein. Taken together, these studies represent a platform for establishing the structural basis of Phafin2 molecular interactions and the role of the two potentially redundant PtdIns(3)P‐binding domains of the protein in endomembrane compartments.     

Thermal versus guanidine-induced unfolding of ubiquitin. An analysis in terms of the contributions from charge-charge interactions to protein stability.

We have characterized the guanidine-induced unfolding of both yeast and bovine ubiquitin at 25 degrees C and in the acidic pH range on the basis of fluorescence and circular dichroism measurements. Unfolding Gibbs energy changes calculated by linear extrapolation from high guanidine unfolding data are found to depend very weakly on pH. A simple explanation for this result involves the two following assumptions: (1) charged atoms of ionizable groups are exposed to the solvent in native ubiquitin (as supported by accessible surface area calculations), and Gibbs energy contributions associated with charge desolvation upon folding (a source of pK shifts) are small; (2) charge-charge interactions (another source of pK shifts upon folding) are screened out in concentrated guanidinium chloride solutions. We have also characterized the thermal unfolding of both proteins using differential scanning calorimetry. Unfolding Gibbs energy changes calculated from the calorimetric data do depend strongly on pH, a result that we attribute to the pH dependence of charge-charge interactions (not eliminated in the absence of guanidine). In fact, we find good agreement between the difference between the two series of experimental unfolding Gibbs energy changes (determined from high guanidine unfolding data by linear extrapolation and from thermal denaturation data in the absence of guanidine) and the theoretical estimates of the contribution from charge-charge interactions to the Gibbs energy change for ubiquitin unfolding obtained by using the solvent-accessibility-corrected Tanford-Kirkwood model, together with the Bashford-Karplus (reduced-set-of-sites) approximation. This contribution is found to be stabilizing at neutral pH, because most charged groups on the native protein interact mainly with groups of the opposite charge, a fact that, together with the absence of large charge-desolvation contributions, may explain the high stability of ubiquitin at neutral pH. In general, our analysis suggests the possibility of enhancing protein thermal stability by adequately redesigning the distribution of solvent-exposed, charged residues on the native protein surface.     

A novel method to determine thermal transition curves of disulfide-containing proteins and their structured folding intermediates

The stability of a protein or of its folding intermediates is frequently characterized by its resistance to chemical and/or thermal denaturation. The folding/unfolding process is generally followed by spectroscopic methods such as absorbance, fluorescence, circular dichroism spectroscopy, etc. Here, we demonstrate a new method, by using HPLC, for determining the thermal unfolding transitions of disulfide-containing proteins and their structured folding intermediates. The thermal transitions of a model protein, ribonuclease A (RNase A), and a recently found unfolding intermediate of onconase (ONC), des [30-75], have been estimated by this method. Finally, the advantages of this method over traditional techniques are discussed by providing specific examples.     

Stability of Ligand-induced Protein Conformation Influences Affinity in Maltose-binding Protein

Our understanding of what determines ligand affinity of proteins is poor, even with high-resolution structures available. Both the non-covalent ligand–protein interactions and the relative free energies of available conformations contribute to the affinity of a protein for a ligand. Distant, non-binding site residues can influence the ligand affinity by altering the free energy difference between a ligand-free and ligand-bound conformation. Our hypothesis is that when different ligands induce distinct ligand-bound conformations, it should be possible to tweak their affinities by changing the free energies of the available conformations. We tested this idea for the maltose-binding protein (MBP) from Escherichia coli. We used single-molecule Förster resonance energy transfer (smFRET) to distinguish several unique ligand-bound conformations of MBP. We engineered mutations, distant from the binding site, to affect the stabilities of different ligand-bound conformations. We show that ligand affinity can indeed be altered in a conformation-dependent manner. Our studies provide a framework for the tuning of ligand affinity, apart from modifying binding site residues.     

Evolutionary stabilization of the gene-3-protein of phage fd reveals the principles that govern the thermodynamic stability of two-domain proteins

The gene-3-protein (G3P) of filamentous phage is essential for their propagation. It consists of three domains. The CT domain anchors G3P in the phage coat, the N2 domain binds to the F pilus of Escherichia coli and thus initiates infection, and the N1 domain continues by interacting with the TolA receptor. Phage are thus only infective when the three domains of G3P are tightly linked, and this requirement is exploited by Proside, an in vitro selection method for proteins with increased stability. In Proside, a repertoire of variants of the protein to be stabilized is inserted between the N2 and the CT domains of G3P. Stabilized variants can be selected because they resist cleavage by a protease and thus maintain the essential linkage between the domains. The method is limited by the proteolytic stability of G3P itself. We improved the stability of G3P by subjecting the phage without a guest protein to rounds of random in vivo mutagenesis and proteolytic Proside selections. Variants of G3P with one to four mutations were selected, and the temperature at which the corresponding phage became accessible for a protease increased in a stepwise manner from 40 degrees C to almost 60 degrees C. The N1-N2 fragments of wild-type gene-3-protein and of the four selected variants were purified and their stabilities towards thermal and denaturant-induced unfolding were determined. In the biphasic transitions of these proteins domain dissociation and unfolding of N2 occur in a concerted reaction in the first step, followed by the independent unfolding of domain N1 in the second step. N2 is thus less stable than N1, and it unfolds when the interactions with N1 are broken. The strongest stabilizations were caused by mutations in domain N2, in particular in its hinge subdomain, which provides many stabilizing interactions between the N1 and N2 domains. These results reveal how the individual domains and their assembly contribute to the overall stability of two-domain proteins and how mutations are optimally placed to improve the stability of such proteins.     

Mechanical unfoldons as building blocks of maltose-binding protein

Identifying independently folding cores or substructures is important for understanding and assaying the structure, function and assembly of large proteins. Here, we suggest mechanical stability as a criterion to identify building blocks of the 366 amino acid maltose-binding protein (MBP). We find that MBP, when pulled at its termini, unfolds via three (meta-) stable unfolding intermediates. Consequently, the MBP structure consists of four structural blocks (unfoldons) that detach sequentially from the folded structure upon force application. We used cysteine cross-link mutations to characterize the four unfoldons structurally. We showed that many MBP constructs composed of those building blocks indeed form stably folded structures in solution. Mechanical unfoldons may provide a new tool for a systematic search for stable substructures of large proteins.     

Ligand binding analysis and screening by chemical denaturation shift

The identification of small molecule ligands is an important first step in drug development, especially drugs that target proteins with no intrinsic activity. Toward this goal, it is important to have access to technologies that are able to measure binding affinities for a large number of potential ligands in a fast and accurate way. Because ligand binding stabilizes the protein structure in a manner dependent on concentration and binding affinity, the magnitude of the protein stabilization effect elicited by binding can be used to identify and characterize ligands. For example, the shift in protein denaturation temperature (T_m shift) has become a popular approach to identify potential ligands. However, T_m shifts cannot be readily transformed into binding affinities, and the ligand rank order obtained at denaturation temperatures (?60 °C) does not necessarily coincide with the rank order at physiological temperature. An alternative approach is the use of chemical denaturation, which can be implemented at any temperature. Chemical denaturation shifts allow accurate determination of binding affinities with a surprisingly wide dynamic range (high micromolar to sub nanomolar) and in situations where binding changes the cooperativity of the unfolding transition. In this article, we develop the basic analytical equations and provide several experimental examples.     

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.     

Stability of a guest protein depends on stability of a host protein in insertional fusion

Insertional fusion between host and guest protein domains has been employed to create multi-domain protein complexes displaying integrated and coupled functionalities. The effects of insertional fusion on the stability of a guest protein are however rather controversial. In the study described here, we examined whether the stability of inserted TEM1 beta-lactamase (BLA), as a guest protein, might be affected by the stability of a maltodextrin-binding protein (MBP), as a host protein. Our results indicate that expression levels and in vitro stability of the BLA domain were significantly higher when inserted into thermophilic MBP from Pyrococcus furiosus (PfMBP) compared to mesophilic MBP from Escherichia coli (EcMBP). Moreover, insertion into PfMBP at selected sites was found to improve thermal stability of the BLA domain without compromise in expression levels and BLA activity. Kinetic stabilization during prolonged thermal denaturation of the BLA domain was not guaranteed by insertion into PfMBP, but rather relied on the insertion sites. Taken together, we provide evidence that (i) the stability of the guest protein depended on the stability of the host protein in insertional fusion and (ii) insertion into PfMBP, at least at selected locations, can serve as a novel method of improving protein thermal stability.     

Thermal unfolding of eosinophil cationic protein/ribonuclease 3: a nonreversible process

Eosinophil cationic protein (ECP)/ribonuclease 3 is a member of the RNase A superfamily involved in inflammatory processes mediated by eosinophils. ECP is bactericidal, helminthotoxic, and cytotoxic to tracheal epithelium cells and to several mammalian cell lines although its RNase activity is low. We studied the thermal stability of ECP by fourth-derivative UV absorbance spectra, circular dichroism, differential scanning calorimetry, and Fourier transform infrared spectroscopy. The T (1/2) values obtained with the different techniques were in very good agreement (T (1/2) approximately 72 degrees C), and the stability was maintained in the pH range between 5 and 7. The ECP calorimetric melting curve showed, in addition to the main transition, a pretransitional conformational change with a T (1/2) of 44 degrees C. Both calorimetric transitions disappeared after successive re-heatings, and the ratio DeltaH versus DeltaH (vH) of 2.2 indicated a significant deviation from the two-state model. It was observed that the thermal unfolding was irreversible. The unfolding process gives rise to changes in the environment of aromatic amino acids that are partially maintained in the refolded protein with the loss of secondary structure and the formation of oligomers. From the thermodynamic analysis of ECP variants, the contribution of specific amino acids, such as Trp10 and the region 115-122, to thermal stability was also determined. The high thermal stability of ECP may contribute to its resistance to degradation when the protein is secreted to the extracellular medium during the immune response.     

Stability diagram and unfolding of a modified cytochrome c: what happens in the transformation regime?

We determined the stability diagram of a modified cytochrome c protein in a glycerol water mixture by measuring the first and the second moment of the fluorescence from the chromophore as a function of temperature and pressure. Temperature and pressure were varied between 273 and 363 K and 0.0001 and 1 GPa, respectively. The shift of the fluorescence maximum showed a characteristic sigmoid-like pattern from which information on the microscopic processes during unfolding is obtained: as the transformation regime is entered, the fluorescence shows a significant blue shift. The conclusion is that water molecules get into contact with the chromophore. They lead to strong electrostatic contributions in the solvent shift, which counteract the red shifting dispersion interactions. Assuming that there are just two relevant states that determine the stability diagram, the complete set of thermodynamic parameters can be determined from the data. However, under certain pressure-temperature conditions the fluorescence pattern is more complicated, pointing toward reentrant transitions and, possibly, to consecutive steps in the unfolding process.