An examination of the involvement of proline peptide isomerization in protein folding.

The proline peptide isomerization model of protein folding predicts that the fraction of denatured polypeptide chains which rapidly fold should be quantitatively related to the numbers of cis and trans proline residues in the folded polypeptide conformation. However, we find that neither the comparative nor the absolute kinetic pattern for folding of the homologous proteins, tuna heart and horse heart ferricytochrome c which differ by one proline residue, is compatible with the quantitative predictions of the proline peptide isomerization model.     

Bioactivity of folding intermediates studied by the recovery of enzymatic activity during refolding

The peptide bond preceding proline residues realizes a cis/trans conformational switch with high switching resistance in native proteins and folding intermediates. Therefore, individual isomers have the potential to differ in bioactivity. However, information about isomer-specific bioactivities is difficult to obtain because of the risk of affecting isomeric distribution by bioactivity assay components.Here we present an approach that allows for the measurement of the recovery of enzymatic activities of wild-type RNase T₁ and RNase T₁ variants during refolding under conditions where the population of enzyme-substrate or enzyme-product complexes is negligible. Recovery of enzymatic activity was continuously monitored within the visible range of the spectrum by addition of a fluorescence-labeled nucleotide substrate to the refolding sample. We found that a nonnative trans conformation at Pro39 renders the RNase T₁ almost completely inactive. A folding intermediate having a nonnative trans conformation at Pro55 shows about 46% of the enzymatic activity referred to the native state. Pro55, in contrast to the active site located Pro39, is situated in a solvent-exposed loop region remote from active-site residues. In both cases, peptidyl prolyl cis/trans isomerases accelerate the regain of nucleolytic activity. Our findings show that even if there is a considerable distance between the site of isomerization and the active site, conformational control of the bioactivity of proteins is likely to occur, and that the surface location of prolyl bonds suffices for the control of buried active sites mediated by peptidyl prolyl cis/trans isomerases.     

Proline isomerization in bovine pancreatic ribonuclease A. 2. Folding conditions

The kinetics of cis-trans isomerization of individual X-Pro peptide groups is used to study the backbone dynamics of bovine pancreatic ribonuclease A (RNase A). We previously developed and validated a fluorescence method for monitoring the cis-trans isomerization of the Tyr92-Pro93 and Asn113-Pro114 peptide groups of RNase A under unfolding conditions [Juminaga, D., Wedemeyer, W. J., and Scheraga, H. A. (1998) Biochemistry 37, 11614-11620]. The essence of this method is to introduce a fluorescent residue (Tyr or Trp) in a position adjacent to the isomerizing proline (if one is not already present) and to eliminate the fluorescence of other such residues adjacent to prolines by mutating them to phenylalanine. Here, we extend this method to observe the cis-trans isomerization of these peptide groups under folding conditions using two site-directed mutants (Y92F and Y115F) of RNase A. Both isomerizations decelerate with increasing concentrations of GdnHCl, with nearly identical m values (1.11 and 1.19 M(-1), respectively) and extrapolated zero-GdnHCl time constants (42 and 32 s, respectively); by contrast, under unfolding conditions, the cis-trans isomerizations of both Pro93 and Pro114 are independent of GdnHCl concentration. Remarkably, the isomerization rates under folding conditions at GdnHCl concentrations above 1 M are significantly slower than those measured under unfolding conditions. The temperature dependence of the Pro114 isomerization under folding conditions is also unusual; whereas Pro93 exhibits an activation energy typical of proline isomerization (19.4 kcal/mol), Pro114 exhibits a sharply reduced activation energy of 5.7 kcal/mol. A structurally plausible model accounts for these results and, in particular, shows that folding conditions strongly accelerate the cis-trans isomerization of both peptide groups to their native cis conformation, suggesting the presence of flickering local structure in their beta-hairpins.     

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.     

Cis-trans isomerization of proline in the peptide (his 105-Val 124) of ribonuclease a containing the primary nucleation site

Conformational energy calculations have been carried out to determine the relative stabilities of the C-terminal sequence 105–124 of ribonuclease A, withcis andtrans forms, respectively, of Asn 113-Pro 114. Thecis form of Pro 114 is the one that occurs in the native protein. This peptide contains the sequence 106–118, which, on the basis of both theoretical and experimental studies, is thought to constitute the primary nucleation site for the folding of ribonuclease A. It is shown that both conformations of the isolated peptide (with Pro 114 in thecis andtrans forms, respectively) are of approximately equal stability. Both forms have similar conformations from residues 105–110 and 118–124, while they differ in the bend region involving residues 111–117. Calculations have also been carried out to deduce the possible low-energy paths for the interconversion between thecis andtrans forms of both Pro 114 and Pro 117. It is shown that there are two low-energy paths (with a minimum activation energy of 16.5 kcal/mole) for the interconversion of Pro 114. Attractive nonbonded interaction energies stabilize the transition state on these paths. Only one relatively low-energy path (with an activation energy of 18 kcal/mole) could be found for the isomerization of Pro 117, which occur in thetrans form in the native protein; in this case, allcis forms have significantly higher energy than thetrans form. These calculations thus show that native-like forms for the isolated peptide can exist with Pro 114 in either thecis or thetrans form and that these forms are readily interconvertible.     

Molecular Determinants of a Native-State Prolyl Isomerization

Prolyl cis/trans isomerizations determine the rates of many protein-folding reactions, and they can serve as molecular switches and timers. The energy required to shift the prolyl cis/trans equilibrium during these processes originates from conformational reactions that are linked structurally and energetically with prolyl isomerization. We used the N2 domain of the gene-3-protein of phage fd to elucidate how such an energetic linkage develops in the course of folding. The Asp160-Pro161 bond at the tip of a β hairpin of N2 is cis in the crystal structure, but in fact, it exists as a mixture of conformers in folded N2. During refolding, about 10 kJ mol^− 1 of conformational energy becomes available for a 75-fold shift of the cis/trans equilibrium constant at Pro161, from 7/93 in the unfolded to 90/10 in the folded form. We combined single- and double-mixing kinetic experiments with a mutational analysis to identify the structural origin of this proline shift energy and to elucidate the molecular path for the transfer of this energy to Pro161. It originates largely, if not entirely, from the two-stranded β sheet at the base of the Pro161 hairpin. The two strands improve their stabilizing interactions when Pro161 is cis, and this stabilization is propagated to Pro161, because the connector peptides between the β strands and Pro161 are native-like folded when Pro161 is cis. In the presence of a trans-Pro161, the connector peptides are locally unfolded, and thus, Pro161 is structurally and energetically uncoupled from the β sheet. Such interrelations between local folding and prolyl isomerization and the potential modulation by prolyl isomerases might also be used to break and reestablish slow communication pathways in proteins.     

Conformational stability and mechanism of folding of ribonuclease T1

Urea and thermal unfolding curves for ribonuclease T1 (RNase T1) were determined by measuring several different physical properties. In all cases, steep, single-step unfolding curves were observed. When these results were analyzed by assuming a two-state folding mechanism, the plots of fraction unfolded protein versus denaturant were coincident. The dependence of the free energy of unfolding, delta G (in kcal/mol), on urea concentration is given by delta G = 5.6 - 1.21 (urea). The parameters characterizing the thermodynamics of unfolding are: midpoint of the thermal unfolding curve, Tm = 48.1 degrees C, enthalpy change at Tm, delta Hm = 97 kcal/mol, and heat capacity change, delta Cp = 1650 cal/mol deg. A single kinetic phase was observed for both the folding and unfolding of RNase T1 in the transition and post-transition regions. However, two slow kinetic phases were observed during folding in the pre-transition region. These two slow phases account for about 90% of the observed amplitude, indicating that a faster kinetic phase is also present. The slow phases probably result from cis-trans isomerization at the 2 proline residues that have a cis configuration in folded RNase T1. These results suggest that RNase T1 folds by a highly cooperative mechanism with no structural intermediates once the proline residues have assumed their correct isomeric configuration. At 25 degrees C, the folded conformation is more stable than the unfolded conformations by 5.6 kcal/mol at pH 7 and by 8.9 kcal/mol at pH 5, which is the pH of maximum stability. At pH 7, the thermodynamic data indicate that the maximum conformational stability of 8.3 kcal/mol will occur at -6 degrees C.     

Subdomain interactions as a determinant in the folding and stability of T4 lysozyme.

The folding of large, multidomain proteins involves the hierarchical assembly of individual domains. It remains unclear whether the stability and folding of small, single-domain proteins occurs through a comparable assembly of small, autonomous folding units. We have investigated the relationship between two subdomains of the protein T4 lysozyme. Thermodynamically, T4 lysozyme behaves as a cooperative unit and the unfolding transition fits a two-state model. The structure of the protein, however, resembles a dumbbell with two potential subdomains: an N-terminal subdomain (residues 13-75), and a C-terminal subdomain (residues 76-164 and 1-12). To investigate the effect of uncoupling these two subdomains within the context of the native protein, we created two circular permutations, both at the subdomain interface (residues 13 and 75). Both variants adopt an active wild-type T4 lysozyme fold. The protein starting with residue 13 is 3 kcal/mol less stable than wild type, whereas the protein beginning at residue 75 is 9 kcal/mol less stable, suggesting that the placement of the termini has a major effect on protein stability while minimally affecting the fold. When isolated as protein fragments, the C-terminal subdomain folds into a marginally stable helical structure, whereas the N-terminal subdomain is predominantly unfolded. ANS fluorescence studies indicate that, at low pH, the C-terminal subdomain adopts a loosely packed acid state. An acid state intermediate is also seen for all of the full-length variants. We propose that this acid state is comprised of an unfolded N-terminal subdomain and a loosely folded C-terminal subdomain.     

Oligomeric States along the Folding Pathways of β2-Microglobulin: Kinetics,Thermodynamics, and Structure

The transition of proteins from their soluble functional state to amyloid fibrils and aggregates is associated with the onset of several human diseases. Protein aggregation often requires some structural reshaping and the subsequent formation of intermolecular contacts. Therefore, the study of the conformation of excited protein states and their ability to form oligomers is of primary importance for understanding the molecular basis of amyloid fibril formation. Here, we investigated the oligomerization processes that occur along the folding of the amyloidogenic human protein β2-microglobulin. The combination of real-time two-dimensional NMR data with real-time small-angle X-ray scattering measurements allowed us to derive thermodynamic and kinetic information on protein oligomerization of different conformational states populated along the folding pathways. In particular, we could demonstrate that a long-lived folding intermediate (I-state) has a higher propensity to oligomerize compared to the native state. Our data agree well with a simple five-state kinetic model that involves only monomeric and dimeric species. The dimers have an elongated shape with the dimerization interface located at the apical side of β2-microglobulin close to Pro32, the residue that has a trans conformation in the I-state and a cis conformation in the native (N) state. Our experimental data suggest that partial unfolding in the apical half of the protein close to Pro32 leads to an excited state conformation with enhanced propensity for oligomerization. This excited state becomes more populated in the transient I-state due to the destabilization of the native conformation by the trans-Pro32 configuration.     

Investigation of the folding pathway of the TEM-1 β-lactamase

The TEM-1 β-lactamase is a globular protein containing 12 proline residues. The folding mechanism of this enzyme was investigated by kinetic and equilibrium experiments with the help of fluorescence spectroscopy and circular dichroism. The equilibrium denaturation of the protein induced by guanidine hydrochloride occurs in two discrete steps, indicating the existence of a thermodynamically stable intermediate state. Thisstate is 5.2 ± 0.4 kcal/mol less stable than the native conformation and 5.7 ± 0.2 kcal/mol more stable than the fully denaturedprotein. This intermediate state exhibits a high content of native secondary structure elements but is devoid of specific tertiary organization; its relation to the “molten globule” is discussed. Refolding kinetic experimentsrevealed the existence of a transient intermediate conformation between thethermodynamically stable intermediate and the native protein. This transient intermediate appears rapidly during the folding reaction. It exhibits a secondary structure content very similar to that of the native protein and has also recovered a significant amount of tertiary organisation. The final refolding step of the TEM-1 β-lactamase, leading to the native enzyme, is dominated by two major slow kinetic phases which probablyreflect a very complex process kinetically limited by proline cis/transisomerization. © 1995 Wiley-Liss, Inc.     

Multiple roles of prolyl residues in structure and folding

To explore the ways that proline residues may influence the conformational options of a polypeptide backbone, we have characterized Pro-->Ala mutants of cellular retinoic acid-binding protein I (CRABP I). While all three Xaa-Pro bonds are in the trans conformation in the native protein and the equilibrium stability of each mutant is similar to that of the parent protein, each has distinct effects on folding and unfolding kinetics. The mutation of Pro105 does not alter the kinetics of folding of CRABP I, which indicates that the flexible loop containing this residue is passive in the folding process. By contrast, replacement of Pro85 by Ala abolishes the observable slow phase of folding, revealing that correct configuration of the 84-85 peptide bond is prerequisite to productive folding. Substitution of Pro39 by Ala yields a protein that folds and unfolds more slowly. Removal of the conformational constraint imposed by the proline ring likely raises the transition state barrier by increasing the entropic cost of narrowing the conformational ensemble. Additionally, the Pro-->Ala mutation removes a helix-termination signal that is important for efficient folding to the native state.     

The role of proline residues in the folding kinetics of the bovine pancreatic trypsin inhibitor derivative RCAM(14–38)

The role of proline residues in the folding of the trypsin inhibitor derivative RCAM(14–38) has been studied by testing for slow-folding species of the unfolded protein, which could result from the introduction of wrong proline isomers after unfolding. The unfolded protein at 25 °C contains chiefly fast-folding (U_F) molecules: they refold with a time constant of 40 milliseconds at pH 6.8 in 1.9 m-guanidinium chloride. At least one minor slow-folding (U_s) species has been found, using fluorescence to monitor refolding. The reaction in which this U_s species is formed after unfolding shows the properties expected for the cis: Irans isomerization of a proline residue. When refolding is monitored by tyrosine absorbance, two minor slow reactions are found. The faster reaction is in the same time range (15 s at 25 °C) as that studied by fluorescence, and the slower reaction is quite slow (200 s at 25 °C). It is not known whether the slower reaction results from a second U_s species. There are four trans proline residues in bovine pancreatic trypsin inhibitor: the proportion of slow-folding molecules (not more than 25% at 25 °C) is smaller than expected if every proline residue can produce a U_s species and if the cis to trans ratio of each residue after unfolding is at least 0.1:0.9.Criteria based on folding kinetics are given for classifying the types of folding reaction shown by unfolded molecules containing a single wrong proline isomer. Levitt (1980) has classified three types of proline residues according to the energy difference (small, intermediate or large) between the native protein and the predicted minimum energy structure containing a wrong proline isomer. He suggests that these three types of proline residues can be recognized by the types of folding reactions they produce. Only type II (intermediate) folding reactions have thus far been characterized by the criteria introduced here. We point out that the type of folding reaction depends also on the folding conditions, and a possible explanation for this effect is given.     

Elimination of a cis-Proline-Containing Loop and Turn Optimization Stabilizes a Protein and Accelerates Its Folding

In the N2 domain of the gene-3-protein of phage fd, two consecutive β-strands are connected by a mobile loop of seven residues (157-163). The stability of this loop is low, and the Asp160-Pro161 bond at its tip shows conformational heterogeneity with 90% being in the cis and 10% in the trans form. The refolding kinetics of N2 are complex because the molecules with cis or trans isomers at Pro161 both fold to native-like conformations, albeit with different rates. We employed consensus design to shorten the seven-residue irregular loop around Pro161 to a four-residue type I′ turn without a proline. This increased the conformational stability of N2 by almost 10 kJ mol^− 1 and abolished the complexity of the folding kinetics. Turn sequences obtained from in vitro selections for increased stability strongly resembled those derived from the consensus design. Two other type I′ turns of N2 could also be stabilized by consensus design. For all three turns, the gain in stability originates from an increase in the rate of refolding. The turns form native-like structures early during refolding and thus stabilize the folding transition state. The crystal structure of the variant with all three stabilized turns confirms that the 157-163 loop was in fact shortened to a type I′ turn and that the other turns maintained their type I′ conformation after sequence optimization.     

Selective inactivation of rabbit muscle glycerophosphate dehydrogenase by diazotized 3-aminopyridine adenine dinucleotide.

Incorporation of a proline analog into collagen polypeptides was studied by incubating matrix-free tendon cells from 17-day-old chick embryos with cis-4-hydroxy-l-proline. Velocity sedimentation of intracellular polypeptides provided further evidence that incorporation of the analog into protein prevented the pro-α- and pro-γ-chains of procollagen from folding into a stable triple-helical conformation. The size of the newly synthesized intracellular and extracellular protein was examined under conditions which prevented proteolysis during processing of the samples. In contrast to previous observations, the results demonstrated that there was little if any intracellular degradation of nonhelical pro-α-and pro-γ-chains containing the proline analog, and a fraction of the nonhelical pro-γ-chains was secreted into the medium without extensive degradation. In further studies, the cells were incubated with ¹⁴C lysine, and the synthesis of glycosylated hydroxylysyl residues was measured in control cells and in cells incubated with cis-4-hydroxy-l-proline. The results demonstrated that the content of glycosylated hydroxylysyl residues in nonhelical pro-γ-chains containing cis-4-hydroxy-l-proline was increased twofold as compared to the triple-helical procollagen in control cells. The results suggested that under control conditions folding into the triplehelical conformation limits the extent of glycosylation of collagen. If folding is prevented or delayed, procollagen polypeptides are more extensively glycosylated.     

Local control of peptide conformation: Stabilization of cis proline peptide bonds by aromatic proline interactions

In the native state of proteins there is a marked tendency for an aromatic amino acid to precede a cis proline. There are also significant differences between the three aromatic amino acids with Tyr exhibiting a noticeably higher propensity than Phe or Trp to precede a cis proline residue. In order to study the role that local interactions play in these conformation preferences, a set of tetrapeptides of the general sequence acetyl-Gly-X-Pro-Gly-carboxamide (GXPG), where X = Tyr, Phe, Trp, Ala, or cyclohexyl alanine, were synthesized and studied by nmr. Analysis of the nmr data shows that none of the peptides adopt a specific backbone structure. Ring current shifts, the equilibrium constants, the Van't Hoff enthalpy, and the measured rate of cis-trans isomerization all indicate that the cis proline conformer is stabilized by favorable interactions between the aromatic ring and the proline residue. Analysis of the side chain conformation of the aromatic residue and analysis of the chemical shifts of the pyrrolidine ring protons shows that the aromatic side chain adopts a preferred conformation in the cis form. The distribution of rotamers and the effect of an aromatic residue on the cis-trans equilibrium indicate that the preferred conformation is populated to approximately 62% for the Phe containing peptide, 67% for the Tyr containing peptide, and between 75 and 80% for the Trp containing peptide. The interaction is unaffected by the addition of 8M urea. These local interactions favor an aromatic residue immediately preceding a cis proline, but they cannot explain the relative propensities for Phe-Pro, Tyr-Pro, and Trp-Pro cis peptide bonds observed in the native state of proteins. In the model peptides the percentage of the cis proline conformer is 21% GYPG while it is 17% for GFPG. This difference is considerably smaller than the almost three to one preponderance observed for cis Tyr-Pro peptide bonds vs cis Phe-Pro peptide bonds in the protein database. © 1998 John Wiley & Sons, Inc. Biopoly 45: 381–394, 1998     

Empirical evaluation of the influence of side chains on the conformational entropy of the polypeptide backbone

Changes in amino acid side chains have long been recognized to alterthe range and distribution of ?, ψ angles found in the main chain of polypeptides. Altering the range and distribution of ?, ψ angles also alters the conformational entropy of the flexible denatured state and may thus stabilize or destabilize it relative to the comparatively conformationally rigid native state. A database of 12,320 residues from 61 nonhomologous, high resolution crystal structures was examined to determine the ?, ψ conformational preferences of each of the 20 amino acids. These observed distributions in the native state of proteins are assumed to also reflect the distributions found in the denatured state. The distributionswere used to approximate the energy surface for each residue, allowing the calculation of relative conformational entropies for each residue relative to glycine. In the most extreme case, replacement of glycine by proline, conformational entropy changes will stabilize the native state relative to the denatured state by ?0.82 ± 0.08 kcal/mol at 20°C. Surprisingly, alanine is found to be the most ordered residue other than proline. This unexpected result is a result of the high percentage of alanines found in helical conformations. This either indicates that the observed distributions in the native state do not reflect the distributions in the denatured state, or that alanine is much more likely to adopt a helical conformation in the denatured state than residues with longer side chains. Among those residues with ?, ψ angles compatible with helix incorporation the percentage of alanines actually in helices is very similar to other residues. This and the consistent ordering of alanine relative to other residues regardless of secondary structure are evidence that ?, ψ distributions in native states reflect those in the denatured states. © 1995 Wiley-Liss, Inc.     

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.     

Proline in alpha-helix: stability and conformation studied by dynamics simulation

Free-energy simulations have been used to estimate the change in the conformational stability of short polyalanine alpha-helices when one of the alanines is replaced by a proline residue. For substituting proline in the middle of the helix the change in free energy of folding (delta delta G degrees) was calculated as 14 kJ/mol (3.4 kcal/mol), in excellent agreement with the one available experimental value. The helix containing proline was found to be strongly kinked; the free energy for reducing the angle of the kink from 40 degrees to 15 degrees was calculated, and found to be small. A tendency to alternate hydrogen bonding schemes was observed in the proline-containing helix. These observations for the oligopeptide agree well with the observation of a range of kink angles (18-35 degrees) and variety of hydrogen bonding schemes, in the rare instances where proline occurs in helices in globular proteins. For substituting proline at the N-terminus of the helix the change in free energy of folding (delta delta G degrees) was calculated as -4 kJ/mol in the first helical position (N1) and +6 kJ/mol in the second helical position (N2). The observed frequent occurrence of proline in position N1 in alpha-helices in proteins therefore has its origin in stability differences of secondary structure. The conclusion reached here that proline may be a better helix former in position N1 than (even) alanine, and thus be a helix initiator may be testable experimentally by measurements of fraction helical conformation of individual residues in oligopeptides of appropriate sequence. The relevance of these results in regards to the frequent occurrence of proline-containing helices in certain membrane proteins is discussed.     

Cis proline mutants of ribonuclease A. I. Thermal stability.

A chemically synthesized gene for ribonuclease A has been expressed in Escherichia coli using a T7 expression system (Studier, F.W., Rosenberg, A.H., Dunn, J.J., & Dubendorff, J.W., 1990, Methods Enzymol. 185, 60-89). The expressed protein, which contains an additional N-terminal methionine residue, has physical and catalytic properties close to those of bovine ribonuclease A. The expressed protein accumulates in inclusion bodies and has scrambled disulfide bonds; the native disulfide bonds are regenerated during purification. Site-directed mutations have been made at each of the two cis proline residues, 93 and 114, and a double mutant has been made. In contrast to results reported for replacement of trans proline residues, replacement of either cis proline is strongly destabilizing. Thermal unfolding experiments on four single mutants give delta Tm approximately equal to 10 degrees C and delta delta G0 (apparent) = 2-3 kcal/mol. The reason is that either the substituted amino acid goes in cis, and cis<==>trans isomerization after unfolding pulls the unfolding equilibrium toward the unfolded state, or else there is a conformational change, which by itself is destabilizing relative to the wild-type conformation, that allows the substituted amino acid to form a trans peptide bond.     

Synthetic biology of proteins: tuning GFPs folding and stability with fluoroproline

Background

Proline residues affect protein folding and stability via cis/trans isomerization of peptide bonds and by the C^γ-exo or -endo puckering of their pyrrolidine rings. Peptide bond conformation as well as puckering propensity can be manipulated by proper choice of ring substituents, e.g. C^γ-fluorination. Synthetic chemistry has routinely exploited ring-substituted proline analogs in order to change, modulate or control folding and stability of peptides.

Methodology/Principal Findings

In order to transmit this synthetic strategy to complex proteins, the ten proline residues of enhanced green fluorescent protein (EGFP) were globally replaced by (4R)- and (4S)-fluoroprolines (FPro). By this approach, we expected to affect the cis/trans peptidyl-proline bond isomerization and pyrrolidine ring puckering, which are responsible for the slow folding of this protein. Expression of both protein variants occurred at levels comparable to the parent protein, but the (4R)-FPro-EGFP resulted in irreversibly unfolded inclusion bodies, whereas the (4S)-FPro-EGFP led to a soluble fluorescent protein. Upon thermal denaturation, refolding of this variant occurs at significantly higher rates than the parent EGFP. Comparative inspection of the X-ray structures of EGFP and (4S)-FPro-EGFP allowed to correlate the significantly improved refolding with the C^γ-endo puckering of the pyrrolidine rings, which is favored by 4S-fluorination, and to lesser extents with the cis/trans isomerization of the prolines.

Conclusions/Significance

We discovered that the folding rates and stability of GFP are affected to a lesser extent by cis/trans isomerization of the proline bonds than by the puckering of pyrrolidine rings. In the C^γ-endo conformation the fluorine atoms are positioned in the structural context of the GFP such that a network of favorable local interactions is established. From these results the combined use of synthetic amino acids along with detailed structural knowledge and existing protein engineering methods can be envisioned as a promising strategy for the design of complex tailor-made proteins and even cellular structures of superior properties compared to the native forms.