Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

Ultrafast signals in protein folding and the polypeptide contracted state

To test the significance of ultrafast protein folding signals (≪1 msec), we studied cytochrome c (Cyt c) and two Cyt c fragments with major C-terminal segments deleted. The fragments remain unfolded under all conditions and so could be used to define the unfolded baselines for protein fluorescence and circular dichroism (CD) as a function of denaturant concentration. When diluted from high to low denaturant in kinetic folding experiments, the fragments readjust to their new baseline values in a “burst phase” within the mixing dead time. The fragment burst phase reflects a contraction of the polypeptide from a more extended unfolded condition at high denaturant to a more contracted unfolded condition in the poorer, low denaturant solvent. Holo Cyt c exhibits fluorescence and CD burst phase signals that are essentially identical to the fragment signals over the whole range of final denaturant concentrations, evidently reflecting the same solvent-dependent, relatively nonspecific contraction and not the formation of a specific folding intermediate. The significance of fast folding signals in Cyt c and other proteins is discussed in relation to the hypothesis of an initial rate-limiting search-nucleation-collapse step in protein folding [Sosnick, T. R., Mayne, L. & Englander, S. W. (1996) Proteins Struct. Funct. Genet. 24, 413–426].

Linking topography of its potential surface with the dynamics of folding of a protein model

The “3-color, 46-bead” model of a folding polypeptide is the vehicle for adapting to proteins a mode of analysis used heretofore for atomic clusters, to relate the topography of the potential surface to the dynamics that lead to formation of selected structures. The analysis is based on sequences of stationary points—successive minima, joined by saddles—that rise monotonically in energy from basin bottoms. Like structure-seeking clusters, the potential surface of the model studied here is staircase-like, rather than sawtooth-like, with highly collective motions required for passage from one minimum to the next. The surface has several deep basins whose minima correspond to very similar structures, but which are separated by high energy barriers.

RNA folding kinetics regulates translation of phage MS2 maturation gene

The gene for the maturation protein of the single-stranded RNA coliphage MS2 is preceded by an untranslated leader of 130 nt, which folds into a cloverleaf, i.e., three stem–loop structures enclosed by a long distance interaction (LDI). This LDI prevents translation because its 3′ moiety contains the Shine–Dalgarno sequence of the maturation gene. Previously, several observations suggested that folding of the cloverleaf is kinetically delayed, providing a time window for ribosomes to access the RNA. Here we present direct evidence for this model. In vitro experiments show that ribosome binding to the maturation gene is faster than refolding of the denatured cloverleaf. This folding delay appears related to special properties of the leader sequence. We have replaced the three stem–loop structures by a single five nt loop. This change does not affect the equilibrium structure of the LDI. Nevertheless, in this construct, the folding delay has virtually disappeared, suggesting that now the RNA folds faster than ribosomes can bind. Perturbation of the cloverleaf by an insertion makes the maturation start permanently accessible. A pseudorevertant that evolved from an infectious clone carrying the insertion had overcome this defect. It showed a wild-type folding delay before closing down the maturation gene. This experiment reveals the biological significance of retarded cloverleaf formation.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Structure and function in rhodopsin: Packing of the helices in the transmembrane domain and folding to a tertiary structure in the intradiscal domain are coupled*

A previous study of the retinitis pigmentosa mutation L125R and two designed mutations at this site, L125A and L125F, showed that these mutations cause partial or total misfolding of the opsins expressed in COS cells from the corresponding mutant opsin genes. We now report on expression and characterization of the opsins from the following retinitis pigmentosa mutants in the transmembrane domain of rhodopsin that correspond to six of the seven helices: G51A and G51V (helix A), G89D (helix B), A164V (helix D), H211P (helix E), P267L and P267R (helix F), and T297R (helix G). All the mutations caused partial misfolding of the opsins as observed by the UV/visible absorption characteristics and by separation of the expressed opsins into fractions that bound 11-cis-retinal to form the corresponding mutant rhodopsins and those that did not bind 11-cis-retinal. Further, all the mutant rhodopsins prepared from the above mutants, except for G51A, showed strikingly abnormal bleaching behavior with abnormal metarhodopsin II photointermediates. The results show that retinitis pigmentosa mutations in every one of the transmembrane helices can cause misfolding of the opsin. Therefore, on the basis of these and previous results, we conclude that defects in the packing of the transmembrane helices resulting from these mutations are relayed to the intradiscal domain, where they cause misfolding of the opsin by inducing the formation of a disulfide bond other than the native Cys-110—Cys-187 disulfide bond. Thus, there is coupling between packing of the helices in the transmembrane domain and folding to a tertiary structure in the intradiscal domain.

Protein folding kinetics exhibit an Arrhenius temperature dependence when corrected for the temperature dependence of protein stability

The anomalous temperature dependence of protein folding has received considerable attention. Here we show that the temperature dependence of the folding of protein L becomes extremely simple when the effects of temperature on protein stability are corrected for; the logarithm of the folding rate is a linear function of 1/T on constant stability contours in the temperature–denaturant plane. This convincingly demonstrates that the anomalous temperature dependence of folding derives from the temperature dependence of the interactions that stabilize proteins, rather than from the super Arrhenius temperature dependence predicted for the configurational diffusion constant on a rough energy landscape. However, because of the limited temperature range accessible to experiment, the results do not rule out models with higher order temperature dependences. The significance of the slope of the stability-corrected Arrhenius plots is discussed.

Characteristic temperatures of folding of a small peptide

We perform a generalized-ensemble simulation of a small peptide taking the interactions among all atoms into account. From this simulation we obtain thermodynamic quantities over a wide range of temperatures. In particular, we show that the folding of a small peptide is a multistage process associated with two characteristic temperatures, the collapse temperature Tθ and the folding temperature Tƒ. Our results give supporting evidence for the energy landscape picture and funnel concept. These ideas were previously developed in the context of studies of simplified protein models, and here are checked in an all-atom Monte Carlo simulation.

Evolution of the folding ability of proteins through functional selection

An evolutionary process is simulated with a simple spin-glass-like model of proteins to examine the origin of folding ability. At each generation, sequences are randomly mutated and subjected to a simulation of the folding process based on the model. According to the frequency of local configurations at the active sites, sequences are selected and passed to the next generation. After a few hundred generations, a sequence capable of folding globally into a native conformation emerges. Moreover, the selected sequence has a distinct energy minimum and an anisotropic funnel on the energy surface, which are the imperative features for fast folding of proteins. The proposed model reveals that the functional selection on the local configurations leads a sequence to fold globally into a conformation at a faster rate.

Folding in vivo of a newly translated yeast cytosolic enzyme is mediated by the SSA class of cytosolic yeast Hsp70 proteins

The nature of chaperone action in the eukaryotic cytosol that assists newly translated cytosolic proteins to reach the native state has remained poorly defined. Actin, tubulin, and Gα transducin are assisted by the cytosolic chaperonin, CCT, but many other proteins, for example, ornithine transcarbamoylase (OTC), a cytosolic homotrimeric enzyme of yeast, do not require CCT action. Here, we observe that yeast cytosolic OTC is assisted to its native state by the SSA class of yeast cytosolic Hsp70 proteins. In vitro, refolding of OTC diluted from denaturant was assisted by crude yeast cytosol and ATP and found to be directed by SSA1/2. In vivo, when OTC was induced in a temperature-sensitive SSA-deficient strain, it exhibited reduced specific activity, and nonnative subunits were detected in the soluble fraction. These findings indicate that, in vivo, the Hsp70 system assists in folding at least some newly translated cytosolic enzymes, most likely functioning in a posttranslational manner.

Folding and aggregation of designed proteins

Protein aggregation is studied by following the simultaneous folding of two designed identical 20-letter amino acid chains within the framework of a lattice model and using Monte Carlo simulations. It is found that protein aggregation is determined by elementary structures (partially folded intermediates) controlled by local contacts among some of the most strongly interacting amino acids and formed at an early stage in the folding process.

Detection of residue contacts in a protein folding intermediate

Protein folding can be described in terms of the development of specific contacts between residues as a highly disordered polypeptide chain converts into the native state. Here we describe an NMR based strategy designed to detect such contacts by observation of nuclear Overhauser effects (NOEs). Experiments with α-lactalbumin reveal the existence of extensive NOEs between aromatic and aliphatic protons in the archetypal molten globule formed by this protein at low pH. Analysis of their time development provides direct evidence for near-native compactness of this state. Through a rapid refolding procedure the NOE intensity can be transferred efficiently into the resolved and assigned spectrum of the native state. This demonstrates the viability of using this approach to map out time-averaged interactions between residues in a partially folded protein.

From residue matching patterns to protein folding topographies: General model and bovine pancreatic trypsin inhibitor

A coarse-grained model for protein-folding dynamics is introduced based on a discretized representation of torsional modes. The model, based on the Ramachandran map of the local torsional potential surface and the class (hydrophobic/polar/neutral) of each residue, recognizes patterns of both torsional conformations and hydrophobic-polar contacts, with tolerance for imperfect patterns. It incorporates empirical rates for formation of secondary and tertiary structure. The method yields a topological representation of the evolving local torsional configuration of the folding protein, modulo the basins of the Ramachandran map. The folding process is modeled as a sequence of transitions from one contact pattern to another, as the torsional patterns evolve. We test the model by applying it to the folding process of bovine pancreatic trypsin inhibitor, obtaining a kinetic description of the transitions between the contact patterns visited by the protein along the dominant folding pathway. The kinetics and detailed balance make it possible to invert the result to obtain a coarse topographic description of the potential energy surface along the dominant folding pathway, in effect to go backward or forward between a topological representation of the chain conformation and a topographical description of the potential energy surface governing the folding process. As a result, the strong structure-seeking character of bovine pancreatic trypsin inhibitor and the principal features of its folding pathway are reproduced in a reasonably quantitative way.

Conformational influences of glycosylation of a peptide: A possible model for the effect of glycosylation on the rate of protein folding

Improved strategies for synthesis make it possible to expand the range of glycopeptides available for detailed conformational studies. The glycopeptide 1 was synthesized using a new solid phase synthesis of carbohydrates and a convergent coupling to peptide followed by deprotection. Its conformational properties were subjected to NMR analysis and compared with a control peptide 2 prepared by conventional solid phase methods. Whereas peptide 2 fails to manifest any appreciable secondary structure, the glycopeptide 1 does show considerable conformational bias suggestive of an equilibrium between an ordered and a random state. The implications of this ordering effect for the larger issue of protein folding are considered.

Efficient folding of proteins with multiple disulfide bonds in the Escherichia coli cytoplasm

Under physiological conditions, the Escherichia coli cytoplasm is maintained in a reduced state that strongly disfavors the formation of stable disulfide bonds in proteins. However, mutants in which the reduction of both thioredoxins and glutathione is impaired (trxB gor mutants) accumulate oxidized, enzymatically active alkaline phosphatase in the cytoplasm. These mutants grow very poorly in the absence of an exogenous reductant and accumulate extragenic suppressors at a high frequency. One such suppressor strain, FA113, grows almost as rapidly as the wild type in the absence of reductant, exhibits slightly faster kinetics of disulfide bond formation, and has fully induced activity of the transcriptional activator, OxyR. FA113 gave substantially higher yields of properly oxidized proteins compared with wild-type or trxB mutant strains. For polypeptides with very complex patterns of disulfide bonds, such as vtPA and the full-length tPA, the amount of active protein was further enhanced up to 15-fold by co-expression of TrxA (thioredoxin 1) mutants with different redox potentials, or 20-fold by the protein disulfide isomerase, DsbC. Remarkably, higher yields of oxidized, biologically active proteins were obtained by expression in the cytoplasm of E. coli FA113 compared with what could be achieved via secretion into the periplasm of a wild-type strain, even under optimized conditions. These results demonstrate that the cytoplasm can be rendered sufficiently oxidizing to allow efficient formation of native disulfide bonds without compromising cell viability.