Protein hydrophobicity and stability support the thermodynamic theory of protein degradation

Rat liver enzymes were used to study the relationship between their in vivo half-lives and their apparent hydrophobicity or their resistance to inactivation by mechanical shaking. The apparent hydrophobicity of these enzymes, measured as the percent of the protein recovered from an octyl-Sepharose column, is correlated with their known half-lives (r = 0.75, P less than 0.01). The presence of specific ligands which are known to increase compactness by impeding unfolding of proteins decreased the apparent hydrophobicity of fructose-1,6-bisphosphatase, glucose-6-phosphate dehydrogenase, glyceraldehyde-3-phosphate dehydrogenase, and pyruvate kinase. Resistance of enzymes to inactivation by mechanical shaking correlated well with their in vivo half-lives (r = 0.90, P less than 0.01). When the shaking experiments were done in the presence of substrates, fructose-1,6-bisphosphatase, glucose-6-phosphate dehydrogenase, glyceraldehyde-3-phosphate dehydrogenase and lactate dehydrogenase were protected from inactivation.     

A New Method to Characterize Hydrophobic Organization of Proteins: Application to Rational Protein Engineering of Barnase

Abstract

We present a new algorithm for characterization of protein spatial structure basing on the molecular hydrophobicity potential approach. The method is illustrated by the analysis of three-dimensional structure of barnase and barnase-barstar complex. Current approach enables identification of amino acid residues situated in unfavorable environment (these residues may be “active” for binding), and to map quantitatively hydrophobic, hydrophilic and unfavorable hydrophobic-hydrophilic intra-and inter-molecular contacts involving backbone and side-chain segments of amino acid residues. Calculation of individual contributions of amino acid residues to such contacts permits identification of structurally-important residues. The contact plots obtained with molecular hydrophobicity potential calculations, provide easy rules to choose sites for mutations, which can increase a strength of intra- or inter-molecular hydrophobic interactions. The unfavorable hydrophobic-hydrophilic contact can be mutated to favorable hydrophobic, and already existing weak hydrophobic contact can be strengthen by increasing hydrophobicity of residues in contact. Basing on the analysis of the contact plots, we suggest several mutations of barnase which are supposed to increase intramolecular hydrophobic interactions, and thus might lead to increased stability of the protein. Part of these mutations was studied previously experimentally, and indeed stabilized barnase. The other of predicted mutations were not studied experimentally yet. Several new mutations of barnase and barstar are also proposed to enhance the hydrophobic interactions on their binding interface.     

Probe into a multi-protein prokaryotic organelle using thermal scanning assay reveals distinct properties of the core and the shell

BackgroundBacterial microcompartments represent the only reported category of prokaryotic organelles that are capable of functioning as independent bioreactors. In this organelle, a biochemical pathway with all the enzyme machinery is encapsulated within an all protein shell. The shell proteins and the enzymes have distinct structural features. It is hypothesized that flat shell proteins align sideways to form extended sheets and, the globular enzymes fill up the central core of the organelle.MethodsUsing differential scanning fluorimetry, we explored the structure and functional alteration of Pdu BMC, involving tertiary or quaternary structures.ResultsOur findings exhibit that these intact BMCs as a whole behave similar to a globular protein with a rich hydrophobic core, which is exposed upon thermal insult. The encapsulated enzymes itself have a strong hydrophobic core, which is in line with the hydrophobic-collapse model of protein folding. The shell proteins, on the other hand, do not have a strong hydrophobic core and show a significant portion of exposed hydrophobic patches.ConclusionWe show for the first time the thermal unfolding profile of the BMC domain proteins and the unique exposure of hydrophobic patches in them might be required for anchoring the enzymes leading to better packaging of the micro-compartments.General significanceThese observations indicate that the genesis of these unique bacterial organelles is driven by the hydrophobic interactions between the shell and the enzymes. Insights from this work will aid in the genetic and biochemical engineering of thermostable efficient enzymatic biomaterials.     

Thermal stability of proteins in the presence of poly(ethylene glycols)

Thermal unfolding of ribonuclease, lysozyme, chymotrypsinogen, and beta-lactoglobulin was studied in the absence or presence of poly(ethylene glycols). The unfolding curves were fitted to a two-state model by a nonlinear least-squares program to obtain values of delta H, delta S, and the melting temperature Tm. A decrease in thermal transition temperature was observed in the presence of poly(ethylene glycol) for all of the protein systems studied. The magnitude of such a decrease depends on the particular protein and the molecular size of poly(ethylene glycol) employed. A linear relation can be established between the magnitude of the decrease in transition temperature and the average hydrophobicity of these proteins; namely, the largest observable decrease is associated with the protein of the highest hydrophobicity. Further analysis of the thermal unfolding data reveals that poly(ethylene glycols) significantly effect the relation between delta H degrees of unfolding and temperature for all the proteins studied. For beta-lactoglobulin, a plot of delta H versus Tm indicates a change in slope from a negative to a positive value, thus implying a change in delta Cp in thermal unfolding caused by the presence of poly(ethylene glycols). Results from solvent-protein interaction studies indicate that at high temperature poly(ethylene glycol) 1000 preferentially interacts with the denatured state of protein but is excluded from the native state at low temperature. These observations are consistent with the fact that poly(ethylene glycols) are hydrophobic in nature and will interact favorably with the hydrophobic side chains exposed upon unfolding; thus, it leads to a lowering of thermal transition temperature.     

Comparison of the structures and stabilities of coiled‐coil proteins containing hexafluoroleucine and t‐butylalanine provides insight into the stabilizing effects of highly fluorinated amino acid side‐chains

Highly fluorinated analogs of hydrophobic amino acids are well known to increase the stability of proteins toward thermal unfolding and chemical denaturation, but there is very little data on the structural consequences of fluorination. We have determined the structures and folding energies of three variants of a de novo designed 4‐helix bundle protein whose hydrophobic cores contain either hexafluoroleucine (hFLeu) or t‐butylalanine (tBAla). Although the buried hydrophobic surface area is the same for all three proteins, the incorporation of tBAla causes a rearrangement of the core packing, resulting in the formation of a destabilizing hydrophobic cavity at the center of the protein. In contrast, incorporation of hFLeu, causes no changes in core packing with respect to the structure of the nonfluorinated parent protein which contains only leucine in the core. These results support the idea that fluorinated residues are especially effective at stabilizing proteins because they closely mimic the shape of the natural residues they replace while increasing buried hydrophobic surface area.     

Protein instability during HIC: hydrogen exchange labeling analysis and a framework for describing mobile and stationary phase effects

Unfolding of marginally stable proteins is a significant factor in commercial application of hydrophobic interaction chromatography (HIC). In this work, hydrogen-deuterium isotope exchange labeling has been used to monitor protein unfolding on HIC media for different stationary phase hydrophobicities and as a function of ammonium sulfate concentration. Circular dichroism and Raman spectroscopy were also used to characterize the structural perturbations experienced by solution phase protein that had been exposed to media and by protein adsorbed on media. As expected, greater instability is seen on chromatographic media with greater apparent hydrophobicity. However, increased salt concentrations also led to more unfolding, despite the well-known stabilizing effect of ammonium sulfate in solution. A thermodynamic framework is proposed to account for the effects of salt on both adsorption and stability during hydrophobic chromatography. Using appropriate estimates of input quantities, analysis with the framework can explain how salt effects on stability in chromatographic systems may contrast with solution stability.     

Hydrophobicity and stability for a family of model proteins

According to the conventional definition, the hydrophobic effect is a result of thermodynamic changes occurring when a nonpolar group dissolves in water and attributable to the fact that water in contact with such a group has special structural and energetic properties. Disagreement now exists as to whether this effect promotes or hinders protein denaturation. Taking the heat capacity change of unfolding as a measure of the hydrophobicity of the protein interior, others have shown that protein stabilities are systematically affected by changes in hydrophobicity. It has been suggested that the observed trends show that hydrophobic hydration is intrinsically a destabilizing factor. Model calculations using known equations for the stability curves and certain simplifying assumptions now show that such regularities provide no evidence for or against this conclusion. All available data can be rationalized if hydrophobic terms are evaluated from models that require a positive hydrophobic contribution to the Gibbs energy of unfolding. The calculations also confirm the recent finding that any set of proteins with denaturation temperatures between about 330 and 380 K that exhibits entropy convergence at about 386 K is thermodynamically required to show enthalpy convergence at approximately the same temperature. © 1993 John Wiley & Sons, Inc.     

Predicting Atomic Details of the Unfolding Pathway for YibK,a Knotted Protein from the SPOUT Superfamily

Abstract

Several protein structures have been reported to contain intricate knots of the polypeptide backbone but the mechanism of the (un)folding process of knotted proteins remains unknown. The members of the SPOUT superfamily of RNA methyltransferases are some of the most intensely studied systems for investigation of the knot formation and function. YibK (whose biochemical function remains unknown) is the representative protein of the SPOUT superfamily. This protein exhibits a deep trefoil knot at the C-terminus.

We conducted an extensive computational analysis of the unfolding process for the monomeric form of YibK. In order to predict the (un)folding pathway of YibK, we have calculated the order of secondary structure disassembly using UNFOLD, and performed thermal unfolding simulations using classical Molecular Dynamics (MD), as well as simulations employing reduced representation of the peptide chain using either MD with the UNRES method or the Monte Carlo (MC) unfolding with the REFINER method.

Results obtained from all methods used in this work are in qualitative agreement. We found that YibK unfolds through four intermediate states. The trefoil knot in YibK disappears at the end of the unfolding process, long after the protein loses its native topology. We observed that the C-terminus leaves the knotting loop folded into a hairpin-like structure, in agreement with the results of coarse-grained simulation reported earlier. We propose that the folding pathway of YibK corresponds to the reversed sequence of events observed in the unfolding pathway elucidated in this study. Thus, we predict that the knot formation is the slowest part of the YibK folding process.     

Underlying Hydrophobic Sequence Periodicity of Protein Tertiary Structure

Abstract

Hydropathy plots or window averages over local stretches of the sequence of residue hydrophobicity have revealed patterns related to various protein tertiary structural features. This has enabled identification of regions of the sequence that are at the surface or within the interior of globular soluble proteins, regions located within the lipid bilayer of transmembrane proteins, portions of the sequence that characterize repeating motifs, as well as motifs that usefully characterize different protein structural families. This, therefore, provides one example of the generally expressed maxim that “sequence determines structure”. On the other hand, a number of previous investigations have shown the rapidly varying values of residue hydrophobicity along the sequence to be distributed approximately randomly. So one might question just how much of the sequence actually determines structure. It is, therefore, of interest to extract that part of this rapidly varying distribution of residue hydrophobicity that is responsible for the longer wavelength variations that correlate with protein tertiary structural features and to determine their prevalence within the entire distribution. This is accomplished by a finite Fourier analysis of the sequence of residue hydrophobicity and of a new measure of residue distance from the protein interior. Calculations are performed on a number of globins, immunoglobulins, cuprodoxins, and papain-like structures. The spectral power of the Fourier amplitudes of the frequencies extracted, whose inverse transforms underlie the windowed values of residue hydrophobicity is shown to be a small fraction of the total power of the hydrophobicity distribution and thereby consistent with a distribution that might appear to be predominantly random. The wide range of sequence identity between proteins having the same fold, all exhibiting similar small fractions of power amplitude that correlate with the longer wavelength inside-to- outside excursions of the amino acid residues, supports the general contention that close sequence identity is an expression of a close evolutionary relationship rather than an expression of structural similarity. Practical implications of the present analysis for protein structure prediction and engineering are also described.     

Solid-phase iodination: in vitro labeling with 125I of proteins bound to nitrocellulose and substituted Sepharose

A solid-phase iodination method is described which employs either nitrocellulose paper, phenyl- and octyl-Sepharose beads, or octyl hydroxylapatite as matrices to adsorb proteins. Nitrocellulose lends itself to cases where denaturation of the iodinated proteins due to the use of chaotropic reagents or strong acids for protein elution can be tolerated. On the other hand, substituted Sepharoses, preferably octyl-Sepharose, should be used when preservation of the biological activity of the iodinated protein molecules is required; immunoglobulins and protein A, for instance, could be recovered as functionally active molecules because they were extracted from the hydrophobic matrices under nondenaturing conditions. Both methods are advantageous if, for example, series of fractions from column chromatographies (including HPLC) are to be iodinated and subsequently analyzed by gel electrophoresis or bioassays. Furthermore, the amount of radioactive waste can be reduced considerably.     

The Unfolding MD Simulations of Cyclophilin: Analyzed by Surface Contact Networks and Their Associated Metrics

Currently, considerable interest exists with regard to the dissociation of close packed aminoacids within proteins, in the course of unfolding, which could result in either wet or dry moltenglobules. The progressive disjuncture of residues constituting the hydrophobic core ofcyclophilin from L. donovani (LdCyp) has been studied during the thermal unfolding of the molecule, by molecular dynamics simulations. LdCyp has been represented as a surface contactnetwork (SCN) based on the surface complementarity (S_m) of interacting residues within themolecular interior. The application of S_m to side chain packing within proteins make it a very sensitive indicator of subtle perturbations in packing, in the thermal unfolding of the protein. Network based metrics have been defined to track the sequential changes in the disintegration ofthe SCN spanning the hydrophobic core of LdCyp and these metrics prove to be highly sensitive compared to traditional metrics in indicating the increased conformational (and dynamical) flexibility in the network. These metrics have been applied to suggest criteria distinguishing DMG, WMG and transition state ensembles and to identify key residues involved in crucial conformational/topological events during the unfolding process.     

Unfolding Studies of <Emphasis Type="Italic">Escherichia coli</Emphasis> Maltodextrin Glucosidase Monitored by Fluorescence Spectroscopy

Equilibrium unfolding of a 69-kDa monomeric Escherichia coli maltodextrin glucosidase (MalZ) was studied using intrinsic and extrinsic fluorescence spectroscopy. The unfolding transition of MalZ followed a three-state process, involving the formation of a stable intermediate state having more exposed hydrophobic surface. It was found that the protein structure can be easily perturbed by low concentration of guanidium hydrochloride (GdnHCl) and, at a GdnHCl concentration of 2 M, MalZ was denatured completely. The active site of the protein also has been proved to be sensitive to a low concentration of GdnHCl since MalZ deactivated at 0.5 M GdnHCl completely. The surface hydrophobicity and ANS-binding site of the protein have been determined to be 150.7 and 0.24, respectively. Perhaps the formation of the stable unfolding intermediate, having higher surface hydrophobicity, may be one of the reasons for aggregation of MalZ and its recognition by chaperonin GroEL during the assisted folding pathway.     

Role of hydrophobic surface proteins in mediating adherence of group B streptococci to epithelial cells.

Determination of the cell-surface hydrophobicity of group B streptococci by hydrophobic interaction chromatography on phenyl-Sepharose revealed that human and bovine group B streptococcal isolates with protein surface antigens, either alone or in combination with polysaccharide antigens, were mainly hydrophobic, whereas those with polysaccharide antigens alone were mainly hydrophilic. Removal of capsular neuraminic acid enhanced, and pronase treatment reduced, surface hydrophobicity. The hydrophobic surface proteins, solubilized by mutanolysin treatment of the bacteria and isolated by hydrophobic interaction chromatography, appeared in SDS-PAGE as numerous protein bands. Staphylococcal carrier cells loaded with antibodies produced against hydrophobic surface proteins agglutinated specifically with hydrophobic group B streptococci. No agglutination reaction was observed with hydrophilic cultures. Hydrophobic group B streptococci adhered to buccal epithelial cells in significantly higher numbers than did hydrophilic cultures. The adherence of group B streptococci to epithelial cells was inhibited in the presence of isolated hydrophobic proteins and in the presence of specific antibodies produced against hydrophobic proteins. The results of this study demonstrate a close relation between the occurrence of type-specific antigens, surface hydrophobicity and the adherence of group B streptococci to epithelial cells.     

Cloning and Sequence Analysis of the Structural Gene for the bc 1 -Type Rieske Iron-Sulfur Protein from Thermus thermophilus HB8

The structural gene encoding the Rieske iron-sulfur protein from Thermus thermophilus HB8 has been cloned and sequenced. The gene encodes a protein of 209 amino acids that begins with a hydrophilic N-terminus followed by a stretch of 21 hydrophobic amino acids that could serve as a transmembrane helix. The remainder of the protein has a hydrophobicity pattern typical of a water-soluble protein. A phylogenetic analysis of 26 Rieske proteins that are part of bc ₁ or b ₆ f complexes shows that they fall into three major groups: eubacterial and mitochondrial, cyanobacterial and plastid, and five highly divergent outliers, including that of Thermus. Although the overall homology with other Rieske proteins is very low, the C-terminal half of the Thermus protein contains the signature sequence CTHLGC-(13X)-CPCH that most likely provides the ligands of the [2Fe-2S] cluster. It is proposed that this region of the protein represents a small domain that folds independently and that the encoding DNA sequence may have been transferred during evolution to several unrelated genes to provide the cluster attachment site to proteins of different origin. The role of individual residues in this domain of the Thermus protein is discussed vis-a-vis the three-dimensional structure of the bovine protein (Iwata et al., 1996 Structure 4, 567–579).     

Influence of pH on the solubility and conformational characteristics of muscle proteins from mullet (Mugil cephalus)

Mullet (Mugil cephalus) muscle homogenates were adjusted to different pH ranging from 2 to12 and the proteins extracted were evaluated for changes in solubility and conformational characteristics viz. surface hydrophobicity and reactive sulphydryl groups. Altering the pH of muscle homogenate to acidic or alkaline increased protein solubility. The hydrophobicity of the proteins increased on exposure to extreme pH indicating unfolding. The reactive sulphydryl groups decreased at acidic and alkaline pH with the lowest at pH 4. When the pH of the muscle homogenates was brought back to the original pH (6.3), the protein solubility was found to decrease. Reactive sulphydryl groups and ANS hydrophobicity of the proteins increased on readjusting the pH resulting in a molten-globule state. The electrophortogram of the samples corresponded well with the observations. Alterations in functional properties of these modified proteins are an area of interest for commercial application.     

Insight into the stability of the hydrophobic binding proteins of Escherichia coli: assessing the proteins for use as biosensors

Spectroscopic methods were used to monitor the unfolding of the leucine specific (LS) and the leucine-isoleucine-valine (LIV) binding proteins. Our studies indicate that ligand-free protein undergoes a simple two-state unfolding, whereas the protein-ligand complex undergoes a three-state unfolding model. Ligand binding causes significant stabilization of both proteins. There is correlation between ligand hydrophobicity and protein stabilization: the most hydrophobic ligand, isoleucine, causes the most significant stabilization of LIV protein. A disulfide bond present in N-domain of both proteins makes a large contribution to the protein stability of these periplasmic binding receptors.     

A Comprehensive Multidimensional-Embedded, One-Dimensional Reaction Coordinate for Protein Unfolding/Folding

The goal of the Dynameomics project is to perform, store, and analyze molecular dynamics simulations of representative proteins, of all known globular folds, in their native state and along their unfolding pathways. To analyze unfolding simulations, the location of the protein along the unfolding reaction coordinate (RXN) must be determined. Properties such as the fraction of native contacts and radius of gyration are often used; however, there is an issue regarding degeneracy with these properties, as native and nonnative species can overlap. Here, we used 15 physical properties of the protein to construct a multidimensional-embedded, one-dimensional RXN coordinate that faithfully captures the complex nature of unfolding. The unfolding RXN coordinates for 188 proteins (1534 simulations and 22.9 μs in explicit water) were calculated. Native, transition, intermediate, and denatured states were readily identified with the use of this RXN coordinate. A global native ensemble based on the native-state properties of the 188 proteins was created. This ensemble was shown to be effective for calculating RXN coordinates for folds outside the initial 188 targets. These RXN coordinates enable, high-throughput assignment of conformational states, which represents an important step in comparing protein properties across fold space as well as characterizing the unfolding of individual proteins.     

Design and Structure of an Equilibrium Protein Folding Intermediate: A Hint into Dynamical Regions of Proteins

Partly unfolded protein conformations close to the native state may play important roles in protein function and in protein misfolding. Structural analyses of such conformations which are essential for their fully physicochemical understanding are complicated by their characteristic low populations at equilibrium. We stabilize here with a single mutation the equilibrium intermediate of apoflavodoxin thermal unfolding and determine its solution structure by NMR. It consists of a large native region identical with that observed in the X-ray structure of the wild-type protein plus an unfolded region. Small-angle X-ray scattering analysis indicates that the calculated ensemble of structures is consistent with the actual degree of expansion of the intermediate. The unfolded region encompasses discontinuous sequence segments that cluster in the 3D structure of the native protein forming the FMN cofactor binding loops and the binding site of a variety of partner proteins. Analysis of the apoflavodoxin inner interfaces reveals that those becoming destabilized in the intermediate are more polar than other inner interfaces of the protein. Natively folded proteins contain hydrophobic cores formed by the packing of hydrophobic surfaces, while natively unfolded proteins are rich in polar residues. The structure of the apoflavodoxin thermal intermediate suggests that the regions of natively folded proteins that are easily responsive to thermal activation may contain cores of intermediate hydrophobicity.