A general rule for the relationship between hydrophobic effect and conformational stability of a protein: stability and structure of a series of hydrophobic mutants of human lysozyme

Previous research has shown that head direction (HD) cells in both the anterior dorsal thalamus (ADN) and the postsubiculum (PoS) in rats discharge in relation to familiar, visual landmarks in the environment. This study assessed whether PoS and ADN HD cells would be similarly responsive to nonvisual or unfamiliar environmental cues. After visual input was eliminated by blindfolding the rats, HD cells maintained direction-specific discharge, but their preferred firing directions became less stable. In addition, rotations of the behavioral apparatus indicated that some nonvisual cues (presumably tactile, olfactory, or both) exerted above chance stimulus control over a cell's preferred firing direction. However, a prominent auditory cue was not effective in exerting stimulus control over a cell's preferred direction. HD cell activity also was assessed after rotation of a novel visual cue exposed to the rat for 1, 3, or 8 min. An 8-min exposure was enough time for a novel visual cue to gain control over a cell's preferred direction, whereas an exposure of 1 or 3 min led to control in only about half the sessions. These latter results indicate that HD cells rely on a rapid learning mechanism to develop associations with landmark cues.     

Stability of hydrophobic nuclei of protein molecules depending of the degree of packing

The values of free energies for the formation of hydrophobic cores with different degree of packing of nonpolar side chains have been calculated. It has been shown that the difference between the values of free energy for tightly packed cores inaccessible to water solvent molecules and the values of free energies for less tightly packed hydrophobic cores whose surface is in contact with the water solvent, is 50-60 kJ/mol and coincides with the value of the enzyme free energy change during its interaction with the substrate.     

Extremely fast folding of a very stable leucine zipper with a strengthened hydrophobic core and lacking electrostatic interactions between helices

The dimer interface of a leucine zipper involves hydrophobic as well as electrostatic interactions between the component helices. Here we ask how hydrophobic effects and electrostatic repulsion balance the rate of folding and thermodynamic stability of a designed dimeric leucine zipper formed by the acidic peptide A that contains four repeating sequence units, (abcdefg)4. The aliphatic a and d residues of peptide A were the same as in the GCN4 leucine zipper but the e and g positions were occupied by Glu, which prevented folding above pH 6 because of electrostatic repulsion. Leucine zipper A2 was formed by protonation of the e and g side chains with a sharp transition midpoint at pH 5.2. Folding could be described by a two-state transition from two unfolded random coil monomers to a coiled coil dimer. There was a linear relationship between the logarithm of the rate constants and the number of repulsive charges on the folded leucine zipper dimer. The same linear relationship applied to the free energy of unfolding and the number of repulsive charges at thermodynamic equilibrium. Fully protonated peptide A folded at a near diffusion-limited rate (kon = 3 x 10(8) M-1 s-1), and the free energy of folding was -55 kJ mol-1 at 25 degrees C. The present work shows that protonation of Glu in positions e and g increases both the folding rate and the stability of the leucine zipper in the absence of any interhelical electrostatic interactions. Protonated Glu is proposed to act like a nonpolar residue and to strengthen the hydrophobic core by folding back toward the core residues in the a and d positions. This effect adds more to the free energy of unfolding and to the rate of folding than maximizing the number of salt bridges across the helix interface in an electrostatically stabilized heterodimeric leucine zipper [Wendt, H., Leder, L., H?rm?, H., Jelesarov, I., Baici, A., and Bosshard, H. R. (1997) Biochemistry 36, 204-213].     

Theory of hydrophobic interactions

An earlier theoretical approach to the hydrophobic interaction, based on statistical mechanical treatments of models for liquid water and for aqueous solutions of hydrocarbons, is summarized here. Experimental verification of the theoretical thermodynamic parameters for hydrophobic interactions, and applications of the theory to some aspects of protein structure, are presented.     

In vivo analyses of interactions between SecE and SecY, core components of the Escherichia coli protein translocation machinery

We have carried out structure-function studies on the cytoplasmic membrane protein, SecE, a component of the Escherichia coli secretion machinery. SecE, along with SecY, form a complex in the cytoplasmic membrane essential for protein translocation. By directed mutagenesis, we altered highly conserved residues of the second cytoplasmic domain (CD2) and of the COOH-terminal periplasmic region (PD2) of SecE. These mutants, as well as previously constructed mutations in the third membrane-spanning segment of SecE (MSS3), were tested for their ability to complement a secE null mutation, for their effects on protein export in vivo, and for their ability to form a stable complex with SecY. Most single mutations at the conserved positions in CD2 caused secretion defects, but had little effect on growth at 37 degrees C. Double mutations in CD2, or the introduction or removal of proline residues, affected growth and protein translocation more severely. Co-immunoprecipitations of SecE and SecY revealed that all mutant proteins, except those altered in PD2, destabilized the SecE-SecY complex. These results suggest that several regions contribute to the formation of a stable SecE-SecY complex but the elimination of a single contact point does not necessarily affect the functionality of the complex.     

EPR mapping of interactions between spin-labeled variants of human carbonic anhydrase II and GroEL: evidence for increased flexibility of the hydrophobic core by the interaction

Human carbonic anhydrase II (HCA II) interacts weakly with GroEL at room temperature. To further investigate this interaction we used electron paramagnetic resonance (EPR) spectroscopy to study HCA II cysteine mutants spin-labeled at selected positions. From our results it is evident that protein-protein interactions can be specifically mapped by site-directed spin-labeling and EPR measurements. HCA II needs to be unfolded to about the same extent as a GuHCl-induced molten-globule intermediate of the enzyme to interact with GroEL. The interaction with GroEL includes interactions with outer parts of the HCA II molecule, such as peripheral beta-strands and the N-terminal domain, which have previously been shown to be rather unstable. As a result of the interaction, the rigid and compact hydrophobic core exhibits higher flexibility than in the molten globule, which is likely to facilitate rearrangements of misfolded structure during the folding process. The degree of binding to GroEL and accompanying inactivation of the enzyme depend on the stability of the HCA II variant, and nonspecific hydrophobic interactions appear to be most important in stabilizing the GroEL-substrate complex.     

Native protein fluctuations: the conformational-motion temperature and the inverse correlation of protein flexibility with protein stability

We study the fluctuations of native proteins by exact enumeration using the HP lattice model. The model fluctuations increase with temperature. We observe a low-temperature point, below which large fluctuations are frozen out. This prediction is consistent with the observation by Tilton et al. [R. F. Tilton, Jr., J. C. Dewan, and G. A. Petsko, Biochemistry 31, 2469 (1992)], that the thermal motions of ribonuclease A increase sharply above about 200 K. We also explore protein "flexibility" as defined by Debye-Waller-like factors and solvent accessibilities of core residues to hydrogen exchange. We find that proteins having greater stability tend to have fewer large fluctuations, and hence lower flexibilities. If flexibility is necessary for enzyme catalysis, this could explain why proteins from thermophilic organisms, which are exceptionally stable, may be catalytically inactive at normal temperatures.     

Temperature, stability, and the hydrophobic interaction

Changes in free energy are normally used to track the effect of temperature on the stability of proteins and hydrophobic interactions. Use of this procedure on the aqueous solubility of hydrocarbons, a standard representation of the hydrophobic effect, leads to the conclusion that the hydrophobic effect increases in strength as the temperature is raised to approximately 140 degrees C. Acceptance of this interpretation leads to a number of far-reaching conclusions that are at variance with the original conception of the hydrophobic effect and add considerably to the complexity of interpretation. There are two legitimate thermodynamic functions that can be used to look at stability as a function of temperature: the standard Gibbs free energy change, deltaG degrees, and deltaG degrees/T. The latter is proportional to the log of the equilibrium constant and is sometimes called the Massieu-Planck function. Arguments are presented for using deltaG degrees/T rather than deltaG degrees for variations in stability with temperature. This makes a considerable difference in the interpretation of the hydrophobic interaction, but makes little change in the stability profile of proteins. Protein unfolding and the aqueous solubility of benzene are given as examples. The contrast between protein unfolding and the hydration of nonpolar molecules provides a rough estimate of the contribution of other factors that stabilize and destabilize protein structure.     

De novo design of the hydrophobic core of ubiquitin

We have previously reported the development and evaluation of a computational program to assist in the design of hydrophobic cores of proteins. In an effort to investigate the role of core packing in protein structure, we have used this program, referred to as Repacking of Cores (ROC), to design several variants of the protein ubiquitin. Nine ubiquitin variants containing from three to eight hydrophobic core mutations were constructed, purified, and characterized in terms of their stability and their ability to adopt a uniquely folded native-like conformation. In general, designed ubiquitin variants are more stable than control variants in which the hydrophobic core was chosen randomly. However, in contrast to previous results with 434 cro, all designs are destabilized relative to the wild-type (WT) protein. This raises the possibility that beta-sheet structures have more stringent packing requirements than alpha-helical proteins. A more striking observation is that all variants, including random controls, adopt fairly well-defined conformations, regardless of their stability. This result supports conclusions from the cro studies that non-core residues contribute significantly to the conformational uniqueness of these proteins while core packing largely affects protein stability and has less impact on the nature or uniqueness of the fold. Concurrent with the above work, we used stability data on the nine ubiquitin variants to evaluate and improve the predictive ability of our core packing algorithm. Additional versions of the program were generated that differ in potential function parameters and sampling of side chain conformers. Reasonable correlations between experimental and predicted stabilities suggest the program will be useful in future studies to design variants with stabilities closer to that of the native protein. Taken together, the present study provides further clarification of the role of specific packing interactions in protein structure and stability, and demonstrates the benefit of using systematic computational methods to predict core packing arrangements for the design of proteins.     

Free energy of burying hydrophobic residues in the interface between protein subunits

We have obtained an experimental estimate of the free energy change associated with variations at the interface between protein subunits, a subject that has raised considerable interest since the concept of accessible surface area was introduced by Lee and Richards [Lee, B. & Richards, F. M. (1971) J. Mol. Biol. 55, 379-400]. We determined by analytical ultracentrifugation the dimer-tetramer equilibrium constant of five single and three double mutants of human Hb. One mutation is at the stationary alpha1 beta1 interface, and all of the others are at the sliding alpha1 beta2 interface where cleavage of the tetramer into dimers and ligand-linked allosteric changes are known to occur. A surprisingly good linear correlation between the change in the free energy of association of the mutants and the change in buried hydrophobic surface area was obtained, after corrections for the energetic cost of losing steric complementarity at the alphabeta dimer interface. The slope yields an interface stabilization free energy of -15 +/- 1.2 cal/mol upon burial of 1 A2 of hydrophobic surface, in very good agreement with the theoretical estimate given by Eisenberg and McLachlan [Eisenberg, D. & McLachlan, A. D. (1986) Nature (London) 319, 199-203].     

Protein stabilization by hydrophobic interactions at the surface

The contribution of the solvent-exposed residue 63 to thermal stability of the thermolysin-like neutral protease of Bacillus stearothermophilus was studied by analyzing the effect of twelve different amino acid substitutions at this position. The thermal stability of the enzyme was increased considerably by introducing Arg, Lys or bulky hydrophobic amino acids. In general, the effects of the mutations showed that hydrophobic contacts in this surface-located region of the protein are a major determinant of thermal stability. This observation contrasts with general concepts concerning the contribution of surface-located residues and surface hydrophobicity to protein stability and indicates new ways for protein stabilization by site-directed mutagenesis.     

Simulating the minimum core for hydrophobic collapse in globular proteins

To investigate the nature of hydrophobic collapse considered to be the driving force in protein folding, we have simulated aqueous solutions of two model hydrophobic solutes, methane and isobutylene. Using a novel methodology for determining contacts, we can precisely follow hydrophobic aggregation as it proceeds through three stages: dispersed, transition, and collapsed. Theoretical modeling of the cluster formation observed by simulation indicates that this aggregation is cooperative and that the simulations favor the formation of a single cluster midway through the transition stage. This defines a minimum solute hydrophobic core volume. We compare this with protein hydrophobic core volumes determined from solved crystal structures. Our analysis shows that the solute core volume roughly estimates the minimum core size required for independent hydrophobic stabilization of a protein and defines a limiting concentration of nonpolar residues that can cause hydrophobic collapse. These results suggest that the physical forces driving aggregation of hydrophobic molecules in water is indeed responsible for protein folding.     

Molecular analysis of the interactions between protein kinase C-epsilon and filamentous actin

Protein kinase C-epsilon (PKC-epsilon) contains a putative actin binding motif that is unique to this individual member of the PKC gene family. We have used deletion mutagenesis to determine whether this hexapeptide motif is required for the physical association of PKC-epsilon and actin. Full-length recombinant PKC-epsilon, but not PKC-betaII, -delta, -eta, or -zeta, bound to filamentous actin in a phorbol ester-dependent manner. Deletion of PKC-epsilon amino acids 222-230, encompassing a putative actin binding motif, completely abrogated this binding activity. When NIH 3T3 cells overexpressing either PKC-epsilon or the deletion mutant of this isozyme were treated with phorbol ester only wild-type PKC-epsilon colocalized with actin in zones of cell adhesion. In binary reactions, it was possible to demonstrate that purified filamentous actin is capable of directly stimulating PKC-epsilon phosphotransferase activity. These and other findings support the hypothesis that a conformationally hidden actin binding motif in the PKC-epsilon sequence becomes exposed upon activation of this isozyme and functions as a dominant localization signal in NIH 3T3 fibroblasts. This protein-protein interaction is sufficient to maintain PKC-epsilon in a catalytically active conformation.     

The hydrophobic repeated domain of the Clostridium cellulovorans cellulose-binding protein (CbpA) has specific interactions with endoglucanases

We overexpressed one of the hydrophobic repeated domains (HBDs) (110 amino acid residues) of the cellulose-binding protein (CbpA) from Clostridium cellulovorans by making a hybrid protein with the Escherichia coli maltose-binding protein (MalE). The HBD was purified to homogeneity, and interactions between the HBD and endoglucanases were analyzed by a novel interaction Western blotting (immunoblotting) method. The HBD had specific interactions with endoglucanases (EngB and EngD) from C. cellulovorans. These results indicated that the HBD was an endoglucanase binding site of CbpA.     

A method for detecting hydrophobic patches on protein surfaces

A method for the detection of hydrophobic patches on the surfaces of protein tertiary structures is presented. It delineates explicit contiguous pieces of surface of arbitrary size and shape that consist solely of carbon and sulphur atoms using a dot representation of the solvent-accessible surface. The technique is also useful in detecting surface segments with other characteristics, such as polar patches. Its potential as a tool in the study of protein-protein interactions and substrate recognition is demonstrated by applying the method to myoglobin, Leu/IIe/Val-binding protein, lipase, lysozyme, azurin, triose phosphate isomerase, carbonic anhydrase, and phosphoglycerate kinase. Only the largest patches, having sizes exceeding random expectation, are deemed meaningful. In addition to well-known hydrophobic patches on these proteins, a number of other patches are found, and their significance is discussed. The method is simple, fast, and robust. The program text is obtainable by anonymous ftp.     

The snRNP core protein SmB and tissue-specific SmN protein are differentially distributed between snRNP particles

The SmN protein is a tissue specific component of the small nuclear ribonucleoprotein particle which is closely related to the ubiquitously expressed SmB protein but is expressed only in the brain and heart. To investigate the function of SmN, its localisation within different snRNP particles was investigated using a range of anti-snRNP monoclonal antibodies. SmN and SmB were found to exhibit different patterns of association with snRNP particles in two cell lines, ND7 and F9 which express SmN. In both cases, SmN was found to be present in the U-2 snRNP but was excluded from the U-1 snRNPs whereas SmB was present in both U-1 and U-2 snRNPs. Data from transfected 3T3 mouse fibroblasts cell lines artificially expressing a low level of SmN also confirm this observation. In contrast, SmN was found to be an integral component of both the U-1 and U-2 snRNPs in both 3T3 cells artificially expressing high levels of SmN and in adult rat brain which has a naturally high level of SmN expression. Taken together, the results suggest that the pre-U1 snRNP particle has a lower affinity for SmN than for SmB. Thus, SmN expressed at low levels incorporates into U2, but SmN expressed at high levels incorporates into both U1 and U2 snRNPs and replaces SmB. The significance of these effects is discussed in terms of the potential role played by SmN in constitutive and alternative splicing pathways in neuronal cells.     

Correlation between protein stability and crystal properties of designed ROP variants

Six variants of the ROP protein, designed with the aim to analyze by X-ray crystallography loop formation and core packing interactions in 4-alpha-helical bundles, have been purified and a search of their crystallization conditions has been carried out. Five mutants yield crystals that are suitable for medium to high resolution X-ray diffraction studies. For all mutants crystal size, sensitivity to X-irradiation and diffraction limit are correlated to their stability as determined by differential scanning calorimetry, in a manner which is not yet understood in detail.     

Characterization of interactions between the anti-apoptotic protein BAG-1 and Hsc70 molecular chaperones

The anti-cell death protein BAG-1 binds to 70-kDa heat shock proteins (Hsp70/Hsc70) and modulates their chaperone activity. Among other facilitory roles, BAG-1 may serve as a nucleotide exchange factor for Hsp70/Hsc70 family proteins and thus represents the first example of a eukaryotic homologue of the bacterial co-chaperone GrpE. In this study, the interactions between BAG-1 and Hsc70 are characterized and compared with the analogous GrpE-DnaK bacterial system. In contrast to GrpE, which binds DnaK as a dimer, BAG-1 binds to Hsc70 as a monomer with a 1:1 stoichiometry. Dynamic light scattering, sedimentation equilibrium, and circular dichroism measurements provided evidence that BAG-1 exists as an elongated, highly helical monomer in solution. Isothermal titration microcalorimetry was used to determine the complex stoichiometry and an equilibrium dissociation constant, KD, of 100 nM. Kinetic analysis using surface plasmon resonance yielded a KD consistent with the calorimetrically determined value. Molecular modeling permitted a comparison of structural features between the functionally homologous BAG-1 and GrpE proteins. These data were used to propose a mechanism for BAG-1 in the regulation of Hsp70/Hsc70 chaperone activity.     

Potentiating interactions between morphogenetic protein and neurotrophic factors in developing neurons

A highly specific antibody against sulfatide, a myelin-associated glycolipid, has been investigated using indirect double immunocytochemistry in rat primary astroglial cultures from cerebral cortex. Sulfatide was expressed in a selected subpopulation of astrocytes (2-3%) and was found to be exclusively located intracellularly. The sulfatide-positive cells appeared in two different morphologies: flat and stellate. Immunolabeling of the astroglial cultures showed that sulfatide always co-existed with GFAP or S-100, and in some cells with GD3 (flat 90% and stellate 50%) or A2B5 (1%) antibody. The sulfatide-positive cells did not bind the O1 antibody, which is used as a marker for oligodendrocytes. Glial cultures from other regions and mixed cultures, with both neurons and glial cells, were examined and showed similar results. Biochemical analysis by TLC-ELISA verified the presence of sulfatide in the astroglial culture and showed decreasing amounts of sulfatide with days in vitro; 0.05 nmol/mg protein at day 10 and 0.01 nmol/mg protein at day 17. This analysis also showed that neither sulpholactosylceramide nor seminolipid was present, each of which also has affinity for the sulfatide antibody. This selective and intracellular expression encourages further identification of the astrocytes expressing sulfatide and the biological role of sulfatide in these cells.