Coexistence of Flexibility and Stability of Proteins: An Equation of State

We consider a recently suggested “equation of state” for natively folded proteins, and verify its validity for a set of about 5800 proteins. The equation is based on a fractal viewpoint of proteins, on a generalization of the Landau-Peierls instability, and on a marginal stability criterion. The latter allows for coexistence of stability and flexibility of proteins, which is required for their proper function. The equation of state relates the protein fractal dimension , its spectral dimension , and the number of amino acids N. Using structural data from the protein data bank (PDB) and the Gaussian network model (GNM), we compute and for the entire set and demonstrate that the equation of state is well obeyed. Addressing the fractal properties and making use of the equation of state may help to engineer biologically inspired catalysts.     

Origin of the anomalous circular dichroism spectra of many apomyoglobin mutants

Several authors have reported that many sperm whale apomyoglobin mutants show anomalous circular dichroism spectra. These mutants have a low molar ellipticity compared to the wild-type protein but in several cases have the same stability of unfolding. A model in which native apomyoglobin is not folded in the same manner as that in other proteins and in which mutants show progressive reductions in their degree of folding has been suggested to explain this phenomenon. However, nuclear magnetic resonance of the native apomyoglobin conformation has shown that this state is folded and compact, raising the possibility that the anomalous circular dichroism spectra could have another explanation. We studied several mutants with anomalous circular dichroism spectra and found that these proteins were all contaminated with nucleic acid that contributed to the ultraviolet absorption and caused uncertainty in the determination of protein concentration. The resulting overestimation of the concentration of apomyoglobin explains the phenomenon of anomalous circular dichroism spectra. We describe a procedure to remove the contaminant nucleic acid which yields accurate protein concentration measurements and provides the normal circular dichroism spectra. Our findings support a well-structured native conformation for apomyoglobin and may also be of the interest to scientists working with the purification of recombinant proteins.     

N- and C-terminal fragments of a globular protein constructed by elongation of modules as a units associated for functional complementation

We have been interested in partially folded proteins with marginal stability and activity, because they have a potential to be mature proteins by artificial evolution. A module is defined as a contiguous peptide chain forming a compact region in a globular protein. Modules may be used as building blocks to create partially folded proteins. Barnase, a ribonuclease consisting of 110 amino acids, has been divided into six modules (M1-M6), four peptide fragments, M12 (1-52), M123 (1-73), M1234 (1-88) and M12345 (1-98), have been constructed by progressive elongation of the modules from the N-terminus. Only M12345 (1-98) had a partially folded conformation, but it lacked detectable RNase activity. A mixture of M12345 (1-98) with M56 (89-110) showed weak but distinct RNase activity. Unfolded M12345 (1-96) was constructed by removal of two residues from the C-terminus of M12345 (1-98). The mixture of M12345 (1-96) with M56 (89-110) also showed RNase activity. Further, the interaction endowed M12345 (1-96) with conformational stability. We propose that N- and C-terminal fragments obtained by successive elongation of modules would interact to be a complex with marginal stability and activity, which would be used for creating a mature complex by artificial evolution.     

NOE data demonstrating a compact unfolded state for an SH3 domain under non-denaturing conditions.

The N-terminal SH3 domain of drk (drkN SH3) is unstable, existing in equilibrium between a folded state (Fexch) and an unfolded state (Uexch) under non-denaturing buffer conditions. Using a15N/2H-labeled sample, long range amide NOEs can be observed in the Uexchstate as a result of reduced relaxation, in some cases correlating protons over 40 residues apart. These long range NOEs disappear upon addition of 2 M guanidinium chloride, demonstrating that there are substantial differences between the Uexchand the guanidine denatured states. Calculations using the long range NOEs of the Uexchstate yield highly compact structures having non-native turns and a non-native buried tryptophan residue. These structures agree with experimental stopped-flow fluorescence data and analytical ultracentrifugation results. Since protein stability depends on the structural and dynamic properties of both the folded and unfolded states, this study provides insights into the stability of the drkN SH3 domain. These results provide the first strong NOE-based evidence for compact unfolded states of proteins and suggest that some unfolded states under physiological conditions have specific interactions leading to compact structures.     

Surface hydrophobic groups, stability, and flip-flopping in lattice proteins

Two-dimensional lattice protein models were studied in two approximations of the conformational equilibrium to elucidate the role of surface hydrophobic groups in their stabilities. We demonstrate that stability of any compactly folded sequence is determined by its ability to "flip-flop" (refold) into alternative compact structures. The degree of stability required for folded sequences determines the average numbers of surface hydrophobic groups in stable lattice structures which are in good agreement with ratios of core to surface hydrophobic groups in real proteins. However, the average destabilization of the native structure per surface hydrophobic group is small (0-0.25 kcal/mol), often disagrees with the free energies derived from the ratios of core to surface hydrophobic groups in the same structures, and has a combinatorial entropic nature independent of the strength of structure stabilizing interactions. This suggests that the free energies derived from the core to surface ratios of hydrophobic groups in real proteins have little to do with folding thermodynamics. On average, sequences with highly stable native structures are the least hydrophobic. The results suggest that in designing novel stable proteins hydrophobic groups on the surface should be avoided to reduce the possibility of flip-flopping. The average stability of highly designable structures is never higher than that of some low designability structures, contrary to the accepted view. In the equilibrium approximation with alternative compact and partially unfolded structures, the requirement of high stability selects a unique 5 x 5 structure formed by only a few sequences, suggesting much stronger sequence selectivity than commonly thought.     

Compaction of ribosomal protein S6 by sucrose occurs only under native conditions

The effect of osmolyte sucrose on the stability and compaction of the folded and unfolded states of ribosomal protein S6 from Thermus thermophilus was analyzed. Confirming previous results obtained with sodium sulfate and trehalose, refolding stopped-flow measurements of S6 show that sucrose favors the conversion of the unfolded state ensemble to a highly compact structure (75% as compact as the folded state). This conversion occurs when the unfolded state is suddenly placed under native conditions and the compact state accumulates in a transient off-folding pathway. This effect of sucrose on the compaction of the unfolded state ensemble is counteracted by guanidinium hydrochloride. The compact state does not accumulate at higher guanidinium concentrations and the unfolded state ensemble does not display increased compaction in the presence of 6 M guanidinium as evaluated by collisional quenching of tryptophan fluorescence. In contrast, accessibility of the tryptophan residue of folded S6 above 1 M sucrose concentration decreased as a result of an increased compaction of the folded state. Unfolding stopped-flow measurements of S6 reflect this increased compaction of the folded state, but the unfolding pathway is not affected by sucrose. Compaction of folded and unfolded S6 induced by sucrose occurs under native conditions indicating that decreased protein conformational entropy significantly contributes to the mechanism of protein stabilization by osmolytes.     

An electrostatic basis for the stability of thermophilic proteins

Two factors provide key contributions to the stability of thermophilic proteins relative to their mesophilic homologues: electrostatic interactions of charged residues in the folded state and the dielectric response of the folded protein. The dielectric response for proteins in a "thermophilic series" globally modulates the thermal stability of its members, with the calculated dielectric constant for the protein increasing from mesophiles to hyperthermophiles. This variability results from differences in the distribution of charged residues on the surface of the protein, in agreement with structural and genetic observations. Furthermore, the contribution of electrostatic interactions to the stability of the folded state is more favorable for thermophilic proteins than for their mesophilic homologues. This leads to the conclusion that electrostatic interactions play an important role in determining the stability of proteins at high temperatures. The interplay between electrostatic interactions and dielectric response also provides further rationalization for the enhanced stability of thermophilic proteins with respect to cold-denaturation. Taken together, the distribution of charged residues and their fluctuations have been shown to be factors in modulating protein stability over the entire range of biologically relevant temperatures.     

Charge-charge interactions influence the denatured state ensemble and contribute to protein stability

Several recent studies have shown that it is possible to increase protein stability by improving electrostatic interactions among charged groups on the surface of the folded protein. However, the stability increases are considerably smaller than predicted by a simple Coulomb's law calculation, and in some cases, a charge reversal on the surface leads to a decrease in stability when an increase was predicted. These results suggest that favorable charge-charge interactions are important in determining the denatured state ensemble, and that the free energy of the denatured state may be decreased more than that of the native state by reversing the charge of a side chain. We suggest that when the hydrophobic and hydrogen bonding interactions that stabilize the folded state are disrupted, the unfolded polypeptide chain rearranges to compact conformations with favorable long-range electrostatic interactions. These charge-charge interactions in the denatured state will reduce the net contribution of electrostatic interactions to protein stability and will help determine the denatured state ensemble. To support this idea, we show that the denatured state ensemble of ribonuclease Sa is considerably more compact at pH 7 where favorable charge-charge interactions are possible than at pH 3, where unfavorable electrostatic repulsion among the positive charges causes an expansion of the denatured state ensemble. Further support is provided by studies of the ionic strength dependence of the stability of charge-reversal mutants of ribonuclease Sa. These results may have important implications for the mechanism of protein folding.     

Global and Local Mobility of Apocalmodulin Monitored through Fast-Field Cycling Relaxometry

Calmodulin (CaM) is a ubiquitous eukaryotic protein with two conformationally independent domains that can bind up to two calcium ions each. In the calcium-bound state, CaM is able to regulate a vast number of cellular activities by binding to a multiplicity of target proteins in different modes. Its versatility has been ascribed to its anomalously high flexibility. The calcium-free form (apoCaM), which is the resting state of CaM in cells, is also able to functionally bind a number of protein targets, but its dynamics has received less attention. At variance with the calcium-bound form, the crystal structure of apoCaM shows a compact organization of the two domains, but NMR measurements could not detect any contact between them, thus indicating the presence of mobility in solution. The mobility of apoCaM is here investigated through protein proton relaxation rate measurements performed with a high-sensitivity fast-field cycling relaxometer. Such measurements provide direct access to the spectral density function and show that 1), the reorientation time is in agreement with a closed form of the protein; but 2), the collective order parameter is much smaller than for other well folded compact proteins, indicating that a remarkably large side-chain mobility must be present.     

The use of osmolytes to facilitate protein NMR spectroscopy

Summary A method of stabilizing folded proteins is described, which allows NMR studies under conditions where a protein would normally be unfolded. This enables stable proteins to be examined at elevated temperatures, or spectra recorded on samples that are insufficiently stable under normal conditions. Up to two molar perdeuterated glycine, a potent osmolyte, can be added to aqueous protein NMR samples without altering the folded three-dimensional structure or function of the protein. However, the stability of the folded form is dramatically increased. This is illustrated for the protein lysozyme at high temperature (348 K) where the structural integrity is destroyed in standard aqueous solution, but is retained in the osmolyte solution. We hope that the technique will be of value to those studying by NMR the structural biology of protein fragments and mutants, which are often of reduced stability compared with the original proteins.To whom correspondence should be addressed.     

Irreversible unfolding of the neutral pH form of influenza hemagglutinin demonstrates that it is not in a metastable state

The thermal denaturation of the proteins of influenza virus has been measured by differential scanning calorimetry in the presence and absence of lipids as a function of scan rate. We have applied theories of irreversible thermodynamics to obtain the activation energy. In the presence of liposomes of dioleoylphosphatidylcholine with the ganglioside, GD(1a), the denaturation temperature of the hemagglutinin protein is lowered. This lowering of thermal stability is also reflected in the temperature dependence of the circular dichroism spectra. Quasi-elastic light scattering confirms that liposomes containing GD(1a) interact with the virus and inhibit the growth in the size of the particle as a function of temperature. Although the virus can fuse with the liposomes at higher temperatures, the enthalpy change for this process is not detectable. Our results also demonstrate that the compact folded structure of the influenza hemagglutinin protein is not a kinetically trapped metastable high-energy form.     

Analysis of the stability of multimeric proteins by effective DeltaG and effective m-values

Analyzing the stability of a multimeric protein is challenging because of the intrinsic difficulty in handling the mathematical model for the folded multimer-unfolded monomer equilibrium. To circumvent this problem, we introduce the concept of effective stability, DeltaGeff (= -RTlnKeff), where Keff is the equilibrium constant expressed in monomer units. Analysis of the denaturant effect on DeltaGeff gives new insight into the stability of multimeric proteins. When a multimeric protein is mostly folded, the dependence of effective stability on denaturant concentration (effective m-value) is simply the m-value of its monomeric unit. However, when the protein is mostly unfolded, its stability depends on denaturant concentration with the m-value of its multimeric form. We also find that the effective m-value at the Cm is a good approximation of the apparent m-value determined by fitting the equilibrium unfolding data from multimeric proteins with a two-state monomer model. Moreover, when the m-value of a monomeric unit is estimated from its size, the effective stability of a multimeric protein can be determined simply from Cm and this estimated m-value. These simple and intuitive approaches will allow a facile analysis of the stability of multimeric proteins. These analyses are also applicable for high-throughput analysis of protein stability on a proteomic scale.     

Selection of stably folded proteins by phage-display with proteolysis.

To facilitate the process of protein design and learn the basic rules that control the structure and stability of proteins, combinatorial methods have been developed to select or screen proteins with desired properties from libraries of mutants. One such method uses phage-display and proteolysis to select stably folded proteins. This method does not rely on specific properties of proteins for selection. Therefore, in principle it can be applied to any protein. Since its first demonstration in 1998, the method has been used to create hyperthermophilic proteins, to evolve novel folded domains from a library generated by combinatorial shuffling of polypeptide segments and to convert a partially unfolded structure to a fully folded protein.     

Structural characterization of partially folded intermediates of apomyoglobin H64F

We present a detailed investigation of unfolded and partially folded states of a mutant apomyoglobin (apoMb) where the distal histidine has been replaced by phenylalanine (H64F). Previous studies have shown that substitution of His64, located in the E helix of the native protein, stabilizes the equilibrium molten globule and native states and leads to an increase in folding rate and a change in the folding pathway. Analysis of changes in chemical shift and in backbone flexibility, detected via [1H]-15N heteronuclear nuclear Overhauser effect measurements, indicates that the phenylalanine substitution has only minor effects on the conformational ensemble in the acid- and urea-unfolded states, but has a substantial effect on the structure, dynamics, and stability of the equilibrium molten globule intermediate formed near pH 4. In H64F apomyoglobin, additional regions of the polypeptide chain are recruited into the compact core of the molten globule. Since the phenylalanine substitution has negligible effect on the unfolded ensemble, its influence on folding rate and stability comes entirely from interactions within the compact folded or partly folded states. Replacement of His64 with Phe leads to favorable hydrophobic packing between the helix E region and the molten globule core and leads to stabilization of helix E secondary structure and overall thermodynamic stabilization of the molten globule. The secondary structure of the equilibrium molten globule parallels that of the burst phase kinetic intermediate; both intermediates contain significant helical structure in regions of the polypeptide that comprise the A, B, E, G, and H helices of the fully folded protein.     

Non‐additive stabilization by halogenated amino acids reveals protein plasticity on a sub‐angstrom scale

It has been a long‐standing goal to understand the structure‐stability relationship of proteins, as optimal stability is essential for protein function and highly desirable for protein therapeutics. Halogenation has emerged as a minimally invasive strategy to probe the physical characteristics of proteins in solution, as well as enhance the structural stabilities of proteins for therapeutic applications. Although advances in synthetic chemistry and genetic code expansion have allowed for the rapid synthesis of proteins with diverse chemical sequences, much remains to be learned regarding the impact of these mutations on their structural integrity. In this contribution, we present a systematic study of three well‐folded model protein systems, in which their structural stabilities are assessed in response to various hydrogen‐to‐halogen atom mutations. Halogenation allows for the perturbation of proteins on a sub‐angstrom scale, offering unprecedented precision of protein engineering. The thermodynamic results from these model systems reveal that in certain cases, proteins can display modest steric tolerance to halogenation, yielding non‐additive consequences to protein stability. The observed sub‐angstrom sensitivity of protein stability highlights the delicate arrangement of a folded protein core structure. The stability data of various halogenated proteins presented herein should also provide guidelines for using halogenation as a strategy to improve the stability of protein therapeutics.     

8-anilino-1-naphthalene sulfonic acid (ANS) induces folding of acid unfolded cytochrome c to molten globule state as a result of electrostatic interactions.

Hydrophobic interaction of 8-anilino-1-naphthalene sulfonic acid (ANS) with proteins is one of the widely used methods for characterizing/detecting partially folded states of proteins. We have carried out a systematic investigation on the effect of ANS, a charged hydrophobic fluorescent dye, on structural properties of acid-unfolded horse heart cytochrome c at pH 2.0 by a combination of optical methods and electrospray ionization mass spectroscopy (ESI MS). ANS was found to induce, a secondary structure similar to native protein and quenching of fluorescence of tryptophan residue, in the acid-unfolded protein. However, the tertiary structure was found to be disrupted thus indicating that ANS stabilizes a molten globule state in acid-unfolded protein. To understand the mechanism of ANS-induced folding of acid-unfolded cytochrome c, comparative ESI MS, soret absorption, and tryptophan fluorescence studies using nile red, a neutral hydrophobic dye, and ANS were carried out. These studies suggested that, at low pH, electrostatic interactions between negatively charged ANS molecules and positively charged amino acid residues present in acid-unfolded cytochrome c are probably responsible for ANS-induced folding of acid-unfolded protein to partially folded compact state or molten globule state. This is the first experimental demonstration of ANS induced folding of unfolded protein and puts to question the usefulness of ANS for characterization/determination of partially folded intermediates of proteins observed under low pH conditions.     

Effect of salt on the urea-unfolded form of barstar probed by m value measurements

To probe for residual structure present in the urea-unfolded form of the small protein barstar, to determine how salt might modulate such structure, and to determine how such structure might affect the stability of the protein, mutant variants that display m values different from that of the wild-type protein have been studied. The mutant proteins were obtained by site-directed mutagenesis at residue positions located on the surface of the folded protein. The m value, which represents the preferential free energy of interaction of urea with the unfolded form in comparison to that with the folded state, was determined from equilibrium urea-induced unfolding curves. Mutant proteins for which the m values were significantly greater than (m(+) mutant forms), significantly smaller than (m(-) mutant forms), or similar to (m(0) mutant forms) the m value determined for the wild-type protein were studied. The unfolded forms of the m(0), m(+) and m(-) mutant proteins represent different components within the unfolded form ensemble, which differ from each other in their solvent-exposed surface areas. Hence, the m value has been used as a measure of residual structure in the unfolded form. To further understand the nature of structures present in the unfolded form ensemble, the effects of the salt KCl on the stabilities of the wild-type and the mutant proteins, as well as on the structures present in the unfolded form ensemble, were also studied. It was found that the m values of the m(0), m(+) and m(-) mutant proteins all converge to the wild-type m value in the presence of KCl. This result indicates that the salt modulates residual structure in the unfolded form by screening electrostatic interactions that maintain compact and expanded components in the unfolded protein ensemble. The use of free energy cycles has allowed the effect of salt on the structure and free energy of the unfolded protein to be related to the stability of the protein.     

Pulse proteolysis: a simple method for quantitative determination of protein stability and ligand binding

Thermodynamic stability is fundamental to the biology of proteins. Information on protein stability is essential for studying protein structure and folding and can also be used indirectly to monitor protein-ligand or protein-protein interactions. While clearly valuable, the experimental determination of a protein's stability typically requires biophysical instrumentation and substantial quantities of purified protein, which has limited the use of this technique as a general laboratory method. We report here a simple new method for determining protein stability by using pulse proteolysis with varying concentrations of denaturant. Pulse proteolysis is designed to digest only the unfolded proteins in an equilibrium mixture of folded and unfolded proteins that relaxes on a time scale longer than the proteolytic pulse. We used this method to study the stabilities of Escherichia coli ribonuclease H and its variants, both in purified form and directly from cell lysates. The DeltaG(unf) degrees values obtained by this technique were in agreement with those determined by traditional methods. We also successfully used this method to monitor the binding of maltose-binding protein to maltose, as well as to rapidly screen cognate ligands for this protein. The simplicity of pulse proteolysis suggests that it is an excellent strategy for the high-throughput determination of protein stability in protein engineering and drug discovery applications.     

Absence of a stable intermediate on the folding pathway of protein A.

The B-domain of protein A has one of the simplest protein topologies, a three-helix bundle. Its folding has been studied as a model for elementary steps in the folding of larger proteins. Earlier studies suggested that folding might occur by way of a helical hairpin intermediate. Equilibrium hydrogen exchange measurements indicate that the C-terminal helical hairpin could be a potential folding intermediate. Kinetic refolding experiments were performed using stopped-flow circular dichroism and NMR hydrogen-deuterium exchange pulse labeling. Folding of the entire molecule is essentially complete within the 6 ms dead time of the quench-flow apparatus, indicating that the intermediate, if formed, progresses rapidly to the final folded state. Site-directed mutagenesis of the isoleucine residue at position 16 was used to generate a variant protein containing tryptophan (the 116 W mutant). The formation of the putative folding intermediate was expected to be favored in this mutant at the expense of the native folded form, due to predicted unfavorable steric interactions of the bulky tryptophan side chain in the folded state. The 116 W mutant refolds completely within the dead time of a stopped-flow fluorescence experiment. No partly folded intermediate could be detected by either kinetic or equilibrium measurements. Studies of peptide fragments suggest that the protein A sequence has an intrinsic propensity to form a helix II/helix III hairpin. However, its stability appears to be marginal (of the order of 1/2 kT) and it could not be an obligatory intermediate on a defined folding pathway. These results explicitly demonstrate that the protein A B domain folds extremely rapidly by an apparent two-state mechanism without formation of stable partly folded intermediates. Similar mechanisms may also be involved in the rapid folding of subdomains of larger proteins to form the compact molten globule intermediates that often accumulate during the folding process.