Thermal stabilities of globular proteins

Statistical thermodynamic theory has recently been developed to account for the stabilities of globular proteins. Here we extend that work to predict the dependence on temperature. Folding is assumed to be driven by solvophobic interactions and opposed by the conformational entropy. The temperature dependence of the solvophobic interaction is taken from the transfer experiments on amino acids by Tanford and Nozaki and on model solutes by Gill and Wads?. One long-standing puzzle has been why proteins denature upon heating, since the solvophobic force to fold strengthens with increasing temperature. This is resolved by the theory, which predicts two first-order phase transitions. "Cold denaturation" is driven principally by the weakening of the solvophobic interaction, but normal denaturation is driven principally by the gain of conformational entropy of the chain. Predictions of the thermodynamic state functions are in reasonable agreement with the calorimetric experiments of Privalov and Khechinashvili. Comparison of the theory with experiments suggests that there may be an additional enthalpic driving force toward folding which is not due to the solvophobic interactions.     

Mutational analysis of a CBM family 5 chitin-binding domain of an alkaline chitinase from Bacillus sp. J813

Chitinase J from alkaliphilic Bacillus sp. J813 comprises a glycoside hydrolase (GH) family 18 catalytic domain (CatD), a fibronectin type III like domain, and a carbohydrate-binding module (CBM) family 5 chitin-binding domain (ChBD). It has been suggested that the ChBD binds to insoluble chitin and enhances its degradation by the CatD. To investigate the roles of two aromatic residues (Trp541 and Trp542), which are exposed on the surface of the ChBD, mutational analysis was performed. Single and double mutations of the two aromatic residues decreased binding and hydrolyzing abilities toward insoluble chitin. This result suggests that the ChBD binds to chitin by hydrophobic interactions via two surface-exposed aromatic residues. However, the double mutant, which has no such aromatic residue, bound to chitin at pH 5.2, probably by electrostatic interactions. Moreover, the ChBD bound to insoluble chitosan by electrostatic interactions.     

Contribution of enthalpy to the energetics of complex formation of aromatic ligands with DNA

The energy contributions of electrostatic, van der Waals interactions, hydrogen bonds, and interactions of charge transfer type to the enthalpy of complex formation of the double-stand DNA with the antitumor antibiotics daunomycin, nogalamycin, and novantron, as well as the mutagens ethidium bromide and proflavine have been calculated. According to the calculations, the van der Waals component (except for nogalamycin) is energetically favorable during complex formation of the antibiotics with DNA, and the contributions of H bonds and electrostatic interactions are unfavorable, with the probability of charge transfer in the complexes being low. It has been shown that the relatively low value of the experimental enthalpy of binding is the sum of components greater in absolute value and different in the sign, which is the cause of large errors in estimating the total enthalpy of complex formation of aromatic ligands with DNA.     

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.     

Mutational Analysis of a CBM Family 5 Chitin-Binding Domain of an Alkaline Chitinase from Bacillus sp. J813

Chitinase J from alkaliphilic Bacillus sp. J813 comprises a glycoside hydrolase (GH) family 18 catalytic domain (CatD), a fibronectin type III like domain, and a carbohydrate-binding module (CBM) family 5 chitin-binding domain (ChBD). It has been suggested that the ChBD binds to insoluble chitin and enhances its degradation by the CatD. To investigate the roles of two aromatic residues (Trp541 and Trp542), which are exposed on the surface of the ChBD, mutational analysis was performed. Single and double mutations of the two aromatic residues decreased binding and hydrolyzing abilities toward insoluble chitin. This result suggests that the ChBD binds to chitin by hydrophobic interactions via two surface-exposed aromatic residues. However, the double mutant, which has no such aromatic residue, bound to chitin at pH 5.2, probably by electrostatic interactions. Moreover, the ChBD bound to insoluble chitosan by electrostatic interactions.     

Halogen bond tunability I: the effects of aromatic fluorine substitution on the strengths of halogen-bonding interactions involving chlorine,bromine, and iodine

In the past several years, halogen bonds have been shown to be relevant in crystal engineering and biomedical applications. One of the reasons for the utility of these types of noncovalent interactions in the development of, for example, pharmaceutical ligands is that their strengths and geometric properties are very tunable. That is, substitution of atoms or chemical groups in the vicinity of a halogen can have a very strong effect on the strength of the halogen bond. In this study we investigate halogen-bonding interactions involving aromatically-bound halogens (Cl, Br, and I) and a carbonyl oxygen. The properties of these halogen bonds are modulated by substitution of aromatic hydrogens with fluorines, which are very electronegative. It is found that these types of substitutions have dramatic effects on the strengths of the halogen bonds, leading to interactions that can be up to 100% stronger. Very good correlations are obtained between the interaction energies and the magnitudes of the positive electrostatic potentials (σ-holes) on the halogens. Interestingly, it is seen that the substitution of fluorines in systems containing smaller halogens results in electrostatic potentials resembling those of systems with larger halogens, with correspondingly stronger interaction energies. It is also shown that aromatic fluorine substitutions affect the optimal geometries of the halogen-bonded complexes, often as the result of secondary interactions.     

Interactions of aromatic residues of proteins with nucleic acids. Fluorescence studies of the binding of oligopeptides containing tryptophan and tyrosine residues to polynucleotides.

The binding of oligopeptides of general structure Lys-X-Lys (where X is an aromatic residue) to several polynucleotides has been studied by fluorescence spectroscopy. Two types of complexes are formed, both involving electrostatic interactions between lysyl residues and phosphate groups as shown by the ionic strength and pH dependence of binding. The fluorescence quantum yield of the first complex is identical with that of the free peptide. The other complex involves a stacking of the nucleic acid bases with the aromatic amino acid whose fluorescence is quenched. Fluorescence data have been quantitatively analyzed according to a model involving these two types of complexes. Association constants and the size of binding sites have been determined. Stacking interactions are favored in single-stranded polynucleotides as compared to double-stranded ones. A short oligopeptide such as Lys-X-Lys is thus able to distinguish between single-stranded and double-stranded nucleic acids. Fluorescence results are compared to those obtained by proton magnetic resonance and circular dichroism.     

Calculation and display of electrostatic potentials

A general methodology is developed for incorporating accurate electrostatic information from ab initio molecular orbital calculations into molecular mechanics calculations. Examples are given of the method applied to simple aromatic organic molecules. A program has been developed for displaying the results of the ab initio calculations on a Silicon Graphics workstation. The technique developed here provides an alternative method for including electrostatic interactions in molecular mechanics calculations and is compared with other methods for determining atomic charges.     

Structural studies of metalloporphyrins. Part XIa: Complexes of water-soluble zinc(II) porphyrins with amino acids: influence of ligand-ligand interactions on the stability of the complexes

The stability constants of a series of complexes of the cationic water-soluble porphyrin ZnTMPyP with various amino acids have been determined by 1H NMR spectroscopy at pH 10.5. The following stability order has been observed: Tyr greater than Phe, Glu greater than Asp greater than Ile greater than Val greater than Gly. These results can be best rationalized by invoking complex stabilization due to ligand-ligand (e.g., stacking or electrostatic) interactions. Evidence for stacking interactions between the porphyrin ring and the aromatic ring of phenylalanine, tyrosine, and tryptophan was further provided by study of the complexation of these amino acids with the free-base porphyrin TMPyPH2. In this case, complexation constants increased in the order: Phe less than Tyr less than Trp. Attempts to form complexes of the amino acids with the anionic porphyrin ZnTCPP proved unsuccessful, indicating that electrostatic interactions play a major role in the stability of the zinc porphyrin-amino acids complexes.     

Halogen bond tunability II: the varying roles of electrostatic and dispersion contributions to attraction in halogen bonds

In a previous study we investigated the effects of aromatic fluorine substitution on the strengths of the halogen bonds in halobenzene…acetone complexes (halo?=?chloro, bromo, and iodo). In this work, we have examined the origins of these halogen bonds (excluding the iodo systems), more specifically, the relative contributions of electrostatic and dispersion forces in these interactions and how these contributions change when halogen σ-holes are modified. These studies have been carried out using density functional symmetry adapted perturbation theory (DFT-SAPT) and through analyses of intermolecular correlation energies and molecular electrostatic potentials. It is found that electrostatic and dispersion contributions to attraction in halogen bonds vary from complex to complex, but are generally quite similar in magnitude. Not surprisingly, increasing the size and positive nature of a halogen’s σ-hole dramatically enhances the strength of the electrostatic component of the halogen bonding interaction. Not so obviously, halogens with larger, more positive σ-holes tend to exhibit weaker dispersion interactions, which is attributable to the lower local polarizabilities of the larger σ-holes.

Probing the specificity of the S1 binding site of subtilisin Carlsberg with boronic acids

The binding properties and limitations of the key S1 site of subtilisin Carlsberg have been probed with boronic acid inhibitors bearing structurally varied substituents ranging from small alkyl to large aromatic groups. The data permit structural features favoring, and disfavoring, good S1 binding to be clarified. In addition, applications of electrostatic energy calculations have identified a hitherto unsuspected region of positive potential in the fundamentally hydrophobic S1 pocket, whose interactions with electronegative substituents of inhibitors can make significant binding contributions.     

Energetics of complex formation of the DNA hairpin structure d(GCGAAGC) with aromatic ligands

The energy contributions of various physical interactions to the total Gibbs energy of complex formation of the biologically important DNA hairpin d(GCGAAGC) with aromatic antitumor antibiotics daunomycin and novantron and the mutagens ethidium and proflavine have been calculated. It has been shown that the relatively small value of the total energy of binding of the ligands to the hairpin is the sum of components great in absolute value and different in sign. The contributions of van der Waals interactions and both intra- and intermolecular hydrogen bonds and bonds with aqueous environment have been studied. According to the calculations, the hydrophobic and van der Waals components are energetically favorable in complex formation of the ligands with the DNA pairpin d(GCGAAGC), whereas the electrostatic (with consideration of hydrogen bonds) and entropic components are unfavorable.     

Energetics of complex formation of the dna hairpin structure d(GCGAAGC) with aromatic ligands

The energy contributions of various physical interactions to the total Gibbs energy of complex formation of the biologically important DNA hairpin d(GCGAAGC) with aromatic antitumor antibiotics daunomycin and novantrone and mutagens ethidium bromide and proflavin have been calculated. It has been shown that the relatively small value of the total energy of binding of the ligands to the hairpin is the sum of components great in absolute value and different in sign. The contributions of van der Waals interactions and both intra- and intermolecular hydrogen bonds and bonds with aqueous environment have been studied. According to the calculations, the hydrophobic and van der Waals components are energetically favorable in complex formation of the ligands with the hairpin d(GCGAAGC), whereas the electrostatic (with consideration of hydrogen bonds) and entropic components are unfavorable.     

Electrostatic contribution to the energy of binding of aromatic ligands with DNA

The method of solution of the nonlinear Poisson-Boltzmann equation was used to calculate electrostatic energy of binding of various aromatic ligands with DNA oligomers of different length. Analysis of the electrostatic contribution was made in terms of a two-step DNA binding process: formation of the intercalation cavity and insertion of the ligand. The total electrostatic energy was also partitioned into components: the energy of atom-atom coulombic interactions and the energy of interaction with surrounding water. The results indicate that electrostatic interactions are, as a whole, unfavorable to the intercalation process and that a correct analysis of structure-energy interrelation for Ligand-DNA interactions should only be accomplished at the level of the components rather than at the level of total electrostatic energy.     

Hydrophobic, aromatic, and electrostatic interactions play a central role in amyloid fibril formation and stability

Amyloid-like fibrous crystals formed by the peptide KFFEAAAKKFFE have been previously characterized and provide an ideal model system to examine the importance of specific interactions by introducing specific substitutions. We find that the removal of any phenylalanine residue completely abrogates assembly ability, while charged residues modulate interactions within the structure resulting in alternative fibrillar morphologies. X-ray fiber diffraction analysis reveals that the essential backbone packing of the peptide molecules is maintained, while small changes accommodate differences in side chain size in the variants. We conclude that even very short peptides are adaptable and add to the growing knowledge regarding amyloid polymorphisms. Additionally, this work impacts on our understanding of the importance of residue composition for amyloidogenic peptides, in particular the roles of electrostatic, aromatic, and hydrophobic interactions in amyloid assembly.     

Lecithin reversed micelle formation and interactions with p,p′-DDA and kelthane. NMR observations

The behavior of the reversed micellar system of lecithin in benzene in the presence of water has been investigated using lanthanide NMR shift reagents and two oxygenated chlorinated hydrocarbon pesticide systems as both reporter molecules and amphiphilic ligands. Increased solubilization of these hydrocarbons in benzene in the presence of lecithin was observed with concomitant deshielding of aromatic proton signals. The experiments indicate that the deshielding may not be due to interactions with lecithin bound phosphate groups but may be forced electrostatic interactions of an intermolecular type. Negligible solubility of these hydrocarbons in benzene was observed in the presence of lecithin in which the anions associated with choline have been partially replaced with the 2,2-bis(4-chlorophenyl) acetate ion.     

阳离子-π相互作用在两种典型蛋白质折叠型中的偏好性

阳离子-π相互作用是一种在阳离子和芳香性体系之间形成的一种作用力。在蛋白质中,带正电荷的氨基酸(Lys、Arg)和芳香族氨基酸(Phe、Tyr、Trp)之间可以形成阳离子-π相互作用。对α/β类蛋白中两种典型折叠类型――单绕和双绕的研究表明:(1)单绕结构中阳离子-π相互作用的分布密度是双绕结构的2.3倍。(2)Arg-Phe组合偏好在双绕中出现,Arg-Tyr组合偏好在单绕中出现。(3)在单绕中除Lys-Phe组合外,其余5种组合的阳离子-π相互作用能量高于双绕的对应组合,其中以Arg-Trp组合的能量最高。(4)在单绕结构中,样本所含氨基酸残基数量和样本中阳离子-π的数量有明显的相关性,在双绕结构中没有发现类似的相关性。(5)在单绕和双绕结构当中,把阳离子-π相互作用能量分解为静电能和范德华力能揭示出静电能与范德华力能之比接近2∶1,静电作用在阳离子-π相互作用中起主要作用。     

Solvent denaturation and stabilization of globular proteins

Statistical thermodynamic theory has recently been developed to account for the stabilities of globular proteins. Here we extend that work to predict the effects of solvents on protein stability. Folding is assumed to be driven by solvophobic interactions and opposed by conformational entropy. The solvent dependence of the solvophobic interactions is taken from transfer experiments of Nozaki and Tanford on amino acids into aqueous solutions of urea or guanidine hydrochloride (GuHCl). On the basis of the assumption of two pathways involving collapse and formation of a core, the theory predicts that increasing denaturant should lead to a two-state denaturation transition (i.e., there is a stable state along each path separated by a free energy barrier). The denaturation midpoint is predicted to occur at higher concentrations of urea than of GuHCl. At neutral pH, the radius of the solvent-denatured state should be much smaller than for a random-flight chain and increase with either denaturant concentration or number of polar residues in the chain. A question of interest is whether free energies of folding should depend linearly on denaturant, as is often assumed. The free energy is predicted to be linear for urea but to have some small curvature for GuHCl. Predicted slopes and exposed areas of the unfolded states are found to be in generally good agreement with experiments. We also discuss stabilizing solvents and compare thermal with solvent denaturation.     

Influence of cation-pi interactions in different folding types of membrane proteins

Cation-pi interactions play an important role to the stability of protein structures. In this work, we analyze the influence of cation-pi interactions in three-dimensional structures of membrane proteins. We found that transmembrane strand (TMS) proteins have more number of cation-pi interactions than transmembrane helical (TMH) proteins. In TMH proteins, both the positively charged residues Lys and Arg equally experience favorable cation-pi interactions whereas in TMS proteins, Arg is more likely than Lys to be in such interactions. There is no relationship between number of cation-pi interactions and number of residues in TMH proteins whereas a good correlation was observed in TMS proteins. The average cation-pi interaction energy for TMH proteins is -16 kcal/mol and that for TMS proteins is -27 kcal/mol. The pair-wise cation-pi interaction energy between aromatic and positively charged residues showed that Lys-Trp energy is stronger in TMS proteins than TMH proteins; Arg-Phe, Arg-Tyr and Lys-Phe have higher energy in TMH proteins than TMS proteins. The decomposition of energies into electrostatic and van der Waals revealed that the contribution from electrostatic energy is twice as that from van der Waals energy in both TMH and TMS proteins. The results obtained in the present study would be helpful to understand the contribution of cation-pi interactions to the stability of membrane proteins.     

Aromatic-aromatic interactions in the zinc finger motif. Analysis of the two-dimensional nuclear magnetic resonance structure of a mutant domain.

The folding and stability of globular proteins are determined by a variety of chemical mechanisms, including hydrogen bonds, salt bridges and the hydrophobic effect. Of particular interest are weakly polar interactions involving aromatic rings, which are proposed to regulate the geometry of closely packed protein interiors. Such interactions reflect the electrostatic contribution of pi-electrons and, unlike van der Waals' interactions and the hydrophobic effect, may, in principle, introduce a directional force in a protein's hydrophobic core. Although the weakly polar hypothesis is supported by a statistical analysis of protein structures, the general importance of such contributions to protein folding and stability is unclear. Here, we show the presence of alternative aromatic-aromatic interactions in the two-dimensional nuclear magnetic resonance structure of a mutant Zn finger. Changes in aromatic packing lead in turn to local and non-local differences between the structures of a wild-type and mutant domain. The results provide insight into the evolution of Zn finger sequences and have implications for understanding how geometric relationships may be chemically encoded in a simple sequence template.