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.     

Increasing temperature accelerates protein unfolding without changing the pathway of unfolding

We have traditionally relied on extremely elevated temperatures (498K, 225 degrees C) to investigate the unfolding process of proteins within the timescale available to molecular dynamics simulations with explicit solvent. However, recent advances in computer hardware have allowed us to extend our thermal denaturation studies to much lower temperatures. Here we describe the results of simulations of chymotrypsin inhibitor 2 at seven temperatures, ranging from 298K to 498K. The simulation lengths vary from 94ns to 20ns, for a total simulation time of 344ns, or 0.34 micros. At 298K, the protein is very stable over the full 50ns simulation. At 348K, corresponding to the experimentally observed melting temperature of CI2, the protein unfolds over the first 25ns, explores partially unfolded conformations for 20ns, and then refolds over the last 35ns. Above its melting temperature, complete thermal denaturation occurs in an activated process. Early unfolding is characterized by sliding or breathing motions in the protein core, leading to an unfolding transition state with a weakened core and some loss of secondary structure. After the unfolding transition, the core contacts are rapidly lost as the protein passes on to the fully denatured ensemble. While the overall character and order of events in the unfolding process are well conserved across temperatures, there are substantial differences in the timescales over which these events take place. We conclude that 498K simulations are suitable for elucidating the details of protein unfolding at a minimum of computational expense.     

Influence of an extrinsic cross-link on the folding pathway of ribonuclease A. Kinetics of folding-unfolding

The kinetics of folding/unfolding of cross-linked Lys7-dinitrophenylene-Lys41-ribonuclease A were studied and compared to those of unmodified ribonuclease A (RNase A) at various concentrations of guanidine hydrochloride. The folding of the denatured cross-linked protein involved one fast-folding species (22 +/- 4%) and two slow-folding species, as observed in unmodified ribonuclease A. Also, a nativelike intermediate, analogous to that reported previously for unmodified ribonuclease A [Cook, K. H., Schmid, F. X., & Baldwin, R. L. (1979) Proc. Natl. Acad. Sci. U.S.A. 76, 6157], has been detected on the folding pathway of cross-linked ribonuclease A. The extrinsic cross-link between Lys7 and Lys41 did not affect the rate constants for the folding kinetics of these three species. The cross-link did, however, significantly affect the rate constant for unfolding of the native protein. The conformation of the protein in the transition state of the unfolding pathway was deduced from an analysis of the kinetic data. It appears that the 41 N-terminal residues are unfolded in the transition state of the unfolding pathway. Thus, the unfolding pathway of RNase A is sequential in that further unfolding (after the transition state) follows the unfolding of the 41 N-terminal residues. Also, the conformation of the 41 N-terminal residues does not play a role in the folding pathway. Presumably, if the cross-link were introduced instead between two other residues that are in the segment(s) involved in the rate-limiting step(s), it could increase the refolding rate constants and possibly the concentration of fast-folding species.     

Local structure in a tryptic fragment of performic acid oxidized ribonuclease A corresponding to a proposed polypeptide chain-folding initiation site detected by tyrosine fluorescence lifetime and proton magnetic resonance measurements

The effects of proline and X-Pro peptide bond conformations on the fluorescence properties of tyrosine in peptides corresponding to parts of a proposed chain-folding initiation site in bovine pancreatic ribonuclease A are examined by time-resolved and steady-state fluorescence spectroscopy. In peptides with Tyr-Pro sequences, the conformational constraints of proline on a preceding residue result in significant fluorescence quenching for both trans and cis peptide bond conformations. Small peptides containing Pro-Tyr sequences, on the other hand, do not exhibit fluorescence quenching compared to Ac-Tyr-NHMe. Studies of fluorescence decay in the tryptic fragment of performic acid oxidized ribonuclease corresponding to residues 105-124 (i.e., O-T-16) demonstrate the presence of at least two environments of the single tyrosine chromophore (in the sequence Asn113-Pro114-Tyr115). In these two (ensemble-averaged) environments, tyrosine has shorter and longer lifetimes, respectively, than in Ac-Tyr-NHMe. The fluorescence heterogeneity in O-T-16 does not correlate with X-Pro cis/trans conformational heterogeneity that can be detected by nuclear magnetic resonance (NMR) spectroscopy. Instead, the fluorescence heterogeneity in O-T-16 arises from the presence of multiple conformations with the same X-Pro peptide bond conformations which interconvert rapidly on the 1H NMR time scale (tau much less than 1 ms) but are distinguishable on the fluorescence lifetime time scale (tau greater than or equal to 1 ns). From comparisons with the tyrosine fluorescence decay of smaller synthetic peptides, it is concluded that the long-lifetime tyrosine fluorescence component of O-T-16 arises from interactions involving residues outside the Asn113-Pro114-Tyr115-Val116-Pro117 sequence, which either stabilize particular local conformations in the vicinity of Tyr115 or act directly to protect Tyr115 from efficient fluorescence quenching. The short-lifetime component of O-T-16 is also observed for the pentapeptide Ac-Asn-Pro-Tyr-Val-Pro-NHMe. The data provide evidence for a nonrandom polypeptide conformation of O-T-16 under conditions of solvent pH and temperature at which the complete disulfide-intact ribonuclease molecule is fully folded. Implications of this work for the interpretation of fluorescence-detected unfolding experiments are discussed.     

The asparagine-stabilized beta-turn of apamin: contribution to structural stability from dynamics simulation and amide hydrogen exchange analysis

Molecular dynamics simulations of bee venom apamin, and an analogue having an Asn to Ala substitution at residue 2 (apamin-N2A), were analyzed to explore the contribution of hydrogen bonds involving Asn2 to local (beta-turn residues N2, C3, K4, A5) and global stability. The wild-type peptide retained a stable conformation during 2.4 ns of simulation at 67 degrees C, with high beta-turn stability characterized by backbone-side chain hydrogen bonds involving beta-turn residues K4 and A5, with the N2 side chain amide carbonyl. The loss of stabilizing interactions involving the N2 side chain resulted in the loss of the beta-turn conformation in the apamin N2A simulations (27 or 67 degrees C). This loss of beta-turn stability propagates throughout the peptide structure, with destabilization of the C-terminal helix connected to the N-terminal region by two disulfide bonds. Backbone stability in a synthetic peptide analogue (apamin-N2A) was characterized by NMR and amide hydrogen exchange measurements. Consistent with the simulations, loss of hydrogen bonds involving the N2 side chain resulted in destabilization of both the N-terminal beta-turn and the C-terminal helix. Amide exchange protection factors in the C-terminal helix were reduced by 9-11-fold in apamin N2A as compared with apamin, corresponding to free energy (deltaDeltaG(uf)) of around 1.5 kcal M(-1) at 20 degrees C. This is equivalent to the contribution of hydrogen bond interactions involving the N2 side chain to the stability of the beta-turn. Together with additional measures of exchange protection factors, the three main contributions to backbone stability in apamin that account for virtually the full thermodynamic stability of the peptide have been quantitated.     

Molecular dynamics simulations of synthetic peptide folding

Because the time scale of protein folding is much greater than that of the widely used simulations of native structures, a detailed report of molecular dynamics simulations of folding has not been available. In this study, we Included the average solvent effect in the potential functions to simplify the calculation of the solvent effect and carried out long molecular dynamics simulations of the alanine-based synthetic peptides at 274 K. From either an extended or a randomly generated conformation, the simulations approached a helix-coil equilibrium in about 3 ns. The multiple minima problem did not prevent helix folding. The calculated helical ratio of Ac-AAQAAAAQAAAAQAAY-NH₂ was 47%, in good agreement with the circular dichroism measurement (about 50%). A helical segment with frayed ends was the most stable conformation, but the hydrophobic interaction favored the compact, distorted helix-turn-helix conformations. The transition between the two types of conformations occurred in a much larger time scale than helix propagation. The transient hydrogen bonds between the glutamine side chain and the backbone carbonyl group could reduce the free energy barrier of helix folding and unfolding. The substitution of a single alanine residue in the middle of the peptide with valine or glycine decreased the average helical ratio significantly, in agreement with experimental observations. © 1996 Wiley-Liss, Inc.     

N-terminal truncation of prion protein affects both formation and conformation of abnormal protease-resistant prion protein generated in vitro

Transmissible spongiform encephalopathy diseases are characterized by conversion of the normal protease-sensitive host prion protein, PrP-sen, to an abnormal protease-resistant form, PrP-res. In the current study, deletions were introduced into the flexible tail of PrP-sen (23) to determine if this region was required for formation of PrP-res in a cell-free assay. PrP-res formation was significantly reduced by deletion of residues 34-94 relative to full-length hamster PrP. Deletion of another nineteen amino acids to residue 113 further reduced the amount of PrP-res formed. Furthermore, the presence of additional proteinase K cleavage sites indicated that deletion to residue 113 generated a protease-resistant product with an altered conformation. Conversion of PrP deletion mutants was also affected by post-translational modifications to PrP-sen. Conversion of unglycosylated PrP-sen appeared to alter both the amount and the conformation of protease-resistant PrP-res produced from N-terminally truncated PrP-sen. The N-terminal region also affected the ability of hamster PrP to block mouse PrP-res formation in scrapie-infected mouse neuroblastoma cells. Thus, regions within the flexible N-terminal tail of PrP influenced interactions required for both generating and disrupting PrP-res formation.     

Structural studies of a folding intermediate of bovine pancreatic ribonuclease A by continuous recycled flow

A new technique, continuous recycled flow (CRF) spectroscopy, has been developed for observing intermediates of any thermally induced, reversible reaction with a half-life of 10 s or longer. The structure can be probed by any spectroscopic method which does not perturb the system. Prolonged signal acquisitions of 8 h for ribonuclease A are possible. CRF was used to investigate the structure of the slow-folding intermediates of chemically intact ribonuclease A (RNase A) during thermal unfolding/folding under acidic conditions. The following conclusions were reached on the basis of the proton nuclear magnetic resonance and far-ultraviolet circular dichroism spectra of a folding intermediate(s): (A) The conformation of the detected folding intermediate(s) is similar to that of the heat-denatured protein. There is only limited formation of new structures. (B) The N-terminal alpha-helix is partially stable under these conditions and is in rapid (less than 10 ms) equilibrium with the denatured conformation. (C) There are long-range interactions between the hydrophobic residues of the N-terminal alpha-helix and the rest of the protein. These interactions persist well above the melting point. (D) An aliphatic methyl group reports on the formation of a new structure(s) that lie(s) outside of the N-terminal region. (E) The structures detected in chemically modified, nonfolding forms of the RNase A are also present in the folding intermediate(s). There are, however, additional interactions that are unique to chemically intact RNase A.     

Dissecting the stability of a beta-hairpin peptide that folds in water: NMR and molecular dynamics analysis of the beta-turn and beta-strand contributions to folding.

NMR studies of the folding and conformational properties of a beta-hairpin peptide, several peptide fragments of the hairpin, and sequence-modified analogues, have enabled the various contributions to beta-hairpin stability in water to be dissected. Temperature and pH-induced unfolding studies indicate that the folding-unfolding equilibrium approximates to a two-state model. The hairpin is highly resistant to denaturation and is still significantly folded in 7 M urea at 298 K. Thermodynamic analysis shows the hairpin to fold in water with a significant change in heat capacity, however, DeltaCp degrees in 7 M urea is reduced. V/Y-->A mutations on one strand of the hairpin reduce folding to <10 %, consistent with a hydrophobic stabilisation model. We show that in a truncated peptide (residues 6-16) lacking the hydrophobic residues on one beta-strand, the type I' Asn-Gly turn in the sequence SINGKK is significantly populated in water in the absence of interstrand hydrophobic contacts. Unrestrained molecular dynamics simulations of unfolding, using an explicit solvation model, show that the conformation of the NG turn persists for longer than the AG analogue, which has a much lower propensity for type I' turn formation from a data base analysis of preferred turns. The origin of the high stability of the Asn-Gly turn is not entirely clear; data base analysis of 66 NG turns, together with molecular dynamics simulations, reveals no participation of the Asn side-chain in turn-stabilising interactions with the peptide backbone. However, hydration analysis of the molecular dynamics simulations reveals a pocket of "high density" water bridging between the Asn side-chain and peptide main-chain that suggests solvent-mediated interactions may play an important role in modulating phi,psi propensities in the NG turn region.     

Structure of the Murine Unglycosylated IgG1 Fc Fragment

A prototypic IgG antibody can be divided into two major structural units: the antigen-binding fragment (Fab) and the Fc fragment that mediates effector functions. The IgG Fc fragment is a homodimer of the two C-terminal domains (C_H2 and C_H3) of the heavy chains. Characteristic of the Fc part is the presence of a sugar moiety at the inner face of the C_H2 domains. The structure of this complex branched oligosaccharide is generally resolved in crystal structures of Fc fragments due to numerous well-defined sugar-protein interactions and a small number of sugar-sugar interactions. This suggested that sugars play an important role in the structure of the Fc fragment. To address this question directly, we determined the crystal structure of the unglycosylated Fc fragment of the murine IgG1 MAK33. The structures of the C_H3 domains of the unglycosylated Fc fragment superimpose perfectly with the structure of the isolated MAK33 C_H3 domain. The unglycosylated C_H2 domains, in contrast, approach each other much more closely compared to known structures of partly deglycosylated Fc fragments with rigid-body motions between 10 and 14 Å, leading to a strongly “closed” conformation of the unglycosylated Fc fragment. The glycosylation sites in the C′E loop and the BC and FG loops are well defined in the unglycosylated C_H2 domain, however, with increased mobility and with a significant displacement of about 4.9 Å for the unglycosylated Asn residue compared to the glycosylated structure. Thus, glycosylation both stabilizes the C′E-loop conformation within the C_H2 domain and also helps to ensure an “open” conformation, as seen upon Fc receptor binding. These structural data provide a rationale for the observation that deglycosylation of antibodies often compromises their ability to bind and activate Fcγ receptors.     

Local secondary structure in ribonuclease A denatured by guanidine · HCl near 1 °C

Recent work has shown that a short α-helix can be stable in water near 1 °C when stabilized by specific interactions between side-chains, while earlier “host-guest” results with random copolymers have shown that a short α-helix is unstable in water at all temperatures in the absence of stabilizing side-chain interactions. As regards the mechanism of protein folding, it is now reasonable on energetic grounds to consider isolated α-helices and β-sheets as the first intermediates on the pathway of protein folding. Proton nuclear magnetic resonance is used here to detect isolated secondary structures in ribonuclease A denatured by guanidine · HCl (GuHCl). Temperatures near 1 °C are used because the low-temperature stability of the C-peptide helix may be a general property of isolated secondary structures in water.Our procedure is to titrate with GuHCl the C2H resonance lines of the four histidine residues of denatured ribonuclease A. Studies of model peptides (C-peptide (lactone) and C-peptide carboxylate, residues 1 to 13 of ribonuclease A; S-peptide, residues 1 to 20) show linear titration curves for the C2H resonance of His12 above 0.5 M-GuHCl, once helix unfolding is complete. Deviations from this line are used to monitor helix formation. The GuHCl titration curves of the other three histidine residues are also linear, once unfolding is complete. The results show that the helix found in C-peptide and S-peptide is also found in denatured ribonuclease A, where it behaves as an isolated helix not stabilized significantly by interactions with other chain segments. Studies of denatured S-protein show that the remaining three His residues, His48, His105 and His119, are involved in structure only below 1 m-GuHCl at 9 °C, pH 1.9. The nature of this structure is not known. The main conclusion from this work is that the His12 helix can be observed as a stable, isolated helix in denatured ribonuclease A near 1 °C, and that none of the other three His residues is involved in a comparably stable local structure. In native ribonuclease A, His12 is within an α-helix and the other three His residues are involved in a 3-stranded β-sheet structure.The helix-coil transition of C-peptide has also been studied for other side-chain resonances by GuHCl titration. Typically, but not always, the titration curves are linear after helix unfolding takes place and resonance lines from different residues of the same amino acid type can be resolved in GuHCl solutions. This is true of the four histidine residues of ribonuclease A although their pK values in 5 m-GuHCl are nearly the same. In C-peptide, the βCH₃ resonance of Ala6 is affected strongly by GuHCl while the lines of Ala4 and Ala5 are shifted only weakly by GuHCl. Evidently the interactions between GuHCl and side-chains in an unfolded peptide depend upon neighboring groups.     

Structural and dynamic effects of alpha-helix deletion in Sso7d: implications for protein thermal stability

Sso7d is a 62-residue protein from the hyperthemophilic archaeon Sulfolobus solfataricus with a denaturation temperature close to 100 degrees C around neutral pH. An engineered form of Sso7d truncated at leucine 54 (L54Delta) is significantly less stable, with a denaturation temperature of 53 degrees C. Molecular dynamics (MD) studies of Sso7d and its truncated form at two different temperatures have been performed. The results of the MD simulations at 300 K indicate that: (1) the flexibility of Sso7d chain at 300 K agrees with that detected from X-ray and NMR structural studies; (2) L54Delta remains stable in the native folded conformation and possesses an overall dynamic behavior similar to that of the parent protein. MD simulations performed at 500 K, 10 ns long, indicate that, while Sso7d is in-silico resistant to high temperature, the truncated variant partially unfolds, revealing the early phases of the thermal unfolding pathway of the protein. Analysis of the trajectories of L54Delta suggests that the unzipping of the N-terminal and C-terminal beta-strands should be the first event of the unfolding pathway, and points out the regions more resistant to thermal unfolding. These findings allow one to understand the role played by specific interactions connecting the two ends of the chain for the high thermal stability of Sso7d, and support recent hypotheses on its folding mechanism emerged from site-directed mutagenesis studies.     

Thermostabilization of Escherichia coli ribonuclease HI by replacing left-handed helical Lys95 with Gly or Asn.

From the systematic replacements of amino acid residues of Escherichia coli ribonuclease HI with those of its thermophilic counterpart, the basic protrusion domain including region 6 (R6) from residues 91 to 95 was found to increase the structural stability of the mutant protein (Kimura, S., Nakamura, H., Hashimoto, T., Oobatake, M., and Kanaya, S. (1992) J. Biol. Chem. 267, 21535-21542). Further mutagenesis concentrating in the R6 region has revealed that replacements of Lys95 at the left-handed structure with Gly or Asn essentially enhances the protein stability. Gly and Asn substitutions stabilize the protein up to 1.9 kcal/mol and 0.9 kcal/mol in the free energy changes of unfolding, respectively. We propose that the amino acid substitution of left-handed non-Gly residue with Gly or Asn residue can be used as one of the general strategies to enhance protein stability, when such a non-Gly residue itself does not seriously contribute to protein stability.     

The role of a beta-bulge in the folding of the beta-hairpin structure in ubiquitin

It is known that the peptide corresponding to the N-terminal beta-hairpin of ubiquitin, U(1-17), can populate the monomeric beta-hairpin conformation in aqueous solution. In this study, we show that the Gly-10 that forms the bulge of the beta-turn in this hairpin is very important to the stability of the hairpin. The deletion of this residue to desG10(1-16) unfolds the structure of the peptide in water. Even under denaturing conditions, this bulge appears to be important in maintaining the residual structure of ubiquitin, which involves tertiary interactions within the sequence 1 to 34 in the denatured state. We surmise that this residual structure functions as one of the nucleation centers in the folding process and is important in stabilizing the transition state. In accordance with this idea, deleting Gly-10 slows down the refolding and unfolding rate by about one half.     

Pressure as a denaturing agent in studies of single-point mutants of an amyloidogenic protein human cystatin c

Recently, we presented a convenient method combining a deuterium‐hydrogen exchange and electrospray mass spectrometry for studying high‐pressure denaturation of proteins (Stefanowicz et al., Biosci Rep 2009; 30:91–99). Here, we present results of pressure‐induced denaturation studies of an amyloidogenic protein—the wild‐type human cystatin C (hCC) and its single‐point mutants, in which Val57 residue from the hinge region was substituted by Asn, Asp or Pro, respectively. The place of mutation and the substituting residues were chosen mainly on a basis of theoretical calculations. Observation of H/D isotopic exchange proceeding during pressure induced unfolding and subsequent refolding allowed us to detect differences in the proteins stability and folding dynamics. On the basis of the obtained results we can conclude that proline residue at the hinge region makes cystatin C structure more flexible and dynamic, what probably facilitates the dimerization process of this hCC variant. Polar asparagine does not influence stability of hCC conformation significantly, whereas charged aspartic acid in 57 position makes the protein structure slightly more prone to unfolding. Our experiments also point out pressure denaturation as a valuable supplementary method in denaturation studies of mutated proteins. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Equilibrium structure and folding of a helix-forming peptide: circular dichroism measurements and replica-exchange molecular dynamics simulations

We have performed experimental measurements and computer simulations of the equilibrium structure and folding of a 21-residue alpha-helical heteropeptide. Far ultraviolet circular dichroism spectroscopy is used to identify the presence of helical structure and to measure the thermal unfolding curve. The observed melting temperature is 296 K, with a folding enthalpy of -11.6 kcal/mol and entropy of -39.6 cal/(mol K). Our simulations involve 45 ns of replica-exchange molecular dynamics of the peptide, using eight replicas at temperatures between 280 and 450 K, and the program CHARMM with a continuum solvent model. In a 30-ns simulation started from a helical structure, conformational equilibrium at all temperatures was reached after 15 ns. This simulation was used to calculate the peptide melting curve, predicting a folding transition with a melting temperature in the 330-350 K range, enthalpy change of -10 kcal/mol, and entropy change of -30 cal/(mol K). The simulation results were also used to analyze the peptide structural fluctuations and the free-energy surface of helix unfolding. In a separate 15-ns replica-exchange molecular dynamics simulation started from the extended structure, the helical conformation was first attained after approximately 2.8 ns, and equilibrium was reached after 10 ns of simulation. These results showed a sequential folding process with a systematic increase in the number of hydrogen bonds until the helical state is reached, and confirmed that the alpha-helical state is the global free-energy minimum for the peptide at low temperatures.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

Charge-charge interactions in the denatured state influence the folding kinetics of ribonuclease Sa

Gaining a better understanding of the denatured state ensemble of proteins is important for understanding protein stability and the mechanism of protein folding. We studied the folding kinetics of ribonuclease Sa (RNase Sa) and a charge-reversal variant (D17R). The refolding kinetics are similar, but the unfolding rate constant is 10-fold greater for the variant. This suggests that charge-charge interactions in the denatured state and the transition state ensembles are more favorable in the variant than in RNase Sa, and shows that charge-charge interactions can influence the kinetics and mechanism of protein folding.     

Effect of urea on peptide conformation in water: molecular dynamics and experimental characterization

Molecular dynamics simulations of a ribonuclease A C-peptide analog and a sequence variant were performed in water at 277 and 300 K and in 8 M urea to clarify the molecular denaturation mechanism induced by urea and the early events in protein unfolding. Spectroscopic characterization of the peptides showed that the C-peptide analog had a high alpha-helical content, which was not the case for the variant. In the simulations, interdependent side-chain interactions were responsible for the high stability of the alpha-helical C-peptide analog in the different solvents. The other peptide displayed alpha-helical unwinding that propagated cooperatively toward the N-terminal. The conformations sampled by the peptides depended on their sequence and on the solvent. The ability of water molecules to form hydrogen bonds to the peptide as well as the hydrogen bond lifetimes increased in the presence of urea, whereas water mobility was reduced near the peptide. Urea accumulated in excess around the peptide, to which it formed long-lived hydrogen bonds. The unfolding mechanisms induced by thermal denaturation and by urea are of a different nature, with urea-aqueous solutions providing a better peptide solvation than pure water. Our results suggest that the effect of urea on the chemical denaturation process involves both the direct and indirect mechanisms.