Context-dependent effects of proline residues on the stability and folding pathway of ubiquitin

Overview of protein folding mechanisms: experimental and theoretical approaches to probing energy landscapes

We present an overview of the current experimental and theoretical approaches to studying protein folding mechanisms, set against current models of the folding energy landscape. We describe how stability and folding kinetics can be determined experimentally and how this data can be interpreted in terms of the characteristic features of various models from the simplest two-state pathway to a multi-state mechanism. We summarize the pros and cons of a range of spectroscopic methods for measuring folding rates and present a theoretical framework, coupled with protein engineering approaches, for elucidating folding mechanisms and structural features of folding transition states. A series of case studies are used to show how experimental kinetic data can be interpreted in the context of non-native interactions, populated intermediates, parallel folding pathways, and sequential transition states. We also show how computational methods now allow transient species of high energy, such as folding transition states, to be modeled on the basis of experimental Φ-value analysis derived from the effects of point mutations on folding kinetics.

Rational design of protein stability: effect of (2S,4R)-4-fluoroproline on the stability and folding pathway of ubiquitin

BACKGROUND

Many strategies have been employed to increase the conformational stability of proteins. The use of 4-substituted proline analogs capable to induce pre-organization in target proteins is an attractive tool to deliver an additional conformational stability without perturbing the overall protein structure. Both, peptides and proteins containing 4-fluorinated proline derivatives can be stabilized by forcing the pyrrolidine ring in its favored puckering conformation. The fluorinated pyrrolidine rings of proline can preferably stabilize either a C(γ)-exo or a C(γ)-endo ring pucker in dependence of proline chirality (4R/4S) in a complex protein structure. To examine whether this rational strategy can be generally used for protein stabilization, we have chosen human ubiquitin as a model protein which contains three proline residues displaying C(γ)-exo puckering.

METHODOLOGY/PRINCIPAL FINDINGS

While (2S,4R)-4-fluoroproline ((4R)-FPro) containing ubiquitinin can be expressed in related auxotrophic Escherichia coli strain, all attempts to incorporate (2S,4S)-4-fluoroproline ((4S)-FPro) failed. Our results indicate that (4R)-FPro is favoring the C(γ)-exo conformation present in the wild type structure and stabilizes the protein structure due to a pre-organization effect. This was confirmed by thermal and guanidinium chloride-induced denaturation profile analyses, where we observed an increase in stability of -4.71 kJ·mol(-1) in the case of (4R)-FPro containing ubiquitin ((4R)-FPro-ub) compared to wild type ubiquitin (wt-ub). Expectedly, activity assays revealed that (4R)-FPro-ub retained the full biological activity compared to wt-ub.

CONCLUSIONS/SIGNIFICANCE

The results fully confirm the general applicability of incorporating fluoroproline derivatives for improving protein stability. In general, a rational design strategy that enforces the natural occurring proline puckering conformation can be used to stabilize the desired target protein.

The action of fish peptide Orpotrin analogs on microcirculation

In order to investigate the relationship between the primary structure of Orpotrin, a vasoactive peptide previously isolated from the freshwater stingray Potamotrygon gr. orbignyi, and its microcirculatory effects, three Orpotrin analogs were synthesized. The analogs have a truncated N-terminal with a His residue deletion and two substituted amino acid residues, where one Nle is substituted for one internal Lys residue and the third analog has a substitution of a Pro for an Ala (Orp-desH(1) , Orp-Nle and Orp-Pro/Ala, respectively). Only Orp-desH(1) could induce a lower vasoconstriction effect compared with the natural Orpotrin, indicating that besides the N-terminal, the positive charge of Lys and the Pro residues located at the center of the amino acid chain is crucial for this vasoconstriction effect. Importantly, the suggestions made with bioactive peptides were based on the molecular modeling and dynamics of peptides, the presence of key amino acids and shared activity in microcirculation, characterized by intravital microscopy. Moreover, this study has demonstrated that even subtle changes in the primary structure of Orpotrin alter the biological effects of this native peptide significantly, which could be of interest for biotechnological applications.

Mutational investigation of protein folding transition states by Phi-value analysis and beyond: lessons from SH3 domain folding

Understanding how proteins adopt their unique native structures requires a complete structural characterization of the rate-limiting transition state(s) along the folding pathway. By definition, transition states are not significantly populated and are only accessible via folding kinetics studies. In this respect, interpreting the kinetic effects of amino acid substitutions (especially to Ala) via Phi-value analysis is the most common method to probe the structure of these transient, yet important states. A critical review of the key assumptions required for rigorous interpretation of Phi values reveals that a multiple substitution strategy in which a position of interest is mutated to a variety of amino acids, and not exclusively to Ala, provides the best means to characterize folding transition states. This approach has proven useful in revealing non-native interactions and (or) conformations in folding transition states. Moreover, by simultaneously examining the folding kinetics of multiple substitutions made at a single surface-exposed position using the Brønsted analysis the backbone conformation in a folding transition state can be investigated. For folding equilibria with exchange rates on the order of milliseconds, the kinetic parameters for Phi-value analysis can be obtained from NMR relaxation dispersion experiments, under fully native conditions, along with a wealth of high-resolution structural information about the states in exchange (native, denatured, and intermediate states that populate the pathway). This additional structural information, which is not readily obtained through stopped-flow based methods, can significantly facilitate the interpretation of Phi values because it often reports on the validity of the assumptions required for a rigorous interpretation of Phi values.

Helix mutations stabilize a late productive intermediate on the folding pathway of ubiquitin

We have investigated the relative placement of rate-limiting energy barriers and the role of productive or obstructive intermediates on the folding pathway of yeast wild-type ubiquitin ( wt-Ub) containing the F45W mutation. To manipulate the folding barriers, we have designed a family of mutants in which stabilizing substitutions have been introduced incrementally on the solvent-exposed surface of the main alpha-helix (residues 23-34), which has a low intrinsic helical propensity in the native sequence. Although the U --> I and I --> N transitions are not clearly delineated in the kinetics of wt-Ub, we show that an intermediate becomes highly populated and more clearly resolved as the predicted stability of the helix increases. The observed acceleration in the rate of folding correlates with helix stability and is consistent with the I-state representing a productive rather than misfolded state. A Leffler analysis of the effects on kinetics of changes in stability within the family of helix mutants results in a biphasic correlation in both the refolding and unfolding rates that suggest a shift from a nucleation-condensation mechanism (weakly stabilized helix) toward a diffusion-collision model (highly stabilized helix). Through the introduction of helix-stabilizing mutations, we are able to engineer a well-resolved I-state on the folding pathway of ubiquitin which is likely to be structurally distinct from that which is only weakly populated on the folding pathway of wild-type ubiquitin.

Surface loop dynamics in adeno-associated virus capsid assembly

The HI loop is a prominent domain on the adeno-associated virus (AAV) capsid surface that extends from each viral protein (VP) subunit overlapping the neighboring fivefold VP. Despite the highly conserved nature of the residues at the fivefold pore, the HI loops surrounding this critical region vary significantly in amino acid sequence between the AAV serotypes. In order to understand the role of this unique capsid domain, we ablated side chain interactions between the HI loop and the underlying EF loop in the neighboring VP subunit by generating a collection of deletion, insertion, and substitution mutants. A mutant lacking the HI loop was unable to assemble particles, while a substitution mutant (10 glycine residues) assembled particles but was unable to package viral genomes. Substitution mutants carrying corresponding regions from AAV1, AAV4, AAV5, and AAV8 yielded (i) particles with titers and infectivity identical to those of AAV2 (AAV2 HI1 and HI8), (ii) particles with a decreased virus titer (1 log) but normal infectivity (HI4), and (iii) particles that synthesized VPs but were unable to assemble into intact capsids (HI5). AAV5 HI is shorter than all other HI loops by one amino acid. Replacing the missing residue (threonine) in AAV2 HI5 resulted in a moderate particle assembly rescue. In addition, we replaced the HI loop with peptides varying in length and amino acid sequence. This region tolerated seven-amino-acid peptide substitutions unless they spanned a conserved phenylalanine at amino acid position 661. Mutation of this highly conserved phenylalanine to a glycine resulted in a modest decrease in virus titer but a substantial decrease (1 log order) in infectivity. Subsequently, confocal studies revealed that AAV2 F661G is incapable of efficiently completing a key step in the infectious pathway nuclear entry, hinting at a possible perturbation of VP1 phospholipase activity. Molecular modeling studies with the F661G mutant suggest that disruption of interactions between F661 and an underlying P373 residue in the EF loop of the neighboring subunit might adversely affect incorporation of the VP1 subunit at the fivefold axis. Western blot analysis confirmed inefficient incorporation of VP1, as well as a proteolytically processed VP1 subunit that could account for the markedly reduced infectivity. In summary, our studies show that the HI loop, while flexible in amino acid sequence, is critical for AAV capsid assembly, proper VP1 subunit incorporation, and viral genome packaging, all of which implies a potential role for this unique surface domain in viral infectivity.

A disease state mutation unfolds the parkin ubiquitin-like domain

E3 ubiquitin ligases are essential enzymes in the ubiquitination pathway responsible for the recognition of specific E2 conjugating enzymes and for transferring ubiquitin to a substrate targeted for degradation. In autosomal recessive juvenile Parkinson's disease, an early onset form of Parkinson's disease, point mutations in the E3 ligase parkin are one of the most commonly observed traits. Parkin is a multidomain E3 ligase that contains an N-terminal ubiquitin-like domain that interacts with, and effects the ubiquitination of, substrates such as cyclin E, p38 and synphilin. In this work we have examined the folding and structure of the parkin ubiquitin-like domain (Ubld) and of the protein with two causative disease mutations (K48A and R42P). Parallel experiments with the protein ubiquitin were done in order to determine if the same mutations were detrimental to the ubiquitin structure and stability. Despite similar folds between the parkin Ubld and ubiquitin, urea unfolding experiments show that the parkin Ubld is surprisingly approximately 10.6 kJ/mol less stable than ubiquitin. The K48A mutation had little effect on the stability of the parkin Ubld or ubiquitin indicating that this mutation contributes to defective protein-protein interactions. In contrast, the single point mutation R42P in parkin's Ubld caused poor expression and degradation of the protein. To avoid these problems, a GB1-Ubld fusion protein was characterized by NMR spectroscopy to show that the R42P mutation causes the complete unfolding of the parkin Ubld. This observation provides a rationale for the more rapid degradation of parkin carrying the R42P mutation in vivo, and its inability to interact with some substrate proteins. Our work provides the first structural and folding insight into the effects of causative mutations within the ubiquitin-like domain in autosomal recessive juvenile Parkinson's disease.

Structure of the cytoplasmic domain of erythrocyte band 3 hereditary spherocytosis variant P327R: band 3 Tuscaloosa

Previous studies have shown that a single P327R point mutation in the cytoplasmic domain of band 3 (cdb3) protein, known as band 3 Tuscaloosa, leads to a reduction in protein 4.2 content of the erythrocyte membrane and hemolytic anemia. Recent studies have shown that this point mutation does not dissociate the cdb3 dimer, nor does it lead to large-scale rearrangement of the protein structure (Bustos, S. P., and Reithmeier, R. A. F. (2006) Biochemistry 45, 1026-1034). To better define the structural changes in cdb3 that lead to the hemolytic anemia phenotype, site-directed spin labeling (SDSL), in combination with continuous wave electron paramagnetic resonance (EPR) and pulsed double electron-electron resonance (DEER) spectroscopies, has been employed in this study to compare the structure of the R327 variant with wild type P327 cdb3. It is confirmed that the P327R mutation does not dissociate the cdb3 dimer, nor does it change the spatial orientation of the two peripheral domains relative to the dimer interface. However, it does affect the packing of the C-terminal end of helix 10 of the dimerization arms in a subpopulation of cdb3 dimers, it leads to spectral changes at some residues in beta-strand 11 and in the N-terminal end of helix10, and it produces measurable spectral changes at other residues that are near the mutation site. The data indicate that the structural changes are subtle and are localized to one surface of the cdb3 dimer. The spectroscopic description of structural features of the P327R variant provides important clues about the location of one potential protein 4.2 binding surface on cdb3 as well as new insight into the structural basis of the membrane destabilization.

Protein folding kinetics provides a context-independent assessment of beta-strand propensity in the Fyn SH3 domain

Structural database-derived propensities for amino acids to adopt particular local protein structures, such as alpha-helix and beta-strand, have long been recognized and effectively exploited for the prediction of protein secondary structure. However, the experimental verification of database-derived propensities using mutagenesis studies has been problematic, especially for beta-strand propensities, because local structural preferences are often confounded by non-local interactions arising from formation of the native tertiary structure. Thus, the overall thermodynamic stability of a protein is not always altered in a predictable manner by changes in local structural propensity at a single position. In this study, we have undertaken an investigation of the relationship between beta-strand propensity and protein folding kinetics. By characterizing the effects of a wide variety of amino acid substitutions at two different beta-strand positions in an SH3 domain, we have found that the observed changes in protein folding rates are very well correlated to beta-strand propensities for almost all of the substitutions examined. In contrast, there is little correlation between propensities and unfolding rates. These data indicate that beta-strand conformation is well formed in the structured portion of the SH3 domain transition state, and that local structure propensity strongly influences the stability of the transition state. Since the transition state is known to be packed more loosely than the native state and likely lacks many of the non-local stabilizing interactions seen in the native state, we suggest that folding kinetics studies may generally provide an effective means for the experimental validation of database-derived local structural propensities.

The Effect of Increasing the Stability of Non-native Interactions on the Folding Landscape of the Bacterial Immunity Protein Im9

How stabilising non-native interactions influence protein folding energy landscapes is currently not well understood: such interactions could speed folding by reducing the conformational search to the native state, or could slow folding by increasing ruggedness. Here, we examine the influence of non-native interactions in the folding process of the bacterial immunity protein Im9, by exploiting our ability to manipulate the stability of the intermediate and rate-limiting transition state (TS) in the folding of this protein by minor alteration of its sequence or changes in solvent conditions. By analysing the properties of these species using Phi-value analysis, and exploration of the structural properties of the TS ensemble using molecular dynamics simulations, we demonstrate the importance of non-native interactions in immunity protein folding and demonstrate that the rate-limiting step involves partial reorganisation of these interactions as the TS ensemble is traversed. Moreover, we show that increasing the contribution to stability made by non-native interactions results in an increase in Phi-values of the TS ensemble without altering its structural properties or solvent-accessible surface area. The data suggest that the immunity proteins fold on multiple, but closely related, micropathways, resulting in a heterogeneous TS ensemble that responds subtly to mutation or changes in the solvent conditions. Thus, altering the relative strength of native and non-native interactions influences the search to the native state by restricting the pathways through the folding energy landscape.

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.

Population of On-pathway Intermediates in the Folding of Ubiquitin

The role that intermediate states play in protein folding is the subject of intense investigation and in the case of ubiquitin has been controversial. We present fluorescence-detected kinetic data derived from single and double mixing stopped-flow experiments to show that the F45W mutant of ubiquitin (WT*), a well-studied single-domain protein and most recently regarded as a simple two-state system, folds via on-pathway intermediates. To account for the discrepancy we observe between equilibrium and kinetic stabilities and m-values, we show that the polypeptide chain undergoes rapid collapse to an intermediate whose presence we infer from a fast lag phase in interrupted refolding experiments. Double-jump kinetic experiments identify two direct folding phases that are not associated with slow isomerisation reactions in the unfolded state. These two phases are explained by kinetic partitioning which allows molecules to reach the native state from the collapsed state via two possible competing routes, which we further examine using two destabilised ubiquitin mutants. Interrupted refolding experiments allow us to observe the formation and decay of an intermediate along one of these pathways. A plausible model for the folding pathway of ubiquitin is presented that demonstrates that obligatory intermediates and/or chain collapse are important events in restricting the conformational search for the native state of ubiquitin.

Engineering Stabilising β-Sheet Interactions into a Conformationally Flexible Region of the Folding Transition State of Ubiquitin

Protein engineering studies suggest that the transition state for the folding of ubiquitin is highly polarised towards the N-terminal part of the sequence and involves a nucleus of residues within the beta-hairpin (residues 1-17) and main alpha-helix (residues 23-34). In contrast, the observation of small phi-values for residues in the C-terminal portion of the sequence (residues 35-76), coupled with a folding topology that results in a much higher contact order, suggests that fast folding of ubiquitin is dependent upon configurational flexibility in the C-terminal part of the polypeptide chain to ensure passage down a relatively smooth folding funnel to the native state. We show that the introduction of a small mini-hairpin motif as an extension of the native 43-50 hairpin stabilises local interactions in the C-terminal part of the sequence, resulting largely in a deceleration of the unfolding kinetics without perturbing the apparent two-state folding mechanism. However, a single-point Leu-->Phe substitution within the engineered hairpin sequence leads to the premature collapse of the denatured ensemble through the stabilisation of non-native interactions and the population of a compact intermediate. Non-linear effects in the kinetic data at low concentrations of denaturant suggest that the collapsed state, which is further stabilised in the presence of cosmotropic salts, may subsequently fold directly to the native state through a "triangular" reaction scheme involving internal rearrangement rather than unfolding and refolding.

Chemical synthesis of the RGD-protein decorsin: Pro→Ala replacement reduces protein thermostability

Decorsin is a 39-residue polypeptide chain, crosslinked by three disulfide bridges, that strongly inhibits platelet aggregation. We report the chemical synthesis and characterization of analogs of decorsin with the aim of investigating the role of proline residues in protein structure, stability and biological activity. Decorsin analogs have been synthesized in which one (P23A and P24A decorsin) or two (P23,24A decorsin) proline residues have been substituted by alanine. The crude synthetic polypeptides were purified by reversed-phase HPLC in their reduced form and allowed to refold oxidatively to their disulfide-crosslinked species. The homogeneity of the synthetic mini-proteins, and also the correct pairing of the three disulfide bridges, were established by a number of analytical criteria, including fingerprinting analysis of the refolded synthetic analogs by using thermolysin and proteinase K as proteolytic enzymes. Replacement of proline by alanine results in a significant and cumulative decrease of the high thermal stability (Tm 74 degrees C) of native decorsin. The mono-substituted analogs display a Tm of 66-67 degrees C, while the double-substituted analog a Tm of 50 degrees C. On the other hand, the overall secondary and tertiary structures were not affected by the Pro-->Ala exchanges, as judged from circular dichroism measurements. Platelet aggregation assays established that the proline substitutions do not impair significantly the biological activity of decorsin. The results of this study clearly indicate that proline residues contribute significantly to the protein thermal stability. Our results are in line with the 'proline rule', previously advanced for explaining the unusual thermal stability of thermophilic enzymes, which usually show an enhanced content of proline residues with respect to their mesophilic counterparts.

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.