Ubiquitin folds via a flip-twist-lock mechanism

To perform specific functional activities, the majority of proteins should fold into their distinct three-dimensional conformations. However, the biologically active conformation of a protein is generally found to be marginally stable than the other conformations that the chain can adopt. How a protein finds its native conformation from its post-synthesis unfolded structure in a complex conformational landscape is the unsolved question that still drives the protein folding community. Here, we report the folding mechanism of a globular protein, ubiquitin, from its chemically denatured state using all-atom molecular dynamics simulations. From the kinetic analysis of the simulated trajectories we show that the folding process can be described by the hydrophobic collapse mechanism, initiated by the “dewetting transition”, and subsequently assisted by the origination of an N-terminal folding nucleus, and finally supported by a native salt-bridge interaction between K11 and E34. We show that ubiquitin folds via an intermediate. Finally, we confirm the presence of “biological water” and explain its role to the folding process.     

Directional structural features of globular proteins

We describe a novel presentation of the conformation of the backbone atoms for proteins of known structure. Given the C^α atom cartesian co-ordinates from X-ray crystallography, a matrix is calculated, where the ijth element of the matrix is the cosine of the angle between the direction of the chain at residue i and the direction of the chain at residue j. These “direction matrices” have distinctive patterns which correspond to α-helix, extended structure, straight or bent segments, “superhelix”, and many other important structural features. We discuss the direction matrices for a number of proteins, and make some generalizations on the basic principles of protein folding.     

Transition state of the rate-limiting step of heat denaturation of Cry3A delta-endotoxin

Heat denaturation of Cry3A delta-endotoxin from Bacillus thuringiensis var. tenebrionis and its 55 kDa fragment was studied by differential scanning microcalorimetry at low pH. Analysis of the calorimetric data has shown that denaturation of Cry3A delta-endotoxin is a nonequilibrium process at heating rates from 0. 125 to 2 K/min. This means that the stability of delta-endotoxin (the apparent temperature of denaturation Tm) under these conditions is under kinetic control rather than under thermodynamic control. It has been shown that heat denaturation of this protein is a one-step kinetic process. The enthalpy of the process and its activation energy were measured as functions of temperature. The data obtained allow confirmation of the fact that the conformation of delta-endotoxin at the transition state only slightly differs from its native conformation with respect to compactness and extent of hydration. The comparison of the activation energy for intact delta-endotoxin and the 55 kDa fragment showed that the transition of the molecule to a transition state does not cause any changes in the conformation of three N-terminal alpha-helices. Complete removal of the N-terminal domain of delta-endotoxin and 40 amino acids from the C-terminus beta-sheet domain III causes an irreversible loss of the tertiary structure. Thus, during protein folding the nucleation core determining protein stability does not involve its three initial alpha-helices but may include the remaining alpha-helices of the N-terminal domain. The functional significance of peculiarities of structure arrangement of the delta-endotoxin molecule is discussed.     

Extended Structures in RNA Folding Intermediates Are Due to Nonnative Interactions Rather than Electrostatic Repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg²⁺ concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.     

Folding, quality control, and secretion of pancreatic ribonuclease in live cells

Although bovine pancreatic RNase is one of the best characterized proteins in respect to structure and in vitro refolding, little is known about its synthesis and maturation in the endoplasmic reticulum (ER) of live cells. We expressed the RNase in live cells and analyzed its folding, quality control, and secretion using pulse-chase analysis and other cell biological techniques. In contrast to the slow in vitro refolding, the protein folded almost instantly after translation and translocation into the ER lumen (t_½ < 3 min). Despite high stability of the native protein, only about half of the RNase reached a secretion competent, monomeric form and was rapidly transported from the rough ER via the Golgi complex (t_½ = 16 min) to the extracellular space (t_½ = 35 min). The rest remained in the ER mainly in the form of dimers and was slowly degraded. The dimers were most likely formed by C-terminal domain swapping since mutation of Asn¹¹³, a residue that stabilizes such dimers, to Ser increased the efficiency of secretion from 59 to 75%. Consistent with stringent ER quality control in vivo, the secreted RNase in the bovine pancreas was mainly monomeric, whereas the enzyme present in the cells also contained 20% dimers. These results suggest that the efficiency of secretion is not only determined by the stability of the native protein but by multiple factors including the stability of secretion-incompetent side products of folding. The presence of N-glycans had little effect on the folding and secretion process.     

Ab initio folding of a trefoil‐fold motif reveals structural similarity with a β‐propeller blade motif

Many protein architectures exhibit evidence of internal rotational symmetry postulated to be the result of gene duplication/fusion events involving a primordial polypeptide motif. A common feature of such structures is a domain‐swapped arrangement at the interface of the N‐ and C‐termini motifs and postulated to provide cooperative interactions that promote folding and stability. De novo designed symmetric protein architectures have demonstrated an ability to accommodate circular permutation of the N‐ and C‐termini in the overall architecture; however, the folding requirement of the primordial motif is poorly understood, and tolerance to circular permutation is essentially unknown. The β‐trefoil protein fold is a threefold‐symmetric architecture where the repeating ~42‐mer “trefoil‐fold” motif assembles via a domain‐swapped arrangement. The trefoil‐fold structure in isolation exposes considerable hydrophobic area that is otherwise buried in the intact β‐trefoil trimeric assembly. The trefoil‐fold sequence is not predicted to adopt the trefoil‐fold architecture in ab initio folding studies; rather, the predicted fold is closely related to a compact “blade” motif from the β‐propeller architecture. Expression of a trefoil‐fold sequence and circular permutants shows that only the wild‐type N‐terminal motif definition yields an intact β‐trefoil trimeric assembly, while permutants yield monomers. The results elucidate the folding requirements of the primordial trefoil‐fold motif, and also suggest that this motif may sample a compact conformation that limits hydrophobic residue exposure, contains key trefoil‐fold structural features, but is more structurally homologous to a β‐propeller blade motif.     

Identification and characterization of structural domains of human ERp57: association with calreticulin requires several domains

The amino acid sequence of ERp57, which functions in the endoplasmic reticulum together with the lectins calreticulin and calnexin to achieve folding of newly synthesized glycoproteins, is highly similar to that of protein disulfide isomerase (PDI), but they have their own distinct roles in protein folding. We have characterized the domain structure of ERp57 by limited proteolysis and N-terminal sequencing and have found it to be similar but not identical to that of PDI. ERp57 had three major protease-sensitive regions, the first of which was located between residues 120 and 150, the second between 201 and 215, and the third between 313 and 341, the data thus being consistent with a four-domain structure abb'a'. Recombinant expression in Escherichia coli was used to verify the domain boundaries. Each single domain and a b'a' double domain could be produced in the form of soluble, folded polypeptides, as verified by circular dichroism spectra and urea gradient gel electrophoresis. When the ability of ERp57 and its a and a' domains to fold denatured RNase A was studied by electrospray mass analyses, ERp57 markedly enhanced the folding rate at early time points, although less effectively than PDI, but was an ineffective catalyst of the overall process. The a and a' domains produced only minor, if any, increases in the folding rate at the early stages and no increase at the late stages. Interaction of the soluble ERp57 domains with the P domain of calreticulin was studied by chemical cross-linking in vitro. None of the single ERp57 domains nor the b'a' double domain could be cross-linked to the P domain, whereas cross-linking was obtained with a hybrid ERpabb'PDIa'c polypeptide but not with ERpabPDIb'a'c, indicating that multiple domains are involved in this protein-protein interaction and that the b' domain of ERp57 cannot be replaced by that of PDI.     

The soluble, periplasmic domain of OmpA folds as an independent unit and displays chaperone activity by reducing the self-association propensity of the unfolded OmpA transmembrane β-barrel

OmpA is one of only a few transmembrane proteins whose folding and stability have been investigated in detail. However, only half of the OmpA mass encodes its transmembrane β-barrel; the remaining sequence is a soluble domain that is localized to the periplasmic side of the outer membrane. To understand how the OmpA periplasmic domain contributes to the stability and folding of the full-length OmpA protein, we cloned, expressed, purified and studied the OmpA periplasmic domain independently of the OmpA transmembrane β-barrel region. Our experiments showed that the OmpA periplasmic domain exists as an independent folding unit with a free energy of folding equal to − 6.2 (± 0.1) kcal mol^-1 at 25 °C. Using circular dichroism, we determined that the OmpA periplasmic domain adopts a mixed alpha/beta secondary structure, a conformation that has previously been used to describe the partially folded non-native state of the full-length OmpA. We further discovered that the OmpA periplasmic domain reduces the self-association propensity of the unfolded barrel domain, but only when covalently attached (in cis). In vitro folding experiments showed that self-association competes with β-barrel folding when allowed to occur before the addition of membranes, and the periplasmic domain enhances the folding efficiency of the full-length protein by reducing its self-association. These results identify a novel chaperone function for the periplasmic domain of OmpA that may be relevant for folding in vivo. We have also extensively investigated the properties of the self-association reaction of unfolded OmpA and found that the transmembrane region must form a critical nucleus comprised of three molecules before undergoing further oligomerization to form large molecular weight species. Finally, we studied the conformation of the unfolded OmpA monomer and found that the folding-competent form of the transmembrane region adopts an expanded conformation, which is in contrast to previous studies that have suggested a collapsed unfolded state.     

A Single Mutation at Residue 25 Populates the Folding Intermediate of E. coli RNase H and Reveals a Highly Dynamic Partially Folded Ensemble

Understanding the nature of partially folded intermediates transiently populated during protein folding is important for understanding both protein folding and misfolding. These ephemeral species, however, often elude direct experimental characterization. The well-characterized protein ribonuclease H (RNase H) from Escherichia coli populates an on-pathway intermediate identified in both bulk studies and single-molecule mechanical unfolding experiments. Here, we set out to trap the transient intermediate of RNase H at equilibrium by selectively destabilizing the region of the protein known to be unfolded in this species. Surprisingly, a single change at Ile25 (I25A) resulted in the equilibrium population of the intermediate under near-native conditions. The intermediate was undetectable in a series of heteronuclear single quantum coherences, revealing the dynamic nature of this partially unfolded form on the timescale of NMR detection. This result is in contrast to studies in which the structures of trapped intermediates are solved by NMR, indicating that they are well packed and native-like. The dynamic nature of the RNase H intermediate may be important for its role as an on-pathway, productive species that promotes efficient folding.     

Role of two proline-containing turns in the folding of porcine ribonuclease

Unfolded ribonuclease (RNase) from porcine pancreas consists of a mixture of fast and slow-refolding species. The equilibrium distribution of these species differs strongly from other homologous RNases, because an additional proline residue is present at position 115 of the porcine protein. The major slow-folding species of porcine RNase contains incorrect proline isomers at Pro93 and at Pro114-Pro115. Both positions are presumably part of beta-turn structures in the native protein, as deduced from the structure of the homologous bovine RNase A. The folding kinetics of these molecules depend strongly on the conditions used. Under unfavorable conditions (near the unfolding transition), refolding is virtually blocked by the presence of the incorrect proline peptide bonds and partially folded intermediates with incorrect isomers could not be detected. As a consequence, folding is very slow under such conditions and the re-isomerization of Pro114-Pro115 is the first and rate-limiting step of folding. Under strongly native conditions (such as in the presence of ammonium sulfate), refolding is much faster. A largely folded intermediate accumulates with the turns around Pro93 and Pro114-Pro115 still in the non-native conformation. These results suggest that incorrect proline isomers strongly influence protein folding and that, under favorable conditions, the polypeptide chain can fold with two beta-turns locked into a non-native conformation. We conclude, therefore, that early formation of correct turn structure is not necessarily required for protein folding. However, the presence of incorrect turns, locked-in by non-native proline isomers, strongly decreases the rate of refolding. Alternative pathways of folding exist. The choice of pathway depends on the number and distribution of incorrect proline isomers and on the folding conditions.     

Solution NMR evidence for a cis Tyr-Ala peptide group in the structure of [Pro93Ala] bovine pancreatic ribonuclease A

Proline peptide group isomerization can result in kinetic barriers in protein folding. In particular, the cis proline peptide conformation at Tyr92-Pro93 of bovine pancreatic ribonuclease A (RNase A) has been proposed to be crucial for chain folding initiation. Mutation of this proline-93 to alanine results in an RNase A molecule, P93A, that exhibits unfolding/refolding kinetics consistent with a cis Tyr92-Ala93 peptide group conformation in the folded structure (Dodge RW, Scheraga HA, 1996, Biochemistry 35:1548-1559). Here, we describe the analysis of backbone proton resonance assignments for P93A together with nuclear Overhauser effect data that provide spectroscopic evidence for a type VI beta-bend conformation with a cis Tyr92-Ala93 peptide group in the folded structure. This is in contrast to the reported X-ray crystal structure of [Pro93Gly]-RNase A (Schultz LW, Hargraves SR, Klink TA, Raines RT, 1998, Protein Sci 7:1620-1625), in which Tyr92-Gly93 forms a type-II beta-bend with a trans peptide group conformation. While a glycine residue at position 93 accommodates a type-II bend (with a positive value of phi93), RNase A molecules with either proline or alanine residues at this position appear to require a cis peptide group with a type-VI beta-bend for proper folding. These results support the view that a cis Pro93 conformation is crucial for proper folding of wild-type RNase A.     

The Folding Trajectory of RNase H Is Dominated by Its Topology and Not Local Stability: A Protein Engineering Study of Variants that Fold via Two-State and Three-State Mechanisms

Proteins can sample a variety of partially folded conformations during the transition between the unfolded and native states. Some proteins never significantly populate these high-energy states and fold by an apparently two-state process. However, many proteins populate detectable, partially folded forms during the folding process. The role of such intermediates is a matter of considerable debate. A single amino acid change can convert Escherichia coli ribonuclease H from a three-state folder that populates a kinetic intermediate to one that folds in an apparent two-state fashion. We have compared the folding trajectories of the three-state RNase H and the two-state RNase H, proteins with the same native-state topology but altered regional stability, using a protein engineering approach. Our data suggest that both versions of RNase H fold through a similar trajectory with similar high-energy conformations. Mutations in the core and the periphery of the protein affect similar aspects of folding for both variants, suggesting a common trajectory with folding of the core region followed by the folding of the periphery. Our results suggest that formation of specific partially folded conformations may be a general feature of protein folding that can promote, rather than hinder, efficient folding.     

Molecular chaperones in cellular protein folding

The discovery of “molecular chaperones” has dramatically changed our concept of cellular protein folding. Rather than folding spontaneously, most newly synthesized polypeptide chains seem to acquire their native conformation in a reaction mediated by these versatile helper proteins. Understanding the structure and function of molecular chaperones is likely to yield useful applications for medicine and biotechnology in the future.     

Crystal Structures of a Group II Chaperonin Reveal the Open and Closed States Associated with the Protein Folding Cycle

Chaperonins are large protein complexes consisting of two stacked multisubunit rings, which open and close in an ATP-dependent manner to create a protected environment for protein folding. Here, we describe the first crystal structure of a group II chaperonin in an open conformation. We have obtained structures of the archaeal chaperonin from Methanococcus maripaludis in both a peptide acceptor (open) state and a protein folding (closed) state. In contrast with group I chaperonins, in which the equatorial domains share a similar conformation between the open and closed states and the largest motions occurs at the intermediate and apical domains, the three domains of the archaeal chaperonin subunit reorient as a single rigid body. The large rotation observed from the open state to the closed state results in a 65% decrease of the folding chamber volume and creates a highly hydrophilic surface inside the cage. These results suggest a completely distinct closing mechanism in the group II chaperonins as compared with the group I chaperonins.     

Folding Simulations of a De Novo Designed Protein with a βαβ Fold

βαβ structural motifs are commonly used building blocks in protein structures containing parallel β-sheets. However, to our knowledge, no stand-alone βαβ structure has been observed in nature to date. Recently, for the first time that we know of, a small protein with an independent βαβ structure (DS119) was successfully designed in our laboratory. To understand the folding mechanism of DS119, in the study described here, we carried out all-atom molecular dynamics and coarse-grained simulations to investigate its folding pathways and energy landscape. From all-atom simulations, we successfully observed the folding event and got a stable folded structure with a minimal root mean-square deviation of 2.6 Å with respect to the NMR structure. The folding process can be described as a fast collapse phase followed by rapid formation of the central helix, and then slow formation of a parallel β-sheet. By using a native-centric Gō-like model, the cooperativity of the system was characterized in terms of the calorimetric criterion, sigmoidal transitions, conformation distribution shifts, and free-energy profiles. DS119 was found to be an incipient downhill folder that folds more cooperatively than a downhill folder, but less cooperatively than a two-state folder. This may reflect the balance between the two structural elements of DS119: the rapidly formed α-helix and the slowly formed parallel β-sheet. Folding times estimated from both the all-atom simulations and the coarse-grained model were at microsecond level, making DS119 another fast folder. Compared to fast folders reported previously, DS119 is, to the best of our knowledge, the first that exhibits a parallel β-sheet.     

The stability of the TIM-barrel domain of a psychrophilic chitinase

Chitinase 60 from the psychrophilic bacterium Moritella marina (MmChi60) is a four-domain protein whose structure revealed flexible hinge regions between the domains, yielding conformations in solution that range from fully extended to compact. The catalytic domain is a shallow-grooved TIM-barrel. Heat-induced denaturation experiments of the wild-type and mutants resulting from the deletions of the two-Ig-like domains and the chitin binding domain reveal calorimetric profiles that are consistent with non-collaborative thermal unfolding of the individual domains, a property that must be associated to the “hinge-regions”. The calorimetric measurements of the (β/α)₈ catalytic domain reveal that the thermal unfolding is a slow-relaxation transition exhibiting a stable, partially structured intermediate state. Circular dichroism provides evidence that the intermediate exhibits features of a molten globule i.e., loss of tertiary structure while maintaining the secondary structural elements of the native. GdnHCl-induced denaturation studies of the TIM-barrel demonstrate an extraordinarily high resistance to the denaturant. Slow-relaxation kinetics characterize the unfolding with equilibration times exceeding six days, a property that is for the first time observed for a psychrophilic TIM barrel. On the other hand, the thermodynamic stability is ΔG=6.75±1.3 kcal/mol, considerably lower than for structural-insertions-containing barrels. The mutant E153Q used for the crystallographic studies of MmChi60 complexes with NAG ligands has a much lower stability than the wild-type.     

Refined models for computer simulation of protein folding. Applications to the study of conserved secondary structure and flexible hinge points during the folding of pancreatic trypsin inhibitor

A new model and parameters are proposed for the computer simulation of protein folding, which satisfy requirements for a fully automatic simulation as discussed in recent critical reviews.The parameters were obtained, refined or checked by empirical observations on proteins of known sequence and conformation, in order to avoid as much as possible theoretical deductions about the nature of the interactions between groups in proteins, which may not be justified by the current status of the art.The major improvement over previous methods is to retain a more realistic and complete representation of the protein backbone, and to alternatively reduce the number of variables by coupling their behaviour. As an example, the method is applied to simulate the folding of pancreatic trypsin inhibitor, and leads to a root-mean-square fit of 6·0 Å with good secondary structure.This also allows a more detailed examination of secondary structure transitions during protein folding than has been possible hitherto. Although, in the simulation discussed most extensively, the advantage of initial statistical predictions is demonstrated, the secondary structure was free to change in the simulation. A simulation from an extended chain is also described, and refinements tested.By observing changes in secondary structure during the simulated folding, it is shown that α-helices and extended chain regions predicted at the outset, or formed early in the simulation, are conserved, and that certain residues are crucial as flexible hinge-points to bring the secondary structure together in order to achieve tertiary packing.In view of recent debate about the importance of glycyl residues as hinge-points, and the danger of imparting glycyl-like backbone behaviour to non-glycyl residues suspected to be hinge-points, it is of considerable interest that the hinge-point residues identified by us are not, in general, glycyl residues. This makes an important distinction between a “reverse turn region”, for which glycine is statistically a strong candidate, and a hinge-point in the protein backbone. It is discussed that reverse turns are locally determined and likely to be fairly stable during the folding process, while hinge-points are determined by tertiary interactions. This distinction, implicit in most papers concerned with statistical methods of secondary structure prediction, has not been made clearly in recent reports of folding simulations.     

Laser Raman scattering of cobramine B, a basic protein from cobra venom

Cobramine B, a small basic protein from cobra venom, is selected as a model for studying the scattering intensity of tyrosyl ring vibrations in the Raman spectra of proteins. All three tyrosines in this protein appear to be “buried” in the interior of the molecule and probably involved in interactions which are similar to those of the three “buried” tyrosines in RNase A when it is dissolved in water. Spectral evidence is presented and discussed. The Raman spectra in the 300–1800 cm^?1 region of cobramine B in the solid and solution are compared quantitatively. Several differences exist between the two spectra and may be interpreted in terms of difference in conformation. In the amide I region, a strong single line was observed at 1672 cm^?1 both in the solid and solution spectra, suggesting that this protein may contain a large fraction of antiparallel-β structure. This is supported by the presence of a line at 1235 cm^?1 in the amide III region, which is also characteristic of β-structure. The resolved peaks at 1254 and 1270 cm^?1 indicate the coexistence of some hydrogen-bonded random-coil and some α-helix with the β-structure.     

Minichaperone (GroEL191-345) mediated folding of MalZ proceeds by binding and release of native and functional intermediates

The isolated apical domain of GroEL consisting of residues 191–345 (known as “minichaperone”) binds and assists the folding of a wide variety of client proteins without GroES and ATP, but the mechanism of its action is still unknown. In order to probe into the matter, we have examined minichaperone-mediated folding of a large aggregation prone protein Maltodextrin-glucosidase (MalZ). The key objective was to identify whether MalZ exists free in solution, or remains bound to, or cycling on and off the minichaperone during the refolding process. When GroES was introduced during refolding process, production of the native MalZ was inhibited. We also observed the same findings with a trap mutant of GroEL, which stably captures a predominantly non-native MalZ released from minichaperone during refolding process, but does not release it. Tryptophan and ANS fluorescence measurements indicated that refolded MalZ has the same structure as the native MalZ, but that its structure when bound to minichaperone is different. Surface plasmon resonance measurements provide an estimate for the equilibrium dissociation constant KD for the MalZ-minichaperone complex of 0.21 ± 0.04 μM, which are significantly higher than for most GroEL clients. This showed that minichaperone interacts loosely with MalZ to allow the protein to change its conformation and fold while bound during the refolding process. These observations suggest that the minichaperone works by carrying out repeated cycles of binding aggregation-prone protein MalZ in a relatively compact conformation and in a partially folded but active state, and releasing them to attempt to fold in solution.