Folding and unfolding kinetics of the proline-to-alanine mutants of bovine pancreatic ribonuclease A

Nonenzymatic Posttranslational Modifications and Peptide Cleavages Observed in Peptide Epimers

Posttranslational modifications (PTMs) play vital roles in cellular homeostasis and are implicated in various pathological conditions. This work uses two ion mobility spectrometry-mass spectrometry (IMS-MS) modalities, drift-tube IMS (DT-IMS) and trapped IMS (TIMS), to characterize three important nonenzymatic PTMs that induce no mass loss: l/d isomerization, aspartate/isoaspartate isomerization, and cis/trans proline isomerization. These PTMs are assessed in a single peptide system, the recently discovered pleurin peptides, Plrn2, from Aplysia californica. We determine that the DT-IMS-MS/MS can capture and locate asparagine deamidation into aspartate and its subsequent isomerization to isoaspartate, a key biomarker for age-related diseases. Additionally, nonenzymatic peptide cleavage via in-source fragmentation is evaluated for differences in the intensities and patterns of fragment peaks between these PTMs. Peptide fragments resulting from in-source fragmentation, preceded by peptide denaturation by liquid chromatography (LC) mobile phase, exhibited cis/trans proline isomerization. Finally, the effects of differing the fragmentation voltage at the source and solution-based denaturation conditions on in-source fragmentation profiles are evaluated, confirming that LC denaturation and in-source fragmentation profoundly impact N-terminal peptide bond cleavages of Plrn2 and the structures of their fragment ions. With that, LC-IMS-MS/MS coupled with in-source fragmentation could be a robust method to identify three important posttranslational modifications: l/d isomerization, Asn-deamidation leading to Asp/IsoAsp isomerization, and cis/trans proline isomerization.

Differentiating between the sequence of structural events on alternative pathways of folding of a heterodimeric protein

Distinguishing between competing pathways of folding of a protein, on the basis of how they differ in their progress of structure acquisition, remains an important challenge in protein folding studies. A previous study had shown that the heterodimeric protein, double chain monellin (dcMN) switches between alternative folding pathways upon a change in guanidine hydrochloride (GdnHCl) concentration. In the current study, the folding of dcMN has been characterized by the pulsed hydrogen exchange (HX) labeling methodology used in conjunction with mass spectrometry. Quantification of the extent to which folding intermediates accumulate and then disappear with time of folding at both low and high GdnHCl concentrations, where the folding pathways are known to be different, shows that the folding mechanism is describable by a triangular three-state mechanism. Structural characterization of the productive folding intermediates populated on the alternative pathways has enabled the pathways to be differentiated on the basis of the progress of structure acquisition that occurs on them. The intermediates on the two pathways differ in the extent to which the α-helix and the rest of the β-sheet have acquired structure that is protective against HX. The major difference is, however, that β2 has not acquired any protective structure in the intermediate formed on one pathway, but it has acquired significant protective structure in the intermediate formed on the alternative pathway. Hence, the sequence of structural events is different on the two alternative pathways.

Weak Ligaments and Sloping Joints: A New Hypothesis for Development of Congenital Atlantoaxial Dislocation and Basilar Invagination

OBJECTIVE

Developmental bony craniovertebral junction (CVJ) anomalies seem to have a genetic basis and also abnormal joint morphology causing atlantoaxial dislocation (AAD) and basilar invagination (BI).

METHODS

DNA extracted polymerase chain reaction single-stranded conformation polymorphism (SSCP) performed for mutation screening of FBN1 gene (n = 50 cases+ 50 age/sex-matched normal; total: 100). Samples with a deviated pattern of bands in SSCP were sequenced to detect the type of variation. Computed tomography (CT) scans of 100 patients (15-45 years old) compared with an equal number of age/sex-matched controls (21.9 ± 8.2 years). Joint parameters studied: sagittal joint inclination (SI), craniocervical tilt (CCT), coronal joint inclination (CI).

RESULTS

Thirty-nine samples (78%) showed sequence variants. Exon 25, 26, 27, and 28 showed variable patterns of DNA bands in SSCP, which on sequencing gives various types of DNA sequence variations in intronic region of the FBN1 gene in 14%, 14%, 6%, and 44% respectively. CT radiology:SI and CCT correlated with both BI and AAD (p < 0.01). The mean SI value in controls: 83.35° ± 8.65°, and in patients with BI and AAD:129° ± 24.05°. Mean CCT in controls: 60.2° ± 9.2°, and in patients with BI and AAD: 86.0° ± 18.1°. Mean CI in controls:110.3° ± 4.23°, and in cases: 125.15° ± 16.4°.

CONCLUSION

The study showed mutations in FBN1 gene (reported in Marfan syndrome). There is also an alteration of joint morphology, correlating with AAD and BI severity. Hence, we propose a double-hit hypothesis: the presence of weak ligaments (due to FB1 gene alterations) and abnormal joint morphology may contribute to AAD and BI.

Role of cis-trans proline isomerization in the function of pathogenic enterobacterial Periplasmic Binding Proteins

Periplasmic Binding Proteins (PBPs) trap nutrients for their internalization into bacteria by ABC transporters. Ligand binding triggers PBP closure by bringing its two domains together like a Venus flytrap. The atomic determinants that control PBP opening and closure for nutrient capture and release are not known, although it is proposed that opening and ligand release occur while in contact with the ABC transporter for concurrent substrate translocation. In this paper we evaluated the effect of the isomerization of a conserved proline, located near the binding site, on the propensity of PBPs to open and close. ArgT/LAO from Salmonella typhimurium and HisJ from Escherichia coli were studied through molecular mechanics at two different temperatures: 300 and 323 K. Eight microseconds were simulated per protein to analyze protein opening and closure in the absence of the ABC transporter. We show that when the studied proline is in trans, closed empty LAO and HisJ can open. In contrast, with the proline in cis, opening transitions were much less frequent and characterized by smaller changes. The proline in trans also renders the open trap prone to close over a ligand. Our data suggest that the isomerization of this conserved proline modulates the PBP mechanism: the proline in trans allows the exploration of conformational space to produce trap opening and closure, while in cis it restricts PBP movement and could limit ligand release until in productive contact with the ABC transporter. This is the first time that a proline isomerization has been related to the control of a large conformational change like the PBP flytrap mechanism.

Stability and folding of amphibian ribonuclease A superfamily members in comparison with mammalian homologues

Comparative studies on homologous proteins can provide knowledge on how limited changes in the primary structure find their expression in large effects on catalytic activity, stability or the folding behavior. For more than half a century, members of the ribonuclease A superfamily have been the subject of a myriad of studies on protein folding and stability. Both the unfolding and refolding kinetics as well as the structure of several folding intermediates of ribonuclease A have been characterized in detail. Moreover, the RNA-degrading activity of these enzymes provides a basis for their cytotoxicity, which renders them potential tumor therapeutics. Because amphibian ribonuclease A homologues evade the human ribonuclease inhibitor, they emerged as particularly promising candidates. Interestingly, the amphibian ribonuclease A homologues investigated to date are more stable than the mammalian homologues. Nevertheless, despite the generation of numerous genetically engineered variants, knowledge of the folding of amphibian ribonuclease A homologues remains rather limited. An exception is onconase, a ribonuclease A homologue from Rana pipiens, which has been characterized in detail. This review summarizes the data on the unfolding and refolding kinetics and pathways, as well on the stability of amphibian ribonuclease A homologues compared with those of ribonuclease A, the best known member of this superfamily.

An intramolecular chaperone inserted in bacteriophage P22 coat protein mediates its chaperonin-independent folding

The bacteriophage P22 coat protein has the common HK97-like fold but with a genetically inserted domain (I-domain). The role of the I-domain, positioned at the outermost surface of the capsid, is unknown. We hypothesize that the I-domain may act as an intramolecular chaperone because the coat protein folds independently, and many folding mutants are localized to the I-domain. The function of the I-domain was investigated by generating the coat protein core without its I-domain and the isolated I-domain. The core coat protein shows a pronounced folding defect. The isolated I-domain folds autonomously and has a high thermodynamic stability and fast folding kinetics in the presence of a peptidyl prolyl isomerase. Thus, the I-domain provides thermodynamic stability to the full-length coat protein so that it can fold reasonably efficiently while still allowing the HK97-like core to retain the flexibility required for conformational switching during procapsid assembly and maturation.

Consequences of proline-to-alanine substitutions for the stability and refolding of onconase

Peptidyl-prolyl isomerization reactions can make for rate-limiting steps in protein folding due to their high activation energy. Onconase, an unusually stable ribonuclease A homologue from the Northern leopard frog, contains four trans proline residues in its native state. During the refolding from its guanidine hydrochloride unfolded state, which includes the formation of a folding intermediate, the slowest of the three phases has earlier been attributed to a cis-to-trans peptidyl-prolyl isomerization reaction. We thus substituted all four proline residues individually by alanine and investigated the effect of the amino acid substitutions on the folding and stability of the onconase variants. All onconase variants proved to adopt a tertiary structure comparable with that of the wild-type protein. Although the slow phase was not eliminated for any of the variants, the P43A substitution resulted in an increase in the rate constant of the fast folding phase, i.e. a faster formation of the folding intermediate. This variant also exhibits a significant increase in thermodynamic stability. As residue 43 belongs to those residues that are protected from hydrogen exchange with the solvent in the folding intermediate, the increase in the rate constant and stability of the P43A variant emphasizes the importance of the intermediate for the folding of onconase.

Protein prosthesis: β-peptides as reverse-turn surrogates

The introduction of non-natural modules could provide unprecedented control over folding/unfolding behavior, conformational stability, and biological function of proteins. Success requires the interrogation of candidate modules in natural contexts. Here, expressed protein ligation is used to replace a reverse turn in bovine pancreatic ribonuclease (RNase A) with a synthetic β-dipeptide: β²-homoalanine-β³-homoalanine. This segment is known to adopt an unnatural reverse-turn conformation that contains a 10-membered ring hydrogen bond, but one with a donor-acceptor pattern opposite to that in the 10-membered rings of natural reverse turns. The RNase A variant has intact enzymatic activity, but unfolds more quickly and has diminished conformational stability relative to native RNase A. These data indicate that hydrogen-bonding pattern merits careful consideration in the selection of beneficial reverse-turn surrogates.

Chain termini cross-talk in the swapping process of bovine pancreatic ribonuclease

3D domain swapping is the process by which two or more protein molecules exchange part of their structure to form intertwined dimers or higher oligomers. Bovine pancreatic ribonuclease (RNase A) is able to swap the N-terminal α-helix (residues 1-13) and/or the C-terminal β-strand (residues 116-124), thus forming a variety of oligomers, including two different dimers. Cis-trans isomerization of the Asn113-Pro114 peptide group was observed when the protein formed the C-terminal swapped dimer. To study the effect of the substitution of Pro114 on the swapping process of RNase A, we have prepared and characterized the P114A monomeric and dimeric variants of the enzyme. In contrast with previous reports, the crystal structure and NMR data on the monomer reveals a mixed cis-trans conformation for the Asn113-Ala114 peptide group, whereas the X-ray structure of the C-terminal swapped dimer of the variant is very close to that of the corresponding dimer of RNase A. The mutation at the C-terminus affects the capability of the N-terminal α-helix to swap and the stability of both dimeric forms. The present results underscore the importance of the hydration shell in determining the cross-talk between the chain termini in the swapping process of RNase A.

Characterization of secondary amide peptide bond isomerization: thermodynamics and kinetics from 2D NMR spectroscopy

Secondary amide cis peptide bonds are of even lower abundance than the cis tertiary amide bonds of prolines, yet they are of biochemical importance. Using 2D NMR exchange spectroscopy (EXSY) we investigated the formation of cis peptide bonds in several oligopeptides: Ac-G-G-G-NH(2) , Ac-I-G-G-NH(2) , Ac-I-G-G-N-NH(2) and its cyclic form: I-G-G-N in dimethylsulfoxide (DMSO). From the NMR studies, using the amide protons as monitors, an occurrence of 0.13-0.23% of cis bonds was obtained at 296 K. The rate constants for the trans to cis conversion determined from 2D EXSY spectroscopy were 4-9 × 10(-3) s(-1) . Multiple minor conformations were detected for most peptide bonds. From their thermodynamic and kinetic properties the cis isomers are distinguished from minor trans isomers that appear because of an adjacent cis peptide bond. Solvent and sequence effects were investigated utilizing N-methylacetamide (NMA) and various peptides, which revealed a unique enthalpy profile in DMSO. The cyclization of a tetrapeptide resulted in greatly lowered cis populations and slower isomerization rates compared to its linear counterpart, further highlighting the impact of structural constraints.

X-ray crystallographic studies of RNase A variants engineered at the most destabilizing positions of the main hydrophobic core: further insight into protein stability

To investigate the structural origin of decreased pressure and temperature stability, the crystal structure of bovine pancreatic ribonuclease A variants V47A, V54A, V57A, I81A, I106A, and V108A was solved at 1.4-2.0 A resolution and compared with the structure of wild-type protein. The introduced mutations had only minor influence on the global structure of ribonuclease A. The structural changes had individual character that depends on the localization of mutated residue, however, they seemed to expand from mutation site to the rest of the structure. Several different parameters have been evaluated to find correlation with decrease of free energy of unfolding DeltaDeltaG(T), and the most significant correlation was found for main cavity volume change. Analysis of the difference distance matrices revealed that the ribonuclease A molecule is organized into five relatively rigid subdomains with individual response to mutation. This behavior consistent with results of unfolding experiments is an intrinsic feature of ribonuclease A that might be surviving remnants of folding intermediates and reflects the dynamic nature of the molecule.

3S-fluoroproline as a probe to monitor proline isomerization during protein folding by 19F-NMR

Variable-temperature inversion transfer NMR is used to determine the kinetic and thermodynamic parameters of cis-trans isomerization of N-Ac-(3R) and (3S)-fluoroproline-OMe.

Asymmetric kinetics of protein structural changes

Thermodynamic and kinetic understanding of structural transformations in proteins is critical to new developments in medicine and biotechnology. These fields often require the design of mechanism-based modulators of protein function. Researchers increasingly consider these structural changes-such as folding/unfolding or shuttling between active and inactive states-within the energy landscape concept that supposes a high-dimensional, rugged conformational surface. The unevenness, or asperity, of this conformational surface results from energetic barriers and kinetic traps. However, for a large number of protein reactions, such as reversible folding/unfolding, the literature only reports simple two-state transitions, which calls into question the use of a more complex energy landscape model. The question is: are these reactions really that simple, or are we misled by a biased experimental approach? In this Account, we argue in favor of the latter possibility. Indeed, the frequently employed temperature-jump method only allows recording protein structure changes in the heating direction. Under those conditions, it might not be possible to detect other kinetic pathways that could have been taken in the cooling direction. Recently, however, we have developed bidirectional pressure- and temperature-jump methods, which can offer new insights. Here, we show the potential of these methods both for studying protein folding/unfolding reactions, taking ribonuclease A as model, and for studying functionally relevant protein conformational changes, using the open/closed allosteric transition of tryptophan synthase. For example, the heating and cooling temperature-jump induced kinetics involved in the folding/unfolding conformational surface of ribonuclease A is illustrated above. In both of our model systems, the kinetic transition states of several reaction steps were path-dependent, i.e. the rates and thermodynamic activation parameters depend on the direction of the applied pressure and temperature perturbation. This asymmetry suggests that proteins cope with external stress by adapting their structure to form different ensembles of conformational substates. These states are distinguished by their activation enthalpy and entropy barriers, which can be strongly negative in the folding direction. Based on our analysis of activation compressibility and heat capacity, hydration and packing defects of the kinetic transition states are also very important for determining the reaction path. We expect that a more generalized use of this experimental approach should allow researchers to obtain greater insight into the mechanisms of physiologically relevant protein structural changes.

Refolding of ribonuclease A monitored by real-time photo-CIDNP NMR spectroscopy

Photo-CIDNP NMR spectroscopy is a powerful method for investigating the solvent accessibility of histidine, tyrosine and tryptophan residues in a protein. When coupled to real-time NMR, this technique allows changes in the environments of these residues to be used as a probe of protein folding. In this paper we describe experiments performed to monitor the refolding of ribonuclease A following dilution from a high concentration of chemical denaturant. These experiments provide a good example of the utility of this technique which provides information that is difficult to obtain by other biophysical methods. Real-time photo-CIDNP measurements yield residue-specific kinetic data pertaining to the folding reaction, interpreted in terms of current knowledge of the folding of bovine pancreatic ribonuclease A.

Effects of tyrosine mutations on the conformational and oxidative folding of ribonuclease a: a comparative study

Ribonuclease A (RNase A) undergoes more rapid conformational folding with its disulfide bonds intact than during oxidative folding from its reduced form. In this study, the effects of the mutants Y92G, Y92A, and Y92L on both the conformational and oxidative folding pathways were examined to determine the role of native interactions in different types of conformational searches for the biologically active structure of a protein. These mutations did not affect the overall conformational folding pathway of RNase A. However, in the mutants Y92G and Y92A, a key structured disulfide-bonded species, des-[65-72], involved in the oxidative folding pathway of RNase A, was destabilized. These results demonstrate the importance of native interactions in the folding process, namely, protection of a native (40-95) disulfide bond by a nearby tyrosyl-prolyl stacking interaction, when disulfide bonds are allowed to undergo SH/S-S reshuffling.

Distinct unfolding and refolding pathways of ribonuclease a revealed by heating and cooling temperature jumps

Heating and cooling temperature jumps (T-jumps) were performed using a newly developed technique to trigger unfolding and refolding of wild-type ribonuclease A and a tryptophan-containing variant (Y115W). From the linear Arrhenius plots of the microscopic folding and unfolding rate constants, activation enthalpy (DeltaH(#)), and activation entropy (DeltaS(#)) were determined to characterize the kinetic transition states (TS) for the unfolding and refolding reactions. The single TS of the wild-type protein was split into three for the Y115W variant. Two of these transition states, TS1 and TS2, characterize a slow kinetic phase, and one, TS3, a fast phase. Heating T-jumps induced protein unfolding via TS2 and TS3; cooling T-jumps induced refolding via TS1 and TS3. The observed speed of the fast phase increased at lower temperature, due to a strongly negative DeltaH(#) of the folding-rate constant. The results are consistent with a path-dependent protein folding/unfolding mechanism. TS1 and TS2 are likely to reflect X-Pro(114) isomerization in the folded and unfolded protein, respectively, and TS3 the local conformational change of the beta-hairpin comprising Trp(115). A very fast protein folding/unfolding phase appears to precede both processes. The path dependence of the observed kinetics is suggestive of a rugged energy protein folding funnel.

Implementation of a k/k(0) method to identify long-range structure in transition states during conformational folding/unfolding of proteins

A previously introduced kinetic-rate constant (k/k(0)) method, where k and k(0) are the folding (unfolding) rate constants in the mutant and the wild-type forms, respectively, of a protein, has been applied to obtain qualitative information about structure in the transition state ensemble (TSE) of bovine pancreatic ribonuclease A (RNase A), which contains four native disulfide bonds. The method compares the folding (unfolding) kinetics of RNase A, with and without a covalent crosslink and tests whether the crosslinked residues are associated in the folding (unfolding) transition state (TS) of the noncrosslinked version. To confirm that the fifth disulfide bond has not introduced a significant structural perturbation, we solved the crystal structure of the V43C-R85C mutant to 1.6 A resolution. Our findings suggest that residues Val43 and Arg85 are not associated, and that residues Ala4 and Val118 may form nonnative contacts, in the folding (unfolding) TSE of RNase A.

Biophysical Characterisation of the Small Ankyrin Repeat Protein Myotrophin

The 118 residue protein myotrophin is composed of four ankyrin repeats that stack linearly to form an elongated, predominantly alpha-helical structure. The protein folds via a two-state mechanism at equilibrium. The free energy change of unfolding in water (DeltaG(U-N)(H(2)O)) is 5.8 kcal.mol(-1). The chevron plot reveals that the folding reaction has a broad energy barrier and that it conforms to a two-state mechanism. The rate of folding in water (k(f)(H(2)O)) of 95 s(-1) is several orders of magnitude slower than the value predicted by topological calculations. Proline mutants were used to show that the minor kinetic phases observed for myotrophin arise from heterogeneity of the ground states due to cis-trans isomerisation of prolyl as well as non-prolyl peptide bonds. Myotrophin is the first example of a naturally occurring ankyrin repeat protein that conforms to an apparent two-state mechanism at equilibrium and under kinetic conditions, making it highly suitable for high resolution protein folding studies.

Conformational Preferences of Non-Prolyl and Prolyl Residues

The conformational study on Ac-Ala-NHMe (the alanine dipeptide) and Ac-Pro-NHMe (the proline dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the differences in the backbone conformational preference and the cis-trans isomerization for the non-prolyl and prolyl residues in the gas phase and in the solutions (chloroform and water). For the alanine and proline dipeptides, with the increase of solvent polarity, the populations of the conformation tC with an intramolecular C(7) hydrogen bond significantly decrease, and those of the polyproline II-like conformation tF and the alpha-helical conformation tA increase, which is in good agreement with the results from circular dichroism and NMR experiments. For both the dipeptides, as the solvent polarity increases, the relative free energy of the cis conformer to the trans conformer decreases and the rotational barrier to the cis-trans isomerization increases. It is found that the cis-trans isomerization proceeds in common through only the clockwise rotation with omega' approximately +120 degrees about the non-prolyl and prolyl peptide bonds in both the gas phase and the solutions. The pertinent distance d(N...H-N(NHMe)) can successfully describe the increase in the rotational barriers for the non-prolyl and prolyl trans-cis isomerization as the solvent polarity increases and the higher barriers for the non-prolyl residue than for the prolyl residue, as seen in experimental and calculated results. By analysis of the contributions to rotational barriers, the cis-trans isomerization for the non-prolyl and prolyl peptide bonds is proven to be entirely enthalpy driven in the gas phase and in the solutions. The calculated cis populations and rotational barriers to the cis-trans isomerization for both the dipeptides in chloroform and/or water accord with the experimental values.

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.

Genetic selection for critical residues in ribonucleases

Homologous mammalian proteins were subjected to an exhaustive search for residues that are critical to their structure/function. Error-prone polymerase chain reactions were used to generate random mutations in the genes of bovine pancreatic ribonuclease (RNase A) and human angiogenin, and a genetic selection based on the intrinsic cytotoxicity of ribonucleolytic activity was used to isolate inactive variants. Twenty-three of the 124 residues in RNase A were found to be intolerant to substitution with at least one particular amino acid. Twenty-nine of the 123 residues in angiogenin were likewise intolerant. In both RNase A and angiogenin, only six residues appeared to be wholly intolerant to substitution: two histidine residues involved in general acid/base catalysis and four cysteine residues that form two disulfide bonds. With few exceptions, the remaining critical residues were buried in the hydrophobic core of the proteins. Most of these residues were found to tolerate only conservative substitutions. The importance of a particular residue as revealed by this genetic selection correlated with its sequence conservation, though several non-conserved residues were found to be critical for protein structure/function. Despite voluminous research on RNase A, the importance of many residues identified herein was unknown, and those can now serve as targets for future work. Moreover, a comparison of the critical residues in RNase A and human angiogenin, which share only 35% amino acid sequence identity, provides a unique perspective on the molecular evolution of the RNase A superfamily, as well as an impetus for applying this methodology to other ribonucleases.

Pressure-jump-induced kinetics reveals a hydration dependent folding/unfolding mechanism of ribonuclease A

Pressure-jump (p-jump)-induced relaxation kinetics was used to explore the energy landscape of protein folding/unfolding of Y115W, a fluorescent variant of ribonuclease A. Pressure-jumps of 40 MPa amplitude (5 ms dead-time) were conducted both to higher (unfolding) and to lower (folding) pressure, in the range from 100 to 500 MPa, between 30 and 50 degrees C. Significant deviations from the expected symmetrical protein relaxation kinetics were observed. Whereas downward p-jumps resulted always in single exponential kinetics, the kinetics induced by upward p-jumps were biphasic in the low pressure range and monophasic at higher pressures. The relative amplitude of the slow phase decreased as a function of both pressure and temperature. At 50 degrees C, only the fast phase remained. These results can be interpreted within the framework of a two-dimensional energy surface containing a pressure- and temperature-dependent barrier between two unfolded states differing in the isomeric state of the Asn-113-Pro-114 bond. Analysis of the activation volume of the fast kinetic phase revealed a temperature-dependent shift of the unfolding transition state to a larger volume. The observed compensation of this effect by glycerol offers an explanation for its protein stabilizing effect.

The crystal structure of the cis-proline to glycine variant (P114G) of ribonuclease A

Replacement of a cis-proline by glycine at position 114 in ribonuclease A leads to a large decrease in thermal stability and simplifies the refolding kinetics. A crystallographic approach was used to determine whether the decrease in thermal stability results from the presence of a cis glycine peptide bond, or from a localized structural rearrangement caused by the isomerization of the mutated cis 114 peptide bond. The structure was solved at 2.0 A resolution and refined to an R-factor of 19.5% and an R(free) of 21.9%. The overall conformation of the protein was similar to that of wild-type ribonuclease A; however, there was a large localized rearrangement of the mutated loop (residues 110-117-a 9.3 A shift of the Calpha atom of residue 114). The peptide bond before Gly114 is in the trans configuration. Interestingly, a large anomalous difference density was found near residue 114, and was attributed to a bound cesium ion present in the crystallization experiment. The trans isomeric configuration of the peptide bond in the folded state of this mutant is consistent with the refolding kinetics previously reported, and the associated protein conformational change provides an explanation for the decreased thermal stability.

Specifically collapsed intermediate in the early stage of the folding of ribonuclease A

Nature of the burst-phase signals of protein folding has been the subject of much debate as to whether the signals represent the formation of early intermediates or the non-specific collapse of unfolded polypeptides. To distinguish the two possibilities, the submillisecond folding dynamics of ribonuclease A (RNase A) was examined, and compared with those of the disulfide bond-ruptured analog of RNase A (r-RNase A). The circular dichroism measurements on RNase A showed the burst-phase signal within 320 micros after the initiation of the folding reaction, which was identical to that observed for r-RNase A. In contrast, the burst phase increase in the extrinsic fluorescence from 1-anilino-8-naphthalene sulfonate (ANS) was observed for RNase A but not for r-RNase A. The kinetic titration experiment of the ANS fluorescence intensity showed the presence of a specific binding site for ANS in the fast-refolding component of RNase A. The small-angle X-ray scattering measurements at approximately 22 ms after initiating the folding reaction demonstrated that the burst phase conformations of the medium and slow-refolding components of RNase A were distinctly smaller than that of r-RNase A. These results indicated the difference in the burst phase conformations of RNase A and r-RNase A. Since r-RNase A is denatured in the physiological solution condition, the burst-phase signal of RNase A was interpreted as the formation of the folding intermediate with specific conformations.

Probing the unfolding region of ribonuclease A by site-directed mutagenesis

Ribonuclease A contains two exposed loop regions, around Ala20 and Asn34. Only the loop around Ala20 is sufficiently flexible even under native conditions to allow cleavage by nonspecific proteases. In contrast, the loop around Asn34 (together with the adjacent beta-sheet around Thr45) is the first region of the ribonuclease A molecule that becomes susceptible to thermolysin and trypsin under unfolding conditions. This second region therefore has been suggested to be involved in early steps of unfolding and was designated as the unfolding region of the ribonuclease A molecule. Consequently, modifications in this region should have a great impact on the unfolding and, thus, on the thermodynamic stability. Also, if the Ala20 loop contributes to the stability of the ribonuclease A molecule, rigidification of this flexible region should stabilize the entire protein molecule. We substituted several residues in both regions without any dramatic effects on the native conformation and catalytic activity. As a result of their remarkably differing stability, the variants fell into two groups carrying the mutations: (a) A20P, S21P, A20P/S21P, S21L, or N34D; (b) L35S, L35A, F46Y, K31A/R33S, L35S/F46Y, L35A/F46Y, or K31A/R33S/F46Y. The first group showed a thermodynamic and kinetic stability similar to wild-type ribonuclease A, whereas both stabilities of the variants in the second group were greatly decreased, suggesting that the decrease in DeltaG can be mainly attributed to an increased unfolding rate. Although rigidification of the Ala20 loop by introduction of proline did not result in stabilization, disturbance of the network of hydrogen bonds and hydrophobic interactions that interlock the proposed unfolding region dramatically destabilized the ribonuclease A molecule.

Individual tyrosine side-chain contributions to circular dichroism of ribonuclease

Experimental and theoretical studies using site-directed mutants of ribonuclease A (RNase A) offer more extensive information on the tyrosine side-chain contributions to the circular dichroism (CD) of the enzyme. Bovine pancreatic RNase A has three exposed tyrosine residues (Tyr73, Tyr76, and Tyr115) and three buried tyrosine residues (Tyr25, Tyr92 and Tyr97). The difference CD spectra between the wild type and the mutants at pH 7.0 (Deltaepsilon(277,wt) - Deltaepsilon(277,mut)) show bands with more negative DeltaDeltaepsilon(277) values for Y73F and Y115F than those for Y25F and Y92F and bands with positive DeltaDeltaepsilon(277) values for Y76F and Y97F. The theoretical calculations are in good semiquantitative agreement for all the mutants. The pH difference spectrum (pH 11.3-7.0) for the wild type shows a negative band at 295 nm and an enhanced positive band at 245 nm. The three mutants at buried tyrosine sites and one mutant at an exposed tyrosine site (Y76F) exhibit pH-difference spectra that are similar to that of the wild type. In contrast, two mutants at exposed tyrosine sites (Y73F and Y115F) exhibit diminished 295-nm negative bands and, instead of positive bands at 245 nm, negative bands are observed. Our results indicate that Tyr73 and Tyr115, two of the exposed tyrosine residues, are the largest contributors to the 277- and 245-nm CD bands of RNaseA, but the buried tyrosine residues and the one remaining exposed residue also contribute to these bands. Disulfide contributions to the 277- and 240-nm bands and the peptide contribution to the 240-nm band are confirmed theoretically.

Proline replacements and the simplification of the complex, parallel channel folding mechanism for the alpha subunit of Trp synthase, a TIM barrel protein

The kinetic folding mechanism for the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli involves four parallel channels whose inter-conversions are controlled by three cis/trans prolyl isomerization reactions (tau(1), tau(2) and tau(3)). A previous mutational analysis of all 19 proline positions, including the unique cis Asp27-Pro28 peptide bond, revealed that the G(3)P28G, P78A or P96A mutations selectively eliminated the fast, tau(1) (ten seconds), folding phase, while the P217M and P261A mutations eliminated the medium, tau(2) (40 seconds) and the slow, tau(3) ( approximately 300 seconds) folding phases, respectively. To further elucidate the role of these proline residues and to simplify the folding mechanism, a series of double and triple mutants were constructed at these critical positions, and comprehensive kinetic and thermodynamic experiments were performed. Although it was not possible to construct a stable system that was free of proline isomerization constraints, a double mutant variant, G(3)P28G/P217M, in which the refolding of more than 90% of the unfolded protein is not limited by proline isomerization reactions was identified. Further, long-range interactions between several of these residues appear to be a crucial part of the cooperative network of structure that stabilizes the TIM barrel motif for alphaTS.

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.

Protein prosthesis: a nonnatural residue accelerates folding and increases stability

Nonnatural residues can endow proteins with desirable properties. Here, replacing a proline residue that has a cis peptide bond in native ribonuclease A with 5,5-dimethyl-l-proline is shown to accelerate protein folding by 6-fold and enhance conformational stability by DeltaTm = 2.8 +/- 0.3 degrees C while having no effect on enzymatic activity. The rational use of this and other prosthetic segments could enable chemotherapeutic proteins to survive longer in vivo or retain activity after oral administration.

Sequential vs. parallel protein-folding mechanisms: experimental tests for complex folding reactions

The recent emphasis on rough energy landscapes for protein folding reactions by theoreticians, and the many observations of complex folding kinetics by experimentalists provide a rationale for a brief literature survey of various empirical approaches for validating the underlying mechanisms. The determination of the folding mechanism is a key step in defining the energy surface on which the folding reactions occurs and in interpreting the effects of amino acid replacements on this reaction. Case studies that illustrate methods for differentiating between sequential and parallel channel folding mechanisms are presented. The ultimate goal of such efforts is to understand how the one-dimensional information contained in the amino acid sequence is rapidly and efficiently translated into three-dimensional structure.

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.

A cis-prolyl peptide bond isomerization dominates the folding of the alpha subunit of Trp synthase, a TIM barrel protein

The cis/trans isomerization of prolyl peptide bonds has been suggested to dominate the folding of the alpha subunit of tryptophan synthase from Escherichia coli (alphaTS). To test the role of the unique cis isomer between Asp27 and Pro28, the folding properties of P28A, P28G and G(3)P28G, a three-glycine insertion mutant between Asp27 and Gly28, were investigated using urea as a denaturant. Circular dichroism analysis demonstrated that none of the mutations perturb the secondary structure significantly, although the aromatic side-chain packing is altered for P28A and P28G. All three mutant proteins inherited the three-state thermodynamic behavior observed in wild-type alphaTS, ensuring that the fundamental features of the energy surface are intact. Kinetic studies showed that neither alanine nor glycine substitutions at Pro28 results in the elimination of any slow-refolding phases. By contrast, the G(3)P28G mutant eliminates the fastest of the slow-refolding phases and one of the two unfolding phases. Double-jump experiments on G(3)P28G confirm the assignment of the missing refolding phase to the isomerization of the Asp27-Pro28 peptide bond. These results imply that the local stability conveyed by the tight, overlapping turns containing the cis peptide bond is sufficient to favor the cis isomer for several non-prolyl residues. The free energy required to drive the isomerization reaction is provided by the formation of the stable intermediate, demonstrating that the acquisition of structure and stability is required to induce subsequent rate-limiting steps in the folding of alphaTS.

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.

Solution NMR evidence for a cis Tyr-Ala peptide group in the structure of [Pro93Ala] bovine pancreatic ribonuclease A

Proline peptide group isomerization can result in kinetic barriers in protein folding. In particular, the cis proline peptide conformation at Tyr92-Pro93 of bovine pancreatic ribonuclease A (RNase A) has been proposed to be crucial for chain folding initiation. Mutation of this proline-93 to alanine results in an RNase A molecule, P93A, that exhibits unfolding/refolding kinetics consistent with a cis Tyr92-Ala93 peptide group conformation in the folded structure (Dodge RW, Scheraga HA, 1996, Biochemistry 35:1548-1559). Here, we describe the analysis of backbone proton resonance assignments for P93A together with nuclear Overhauser effect data that provide spectroscopic evidence for a type VI beta-bend conformation with a cis Tyr92-Ala93 peptide group in the folded structure. This is in contrast to the reported X-ray crystal structure of [Pro93Gly]-RNase A (Schultz LW, Hargraves SR, Klink TA, Raines RT, 1998, Protein Sci 7:1620-1625), in which Tyr92-Gly93 forms a type-II beta-bend with a trans peptide group conformation. While a glycine residue at position 93 accommodates a type-II bend (with a positive value of phi93), RNase A molecules with either proline or alanine residues at this position appear to require a cis peptide group with a type-VI beta-bend for proper folding. These results support the view that a cis Pro93 conformation is crucial for proper folding of wild-type RNase A.

Pressure versus heat-induced unfolding of ribonuclease A: the case of hydrophobic interactions within a chain-folding initiation site

To investigate the characteristics of the postulated carboxy terminal chain-folding initiation site in bovine pancreatic ribonuclease A (RNase A) (residues 106-118), important in the early stages of the folding pathway, we have engineered by site-directed mutagenesis a set of 14 predominantly conservative hydrophobic variants of the protein. The stability of each variant has been compared by pressure and temperature-induced unfolding, monitored by fourth derivative UV absorbance spectroscopy. Apparently simple two-state, reversible unfolding transitions are observed, suggesting that the disruption of tertiary structure of each protein at high pressure or temperature is strongly cooperative. Within the limits of the technique, we are unable to detect significant differences between the two processes of denaturation. Both steady-state kinetic parameters for the enzyme reaction and UV CD spectra of each RNase A variant indicate that truncation of hydrophobic side chains in this region has, in general, little or no effect on the native structure and function of the enzyme. Furthermore, the decreases in free energy of unfolding upon pressure and thermal denaturation of all the variants, particularly those modified at residues 106 and 108, suggest that the hydrophobic residues and side chain packing interactions of this region play an important role in maintaining the conformational stability of RNase A. We also demonstrate the potential of Tyr115 replacement by Trp as a non-destabilizing fluorescence probe of conformational changes local to the region.