The rate-limiting step in the folding of the cis-Pro167Thr mutant of TEM-1 β-lactamase is the trans to cis isomerization of a non-proline peptide bond

The stability and kinetics of unfolding and refolding of the P167T mutant of the TEM-1 β-lactamase have been investigated as a function of guanidine hydrochloride concentration. The activity of the mutant enzyme was not significantly modified, which strongly suggests that the Glu166–Thr167 peptide bond, like the Glu166–Pro167, is cis. The mutation, however, led to a significant decrease in the stability of the native state relative to both the thermodynamically stable intermediate and the fully unfolded state of the protein. In contrast to the two slower phases seen in the refolding of the wild-type enzyme, only one phase was detected in the refolding of the mutant, indicating a determining role of proline 167 in the kinetics of folding of the wild-type enzyme. The former phases are replaced by rapid refolding when the enzyme is unfolded for short periods of time, but the latter is independent of the time of unfolding. The monophasic refolding reaction of the mutant is proposed to reflect mainly the trans→cis isomerization of the Glu166–Thr167 peptide bond. © 1996 John Wiley & Sons, Inc.     

Parallel channels and rate-limiting steps in complex protein folding reactions: prolyl isomerization and the alpha subunit of Trp synthase,a TIM barrel protein

A kinetic folding mechanism for the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli, involving four parallel channels with multiple native, intermediate and unfolded forms, has recently been proposed. The hypothesis that cis/trans isomerization of several Xaa-Pro peptide bonds is the source of the multiple folding channels was tested by measuring the sensitivity of the three rate-limiting phases (tau(1), tau(2), tau(3)) to catalysis by cyclophilin, a peptidyl-prolyl isomerase. Although the absence of catalysis for the tau(1) (fast) phase leaves its assignment ambiguous, our previous mutational analysis demonstrated its connection to the unique cis peptide bond preceding proline 28. The acceleration of the tau(2) (medium) and tau(3) (slow) refolding phases by cyclophilin demonstrated that cis/trans prolyl isomerization is also the source of these phases. A collection of proline mutants, which covered all of the remaining 18 trans proline residues of alphaTS, was constructed to obtain specific assignments for these phases. Almost all of the mutant proteins retained the complex equilibrium and kinetic folding properties of wild-type alphaTS; only the P217A, P217G and P261A mutations caused significant changes in the equilibrium free energy surface. Both the P78A and P96A mutations selectively eliminated the tau(1) folding phase, while the P217M and P261A mutations eliminated the tau(2) and tau(3) folding phases, respectively. The redundant assignment of the tau(1) phase to Pro28, Pro78 and Pro96 may reflect their mutual interactions in non-random structure in the unfolded state. The non-native cis isomers for Pro217 and Pro261 may destabilize an autonomous C-terminal folding unit, thereby giving rise to kinetically distinct unfolded forms. The nature of the preceding amino acid, the solvent exposure, or the participation in specific elements of secondary structure in the native state, in general, are not determinative of the proline residues whose isomerization reactions can limit folding.     

Proline peptide isomerization and protein folding

The unfolding-refolding of proteins is a cooperative process and, as judged by equilibrium properties, occurs in one step involving the native,N, and the unfoldedU, conformational states. Kinetic studies have shown that the denatured protein exists as a mixture of slow-(U)_Sand fast-(U)_Frefolding forms produced by proline peptidecis-trans isomerization. Proline residues inU _Fare in the same configuration as in the native protein while they are in non-native configuration inU _S. For protein folding to occur quicklyU _Smust be converted intoU _F. The fact that the equilibrium and kinetic properties of are the same as those found for prolinecis-trans isomerization taken together with the absence of slow phase in the kinetics of refolding of a protein devoid of proline, support this view. However, the absence of a linear correlation between half-time of reactivation of denatured enzymes and their proline-contents, as well as the dissimilarities in the kinetic properties of in unfolding and refolding experiments are not consistent with the model. Conformational energy calculation and experimental results on refolding of proteins suggest that some proline residues are non-essential. They will not block protein folding even in wrong isomeric form. The native-like folded structure with incorrect proline isomers will serve as intermediate state(s) in which these prolines will more readily isomerize to the correct isomeric form. The picture becomes more complex when one considers the consequence ofcis-trans isomerism of non-proline residues on protein folding.     

Structure of a partially unfolded form of Escherichia coli dihydrofolate reductase provides insight into its folding pathway

Proteins frequently fold via folding intermediates that correspond to local minima on the conformational energy landscape. Probing the structure of the partially unfolded forms in equilibrium under native conditions can provide insight into the properties of folding intermediates. To elucidate the structures of folding intermediates of Escherichia coli dihydrofolate reductase (DHFR), we investigated transient partial unfolding of DHFR under native conditions. We probed the structure of a high‐energy conformation susceptible to proteolysis (cleavable form) using native‐state proteolysis. The free energy for unfolding to the cleavable form is clearly less than that for global unfolding. The dependence of the free energy on urea concentration (m‐value) also confirmed that the cleavable form is a partially unfolded form. By assessing the effect of mutations on the stability of the partially unfolded form, we found that native contacts in a hydrophobic cluster formed by the F‐G and Met‐20 loops on one face of the central β‐sheet are mostly lost in the partially unfolded form. Also, the folded region of the partially unfolded form is likely to have some degree of structural heterogeneity. The structure of the partially unfolded form is fully consistent with spectroscopic properties of the near‐native kinetic intermediate observed in previous folding studies of DHFR. The findings suggest that the last step of the folding of DHFR involves organization in the structure of two large loops, the F‐G and Met‐20 loops, which is coupled with compaction of the rest of the protein.     

Characterization of secondary amide peptide bond isomerization: thermodynamics and kinetics from 2D NMR spectroscopy

Secondary amide cis peptide bonds are of even lower abundance than the cis tertiary amide bonds of prolines, yet they are of biochemical importance. Using 2D NMR exchange spectroscopy (EXSY) we investigated the formation of cis peptide bonds in several oligopeptides: Ac-G-G-G-NH(2) , Ac-I-G-G-NH(2) , Ac-I-G-G-N-NH(2) and its cyclic form: I-G-G-N in dimethylsulfoxide (DMSO). From the NMR studies, using the amide protons as monitors, an occurrence of 0.13-0.23% of cis bonds was obtained at 296 K. The rate constants for the trans to cis conversion determined from 2D EXSY spectroscopy were 4-9 × 10(-3) s(-1) . Multiple minor conformations were detected for most peptide bonds. From their thermodynamic and kinetic properties the cis isomers are distinguished from minor trans isomers that appear because of an adjacent cis peptide bond. Solvent and sequence effects were investigated utilizing N-methylacetamide (NMA) and various peptides, which revealed a unique enthalpy profile in DMSO. The cyclization of a tetrapeptide resulted in greatly lowered cis populations and slower isomerization rates compared to its linear counterpart, further highlighting the impact of structural constraints.     

Structural mechanism governing cis and trans isomeric states and an intramolecular switch for cis/trans isomerization of a non-proline peptide bond observed in crystal structures of scorpion toxins

Non-proline cis peptide bonds have been observed in numerous protein crystal structures even though the energetic barrier to this conformation is significant and no non-prolyl-cis/trans-isomerase has been identified to date. While some external factors, such as metal binding or co-factor interaction, have been identified that appear to induce cis/trans isomerization of non-proline peptide bonds, the intrinsic structural basis for their existence and the mechanism governing cis/trans isomerization in proteins remains poorly understood. Here, we report the crystal structure of a newly isolated neurotoxin, the scorpion alpha-like toxin Buthus martensii Karsch (BmK) M7, at 1.4A resolution. BmK M7 crystallizes as a dimer in which the identical non-proline peptide bond between residues 9 and 10 exists either in the cis conformation or as a mixture of cis and trans conformations in either monomer. We also determined the crystal structures of several mutants of BmK M1, a representative scorpion alpha-like toxin that contains an identical non-proline cis peptide bond as that observed in BmK M7, in which residues within or neighboring the cis peptide bond were altered. Substitution of an aspartic acid residue for lysine at residue 8 in the BmK M1 (K8D) mutant converted the cis form of the non-proline peptide bond 9-10 into the trans form, revealing an intramolecular switch for cis-to-trans isomerization. Cis/trans interconversion of the switch residue at position 8 appears to be sequence-dependent as the peptide bond between residues 9 and 10 retains its wild-type cis conformation in the BmK M1 (K8Q) mutant structure. The structural interconversion of the isomeric states of the BmK M1 non-proline cis peptide bond may relate to the conversion of the scorpion alpha-toxins subgroups.     

A cis-prolyl peptide bond isomerization dominates the folding of the alpha subunit of Trp synthase,a TIM barrel protein

The cis/trans isomerization of prolyl peptide bonds has been suggested to dominate the folding of the alpha subunit of tryptophan synthase from Escherichia coli (alphaTS). To test the role of the unique cis isomer between Asp27 and Pro28, the folding properties of P28A, P28G and G(3)P28G, a three-glycine insertion mutant between Asp27 and Gly28, were investigated using urea as a denaturant. Circular dichroism analysis demonstrated that none of the mutations perturb the secondary structure significantly, although the aromatic side-chain packing is altered for P28A and P28G. All three mutant proteins inherited the three-state thermodynamic behavior observed in wild-type alphaTS, ensuring that the fundamental features of the energy surface are intact. Kinetic studies showed that neither alanine nor glycine substitutions at Pro28 results in the elimination of any slow-refolding phases. By contrast, the G(3)P28G mutant eliminates the fastest of the slow-refolding phases and one of the two unfolding phases. Double-jump experiments on G(3)P28G confirm the assignment of the missing refolding phase to the isomerization of the Asp27-Pro28 peptide bond. These results imply that the local stability conveyed by the tight, overlapping turns containing the cis peptide bond is sufficient to favor the cis isomer for several non-prolyl residues. The free energy required to drive the isomerization reaction is provided by the formation of the stable intermediate, demonstrating that the acquisition of structure and stability is required to induce subsequent rate-limiting steps in the folding of alphaTS.     

A hydrophobic cluster forms early in the folding of dihydrofolate reductase

The rapid kinetic phase that leads from unfolded species to transient folding intermediates in dihydrofolate reductase from Escherichia coli was examined by site-directed mutagenesis and by physicochemical means. The absence of this fluorescence-detected phase in the refolding of the Trp-74Phe mutant protein strongly implies that this early phase in refolding can be assigned to just one of the five Trp residues in the protein, Trp-74. In addition, water-soluble fluorescence quenching agents, iodide and cesium, have a much less significant effect on this early step in refolding than on the slower phases that lead to native and native-like conformers. These and other data imply that an important early event in the folding of dihydrofolate reductase is the formation of a hydrophobic cluster which protects Trp-74 from solvent.     

Breaking non-native hydrophobic clusters is the rate-limiting step in the folding of an alanine-based peptide

The formation mechanism of an alanine-based peptide has been studied by all-atom molecular dynamics simulations with a recently developed all-atom point-charge force field and the Generalize Born continuum solvent model at an effective salt concentration of 0.2M. Thirty-two simulations were conducted. Each simulation was performed for 100 ns. A surprisingly complex folding process was observed. The development of the helical content can be divided into three phases with time constants of 0.06-0.08, 1.4-2.3, and 12-13 ns, respectively. Helices initiate extreme rapidly in the first phase similar to that estimated from explicit solvent simulations. Hydrophobic collapse also takes place in this phase. A folding intermediate state develops in the second phase and is unfolded to allow the peptide to reach the transition state in the third phase. The folding intermediate states are characterized by the two-turn short helices and the transition states are helix-turn-helix motifs-both of which are stabilized by hydrophobic clusters. The equilibrium helical content, calculated by both the main-chain Phi-Psi torsion angles and the main-chain hydrogen bonds, is 64-66%, which is in remarkable agreement with experiments. After corrected for the solvent viscosity effect, an extrapolated folding time of 16-20 ns is obtained that is in qualitative agreement with experiments. Contrary to the prevailing opinion, neither initiation nor growth of the helix is the rate-limiting step. Instead, the rate-limiting step for this peptide is breaking the non-native hydrophobic clusters in order to reach the transition state. The implication to the folding mechanisms of proteins is also discussed.     

Intramolecular catalysis of a proline isomerization reaction in the folding of dihydrofolate reductase.

The cis/trans isomerization of the peptide bond preceding proline residues in proteins can limit the rate at which a protein folds to its native conformation. Mutagenic analyses of dihydrofolate reductase (DHFR) from Escherichia coli show that this isomerization reaction can be intramolecularly catalyzed by a side chain from an amino acid which is distant in sequence but adjacent in the native conformation. The guanidinium NH2 nitrogen of Arg 44 forms one hydrogen bond to the imide nitrogen and a second to the carbonyl oxygen of Pro 66 in wild-type DHFR. Replacement of Arg 44 with Leu results in a change of the nature of the two slow steps in refolding from being limited by the acquisition of secondary and/or tertiary structure to being limited by isomerization. The simultaneous replacement of Pro 66 with Ala (i.e., the Leu 44/Ala 66 double mutant) eliminates this isomerization reaction and once again makes protein folding the limiting process. Apparently, one or both of the hydrogen bonds between Arg 44 and Pro 66 accelerate the isomerization of the Gln 65-Pro 66 peptide bond. The replacement of Arg 44 with Leu affects the kinetics of the slow folding reactions in a fashion which indicates that the crucial hydrogen bonds form in the transition states for the rate-limiting steps in folding.     

Time-resolved FTIR difference spectroscopy as tool for investigating refolding reactions of ribonuclease T1 synchronized with trans --> cis prolyl isomerization

The structurally well-characterized enzyme ribonuclease T1 was used as a model protein to further evaluate time-resolved Fourier transform IR difference spectroscopy in conjunction with temperature-jump techniques as a useful detection technique for protein folding studies. Compared to the wild-type protein, it was confirmed that the lack of one cis-proline bond at position 55 of the S54G/P55N variant is sufficient to significantly simplify and accelerate the refolding process. This result was sustained by the characterization of the early refolding events that occurred within the experimental dead time.     

Dynamics of hydrogen bond desolvation in protein folding

As proteins fold, a progressive structuring, immobilization and eventual exclusion of water surrounding backbone hydrogen bonds takes place. This process turns hydrogen bonds into major determinants of the folding pathway and compensates for the penalty of desolvation of the backbone polar groups. Taken as an average over all hydrogen bonds in a native fold, this extent of protection is found to be nearly ubiquitous. It is dynamically crucial, determining a constraint in the long-time limit behavior of coarse-grained ab initio simulations. Furthermore, an examination of one of the longest available (1micros) all-atom simulations with explicit solvent reveals that this average extent of protection is a constant of motion for the folding trajectory. We propose how such a stabilization is best achieved by clustering five hydrophobes around the backbone hydrogen bonds, an arrangement that yields the optimal stabilization. Our results support and clarify the view that hydrophobic surface burial should be commensurate with hydrogen-bond formation and enable us to define a basic desolvation motif inherent to structure and folding dynamics.     

Probing minimal independent folding units in dihydrofolate reductase by molecular dissection.

Molecular dissection was employed to identify minimal independent folding units in dihydrofolate reductase (DHFR) from Escherichia coli. Eight overlapping fragments of DHFR, spanning the entire sequence and ranging in size from 36 to 123 amino acids, were constructed by chemical cleavage. These fragments were designed to examine the effect of tethering multiple elements of secondary structure on folding and to test if the secondary structural domains represent autonomous folding units. CD and fluorescence spectroscopy demonstrated that six fragments containing up to a total of seven alpha-helices or beta-strands and, in three cases, the adenine binding domain (residues 37-86), are largely disordered. A stoichiometric mixture of the two fragments comprising the large discontinuous domain, 1-36 and 87-159, also showed no evidence for folding beyond that observed for the isolated fragments. A fragment containing residues 1-107 appears to have secondary and tertiary structure; however, spontaneous self-association made it impossible to determine if this structure solely reflects the behavior of the monomeric form. In contrast, a monomeric fragment spanning residues 37-159 possesses significant secondary and tertiary structure. The urea-induced unfolding of fragment 37-159 in the presence of 0.5 M ammonium sulfate was found to be a well-defined, two-state process. The observation that fragment 37-159 can adopt a stable native fold with unique, aromatic side-chain packing is quite striking because residues 1-36 form an integral part of the structural core of the full-length protein.     

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.     

Green‐lighting green fluorescent protein: Faster and more efficient folding by eliminating a cis–trans peptide isomerization event

Wild‐type green fluorescent protein (GFP) folds on a time scale of minutes. The slow step in folding is a cis–trans peptide bond isomerization. The only conserved cis‐peptide bond in the native GFP structure, at P89, was remodeled by the insertion of two residues, followed by iterative energy minimization and side chain design. The engineered GFP was synthesized and found to fold faster and more efficiently than its template protein, recovering 50% more of its fluorescence upon refolding. The slow phase of folding is faster and smaller in amplitude, and hysteresis in refolding has been eliminated. The elimination of a previously reported kinetically trapped state in refolding suggests that X‐P89 is trans in the trapped state. A 2.55 Å resolution crystal structure revealed that the new variant contains only trans‐peptide bonds, as designed. This is the first instance of a computationally remodeled fluorescent protein that folds faster and more efficiently than wild type.     

Microsecond subdomain folding in dihydrofolate reductase

The characterization of microsecond dynamics in the folding of multisubdomain proteins has been a major challenge in understanding their often complex folding mechanisms. Using a continuous-flow mixing device coupled with fluorescence lifetime detection, we report the microsecond folding dynamics of dihydrofolate reductase (DHFR), a two-subdomain α/β/α sandwich protein known to begin folding in this time range. The global dimensions of early intermediates were monitored by Förster resonance energy transfer, and the dynamic properties of the local Trp environments were monitored by fluorescence lifetime detection. We found that substantial collapse occurs in both the locally connected adenosine binding subdomain and the discontinuous loop subdomain within 35 μs of initiation of folding from the urea unfolded state. During the fastest observable ∼ 550 μs phase, the discontinuous loop subdomain further contracts, concomitant with the burial of Trp residue(s), as both subdomains achieve a similar degree of compactness. Taken together with previous studies in the millisecond time range, a hierarchical assembly of DHFR—in which each subdomain independently folds, subsequently docks, and then anneals into the native conformation after an initial heterogeneous global collapse—emerges. The progressive acquisition of structure, beginning with a continuously connected subdomain and spreading to distal regions, shows that chain entropy is a significant organizing principle in the folding of multisubdomain proteins and single-domain proteins. Subdomain folding also provides a rationale for the complex kinetics often observed.     

Peptide fragment studies on the folding elements of dihydrofolate reductase from Escherichia coli

One of the necessary conditions for a protein to be foldable is the presence of a complete set of “folding elements” (FEs) that are short, contiguous peptide segments distributed over an amino acid sequence. The FE‐assembly model of protein folding has been proposed, in which the FEs play a role in guiding structure formation through FE–FE interactions early in folding. However, two major issues remain to be clarified regarding the roles of the FEs in determining protein foldability. Are the FEs AFUs that can form nativelike structures in isolation? Is the presence of only the FEs without mutual connections a sufficient condition for a protein to be foldable? Here, we address these questions using peptide fragments corresponding to the FEs of DHFR from Escherichia coli. We show by CD measurement that the FE peptides are unfolded under the native conditions, and some of them have the propensities toward non‐native helices. MD simulations also show the non‐native helical propensities of the peptides, and the helix contents estimated from the simulations are well correlated with those estimated from the CD in TFE. Thus, the FEs of DHFR are not AFUs, suggesting the importance of the FEs in nonlocal interactions. We also show that equimolar mixtures of the FE peptides do not induce any structural formation. Therefore, mutual connections between the FEs, which should strengthen the nonlocal FE–FE interactions, are also one of the necessary conditions for a protein to be foldable. Proteins 2006. © 2005 Wiley‐Liss, Inc.     

Osmolyte effects on kinetics of FKBP12 C22A folding coupled with prolyl isomerization

Unfolding and refolding kinetics of human FKBP12 C22A were monitored by fluorescence emission over a wide range of urea concentration in the presence and absence of protecting osmolytes glycerol, proline, sarcosine and trimethylamine-N-oxide (TMAO). Unfolding is well described by a mono-exponential process, while refolding required a minimum of two exponentials for an adequate fit throughout the urea concentration range considered. The bi-exponential behavior resulted from complex coupling between protein folding, and prolyl isomerization in the denatured state in which the urea-dependent rate constant for folding was greater than, equal to, and less than the rate constants for prolyl isomerization within the urea concentration range of zero to five molar. Amplitudes and the observed folding and unfolding rate constants were fitted to a reversible three-state model composed of two sequential steps involving the native state and a folding-competent denatured species thermodynamically linked to a folding-incompetent denatured species. Excellent agreement between thermodynamic parameters for FKBP12 C22A folding calculated from the kinetic parameters and those obtained directly from equilibrium denaturation assays provides strong support for the applicability of the mechanism, and provides evidence that FKBP12 C22A folding/unfolding is two-state, with prolyl isomer heterogeneity in the denatured ensemble. Despite the chemical diversity of the protecting osmolytes, they all exhibit the same kinetic behavior of increasing the rate constant of folding and decreasing the rate constant for unfolding. Osmolyte effects on folding/unfolding kinetics are readily explained in terms of principles established in understanding osmolyte effects on protein stability. These principles involve the osmophobic effect, which raises the Gibbs energy of the denatured state due to exposure of peptide backbone, thereby increasing the folding rate. This effect also plays a key role in decreasing the unfolding rate when, as is often the case, the activated complex exposes more backbone than is exposed in the native state.     

Probing the interactions between the folding elements early in the folding of Escherichia coli dihydrofolate reductase by systematic sequence perturbation analysis

One of the necessary conditions for a protein to be foldable is the presence of a complete set of folding elements (FEs) that are short contiguous peptide segments distributed over an amino acid sequence. Previous studies indicated the FE assembly model of protein folding, in which the FEs interact with each other and coalesce to form an intermediate(s) early in the folding reaction. This suggests that a clue to the understanding of the determinants of protein foldability can be found by investigating how the FEs interact with each other early in the folding and thereby elucidating roles of the FEs in protein folding. To reveal the formation process of FE-FE interactions, we studied the early folding events of Escherichia coli dihydrofolate reductase (DHFR) utilizing systematic sequence perturbation analysis. Here, systematic single amino acid substitutions were introduced inside of the FEs (W30X in FE2, V40X in FE3, N59X in FE4, and I155X in FE10; X refers to various amino acid residues), and their kinetic refolding reactions were measured by stopped-flow circular dichroism and fluorescence. We show that the interactions around Trp30 and Ile155 are formed in the burst phase intermediate, while those around Val40 and Asn59 are formed in the transition state of the subsequent folding phase (tau5-phase) and in much later processes, respectively. These and previous results suggest that FE2 and FE10, and also FE1 and FE7, involved in the loop subdomain of DHFR, interact with each other within a millisecond time range, while the stable FE3-FE4 interactions are formed in the later processes. This may highlight the important roles of the FEs mainly inside of the loop subdomain in formation of the burst phase intermediate having a hydrophobic cluster and native-like overall topology and in acquisition of the foldability of DHFR.     

Correlation of folding kinetics with the number and isomerization states of prolines in three homologous proteins of the RNase family

Several studies attribute the slower phases in protein folding to prolyl isomerizations, and several others do not. A correlation exists between the number of prolines in a protein and the complexity of the mechanism with which it folds. In this study, we have demonstrated a direct correlation between the number of cis-prolyl bonds in a native protein and the complexity with which it folds via slower phases by studying the folding of three structurally homologous proteins of the ribonuclease family, namely RNase A, onconase and angiogenin, which differ in the number and isomerization states of their proline residues.