Conversion of a c type cytochrome to a b type that spontaneously forms in vitro from apo protein and heme: Implications for c type cytochrome biogenesis and folding

Cytochrome c552 from Hydrogenobacter thermophilus, a thermophilic bacterium, has been converted into a b type cytochrome, after mutagenesis of both heme-binding cysteines to alanine and expression in the cytoplasm of Escherichia coli. The b type variant is less stable, with the guanidine hydrochloride unfolding midpoint occurring at a concentration 2 M lower than for the wild-type protein. The reduction potential is 75 mV lower than that of the recombinant wild-type protein. The heme can be removed from the b type variant, thus generating an apo protein that has, according to circular dichroism spectroscopy, an α-helical content different from that of the holo b type protein. The latter is readily reformed in vitro by addition of heme to the apo protein. This reforming suggests that previously observed assembly of cytochrome c552, which has the typical class I cytochrome c fold, in the E. coli cytoplasm is a consequence of spontaneous thioether bond formation after binding of heme to a prefolded polypeptide. These observations have implications for the general problem of c type cytochrome biogenesis.

Single-molecule protein folding: Diffusion fluorescence resonance energy transfer studies of the denaturation of chymotrypsin inhibitor 2

We report single-molecule folding studies of a small, single-domain protein, chymotrypsin inhibitor 2 (CI2). CI2 is an excellent model system for protein folding studies and has been extensively studied, both experimentally (at the ensemble level) and theoretically. Conformationally assisted ligation methodology was used to synthesize the proteins and site-specifically label them with donor and acceptor dyes. Folded and denatured subpopulations were observed by fluorescence resonance energy transfer (FRET) measurements on freely diffusing single protein molecules. Properties of these subpopulations were directly monitored as a function of guanidinium chloride concentration. It is shown that new information about different aspects of the protein folding reaction can be extracted from such subpopulation properties. Shifts in the mean transfer efficiencies are discussed, FRET efficiency distributions are translated into potentials, and denaturation curves are directly plotted from the areas of the FRET peaks. Changes in stability caused by mutation also are measured by comparing pseudo wild-type CI2 with a destabilized mutant (K17G). Current limitations and future possibilities and prospects for single-pair FRET protein folding investigations are discussed.

Folding and activity of circularly permuted forms of a polytopic membrane protein

The transmembrane subunit of the Glc transporter (IICBGlc), which mediates uptake and concomitant phosphorylation of glucose, spans the membrane eight times. Variants of IICBGlc with the native N and C termini joined and new N and C termini in the periplasmic and cytoplasmic surface loops were expressed in Escherichia coli. In vivo transport/in vitro phosphotransferase activities of the circularly permuted variants with the termini in the periplasmic loops 1 to 4 were 35/58, 32/37, 0/3, and 0/0% of wild type, respectively. The activities of the variants with the termini in the cytoplasmic loops 1 to 3 were 0/25, 0/4 and 24/70, respectively. Fusion of alkaline phosphatase to the periplasmic C termini stabilized membrane integration and increased uptake and/or phosphorylation activities. These results suggest that internal signal anchor and stop transfer sequences can function as N-terminal signal sequences in a circularly permuted α-helical bundle protein and that the orientation of transmembrane segments is determined by the amino acid sequence and not by the sequential appearance during translation. Of the four IICBGlc variants with new termini in periplasmic loops, only the one with the discontinuity in loop 4 is inactive. The sequences of loop 4 and of the adjacent TM7 and TM8 are conserved in all phosphoenolpyruvate-dependent carbohydrate:phosphotransferase system transporters of the glucose family.

Coupling between protein folding and allostery in the GroE chaperonin system

GroEL is an allosteric protein that facilitates protein folding in an ATP-dependent manner. Herein, the relationship between cooperative ATP binding by GroEL and the kinetics of GroE-assisted folding of two substrates with different GroES dependence, mouse dihydrofolate reductase (mDHFR) and mitochondrial malate dehydrogenase, is examined by using cooperativity mutants of GroEL. Strong intra-ring positive cooperativity in ATP binding by GroEL decreases the rate of GroEL-assisted mDHFR folding owing to a slow rate of the ATP-induced transition from the protein-acceptor state to the protein-release state. Inter-ring negative cooperativity in ATP binding by GroEL is found to affect the kinetic partitioning of mDHFR, but not of mitochondrial malate dehydrogenase, between folding in solution and folding in the cavity underneath GroES. Our results show that protein folding by this “two-stroke motor” is coupled to cooperative ATP binding.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Denaturant-induced movement of the transition state of protein folding revealed by high-pressure stopped-flow measurements

The small all-β protein tendamistat folds and unfolds with two-state kinetics. We determined the volume changes associated with the folding process by performing kinetic and equilibrium measurements at variable pressure between 0.1 and 100 MPa (1 to 1,000 bar). GdmCl-induced equilibrium unfolding transitions reveal that the volume of the native state is increased by 41.4 ± 2.0 cm3/mol relative to the unfolded state. This value is virtually independent of denaturant concentration. The use of a high-pressure stopped-flow instrument enabled us to measure the activation volumes for the refolding (ΔVf0‡) and unfolding reaction (ΔVu0‡) over a broad range of GdmCl concentrations. The volume of the transition state is 60% native-like (ΔVf0‡ = 25.0 ± 1.2 cm3/mol) in the absence of denaturant, indicating partial solvent accessibility of the core residues. The volume of the transition state increases linearly with denaturant concentration and exceeds the volume of the native state above 6 M GdmCl. This result argues for a largely desolvated transition state with packing deficiencies at high denaturant concentrations and shows that the structure of the transition state depends strongly on the experimental conditions.

Structure and function in rhodopsin: Further elucidation of the role of the intradiscal cysteines, Cys-110, -185, and -187, in rhodopsin folding and function

The disulfide bond between Cys-110 and Cys-187 in the intradiscal domain is required for correct folding in vivo and function of mammalian rhodopsin. Misfolding in rhodopsin, characterized by the loss of ability to bind 11-cis-retinal, has been shown to be caused by an intradiscal disulfide bond different from the above native disulfide bond. Further, naturally occurring single mutations of the intradiscal cysteines (C110F, C110Y, and C187Y) are associated with retinitis pigmentosa (RP). To elucidate further the role of every one of the three intradiscal cysteines, mutants containing single-cysteine replacements by alanine residues and the above three RP mutants have been studied. We find that C110A, C110F, and C110Y all form a disulfide bond between C185 and C187 and cause loss of retinal binding. C185A allows the formation of a C110–C187 disulfide bond, with wild-type-like rhodopsin phenotype. C187A forms a disulfide bond between C110 and C185 and binds retinal, and the pigment formed has markedly altered bleaching behavior. However, the opsin from the RP mutant C187Y forms no rhodopsin chromophore.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Stretching lattice models of protein folding

A new class of experiments that probe folding of individual protein domains uses mechanical stretching to cause the transition. We show how stretching forces can be incorporated in lattice models of folding. For fast folding proteins, the analysis suggests a complex relation between the force dependence and the reaction coordinate for folding.

Observation of strange kinetics in protein folding

Highly nonexponential folding kinetics in aqueous solution have been observed during temperature jump-induced refolding of two proteins, yeast phosphoglycerate kinase and a ubiquitin mutant. The observations are most easily interpreted in terms of downhill folding, which posits a heterogeneous ensemble of structures en route to the folded state. The data are also reconciled with exponential kinetics measured under different experimental conditions and with titration experiments indicating cooperative folding.

Structural characterization of an engineered tandem repeat contrasts the importance of context and sequence in protein folding

To test a different approach to understanding the relationship between the sequence of part of a protein and its conformation in the overall folded structure, the amino acid sequence corresponding to an α-helix of T4 lysozyme was duplicated in tandem. The presence of such a sequence repeat provides the protein with “choices” during folding. The mutant protein folds with almost wild-type stability, is active, and crystallizes in two different space groups, one isomorphous with wild type and the other with two molecules in the asymmetric unit. The fold of the mutant is essentially the same in all cases, showing that the inserted segment has a well-defined structure. More than half of the inserted residues are themselves helical and extend the helix present in the wild-type protein. Participation of additional duplicated residues in this helix would have required major disruption of the parent structure. The results clearly show that the residues within the duplicated sequence tend to maintain a helical conformation even though the packing interactions with the remainder of the protein are different from those of the original helix. It supports the hypothesis that the structures of individual α-helices are determined predominantly by the nature of the amino acids within the helix, rather than the structural environment provided by the rest of the protein.

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.

Temperature and length scale dependence of hydrophobic effects and their possible implications for protein folding

The Lum–Chandler–Weeks theory of hydrophobicity [Lum, K., Chandler, D. & Weeks, J. D. (1999) J. Phys. Chem. 103, 4570–4577] is applied to treat the temperature dependence of hydrophobic solvation in water. The application illustrates how the temperature dependence for hydrophobic surfaces extending less than 1 nm differs significantly from that for surfaces extending more than 1 nm. The latter is the result of water depletion, a collective effect, that appears at length scales of 1 nm and larger. Because of the contrasting behaviors at small and large length scales, hydrophobicity by itself can explain the variable behavior of entropies of protein folding.

Folding simulations of a three-stranded antiparallel β-sheet peptide

Protein folding is a grand challenge of the postgenomic era. In this paper, 58 folding events sampled during 47 molecular dynamics trajectories for a total simulation time of more than 4 μs provide an atomic detail picture of the folding of a 20-residue synthetic peptide with a stable three-stranded antiparallel β-sheet fold. The simulations successfully reproduce the NMR solution conformation, irrespective of the starting structure. The sampling of the conformational space is sufficient to determine the free energy surface and localize the minima and transition states. The statistically predominant folding pathway involves the formation of contacts between strands 2 and 3, starting with the side chains close to the turn, followed by association of the N-terminal strand onto the preformed 2–3 β-hairpin. The folding mechanism presented here, formation of a β-hairpin followed by consolidation, is in agreement with a computational study of the free energy surface of another synthetic three-stranded antiparallel β-sheet by Bursulaya and Brooks [(1999) J. Am. Chem. Soc. 121, 9947–9951]. Hence, it might hold in general for antiparallel β-sheets with short turns.

Absence of stable intermediates on the folding pathway of barnase

Barnase is one of the few protein models that has been studied extensively for protein folding. Previous studies led to the conclusion that barnase folds through a very stable submillisecond intermediate (≈3 kcal/mol). The structure of this intermediate was characterized intensively by using a protein engineering approach. This intermediate has now been reexamined with three direct and independent methods. (i) Hydrogen exchange experiments show very small protection factors (≈2) for the putative intermediate, indicating a stability of ≈0.0 kcal/mol. (ii) Denaturant-dependent unfolding of the putative intermediate is noncooperative and indicates a stability less than 0.0 kcal/mol. (iii) The logarithm of the unfolding rate constant of native barnase vs. denaturant concentrations is not linear. Together with the measured rate (“I” to N), this nonlinear behavior accounts for almost all of the protein stability, leaving only about 0.3 kcal/mol that could be attributed to the rapidly formed intermediate. Other observations previously interpreted to support the presence of an intermediate are now known to have alternative explanations. These results cast doubts on the previous conclusions on the nature of the early folding state in barnase and therefore should have important implications in understanding the early folding events of barnase and other proteins in general.