Energetics of protein folding

The energetics of protein folding determine the 3D structure of a folded protein. Knowledge of the energetics is needed to predict the 3D structure from the amino acid sequence or to modify the structure by protein engineering. Recent developments are discussed: major factors are reviewed and auxiliary factors are discussed briefly. Major factors include the hydrophobic factor (burial of non-polar surface area) and van der Waals interactions together with peptide hydrogen bonds and peptide solvation. The long-standing model for the hydrophobic factor (free energy change proportional to buried non-polar surface area) is contrasted with the packing-desolvation model and the approximate nature of the proportionality between free energy and apolar surface area is discussed. Recent energetic studies of forming peptide hydrogen bonds (gas phase) are reviewed together with studies of peptide solvation in solution. Closer agreement is achieved between the 1995 values for protein unfolding enthalpies in vacuum given by Lazaridis-Archontis-Karplus and Makhatadze-Privalov when the solvation enthalpy of the peptide group is taken from electrostatic calculations. Auxiliary factors in folding energetics include salt bridges and side-chain hydrogen bonds, disulfide bridges, and propensities to form alpha-helices and beta-structure. Backbone conformational entropy is a major energetic factor which is discussed only briefly for lack of knowledge.     

Energetics of membrane protein folding and stability

The critical role of membrane proteins in a myriad of biological and physiological functions has spawned numerous investigations over the past several decades with the long-term goal of identifying the molecular origins and energetic forces that stabilize these proteins within the membrane. Parallel structural and thermodynamics studies on several systems have provided significant insight regarding the driving forces governing folding, assembly, insertion, and translocation of membrane proteins. The present review surveys families of membrane-associated proteins including α-helical and β-barrel structures, viral surface receptors, and pore-forming toxins, citing representative proteins within each of these classes for further scrutiny in terms of structure-function relationships and global conformational stability. This overview presents seminal findings from pioneering studies on the energetics of membrane protein folding and stability to modern techniques that are exploiting the use of molecular genetics and single molecule studies. An overall consensus regarding the molecular origins of membrane protein stability is that a number of intrinsic properties resemble features of soluble proteins, yet there are distinct energetic differences arising from specific intra- and intermolecular interactions within the membrane. The combined efforts from structural, energetics, and dynamics approaches offer unique insights and improve our fundamental understanding of the driving forces dictating membrane protein folding and stability.     

Energetics of multihelix interactions in protein folding: application to myoglobin

Conformational-energy calculations have been carried out in order to determine favorable packing arrangements within a group of α-helices. The influence of side chains and of the number of interacting α-helices on the mode of packing was analyzed. In this work, our earlier methods for computing the packing energy of a pair of α-helices [Chou, K.-C., Némethy, G. & Scheraga, H. A. (1984) J. Am. Chem. Soc. 106 , 3161–3170] have been extended to treat the interactions among several helices. Also, new algorithms allow the matching of standard peptide geometry to x-ray coordinates of helical complexes and the analysis of interrelations between several helices. As a specific test case, the packing of three neighboring α-helices, viz., the A, G, and H helices of sperm whale myoglobin, was considered. Minimum-energy arrangements were computed for the separate A-H and the G-H α-helix pairs as well as for the A-G-H three-helix complex. For the packing of the nearly antiparallel G and H α-helices, the same optimal structure was obtained in two- and three-helix complexes, indicating that a single packing arrangement is specifically favored by interhelix interactions. For the pair of nearly perpendicular A and H α-helices, interactions are less specific, so that there is no unique optimal structure in the two-helix complex; in the three-helix complex, however, a specific mode of packing is favored even for the A-H pair. This result indicates that the presence of other nearby α-helices can influence the packing of a given α-helix pair. The computed arrangement of the A-G-H complex is very close to that of the crystallographically determined structure. These results can be used to make deductions about the likely sequence of events in protein folding, where, in this particular case, it appears that the G-H helix pair may form first and then induce proper orientation of the A helix.     

Physical principles of protein structure and protein folding

Physical principles determining the protein structure and protein folding are reviewed: (i) the molecular theory of protein secondary structure and the method of its prediction based on this theory; (ii) the existence of a limited set of thermodynamically favourable folding patterns of α- and β-regions in a compact globule which does not depend on the details of the amino acid sequence; (iii) the moderns approaches to the prediction of the folding patterns of α- and β-regions in concrete proteins; (iv) experimental approaches to the mechanism of protein folding. The review reflects theoretical and experimental works of the author and his collaborators as well as those of other groups.     

Energetics of folding subtilisin BPN'.

Subtilisin is an unusual example of a monomeric protein with a substantial kinetic barrier to folding and unfolding. Here we document for the first time the in vitro folding of the mature form of subtilisin. Subtilisin was modified by site-directed mutagenesis to be proteolytically inactive, allowing the impediments to folding to be systematically examined. First, the thermodynamics and kinetics of calcium binding to the high-affinity calcium A-site have been measured by microcalorimetry and fluorescence spectroscopy. Binding is an enthalpically driven process with an association constant (Ka) equal to 7 x 10(6) M-1. Furthermore, the kinetic barrier to calcium removal from the A-site (23 kcal/mol) is substantially larger than the standard free energy of binding (9.3 kcal/mol). The kinetics of calcium dissociation from subtilisin (e.g., in excess EDTA) are accordingly very slow (t1/2 = 1.3 h at 25 degrees C). Second, to measure the kinetics of subtilisin folding independent of calcium binding, the high-affinity calcium binding site was deleted from the protein. At low ionic strength (I = 0.01) refolding of this mutant requires several days. The folding rate is accelerated almost 100-fold by a 10-fold increase in ionic strength, indicating that part of the free energy of activation may be electrostatic. At relatively high ionic strength (I = 0.5) refolding of the mutant subtilisin is complete in less than 1 h at 25 degrees C. We suggest that part of the electrostatic contribution to the activation free energy for folding subtilisin is related to the highly charged region of the protein comprising the weak ion binding site (site B).(ABSTRACT TRUNCATED AT 250 WORDS)     

Evolutionary aspects of protein structure and folding

Four basic stages of evolution of protein structure are described based on recent work of the authors targeted specifically on reconstruction of the earliest events in the protein evolution. According to this reconstruction, the initial stage of short peptides of, probably, only few amino-acid residues had been followed by formation of closed loops of the size 25-30 residues, which corresponds to the polymer-statistically optimal ring closure size for mixed polypeptide chains. The next stage involved fusion of the respective small linear genes and formation of protein structures consisting of several closed loops of the nearly standard size, up to 4-6 loops (100-200 amino acid residues) in a typical protein fold. The last, modern stage began with combinatorial fusion of the presumably circular 300-600 bp DNA units and, accordingly, formation of multidomain proteins.     

Evolutionary aspects of protein structure and folding

The traditional reconstruction of molecular events of the past based on sequence conservation becomes very vague beyond one to two billion years ago. There are certain molecular features, however, such as polymer flexibility and loop closure, that are conserved merely because of their physical nature. This allows one to penetrate the earliest stages of protein evolution.     

Crystal structure of a defective folding protein

Maltose-binding protein (MBP or MalE) of Escherichia coli is the periplasmic receptor of the maltose transport system. MalE31, a defective folding mutant of MalE carrying sequence changes Gly 32-->Asp and Ile 33-->Pro, is either degraded or forms inclusion bodies following its export to the periplasmic compartment. We have shown previously that overexpression of FkpA, a heat-shock periplasmic peptidyl-prolyl isomerase with chaperone activity, suppresses MalE31 misfolding. Here, we have exploited this property to characterize the maltose transport activity of MalE31 in whole cells. MalE31 displays defective transport behavior, even though it retains maltose-binding activity comparable with that of the wild-type protein. Because the mutated residues are in a region on the surface of MalE not identified previously as important for maltose transport, we have solved the crystal structure of MalE31 in the maltose-bound state in order to characterize the effects of these changes. The structure was determined by molecular replacement methods and refined to 1.85 A resolution. The conformation of MalE31 closely resembles that of wild-type MalE, with very small displacements of the mutated residues located in the loop connecting the first alpha-helix to the first beta-strand. The structural and functional characterization provides experimental evidence that MalE31 can attain a wild-type folded conformation, and suggest that the mutated sites are probably involved in the interactions with the membrane components of the maltose transport system.     

Energetics and dynamics of SNAREpin folding across lipid bilayers

Membrane fusion occurs when SNAREpins fold up between lipid bilayers. How much energy is generated during SNAREpin folding and how this energy is coupled to the fusion of apposing membranes is unknown. We have used a surface forces apparatus to determine the energetics and dynamics of SNAREpin formation and characterize the different intermediate structures sampled by cognate SNAREs in the course of their assembly. The interaction energy-versus-distance profiles of assembling SNAREpins reveal that SNARE motifs begin to interact when the membranes are 8 nm apart. Even after very close approach of the bilayers (approximately 2-4 nm), the SNAREpins remain partly unstructured in their membrane-proximal region. The energy stabilizing a single SNAREpin in this configuration (35 k(B)T) corresponds closely with the energy needed to fuse outer but not inner leaflets (hemifusion) of pure lipid bilayers (40-50 k(B)T).     

The relation between mRNA folding and protein structure

About 200 mRNA sequences of Escherichia coli and human with matching protein secondary structure data were studied. The mRNA folding for each native sequence and for corresponding randomized sequences was calculated through free energy minimization. We have found that the folding energy of mRNA segments in different protein secondary structures is significantly different. The average Z score is more negative for regular secondary structure (alpha-helix and beta-strand) than that for coil. This suggests that the codon choice in native mRNA sequence coding for protein regular structure contributes more to the mRNA folding stability.     

Perfect temperature for protein structure prediction and folding

We have investigated the influence of the “noise” of inevitable errors in energetic parameters on-protein structure prediction. Because of this noise, only a part of all the interactions operating in a protein chain can be taken into account, and therefore a search for the energy minimum becomes inadequate for protein structure prediction. One can rather rely on statistical mechanics: a calculation carried out at a temperature T_* somewhat below that of protein melting gives the best possible, though always approximate prediction. The early stages of protein folding also “take into account” only a part of all the interactions; consequently, the same temperature T_* is favorable for the self-organization of native-like intermediates in protein folding. © 1995 Wiley-Liss, Inc.     

Energetics of folding of a lysozyme beta-bend.

The forces directing the “β-fold” at residues 52–59 in hen egg-white lysozyme have been explored by theoretical conformational analysis, which includes solvent interaction. It is shown that, whereas the conformation is in its most favorable free-energy state for a folded form, the fold is actually a destabilizing influence which is overcome only by long range interactions. The concept is introduced that nucleation of the tertiary structure initiates the folding process which is localized by the specific sequence. Thus, long range forces “drive” the fold and short range forces “localize” it.     

Evolutionary computer programming of protein folding and structure predictions

In order to understand the mechanism of protein folding and to assist the rational de-novo design of fast-folding, non-aggregating and stable artificial enzymes it is very helpful to be able to simulate protein folding reactions and to predict the structures of proteins and other biomacromolecules. Here, we use a method of computer programming called "evolutionary computer programming" in which a program evolves depending on the evolutionary pressure exerted on the program. In the case of the presented application of this method on a computer program for folding simulations, the evolutionary pressure exerted was towards faster finding deep minima in the energy landscape of protein folding. Already after 20 evolution steps, the evolved program was able to find deep minima in the energy landscape more than 10 times faster than the original program prior to the evolution process.     

Improved photo-CIDNP methods for studying protein structure and folding

Two new techniques offering considerable improvements in the quality of ¹H photo-CIDNP spectra of proteins are demonstrated. Both focus on the problem of progressive photo-degradation of the flavin dye used to generate polarization in exposed tryptophan, tyrosine and histidine side-chains. One approach uses rapid addition and removal of protein/flavin solution between light flashes to mix the NMR sample and introduce fresh dye into the laser-irradiated region. The other involves chemical oxidation of photo-reduced flavin by the addition of hydrogen peroxide. In both cases a larger number of scans can be accumulated before the flavin is exhausted than would otherwise be possible. The techniques are demonstrated by 600 MHz CIDNP-NOESY spectroscopy of bovine holo--lactalbumin, and by real-time CIDNP observation of the refolding of bovine apo--lactalbumin following rapid dilution from a high concentration of chemical denaturant.     

Dynamics, energetics, and structure in protein folding

For many decades, protein folding experimentalists have worked with no information about the time scales of relevant protein folding motions and without methods for estimating the height of folding barriers. Protein folding experiments have been interpreted using chemical models in which the folding process is characterized as a series of equilibria between two or more distinct states that interconvert with activated kinetics. Accordingly, the information to be extracted from experiments was circumscribed to apparent equilibrium constants and relative folding rates. Recent developments are changing this situation dramatically. The combination of fast-folding experiments with the development of analytical methods more closely connected to physical theory reveals that folding barriers in native conditions range from minimally high (approximately 14RT for the very slow folder AcP) to nonexistent. While slow-folding (i.e., > or = 1 ms) single-domain proteins are expected to fold in a two-state fashion, microsecond-folding proteins should exhibit complex behavior arising from crossing marginal or negligible folding barriers. This realization opens a realm of exciting opportunities for experimentalists. The free energy surface of a protein with a marginal (or no) barrier can be mapped using equilibrium experiments, which could resolve energetic factors from structural factors in folding. Kinetic experiments on these proteins provide the unique opportunity to measure folding dynamics directly. Furthermore, the complex distributions of time-dependent folding behaviors expected for these proteins might be accessible to single-molecule measurements. Here, we discuss some of these recent developments in protein folding, emphasizing aspects that can serve as a guide for experimentalists interested in exploiting this new avenue of research.     

Computer simulations of membrane protein folding: structure and dynamics

A lattice model of membrane proteins with a composite energy function is proposed to study their folding dynamics and native structures using Monte Carlo simulations. This model successfully predicts the seven helix bundle structure of sensory rhodopsin I by practicing a three-stage folding. Folding dynamics of a transmembrane segment into a helix is further investigated by varying the cooperativity in the formation of alpha helices for both random folding and assisted folding. The chain length dependence of the folding time of a hydrophobic segment to a helical state is studied for both free and anchored chains. An unusual length dependence in the folding time of anchored chains is observed.     

The consistency principle in protein structure and pathways of folding

• 1.i) It is pointed out that various energy terms contributing to stabilize the native state of globular proteins are consistent in the first approximation with each other in the native state. This means that each energy term is individually minimized at the minimum point of the total energy. I proposed (1) to call this fact “the consistency principle in protein structure.”
• 2.ii) The fair success of various methods of prediction of the secondary structures in globular proteins from their amino acid sequence is often interpreted as indicating the dominance of the short-range interactions in determining the local structures of the polypeptide chains. Partly from such a point of view, the hierarchic condensation model has been popular for the process of protein folding. However the consistency principle indicates that the short-range interactions are just one type of intramolecular interaction which contributes to stabilization of the native structure together with other mutually consistent types of intramolecular interactions. Therefore the hierarchic condensation model is not necessarily a unique model of protein folding.
• 3.iii) Roles of a possible nonspecific globular state, stabilized by nonspecific long-range intramolecular interactions, in the folding process are discussed. It is expected that this nonspecific globular state is observed either as an equilibrium or a kinetic intermediate state between the unfolded and the folded native states. Observation as a kinetic intermediate state is expected to occur in experiments done under strongly refolding conditions. In this case the polypeptide chain in the unfolded state collapses into a nonspecific globule by the action of nonspecific long-range intramolecular interactions. Two possible mechanisms of the transition from the nonspecific globular state to the specific native folded state are discussed.
• 4.iv) In an experiment done under weakly refolding conditions, folding is expected to occur according to the embryo-nucleus model. This model is a refined version of the hierarchic condensation model. Refinement is done by taking into account the fact that the intermediate structures assumed in the hierarchic condensation model are unstable against both the native folded state and the unfolded state. A nucleus is an ordered structure of a certain size. Ordered structures of a size larger than a nucleus tend to fold further to become the native specific globule. Ordered structures of a size smaller than a nucleus tend to unfold. Embryos are intrinsically unstable ordered structures smaller than a nucleus. Folding occurs when embryos grow in size to become a nucleus. The intrinsic instability of embryos is the built-in mechanism to overcome the low resolving power of the short-range interactions in determining local conformations of the polypeptide chain.