Two-dimensional differential scanning calorimetry: simultaneous resolution of intrinsic protein structural energetics and ligand binding interactions by global linkage analysis

Extrapolating differential scanning calorimetry data for monoclonal antibodies to low temperatures

For irreversible denaturation transitions such as those exhibited by monoclonal antibodies, differential scanning calorimetry provides the denaturation temperature, Tm, the rate of denaturation at Tm, and the activation energy at Tm. These three quantities are essential but not sufficient for an accurate extrapolation of the rate of denaturation to temperatures of 25 °C and below. We have observed that the activation energy is not constant but temperature dependent due to the existence of an activation heat capacity, Cp,a. It is shown in this paper that a model that incorporates Cp,a is able to account for previous observations like, for example, that increasing the Tm does not always improve the stability at low temperatures; that some antibodies exhibit lower stabilities at 5 °C than at 25 °C; or that low temperature stabilities do not follow the rank order derived from Tm values. Most importantly, the activation heat capacity model is able to reproduce time dependent stabilities measured by size exclusion chromatography at low temperatures.

The effect of denaturants on protein thermal stability analyzed through a theoretical model considering multiple binding sites

A novel mathematical development applied to protein ligand binding thermodynamics is proposed, which allows the simulation, and therefore the analysis of the effects of multiple and independent binding sites to the Native and/or Unfolded protein conformations, with different binding constant values. Protein stability is affected when it binds to a small number of high affinity ligands or to a high number of low affinity ligands. Differential scanning calorimetry (DSC) measures released or absorbed energy of thermally induced structural transitions of biomolecules. This paper presents the general theoretical development for the analysis of thermograms of proteins obtained for n-ligands bound to the native protein and m-ligands bound to their unfolded form. In particular, the effect of ligands with low affinity and with a high number of binding sites (n and/or m > 50) is analyzed. If the interaction with the native form of the protein is the one that predominates, they are considered stabilizers and if the binding with the unfolded species predominates, it is expected a destabilizing effect. The formalism presented here can be adapted to fitting routines in order to simultaneously obtain the unfolding energy and ligand binding energy of the protein. The effect of guanidinium chloride on bovine serum albumin thermal stability, was successfully analyzed with the model considering low number of middle affinity binding sites to the native state and a high number of weak binding sites to the unfolded state.

Intrinsic affinity of protein - ligand binding by differential scanning calorimetry

Differential scanning calorimetry (DSC) determines the enthalpy change upon protein unfolding and the melting temperature of the protein. Performing DSC of a protein in the presence of increasing concentrations of specifically-binding ligand yields a series of curves that can be fit to obtain the protein-ligand dissociation constant as done in the fluorescence-based thermal shift assay (FTSA, ThermoFluor, DSF). The enthalpy of unfolding, as directly determined by DSC, helps improving the precision of the fit. If the ligand binding is linked to protonation reactions, the intrinsic binding constant can be determined by performing the affinity determination at a series of pH values. Here, the intrinsic, pH-independent, affinity of acetazolamide binding to carbonic anhydrase (CA) II was determined. A series of high-affinity ligands binding to CAIX, an anticancer drug target, and CAII showed recognition and selectivity for the anticancer isozyme. Performing the DSC experiment in buffers of highly different enthalpies of protonation enabled to observe the ligand unbinding-linked protonation reactions and estimate the intrinsic enthalpy of binding. The heat capacity of combined unfolding and unbinding was determined by varying the ligand concentrations. Taken together, these parameters provided a detailed thermodynamic picture of the linked ligand binding and protein unfolding process.

Thermodynamic and biophysical study of fatty acid effector binding to soybean lipoxygenase: implications for allostery driven by helix α2 dynamics

Previous comparative kinetic isotope effects have inferred an allosteric site for fatty acids and their derivatives that modulates substrate selectivity in 15-lipoxygenases. Hydrogen-deuterium exchange also previously revealed regionally defined enhanced protein flexibility, centred at helix α2 - a gate to the substrate entrance. Direct evidence for allosteric binding and a complete understanding of its mechanism remains elusive. In this study, we examine the binding thermodynamics of the fatty acid mimic, oleyl sulfate (OS), with the monomeric model plant 15-LOX, soybean lipoxygenase (SLO), using isothermal titration calorimetry. Dynamic light scattering and differential scanning calorimetry rule out OS-induced oligomerization or structural changes. These data provide evidence that the fatty acid allosteric regulation of SLO is controlled by the dynamics of helix α2.

Reversibility and irreversibility in the temperature denaturation of monoclonal antibodies

There have been numerous studies of the temperature denaturation of monoclonal antibodies (mAbs) using differential scanning calorimetry (DSC). In general, mAbs are characterized by complex temperature denaturation transitions in which the various domains (CH2, CH3, Fab) give rise to different peaks in the heat capacity function. The complexity and overall irreversibility of the temperature denaturation transition is well known and has limited the number of publications with an in-depth analysis of the data. Here we report that the temperature denaturation of the CH2 domain is reversible and only becomes irreversible after denaturation of the Fab domain, which is intrinsically irreversible. For these studies we have used the HIV neutralizing monoclonal antibody 17b. To account for the experimental heat capacity function, a mixed denaturation model that combines multiple reversible and irreversible transitions has been developed. This model accounts well for the DSC data and for the pH dependence of the heat capacity function of 17b and other monoclonal antibodies for which data is available in the literature. It is expected that a more detailed analysis of the stability of monoclonal antibodies will contribute to the development of better approaches to understand and optimize the structural viability of these therapeutic macromolecules.

Applied Biophysical Methods in Fragment-Based Drug Discovery

Fragment-based drug discovery (FBDD) has come of age in the last decade with the FDA approval of four fragment-derived drugs. Biophysical methods are at the heart of hit discovery and validation in FBDD campaigns. The three most commonly used methods, thermal shift, surface plasmon resonance, and nuclear magnetic resonance, can be daunting for the novice user. We aim here to provide the nonexpert user of these methods with a summary of problems and challenges that might be faced, but also highlight the potential gains that each method can contribute to an FBDD project. While our view on FBDD is slightly biased toward enabling structure-guided drug discovery, most of the points we address in this review are also valid for non-structure-focused FBDD.

Evaluation of the binding properties of drugs to albumin from DSC thermograms

There is a demand in rapid and robust methods to determine the affinity of drugs to receptors, enzymes, and transport proteins. Differential scanning calorimetry (DSC) is a common method to prove the existence of ligand-protein binding from the shift of denaturation peak, but it is rarely used to obtain the binding constant values. The work is aimed to prove that the DSC experiments can be a source of reliable values of the binding constants and information on the stoichiometry of drug-albumin binding. DSC thermograms of bovine serum albumin denaturation in the presence of several drugs with different affinity and stoichiometry of binding to albumin: naproxen, warfarin, ibuprofen, and isoniazid were recorded. The dependences of the denaturation peak maximum temperature and area on the molar drug/protein ratio, which varied from 0 to 100, were considered. With the help of numerical modeling of the DSC curves, these dependences were predicted using the binding parameters determined in independent experiments and a simple two-state model of denaturation. The DSC data at relatively small concentrations of ligands are in good agreement with the calculation results. The deviations from the model predictions at high ligand concentrations in the cases of naproxen and ibuprofen indicate that albumin is able to bind several additional molecules of these drugs with its low-affinity sites. The fit was improved by using a sequential binding model with two binding constants K₁ = 1.0 × 10⁷ and K₂ = 1.0 × 10⁴ for naproxen and a cooperative binding model for ibuprofen. The stoichiometry of drug-albumin complexes fully saturated with drug ligand was calculated from the dependence of the denaturation temperature on the drug concentration. In the case of isoniazid, DSC thermograms indicated very weak binding to albumin.

Homogeneous Assay for Target Engagement Utilizing Bioluminescent Thermal Shift

Protein thermal shift assays (TSAs) provide a means for characterizing target engagement through ligand-induced thermal stabilization. Although these assays are widely utilized for screening libraries and validating hits in drug discovery programs, they can impose encumbering operational requirements, such as the availability of purified proteins or selective antibodies. Appending the target protein with a small luciferase (NanoLuc) allows coupling of thermal denaturation with luminescent output, providing a rapid and sensitive means for assessing target engagement in compositionally complex environments such as permeabilized cells. The intrinsic thermal stability of NanoLuc is greater than mammalian proteins, and our results indicate that the appended luciferase does not alter thermal denaturation of the target protein. We have successfully applied the NanoLuc luciferase thermal shift assay (NaLTSA) to several clinically relevant protein families, including kinases, bromodomains, and histone deacetylases. We have also demonstrated the suitability of this assay method for library screening and compound profiling.

Direct and indirect approaches to identify drug modes of action

Phenotypic assays are becoming increasingly more common among drug discovery practices, expanding drug target diversity as lead compounds identified through such screens are not limited to known targets. While increasing diversity is beneficial to the drug discovery process and the fight against disease, the unknown modes of action of new lead compounds can hamper drug discovery as, in most cases, the process of lead compound optimization is made difficult due to the unknown nature of the target; blindly changing substituents can prove fruitless due to the inexhaustible number of potential combinations, and it is therefore desirable to rapidly identify the targets of lead compounds developed through phenotypic screening. In addition, leads identified through target-based screening often have off-target effects that contribute towards drug toxicity, and by identifying those secondary targets, the drugs can be improved. However, the identification of a leads mode of action is far from trivial and now represents a major bottleneck in the drug discovery pipeline. This review looks at some of the recent developments in the identification of drug modes of action, focusing on phenotype-based methods using metabolomics, proteomics, transcriptomics, and genomics to detect changes in phenotype in response to the presence of the drug, and affinity-based methods using modified/unmodified drug as bait to capture and identify targets. © 2017 IUBMB Life, 70(1):9-22, 2018.

A unified framework based on the binding polynomial for characterizing biological systems by isothermal titration calorimetry

Isothermal titration calorimetry (ITC) has become the gold-standard technique for studying binding processes due to its high precision and sensitivity, as well as its capability for the simultaneous determination of the association equilibrium constant, the binding enthalpy and the binding stoichiometry. The current widespread use of ITC for biological systems has been facilitated by technical advances and the availability of commercial calorimeters. However, the complexity of data analysis for non-standard models is one of the most significant drawbacks in ITC. Many models for studying macromolecular interactions can be found in the literature, but it looks like each biological system requires specific modeling and data analysis approaches. The aim of this article is to solve this lack of unity and provide a unified methodological framework for studying binding interactions by ITC that can be applied to any experimental system. The apparent complexity of this methodology, based on the binding polynomial, is overcome by its easy generalization to complex systems.

Affinity assays for decrypting protein targets of unknown function

The application of label-independent biophysical assays (affinity assays) to carefully selected libraries of biochemicals, natural products and known bioactive molecules generates a powerful tool to identify ligands indicative of protein function. Detection of ligand binding does not require the measurement of enzyme activity, thus every cofactor, ion, reactant or product might be identified separately. One such assay, ThermoFluor(®), monitors ligand perturbation of protein thermostability. In conjunction with a functional probe library (FPL), this assay has been used to decrypt the function of orphan proteins.:

Chip-based nanoelectrospray mass spectrometry for protein characterization

In the last several years, significant progress has been made in the development of microfluidic-based analytical technologies for proteomic and drug discovery applications. Chip-based nanoelectrospray coupled to a mass spectrometer detector is one of the recently developed analytical microscale technologies. This technology offers unique advantages for automated nanoelectrospray including reduced sample consumption, improved detection sensitivity and enhanced data quality for proteomic studies. This review presents an overview and introduction of recent developments in chip devices coupled to electrospray mass spectrometers including the development of the automated nanoelectrospray ionization chip device for protein characterization. Applications using automated chip-based nanoelectrospray ionization technology in proteomic and bioanalytical studies are also extensively reviewed in the fields of high-throughput protein identification, protein post-translational modification studies, top-down proteomics, biomarker screening by pattern recognition, noncovalent protein-ligand binding for drug discovery and lipid analysis. Additionally, future trends in chip-based nanoelectrospray technology are discussed.

Target identification strategies in plant chemical biology

The current needs to understand gene function in plant biology increasingly require more dynamic and conditional approaches opposed to classic genetic strategies. Gene redundancy and lethality can substantially complicate research, which might be solved by applying a chemical genetics approach. Now understood as the study of small molecules and their effect on biological systems with subsequent target identification, chemical genetics is a fast developing field with a strong history in pharmaceutical research and drug discovery. In plant biology however, chemical genetics is still largely in the starting blocks, with most studies relying on forward genetics and phenotypic analysis for target identification, whereas studies including direct target identification are limited. Here, we provide an overview of recent advances in chemical genetics in plant biology with a focus on target identification. Furthermore, we discuss different strategies for direct target identification and the possibilities and challenges for plant biology.

Two-state dynamics of the SH3-SH2 tandem of Abl kinase and the allosteric role of the N-cap

The regulation and localization of signaling enzymes is often mediated by accessory modular domains, which frequently function in tandems. The ability of these tandems to adopt multiple conformations is as important for proper regulation as the individual domain specificity. A paradigmatic example is Abl, a ubiquitous tyrosine kinase of significant pharmacological interest. SH3 and SH2 domains inhibit Abl by assembling onto the catalytic domain, allosterically clamping it in an inactive state. We investigate the dynamics of this SH3-SH2 tandem, using microsecond all-atom simulations and differential scanning calorimetry. Our results indicate that the Abl tandem is a two-state switch, alternating between the conformation observed in the structure of the autoinhibited enzyme and another configuration that is consistent with existing scattering data for an activated form. Intriguingly, we find that the latter is the most probable when the tandem is disengaged from the catalytic domain. Nevertheless, an amino acid stretch preceding the SH3 domain, the so-called N-cap, reshapes the free-energy landscape of the tandem and favors the interaction of this domain with the SH2-kinase linker, an intermediate step necessary for assembly of the autoinhibited complex. This allosteric effect arises from interactions between N-cap and the SH2 domain and SH3-SH2 connector, which involve a phosphorylation site. We also show that the SH3-SH2 connector plays a determinant role in the assembly equilibrium of Abl, because mutations thereof hinder the engagement of the SH2-kinase linker. These results provide a thermodynamic rationale for the involvement of N-cap and SH3-SH2 connector in Abl regulation and expand our understanding of the principles of modular domain organization.

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 (Tm shift) has become a popular approach to identify potential ligands. However, Tm 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.

Improved accuracy of low affinity protein-ligand equilibrium dissociation constants directly determined by electrospray ionization mass spectrometry

There is continued interest in the determination by ESI-MS of equilibrium dissociation constants (K(D)) that accurately reflect the affinity of a protein-ligand complex in solution. Issues in the measurement of K(D) are compounded in the case of low affinity complexes. Here we present a K(D) measurement method and corresponding mathematical model dealing with both gas-phase dissociation (GPD) and aggregation. To this end, a rational mathematical correction of GPD (f(sat)) is combined with the development of an experimental protocol to deal with gas-phase aggregation. A guide to apply the method to noncovalent protein-ligand systems according to their kinetic behavior is provided. The approach is validated by comparing the K(D) values determined by this method with in-solution K(D) literature values. The influence of the type of molecular interactions and instrumental setup on f(sat) is examined as a first step towards a fine dissection of factors affecting GPD. The method can be reliably applied to a wide array of low affinity systems without the need for a reference ligand or protein.

Unbiased binding assays for discovering small-molecule probes and drugs

2011 marks the 10-year anniversary of milestone manuscripts describing drafts of the human genome sequence. Over the past decade, a number of new proteins have been linked to disease-many of which fall into classes that have been historically considered challenging from the perspective of drug discovery. Several of these newly associated proteins lack structural information or strong annotation with regard to function, making development of conventional in vitro functional assays difficult. A recent resurgence in the popularity of simple small molecule binding assays has led to new approaches that do not require knowledge of protein structure or function in advance. Here we briefly review selected methods for executing binding assays that have been used successfully to discover small-molecule probes or drug candidates.

Conformational stability of hepatitis C virus NS3 protease

The hepatitis C virus NS3 protease is responsible for the processing of the nonstructural region of viral precursor polyprotein in infected hepatic cells. NS3 has been considered a target for drug discovery for a long time. NS3 is a zinc-dependent serine protease. However, the zinc ion is not involved in the catalytic mechanism, because it is bound far away from the active site. Thus, zinc is essential for the structural integrity of the protein and it is considered to have a structural role. The first thermodynamic study on the conformational equilibrium and stability of NS3 and the effect of zinc on such equilibrium is presented here. In agreement with a previous calorimetric study on the binding of zinc to NS3, the global unfolding heat capacity is dominated by the zinc dissociation step, suggesting that the binding of zinc induces a significant structural rearrangement of the protein. In addition, contrary to other homologous zinc-dependent proteases, the zinc-free NS3 protease is not completely unstructured. It is apparent that the conformational landscape of hepatitis C virus NS3 protease is fairly complex due to its intrinsic plasticity, and to the interactions with its different effectors (zinc and the accessory viral protein NS4A) and their modulation of the population of the different conformational states.

The N-terminal domain of the enzyme I is a monomeric well-folded protein with a low conformational stability and residual structure in the unfolded state

The bacterial phosphoenolpyruvate-dependent sugar phosphotransferase system is a multiprotein complex that phosphorylates and, concomitantly, transports carbohydrates across the membrane into the cell. The first protein of the cascade is a multidomain protein so-called enzyme I (EI). The N-terminal domain of EI from Streptomyces coelicolor, EIN(sc), responsible for the binding to the second protein in the cascade (the histidine phosphocarrier, HPr), was cloned and successfully expressed and purified. We have previously shown that EI(sc) binds to HPr(sc) with smaller affinity than other members of the EI and HPr families [Hurtado-Gómez et al. (2008) Biophys. J., 95, 1336-1348]. We think that the study of the isolated binding HPr(sc) domain, that is EIN(sc), could shed light on the small affinity value measured. Therefore, in this work we present a detailed description of the structural features of the EIN domain, as a first step towards a complete characterization of the molecular recognition process between the two proteins. We show that EIN(sc) is a folded protein, with alpha-helix and beta-sheet structures and also random-coil conformations, as shown by circular dichroism (CD), FTIR and NMR spectroscopies. The acquisition of secondary and tertiary structures, and the burial of hydrophobic regions, occurred concomitantly at acidic pHs, but at very low pH, the domain acquired a molten-globule conformation. The EIN(sc) protein was not very stable, with an apparent conformational free energy change upon unfolding, DeltaG, of 4.1 +/- 0.4 kcal mol(-1), which was pH independent in the range explored (from pH 6.0 to 8.5). The thermal denaturation midpoint, which was also pH invariant, was similar to that measured in the isolated intact EI(sc). Although EIN(sc) shows thermal- and chemical denaturations that seems to follow a two-state mechanism, there is evidence of residual structure in the chemical and thermally unfolded states, as indicated by differential scanning calorimetry and CD measurements.

Relationship between the wavelength maximum of a protein and the temperature dependence of its intrinsic tryptophan fluorescence intensity

Proper determination of the temperature dependence of intrinsic tryptophan fluorescence intensity in native and denatured states is an essential prerequisite for extracting the free energy of protein unfolding from the thermal denaturation profile. The most common method employed in determining the temperature dependence of these conformations is through the determination of slopes of pre- and post-transition baselines. However, simulations of protein unfolding profiles suggest that this method does not work for marginally stable proteins. We show herein that the temperature dependence of the fluorescence intensity of N-acetyl tryptophanamide (NATA) in organic solvents and water may be used to represent the temperature dependence of the fluorescence intensity of tryptophan in native and denatured conformations of a protein, respectively. The wavelength of the emission maximum, lambda (max), of N-acetyl tryptophanamide (NATA) in a particular solvent or tryptophan in proteins is related to the temperature dependence (m) of its fluorescence intensity by the equation: m (K(-1)) = (-0.000299 +/- 2.2 x 10(-5) K(-1) nm(-1)) x lambda (max) (nm) + (0.0919 +/- 0.0025 K(-1)).

Discovery of specific flavodoxin inhibitors as potential therapeutic agents against Helicobacter pylori infection

Helicobacter pylori establishes life-long infections in the gastric mucosa of over 1 billion people worldwide. In many cases, without specific antimicrobial intervention, H. pylori infected individuals will develop type B gastritis, chronic peptic ulcers and, more rarely, gastric neoplasias. Conventional antimicrobial therapy has been complicated by dramatic increases in resistance to macrolides, metronidazole and fluoroquinolones. Here, we report the development of novel therapeutics that specifically target the unique flavodoxin component of an essential metabolic pathway of H. pylori. With the use of high-throughput screening methodology, we have tested 10,000 chemicals and have identified 29 compounds that bind flavodoxin, four of which interrupted in vitro electron transfer to flavodoxin physiological partners. Three of these compounds are bactericidal and promisingly selective for H. pylori. The minimal inhibitory concentrations of two of them are 10 times lower than their minimal cytotoxic concentrations for HeLa cells. Importantly, neither of the four inhibitors is toxic for mice after administration of 1-10 mg kg(-1) doses twice a day for 5 days. Enzymatic, thermodynamic and structural characterization of the inhibitor-flavodoxin complexes suggests these compounds could act by modifying the redox potentials of flavodoxin. These newly discovered inhibitors represent promising selective leads against the different diseases associated to H. pylori infection.

Comparison of Two ESI MS Based H/D Exchange Methods for Extracting Protein Folding Energies

In this report, the model proteins staphylococcal nuclease and ubiquitin were used to test the applicability of two new hydrogen/deuterium exchange (HX) electrospray ionization mass spectrometry (ESI-MS) methods for estimating protein folding energies. Both methods use the H/D exchange of globally protected amide protons (amide protons which are buried in the hydrophobic core) to elucidate protein folding energies. One method is a kinetic-based method and the other is equilibrium-based. The first method, the HX ESI-MS kinetic-based approach is conceptually identical to SUPREX (stability of unpurified proteins from rates of H/D exchange) method but is based on ESI-MS rather than MALDI-MS (matrix assisted laser desorption mass spectrometry). This method employs the time-dependence of H/D exchange using various denaturant concentrations to extract folding energies. Like SUPREX, this approach requires the assumption of EX2 exchange kinetics. The second method, which we call a protein equilibrium population snapshot (PEPS) by HX ESI-MS uses data collected only for a single time point (usually the shortest possible) to obtain a snapshot of the open and closed populations of the protein. The PEPS approach requires few assumptions in the derivation of the equations used for calculation of the folding energies. The extraction of folding energies from mass spectral data is simple and straightforward. The PEPS method is applicable for proteins that follow either EX1 or EX2 HX mechanisms. In our experiments the kinetic-based method produced less accurate ΔG(H(2)O) and m(GdHCl) values for wild-type staphylococcal nuclease and mutants undergoing H/D exchange by EX1, as would be expected. Better results were obtained for ubiquitin which undergoes HX by an EX2 mechanism. Using the PEPS method we obtained ΔG(H(2)O) and m(GdHCl) values that were in good agreement with literature values for both staphylococcal nuclease (EX1) and ubiquitin (EX2). We also show that the observation of straight lines in linear extrapolation method (LEM) plots is not a reliable indicator of the validity of the data obtained using the LEM approach.

Computing protein stabilities from their chain lengths

New amino acid sequences of proteins are being learned at a rapid rate, thanks to modern genomics. The native structures and functions of those proteins can often be inferred using bioinformatics methods. We show here that it is also possible to infer the stabilities and thermal folding properties of proteins, given only simple genomics information: the chain length and the numbers of charged side chains. In particular, our model predicts DeltaH(T), DeltaS(T), DeltaC(p), and DeltaF(T)--the folding enthalpy, entropy, heat capacity, and free energy--as functions of temperature T; the denaturant m values in guanidine and urea; the pH-temperature-salt phase diagrams, and the energy of confinement F(s) of the protein inside a cavity of radius s. All combinations of these phase equilibria can also then be computed from that information. As one illustration, we compute the pH and salt conditions that would denature a protein inside a small confined cavity. Because the model is analytical, it is computationally efficient enough that it could be used to automatically annotate whole proteomes with protein stability information.

Mechanism of low density lipoprotein (LDL) release in the endosome: implications of the stability and Ca2+ affinity of the fifth binding module of the LDL receptor

Uptake of low density lipoproteins (LDL) by their receptor, LDLR, is the primary mechanism by which cells incorporate cholesterol from plasma. Mutations in LDLR lead to familial hypercholesterolemia, a common disease affecting 1 in 500 of the human population. LDLR is a modular protein that uses several small repeats to bind LDL. The repeats contain around 40 residues, including three disulfide bonds and a calcium ion. Repeat 5 (LR5) is critical for LDL and beta-migrating very low density lipoprotein binding. Based on the crystal structure of LDLR at endosomal pH (but close to extracellular calcium concentration), LR5 has been proposed to bind to the epidermal growth factor (EGF) precursor domain of LDLR in the endosome, thus releasing the LDL particles previously bound in extracellular conditions. We report here the conformational stability of LR5 as a function of temperature and calcium concentration under both extracellular and endosomal pH conditions. The repeat was very stable when it bore a bound calcium ion but was severely destabilized in the absence of calcium and even further destabilized at acidic versus neutral pH. The temperature and calcium concentration dependence of LR5 stability clearly indicate that under endosomal conditions the unfolded conformation of the repeat is largely dominant. We thus propose a new mechanism for LDL release in the endosome in which calcium depletion and decreased stability at acidic pH drives LR5 unfolding, which triggers LDL release from the receptor. Subsequent binding of LR5 to the EGF precursor domain, if it takes place at low calcium concentrations, would contribute to a further shifting of the equilibrium toward dissociation.

Molten globule and native state ensemble of Helicobacter pylori flavodoxin: can crowding, osmolytes or cofactors stabilize the native conformation relative to the molten globule?

Partly unfolded protein conformations close in energy to the native state may be involved in protein functioning and also be related to folding diseases, but yet their structure and energetics are poorly understood. One such conformation, the monomeric and well-behaved molten globule of Helicobacter pylori apoflavodoxin, is here investigated to provide, in a wide pH interval, a complete thermodynamic description of its unfolding equilibrium and the equilibrium linking molten globule and native state. All thermodynamic and molecular properties of the molten globule here analyzed are characteristic of a partly unfolded conformation, and their differences with those of the native state are typically quantitative rather than qualitative. The stability data depict a native state ensemble where the relative populations of the different intermediates are strongly modulated by pH. Whereas the molten globule is dominant at pH 2.0, at neutral pH it is just the least stable of three partly unfolded intermediates populated by this protein. It is of interest that the energy rank of these intermediates at pH 7.0 is consistent with their likelihood to overcome the native state and become the more stable conformation when the native state protein is subjected to heat or mutation stress. Given the small volume difference between molten globule and native state, neither crowding agents nor osmolytes can drive the molten globule back to the native state. This observation, which is in qualitative accord with predictions of simple excluded volume theory, indicates that molecular crowding in vivo is not an effective mechanism to minimize partial unfolding events leading to equilibrium intermediates.

The flavodoxin from Helicobacter pylori: structural determinants of thermostability and FMN cofactor binding

Flavodoxin has been recently recognized as an essential protein for a number of pathogenic bacteria including Helicobacter pylori, where it has been proposed to constitute a target for antibacterial drug development. One way we are exploring to screen for novel inhibitory compounds is to perform thermal upshift assays, for which a detailed knowledge of protein thermostability and cofactor binding properties is of great help. However, very little is known on the stability and ligand binding properties of H. pylori flavodoxin, and its peculiar FMN binding site together with the variety of behaviors observed within the flavodoxin family preclude extrapolations. We have thus performed a detailed experimental and computational analysis of the thermostability and cofactor binding energetics of H. pylori flavodoxin, and we have found that the thermal unfolding equilibrium is more complex that any other previously described for flavodoxins as it involves the accumulation of two distinct equilibrium intermediates. Fortunately the entire stability and binding data can be satisfactorily fitted to a model, summarized in a simple phase diagram, where the cofactor only binds to the native state. On the other hand, we show how variability of thermal unfolding behavior within the flavodoxin family can be predicted using structure-energetics relationships implemented in the COREX algorithm. The different distribution and ranges of local stabilities of the Anabaena and H. pylori apoflavodoxins explain the essential experimental differences observed: much lower Tm1, greater resistance to global unfolding, and more pronounced cold denaturation in H. pylori. Finally, a new strategy is proposed to identify using COREX structural characteristics of equilibrium intermediate states populated during protein unfolding.

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.

Differential scanning calorimetry as a tool to estimate binding parameters in multiligand binding proteins

The stability of proteins and their interactions with other molecules is a topic of special interest in biochemistry because many cellular processes depend on that. New methods and approaches are constantly developed to elucidate the energetics of biomolecular recognition. In this sense, the application of the theory of macromolecular unfolding linked to ligand binding to differential scanning calorimetry (DSC) has proved to be a useful tool to simultaneously characterize the energetics of unfolding and binding. Although the general theory is well known, the applicability of DSC to study the interaction of biomolecules is not common. In the current work, we estimated the binding parameters of 8-anilinonaphthalene-1-sulfonic acid to human serum albumin using DSC. This model system was chosen due to both the complex stoichiometry and the moderate binding constants. From DSC curves acquired at different ligand concentrations, we obtained the number of bound ligands, the binding constants, and the binding enthalpy for each independent binding site. Compared with those parameters determined by titration calorimetry, the results highlight the potentiality of DSC to estimate binding parameters in multiligand binding proteins.

Molecular basis for increased susceptibility of isolates with atazanavir resistance-conferring substitution I50L to other protease inhibitors

Protease inhibitors (PIs) are highly effective drugs against the human immunodeficiency virus (HIV), yet long-term therapeutic use is limited by emergence of HIV type 1 (HIV-1) protease substitutions that confer cross-resistance to multiple protease inhibitor drugs. Atazanavir is a highly potent HIV protease inhibitor with a distinct resistance profile that includes effectiveness against most HIV-1 isolates resistant to one or two PIs. The signature resistance substitution for atazanavir is I50L, and it is frequently (53%) accompanied by a compensatory A71V substitution that helps restore viability and increases atazanavir resistance levels. We measured the binding affinities of wild-type (WT) and I50L/A71V HIV-1 proteases to atazanavir and other currently approved PIs (ritonavir, lopinavir, saquinavir, nelfinavir, indinavir, and amprenavir) by isothermal titration calorimetry. Remarkably, we find that all of the PIs have 2- to 10-fold increased affinities for I50L/A71V protease, except for atazanavir. The results are also manifested by thermal stability measures of affinity for WT and I50L/A71V proteases. Additional biophysical and enzyme kinetics experiments show I50L/A71V protease is a stable enzyme with catalytic activity that is slightly reduced (34%) relative to the WT. Computational modeling reveals that the unique resistance phenotype of I50L/A71V protease likely originates from bulky tert-butyl groups at P2 and P2' (specific to atazanavir) that sterically clash with methyl groups on residue L50. The results of this study provide a molecular understanding of the novel hypersusceptibility of atazanavir-resistant I50L/A71V-containing clinical isolates to other currently approved PIs.

Ligand-induced thermostability in proteins: thermodynamic analysis of ANS-albumin interaction

A comparative thermodynamic study of the interaction of anilinonaphthalene sulfonate (ANS) derivatives with bovine serum albumin (BSA) was performed by using differential scanning calorimetry (DSC) and isothermal titration calorimetry (ITC). The chemically related ligands, 1,8-ANS and 2,6-ANS, present a similar affinity for BSA with different binding energetics. The analysis of the binding driving forces suggests that not only hydrophobic effect but also electrostatic interactions are relevant, even though they have been extensively used as probes for non-polar domains in proteins. Ligand association leads to an increase in protein thermostability, indicating that both dyes interact mainly with native BSA. ITC data show that 1,8-ANS and 2,6-ANS have a moderate affinity for BSA, with an association constant of around 1-9x10(5) M(-1) for the high-affinity site. Ligand binding is disfavoured by conformational entropy. The theoretical model used to simulate DSC data satisfactorily reproduces experimental thermograms, validating this approach as one which provides new insights into the interaction between one or more ligands with a protein. By comparison with 1,8-ANS, 2,6-ANS appears as a more "inert" probe to assess processes which involve conformational changes in proteins.

Isothermal titration calorimetry

Isothermal titration calorimetry is an ideal technique for measuring biological binding interactions. It does not rely on the presence of chromophores or fluorophores, nor does it require an enzymatic assay. Because the technique relies only on the detection of a heat effect upon binding, it can be used to measure the binding constant, K, the enthalpy of binding, DeltaH degrees and the stoichiometry, or number of binding sites, n. This chapter describes instrumentation, experimental design, and the theoretical underpinnings necessary to run and analyze a calorimetric binding experiment.

Ligand-linked stability of mutants of the C-domain of calmodulin

There is a necessary energetic linkage between ligand binding and stability in biological molecules. The critical glutamate in Site 4 was mutated to create two mutants of the C-domain of calmodulin yielding E140D and E140Q. These proteins were stably folded in the absence of calcium, but had dramatically impaired binding of calcium. We determined the stability of the mutant proteins in the absence and presence of calcium using urea-induced unfolding monitored by circular dichroism (CD) spectroscopy. These calcium-dependent unfolding curves were fit to models that allowed for linkage of stability to binding of a single calcium ion to the native and unfolded states. Simultaneous analysis of the unfolding profiles for each mutant yielded estimates for calcium-binding constants that were consistent with results from direct titrations monitored by fluorescence. Binding to the unfolded state was not an important energetic contributor to the ligand-linked stability of these mutants.

Thermal Denaturation: A Method to Rank Slow Binding, High-Affinity P38α MAP Kinase Inhibitors

It has been reported that the diaryl urea class of p38alpha inhibitors binds to p38 map kinase with both high affinity and slow binding kinetics (Pargellis et al. Nat. Struct. Biol. 2002, 9, 268-272). The slow binding kinetics of this class of inhibitors is believed to be the result of binding to an allosteric pocket adjacent to the p38alpha active site. The use of traditional kinetic and equilibrium methods to measure the binding affinity of this class of compounds has created many challenges for determination of structure-activity relationships (SAR). The thermal denaturation method provides a means of measuring high-affinity interactions. In this paper, the method of thermal denaturation will be described as it has been applied to the diaryl urea class of p38 map kinase inhibitors.

Quantitative determination of noncovalent binding interactions using automated nanoelectrospray mass spectrometry

Electrospray ionization mass spectrometry (ESI-MS) has proven to be an extremely powerful tool for studying biomolecular structures and noncovalent interactions. Here we report a method using a fully automated, chip-based nanoESI-MS system to determine the dissociation constants (Kd) for the complexes of two different proteins with their ligands. The automated nanoelectrospray system, consisting of the NanoMate and ESI chip, serves functionally as a combination of autosampler and nanoelectrospray ionization source. This system provides all the advantages of conventional nanoelectrospray plus automated, high-throughput analyses without carryover. The automated nanoESI system was used to investigate quantitative noncovalent interactions between ribonuclease A (RNase A) and cytidylic acid ligands (2'-CMP, CTP), a well-characterized model protein-ligand complex, and between an inactive endocellulase mutant (Thermobifida fusca Cel6A D117Acd) and four oligosaccharide ligands (cellotriose, cellotetraose, cellopentaose, cellohexaose). Both titration and competitive binding approaches were performed prior to automated nanoESI-MS analysis with a Q-TOF mass spectrometer. Dissociation constants for each complex were calculated from the sum of ion peak areas of free and complexed proteins during the titration and competition experiments. The measured Kd values for the RNase A-CMP and Cel6A D117Acd-G3 complexes were found to be in excellent agreement with the available published values obtained by standard spectroscopic titration techniques. To our knowledge, this is the first report of using an ESI-MS approach to study the interactions between a cellulase and oligosaccharides. The results provide new insights for understanding the nature of cellulase-cellulose interactions.

Protein stability induced by ligand binding correlates with changes in protein flexibility

The interaction between ligands and proteins usually induces changes in protein thermal stability with modifications in the midpoint denaturation temperature, enthalpy of unfolding, and heat capacity. These modifications are due to the coupling of unfolding with binding equilibrium. Furthermore, they can be attained by changes in protein structure and conformational flexibility induced by ligand interaction. To study these effects we have used bovine serum albumin (BSA) interacting with three different anilinonaphthalene sulfonate derivatives (ANS). These ligands have different effects on protein stability, conformation, and dynamics. Protein stability was studied by differential scanning calorimetry and fluorescence spectroscopy, whereas conformational changes were detected by circular dichroism and infrared spectroscopy including kinetics of hydrogen/deuterium exchange. The order of calorimetric midpoint of denaturation was: 1,8-ANS-BSA > 2,6-ANS-BSA > free BSA >> (nondetected) bis-ANS-BSA. Both 1,8-ANS and 2,6-ANS did not substantially modify the secondary structure of BSA, whereas bis-ANS induced a distorted alpha-helix conformation with an increase of disordered structure. Protein flexibility followed the order: 1,8-ANS-BSA < 2,6-ANS-BSA < free BSA << bis-ANS-BSA, indicating a clear correlation between stability and conformational flexibility. The structure induced by an excess of bis-ANS to BSA is compatible with a molten globule-like state. Within the context of the binding landscape model, we have distinguished five conformers (identified by subscript): BSA(1,8-ANS), BSA(2,6-ANS), BSA(free), BSA(bis-ANS), and BSA(unfolded) among the large number of possible states of the conformational dynamic ensemble. The relative population of each distinguishable conformer depends on the type and concentration of ligand and the temperature of the system.

Phase diagrams: a graphical representation of linkage relations

It is shown here that phase diagrams of ligand-binding biological macromolecules provide a powerful tool for the analysis of reaction mechanisms. The present study provides simple rules for the construction and interpretation of such phase diagrams. We give examples for the derivation of reaction schemes for macromolecules that can bind two different kinds of ligands. By sampling one dimension of a phase diagram it is possible to reconstruct the second dimension, including the correct stoichiometry, positive and negative linkage between the ligands and equilibrium binding constants for the complete series of reactions. The discussion is generalised to temperature and pressure-dependent phase diagrams. To exemplify the new diagram method we analyse the pH-dependent binding of trans-beta-indole acrylic acid to apo-Trp repressor, the pH-dependent thermal denaturation of alpha-chymotrypsinogen A, calcium binding and denaturation of annexin I, high affinity zinc binding to a metallo-beta-lactamase and high-pressure and temperature denaturation of RNase A and staphylococcal nuclease.

A general mass spectrometry-based assay for the quantitation of protein-ligand binding interactions in solution

A new method that utilizes matrix-assisted laser desorption/ionization (MALDI) mass spectrometry and exploits the hydrogen/deuterium (H/D) exchange properties of proteins was developed for measuring the thermodynamic properties of protein-ligand complexes in solution. Dissociation constants (Kd values) determined by the method for five model protein-ligand complexes that included those with small molecules, nucleic acids, peptides, and other proteins were generally in good agreement with Kd values measured by conventional methods. Important experimental advantages of the described method over existing methods include: the ability to make measurements in a high-throughput and automated fashion, the ability to make measurements using only picomole quantitities of protein, and the ability to analyze either purified or unpurified protein-ligand complexes.

The linkage between protein folding and functional cooperativity: two sides of the same coin?

During the course of their biological function, proteins undergo different types of structural rearrangements ranging from local to large-scale conformational changes. These changes are usually triggered by their interactions with small-molecular-weight ligands or other macromolecules. Because binding interactions occur at specific sites and involve only a small number of residues, a chain of cooperative interactions is necessary for the propagation of binding signals to distal locations within the protein structure. This process requires an uneven structural distribution of protein stability and cooperativity as revealed by NMR-detected hydrogen/deuterium exchange experiments under native conditions. The distribution of stabilizing interactions does not only provide the architectural foundation to the three-dimensional structure of a protein, but it also provides the required framework for functional cooperativity. In this review, the statistical thermodynamic linkage between protein stability, functional cooperativity, and ligand binding is discussed.

The heat capacity paradox of ligand binding proteins: reconciling the microscopic and macroscopic world

Differential scanning microcalorimetry (DSC) is a superb method for the analysis of protein energetics. However, the relative simplicity of application has led astray many to assume that a proper analysis of the data was possible without a sound knowledge of the underlying statistical thermodynamic principles. In this study, the question is addressed of how to calculate properly the heat capacity signal of a protein in the presence of high affinity ligands. It is shown that the signal corresponds neither to grand canonic nor to canonic heat capacity. Statistical thermodynamic model calculations result only in the observed macroscopic heat capacity signal, if the protein in the calorimetric cell is assumed to form a grand canonic ensemble (T, p, mu controlled) which is, however, heated under constraints typical for a canonic ensemble (T, p, N controlled). As a consequence, the microscopic statistical thermodynamic heat capacity must be carefully distinguished from the macroscopically observable thermodynamic heat capacity in those cases where proteins unfold in the presence of high affinity ligands.

Thermodynamic effects of active-site ligands on the reversible, partial unfolding of dodecameric glutamine synthetase from Escherichia coli: calorimetric studies

Dodecameric glutamine synthetase (GS) from Escherichia coli undergoes reversible, thermally induced partial unfolding without subunit dissociation. A single endotherm for Mn.GS (+/- active-site ligands) in the presence of 1 mM free Mn2+ and 100 mM KCl at pH 7 is observed by differential scanning calorimetry (DSC). Previous deconvolutions of DSC data for Mn.GS showed only two two-state transitions (with similar tm values; 51.6 +/- 2 degrees C), and indicated that cooperative interactions link partial unfolding reactions of all subunits within the Mn.enzyme dodecamer [Ginsburg, A., & Zolkiewski, M. (1991) Biochemistry 30, 9421]. A net uptake of 8.0 equiv of H+ by Mn.GS occurs during partial unfolding, as determined in the present DSC experiments conducted with four buffers having different heats of protonation at 50 degrees C. These data gave a value of 176 +/- 12 kcal (mol of dodecamer)-1 for delta Hcal corrected for buffer protonation. L-Glutamine and L-Met-(SR)-sulfoximine stabilize the Mn.GS dodecamer through the free energies of ligand binding, and these were shown to be partially and totally released, respectively, from the 12 active sites at high temperature. Ligand effects on Tm values from DSC were similar to those from spectral measurements of Trp and Tyr exposures in two subunit domains. Effects of varying [ADP] on DSC profiles of Mn.GS were complex; Tm is increased by low [ADP] and decreased by > 100 microM free ADP. This is due to the exposure of an additional low-affinity ADP binding site per GS subunit at high temperature with log K1' = 4.3 and log K2' = 3.6 at 60 degrees C relative to log K' = 5.5 for ADP binding at 30 degrees C, as determined by isothermal calorimetric and fluorescence titrations. Moreover, delta Hcal at > 27% saturation with ADP (corrected for ADP binding/dissociation) is approximately 80-100 kcal/mol more than in the absence of ligands. Changes in domain interactions could result from ADP bridging subunit contacts in the dodecamer. Each of the active-site ligands investigated here produces different effects on DSC profiles without uncoupling the extremely cooperative, partial unfolding reactions in the Mn.GS dodecamer.

Determination of RNase A/2'-cytidine monophosphate binding affinity and enthalpy by a global fit of thermal unfolding curves

A spectropolarimeter was used to measure the thermal response curves of RNase A in the presence and absence of the ligand cytidine 2'-monophosphate. A coupled equilibrium model was used to describe the dissociation of the protein-ligand complex (NL<==>N + L) and the thermal unfolding of the free protein (N<==>U). The unfolding curves of the protein in the presence of several different concentrations of ligand were fit to this coupled equilibrium model using global linkage analysis. The best-fitted values for the thermal unfolding of the apo-protein were 60.9 +/- 0.2 degrees C (T(m)) and 105.5 +/- 1.4 kcal/mol (DeltaH), while the fitted values for the dissociation of the protein-ligand complex were 1.6 +/- 0.4 microM (K(D)) and 18.7 +/- 1.0 kcal mol(-1) (DeltaH(L)). These values were in excellent agreement with values obtained by other methods. The sensitivity of the data fitting was compared using linear or quadratic temperature dependence for the response curves of the free ligand (L), native apo-protein (N), native ligand-bound protein (NL), and unfolded apo-protein (U). There was no significant improvement in the precision of the fitted data when the temperature-dependent response for each species (N, L, NL, and U) was expressed as quadratic functions of temperature.