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.     

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.     

Kinetic coupling of folding and prolyl isomerization of beta2-microglobulin studied by mutational analysis

β₂-Microglobulin (β2-m), a protein responsible for dialysis-related amyloidosis, adopts a typical immunoglobulin domain fold with the N-terminal peptide bond of Pro32 in a cis isomer. The refolding of β2-m is limited by the slow trans-to-cis isomerization of Pro32, implying that intermediates with a non-native trans-Pro32 isomer are precursors for the formation of amyloid fibrils. To obtain further insight into the Pro-limited folding of β2-m, we studied the Gdn-HCl-dependent unfolding/refolding kinetics using two mutants (W39 and P32V β2-ms) as well as the wild-type β2-m. W39 β2-m is a triple mutant in which both of the authentic Trp residues (Trp60 and Trp95) are replaced by Phe and a buried Trp common to other immunoglobulin domains is introduced at the position of Leu39 (i.e., L39W/W60F/W95F). W39 β2-m exhibits a dramatic quenching of fluorescence upon folding, enabling a detailed analysis of Pro-limited unfolding/refolding. On the other hand, P32V β2-m is a mutant in which Pro32 is replaced by Val, useful for probing the kinetic role of the trans-to-cis isomerization of Pro32. A comparative analysis of the unfolding/refolding kinetics of these mutants including three types of double-jump experiments revealed the prolyl isomerization to be coupled with the conformational transitions, leading to apparently unusual kinetics, particularly for the unfolding. We suggest that careful consideration of the kinetic coupling of unfolding/refolding and prolyl isomerization, which has tended to be neglected in recent studies, is essential for clarifying the mechanism of protein folding and, moreover, its biological significance.     

Energetic coupling between native-state prolyl isomerization and conformational protein folding

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.     

Energetic Communication between Functional Sites of the Gene-3-Protein during Infection by Phage fd

To initiate infection of Escherichia coli, phage fd uses its gene-3-protein (G3P) to bind first to an F pilus and then to the TolA protein at the cell surface. G3P is normally auto-inhibited because a tight interaction between the two N-terminal domains N1 and N2 buries the TolA binding site. Binding of N2 to the pilus activates G3P by initiating long-range conformational changes that are relayed to the domain interface and to a proline timer. We discovered that the 23–28 loop of the N1 domain is critical for propagating these conformational signals. The analysis of the stability and the folding dynamics of G3P variants with a shortened loop combined with TolA interaction studies and phage infection experiments reveal how the contact between the N2 domain and the 23–28 loop of N1 is energetically linked with the interdomain region and the proline timer and how it affects phage infectivity. Our results illustrate how conformational transitions and prolyl cis/trans isomerization can be coupled energetically and how conformational signals to and from prolines can be propagated over long distances in proteins.     

Native-State Heterogeneity of β2-Microglobulin as Revealed by Kinetic Folding and Real-Time NMR Experiments

The kinetic folding of β₂-microglobulin from the acid-denatured state was investigated by interrupted-unfolding and interrupted-refolding experiments using stopped-flow double-jump techniques. In the interrupted unfolding, we first unfolded the protein by a pH jump from pH 7.5 to pH 2.0, and the kinetic refolding assay was carried out by the reverse pH jump by monitoring tryptophan fluorescence. Similarly, in the interrupted refolding, we first refolded the protein by a pH jump from pH 2.0 to pH 7.5 and used a guanidine hydrochloride (GdnHCl) concentration jump as well as the reverse pH jump as unfolding assays. Based on these experiments, the folding is represented by a parallel-pathway model, in which the molecule with the correct Pro32 cis isomer refolds rapidly with a rate constant of 5–6 s^? 1, while the molecule with the Pro32 trans isomer refolds more slowly (pH 7.5 and 25 °C). At the last step of folding, the native-like trans conformer produced on the latter pathway isomerizes very slowly (0.001–0.002 s^? 1) into the native cis conformer. In the GdnHCl-induced unfolding assays in the interrupted refolding, the native-like trans conformer unfolded remarkably faster than the native cis conformer, and the direct GdnHCl-induced unfolding was also biphasic, indicating that the native-like trans conformer is populated at a significant level under the native condition. The one-dimensional NMR and the real-time NMR experiments of refolding further indicated that the population of the trans conformer increases up to 7–9% under a more physiological condition (pH 7.5 and 37 °C).     

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.     

Effect of proline residues on protein folding

Conformational energy calculations have been used to study the role of the proline residues in the folding of bovine pancreatic trypsin inhibitor. In the calculation, each of the four proline residues of this small protein is forced from the trans to cis peptide isomer while still part of the native folded structure. The cis proline residue can always be accommodated by small changes of the native conformation (< 1 Å root-mean-square deviation). For three of the four proline residues, Pro2, Pro9 and Pro 13, being in the cis form is calculated to destabilize the folded conformation by less than 11 kcal/mol, suggesting that rapid folding to a stable native-like conformation can occur with either isomeric form. For one of these three, Pro13, the destabilization is only 1 kcal/mol, suggesting the existence of an alternative folded native conformation with Pro13 cis. The fourth proline residue, Pro8, is calculated to destabilize the native conformation by so much (33 kcal/mol) that it will block folding in the manner proposed by Brandts et al. (1975).     

Coevolution of Antibody Stability and Vκ CDR-L3 Canonical Structure

Antibodies recognize antigens through six hypervariable loops, five of which have a limited set of conformations known as canonical structures. For κ light chains, the majority of CDR-L3 [the third hypervariable loop of the light chain variable domain (V_L)] adopts the type 1 canonical structure (CS1), with a cis-proline at position 95. Here, we present the design and structural studies of the monoclonal antibody mAb15 and related mutants that contained a series of progressively germline mutations only in the heavy chain variable domain (V_H) that ultimately led to an increase of more than 11 °C in the melting temperature (T_m) of the antigen-binding fragment (Fab). The all-trans CDR-L3 structure in the wild type is significantly different from any known CDR-L3 canonical structures. In the thermally stable mutants, the L94^L-S95^L peptide bond adopts an energetically unfavorable non-X-proline cis conformation, but the overall CDR-L3 loop converted to CS1. The stabilized V_H appears to function as a specific molecular chaperone that facilitated the trans-cis isomerization of S95^L. Thus, it is plausible that proline is the evolutionary choice to maintain overall structure and stability for V_L. These results provide new insights into the evolution of CS1 and suggest a potential molecular switch mechanism at position 95 that links CDR-L3 structural diversity and antibody stability and will have implications for antibody engineering.     

Exclusion of the native alpha-helix from the amyloid fibrils of a mixed alpha/beta protein

Members of the cystatin superfamily are involved in an inherited form of cerebral amyloid angiopathy and readily form amyloid fibrils in vitro. We have determined the structured core of human stefin B (cystatin B) amyloid fibrils using quenched hydrogen exchange and NMR. The core contains residues from four of the five strands of the native β-sheet, delimited by unprotected loop regions analogous to those of the native monomeric structure. However, non-native features are also apparent, the most striking of which is the exclusion of the native α-helix. Before forming amyloid in vitro, cystatins dimerise via 3D domain swapping, and assemble into tetramers with trans to cis isomerism of a conserved proline. In the fibril, the hinge loop that forms an extended β-structure in the dimer remains protected, consistent with the domain-swapping interface being maintained. However, the fibril data are not compatible with a simple 3D domain-swapping model for amyloid formation, and the displacement of the helix points to alternative packing arrangements of native-like β-structure, in which proline isomerism is important in preventing steric clashing.     

Real-time NMR studies on a transient folding intermediate of barstar.

The refolding of barstar, the intracellular inhibitor of barnase, is dominated by the slow formation of a cis peptidyl prolyl bond in the native protein. The triple mutant C40/82A P27A in which two cysteine residues and one trans proline were replaced by alanine was used as model system to investigate the kinetics and structural consequences of the trans/cis interconversion of Pro48. One- and two-dimensional real-time NMR spectroscopy was used to follow the trans/cis interconversion after folding was initiated by rapid dilution of the urea denatured protein. Series of 1H, 15N HSQC spectra acquired with and without the addition of peptidyl prolyl isomerase unambiguously revealed the accumulation of a transient trans-Pro48 intermediate within the dead time of the experiment. Subtle chemical shift differences between the native state and the intermediate spectra indicate that the intermediate is predominantly native-like with a local rearrangement in the Pro48 loop and in the beta-sheet region including residues Tyr47, Ala82, Thr85, and Val50.     

Structure of H/ACA RNP protein Nhp2p reveals cis/trans isomerization of a conserved proline at the RNA and Nop10 binding interface

H/ACA small nucleolar and Cajal body ribonucleoproteins (RNPs) function in site-specific pseudouridylation of eukaryotic rRNA and snRNA, rRNA processing, and vertebrate telomerase biogenesis. Nhp2, one of four essential protein components of eukaryotic H/ACA RNPs, forms a core trimer with the pseudouridylase Cbf5 and Nop10 that binds to H/ACA RNAs specifically. Crystal structures of archaeal H/ACA RNPs have revealed how the protein components interact with each other and with the H/ACA RNA. However, in place of Nhp2p, archaeal H/ACA RNPs contain L7Ae, which binds specifically to an RNA K-loop motif absent from eukaryotic H/ACA RNPs, while Nhp2 binds a broader range of RNA structures. We report solution NMR studies of Saccharomyces cerevisiae Nhp2 (Nhp2p), which reveal that Nhp2p exhibits two major conformations in solution due to cis/trans isomerization of the evolutionarily conserved Pro83. The equivalent proline is in the cis conformation in all reported structures of L7Ae and other homologous proteins. Nhp2p has the expected α-β-α fold, but the solution structures of the major conformation of Nhp2p with trans Pro83 and of Nhp2p-S82W with cis Pro83 reveal that Pro83 cis/trans isomerization affects the positions of numerous residues at the Nop10 and RNA binding interface. An S82W substitution, which stabilizes the cis conformation, also stabilizes the association of Nhp2p with H/ACA snoRNPs expressed in vivo. We propose that Pro83 plays a key role in the assembly of the eukaryotic H/ACA RNP, with the cis conformation locking in a stable Cbf5-Nop10-Nhp2 ternary complex and positioning the protein backbone to interact with the H/ACA RNA.     

Consequences of Domain Insertion on the Stability and Folding Mechanism of a Protein

SlyD, the sensitive-to-lysis protein from Escherichia coli, consists of two domains. They are not arranged successively along the protein chain, but one domain, the “insert-in-flap” (IF) domain, is inserted internally as a guest into a surface loop of the host domain, which is a prolyl isomerase of the FK506 binding protein (FKBP) type. We used SlyD as a model to elucidate how such a domain insertion affects the stability and folding mechanism of the host and the guest domain. For these studies, the two-domain protein was compared with a single-domain variant SlyDΔIF, SlyD* without the chaperone domain (residues 1-69 and 130-165) in which the IF domain was removed and replaced by a short loop, as present in human FKBP12. Equilibrium unfolding and folding kinetics followed an apparent two-state mechanism in the absence and in the presence of the IF domain. The inserted domain decreased, however, the stability of the host domain in the transition region and decelerated its refolding reaction by about 10-fold. This originates from the interruption of the chain connectivity by the IF domain and its inherent instability. To monitor folding processes in this domain selectively, a Trp residue was introduced as fluorescent probe. Kinetic double-mixing experiments revealed that, in intact SlyD, the IF domain folds and unfolds about 1000-fold more rapidly than the FKBP domain, and that it is strongly stabilized when linked with the folded FKBP domain. The unfolding limbs of the kinetic chevrons of SlyD show a strong downward curvature. This deviation from linearity is not caused by a transition-state movement, as often assumed, but by the accumulation of a silent unfolding intermediate at high denaturant concentrations. In this kinetic intermediate, the FKBP domain is still folded, whereas the IF domain is already unfolded.     

Kinetic Control in Protein Folding for Light Chain Amyloidosis and the Differential Effects of Somatic Mutations

Light chain amyloidosis is a devastating disease where immunoglobulin light chains form amyloid fibrils, resulting in organ dysfunction and death. Previous studies have shown a direct correlation between the protein thermodynamic stability and the propensity for amyloid formation for some proteins involved in light chain amyloidosis. Here we investigate the effect of somatic mutations on protein stability and in vitro fibril formation of single and double restorative mutants of the protein AL-103 compared to the wild-type germline control protein. A scan rate dependence and hysteresis in the thermal unfolding and refolding was observed for all proteins. This indicates that the unfolding/refolding reaction is kinetically determined with different kinetic constants for unfolding and refolding even though the process remains experimentally reversible. Our structural analysis of AL-103 and AL-103 delP95aIns suggests a kinetic coupling of the unfolding/refolding process with cis–trans prolyl isomerization. Our data reveal that the deletion of proline 95a (AL-103 delP95aIns), which removes the trans–cis di-proline motif present in the patient protein AL-103, results in a dramatic increment in the thermodynamic stability and a significant delay in fibril formation kinetics with respect to AL-103. Fibril formation is pH dependent; all proteins form fibrils at pH 2; reactions become slower and more stochastic as the pH increases up to pH 7. Based on these results, we propose that, in addition to thermodynamic stability, kinetic stability (possibly influenced by the presence of cis proline 95a) plays a major role in the AL-103 amyloid fibril formation process.     

OmpA can form folded and unfolded oligomers

The monomeric outer membrane protein OmpA from Escherichia coli has long served as a model protein for studying the folding and membrane insertion of β-barrel membrane proteins. Here we report that when OmpA is refolded in limiting amounts of surfactant (close to the cmc), it has a high propensity to form folded and unfolded oligomers. The oligomers exist both in a folded and (partially) unfolded form which both dissociate under denaturing conditions. Oligomerization does not require the involvement of the periplasmic domain and is not strongly affected by ionic strength. The folded dimers can be isolated and show native-like secondary structure; they are resistant to proteolytic attack and do not dissociate in high surfactant concentrations, indicating high kinetic stability once formed. Remarkably, OmpA also forms significant amounts of higher order structures when refolding in the presence of lipid vesicles. We suggest that oligomerization occurs by domain swapping favored by the high local concentration of OmpA molecules congregating on the same micelle or vesicle. In this model, the unfolded oligomer is stabilized by a small number of intermolecular β-strand contacts and subsequently folds to a more stable state where these intermolecular contacts are consolidated in a native-like fashion by contacts between complementary β-strands from different molecules. Our model is supported by the ability of complementary fragments to associate with each other in vitro. Oligomerization is probably avoided in the cell by the presence of cellular chaperones which maintain the protein in a monomeric state.     

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.     

Structure of a Double Transmembrane Fragment of a G-Protein-Coupled Receptor in Micelles

The structure and dynamic properties of an 80-residue fragment of Ste2p, the G-protein-coupled receptor for α-factor of Saccharomyces cerevisiae, was studied in LPPG micelles with the use of solution NMR spectroscopy. The fragment Ste2p(G31-T110) (TM1-TM2) consisted of 19 residues from the N-terminal domain, the first TM helix (TM1), the first cytoplasmic loop, the second TM helix (TM2), and seven residues from the first extracellular loop. Multidimensional NMR experiments on [¹⁵N], [¹⁵N, ¹³C], [¹⁵N, ¹³C, ²H]-labeled TM1-TM2 and on protein fragments selectively labeled at specific amino acid residues or protonated at selected methyl groups resulted in >95% assignment of backbone and side-chain nuclei. The NMR investigation revealed the secondary structure of specific residues of TM1-TM2. TALOS constraints and NOE connectivities were used to calculate a structure for TM1-TM2 that was highlighted by the presence of three α-helices encompassing residues 39-47, 49-72, and 80-103, with higher flexibility around the internal Arg⁵⁸ site of TM1. RMSD values of individually superimposed helical segments 39-47, 49-72, and 80-103 were 0.25 ± 0.10 Å, 0.40 ± 0.13 Å, and 0.57 ± 0.19 Å, respectively. Several long-range interhelical connectivities supported the folding of TM1-TM2 into a tertiary structure typified by a crossed helix that splays apart toward the extracellular regions and contains considerable flexibility in the G⁵⁶VRSG⁶⁰ region. ¹⁵N-relaxation and hydrogen-deuterium exchange data support a stable fold for the TM parts of TM1-TM2, whereas the solvent-exposed segments are more flexible. The NMR structure is consistent with the results of biochemical experiments that identified the ligand-binding site within this region of the receptor.     

Crystal Structure Determination and Functional Characterization of the Metallochaperone SlyD from Thermus thermophilus

SlyD (sensitive to lysis D; product of the slyD gene) is a prolyl isomerase [peptidyl-prolyl cis/trans isomerase (PPIase)] of the FK506 binding protein (FKBP) type with chaperone properties. X-ray structures derived from three different crystal forms reveal that SlyD from Thermus thermophilus consists of two domains representing two functional units. PPIase activity is located in a typical FKBP domain, whereas chaperone function is associated with the autonomously folded insert-in-flap (IF) domain. The two isolated domains are stable and functional in solution, but the presence of the IF domain increases the PPIase catalytic efficiency of the FKBP domain by 2 orders of magnitude, suggesting that the two domains act synergistically to assist the folding of polypeptide chains. The substrate binding surface of SlyD from T. thermophilus was mapped by NMR chemical shift perturbations to hydrophobic residues of the IF domain, which exhibits significantly reduced thermodynamic stability according to NMR hydrogen/deuterium exchange and fluorescence equilibrium transition experiments. Based on structural homologies, we hypothesize that this is due to the absence of a stabilizing β-strand, suggesting in turn a mechanism for chaperone activity by ‘donor-strand complementation.’ Furthermore, we identified a conserved metal (Ni²⁺) binding site at the C-terminal SlyD-specific helical appendix of the FKBP domain, which may play a role in metalloprotein assembly.     

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.     

The equilibrium unfolding intermediate observed at pH 4 and its relationship with the kinetic folding intermediates in green fluorescent protein

Folding mechanisms of a variant of green fluorescent protein (F99S/M153T/V163A) were investigated by a wide variety of spectroscopic techniques. Equilibrium measurements on acid-induced denaturation of the protein monitored by chromophore and tryptophan fluorescence and small-angle X-ray scattering revealed that this protein accumulates at least two equilibrium intermediates, a native-like intermediate and an unfolding intermediate, the latter of which exhibits the characteristics of the molten globule state under moderately denaturing conditions at pH 4. To elucidate the role of the equilibrium unfolding intermediate in folding, a series of kinetic refolding experiments with various combinations of initial and final pH values, including pH 7.5 (the native condition), pH 4.0 (the moderately denaturing condition where the unfolding intermediate is accumulated), and pH 2.0 (the acid-denaturing condition) were carried out by monitoring chromophore and tryptophan fluorescence. Kinetic on-pathway intermediates were accumulated during the folding on the refolding reaction from pH 2.0 to 7.5. However, the signal change corresponding to the conversion from the acid-denatured to the kinetic intermediate states was significantly reduced on the refolding reaction from pH 4.0 to pH 7.5, whereas only the signal change corresponding to the above conversion was observed on the refolding reaction from pH 2.0 to pH 4.0. These results indicate that the equilibrium unfolding intermediate is composed of an ensemble of the folding intermediate species accumulated during the folding reaction, and thus support a hierarchical model of protein folding.