Effect of multiple prolyl isomerization reactions on the stability and folding kinetics of the notch ankyrin domain: experiment and theory

Studies on the folding kinetics of the Notch ankyrin domain have demonstrated that the major refolding phase is slow, the minor refolding phase is limited by the isomerization of prolyl peptide bonds, and that unfolding is multiexponential. Here, we explore the relationship between prolyl isomerization and folding heterogeneity using a combination of experiment and simulation. Proline residues were replaced with alanine, both singly and in various combinations. These destabilizing substitutions combine to eliminate the minor refolding phase, although unfolding heterogeneity persists even when all seven proline residues are replaced. To test whether prolyl isomerization influences the major refolding phase, we modeled folding and prolyl isomerization as a system of sequential reactions. Simulations that use rate constants of the major folding phase of the Notch ankyrin domain to represent intrinsic folding indicate that even with seven prolyl isomerization reactions, only two significant phases should be observed, and that the fast observed phase provides a good approximation of the intrinsic folding in the absence of prolyl isomerization. These results indicate that the major refolding phase of the Notch ankyrin domain reflects an intrinsically slow folding transition, rather than coupling of fast folding events with slow prolyl isomerization steps. This is consistent with the observation that the single observed refolding phase of a construct in which all proline residues are replaced remains slow. Finally, the simulation fails to produce a second unfolding phase at high urea concentrations, indicating that prolyl isomerization does not play a role in the three-state mechanism that leads to this heterogeneity.     

Metastable, partially folded states in the productive folding and in the misfolding and amyloid aggregation of proteins

Understanding the energetic and structural basis of protein folding in a physiological context may represent an important step toward the elucidation of protein misfolding and aggregation events that take place in several pathological states. In particular, investigation of the structure and thermodynamic properties of partially folded intermediate states involved in productive folding or in misfolding/aggregation may provide insight into these processes and suggest novel approaches to prevent misfolding in living organisms. This goal, however, has remained elusive, because such intermediates are often transient and correspond to metastable states that are little populated under physiological conditions. Characterization of these states requires their stabilization by means of manipulation of the experimental conditions, involving changes in temperature, pH, or addition of different types of denaturants. In the past few years, hydrostatic pressure has been increasingly used as a thermodynamic variable in the study of both protein folding and misfolding/aggregation transitions. Compared with other chemical or physical denaturing agents, a unique feature of pressure is its ability to induce subtle changes in protein conformation, allowing the stabilization of partially folded states that are usually not significantly populated under more drastic conditions. Much of the recent work in this field has focused on the characterization of folding intermediates, because they seem to be involved in a variety of disease-causing protein misfolding and aggregation reactions. Here, we review recent examples of the use of hydrostatic pressure as a tool to gain insight into the forces and energetics governing the productive folding or the misfolding and amyloid aggregation of proteions.     

Function of phosducin-like proteins in G protein signaling and chaperone-assisted protein folding

Members of the phosducin gene family were initially proposed to act as down-regulators of G protein signaling by binding G protein βγ dimers (Gβγ) and inhibiting their ability to interact with G protein subunits (G) and effectors. However, recent findings have over-turned this hypothesis by showing that most members of the phosducin family act as co-chaperones with the cytosolic chaperonin complex (CCT) to assist in the folding of a variety of proteins from their nascent polypeptides. In fact rather than inhibiting G protein pathways, phosducin-like protein 1 (PhLP1) has been shown to be essential for G protein signaling by catalyzing the folding and assembly of the Gβγ dimer. PhLP2 and PhLP3 have no role in G protein signaling, but they appear to assist in the folding of proteins essential in regulating cell cycle progression as well as actin and tubulin. Phosducin itself is the only family member that does not participate with CCT in protein folding, but it is believed to have a specific role in visual signal transduction to chaperone Gβγ subunits as they translocate to and from the outer and inner segments of photoreceptor cells during light-adaptation.     

Relationships between the folding rate constant and the topological parameters of small two-state proteins based on general random walk model

In this paper, we propose an analytically tractable model of protein folding based on one-dimensional general random walk. A second-order differential equation for the mean folding time of a single protein is constructed which can be used to derive the observed relationship between the folding rate constant and the number of native contacts. The parameters appearing in the model can be determined by fitting the theoretical prediction to the experimental result. In addition, taking into account the fact that the number of native contacts is almost proportional to the relative contact order, we can also explain the observed relationship between the folding rate constant and the relative contact order.     

The coordination of the isomerization of a conserved non-prolyl cis peptide bond with the rate-limiting steps in the folding of dihydrofolate reductase

The propensity for peptide bonds to adopt the trans configuration in native and unfolded proteins, and the relatively slow rates of cis-trans isomerization reactions, imply that the formation of cis peptide bonds in native conformations are likely to limit folding reactions. The role of the conserved cis Gly95-Gly96 peptide bond in dihydrofolate reductase (DHFR) from Escherichia coli was examined by replacing Gly95 with alanine. The introduction of a beta carbon at position 95 is expected to increase the propensity for the trans isomer and perturb the isomerization reaction required to reach the native conformation. Although G95A DHFR is 1.30 kcal mol(-1) less stable than the wild-type protein, it adopts a well-folded structure that can be chemically denatured in a cooperative fashion. The mutant protein also retains the complex refolding kinetic pattern attributed to a parallel-channel mechanism in wild-type DHFR. The spectroscopic response upon refolding monitored by Trp fluorescence and the absence of a Trp/Trp exciton coupling apparent in the far-UV CD spectrum of the wild-type protein, however, indicated that the tertiary structure of the folded state for G95A DHFR is altered. The addition of methotrexate (MTX), a tight-binding inhibitor, to folded G95A DHFR restored the exciton coupling and the fluorescence properties through five slow kinetic events whose relaxation times are independent of the ligand and the denaturant concentrations. The results were interpreted to mean that MTX-binding drives the formation of the cis isomer of the peptide bond between Ala95 and Gly96 in five compact and stable but not wild-type-like conformations that contain the trans isomer. Folding studies in the presence of MTX for both wild-type and G95A DHFR support the notion that the cis peptide bond between Gly95 and Gly96 in the wild-type protein forms during four parallel rate-limiting steps, which are primarily controlled by folding reactions, and lead directly to a set of native, or native-like, conformers. The isomerization of the cis peptide bond is not a source of the parallel channels that characterize the complex folding mechanism for DHFR.     

Structural analysis of the rate-limiting transition states in the folding of Im7 and Im9: similarities and differences in the folding of homologous proteins

The bacterial immunity proteins Im7 and Im9 fold with mechanisms of different kinetic complexity. Whilst Im9 folds in a two-state transition at pH 7.0 and 10 degrees C, Im7 populates an on-pathway intermediate under these conditions. In order to assess the role of sequence versus topology in the folding of these proteins, and to analyse the effect of populating an intermediate on the landscape for folding, we have determined the conformational properties of the rate-limiting transition state for Im9 folding/unfolding using Phi(F)-value analysis and have compared the results with similar data obtained previously for Im7. The data show that the rate-limiting transition states for Im9 and Im7 folding/unfolding are similar: both are compact (beta(T)=0.94 and 0.89, respectively) and contain three of the four native helices docked around a specific hydrophobic core. Significant differences are observed, however, in the magnitude of the Phi(F)-values obtained for the two proteins. Of the 20 residues studied in both proteins, ten have Phi(F)-values in Im7 that exceed those in Im9 by more than 0.2, and of these five differ by more than 0.4. The data suggest that the population of an intermediate in Im7 results in folding via a transition state ensemble that is conformationally restricted relative to that of Im9. The data are consistent with the view that topology is an important determinant of folding. Importantly, however, they also demonstrate that while the folding transition state may be conserved in homologous proteins that fold with two and three-state kinetics, the population of an intermediate can have a significant effect on the breadth of the transition state ensemble.     

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.     

Role of the Cys 2-Cys 10 disulfide bond for the structure, stability, and folding kinetics of ribonuclease T1.

The Cys 2-Cys 10 disulfide bond in ribonuclease T1 was broken by substituting Cys 2 and Cys 10 by Ser and Asn, respectively, as present in ribonuclease F1. This C2S/C10N variant resembles the wild-type protein in structure and in catalytic activity. Minor structural changes were observed by 2-dimensional NMR in the local environment of the substituted amino acids only. The thermodynamic stability of ribonuclease T1 is strongly reduced by breaking the Cys 2-Cys 10 bond, and the free energy of denaturation is decreased by about 10 kJ/mol. The folding mechanism is not affected, and the trans to cis isomerizations of Pro 39 and Pro 55 are still the rate-limiting steps of the folding process. The differences in the time courses of unfolding and refolding are correlated with the decrease in stability: the folding kinetics of the wild-type protein and the C2S/C10N variant become indistinguishable when they are compared under conditions of identical stability. Apparently, the Cys 2-Cys 10 disulfide bond is important for the stability but not for the folding mechanism of ribonuclease T1. The breaking of this bond has the same effect on stability and folding kinetics as adding 1 M guanidinium chloride to the wild-type protein.     

Analysis of the differences in the folding mechanisms of c-type lysozymes based on contact maps constructed with interresidue average distances

A method for analyzing differences in the folding mechanisms of proteins in the same family is presented. Using only information from the amino acid sequences, contact maps derived from the interresidue average distances are employed. These maps, referred to as average distance maps (ADM), are applied to the folding of c-type lysozymes. The results reveal that the ADMs of these lysozymes reflect the differences in the detailed folding mechanisms. Further possible applications of the present method are also discussed.     

Conformational analysis of polypeptides and proteins for the study of protein folding, molecular recognition, and molecular design

Conformational energy calculations provide an understanding as to how interatomic interactions lead to the three-dimensional structures of polypeptides and proteins, and how these molecules interact with other molecules. Illustrative results of such calculations pertain to model systems (-helices and -sheets, and interactions between them), to various open-chain and cyclic peptides, to fibrous proteins, to globular proteins, and to enzyme-substrate complexes. In most cases, the validity of the computations is established by experimental tests of the predicted structures.This article was presented during the proceedings of the International Conference on Macromolecular Structure and Function, held at the National Defence Medical College, Tokorozawa, Japan, December 1985. This paper first appeared in the Israel Journal of Chemistry, Vol. 27, 1986.     

Hidden intermediates and levinthal paradox in the folding of small proteins

It has long been suggested that existence of partially folded intermediates may be essential for proteins to fold in a biologically meaningful time scale. Although partially folded intermediates have been commonly observed in larger proteins, they are generally not detectable in the kinetic folding of smaller proteins (approximately 100 amino acids or less). Recent native-state hydrogen exchange studies suggest that partially folded intermediates may exist behind the rate-limiting transition state in small proteins and evade detection by conventional kinetic methods.     

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

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

The contribution of proline residues to protein stability is associated with isomerization equilibrium in both unfolded and folded states

It has long been understood that the proline residue has lower configurational entropy than any other amino acid residue due to pyrrolidine ring hindrance. The peptide bond between proline and its preceding amino acid (Xaa-Pro) typically exists as a mixture of cis- and trans-isomers in the unfolded protein. Cis–trans isomerization of Xaa-Pro peptide bonds are infrequent, but still occur in folded proteins. Therefore, the effects of the cis–trans isomerization equilibrium in both unfolded and folded states should be taken into account when estimating the stability contribution of a specific proline residue. In order to study the stability contribution of the four proline residues to the hyperthermophilic protein Ssh10b, in this work, we expressed and purified a series of Pro→Ala mutants of Ssh10b, and performed correlative unfolding experiments in detail. We proposed a new unfolding model including proline isomerization. The model predicts that the contribution of a proline residue to protein stability is associated with the thermodynamic equilibrium between cis- and trans-isomers both in the unfolded and folded states, agreeing well with the experimental results.     

Transient intermediary states with high and low folding probabilities in the apparent two-state folding equilibrium of ACBP at low pH

Measurements of the stability as a function of pH for the acyl-coenzyme A binding protein (ACBP) has shown a significant difference in the pH transition midpoint measured by NMR spectroscopy at pH 3.12 and the transition midpoint measured at pH 2.92 and 2.97 by circular dichroism and by fluorescence spectroscopy, respectively. A similar behavior has not been observed in other proteins. It is suggested that these differences arise because the population of the unfolded molecules still contains significant amounts of native like secondary and tertiary structure. NMR spectroscopy measures the concentration of the two components of the folding unfolding equilibrium individually, whereas circular dichroism and fluorescence measure the concentration of the conformations of the light-absorbing chromophores present in both the folded and the unfolded molecules. In the narrow pH range, nascent structure can be detected as the average amount of secondary structure per unfolded molecule and hydrophobic interactions in the population of unfolded molecules. These structures are not observable immediately by NMR spectroscopy; however, a chemical shift analysis of the peptide backbone (13)C chemical shift indicates strongly the existence of short-lived and transient helical structures at pH 2.3. Magnetization transfer studies have been applied to study the equilibrium between folded and unfolded ACBP near the pH transition point measured by NMR. This study has shown that there are two categories of subpopulations in the population of unfolded ACBP. One for which magnetization can be transferred to the folded form during the folding process, and one for which transfer is not observed. The molecules of the latter population of unfolded protein apparently, do not fold within the time-frame of the magnetization transfer experiment. This result suggests the existence of a subpopulation of the acid-unfolded protein molecules with a high propensity for folding. It is suggested that in this subpopulation, a particular set of native like interactions in the peptide backbone and between side-chains in the peptide chain have to be formed.     

Analysis of residuals from enzyme kinetic and protein folding experiments in the presence of correlated experimental noise

Experimental data from continuous enzyme assays or protein folding experiments often contain hundreds, or even thousands, of densely spaced data points. When the sampling interval is extremely short, the experimental data points might not be statistically independent. The resulting neighborhood correlation invalidates important theoretical assumptions of nonlinear regression analysis. As a consequence, certain goodness-of-fit criteria, such as the runs-of-signs test and the autocorrelation function, might indicate a systematic lack of fit even if the experiment does agree very well with the underlying theoretical model. A solution to this problem is to analyze only a subset of the residuals of fit, such that any excessive neighborhood correlation is eliminated. Substrate kinetics of the HIV protease and the unfolding kinetics of UMP/CMP kinase, a globular protein from Dictyostelium discoideum, serve as two illustrative examples. A suitable data-reduction algorithm has been incorporated into software DYNAFIT [P. Kuzmi?, Anal. Biochem. 237 (1996) 260-273], freely available to all academic researchers from http://www.biokin.com.     

An investigation into the N- and C-capping effects of glycine in cavitand-based four-helix bundle proteins

The capping efficiency of glycine on cavitand-based synthetic four-helix bundles was investigated. Glycine, a common C-capping amino acid, has always been included as a C-terminal residue in our de novo peptides, although the exact contribution of the glyince cap to the overall stability and structure of the caviteins had not previously been examined. The uncapped proteins were found to be less helical according to their CD spectra. In addition, the H/D exchange experiments suggested that the uncapped caviteins were more conformationally flexible. Capped and uncapped caviteins exhibited similar values of unfolding. Overall, it can be concluded that glycine caps are useful, as they reduce helical unravelling and enhance helicity, and thus, glycine will be included as a C-terminal residue in future de novo peptide sequences.     

Comparison of transition states obtained upon modeling of unfolding of immunoglobulin-binding domains of proteins L and G caused by external action with transition states obtained in the absence of force probed by experiments

We have studied the extent of coincidence of the pathway of unfolding of protein globules upon experimental modeling of protein unfolding caused by external actions and denaturants. To this end, we compared experimental Φ-values reported in the literature and Φ-values obtained by us upon modeling of unfolding of immunoglobulin-binding domains of proteins L and G caused by external actions at a constant rate. A comparison of the results of calculation with the experimental data shows that the folding pathways for protein L coincide, while those for protein G do not coincide despite structural similarity of these proteins.     

Characterizing sequence variation in the VP1 capsid proteins of foot and mouth disease virus (serotype 0) with respect to virion structure

The VP1 capsid protein of foot and mouth disease virus (FMDV) is highly polymorphic and contains several of the major immunogenic sites important to effective antibody neutralization and subsequent viral clearance by the immune system. Whether this high level of polymorphism is of adaptive value to the virus remains unknown. In this study we examined sequence data from a set of 55 isolates in order to establish the nature of selective pressures acting on this gene. Using the known molecular structure of VP1, the rates and ratios of different types of nonsynonymous and synonymous changes were compared between different parts of the protein. All parts of the protein are subject to purifying selection, but this is greatest amongst those amino acid residues within β-strands and is significantly reduced at residues exposed on the capsid surface, which include those residues demonstrated by previous mutational analyses to permit the virus to escape from monoclonal antibody binding. The ratios of nonsynonymous substitution resulting in various forms of physicochemically radical and conserved amino acid change were shown to be largely equal throughout these different parts of the protein. There was a consistently higher level of nonsynonymous and charge radical sites in those regions of the gene coding for residues exposed on the outer surface of the capsid and a marked difference in the use of amino acids between surface and nonsurface regions of the protein. However, the analysis is consistent with the hypothesis that the observed sequence variation arises where it is least likely to be disruptive to the higher-order structure of the protein and is not necessarily due to positive Darwinian selection. Received: 8 March 1997 / Accepted: 12 August 1997