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.

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.

The quorum-sensing transcriptional regulator TraR requires its cognate signaling ligand for protein folding, protease resistance, and dimerization

Complexes between the quorum-sensing regulator TraR and its inducing ligand autoinducer (AAI) are soluble in Escherichia coli, whereas apo-TraR is almost completely insoluble. Here we show that the lack of soluble TraR is due in large part to rapid proteolysis, inasmuch as apo-TraR accumulated to high levels in an E. coli strain deficient in Clp and Lon proteases. In pulse labeling experiments, AAI protected TraR against proteolysis only when it was added before the radiolabel. This observation indicates that TraR proteins can productively bind AAI only during their own synthesis on polysomes, whereas fully synthesized apo-TraR proteins are not functional AAI receptors. Purified apo-TraR was rapidly degraded by trypsin to oligopeptides, whereas TraR–AAI complexes were more resistant to trypsin and were cleaved at discrete interdomain linkers, indicating that TraR requires AAI to attain its mature tertiary structure. TraR–AAI complexes eluted from a gel filtration column as dimers and bound DNA as dimers. In contrast, apo-TraR was monomeric, and incubation with AAI under a variety of conditions did not cause dimerization. We conclude that AAI is critical for the folding of nascent TraR protein into its mature tertiary structure and that full-length apo-TraR cannot productively bind AAI and is consequently targeted for rapid proteolysis.

Dissection of the ion-induced folding of the hammerhead ribozyme using 19F NMR

We have used 19F NMR to analyze the metal ion-induced folding of the hammerhead ribozyme by selective incorporation of 5fluorouridine. We have studied the chemical shift and linewidths of 19F resonances of 5-fluorouridine at the 4 and 7 positions in the ribozyme core as a function of added Mg2+. The data fit well to a simple two-state model whereby the formation of domain 1 is induced by the noncooperative binding of Mg2+ with an association constant in the range of 100 to 500 M−1, depending on the concentration of monovalent ions present. The results are in excellent agreement with data reporting on changes in the global shape of the ribozyme. However, the NMR experiments exploit reporters located in the center of the RNA sections undergoing the folding transitions, thereby allowing the assignment of specific nucleotides to the separate stages. The results define the folding pathway at high resolution and provide a time scale for the first transition in the millisecond range.

Effects of inefficient cleavage of the signal sequence of HIV-1 gp 120 on its association with calnexin, folding, and intracellular transport.

The HIV-1 envelope glycoprotein gp120 displays inefficient intracellular transport, which is caused by its retention in the endoplasmic reticulum. Coexpression in insect cells (Sf9) of HIV-1 gp120 with calnexin has shown that their interaction was modulated by the signal sequence of HIV-1 gp120. gp120, with its natural signal sequence, showed a prolonged association with calnexin with a t1/2 of greater than 20 min. Replacement of the natural signal sequence with the signal sequence from mellitin led to a decreased time of association of gp120 with calnexin (t1/2 < 10 min). These different times of calnexin association coincided both with the folding of gp120 as measured by the ability of bind CD4 and with endoplasmic reticulum to Golgi transport as analyzed by the acquisition of partial endoglycosidase H resistance. Using a monospecific antibody to the HIV-1 gp120 natural signal peptide, we showed that calnexin associated with N-glycosylated but uncleaved gp120. Only after dissociation from calnexin was gp120 cleaved, but very inefficiently. Only the small proportion of signal-cleaved gp120 molecules acquired transport competence and were secreted. This is the first report demonstrating the effect of the signal sequence on calnexin association.

Structure and function in rhodopsin: correct folding and misfolding in point mutants at and in proximity to the site of the retinitis pigmentosa mutation Leu-125-->Arg in the transmembrane helix C.

L125R is a mutation in the transmembrane helix C of rhodopsin that is associated with autosomal dominant retinitis pigmentosa. To probe the orientation of the helix and its packing in the transmembrane domain, we have prepared and studied the mutations E122R, I123R, A124R, S127R, L125F, and L125A at, and in proximity to, the above mutation site. Like L125R, the opsin expressed in COS-1 cells from E122R did not bind 11-cis-retinal, whereas those from I123R and S127R formed the rhodopsin chromophore partially. A124R opsin formed the rhodopsin chromophore (lambda max 495 nm) in the dark, but the metarhodopsin II formed on illumination decayed about 6.5 times faster than that of the wild type and was defective in transducin activation. The mutant opsins from L125F and L125A bound 11-cis-retinal only partially, and in both cases, the mixtures of the proteins produced were separated into retinal-binding and non-retinal-binding (misfolded) fractions. The purified mutant rhodopsin from L125F showed lambda max at 500 nm, whereas that from L125A showed lambda max at 503 nm. The mutant rhodopsin L125F showed abnormal bleaching behavior and both mutants on illumination showed destabilized metarhodopsin II species and reduced transducin activation. Because previous results have indicated that misfolding in rhodopsin is due to the formation of a disulfide bond other than the normal disulfide bond between Cys-110 and Cys-187 in the intradiscal domain, we conclude from the misfolding in mutants L125F and L125A that the folding in vivo in the transmembrane domain is coupled to that in the intradiscal domain.

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.

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.

Brownian dynamics simulations of protein folding: Access to milliseconds time scale and beyond

Protein folding occurs on a time scale ranging from milliseconds to minutes for a majority of proteins. Computer simulation of protein folding, from a random configuration to the native structure, is nontrivial owing to the large disparity between the simulation and folding time scales. As an effort to overcome this limitation, simple models with idealized protein subdomains, e.g., the diffusion–collision model of Karplus and Weaver, have gained some popularity. We present here new results for the folding of a four-helix bundle within the framework of the diffusion–collision model. Even with such simplifying assumptions, a direct application of standard Brownian dynamics methods would consume 10,000 processor-years on current supercomputers. We circumvent this difficulty by invoking a special Brownian dynamics simulation. The method features the calculation of the mean passage time of an event from the flux overpopulation method and the sampling of events that lead to productive collisions even if their probability is extremely small (because of large free-energy barriers that separate them from the higher probability events). Using these developments, we demonstrate that a coarse-grained model of the four-helix bundle can be simulated in several days on current supercomputers. Furthermore, such simulations yield folding times that are in the range of time scales observed in experiments.

Definition of amide protection factors for early kinetic intermediates in protein folding

Hydrogen–deuterium exchange experiments have been used previously to investigate the structures of well defined states of a given protein. These include the native state, the unfolded state, and any intermediates that can be stably populated at equilibrium. More recently, the hydrogen–deuterium exchange technique has been applied in kinetic labeling experiments to probe the structures of transiently formed intermediates on the kinetic folding pathway of a given protein. From these equilibrium and nonequilibrium studies, protection factors are usually obtained. These protection factors are defined as the ratio of the rate of exchange of a given backbone amide when it is in a fully solvent-exposed state (usually obtained from model peptides) to the rate of exchange of that amide in some state of the protein or in some intermediate on the folding pathway of the protein. This definition is straightforward for the case of equilibrium studies; however, it is less clear-cut for the case of transient kinetic intermediates. To clarify the concept for the case of burst-phase intermediates, we have introduced and mathematically defined two different types of protection factors: one is Pstruc, which is more related to the structure of the intermediate, and the other is Papp, which is more related to the stability of the intermediate. Kinetic hydrogen–deuterium exchange data from disulfide-intact ribonuclease A and from cytochrome c are discussed to explain the use and implications of these two definitions.

Structure and function of the Bacillus hybrid enzyme GluXyn-1: Native-like jellyroll fold preserved after insertion of autonomous globular domain

The 1,3–1,4-β-glucanase from Bacillus macerans (wtGLU) and the 1,4-β-xylanase from Bacillus subtilis (wtXYN) are both single-domain jellyroll proteins catalyzing similar enzymatic reactions. In the fusion protein GluXyn-1, the two proteins are joined by insertion of the entire XYN domain into a surface loop of cpMAC-57, a circularly permuted variant of wtGLU. GluXyn-1 was generated by protein engineering methods, produced in Escherichia coli and shown to fold spontaneously and have both enzymatic activities at wild-type level. The crystal structure of GluXyn-1 was determined at 2.1 Å resolution and refined to R = 17.7% and R(free) = 22.4%. It shows nearly ideal, native-like folding of both protein domains and a small, but significant hinge bending between the domains. The active sites are independent and accessible explaining the observed enzymatic activity. Because in GluXyn-1 the complete XYN domain is inserted into the compact folding unit of GLU, the wild-type-like activity and tertiary structure of the latter proves that the folding process of GLU does not depend on intramolecular interactions that are short-ranged in the sequence. Insertion fusions of the GluXyn-1 type may prove to be an easy route toward more stable bifunctional proteins in which the two parts are more closely associated than in linear end-to-end protein fusions.

GroEL-GroES-mediated protein folding requires an intact central cavity

The chaperonin GroEL is an oligomeric double ring structure that, together with the cochaperonin GroES, assists protein folding. Biochemical analyses indicate that folding occurs in a cis ternary complex in which substrate is sequestered within the GroEL central cavity underneath GroES. Recently, however, studies of GroEL “minichaperones” containing only the apical substrate binding subdomain have questioned the functional importance of substrate encapsulation within GroEL-GroES complexes. Minichaperones were reported to assist folding despite the fact that they are monomeric and therefore cannot form a central cavity. Here we compare directly the folding activity of minichaperones with that of the full GroEL-GroES system. In agreement with earlier studies, minichaperones assist folding of some proteins. However, this effect is observed only under conditions where substantial spontaneous folding is also observed and is indistinguishable from that resulting from addition of the nonchaperone protein α-casein. By contrast, the full GroE system efficiently promotes folding of several substrates under conditions where essentially no spontaneous folding is observed. These data argue that the full GroEL folding activity requires the intact GroEL-GroES complex, and in light of previous studies, underscore the importance of substrate encapsulation for providing a folding environment distinct from the bulk solution.

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.

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.

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.