Molecular basis of co-operativity in protein folding. III. Structural identification of cooperative folding units and folding intermediates

Protective effects of combining monoclonal antibodies and vaccines against the Plasmodium falciparum circumsporozoite protein

Combinations of monoclonal antibodies (mAbs) against different epitopes on the same antigen synergistically neutralize many viruses. However, there are limited studies assessing whether combining human mAbs against distinct regions of the Plasmodium falciparum (Pf) circumsporozoite protein (CSP) enhances in vivo protection against malaria compared to each mAb alone or whether passive transfer of PfCSP mAbs would improve protection following vaccination against PfCSP. Here, we isolated a panel of human mAbs against the subdominant C-terminal domain of PfCSP (C-CSP) from a volunteer immunized with radiation-attenuated Pf sporozoites. These C-CSP-specific mAbs had limited binding to sporozoites in vitro that was increased by combination with neutralizing human "repeat" mAbs against the NPDP/NVDP/NANP tetrapeptides in the central repeat region of PfCSP. Nevertheless, passive transfer of repeat- and C-CSP-specific mAb combinations did not provide enhanced protection against in vivo sporozoite challenge compared to repeat mAbs alone. Furthermore, combining potent repeat-specific mAbs (CIS43, L9, and 317) that respectively target the three tetrapeptides (NPDP/NVDP/NANP) did not provide additional protection against in vivo sporozoite challenge. However, administration of either CIS43, L9, or 317 (but not C-CSP-specific mAbs) to mice that had been immunized with R21, a PfCSP-based virus-like particle vaccine that induces polyclonal antibodies against the repeat region and C-CSP, provided enhanced protection against sporozoite challenge when compared to vaccine or mAbs alone. Collectively, this study shows that while combining mAbs against the repeat and C-terminal regions of PfCSP provide no additional protection in vivo, repeat mAbs do provide increased protection when combined with vaccine-induced polyclonal antibodies. These data should inform the implementation of PfCSP human mAbs alone or following vaccination to prevent malaria infection.

Structural and Energetic Characterization of the Denatured State from the Perspectives of Peptides, the Coil Library, and Intrinsically Disordered Proteins

The α and polyproline II (PPII) basins are the two most populated regions of the Ramachandran map when constructed from the protein coil library, a widely used denatured state model built from the segments of irregular structure found in the Protein Data Bank. This indicates the α and PPII conformations are dominant components of the ensembles of denatured structures that exist in solution for biological proteins, an observation supported in part by structural studies of short, and thus unfolded, peptides. Although intrinsic conformational propensities have been determined experimentally for the common amino acids in short peptides, and estimated from surveys of the protein coil library, the ability of these intrinsic conformational propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs), an increasingly important class of proteins in cell function, has thus far proven elusive to establish. Recently, we demonstrated that the sequence dependence of the mean hydrodynamic size of IDPs in water and the impact of heat on the coil dimensions, provide access to both the sequence dependence and thermodynamic energies that are associated with biases for the α and PPII backbone conformations. Here, we compare results from peptide-based studies of intrinsic conformational propensities and surveys of the protein coil library to those of the sequence-based analysis of heat effects on IDP hydrodynamic size, showing that a common structural and thermodynamic description of the protein denatured state is obtained.

Binding Thermodynamics to Intrinsically Disordered Protein Domains

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

Impact of Heat on Coil Hydrodynamic Size Yields the Energetics of Denatured State Conformational Bias

Conformational equilibria in the protein denatured state have key roles regulating folding, stability, and function. The extent of conformational bias in the protein denatured state under folding conditions, however, has thus far proven elusive to quantify, particularly with regard to its sequence dependence and energetic character. To better understand the structural preferences of the denatured state, we analyzed both the sequence dependence to the mean hydrodynamic size of disordered proteins in water and the impact of heat on the coil dimensions, showing that the sequence dependence and thermodynamic energies associated with intrinsic biases for the α and polyproline II (PPII) backbone conformations can be obtained. Experiments that evaluate how the hydrodynamic size changes with compositional changes in the protein reveal amino acid specific preferences for PPII that are in good quantitative agreement with calorimetry-measured values from unfolded peptides and those inferred by survey of the protein coil library. At temperatures above 25 °C, the denatured state follows the predictions of a PPII-dominant ensemble. Heat effects on coil hydrodynamic size indicate the α bias is comparable to the PPII bias at cold temperatures. Though historically thought to give poor resolution to structural details, the hydrodynamic size of the unfolded state is found to be an effective reporter on the extent of the biases for the α and PPII backbone conformations.

Sequence Reversal Prevents Chain Collapse and Yields Heat-Sensitive Intrinsic Disorder

Sequence patterns of charge, hydrophobicity, hydrogen bonding, and other amino acid physicochemical properties contribute to mechanisms of protein folding, but how sequence composition and patterns influence the conformational dynamics of the denatured state ensemble is not fully understood. To investigate structure-sequence relationships in the denatured state, we reversed the sequence of staphylococcal nuclease and characterized its structure, thermodynamic character, and hydrodynamic radius using circular dichroism spectroscopy, dynamic light scattering, analytical ultracentrifugation, and size-exclusion chromatography as a function of temperature. The macromolecular size of "Retro-nuclease" is highly expanded in solution with characteristics similar to biological intrinsically disordered proteins. In contradistinction to a disordered state, Retro-nuclease exhibits a broad sigmoid transition of its hydrodynamic dimensions as temperature is increased, indicating a thermodynamically controlled compaction. Counterintuitively, the magnitude of these temperature-induced hydrodynamic changes exceed that observed from thermal denaturation of folded unaltered staphylococcal nuclease. Undetectable by calorimetry and intrinsic tryptophan fluorescence, the lack of heat capacity or fluorescence changes throughout the thermal transition indicate canonical hydrophobic collapse did not drive the Retro-nuclease structural transitions. Temperature-dependent circular dichroism spectroscopy performed on Retro-nuclease and computer simulations correlate to temperature sensitivity in the intrinsic sampling of backbone conformations for polyproline II and α-helix. The experimental results indicate a role for sequence direction in mediating the collapse of the polypeptide chain, whereas the simulation trends illustrate the generality of the observed heat effects on disordered protein structure.

Mixed Macromolecular Crowding: A Protein and Solvent Perspective

In the living cell, biomolecules perform their respective functions in the presence of not only one type of macromolecules but rather in the presence of various macromolecules with different shapes and sizes. In this study, we have investigated the effects of five single macromolecular crowding agents, Dextran 6, Dextran 40, Dextran 70, Ficoll 70, and PEG 8000 and their binary mixtures on the modulation in the domain separation of human serum albumin using a Förster resonance energy transfer-based approach and the translational mobility of a small fluorescent probe fluorescein isothiocyanate (FITC) using fluorescence correlation spectroscopy (FCS). Our observations suggest that mixed crowding induces greater cooperativity in the domain movement as compared to the components of the mixtures. Thermodynamic analyses of the same provide evidence of crossovers from enthalpy-based interactions to effects dominated by hard-sphere potential. When compared with those obtained for individual crowders, both domain movements and FITC diffusion studies show significant deviations from ideality, with an ideal solution being considered to be that arising from the sum of the contributions of those obtained in the presence of individual crowding agents. Considering the fact that domain movements are local (on the order of a few angstroms) in nature while translational movements span much larger lengthscales, our results imply that the observed deviation from simple additivity exists at several possible levels or lengthscales in such mixtures. Moreover, the nature and the type of deviation not only depend on the identities of the components of the crowder mixtures but are also influenced by the particular face of the serum protein (either the domain I-II or the domain II-III face) that the crowders interact with, thus providing further insights into the possible existence of microheterogeneities in such solutions.

Structural Stability of the Coiled-Coil Domain of Tumor Susceptibility Gene (TSG)-101

The tumor susceptibility gene-101 coiled coil domain (TSG101cc) is an integral component of the endosomal maturation machinery and cytokinesis, and also interacts with several transcription factors. The TSG101cc has been crystallized as a homotetramer but is known to interact with two of its binding partners as a heterotrimer. To investigate this apparent discrepancy, we examined the solution thermodynamics of the TSG101cc. Here, we use circular dichroism, differential scanning calorimetry, analytical ultracentrifugation, fluorescence, and structural thermodynamic analysis to investigate the structural stability and the unfolding of the TSG101cc. We demonstrate that TSG101cc exists in solution primarily as a tetramer, which unfolds in a two-state manner. Surprisingly, no homodimeric or homotrimeric species were detected. Structural thermodynamic analysis of the homotetrameric structure and comparison with known oligomeric coiled-coils suggests that the TSG101cc homotetramer is comparatively unstable on a per residue basis. Furthermore, the homotrimeric coiled-coil is predicted to be much less stable than the functional heterotrimeric coiled-coil in the endosomal sorting complex required for transport 1 (ESCRT1). These results support a model whereby the tetramer–monomer equilibrium of TSG101 serves as the cellular reservoir of TSG101, which is effectively outcompeted when its binding partners are present and the heteroternary complex can form.

Temperature stability of proteins: Analysis of irreversible denaturation using isothermal calorimetry

The structural stability of proteins has been traditionally studied under conditions in which the folding/unfolding reaction is reversible, since thermodynamic parameters can only be determined under these conditions. Achieving reversibility conditions in temperature stability experiments has often required performing the experiments at acidic pH or other nonphysiological solvent conditions. With the rapid development of protein drugs, the fastest growing segment in the pharmaceutical industry, the need to evaluate protein stability under formulation conditions has acquired renewed urgency. Under formulation conditions and the required high protein concentration (∼100 mg/mL), protein denaturation is irreversible and frequently coupled to aggregation and precipitation. In this article, we examine the thermal denaturation of hen egg white lysozyme (HEWL) under irreversible conditions and concentrations up to 100 mg/mL using several techniques, especially isothermal calorimetry which has been used to measure the enthalpy and kinetics of the unfolding and aggregation/precipitation at 12°C below the transition temperature measured by DSC. At those temperatures the rate of irreversible protein denaturation and aggregation of HEWL is measured to be on the order of 1 day^-1 . Isothermal calorimetry appears a suitable technique to identify buffer formulation conditions that maximize the long term stability of protein drugs.

Intrinsic α helix propensities compact hydrodynamic radii in intrinsically disordered proteins

Proteins that lack tertiary stability under normal conditions, known as intrinsically disordered, exhibit a wide range of biological activities. Molecular descriptions for the biology of intrinsically disordered proteins (IDPs) consequently rely on disordered structural models, which in turn require experiments that assess the origins to structural features observed. For example, while hydrodynamic size is mostly insensitive to sequence composition in chemically denatured proteins, IDPs show strong sequence-specific effects in the hydrodynamic radius (R_h ) when measured under normal conditions. To investigate sequence-modulation of IDP R_h , disordered ensembles generated by a hard sphere collision model modified with a structure-based parameterization of the solution energetics were used to parse the contributions of net charge, main chain dihedral angle bias, and excluded volume on hydrodynamic size. Ensembles for polypeptides 10-35 residues in length were then used to establish power-law scaling relationships for comparison to experimental R_h from 26 IDPs. Results showed the expected outcomes of increased hydrodynamic size from increases in excluded volume and net charge, and compaction from chain-solvent interactions. Chain bias representing intrinsic preferences for α helix and polyproline II (PP_II ), however, modulated R_h with intricate dependence on the simulated propensities. PP_II propensities at levels expected in IDPs correlated with heightened R_h sensitivity to even weak α helix propensities, indicating bias for common (φ, ψ) are important determinants of hydrodynamic size. Moreover, data show that IDP R_h can be predicted from sequence with good accuracy from a small set of physicochemical properties, namely intrinsic conformational propensities and net charge. Proteins 2017; 85:296-311. © 2016 Wiley Periodicals, Inc.

Hydrodynamic Radii of Intrinsically Disordered Proteins Determined from Experimental Polyproline II Propensities

The properties of disordered proteins are thought to depend on intrinsic conformational propensities for polyproline II (PP_II) structure. While intrinsic PP_II propensities have been measured for the common biological amino acids in short peptides, the ability of these experimentally determined propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs) has not been established. Presented here are results from molecular simulations of disordered proteins showing that the hydrodynamic radius (R_h) can be predicted from experimental PP_II propensities with good agreement, even when charge-based considerations are omitted. The simulations demonstrate that R_h and chain propensity for PP_II structure are linked via a simple power-law scaling relationship, which was tested using the experimental R_h of 22 IDPs covering a wide range of peptide lengths, net charge, and sequence composition. Charge effects on R_h were found to be generally weak when compared to PP_II effects on R_h. Results from this study indicate that the hydrodynamic dimensions of IDPs are evidence of considerable sequence-dependent backbone propensities for PP_II structure that qualitatively, if not quantitatively, match conformational propensities measured in peptides.

Molecular models of disordered protein structures are needed to elucidate the functional mechanisms of intrinsically disordered proteins, a class of proteins implicated in many disease pathologies and human health issues. Several studies have measured intrinsic conformational propensities for polyproline II helix, a key structural motif of disordered proteins, in short peptides. Whether or not these experimental polyproline II propensities, which vary by amino acid type, reproduce structural behavior in intrinsically disordered proteins has yet to be demonstrated. Presented here are simulation results showing that polyproline II propensities from short peptides accurately describe sequence-dependent variability in the hydrodynamic dimensions of intrinsically disordered proteins. Good agreement was observed from a simple molecular model even when charge-based considerations were ignored, predicting that global organization of disordered protein structure is strongly dependent on intrinsic conformational propensities and, for many intrinsically disordered proteins, modulated only weakly by coulombic effects.

Protein receptor-independent plasma membrane remodeling by HAMLET: a tumoricidal protein-lipid complex

A central tenet of signal transduction in eukaryotic cells is that extra-cellular ligands activate specific cell surface receptors, which orchestrate downstream responses. This ‘’protein-centric” view is increasingly challenged by evidence for the involvement of specialized membrane domains in signal transduction. Here, we propose that membrane perturbation may serve as an alternative mechanism to activate a conserved cell-death program in cancer cells. This view emerges from the extraordinary manner in which HAMLET (Human Alpha-lactalbumin Made LEthal to Tumor cells) kills a wide range of tumor cells in vitro and demonstrates therapeutic efficacy and selectivity in cancer models and clinical studies. We identify a ‘’receptor independent” transformation of vesicular motifs in model membranes, which is paralleled by gross remodeling of tumor cell membranes. Furthermore, we find that HAMLET accumulates within these de novo membrane conformations and define membrane blebs as cellular compartments for direct interactions of HAMLET with essential target proteins such as the Ras family of GTPases. Finally, we demonstrate lower sensitivity of healthy cell membranes to HAMLET challenge. These features suggest that HAMLET-induced curvature-dependent membrane conformations serve as surrogate receptors for initiating signal transduction cascades, ultimately leading to cell death.

Structural and functional stabilization of protein entities: state-of-the-art

Within the context of biomedicine and pharmaceutical sciences, the issue of (therapeutic) protein stabilization assumes particular relevance. Stabilization of protein and protein-like molecules translates into preservation of both structure and functionality during storage and/or targeting, and such stabilization is mostly attained through establishment of a thermodynamic equilibrium with the (micro)environment. The basic thermodynamic principles that govern protein structural transitions and the interactions of the protein molecule with its (micro)environment are, therefore, tackled in a systematic fashion. Highlights are given to the major classes of (bio)therapeutic molecules, viz. enzymes, recombinant proteins, (macro)peptides, (monoclonal) antibodies and bacteriophages. Modification of the microenvironment of the biomolecule via multipoint covalent attachment onto a solid surface followed by hydrophilic polymer co-immobilization, or physical containment within nanocarriers, are some of the (latest) strategies discussed aiming at full structural and functional stabilization of said biomolecules.

Relationship between protein stability and functional activity in the presence of macromolecular crowding agents alone and in mixture: An insight into stability-activity trade-off

The cellular environment is crowded with different kinds of molecules with varying sizes, shapes and compositions. Most of the experiments studying the nature and behaviour of a protein have been done on the isolated protein in dilute buffer solutions which actually do not imitate the in vivo situation. To understand the consequences of such crowded environment, we investigated the effect of macromolecular crowding on the stability and activity of hen egg white lysozyme. Two crowding agents, dextran 70 and ficoll 70 which have different shapes and composition, have been employed in this study. To mimic the cellular condition from physiological point of view, the effect of mixtures of both the crowding agents has been also studied. The results indicate that owing to volume exclusion, lysozyme is stabilized while its activity decays with the increasing concentration of both the crowders elucidating the hypothesis of stability-activity trade-off. Mixed macromolecular crowding exerts greater effect than the sum of constituent crowding agents (dextran 70 and ficoll 70).

Temperature effects on the hydrodynamic radius of the intrinsically disordered N-terminal region of the p53 protein

Intrinsically disordered proteins (IDPs) are often characterized in terms of the hydrodynamic radius, Rh . The Rh of IDPs are known to depend on fractional proline content and net charge, where increased numbers of proline residues and increased net charge cause larger Rh . Though sequence and charge effects on the Rh of IDPs have been studied, the temperature sensitivity has been noted only briefly. Reported here are Rh measurements in the temperature range of 5-75°C for the intrinsically disordered N-terminal region of the p53 protein, p53(1-93). Of note, the Rh of this protein fragment was highly sensitive to temperature, decreasing from 35 Å at 5°C to 26 Å at 75°C. Computer generated simulations of conformationally dynamic and disordered polypeptide chains were performed to provide a hypothesis for the heat-induced compaction of p53(1-93) structure, which was opposite to the heat-induced increase in Rh observed for a model folded protein. The simulations demonstrated that heat caused Rh to trend toward statistical coil values for both proteins, indicating that the effects of heat on p53(1-93) structure could be interpreted as thermal denaturation. The simulation data also predicted that proline content contributed minimally to the native Rh of p53(1-93), which was confirmed by measuring Rh for a substitution variant that had all 22 proline residues changed for glycine.

Thermodynamic dissection of large-scale domain motions coupled with ligand binding of enzyme I

Domain motions are central to the biological functions of many proteins. The energetics of the motions, however, is often difficult to characterize when motions are coupled with the ligand binding. Here, we determined the thermodynamic parameters of individual domain motions and ligand binding of enzyme I (EI) using strategic domain-deletion mutants that selectively removed particular motions. Upon ligand binding, EI employs two large-scale domain motions, the hinge motion and the swivel motion, to switch between conformational states of distinct domain-domain orientations. Calorimetric analysis of the EI mutants separated the free energy changes of the binding and motions, demonstrating that the unfavorable hinge motion (ΔG = 1.5 kcal mol(-1)) was driven by the favorable swivel motion (ΔG = -5.2 kcal mol(-1)). The large free energy differences could be explained by the physicochemical nature of the domain interfaces associated with the motions; the hinge motion employed much narrower interface than the swivel motion without any hydrogen bonds or salt bridges. The small heat capacity further suggested that the packing of the domain interfaces associated with the hinge motion was less compact than that commonly observed in proteins. Lastly, thermodynamic analysis of phosphorylated EI suggests that the domain motions are regulated by the ligand binding and the phosphorylation states. Taken together, the thermodynamic dissection approach illustrates how multiple motions and ligand binding are energetically connected during the functional cycle of EI.

An intramolecular interaction involving the N terminus of a streptococcal adhesin affects its conformation and adhesive function

BACKGROUND

P1 is an adhesin on the surface of Streptococcus mutans.

RESULTS

Destroying the high affinity interaction between the N and C termini of S. mutans P1 creates a non-adherent phenotype.

CONCLUSION

The N terminus facilitates proper folding, function, and stability within recombinant P1.

SIGNIFICANCE

The relationship between folding, maturation, and cell surface assembly is critical to understanding the P1 mechanism of action. The adhesin P1 is localized on the surface of the oral pathogen Streptococcus mutans and facilitates an interaction with the glycoprotein complex salivary agglutinin that is comprised primarily of the scavenger receptor gp340. Recent crystal structures of P1 display an unusual structure in which the protein folds back upon itself to form an elongated hybrid helical stalk with a globular head at the apex and a globular C-terminal region at the base. The N terminus of P1 has not yet been characterized. In this report we describe the contribution of an interaction between the N-terminal and C-terminal portions of the protein that is required for proper function of P1 on the surface of S. mutans. Utilizing recombinant N-terminal and C-terminal fragments, we employed isothermal titration calorimetry and native gel electrophoresis to demonstrate that these fragments form a high affinity and stable complex in solution. Furthermore, circular dichroism and surface plasmon resonance measurements indicated that the N-terminal fragment contributes to the folding and increases the functionality of the C-terminal fragment in trans. Finally, we utilized circular dichroism, surface plasmon resonance, and differential scanning calorimetry to show that an N-terminal 106-amino acid segment within P1 contributes to the proper folding and function of the full-length recombinant molecule and increases the stability of its elongated hybrid helical stalk.

Discrimination of ligands with different flexibilities resulting from the plasticity of the binding site in tubulin

Tubulin, an α,β heterodimer, has four distinct ligand binding sites (for paclitaxel, peloruside/laulimalide, vinca, and colchicine). The site where colchicine binds is a promising drug target for arresting cell division and has been observed to accommodate compounds that are structurally diverse but possess comparable affinity. This investigation, using two such structurally different ligands as probes (one being colchicine itself and another, TN16), aims to provide insight into the origin of this diverse acceptability to provide a better perspective for the design of novel therapeutic molecules. Thermodynamic measurements reveal interesting interplay between entropy and enthalpy. Although both these parameters are favourable for TN16 binding (ΔH < 0, ΔS > 0), but the magnitude of entropy has the determining role for colchicine binding as its enthalpic component is destabilizing (ΔH > 0, ΔS > 0). Molecular dynamics simulation provides atomistic insight into the mechanism, pointing to the inherent flexibility of the binding pocket that can drastically change its shape depending on the ligand that it accepts. Simulation shows that in the complexed states both the ligands have freedom to move within the binding pocket; colchicine can switch its interactions like a "flying trapeze", whereas TN16 rocks like a "swing cradle", both benefiting entropically, although in two different ways. Additionally, the experimental results with respect to the role of solvation entropy correlate well with the computed difference in the hydration: water molecules associated with the ligands are released upon complexation. The complementary role of van der Waals packing versus flexibility controls the entropy-enthalpy modulations. This analysis provides lessons for the design of new ligands that should balance between the "better fit" and "flexibility"', instead of focusing only on the receptor-ligand interactions.

Thermodynamics of the GTP-GDP-operated conformational switch of selenocysteine-specific translation factor SelB

SelB is a specialized translation factor that binds GTP and GDP and delivers selenocysteyl-tRNA (Sec-tRNA(Sec)) to the ribosome. By analogy to elongation factor Tu (EF-Tu), SelB is expected to control the delivery and release of Sec-tRNA(Sec) to the ribosome by the structural switch between GTP- and GDP-bound conformations. However, crystal structures of SelB suggested a similar domain arrangement in the apo form and GDP- and GTP-bound forms of the factor, raising the question of how SelB can fulfill its delivery function. Here, we studied the thermodynamics of guanine nucleotide binding to SelB by isothermal titration calorimetry in the temperature range between 10 and 25 °C using GTP, GDP, and two nonhydrolyzable GTP analogs, guanosine 5'-O-(γ-thio)triphosphate (GTPγS) and guanosine 5'-(β,γ-imido)-triphosphate (GDPNP). The binding of SelB to either guanine nucleotide is characterized by a large heat capacity change (-621, -467, -235, and -275 cal × mol(-1) × K(-1), with GTP, GTPγS, GDPNP, and GDP, respectively), associated with compensatory changes in binding entropy and enthalpy. Changes in heat capacity indicate a large decrease of the solvent-accessible surface area in SelB, amounting to 43 or 32 amino acids buried upon binding of GTP or GTPγS, respectively, and 15-19 amino acids upon binding GDP or GDPNP. The similarity of the GTP and GDP forms in the crystal structures can be attributed to the use of GDPNP, which appears to induce a structure of SelB that is more similar to the GDP than to the GTP-bound form.

A thermodynamic definition of protein domains

Protein domains are conspicuous structural units in globular proteins, and their identification has been a topic of intense biochemical interest dating back to the earliest crystal structures. Numerous disparate domain identification algorithms have been proposed, all involving some combination of visual intuition and/or structure-based decomposition. Instead, we present a rigorous, thermodynamically-based approach that redefines domains as cooperative chain segments. In greater detail, most small proteins fold with high cooperativity, meaning that the equilibrium population is dominated by completely folded and completely unfolded molecules, with a negligible subpopulation of partially folded intermediates. Here, we redefine structural domains in thermodynamic terms as cooperative folding units, based on m-values, which measure the cooperativity of a protein or its substructures. In our analysis, a domain is equated to a contiguous segment of the folded protein whose m-value is largely unaffected when that segment is excised from its parent structure. Defined in this way, a domain is a self-contained cooperative unit; i.e., its cooperativity depends primarily upon intrasegment interactions, not intersegment interactions. Implementing this concept computationally, the domains in a large representative set of proteins were identified; all exhibit consistency with experimental findings. Specifically, our domain divisions correspond to the experimentally determined equilibrium folding intermediates in a set of nine proteins. The approach was also proofed against a representative set of 71 additional proteins, again with confirmatory results. Our reframed interpretation of a protein domain transforms an indeterminate structural phenomenon into a quantifiable molecular property grounded in solution thermodynamics.

Mutational analysis of m-values as a strategy to identify cold-resistant substructures of the protein ensemble

Characterizing the native ensemble of protein is an important yet difficult objective of structural biology. The structural dynamics of protein macromolecules play key roles in biological function, but the short lifetimes and low population of near-native states of the protein ensemble limit their ability to be studied directly. In part to address such issues, it was shown recently that the cooperative substructures that populate a protein ensemble could be ascertained by NMR methods performed at very cold temperatures. What is presented here is an argument that these same substructures can also be determined by denaturant-induced unfolding studies performed on protein at room temperature. Data supporting this argument are given for Staphylococcal nuclease, chymotrypsin inhibitor 2, and ubiquitin. The observation of an agreement between the thermodynamics of the protein ensemble simulated under very cold temperatures to the apparent sensitivity of the ensemble to chemical denaturants at room temperature also suggests that the overall structural-thermodynamic character of an ensemble is surprisingly robust and preserved even in the presence of strong denaturing conditions.

Strategies for the thermodynamic characterization of linked binding/local folding reactions within the native state application to the LID domain of adenylate kinase from Escherichia coli

Conformational fluctuations in proteins have emerged as an important aspect of biological function, having been linked to processes ranging from molecular recognition and catalysis to allostery and signal transduction. In spite of the realization of their importance, however, the connections between fluctuations and function have largely been empirical, even when they have been quantitative. Part of the problem in understanding the role of fluctuations in function is the fact that the mere existence of fluctuations complicates the interpretation of classic mutagenesis approaches. Namely, mutagenesis, which is typically targeted to an internal position (to elicit an effect), will change the fluctuations as well as the structure of the native state. Decoupling these effects is essential to an unambiguous understanding of the role of fluctuations in function. Here, we use a mutation strategy that targets surface-exposed sites in flexible parts of the molecule for mutation to glycine. Such mutations leave the ground-state structure unaffected. As a result, we can assess the nature of the fluctuations, develop a quantitative model relating fluctuations to function (in this case, molecular recognition), and unambiguously resolve the probabilities of the fluctuating states. We show that when this approach is applied to Escherichia coli adenylate kinase (AK), unique thermodynamic and structural insights are obtained, even when classic mutagenesis approaches targeted to the same region yield ambiguous results.

A thermodynamic analysis of fibrillar polymorphism

We explore the thermodynamic properties of three different fibrils of the peptide hormone glucagon, formed under different salt conditions (glycine, sulfate and NaCl, respectively), and differing considerably in compactness. The three fibrils display a large variation in the specific heat capacity DeltaC(p) determined by isothermal titration calorimetry. Sulfate fibrils show a negative DeltaC(p) expected from a folding reaction, while the DeltaC(p) for glycine fibrils is essentially zero. NaCl fibrils, which are less stable than the other fibrils, have a large and positive C(p). The predicted change in solvent accessible area is not a useful predictor of fibrillar DeltaC(p) unlike the case for globular proteins. We speculate that strong backbone interactions may lead to the unfavorable burial of polar side residues, water and/or charged groups which all can have major influence on the change in C(p). These results highlight differences in the driving forces of native folding and fibril formation.

Thermodynamics of Protein–Ligand Interactions: History, Presence, and Future Aspects

The understanding of molecular recognition processes of small ligands and biological macromolecules requires a complete characterization of the binding energetics and correlation of thermodynamic data with interacting structures involved. A quantitative description of the forces that govern molecular associations requires determination of changes of all thermodynamic parameters, including free energy of binding (deltaG), enthalpy (deltaH), and entropy (deltaS) of binding and the heat capacity change (deltaCp). A close insight into the binding process is of significant and practical interest, since it provides the fundamental know-how for development of structure-based molecular design-strategies. The only direct method to measure the heat change during complex formation at constant temperature is provided by isothermal titration calorimetry (ITC). With this method one binding partner is titrated into a solution containing the interaction partner, thereby generating or absorbing heat. This heat is the direct observable that can be quantified by the calorimeter. The use of ITC has been limited due to the lack of sensitivity, but recent developments in instrument design permit to measure heat effects generated by nanomol (typically 10-100) amounts of reactants. ITC has emerged as the primary tool for characterizing interactions in terms of thermodynamic parameters. Because heat changes occur in almost all chemical and biochemical processes, ITC can be used for numerous applications, e.g., binding studies of antibody-antigen, protein-peptide, protein-protein, enzyme-inhibitor or enzyme-substrate, carbohydrate-protein, DNA-protein (and many more) interactions as well as enzyme kinetics. Under appropriate conditions data analysis from a single experiment yields deltaH, K(B), the stoichiometry (n), deltaG and deltaS of binding. Moreover, ITC experiments performed at different temperatures yield the heat capacity change (deltaCp). The informational content of thermodynamic data is large, and it has been shown that it plays an important role in the elucidation of binding mechanisms and, through the link to structural data, also in rational drug design. In this review we will present a comprehensive overview to ITC by giving some historical background to calorimetry, outline some critical experimental and data analysis aspects, discuss the latest developments, and give three recent examples of studies published with respect to macromolecule-ligand interactions that have utilized ITC technology.

Characterizing the role of ensemble modulation in mutation-induced changes in binding affinity

Protein conformational fluctuations are key contributors to biological function, mediating important processes such as enzyme catalysis, molecular recognition, and allosteric signaling. To better understand the role of conformational fluctuations in substrate/ligand recognition, we analyzed, experimentally and computationally, the binding reaction between an SH3 domain and the recognition peptide of its partner protein. The fluctuations in this SH3 domain were enumerated by using an algorithm based on the hard sphere collision model, and the binding energetics resulting from these fluctuations were calculated using a structure-based energy function parametrized to solvent accessible surface areas. Surprisingly, this simple model reproduced the effects of mutations on the experimentally determined SH3 binding energetics, within the uncertainties of the measurements, indicating that conformational fluctuations in SH3, and in particular the RT loop region, are structurally diverse and are well-approximated by the randomly configured states. The mutated positions in SH3 were distant to the binding site and involved Ala and Gly substitutions of solvent exposed positions in the RT loop. To characterize these fluctuations, we applied principal coordinate analysis to the computed ensembles, uncovering the principal modes of conformational variation. It is shown that the observed differences in binding affinity between each mutant, and thus the apparent coupling between the mutated sites, can be described in terms of the changes in these principal modes. These results indicate that dynamic loops in proteins can populate a broad conformational ensemble and that a quantitative understanding of molecular recognition requires consideration of the entire distribution of states.

Activation of leukemia-associated RhoGEF by Galpha13 with significant conformational rearrangements in the interface

The transient protein-protein interactions induced by guanine nucleotide-dependent conformational changes of G proteins play central roles in G protein-coupled receptor-mediated signaling systems. Leukemia-associated RhoGEF (LARG), a guanine nucleotide exchange factor for Rho, contains an RGS homology (RH) domain and Dbl homology/pleckstrin homology (DH/PH) domains and acts both as a GTPase-activating protein (GAP) and an effector for Galpha(13). However, the molecular mechanism of LARG activation upon Galpha(13) binding is not yet well understood. In this study, we analyzed the Galpha(13)-LARG interaction using cellular and biochemical methods, including a surface plasmon resonance (SPR) analysis. The results obtained using various LARG fragments demonstrated that active Galpha(13) interacts with LARG through the RH domain, DH/PH domains, and C-terminal region. However, an alanine substitution at the RH domain contact position in Galpha(13) resulted in a large decrease in affinity. Thermodynamic analysis revealed that binding of Galpha(13) proceeds with a large negative heat capacity change (DeltaCp degrees ), accompanied by a positive entropy change (DeltaS degrees ). These results likely indicate that the binding of Galpha(13) with the RH domain triggers conformational rearrangements between Galpha(13) and LARG burying an exposed hydrophobic surface to create a large complementary interface, which facilitates complex formation through both GAP and effector interfaces, and activates the RhoGEF. We propose that LARG activation is regulated by an induced-fit mechanism through the GAP interface of Galpha(13).

The solution structure of BMPR-IA reveals a local disorder-to-order transition upon BMP-2 binding

The structure of the extracellular domain of BMP receptor IA was determined in solution by NMR spectroscopy and compared to its structure when bound to its ligand BMP-2. While most parts of the secondary structure are highly conserved between the bound and unbound forms, large conformational rearrangements can be observed in the beta4beta5 loop of BMPR-IA, which is in contact with BMP-2 and harbors the main binding determinants for the BMPR-IA-BMP-2 interaction. In its unbound form, helix alpha1 in BMPR-IA, which is in the center of the binding epitope for BMP-2, is missing. Since BMP-2 also shows conformational changes in the type I receptor epitope upon binding to BMPR-IA, both binding partners pass through an induced fit mechanism to adapt their binding interfaces to a given interaction surface. The inherent flexibility of both partners possibly explains the promiscuous ligand-receptor interaction observed in the BMP protein superfamily.

Evidence that the bZIP domains of the Jun transcription factor bind to DNA as monomers prior to folding and homodimerization

The Jun oncoprotein belongs to the AP1 family of transcription factors that is collectively engaged in diverse cellular processes by virtue of their ability to bind to the promoters of a wide spectrum of genes in a DNA sequence-dependent manner. Here, using isothermal titration calorimetry, we report detailed thermodynamics of the binding of bZIP domain of Jun to synthetic dsDNA oligos containing the TRE and CRE consensus promoter elements. Our data suggest that binding of Jun to both sites occurs with indistinguishable affinities but with distinct thermodynamic signatures comprised of favorable enthalpic contributions accompanied by entropic penalty at physiological temperatures. Furthermore, anomalously large negative heat capacity changes observed provoke a model in which Jun loads onto DNA as unfolded monomers coupled with subsequent folding and homodimerization upon association. Taken together, our data provide novel insights into the energetics of a key protein-DNA interaction pertinent to cellular signaling and cancer. Our study underscores the notion that the folding and dimerization of transcription factors upon association with DNA may be a more general mechanism employed in protein-DNA interactions and that the conventional school of thought may need to be re-evaluated.

Sequence specific and high affinity recognition of 5'-ACGCGT-3' by rationally designed pyrrole-imidazole H-pin polyamides: thermodynamic and structural studies

Imidazole (Im) and Pyrrole (Py)-containing polyamides that can form stacked dimers can be programmed to target specific sequences in the minor groove of DNA and control gene expression. Even though various designs of polyamides have been thoroughly investigated for DNA sequence recognition, the use of H-pin polyamides (covalently cross-linked polyamides) has not received as much attention. Therefore, experiments were designed to systematically investigate the DNA recognition properties of two symmetrical H-pin polyamides composed of PyImPyIm (5) or f-ImPyIm (3e, f=formamido) tethered with an ethylene glycol linker. These compounds were created to recognize the cognate 5'-ACGCGT-3' through an overlapped and staggered binding motif, respectively. Results from DNaseI footprinting, thermal denaturation, circular dichroism, surface plasmon resonance and isothermal titration microcalorimetry studies demonstrated that both H-pin polyamides bound with higher affinity than their respective monomers. The binding affinity of formamido-containing H-pin 3e was more than a hundred times greater than that for the tetraamide H-pin 5, demonstrating the importance of having a formamido group and the staggered motif in enhancing affinity. However, compared to H-pin 3e, tetraamide H-pin 5 demonstrated superior binding preference for the cognate sequence over its non-cognates, ACCGGT and AAATTT. Data from SPR experiments yielded binding constants of 1.6x10(8)M(-1) and 2.0x10(10)M(-1) for PyImPyIm H-pin 5 and f-ImPyIm H-pin 3e, respectively. Both H-pins bound with significantly higher affinity (ca. 100-fold) than their corresponding unlinked PyImPyIm 4 and f-ImPyIm 2 counterparts. ITC analyses revealed modest enthalpies of reactions at 298 K (DeltaH of -3.3 and -1.0 kcal mol(-1) for 5 and 3e, respectively), indicating these were entropic-driven interactions. The heat capacities (DeltaC(p)) were determined to be -116 and -499 cal mol(-1)K(-1), respectively. These results are in general agreement with DeltaC(p) values determined from changes in the solvent accessible surface areas using complexes of the H-pins bound to (5'-CCACGCGTGG)(2). According to the models, the H-pins fit snugly in the minor groove and the linker comfortably holds both polyamide portions in place, with the oxygen atoms pointing into the solvent. In summary, the H-pin polyamide provides an important molecular design motif for the discovery of future generations of programmable small molecules capable of binding to target DNA sequences with high affinity and selectivity.

Exploring the impact of polyproline II (PII) conformational bias on the binding of peptides to the SEM-5 SH3 domain

The left-handed polyproline II helical structure (P(II)) is observed to be a dominant conformation in the disordered states of protein and small polypeptide chains, even when no prolines are present in the sequence. Recently, in work by Ferreon and Hilser, the energetics associated with Ala and Gly substitutions at a surface exposed proline site were determined calorimetrically by measuring the binding energetics of Sos peptide variants to the C-terminal Src Homology 3 domain of SEM-5. The results were interpreted as a significant conformational bias toward the bound conformation (i.e., P(II)), even when the ligand is unbound. That study was not able to determine, however, whether the conformational bias of the peptides could be explained in terms other than that of a P(II) preference. Here, we test, using a computer algorithm based on the hard sphere collision (HSC) model, the notion of whether a bias in the unbound states of the peptide ligands is specific for the P(II) conformation, or if a bias to any other region of (phi, psi) space can also result in the same observed binding energetics. The results of these computer simulations indicate that, of the regions of (phi, psi) modeled for bias in the small peptides, only the bias to the P(II) conformation, and at rates of bias similar to the experimentally observed rates, quantitatively reproduced the experimental binding energetics.

Estimating free-energy barrier heights for an ultrafast folding protein from calorimetric and kinetic data

Differential scanning calorimetry was used to measure the temperature dependence of the absolute heat capacity of the 35-residue subdomain of the villin headpiece, a protein that folds in 5 mus and is therefore assumed to have a small free-energy barrier separating folded and unfolded states. To obtain an estimate of the barrier height from the calorimetric data, two models, a variable-barrier model and an Ising-like model, were used to fit the heat capacity in excess of the folded state over the temperature range 15-125 degrees C. The variable-barrier model is based on an empirical mathematical form for the density of states, with four adjustable parameters and the enthalpy (H) as a reaction coordinate. The Ising-like model is based on the inter-residue contact map of the X-ray structure with exact enumeration of approximately 10(5) possible conformations, with two adjustable parameters in the partition function, and either the fraction of native contacts (Q) or the number of ordered residues (P) as reaction coordinates. The variable-barrier model provides an excellent fit to the data and yields a barrier height at the folding temperature ranging from 0.4 to 1.1 kcal mol(-1), while the Ising-like model provides a less good fit and yields barrier heights of 2.3 +/- 0.1 kcal mol(-1) and 2.1 +/- 0.1 kcal mol(-1) for the Q and P reaction coordinates, respectively. In both models, the barrier to folding increases with increasing temperature. Assuming a sufficiently large activation energy for diffusion on the free-energy surfaces, both models are consistent with the observation of a temperature-independent folding rate in previously published laser temperature-jump experiments. Analysis of this kinetic data, using an approximate form for the pre-exponential factor of Kramers theory and the 70 ns relaxation time for the fast phase that precedes the unfolding/refolding relaxation to determine the diffusion coefficient, results in a barrier height of 1.6 +/- 0.3 kcal mol-1 for an unspecified reaction coordinate. Although no independent test of the validity of the H, Q, or P reaction coordinates is given, the barrier-height estimates obtained with the three reaction coordinates are in quite good agreement with the value derived from a Kramers analysis of the kinetics that makes no assumptions about the reaction coordinate. However, the higher estimates obtained using Q or P appear more consistent with the finding of barrier-crossing kinetics of a villin mutant that folds in 700 ns, corresponding to a 1.3 kcal mol-1 reduction in the folding barrier relative to wild-type. All of the results suggest that the free-energy barrier to folding is sufficiently low that it should be possible to engineer this protein or find solution conditions that would eliminate the barrier to create the "downhill" folding scenario of Wolynes and Onuchic.

Modifying the N-terminus of polyamides: PyImPyIm has improved sequence specificity over f-ImPyIm

Seven N-terminus modified derivatives of a previously published minor-groove binding polyamide (f-ImPyIm, 1) were synthesized and the biochemical and biophysical chemistry evaluated. These compounds were synthesized with the aim of attaining a higher level of sequence selectivity over f-ImPyIm (1), a previously published strong minor-groove binder. Two compounds possessing a furan or a benzofuran moiety at the N-terminus showed a footprint of 0.5microM at the cognate ACGCGT site (determined by DNase I footprinting); however, the specificity of these compounds was not improved. In contrast, PyImPyIm (4) produced a footprint of 0.5microM but showed a superior specificity using the same technique. When evaluated by thermal melting experiments and circular dichroism using ACGCGT and the non-cognate AAATTT sequence, all compounds were shown to bind in the minor-groove of DNA and stabilize the cognate sequence much better than the non-cognate (except for the non-amido-compound that did not bind either sequence, as expected). PyImPyIm (4) was interesting as the DeltaT(m) for this compound was only 4 degrees C but the footprint was very selective. No binding was observed for this compound with a third DNA (non-cognate, ACCGGT). ITC studies on compound 4 showed exothermic binding with ACGCGT and no heat change was observed for titrating the compound to the other two DNA sequences. The heat capacity (DeltaC(p)) of the PIPI/ACGCGT complex calculated from the hydrophobic interactions and SASA calculations was comparable to the experimental value obtained from ITC (-146calmol(-1)K(-1)). SPR results provided confirmation of the sequence specificity of PyImPyIm (4), with a K(eq) value determined to be 7.1x10(6) M(-1) for the cognate sequence and no observable binding to AAATTT and ACCGGT. Molecular dynamic simulations affirmed that PyImPyIm (4) binds as a dimer in an overlapped conformation, and it fits snugly in the minor-groove of the ACGCGT oligonucleotide. PyImPyIm (4) is an especially interesting molecule, because although the binding affinity is slightly reduced, the specificity with respect to f-ImPyIm (1) is significantly improved.

Coupling of folding and DNA-binding in the bZIP domains of Jun-Fos heterodimeric transcription factor

In response to mitogenic stimuli, the heterodimeric transcription factor Jun-Fos binds to the promoters of a diverse array of genes involved in critical cellular responses such as cell growth and proliferation, cell cycle regulation, embryogenic development and cancer. In so doing, Jun-Fos heterodimer regulates gene expression central to physiology and pathology of the cell in a specific and timely manner. Here, using the technique of isothermal titration calorimetry (ITC), we report detailed thermodynamics of the bZIP domains of Jun-Fos heterodimer to synthetic dsDNA oligos containing the TRE and CRE consensus promoter elements. Our data suggest that binding of the bZIP domains to both TRE and CRE is under enthalpic control and accompanied by entropic penalty at physiological temperatures. Although the bZIP domains bind to both TRE and CRE with very similar affinities, the enthalpic contributions to the free energy of binding to CRE are more favorable than TRE, while the entropic penalty to the free energy of binding to TRE is smaller than CRE. Despite such differences in their thermodynamic signatures, enthalpy and entropy of binding of the bZIP domains to both TRE and CRE are highly temperature-dependent and largely compensate each other resulting in negligible effect of temperature on the free energy of binding. From the plot of enthalpy change versus temperature, the magnitude of heat capacity change determined is much larger than that expected from the direct association of bZIP domains with DNA. This observation is interpreted to suggest that the basic regions in the bZIP domains are largely unstructured in the absence of DNA and only become structured upon interaction with DNA in a coupled folding and binding manner. Our new findings are rationalized in the context of 3D structural models of bZIP domains of Jun-Fos heterodimer in complex with dsDNA oligos containing the TRE and CRE consensus sequences. Taken together, our study demonstrates that enthalpy is the major driving force for a key protein-DNA interaction pertinent to cellular signaling and that protein-DNA interactions with similar binding affinities may be accompanied by differential thermodynamic signatures. Our data corroborate the notion that the DNA-induced protein structural changes are a general feature of the bZIP family of transcription factors.

Ligand effects on the protein ensemble: unifying the descriptions of ligand binding, local conformational fluctuations, and protein stability

Detailed description of the structural and physical basis of allostery, cooperativity, and other manifestations of long-range communication between binding sites in proteins remains elusive. Here we describe an ensemble-based structural-thermodynamic model capable of treating explicitly the coupling between ligand binding reactions, local fluctuations in structure, and global conformational transitions. The H(+) binding reactions of staphylococcal nuclease and the effects of pH on its stability were used to illustrate the properties of proteins that can be described quantitatively with this model. Each microstate in the native ensemble was modeled to have dual structural character; some regions were treated as folded and retained the same atomic geometry as in the crystallographic structure while other regions were treated thermodynamically as if they were unfolded. Two sets of pK(a) values were used to describe the affinity of each H(+) binding site. One set, calculated with a standard continuum electrostatics method, describes H(+) binding to sites in folded parts of the protein. A second set of pK(a) values, obtained from model compounds in water, was used to describe H(+) binding to sites in unfolded regions. An empirical free energy function, parameterized to reproduce folding thermodynamics measured by differential scanning calorimetry, was used to calculate the probability of each microstate. The effects of pH on the distribution of microstates were determined by the H(+) binding properties of each microstate. The validity of the calculations was established by comparison with a number of different experimental observables.

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.

Functional residues serve a dominant role in mediating the cooperativity of the protein ensemble

Conformational fluctuations in proteins have emerged as a potentially important aspect of biological function, although the precise relationship and the implications have yet to be fully explored. Numerous studies have reported that the binding of ligand can influence fluctuations. However, the role of the binding site in mediating these fluctuations is not known. Of particular interest is whether in addition to serving as structural scaffolds for recognition and catalysis, active-site residues may also play a role in modulating the cooperative network. To address this question, we employ an experimentally validated ensemble-based description of proteins to elucidate the extent to which perturbations at different sites can influence the cooperative network in the protein. Applying this method to a database of test proteins, it is found statistically that binding sites are located in regions most able to affect the cooperative network, even for cooperative interactions between residues distant to the binding sites. This indicates that the conformational manifold under native conditions is determined by the network of cooperative interactions within the protein and suggests that proteins have evolved to use these conformational fluctuations in carrying out their functions. Furthermore, because the energetic coupling pattern calculated for each protein is robust and relatively insensitive to sequence, these studies further suggest that binding sites evolved in regions of the protein that are inherently poised to take advantage of the fluctuations in the native structure.

Thermodynamics of proteins in unusual environments

Some aspects of protein thermodynamics in unconventional environments are addressed and discussed. Aqueous medium, especially dilute solution is the 'usual' ambient, which mediates all the interactions between protein and nearby molecules. When the water content is low, the surroundings may be considered 'unusual', exerting new stresses on the protein molecule and demanding different responses and property changes. The unusual systems considered in this article are low-water protein environments, including nearly dry state powders, organic solvent dispersions and reverse micelles' inclusions. The changes of hydration experienced by the protein after immobilization on solid supports are emphasized with respect to the free bulk solution state. Finally, the aqueous medium altered by water connectivity perturbing agents (polysaccharides) or in macromolecular crowding conditions (in the presence of polyols) is also considered as highly not ideal protein environments. The different responses elicited by the protein under the stress induced by drastic surrounding alterations may give insights for the controlled exploitation of the protein's biological and thermodynamic properties.

Binding-linked protonation of a DNA minor-groove agent

The energetics for binding of a diphenyl diamidine antitrypanosomal agent CGP 40215A to DNA have been studied by spectroscopy, isothermal titration calorimetry, and surface plasmon resonance biosensor methods. Both amidines are positively charged under experimental conditions, but the linking group for the two phenyl amidines has a pK(a) of 6.3 that is susceptible to a protonation process. Spectroscopic studies indicate an increase of 2.7 pK(a) units in the linking group when the compound binds to an A/T minor-groove site. Calorimetric titrations in different buffers and pH conditions support the proton-linkage process and are in a good agreement with spectroscopic titrations. The two methods established a proton-uptake profile as a function of pH. The exothermic enthalpy of complex formation varies with different pH conditions. The observed binding enthalpy increases as a function of temperature indicating a negative heat capacity change that is typical for DNA minor-groove binders. Solvent accessible surface area calculations suggest that surface burial accounts for about one-half of the observed intrinsic negative heat capacity change. Biosensor and calorimetric experiments indicate that the binding affinities vary with pH values and salt concentrations due to protonation and electrostatic interactions. The surface plasmon resonance binding studies indicate that the charge density per phosphate in DNA hairpins is smaller than that in polymers. Energetic contributions from different factors were also estimated for the ligand/DNA complex.

Loosely packed papain prosegment displays inhibitory activity

Most protease prosegments are co-synthesized at the N-termini of cysteine proteases and are involved in folding assistance, inhibition, and activation of their mature enzymes. By using circular dichroism, UV-difference and fluorescence spectroscopies, we studied the thermal unfolding of papain prosegment. The transition seems to be two-state and reversible, with an unfolded state prone to aggregation. Unfolding thermodynamic parameters obtained show low values both for deltaH(Tm) and deltaCp(U), indicative of a loosely packed three-dimensional conformation for the prosegment at near-neutral pH conditions. In spite of these results, fluorescence experiments demonstrate that papain prosegment is able to recognize and inhibit its cognate protease. An acid medium induces a molten globule-like state without intermediates, which in turn undergoes an irreversible thermal unfolding. Our results suggest that papain prosegment has a high degree of conformational flexibility, with the ability to form not only a molten globule-like structure in activating conditions, but also requiring an induced fit in order to be functional as inhibitor.

Local conformational fluctuations can modulate the coupling between proton binding and global structural transitions in proteins

Local conformational fluctuations in proteins can affect the coupling between ligand binding and global structural transitions. This finding was established by monitoring quantitatively how the population distribution in the ensemble of microstates of staphylococcal nuclease was affected by proton binding. Analysis of acid unfolding and proton-binding data with an ensemble-based model suggests that local fluctuations: (i) can be effective modulators of ligand-binding affinities, (ii) are important determinants of the cooperativity of ligand-driven global structural transitions, and (iii) are well represented thermodynamically as local unfolding processes. These studies illustrate how an ensemble-based description of proteins can be used to describe quantitatively the interdependence of local conformational fluctuations, ligand-binding processes, and global structural transitions. This level of understanding of the relationship between conformation, energy, and dynamics is required for a detailed mechanistic understanding of allostery, cooperativity, and other complex functional and regulatory properties of macromolecules.

Counteracting Osmolyte Trimethylamine N-Oxide Destabilizes Proteins at pH below Its pK

Earlier studies have reported that trimethylamine N-oxide (TMAO), a naturally occurring osmolyte, is a universal stabilizer of proteins because it folds unstructured proteins and counteracts the deleterious effects of urea and salts on the structure and function of proteins. This conclusion has been reached from the studies of the effect of TMAO on proteins in the pH range 6.0-8.0. In this pH range TMAO is almost neutral (zwitterionic form), for it has a pK(a) of 4.66 +/- 0.10. We have asked the question of whether the effect of TMAO on protein stability is pH-dependent. To answer this question we have carried out thermal denaturation studies of lysozyme, ribonuclease-A, and apo-alpha-lactalbumin in the presence of various TMAO concentrations at different pH values above and below the pK(a) of TMAO. The main conclusion of this study is that near room temperature TMAO destabilizes proteins at pH values below its pK(a), whereas it stabilizes proteins at pH values above its pK(a). This conclusion was reached by determining the T(m) (midpoint of denaturation), delta H(m) (denaturational enthalpy change at T(m)), delta C(p) (constant pressure heat capacity change), and delta G(D) degrees (denaturational Gibbs energy change at 25 degrees C) of proteins in the presence of different TMAO concentrations. Other conclusions of this study are that T(m) and delta G(D) degrees depend on TMAO concentration at each pH value and that delta H(m) and the delta C(p) are not significantly changed in presence of TMAO.

The salt-dependence of a protein-ligand interaction: ion-protein binding energetics

Using the binding of a nucleotide inhibitor (guanosine-3'-monophosphate) to a ribonuclease (ribonuclease Sa) as a model system, we show that the salt-dependence of the interaction arises due to specific ion binding at the site of nucleotide binding. The presence of specific ion-protein binding is concluded from a combination of differential scanning calorimetry and NMR data. Isothermal titration calorimetry data are then fit to determine the energetic profile (enthalpy, entropy, and heat capacity) for both the ion-protein and nucleotide-protein interactions. The results provide insight into the energetics of charge-charge interactions, and have implications for the interpretation of an observed salt-dependence. Further, the presence of specific ion-binding leads to a system behavior as a function of temperature that is drastically different from that predicted from Poisson-Boltzmann calculations.

Protein aggregation and its inhibition in biopharmaceutics

Protein aggregation is arguably the most common and troubling manifestation of protein instability, encountered in almost all stages of protein drug development. Protein aggregation, along with other physical and/or chemical instabilities of proteins, remains to be one of the major road barriers hindering rapid commercialization of potential protein drug candidates. Although a variety of methods have been used/designed to prevent/inhibit protein aggregation, the end results are often unsatisfactory for many proteins. The limited success is partly due to our lack of a clear understanding of the protein aggregation process. This article intends to discuss protein aggregation and its related mechanisms, methods characterizing protein aggregation, factors affecting protein aggregation, and possible venues in aggregation prevention/inhibition in various stages of protein drug development.