The meaning of hydrophobicity

Physicochemical descriptors to discriminate protein-protein interactions in permanent and transient complexes selected by means of machine learning algorithms

Analyzing protein-protein interactions at the atomic level is critical for our understanding of the principles governing the interactions involved in protein-protein recognition. For this purpose, descriptors explaining the nature of different protein-protein complexes are desirable. In this work, the authors introduced Epic Protein Interface Classification as a framework handling the preparation, processing, and analysis of protein-protein complexes for classification with machine learning algorithms. We applied four different machine learning algorithms: Support Vector Machines, C4.5 Decision Trees, K Nearest Neighbors, and Naïve Bayes algorithm in combination with three feature selection methods, Filter (Relief F), Wrapper, and Genetic Algorithms, to extract discriminating features from the protein-protein complexes. To compare protein-protein complexes to each other, the authors represented the physicochemical characteristics of their interfaces in four different ways, using two different atomic contact vectors, DrugScore pair potential vectors and SFCscore descriptor vectors. We classified two different datasets: (A) 172 protein-protein complexes comprising 96 monomers, forming contacts enforced by the crystallographic packing environment (crystal contacts), and 76 biologically functional homodimer complexes; (B) 345 protein-protein complexes containing 147 permanent complexes and 198 transient complexes. We were able to classify up to 94.8% of the packing enforced/functional and up to 93.6% of the permanent/transient complexes correctly. Furthermore, we were able to extract relevant features from the different protein-protein complexes and introduce an approach for scoring the importance of the extracted features.

Buccal bioadhesive drug delivery--a promising option for orally less efficient drugs

Rapid developments in the field of molecular biology and gene technology resulted in generation of many macromolecular drugs including peptides, proteins, polysaccharides and nucleic acids in great number possessing superior pharmacological efficacy with site specificity and devoid of untoward and toxic effects. However, the main impediment for the oral delivery of these drugs as potential therapeutic agents is their extensive presystemic metabolism, instability in acidic environment resulting into inadequate and erratic oral absorption. Parenteral route of administration is the only established route that overcomes all these drawbacks associated with these orally less/inefficient drugs. But, these formulations are costly, have least patient compliance, require repeated administration, in addition to the other hazardous effects associated with this route. Over the last few decades' pharmaceutical scientists throughout the world are trying to explore transdermal and transmucosal routes as an alternative to injections. Among the various transmucosal sites available, mucosa of the buccal cavity was found to be the most convenient and easily accessible site for the delivery of therapeutic agents for both local and systemic delivery as retentive dosage forms, because it has expanse of smooth muscle which is relatively immobile, abundant vascularization, rapid recovery time after exposure to stress and the near absence of langerhans cells. Direct access to the systemic circulation through the internal jugular vein bypasses drugs from the hepatic first pass metabolism leading to high bioavailability. Further, these dosage forms are self-administrable, cheap and have superior patient compliance. Developing a dosage form with the optimum pharmacokinetics is a promising area for continued research as it is enormously important and intellectually challenging. With the right dosage form design, local environment of the mucosa can be controlled and manipulated in order to optimize the rate of drug dissolution and permeation. A rational approach to dosage form design requires a complete understanding of the physicochemical and biopharmaceutical properties of the drug and excipients. Advances in experimental and computational methodologies will be helpful in shortening the processing time from formulation design to clinical use. This paper aims to review the developments in the buccal adhesive drug delivery systems to provide basic principles to the young scientists, which will be useful to circumvent the difficulties associated with the formulation design.

The crystal structure of a hyperthermostable subfamily II isocitrate dehydrogenase from Thermotoga maritima

Isocitrate dehydrogenase (IDH) from the hyperthermophile Thermotoga maritima (TmIDH) catalyses NADP+- and metal-dependent oxidative decarboxylation of isocitrate to alpha-ketoglutarate. It belongs to the beta-decarboxylating dehydrogenase family and is the only hyperthermostable IDH identified within subfamily II. Furthermore, it is the only IDH that has been characterized as both dimeric and tetrameric in solution. We solved the crystal structure of the dimeric apo form of TmIDH at 2.2 A. The R-factor of the refined model was 18.5% (R(free) 22.4%). The conformation of the TmIDH structure was open and showed a domain rotation of 25-30 degrees compared with closed IDHs. The separate domains were found to be homologous to those of the mesophilic mammalian IDHs of subfamily II and were subjected to a comparative analysis in order to find differences that could explain the large difference in thermostability. Mutational studies revealed that stabilization of the N- and C-termini via long-range electrostatic interactions were important for the higher thermostability of TmIDH. Moreover, the number of intra- and intersubunit ion pairs was higher and the ionic networks were larger compared with the mesophilic IDHs. Other factors likely to confer higher stability in TmIDH were a less hydrophobic and more charged accessible surface, a more hydrophobic subunit interface, more hydrogen bonds per residue and a few loop deletions. The residues responsible for the binding of isocitrate and NADP+ were found to be highly conserved between TmIDH and the mammalian IDHs and it is likely that the reaction mechanism is the same.

Genetic algorithms as a tool for helix design – computational and experimental studies on prion protein helix 1

Evolutionary computing is a general optimization mechanism successfully implemented for a variety of numeric problems in a variety of fields, including structural biology. We here present an evolutionary approach to optimize helix stability in peptides and proteins employing the AGADIR energy function for helix stability as scoring function. With the ability to apply masks determining positions, which are to remain constant or fixed to a certain class of amino acids, our algorithm is capable of developing stable helical scaffolds containing a wide variety of structural and functional amino acid patterns. The algorithm showed good convergence behaviour in all tested cases and can be parameterized in a wide variety of ways. We have applied our algorithm for the optimization of the stability of prion protein helix 1, a structural element of the prion protein which is thought to play a crucial role in the conformational transition from the cellular to the pathogenic form of the prion protein, and which therefore poses an interesting target for pharmacological as well as genetic engineering approaches to counter the as of yet uncurable prion diseases. NMR spectroscopic investigations of selected stabilizing and destabilizing mutations found by our algorithm could demonstrate its ability to create stabilized variants of secondary structure elements.

Backbone nuclear relaxation characteristics and calorimetric investigation of the human Grb7-SH2/erbB2 peptide complex

Grb7 is a member of the Grb7 family of proteins, which also includes Grb10 and Grb14. All three proteins have been found to be overexpressed in certain cancers and cancer cell lines. In particular, Grb7 (along with the receptor tyrosine kinase erbB2) is overexpressed in 20%-30% of breast cancers. Grb7 binds to erbB2 and may be involved in cell signaling pathways that promote the formation of metastases and inflammatory responses. In a prior study, we reported the solution structure of the Grb7-SH2/erbB2 peptide complex. In this study, T(1), T(2), and steady-state NOE measurements were performed on the Grb7-SH2 domain, and the backbone relaxation behavior of the domain is discussed with respect to the potential function of an insert region present in all three members of this protein family. Isothermal titration calorimetry (ITC) studies were completed measuring the thermodynamic parameters of the binding of a 10-residue phosphorylated peptide representative of erbB2 to the SH2 domain. These measurements are compared to calorimetric studies performed on other SH2 domain/phosphorylated peptide complexes available in the literature.

Long dynamics simulations of proteins using atomistic force fields and a continuum representation of solvent effects: calculation of structural and dynamic properties

Long dynamics simulations were carried out on the B1 immunoglobulin-binding domain of streptococcal protein G (ProtG) and bovine pancreatic trypsin inhibitor (BPTI) using atomistic descriptions of the proteins and a continuum representation of solvent effects. To mimic frictional and random collision effects, Langevin dynamics (LD) were used. The main goal of the calculations was to explore the stability of tens-of-nanosecond trajectories as generated by this molecular mechanics approximation and to analyze in detail structural and dynamical properties. Conformational fluctuations, order parameters, cross correlation matrices, residue solvent accessibilities, pKa values of titratable groups, and hydrogen-bonding (HB) patterns were calculated from all of the trajectories and compared with available experimental data. The simulations comprised over 40 ns per trajectory for ProtG and over 30 ns per trajectory for BPTI. For comparison, explicit water molecular dynamics simulations (EW/MD) of 3 ns and 4 ns, respectively, were also carried out. Two continuum simulations were performed on each protein using the CHARMM program, one with the all-atom PAR22 representation of the protein force field (here referred to as PAR22/LD simulations) and the other with the modifications introduced by the recently developed CMAP potential (CMAP/LD simulations). The explicit solvent simulations were performed with PAR22 only. Solvent effects are described by a continuum model based on screened Coulomb potentials (SCP) reported earlier, i.e., the SCP-based implicit solvent model (SCP-ISM). For ProtG, both the PAR22/LD and the CMAP/LD 40-ns trajectories were stable, yielding C(alpha) root mean square deviations (RMSD) of about 1.0 and 0.8 A respectively along the entire simulation time, compared to 0.8 A for the EW/MD simulation. For BPTI, only the CMAP/LD trajectory was stable for the entire 30-ns simulation, with a C(alpha) RMSD of approximately 1.4 A, while the PAR22/LD trajectory became unstable early in the simulation, reaching a C(alpha) RMSD of about 2.7 A and remaining at this value until the end of the simulation; the C(alpha) RMSD of the EW/MD simulation was about 1.5 A. The source of the instabilities of the BPTI trajectories in the PAR22/LD simulations was explored by an analysis of the backbone torsion angles. To further validate the findings from this analysis of BPTI, a 35-ns SCP-ISM simulation of Ubiquitin (Ubq) was carried out. For this protein, the CMAP/LD simulation was stable for the entire simulation time (C(alpha) RMSD of approximately 1.0 A), while the PAR22/LD trajectory showed a trend similar to that in BPTI, reaching a C(alpha) RMSD of approximately 1.5 A at 7 ns. All the calculated properties were found to be in agreement with the corresponding experimental values, although local deviations were also observed. HB patterns were also well reproduced by all the continuum solvent simulations with the exception of solvent-exposed side chain-side chain (sc-sc) HB in ProtG, where several of the HB interactions observed in the crystal structure and in the EW/MD simulation were lost. The overall analysis reported in this work suggests that the combination of an atomistic representation of a protein with a CMAP/CHARMM force field and a continuum representation of solvent effects such as the SCP-ISM provides a good description of structural and dynamic properties obtained from long computer simulations. Although the SCP-ISM simulations (CMAP/LD) reported here were shown to be stable and the properties well reproduced, further refinement is needed to attain a level of accuracy suitable for more challenging biological applications, particularly the study of protein-protein interactions.

Thermodynamic anomalies in a lattice model of water: Solvation properties

We investigate a lattice-fluid model of water, defined on a three-dimensional body-centered-cubic lattice. Model molecules possess a tetrahedral symmetry, with four equivalent bonding arms. The model is similar to the one proposed by Roberts and Debenedetti [J. Chem. Phys. 105, 658 (1996)], simplified by removing distinction between "donors" and "acceptors." We focus on the solvation properties, mainly as far as an ideally inert (hydrophobic) solute is concerned. As in our previous analysis, devoted to neat water [J. Chem. Phys. 121, 11856 (2004)], we make use of a generalized first-order approximation on a tetrahedral cluster. We show that the model exhibits quite a coherent picture of water thermodynamics, reproducing qualitatively several anomalous properties observed both in pure water and in solutions of hydrophobic solutes. As far as supercooled liquid water is concerned, the model is consistent with the second critical-point scenario.

Modeling Water, the Hydrophobic Effect, and Ion Solvation

Water plays a central role in the structures and properties of biomolecules--proteins, nucleic acids, and membranes--and in their interactions with ligands and drugs. Over the past half century, our understanding of water has been advanced significantly owing to theoretical and computational modeling. However, like the blind men and the elephant, different models describe different aspects of water's behavior. The trend in water modeling has been toward finer-scale properties and increasing structural detail, at increasing computational expense. Recently, our labs and others have moved in the opposite direction, toward simpler physical models, focusing on more global properties-water's thermodynamics, phase diagram, and solvation properties, for example-and toward less computational expense. Simplified models can guide a better understanding of water in ways that complement what we learn from more complex models. One ultimate goal is more tractable models for computer simulations of biomolecules. This review gives a perspective from simple models on how the physical properties of water-as a pure liquid and as a solvent-derive from the geometric and hydrogen bonding properties of water.

Hydration of an apolar solute in a two-dimensional waterlike lattice fluid

In a previous work, we investigated a two-dimensional lattice-fluid model, displaying some waterlike thermodynamic anomalies. The model, defined on a triangular lattice, is now extended to aqueous solutions with apolar species. Water molecules are of the "Mercedes Benz" type, i.e., they possess a D3 (equilateral triangle) symmetry, with three equivalent bonding arms. Bond formation depends both on orientation and local density. The insertion of inert molecules displays typical signatures of hydrophobic hydration: large positive transfer free energy, large negative transfer entropy (at low temperature), strong temperature dependence of the transfer enthalpy and entropy, i.e., large (positive) transfer heat capacity. Model properties are derived by a generalized first order approximation on a triangle cluster.

A new amphipathy scale

Two new amphipathy scales elaborated from molecular dynamics data are presented. Their applications contribute for the identification of the hydrophobic or hydrophilic regions in proteins solely from the primary structure. The new amphipathy coefficients (AC) reflect the side chain/solvent molecules configurational energies. A polar (water) and an apolar solvent, CCl4, were used resulting in the two ACwater and ACCCl4 scales. These solvents were chosen to simulate the aqueous phases and the transmembrane ambients of cellular membranes where the membrane proteins act. The new amphipathy scales were compared with some previous scales determined by different methods, which were also compared between them, indicating more than 90% of the correlation coefficients are less than 0.9: the scales are strictly dependent on the methodologies used in their determination. The ACCCl4 scale is related with the size of side chain amino acids while ACwater is related with the hydrophobicity of side chain amino acids. The quality of the scales was confirmed by an example of application where ACwater was able to identify correctly the transmembrane, hydrophobic regions of a membrane protein. These results also indicate that water is an important factor responsible for the tertiary structure of membrane proteins.

Evaluation of methods for measuring amino acid hydrophobicities and interactions

The concept of hydrophobicity has been addressed by researchers in all aspects of science, particularly in the fields of biology and chemistry. Over the past several decades, the study of the hydrophobicity of biomolecules, particularly amino acids has resulted in the development of a variety of hydrophobicity scales. In this review, we discuss the various methods of measuring amino acid hydrophobicity and provide explanations for the wide range of rankings that exist among these published scales. A discussion of the literature on amino acid interactions is also presented. Only a surprisingly small number of papers exist in this rather important area of research; measuring pairwise amino acid interactions will aid in understanding structural aspects of proteins.

Equilibrium and kinetic studies on folding of canine milk lysozyme

Thermal and chemical unfolding studies of the calcium-binding canine lysozyme (CL) by fluorescence and circular dichroism spectroscopy show that, upon unfolding in the absence of calcium ions, a very stable equilibrium intermediate state is formed. At room temperature and pH 7.5, for example, a stable molten globule state is attained in 3 M GdnHCl. The existence of such a pure and stable intermediate state allowed us to extend classical stopped-flow fluorescence measurements that describe the transition from the native to the unfolded form, with kinetic experiments that monitor separately the transition from the unfolded to the intermediate state and from the intermediate to the native state, respectively. The overall refolding kinetics of apo-canine lysozyme are characterized by a significant drop in the fluorescence intensity during the dead time, followed by a monoexponential increase of the fluorescence with k = 3.6 s(-1). Furthermore, the results show that, unlike its drastic effect on the stability, Ca(2+)-binding only marginally affects the refolding kinetics. During the refolding process of apo-CL non-native interactions, comparable to those observed in hen egg white lysozyme, are revealed by a substantial quenching of tryptophan fluorescence. The dissection of the refolding process in two distinct steps shows that these non-native interactions only occur in the final stage of the refolding process in which the two domains match to form the native conformation.

The osmophobic effect: natural selection of a thermodynamic force in protein folding

Intracellular organic osmolytes are present in certain organisms adapted to harsh environments and these osmolytes protect intracellular macromolecules against the denaturing environmental stress. In natural selection of organic osmolytes as protein stabilizers, it appears that the osmolyte property selected for is the unfavorable interaction between the osmolyte and the peptide backbone, a solvophobic thermodynamic force that we call the osmophobic effect. Because the peptide backbone is highly exposed to osmolyte in the denatured state, the osmophobic effect preferentially raises the free energy of the denatured state, shifting the equilibrium in favor of the native state. By focusing the solvophobic force on the denatured state, the native state is left free to function relatively unfettered by the presence of osmolyte. The osmophobic effect is a newly uncovered thermodynamic force in nature that complements the well-recognized hydrophobic interactions, hydrogen bonding, electrostatic and dispersion forces that drive protein folding. In organisms whose survival depends on the intracellular presence of osmolytes that can counteract denaturing stresses, the osmophobic effect is as fundamental to protein folding as these well-recognized forces.

Hydrophobic interaction chromatography of proteins

In this article, an overview of hydrophobic interaction chromatography (HIC) of proteins is given. After a brief description of protein hydrophobicity and hydrophobic interactions, we present the different proposed theories for the retention mechanism of proteins in HIC. Additionally, the main parameters to consider for the optimization of fractionation processes by HIC and the stationary phases available were described. Selected examples of protein fractionation by HIC are also presented.

Quantitative comparison of the ability of hydropathy scales to recognize surface beta-strands in proteins

Methods based on the use of hydropathy scales have been used widely to ascertain the secondary structures of proteins. However, over 100 such scales have been reported in the literature, and which of these is the most successful in terms of the prediction rate of the correct structure is not clear. This article, therefore, reports a comprehensive analysis of the relative success of hydropathy scales to locate beta-strands on the surfaces of proteins. The technique we used is based on the technique proposed by Fraser and Parry, but it includes a modification that allows a higher rate of successful prediction and a lower rate of overprediction. We used as a basis for assessing the predictions a database of sequence-unique structures that we previously established. Proteins 2001;42:243-255.

Mechanism of protein folding

The high structural resolution of the main transition states for the formation of native structure for the six small proteins of which Phi-values for a large set of mutants have become available, barstar, barnase, chymotrypsin inhibitor 2, Arc repressor, the src SH3 domain, and a tetrameric p53 domain reveals that for the first 5 of these proteins: (1) Residues that belong to regular secondary structure have a significantly larger average fraction of native structural consolidation than residues in loops; (2) on the other hand, secondary and tertiary structures have built up to the same degree, or at least a high degree, but nonuniformly distributed over the molecule; (3) the most consolidated parts of each protein molecule in the transition state cluster together, and these clusters contain a significantly higher percentage of residues that belong to regular secondary structure than the rest of the molecule. These observations further reconcile the framework model with the nucleation-condensation mechanism for folding: The amazing speed of protein folding can be understood as caused by the catalytic effect of the formation of clusters of residues which have particularly high preferences for the early formation of regular secondary structure in the presence of significant amounts of tertiary structure interactions.

Reversible denaturation of oligomeric human chaperonin 10: denatured state depends on chemical denaturant

Chaperonins cpn60/cpn10 (GroEL/GroES in Escherichia coli) assist folding of nonnative polypeptides. Folding of the chaperonins themselves is distinct in that it entails assembly of a sevenfold symmetrical structure. We have characterized denaturation and renaturation of the recombinant human chaperonin 10 (cpn10), which forms a heptamer. Denaturation induced by chemical denaturants urea and guanidine hydrochloride (GuHCl) as well as by heat was monitored by tyrosine fluorescence, far-ultraviolet circular dichroism, and cross-linking; all denaturation reactions were reversible. GuHCl-induced denaturation was found to be cpn10 concentration dependent, in accord with a native heptamer to denatured monomer transition. In contrast, urea-induced denaturation was not cpn10 concentration dependent, suggesting that under these conditions cpn10 heptamers denature without dissociation. There were no indications of equilibrium intermediates, such as folded monomers, in either denaturant. The different cpn10 denatured states observed in high [GuHCl] and high [urea] were supported by cross-linking experiments. Thermal denaturation revealed that monomer and heptamer reactions display the same enthalpy change (per monomer), whereas the entropy-increase is significantly larger for the heptamer. A thermodynamic cycle for oligomeric cpn10, combining chemical denaturation with the dissociation constant in absence of denaturant, shows that dissociated monomers are only marginally stable (3 kJ/mol). The thermodynamics for co-chaperonin stability appears conserved; therefore, instability of the monomer could be necessary to specify the native heptameric structure.

Guidelines for membrane protein engineering derived from de novo designed model peptides

Notwithstanding great advances in the engineering and structural analysis of globular proteins, relatively limited success has been achieved with membrane proteins--due largely to their intrinsic high insolubility and the concomitant difficulty in obtaining crystals. Progress with de novo synthesis of model membrane-interactive peptides presents an opportunity to construct simpler peptides with definable structures, and permits one to approach an understanding of the properties of the membrane proteins themselves. In the present article, we review how our laboratory and others have used peptide approaches to assess the detailed interactions of peptides with membranes, and primary folding at membrane surfaces and in membranes. Structural studies of model peptides identified the existence of a "threshold hydrophobicity," which controls spontaneous peptide insertion into membranes. Related studies of the relative helicity of peptides in organic media such as n-butanol indicate that the helical propensity of individual residues--not simply their hydrophobicity--may dictate the conformations of peptides in membranes. The overall experimental results provide fundamental guidelines for membrane protein engineering.

Stability and stabilization of globular proteins in solution

Proteins are multifunctional: their amino acid sequences simultaneously determine folding, function and turnover. Correspondingly, evolution selected for compromises between rigidity (stability) and flexibility (folding/function/degradation), to the result that generally the free energy of stabilization of globular proteins in solution is the equivalent to only a few weak intermolecular interactions. Additional increments may come from extrinsic factors such as ligands or specific compatible solutes. Apart from the enthalpic effects, entropy may play a role by reducing the flexibility (cystine bridges, increased proline content), or by water release from residues buried upon folding and association. Additional quaternary interactions and closer packing are typical characteristics of proteins from thermophiles. In halophiles, protein stability and function are maintained by increased ion binding and glutamic acid content, both allowing the protein inventory to compete for water at high salt. Acidophiles and alkalophiles show neutral intracellular pH; proteins facing the outside extremes of pH possess anomalously high contents in ionizable amino acids. Global comparisons of the amino acid compositions and sequences of proteins from mesophiles and extremophiles did not result in general rules of protein stabilization, even after including complete genome sequences into the search. Obviously, proteins are individuals that optimize internal packing and external solvent interactions by very different mechanisms, each protein in its own way. Strategies deduced from specific ultrastable proteins allow stabilizing point mutations to be predicted.

Surface Heterogeneity of Passively Oxidized Silicon Carbide Particles: Hydrophobic-Hydrophilic Partition

The surface of silicon carbide (SiC) particles previously subjected to passive oxidation has been characterized using various techniques such as adsorption of positively and negatively charged surfactants in aqueous solution, immersion microcalorimetry in three probe liquids, and flow microcalorimetry in organic media. Adsorption data show that around one quarter of the surface appears negatively charged and thus hydrophilic, while three quarters appear uncharged and hydrophobic. This is attributed to dissociation of silanols groups. Immersion calorimetry in liquids having well-defined polar and nonpolar components of surface energy shows that the Lifshitz-van der Waals component of SiC is very high and that the acid and basic components are weak. The experimental results appear to be consistent with both computations of surface energy using Lifshitz theory and experimental data previously obtained with other minerals. The three indexes are discussed and it is argued that they represent different terms of the solid surface energy. Copyright 2000 Academic Press.

Molecular recognition and binding of thermal hysteresis proteins to ice

Molecular recognition and binding are two very important processes in virtually all biological and chemical processes. An extremely interesting system involving recognition and binding is that of thermal hysteresis proteins at the ice-water interface. These proteins are of great scientific interest because of their antifreeze activity. Certain fish, insects and plants living in cold weather regions are known to generate these proteins for survival. A detailed molecular understanding of how these proteins work could assist in developing synthetic analogs for use in industry. Although the shapes of these proteins vary from completely alpha-helical to globular, they perform the same function. It is the shapes of these proteins that control their recognition and binding to a specific face of ice. Thermal hysteresis proteins modify the morphology of the ice crystal, thereby depressing the freezing point. Currently there are three hypotheses proposed with respect to the antifreeze activity of thermal hysteresis proteins. From structure-function experiments, ice etching experiments, X-ray structures and computer modeling at the ice-vacuum interface, the first recognition and binding hypothesis was proposed and stated that a lattice match of the ice oxygens with hydrogen-bonding groups on the proteins was important. Additional mutagenesis experiments and computer simulations have lead to the second hypothesis, which asserted that the hydrophobic portion of the amphiphilic helix of the type I thermal hysteresis proteins accumulates at the ice-water interface. A third hypothesis, also based on mutagenesis experiments and computer simulations, suggests that the thermal hysteresis proteins accumulate in the ice-water interface and actually influence the specific ice plane to which the thermal hysteresis protein ultimately binds. The first two hypotheses emphasize the aspect of the protein 'binding or accumulating' to specific faces of ice, while the third suggests that the protein assists in the development of the binding site. Our modeling and analysis supports the third hypothesis, however, the first two cannot be completely ruled out at this time. The objective of this paper is to review the computational and experimental efforts during the past 20 years to elucidate the recognition and binding of thermal hysteresis proteins at the ice-vacuum and ice-water interface.

The maximal affinity of ligands

We explore the question of what are the best ligands for macromolecular targets. A survey of experimental data on a large number of the strongest-binding ligands indicates that the free energy of binding increases with the number of nonhydrogen atoms with an initial slope of approximately -1.5 kcal/mol (1 cal = 4.18 J) per atom. For ligands that contain more than 15 nonhydrogen atoms, the free energy of binding increases very little with relative molecular mass. This nonlinearity is largely ascribed to nonthermodynamic factors. An analysis of the dominant interactions suggests that van der Waals interactions and hydrophobic effects provide a reasonable basis for understanding binding affinities across the entire set of ligands. Interesting outliers that bind unusually strongly on a per atom basis include metal ions, covalently attached ligands, and a few well known complexes such as biotin-avidin.

Fluoroalcohols as structure modifiers in peptides and proteins: hexafluoroacetone hydrate stabilizes a helical conformation of melittin at low pH

The effect of hexafluoroacetone hydrate (HFA) on the structure of the honey bee venom peptide melittin has been investigated. In aqueous solution at low pH melittin is predominantly unstructured. Addition of HFA at pH approximately 2.0 induces a structural transition from the unstructured state to a predominantly helical conformation as suggested by intense diagnostic far UV CD bands. The structural transition is highly cooperative and complete at 3.6 M (50% v/v) HFA. A similar structural transition is also observed in 2,2,2 trifluoroethanol which is complete only at a cosolvent concentration of approximately 8 M. Temperature dependent CD experiments support a 'cold denaturation' of melittin at low concentrations of HFA, suggesting that selective solvation of peptide by HFA is mediated by hydrophobic interactions. NMR studies in 3.6 M HFA establish a well-defined helical structure of melittin at low pH, as suggested by the presence of strong NH/NHi+1 NOEs throughout the sequence, along with many medium range helical NOEs. Structure calculations using NOE-driven distance constraints reveal a well-ordered helical fold with a relatively flexible segment around residues T10-G11-T12. The helical structure of melittin obtained at 3.6 M HFA at low pH is similar to those determined in methanolic solution and perdeuterated dodecylphosphocholine micelles. HFA as a cosolvent facilitates helix formation even in the highly charged C-terminal segment.

Thermodynamics of Interaction between Some Cellulose Ethers and SDS by Titration Microcalorimetry

The interaction between certain nonionic cellulose ethers (ethyl hydroxyethyl cellulose and hydroxypropyl methyl cellulose) and sodium dodecyl sulphate (SDS) has been investigated using isothermal titration microcalorimetry at temperatures between 25-50 degrees C. The observed heat flow curves have been interpreted in terms of a plausible mechanism of the interaction of the substituent groups with SDS monomers and clusters. The data have been related to changes occuring in the system at the macro- and microscopic levels with the addition of surfactants and with temperature. The process consists predominantly of polymer-surfactant interactions initially and surfactant-surfactant interactions at the later stages. A phenomenological model of the cooperative interaction (adsorption) process has been derived, and earlier published equilibrium binding data have been used to recover binding constants and Gibbs energy changes for this process. The adsorption enthalpies and entropies have been recovered along with the heat capacity change. The enthalpic cost of confining the nonpolar regions of the polymers in surfactant clusters is high, but the entropy gain from release of hydration shell water molecules as well as increased freedom of movement of these nonpolar regions in the clusters gives the process a strong entropic driving force. The process is entropy-driven initially and converts to being both enthalpy and entropy-driven at high SDS concentrations. An enthalpy-entropy compensation behavior is seen. Strongly negative heat capacity changes have been obtained resulting from the transfer of nonpolar groups from aqueous into nonpolar environments, as well as a reduction of conformational domains that the chains can populate. Changes in these two components cause the heat capacity change to become less negative at the higher binding levels. The system can be classified as exhibiting nonclassical hydrophobic binding at the later stages of binding. Copyright 1999 Academic Press.

Hydrophobicity regained

A widespread practice is to use free energies of transfer between organic solvents and water (delta G0transfer to define hydrophobicity scales for the amino acid side chains. A comparison of four delta G0transfer scales reveals that the values for hydrogen-bonding side chains are highly dependent on the non-aqueous environment. This property of polar side chains violates the assumptions underlying the paradigm of equating delta G0transfer with hydrophobicity or even with a generic solvation energy that is directly relevant to protein stability and ligand binding energetics. This simple regaining of the original concept of hydrophobicity reveals a flaw in approaches that use delta G0transfer values to derive generic estimates of the energetics of the burial of polar groups, and allows the introduction of a "pure" hydrophobicity scale for the amino acid residues.

Dissecting the energetics of a protein-protein interaction: the binding of ovomucoid third domain to elastase

An understanding of the structural basis for protein-protein interactions, and molecular recognition in general, requires complete characterization of binding energetics. Not only does this include quantification of the changes that occur in all of the thermodynamic parameters upon binding, including the enthalpy, entropy and heat capacity, but a description of how these changes are modulated by environmental conditions, most notably pH. Here, we have investigated the binding of turkey ovomucoid third domain (OMTKY3), a potent serine protease inhibitor, to the serine protease porcine pancreatic elastase (PPE) using isothermal titration calorimetry and structure-based thermodynamic calculations. We find that near neutral pH the binding energetics are influenced by a shift in the pKa of an ionizable group, most likely histidine 57 in the protease active site. Consequently, the observed binding energetics are strongly dependent upon solution conditions. Through a global analysis, the intrinsic energetics of binding have been determined, as have those associated with the pKa shift. The protonation energetics suggest that the drop in pKa is largely due to desolvation of the histidine residue. The resulting deprotonation is necessary for the enzymatic function of elastase. Intrinsically, at 25 degrees C the binding of OMTKY3 to PPE is characterized by an almost negligible enthalpy change, a large positive entropy change, and a large negative heat capacity change. These parameters are consistent with a model of the OMTKY3-PPE complex, which shows a large and significantly apolar protein-protein interface. Thermodynamic calculations based upon changes that occur in polar and apolar solvent-accessible surface area are in very good agreement with the measured intrinsic binding energetics.

Titration calorimetry of anesthetic-protein interaction: negative enthalpy of binding and anesthetic potency

Anesthetic potency increases at lower temperatures. In contrast, the transfer enthalpy of volatile anesthetics from water to macromolecules is usually positive. The transfer decreases at lower temperature. It was proposed that a few selective proteins bind volatile anesthetics with negative delta H, and these proteins are involved in signal transduction. There has been no report on direct estimation of binding delta H of anesthetics to proteins. This study used isothermal titration calorimetry to analyze chloroform binding to bovine serum albumin. The calorimetrically measured delta H cal was -10.37 kJ.mol-1. Thus the negative delta H of anesthetic binding is not limited to signal transduction proteins. The binding was saturable following Fermi-Dirac statistics and is characterized by the Langmuir adsorption isotherms, which is interfacial. The high-affinity association constant, K, was 2150 +/- 132 M-1 (KD = 0.47 mM) with the maximum binding number, Bmax = 3.7 +/- 0.2. The low-affinity K was 189 +/- 3.8 M-1 (KD = 5.29 mM), with a Bmax of 13.2 +/- 0.3. Anesthetic potency is a function of the activity of anesthetic molecules, not the concentration. Because the sign of delta H determines the temperature dependence of distribution of anesthetic molecules, it is irrelevant to the temperature dependence of anesthetic potency.

Solvation: how to obtain microscopic energies from partitioning and solvation experiments

Oil-water partitioning, solubilities, and vapor pressure experiments on small-molecule compounds are often used as models to obtain energies for biomolecular modeling. For example, measured partition coefficients, K, are often inserted into the formula -RT in K to obtain quantities thought to represent microscopic contact interaction free energies. We review evidence here that this procedure does not always give microscopically meaningful free energies. Some partitioning processes, particularly involving polymeric solvents such as octanol or hexadecane, are governed not only by translational entropies and contact interactions, but also by free energies resulting from changes in the conformations of the polymer chains upon solute insertion. The Flory-Huggins theory is more suitable for these situations than is the classical approach. We discuss the physical bases for both approaches.

Evaluation of linked protonation effects in protein binding reactions using isothermal titration calorimetry

A theoretical development in the evaluation of proton linkage in protein binding reactions by isothermal titration calorimetry (ITC) is presented. For a system in which binding is linked to protonation of an ionizable group on a protein, we show that by performing experiments as a function of pH in buffers with varying ionization enthalpy, one can determine the pK(a)'s of the group responsible for the proton linkage in the free and the liganded states, the protonation enthalpy for this group in these states, as well as the intrinsic energetics for ligand binding (delta H(o), delta S(o), and delta C(p)). Determination of intrinsic energetics in this fashion allows for comparison with energetics calculated empirically from structural information. It is shown that in addition to variation of the ligand binding constant with pH, the observed binding enthalpy and heat capacity change can undergo extreme deviations from their intrinsic values, depending upon pH and buffer conditions.

Energetics of hydrogen bonding in proteins: a model compound study

Differences in the energetics of amide-amide and amide-hydroxyl hydrogen bonds in proteins have been explored from the effect of hydroxyl groups on the structure and dissolution energetics of a series of crystalline cyclic dipeptides. The calorimetrically determined energetics are interpreted in light of the crystal structures of the studied compounds. Our results indicate that the amide-amide and amide-hydroxyl hydrogen bonds both provide considerable enthalpic stability, but that the amide-amide hydrogen bond is about twice that of the amide-hydroxyl. Additionally, the interaction of the hydroxyl group with water is seen most readily in its contributions to entropy and heat capacity changes. Surprisingly, the hydroxyl group shows weakly hydrophobic behavior in terms of these contributions. These results can be used to understand the effects of mutations on the stability of globular proteins.

Characterization of the sources of protein-ligand affinity: 1-sulfonato-8-(1')anilinonaphthalene binding to intestinal fatty acid binding protein

1-Sulfonato-8-(1')anilinonaphthalene (1,8-ANS) was employed as a fluorescent probe of the fatty acid binding site of recombinant rat intestinal fatty acid binding protein (1-FABP). The enhancement of fluorescence upon binding allowed direct determination of binding affinity by fluorescence titration experiments, and measurement of the effects on that affinity of temperature, pH, and ionic strength. Solvent isotope effects were also determined. These data were compared to results from isothermal titration calorimetry. We obtained values for the enthalpy and entropy of this interaction at a variety of temperatures, and hence determined the change in heat capacity of the system consequent upon binding. The ANS-1-FABP is enthalpically driven; above approximately 14 degrees C it is entropically opposed, but below this temperature the entropy makes a positive contribution to the binding. The changes we observe in both enthalpy and entropy of binding with temperature can be derived from the change in heat capacity upon binding by integration, which demonstrates the internal consistency of our results. Bound ANS is displaced by fatty acids and can itself displace fatty acids bound to I-FABP. The binding site for ANS appears to be inside the solvent-containing cavity observed in the x-ray crystal structure, the same cavity occupied by fatty acid. From the fluorescence spectrum and from an inversion of the Debye-Hueckel formula for the activity coefficients as a function of added salt, we inferred that this cavity is fairly polar in character, which is in keeping with inferences drawn from the x-ray structure. The binding affinity of ANS is considered to be a consequence of both electrostatic and conditional hydrophobic effects. We speculate that the observed change in heat capacity is produced mainly by the displacement of strongly hydrogen-bonded waters from the protein cavity.

Influence of salts on protein interactions at interfaces of amphiphilic polymers and adsorbents

The protein-binding capacity of two different amphiphilic adsorbents was investigated to determine the effect of solvent additives on the binding of proteins in hydrophobic-interaction chromatography. There was no simple correlation between binding capacity and the lyotropic series such as those suggested by the two different theories proposed by Arakawa and Narhi and Melander and Horváth. Proteins are known to be dynamic flexible objects which continuously undergo changes in conformation and which may well be influenced by chaotropic salts. Are conformational changes of proteins at interfaces an important parameter involved in protein interactions with amphiphilic polymers and adsorbents? In an attempt to answer this question, the reactivity of the thiol group in human serum albumin (HSA) toward N-ethyl-3-(2-pyridyldisulfanyl)propionamide dextran was used as a model system to evaluate its correlation with the lyotropic series. The results indicate that the thiol-disulfide exchange reaction at interfaces of amphiphilic polymers is influenced by the type of salt used.