Development of enzymatic activity during protein folding. Detection of a spectroscopically silent native-like intermediate of muscle acylphosphatase

Structural and thermodynamic investigations on the aggregation and folding of acylphosphatase by molecular dynamics simulations and solvation free energy analysis

Protein engineering method to study the mutation effects on muscle acylphosphatase (AcP) has been actively applied to describe kinetics and thermodynamics associated with AcP aggregation as well as folding processes. Despite the extensive mutation experiments, the molecular origin and the structural motifs for aggregation and folding kinetics as well as thermodynamics of AcP have not been rationalized at the atomic resolution. To this end, we have investigated the mutation effects on the structures and thermodynamics for the aggregation and folding of AcP by using the combination of fully atomistic, explicit-water molecular dynamics simulations, and three-dimensional reference interaction site model theory. The results indicate that the A30G mutant with the fastest experimental aggregation rate displays considerably decreased α1-helical contents as well as disrupted hydrophobic core compared to the wild-type AcP. Increased solvation free energy as well as hydrophobicity upon A30G mutation is achieved due to the dehydration of hydrophilic side chains in the disrupted α1-helix region of A30G. In contrast, the Y91Q mutant with the slowest aggregation rate shows a non-native H-bonding network spanning the mutation site to hydrophobic core and α1-helix region, which rigidifies the native state protein conformation with the enhanced α1-helicity. Furthermore, Y91Q exhibits decreased solvation free energy and hydrophobicity compared to wild type due to more exposed and solvated hydrophilic side chains in the α1-region. On the other hand, the experimentally observed slower folding rates in both mutants are accompanied by decreased helicity in α2-helix upon mutation. We here provide the atomic-level structures and thermodynamic quantities of AcP mutants and rationalize the structural origin for the changes that occur upon introduction of those mutations along the AcP aggregation and folding processes.

Biological function in a non-native partially folded state of a protein

As structural flexibility is known to be required for enzyme catalysis and pattern recognition and a significant fraction of eukaryotic proteins appear to be unfolded or contain unstructured regions, biological activity of conformational states distinct from fully folded structures could be more common than previously thought. By applying a procedure that allows the recovery of enzymatic activity to be monitored in real time, we show that a non-native state populated transiently during folding of the acylphosphatase from Sulfolobus solfataricus is enzymatically active. The structural characterization of this partially folded state reveals that enzymatic activity is possible even if the catalytic site is structurally heterogeneous, whereas the remainder of the structure acts as a scaffold. These results extend the spectrum of biological functions carried out in the absence of a folded state to include enzyme catalysis.

The Oxidative Folding Rate of Bovine Pancreatic Ribonuclease Is Enhanced by a Covalently Attached Oligosaccharide

Bovine pancreatic ribonuclease B (RNase B) differs from RNase A by the presence of an oligosaccharide moiety covalently attached to Asn 34. Oxidative folding studies of RNase B were carried out at different temperatures using DTT(ox) as the oxidizing agent, and the results were compared with those for RNase A. The oxidative folding rates of RNase B are between 1.7 and 1.3 times faster than those of RNase A at the temperatures that were investigated. The folding pathways of RNase B were determined to be similar to those of RNase A in that two structured intermediates, each lacking one native disulfide bond, were found to populate the regeneration pathways at 25 degrees C and pH 8.3. The thermodynamic stabilities of these two glycosylated intermediates, and their rates of formation from their unstructured precursors in the rate-determining step, were found to be higher than those of their unglycosylated counterparts from RNase A. Thus, the underlying cause for the faster rate of oxidative regeneration of native RNase B appears to be both thermodynamic and kinetic due to the higher stability, and faster rate of formation, of the intermediates of RNase B compared to those of RNase A.

Reversal of protein aggregation provides evidence for multiple aggregated States

Observations that prefibrillar aggregates from different amyloidogenic proteins can be solubilised under some conditions have raised questions as to the generality of this phenomenon and the nature of the factors that influence it. By studying aggregates formed from human muscle acylphosphatase (AcP) under mild denaturing conditions, and by using a battery of techniques, we demonstrate that disaggregation is possible under conditions close to physiological where the protein is stable in its native state. In the presence of 25% (v/v) trifluoroethanol (TFE) AcP undergoes partial unfolding and globular aggregates (60-200 nm in diameter) that can assemble further into clusters (400-800 nm in diameter) develop progressively. Yet larger superstructures (>5 microm) are formed when the concentration of the globular aggregates exceeds a critical concentration. After diluting the sample to give a solution containing 5% TFE, the fraction of partially unfolded monomeric protein refolds very rapidly, with a rate constant of approximately 1s(-1). The 60-200 nm globular aggregates disaggregate with an apparent rate constant of approximately 2.5 x 10(-3)s(-1) while the 400-800 nm clusters disassembly more slowly with a rate constant of approximately 3.1 x 10(-4)s(-1). The larger (>5 microm) superstructures are not disrupted under the conditions used here. These results suggest that amyloid formation occurs in discrete steps whose reversibility is increasingly difficult, and dependent on the size of the aggregates, and that disaggregation experiments can provide a powerful method of detecting different types of species within the complex process of aggregation. In addition, our work suggests that destabilization of amyloid aggregates resulting in the conversion of misfolded proteins back to their native states could be an important factor in both the onset and treatment of diseases associated with protein aggregation.

The folding of spectrin domains I: wild-type domains have the same stability but very different kinetic properties

The study of proteins with the same architecture, but different sequence has proven to be a valuable tool in the protein folding field. As a prelude to studies on the folding mechanism of spectrin domains we present the kinetic characterisation of the wild-type forms of the 15th, 16th, and 17th domains of chicken brain alpha-spectrin (referred to as R15, R16 and R17, respectively). We show that the proteins all behave in a two-state manner, with different kinetic properties. The folding rate varies remarkably between different members, with a 5000-fold variation in folding rate and 3000-fold variation in unfolding rate seen for proteins differing only 1 kcal mol(-1) in stability. We show clear evidence for significant complexity in the energy landscape of R16, which shows a change in amplitude outside the stopped-flow timescale and curvature in the unfolding arm of the chevron plot. The accompanying paper describes the characterisation of the folding pathway of this domain.

Catalysis of proline-directed protein phosphorylation by peptidyl-prolyl cis/trans isomerases

Proline-directed protein phosphorylation was shown to depend on the capacity of the targeted Ser(Thr)-Pro bond to exhibit conformational polymorphism. The cis/trans isomer specificity underlying ERK2-catalyzed phosphate transfer leads to a complete discrimination of the cis Ser(Thr)-Pro conformer of oligopeptide substrates. We investigated in vitro the ERK2-catalyzed phosphorylation of Aspergillus oryzae RNase T1 containing two Ser-Pro bonds both of which share high stabilization energy in their respective native state conformation, the cis Ser54-Pro and the trans Ser72-Pro moiety. Despite trans isomer specificity of ERK2, a doubly phosphorylated RNase T1 was found as the final reaction product. Similarly, the RNase T1 S54G/P55N and RNase T1 P73V variants, which retain the prolyl bond conformations of the RNase T1-wt, were both monophosphorylated with a catalytic efficiency kcat/KM of 425 M(-1) s(-1) and 1228 M(-1) s(-1), respectively. However, initial phosphorylation rates did not depend linearly on the ERK2 concentration. The phosphorylation rate of the resulting plateau region at high ERK2 concentrations can be increased up to threefold for the RNase T1 P73V variant in the presence of the peptidyl-prolyl cis/trans isomerase Cyclophilin 18, indicating a conformational interconversion as the rate limiting step in the catalyzed phosphate group transfer. Using peptidyl-prolyl cis/trans isomerases with different substrate specificity, we identified a native state conformational equilibrium of the Ser54-Pro bond with the minor trans Ser54-Pro bond as the phosphorylation-sensitive moiety. This technique can therefore be used for a determination of the ratio and the interconversion rates of prolyl bond isomers in the native state of proteins.

Integrating folding kinetics and protein function: biphasic kinetics and dual binding specificity in a WW domain

Because of the association of beta-sheet formation with the initiation and propagation of amyloid diseases, model systems have been sought to further our understanding of this process. WW domains have been proposed as one such model system. Whereas the folding of the WW domains from human Yes-associated protein (YAP) and Pin have been shown to obey single-exponential kinetics, the folding of the WW domain from formin-binding protein (FBP) 28 has been shown to proceed via biphasic kinetics. From an analysis of free-energy landscapes from atomic-level molecular dynamics simulations, the biphasic folding kinetics observed in the FBP WW domain may be traced to the ability of this WW domain to adopt two slightly different forms of packing in its hydrophobic core. This conformational change is propagated along the peptide backbone and affects the position of a tryptophan residue shown in other WW domains to play a key role in binding. The WW domains of Pin and YAP do not support more than one type of packing each, leading to monophasic folding kinetics. The ability of the FBP WW domain to assume two different types of packing may, in turn, explain the capacity of this WW domain to bind two classes of ligand, a property that is not shared by other WW domains. These findings lead to the hypothesis that lability with respect to conformations separated by an observable barrier as a requirement for function is incompatible with the ability of a protein to fold via single-exponential kinetics.

Regulation of peptide bond cis/trans isomerization by enzyme catalysis and its implication in physiological processes

In some cases, the slow rotational movement underlying peptide bond cis/trans isomerizations is found to control the biological activity of proteins. Peptide bond cis/trans isomerases as cyclophilins, Fk506-binding proteins, parvulins, and bacterial hsp70 generally assist in the interconversion of the polypeptide substrate cis/trans isomers, and rate acceleration is the dominating mechanism of action in cells. We present evidence disputing the hypothesis that some of the molecular properties of these proteins play an auxiliary role in enzyme function.

The kinetics of side chain stabilization during protein folding

The rate of stabilization of side chains during protein folding has never been carefully studied. Recent developments in labeling proteins with (19)F-labeled amino acids coupled with real-time NMR measurements have allowed such measurements to be made. This paper describes the application of this method to the study of several proteins using 6-(19)F-tryptophan as the reporting group. It is found that these side chains adopt their final stable state at the last stages of the folding process and that the stabilization of side chains into their final conformation is a highly cooperative process. It is also possible to show the presence of intermediates in which the side chains are not correctly packed. The technique should be applicable to many systems.

Comparison of the folding processes of distantly related proteins. Importance of hydrophobic content in folding

The N-terminal domain of HypF from Escherichia coli (HypF-N) is a 91 residue protein module sharing the same folding topology and a significant sequence identity with two extensively studied human proteins, muscle and common-type acylphosphatases (mAcP and ctAcP). With the aim of learning fundamental aspects of protein folding from the close comparison of so similar proteins, the folding process of HypF-N has been studied using stopped-flow fluorescence. While mAcP and ctAcP fold in a two-state fashion, HypF-N was found to collapse into a partially folded intermediate before reaching the fully folded conformation. Formation of a burst-phase intermediate is indicated by the roll over in the Chevron plot at low urea concentrations and by the large jump of intrinsic and 8-anilino-1-naphtalenesulphonic acid-derived fluorescence immediately after removal of denaturant. Furthermore, HypF-N was found to fold rapidly with a rate constant that is approximately two and three orders of magnitudes faster than ctAcP and mAcP, respectively. Differences between the bacterial protein and the two human counterparts were also found as to the involvement of proline isomerism in their respective folding processes. The results clearly indicate that features that are often thought to be relevant in protein folding are not highly conserved in the evolution of the acylphosphatase superfamily. The large difference in folding rate between mAcP and HypF-N cannot be entirely accounted for by the difference in relative contact order or related topological metrics. The analysis shows that the higher folding rate of HypF-N is in part due to the relatively high hydrophobic content of this protein. This conclusion, which is also supported by the highly significant correlation found between folding rate and hydrophobic content within a group of proteins displaying the topology of HypF-N and AcPs, suggests that the average hydrophobicity of a protein sequence is an important determinant of its folding rate.

Characterization of single-tryptophan mutants of histidine-containing phosphocarrier protein: evidence for local rearrangements during folding from high concentrations of denaturant

We have used site-directed mutagenesis in combination with a battery of biophysical techniques to probe the stability and folding behavior of a small globular protein, the histidine-containing phosphocarrier protein (HPr). Specifically, the four phenylalanine residues (2, 22, 29, and 48) of the wild-type protein were individually replaced by single tryptophans, thus introducing site-specific probes for monitoring the behavior of the protein. The folding of the tryptophan mutants was investigated by NMR, DSC, CD, intrinsic fluorescence, fluorescence anisotropy, and fluorescence quenching. The heat-induced denaturation of all four mutants, and the GdnHCl-induced unfolding curves of F2W, F29W, and F48W, can be fitted adequately to a two-state model, in agreement with the observations for the wild-type protein. The GdnHCl unfolding transitions of F22W, however, showed the accumulation of an intermediate state at low concentrations of denaturant. Kinetic refolding studies of F2W, F29W, and F48W showed a major single phase, independent of the probe used (CD, fluorescence, and fluorescence anisotropy) and similar to that of the wild-type protein. In contrast, F22W showed two phases in the fluorescence experiments corresponding to the two phases previously observed in ANS binding studies of the wild-type protein [Van Nuland et al. (1998) Biochemistry 37, 622-637]. Residue 22 was found from NMR studies to be part of the binding interface on HPr for ANS. These observations indicate that the second slow phase reflects a local, rather than a global, rearrangement from a well-structured highly nativelike intermediate state to the fully folded native state that has less hydrophobic surface exposed to the solvent. The detection of the second slow phase by the use of selective labeling of different regions of the protein with fluorophores illustrates the need for an integrated approach in order to understand the intricate details of the folding reactions of even the simplest proteins.

Kinetic mechanism and catalysis of a native-state prolyl isomerization reaction

Folding of tendamistat is a rapid two-state process for the majority of the unfolded molecules. In fluorescence-monitored refolding kinetics about 8% of the unfolded molecules fold slowly (lambda=0.083s(-1)), limited by peptidyl-prolyl cis-trans isomerization. This is significantly less than expected from the presence of three trans prolyl-peptide bonds in the native state. In interrupted refolding experiments we detected an additional very slow folding reaction (lambda=0.008s(-1) at pH 2) with an amplitude of about 12%. This reaction is caused by the interconversion of a highly structured intermediate to native tendamistat. The intermediate has essentially native spectroscopic properties and about 2% of it remain populated in equilibrium after folding is complete. Catalysis by human cyclophilin 18 identifies this very slow reaction as a prolyl isomerization reaction. This shows that prolyl-isomerases are able to efficiently catalyze native state isomerization reactions, which allows them to influence biologically important regulatory conformational transitions. Folding kinetics of the proline variants P7A, P9A, P50A and P7A/P9A show that the very slow reaction is due to isomerization of the Glu6-Pro7 and Ala8-Pro9 peptide bonds, which are located in a region that makes strong backbone and side-chain interactions to both beta-sheets. In the P50A variant the very slow isomerization reaction is still present but native state heterogeneity is not observed any more, indicating a long-range destabilizing effect on the alternative native state relative to N. These results enable us to include all prolyl and non-prolyl peptide bond isomerization reactions in the folding mechanism of tendamistat and to characterize the kinetic mechanism and the energetics of a native-state prolyl isomerization reaction.

The involvement of mammalian and plant FK506-binding proteins (FKBPs) in development

The FK506-binding proteins (FKBPs) are peptidyl prolyl cis/trans isomerases and the information gathered in the last 10 years reveals their involvement in diverse biological systems affecting the function and structure of target proteins. Members of the FKBP family were shown to be growth-regulated and participate in signal transduction. In this review we have chosen to focus on a few examples of the mammalian and plant systems in which members of the FKBP family have been demonstrated to affect the function of proteins or development. The technologies that enable production of knockout mice, Arabidopsis mutants and overexpression in transgenic organisms have revealed the contribution of FKBP to development in higher eukaryotes. It appears that members of the FKBP family have conserved some of their basic functions in the animal and plant kingdom, whereas other functions became unique. Studies that will take advantage of the full genome sequence available for Arabidopsis and the human genome, DNA chip technologies and the use of transgenic complementation system will contribute to the elucidation of the molecular mechanism and biological function of FKBPs.

Folding and aggregation are selectively influenced by the conformational preferences of the alpha-helices of muscle acylphosphatase

The native state of human muscle acylphosphatase (AcP) presents two alpha-helices. In this study we have investigated folding and aggregation of a number of protein variants having mutations aimed at changing the propensity of these helical regions. Equilibrium and kinetic measurements of folding indicate that only helix-2, spanning residues 55-67, is largely stabilized in the transition state for folding therefore playing a relevant role in this process. On the contrary, the aggregation rate appears to vary only for the variants in which the propensity of the region corresponding to helix-1, spanning residues 22-32, is changed. Mutations that stabilize the first helix slow down the aggregation process while those that destabilize it increase the aggregation rate. AcP variants with the first helix destabilized aggregate with rates increased to different extents depending on whether the introduced mutations also alter the propensity to form beta-sheet structure. The fact that the first alpha-helix is important for aggregation and the second helix is important for folding indicates that these processes are highly specific. This partitioning does not reflect the difference in intrinsic alpha-helical propensities of the two helices, because helix-1 is the one presenting the highest propensity. Both processes of folding and aggregation do not therefore initiate from regions that have simply secondary structure propensities favorable for such processes. The identification of the regions involved in aggregation and the understanding of the factors that promote such a process are of fundamental importance to elucidate the principles by which proteins have evolved and for successful protein design.

Direct measurement indicates a slow cis/trans isomerization at the secondary amide peptide bond of glycylglycine

Spectral differences between the cis and the trans isomer of a secondary amide peptide bond were used to follow the time course of the cis/trans isomerization of Gly-Gly, Gly-Ala, Ala-Gly, and Ala-Ala dipeptides in the UV/vis region at 220 nm. Isomerization rates and Eyring activation energies were calculated from pH- and LiCl-mediated solvent jump experiments. Rate constants were found to be in a narrow range of 0.29 to 0.64 s(-)(1) for the zwitterionic dipeptides at 25 degrees C. The isomerization rate is about 2-fold higher for the monoionic forms of Gly-Gly. The zwitterionic Gly-Gly has an activation enthalpy DeltaH() of 71.6 +/- 4.9 kJ mol(-)(1) that is in the range of the rotational barriers of aromatic side chain dipeptides that have been measured by (1)H NMR magnetization transfer experiments. Late stages of protein backbone rearrangements often involve crossing the energy barrier for rotational isomerization of imidic peptide bonds. Our findings are consistent with the idea that a wide range of secondary amide peptide bonds are also able to induce slow rate-limiting steps in protein restructuring.

Detection of two partially structured species in the folding process of the amyloidogenic protein beta 2-microglobulin

beta 2-Microglobulin is a small, major histocompatibility complex class I-associated protein that undergoes aggregation and accumulates as amyloid deposits in human tissues as a consequence of long-term haemodialysis. The folding process of this amyloidogenic protein has been studied in vitro by diluting the guanidine hydrochloride-denatured protein in refolding buffer at pH 7.4 and monitoring the folding process by means of a number of spectroscopic probes that allow the native structure of the protein to be detected as it develops. These techniques include fluorescence spectroscopy, far and near-UV circular dichroism, 8-anilino-1-naphthalenesulfonic acid binding and double jump assays. All spectroscopic probes indicate that a significant amount of structure forms within the dead-time of stopped-flow measurements (<5 ms). The folding reaction goes to completion through a fast phase followed by a slow phase, whose rate constants are ca 5.1 and 0.0030 s(-1) in water, respectively. Unfolding-folding double jump experiments, together with the use of peptidyl prolyl isomerase, reveal that the slow phase of folding of beta 2-microglobulin is not fundamentally determined by cis/trans isomerisation of X-Pro peptide bonds. Other folding-unfolding double jump experiments also suggest that the fast and slow phases of folding are not related to independent folding of different populations of protein molecules. Rather, we provide evidence for a sequential mechanism of folding where denatured beta 2-microglobulin collapses to an ensemble of partially folded conformations (I(1)) which fold subsequently to a more highly structured species (I(2)) and, finally, attain the native state. The partially folded species I(2) appears to be closely similar to previously studied amyloidogenic forms of beta 2-microglobulin, such as those adopted by the protein at mildly acid pH values and by a variant with six residues deleted at the N terminus. Since amyloid formation in vivo originates from partial denaturation of beta 2-microglobulin under conditions favouring the folding process, the long-lived, partially structured species detected here might be significantly populated under some physiological conditions and hence might play an important role in the process of amyloid formation.

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.

Initial denaturing conditions influence the slow folding phase of acylphosphatase associated with proline isomerization

The folding kinetics of human common-type acylphosphatase (cAcP) from its urea- and TFE-denatured states have been determined by stopped-flow fluorescence techniques. The refolding reaction from the highly unfolded state formed in urea is characterized by double exponential behavior that includes a slow phase associated with isomerism of the Gly53-Pro54 peptide bond. However, this slow phase is absent when refolding is initiated by dilution of the highly a-helical denatured state formed in the presence of 40% trifluoroethanol (TFE). NMR studies of a peptide fragment corresponding to residues Gly53-Gly69 of cAcP indicate that only the native-like trans isomer of the Gly-Pro peptide bond is significantly populated in the presence of TFE, whereas both the cis and trans isomers are found in an approximately 1:9 ratio for the peptide bond in aqueous solution. Molecular modeling studies in conjunction with NMR experiments suggest that the trans isomer of the Gly53-Pro54 peptide bond is stabilized in TFE by the formation of a nonnative-like hydrogen bond between the CO group of Gly53 and the NH group of Lys57. These results therefore reveal that a specific nonnative interaction in the denatured state can increase significantly the overall efficiency of refolding.

Stabilisation of alpha-helices by site-directed mutagenesis reveals the importance of secondary structure in the transition state for acylphosphatase folding

The effects of stabilising mutations on the folding process of common-type acylphosphatase have been investigated. The mutations were designed to increase the helical propensity of the regions of the polypeptide chain corresponding to the two alpha-helices of the native protein. Various synthetic peptides incorporating the designed mutations were produced and their helical content estimated by circular dichroism. The most substantial increase in helical content is found for the peptide carrying five mutations in the second alpha-helix. Acylphosphatase variants containing the corresponding mutations display, to different extents, enhanced conformational stabilities as indicated by equilibrium urea denaturation experiments monitored by changes of intrinsic fluorescence. All the protein variants studied here refold with apparent two-state kinetics. Mutations in the first alpha-helix are responsible for a small increase in the refolding rate, accompanied by a marked decrease in the unfolding rate. On the other hand, multiple mutations in the second helix result in a considerable increase in the refolding rate without any significant effect on the unfolding rate. Addition of trifluoroethanol was found to accelerate the folding of the acylphosphatase variants, the extent of the acceleration being inversely proportional to the intrinsic rate of folding of the corresponding mutant. The trifluoroethanol-induced acceleration is far less marked for those variants whose alpha-helical structure is efficiently stabilised by amino acid replacements. This observation suggests that trifluoroethanol acts in a similar manner to the stabilising mutations in promoting native-like secondary structure. Analysis of the kinetic data indicates that the second helix is fully consolidated in the transition state for folding of acylphosphatase, whereas the first helix is only partially formed. These data suggest that the second helix is an important element in the folding process of the protein.

Enzymes that catalyse the restructuring of proteins

The large enzyme families of protein disulfide isomerases and peptidyl prolyl cis/trans isomerases have been shown to assist polypeptide restructuring. Various folding states of polypeptides may serve as substrates of the catalysed reaction. Our understanding of the cellular function of these enzymes is increasing as a result of the availability of more specific inhibitors, the discovery of natural substrates and the use of genetically modified organisms. Further highlights of these studies include insights into the three-dimensional structures of enzyme-ligand complexes, as well as into the mechanism of slow folding phases on the atomic level.