Barriers to protein folding: formation of buried polar interactions is a slow step in acquisition of structure.

In the MYL mutant of the Arc repressor dimer, sets of partially buried salt-bridge and hydrogen-bond interactions mediated by Arg-31, Glu-36, and Arg-40 in each subunit are replaced by hydrophobic interactions between Met-31, Tyr-36, and Leu-40. The MYL refolding/dimerization reaction differs from that of wild type in being 10- to 1250-fold faster, having an earlier transition state, and depending upon viscosity but not ionic strength. Formation of the wild-type salt bridges in a hydrophobic environment clearly imposes a kinetic barrier to folding, which can be lowered by high salt concentrations. The changes in the position of the transition state and viscosity dependence can be explained if denatured monomers interact to form a partially folded dimeric intermediate, which then continues folding to form the native dimer. The second step is postulated to be rate limiting for wild type. Replacing the salt bridge with hydrophobic interactions lowers this barrier for MYL. This makes the first kinetic barrier rate limiting for MYL refolding and creates a downhill free-energy landscape in which most molecules which reach the intermediate state continue to form native dimers.

The protein-folding activity of chaperonins correlates with the symmetric GroEL14(GroES7)2 heterooligomer.

Chaperonins GroEL and GroES form, in the presence of ATP, two types of heterooligomers in solution: an asymmetric GroEL14GroES7 "bullet"-shaped particle and a symmetric GroEL14(GroES7)2 "football"-shaped particle. Under limiting concentrations of ATP or GroES, excess ADP, or in the presence of 5'-adenylyl imidodiphosphate, a correlation is seen between protein folding and the amount of symmetric GroEL14(GroES7)2 particles in a chaperonin solution, as detected by electron microscopy or by chemical crosslinking. Kinetic analysis suggests that protein folding is more efficient when carried out by a chaperonin solution populated with a majority of symmetric GroEL14(GroES7)2 particles than by a majority of asymmetric GroEL14GroES7 particles. The symmetric heterooligomer behaves as a highly efficient intermediate of the chaperonin protein folding cycle in vitro.

Optimization of rates of protein folding: the nucleation-condensation mechanism and its implications.

Small, single-module proteins that fold in a single cooperative step may be paradigms for understanding early events in protein-folding pathways generally. Recent experimental studies of the 64-residue chymotrypsin inhibitor 2 (CI2) support a nucleation mechanism for folding, as do some computer stimulations. CI2 has a nucleation site that develops only in the transition state for folding. The nucleus is composed of a set of adjacent residues (an alpha-helix), stabilized by long-range interactions that are formed as the rest of the protein collapses around it. A simple analysis of the optimization of the rate of protein folding predicts that rates are highest when the denatured state has little residual structure under physiological conditions and no intermediates accumulate. This implies that any potential nucleation site that is composed mainly of adjacent residues should be just weakly populated in the denatured state and become structured only in a high-energy intermediate or transition state when it is stabilized by interactions elsewhere in the protein. Hierarchical mechanisms of folding in which stable elements of structure accrete are unfavorable. The nucleation-condensation mechanism of CI2 fulfills the criteria for fast folding. On the other hand, stable intermediates do form in the folding of more complex proteins, and this may be an unavoidable consequence of increasing size and nucleation at more than one site.

Submillisecond events in protein folding.

The pathway of protein folding is now being analyzed at the resolution of individual residues by kinetic measurements on suitably engineered mutants. The kinetic methods generally employed for studying folding are typically limited to the time range of > or = 1 ms because the folding of denatured proteins is usually initiated by mixing them with buffers that favor folding, and the dead time of rapid mixing experiments is about a millisecond. We now show that the study of protein folding may be extended to the microsecond time region by using temperature-jump measurements on the cold-unfolded state of a suitable protein. We are able to detect early events in the folding of mutants of barstar, the polypeptide inhibitor of barnase. A preliminary characterization of the fast phase from spectroscopic and phi-value analysis indicates that it is a transition between two relatively solvent-exposed states with little consolidation of structure.

G protein activation kinetics and spillover of gamma-aminobutyric acid may account for differences between inhibitory responses in the hippocampus and thalamus.

We have developed a model of gamma-aminobutyric acid (GABA)ergic synaptic transmission mediated by GABAA and GABAB receptors, including cooperativity in the guanine nucleotide binding protein (G protein) cascade mediating the activation of K+ channels by GABAB receptors. If the binding of several G proteins is needed to activate the K+ channels, then only a prolonged activation of GABAB receptors evoked detectable currents. This could occur if strong stimuli evoked release in adjacent terminals and the spillover resulted in prolonged activation of the receptors, leading to inhibitory responses similar to those observed in hippocampal slices. The same model also reproduced thalamic GABAB responses to high-frequency bursts of stimuli. In this case, prolonged activation of the receptors was due to high-frequency release conditions. This model provides insights into the function of GABAB receptors in normal and epileptic discharges.

Prediction of protein folding class using global description of amino acid sequence.

We present a method for predicting protein folding class based on global protein chain description and a voting process. Selection of the best descriptors was achieved by a computer-simulated neural network trained on a data base consisting of 83 folding classes. Protein-chain descriptors include overall composition, transition, and distribution of amino acid attributes, such as relative hydrophobicity, predicted secondary structure, and predicted solvent exposure. Cross-validation testing was performed on 15 of the largest classes. The test shows that proteins were assigned to the correct class (correct positive prediction) with an average accuracy of 71.7%, whereas the inverse prediction of proteins as not belonging to a particular class (correct negative prediction) was 90-95% accurate. When tested on 254 structures used in this study, the top two predictions contained the correct class in 91% of the cases.

Cyclophilin catalyzes protein folding in yeast mitochondria.

Cyclophilins are a family of ubiquitous proteins that are the intracellular target of the immunosuppressant drug cyclosporin A. Although cyclophilins catalyze peptidylprolyl cis-trans isomerization in vitro, it has remained open whether they also perform this function in vivo. Here we show that Cpr3p, a cyclophilin in the matrix of yeast mitochondria, accelerates the refolding of a fusion protein that was synthesized in a reticulocyte lysate and imported into the matrix of isolated yeast mitochondria. The fusion protein consisted of the matrix-targeting sequence of subunit 9 of F1F0-ATPase fused to mouse dihydrofolate reductase. Refolding of the dihydrofolate reductase moiety in the matrix was monitored by acquisition of resistance to proteinase K. The rate of refolding was reduced by a factor of 2-6 by 2.5 microM cyclosporin A. This reduced rate of folding was also observed with mitochondria lacking Cpr3p. In these mitochondria, protein folding was insensitive to cyclosporin A. The rate of protein import was not affected by cyclosporin A or by deletion of Cpr3p.

The folding of GroEL-bound barnase as a model for chaperonin-mediated protein folding.

We have analyzed the pathway of folding of barnase bound to GroEL to resolve the controversy of whether proteins can fold while bound to chaperonins (GroEL or Cpn60) or fold only after their release into solution. Four phases in the folding were detected by rapid-reaction kinetic measurements of the intrinsic fluorescence of both wild type and barnase mutants. The phases were assigned from their rate laws, sensitivity to mutations, and correspondence to regain of catalytic activity. At high ratios of denatured barnase to GroEL, 4 mol of barnase rapidly bind per 14-mer of GroEL. At high ratios of GroEL to barnase, 1 mol of barnase binds with a rate constant of 3.5 x 10(7) s-1.M-1. This molecule then refolds with a low rate constant that changes on mutation in parallel with the rate constant for the folding in solution. This rate constant corresponds to the regain of the overall catalytic activity of barnase and increases 15-fold on the addition of ATP to a physiologically relevant value of approximately 0.4 s-1. The multiply bound molecules of barnase that are present at high ratios of GroEL to barnase fold with a rate constant that is also sensitive to mutation but is 10 times higher. If the 110-residue barnase can fold when bound to GroEL and many moles can bind simultaneously, then smaller parts of large proteins should be able to fold while bound.

Toward an outline of the topography of a realistic protein-folding funnel.

Experimental information on the structure and dynamics of molten globules gives estimates for the energy landscape's characteristics for folding highly helical proteins, when supplemented by a theory of the helix-coil transition in collapsed heteropolymers. A law of corresponding states relating simulations on small lattice models to real proteins possessing many more degrees of freedom results. This correspondence reveals parallels between "minimalist" lattice results and recent experimental results for the degree of native character of the folding transition state and molten globule and also pinpoints the needs of further experiments.

A specific transition state for S-peptide combining with folded S-protein and then refolding

We measured the folding and unfolding kinetics of mutants for a simple protein folding reaction to characterize the structure of the transition state. Fluorescently labeled S-peptide analogues combine with S-protein to form ribonuclease S analogues: initially, S-peptide is disordered whereas S-protein is folded. The fluorescent probe provides a convenient spectroscopic probe for the reaction. The association rate constant, kon, and the dissociation rate constant, koff, were both determined for two sets of mutants. The dissociation rate constant is measured by adding an excess of unlabeled S-peptide analogue to a labeled complex (RNaseS*). This strategy allows kon and koff to be measured under identical conditions so that microscopic reversibility applies and the transition state is the same for unfolding and refolding. The first set of mutants tests the role of the α-helix in the transition state. Solvent-exposed residues Ala-6 and Gln-11 in the α-helix of native RNaseS were replaced by the helix destabilizing residues glycine or proline. A plot of log kon vs. log Kd for this series of mutants is linear over a very wide range, with a slope of −0.3, indicating that almost all of the molecules fold via a transition state involving the helix. A second set of mutants tests the role of side chains in the transition state. Three side chains were investigated: Phe-8, His-12, and Met-13, which are known to be important for binding S-peptide to S-protein and which also contribute strongly to the stability of RNaseS*. Only the side chain of Phe-8 contributes significantly, however, to the stability of the transition state. The results provide a remarkably clear description of a folding transition state.

Magnesium-dependent folding of self-splicing RNA: Exploring the link between cooperativity, thermodynamics, and kinetics

Folding of the Tetrahymena self-splicing RNA into its active conformation involves a set of discrete intermediate states. The Mg2+-dependent equilibrium transition from the intermediates to the native structure is more cooperative than the formation of the intermediates from the unfolded states. We show that the degree of cooperativity is linked to the free energy of each transition and that the rate of the slow transition from the intermediates to the native state decreases exponentially with increasing Mg2+ concentration. Monovalent salts, which stabilize the folded RNA nonspecifically, induce states that fold in less than 30 s after Mg2+ is added to the RNA. A simple model is proposed that predicts the folding kinetics from the Mg2+-dependent change in the relative stabilities of the intermediate and native states.

A statistical mechanical model for β-hairpin kinetics

Understanding the mechanism of protein secondary structure formation is an essential part of the protein-folding puzzle. Here, we describe a simple statistical mechanical model for the formation of a β-hairpin, the minimal structural element of the antiparallel β-pleated sheet. The model accurately describes the thermodynamic and kinetic behavior of a 16-residue, β-hairpin-forming peptide, successfully explaining its two-state behavior and apparent negative activation energy for folding. The model classifies structures according to their backbone conformation, defined by 15 pairs of dihedral angles, and is further simplified by considering only the 120 structures with contiguous stretches of native pairs of backbone dihedral angles. This single sequence approximation is tested by comparison with a more complete model that includes the 215 possible conformations and 15 × 215 possible kinetic transitions. Finally, we use the model to predict the equilibrium unfolding curves and kinetics for several variants of the β-hairpin peptide.

Folding funnels and frustration in off-lattice minimalist protein landscapes

A full quantitative understanding of the protein folding problem is now becoming possible with the help of the energy landscape theory and the protein folding funnel concept. Good folding sequences have a landscape that resembles a rough funnel where the energy bias towards the native state is larger than its ruggedness. Such a landscape leads not only to fast folding and stable native conformations but, more importantly, to sequences that are robust to variations in the protein environment and to sequence mutations. In this paper, an off-lattice model of sequences that fold into a β-barrel native structure is used to describe a framework that can quantitatively distinguish good and bad folders. The two sequences analyzed have the same native structure, but one of them is minimally frustrated whereas the other one exhibits a high degree of frustration.

Cytochrome c′ folding triggered by electron transfer: Fast and slow formation of four-helix bundles

Reduced (FeII) Rhodopseudomonas palustris cytochrome c′ (Cyt c′) is more stable toward unfolding ([GuHCl]1/2 = 2.9(1) M) than the oxidized (FeIII) protein ([GuHCl]1/2 = 1.9(1) M). The difference in folding free energies (ΔΔGf° = 70 meV) is less than half of the difference in reduction potentials of the folded protein (100 mV vs. NHE) and a free heme in aqueous solution (≈−150 mV). The spectroscopic features of unfolded FeII–Cyt c′ indicate a low-spin heme that is axially coordinated to methionine sulfur (Met-15 or Met-25). Time-resolved absorption measurements after CO photodissociation from unfolded FeII(CO)–Cyt c′ confirm that methionine can bind to the ferroheme on the microsecond time scale [kobs = 5(2) × 104 s−1]. Protein folding was initiated by photoreduction (two-photon laser excitation of NADH) of unfolded FeIII–Cyt c′ ([GuHCl] = 2.02–2.54 M). Folding kinetics monitored by heme absorption span a wide time range and are highly heterogeneous; there are fast-folding (≈103 s−1), intermediate-folding (102–101 s−1), and slow-folding (10−1 s−1) populations, with the last two likely containing methionine-ligated (Met-15 or Met-25) ferrohemes. Kinetics after photoreduction of unfolded FeIII–Cyt c′ in the presence of CO are attributable to CO binding [1.4(6) × 103 s−1] and FeII(CO)–Cyt c′ folding [2.8(9) s−1] processes; stopped-flow triggered folding of FeIII–Cyt c′ (which does not contain a protein-derived sixth ligand) is adequately described by a single kinetics phase with an estimated folding time constant of ≈4 ms [ΔGf° = −33(3) kJ mol−1] at zero denaturant.

Configurational diffusion down a folding funnel describes the dynamics of DNA hairpins

Elucidating the mechanism of folding of polynucleotides depends on accurate estimates of free energy surfaces and a quantitative description of the kinetics of structure formation. Here, the kinetics of hairpin formation in single-stranded DNA are measured after a laser temperature jump. The kinetics are modeled as configurational diffusion on a free energy surface obtained from a statistical mechanical description of equilibrium melting profiles. The effective diffusion coefficient is found to be strongly temperature-dependent in the nucleation step as a result of formation of misfolded loops that do not lead to subsequent zipping. This simple system exhibits many of the features predicted from theoretical studies of protein folding, including a funnel-like energy surface with many folding pathways, trapping in misfolded conformations, and non-Arrhenius folding rates.