Substrate flow in catalases deduced from the crystal structures of active site variants of HPII from Escherichia coli.

The active site of heme catalases is buried deep inside a structurally highly conserved homotetramer. Channels leading to the active site have been identified as potential routes for substrate flow and product release, although evidence in support of this model is limited. To investigate further the role of protein structure and molecular channels in catalysis, the crystal structures of four active site variants of catalase HPII from Escherichia coli (His128Ala, His128Asn, Asn201Ala, and Asn201His) have been determined at approximately 2.0-A resolution. The solvent organization shows major rearrangements with respect to native HPII, not only in the vicinity of the replaced residues but also in the main molecular channel leading to the heme distal pocket. In the two inactive His128 variants, continuous chains of hydrogen bonded water molecules extend from the molecular surface to the heme distal pocket filling the main channel. The differences in continuity of solvent molecules between the native and variant structures illustrate how sensitive the solvent matrix is to subtle changes in structure. It is hypothesized that the slightly larger H(2)O(2) passing through the channel of the native enzyme will promote the formation of a continuous chain of solvent and peroxide. The structure of the His128Asn variant complexed with hydrogen peroxide has also been determined at 2.3-A resolution, revealing the existence of hydrogen peroxide binding sites both in the heme distal pocket and in the main channel. Unexpectedly, the largest changes in protein structure resulting from peroxide binding are clustered on the heme proximal side and mainly involve residues in only two subunits, leading to a departure from the 222-point group symmetry of the native enzyme. An active role for channels in the selective flow of substrates through the catalase molecule is proposed as an integral feature of the catalytic mechanism. The Asn201His variant of HPII was found to contain unoxidized heme b in combination with the proximal side His-Tyr bond suggesting that the mechanistic pathways of the two reactions can be uncoupled.     

Structural implications of drug-resistant mutants of HIV-1 protease: high-resolution crystal structures of the mutant protease/substrate analogue complexes

Emergence of drug-resistant mutants of HIV-1 protease is an ongoing problem in the fight against AIDS. The mechanisms governing resistance are both complex and varied. We have determined crystal structures of HIV-1 protease mutants, D30N, K45I, N88D, and L90M complexed with peptide inhibitor analogues of CA-p2 and p2-NC cleavage sites in the Gag-pol precursor in order to study the structural mechanisms underlying resistance. The structures were determined at 1.55-1.9-A resolution and compared with the wild-type structure. The conformational disorder seen for most of the hydrophobic side-chains around the inhibitor binding site indicates flexibility of binding. Eight water molecules are conserved in all 9 structures; their location suggests that they are important for catalysis as well as structural stability. Structural differences among the mutants were analyzed in relation to the observed changes in protease activity and stability. Mutant L90M shows steric contacts with the catalytic Asp25 that could destabilize the catalytic loop at the dimer interface, leading to its observed decreased dimer stability and activity. Mutant K45I reduces the mobility of the flap and the inhibitor and contributes to an enhancement in structural stability and activity. The side-chain variations at residue 30 relative to wild-type are the largest in D30N and the changes are consistent with the altered activity observed with peptide substrates. Polar interactions in D30N are maintained, in agreement with the observed urea sensitivity. The side-chains of D30N and N88D are linked through a water molecule suggesting correlated changes at the two sites, as seen with clinical inhibitors. Structural changes seen in N88D are small; however, water molecules that mediate interactions between Asn88 and Thr74/Thr31/Asp30 in other complexes are missing in N88D.     

Model for the three-dimensional structure of vitronectin: predictions for the multi-domain protein from threading and docking.

The structure of vitronectin, an adhesive protein that circulates in high concentrations in human plasma, was predicted through a combination of computational methods and experimental approaches. Fold recognition and sequence-structure alignment were performed using the threading program PROSPECT for each of three structural domains, i.e., the N-terminal somatomedin B domain (residues 1-53), the central region that folds into a four-bladed beta-propeller domain (residues 131-342), and the C-terminal heparin-binding domain (residues 347-459). The atomic structure of each domain was generated using MODELLER, based on the alignment obtained from threading. Docking experiments between the central and C-terminal domains were conducted using the program GRAMM, with limits on the degrees of freedom from a known inter-domain disulfide bridge. The docked structure has a large inter-domain contact surface and defines a putative heparin-binding groove at the inter-domain interface. We also docked heparin together with the combined structure of the central and C-terminal domains, using GRAMM. The predictions from the threading and docking experiments are consistent with experimental data on purified plasma vitronectin pertaining to protease sensitivity, ligand-binding sites, and buried cysteines.     

Model of three-dimensional structure of vitamin D receptor and its binding mechanism with 1alpha,25-dihydroxyvitamin D(3).

Comparative modeling of the vitamin D receptor three-dimensional structure and computational docking of 1alpha,25-dihydroxyvitamin D(3) into the putative binding pocket of the two deletion mutant receptors: (207-423) and (120-422, Delta [164-207]) are reported and evaluated in the context of extensive mutagenic analysis and crystal structure of holo hVDR deletion protein published recently. The obtained molecular model agrees well with the experimentally determined structure. Six different conformers of 1alpha,25-dihydroxyvitamin D(3) were used to study flexible docking to the receptor. On the basis of values of conformational energy of various complexes and their consistency with functional activity, it appears that 1alpha,25-dihydroxyvitamin D(3) binds the receptor in its 6-s-trans form. The two lowest energy complexes obtained from docking the hormone into the deletion protein (207-423) differ in conformation of ring A and orientation of the ligand molecule in the VDR pocket. 1alpha,25-Dihydroxyvitamin D(3) possessing the A-ring conformation with axially oriented 1alpha-hydroxy group binds receptor with its 25-hydroxy substituent oriented toward the center of the receptor cavity, whereas ligand possessing equatorial conformation of 1alpha-hydroxy enters the pocket with A ring directed inward. The latter conformation and orientation of the ligand is consistent with the crystal structure of hVDR deletion mutant (118-425, Delta [165-215]). The lattice model of rVDR (120-422, Delta [164-207]) shows excellent agreement with the crystal structure of the hVDR mutant. The complex obtained from docking the hormone into the receptor has lower energy than complexes for which homology modeling was used. Thus, a simple model of vitamin D receptor with the first two helices deleted can be potentially useful for designing a general structure of ligand, whereas the advanced lattice model is suitable for examining binding sites in the pocket.     

Structure of a legume lectin from the bark of Robinia pseudoacacia and its complex with N-acetylgalactosamine

The structure of the bark lectin RPbAI (isoform A4) from Robinia pseudoacacia has been determined by protein crystallography both in the free form and complexed with N-acetylgalactosamine. The free form is refined at 1.80 A resolution to an R-factor of 18.9% whereas the complexed structure has an R-factor of 19.7% at 2.05 A resolution. Both structures are compared to each other and to other available legume lectin structures. The polypeptide chains of the two structures exhibit the characteristic legume lectin tertiary fold. The quaternary structure resembles that of the Phaseolus vulgaris lectin, the soybean agglutinin, and the Dolichos biflorus lectin, but displays some unique features leading to the extreme stability of this lectin.     

Impact of an easily reducible disulfide bond on the oxidative folding rate of multi-disulfide-containing proteins.

The burial of native disulfide bonds, formed within stable structure in the regeneration of multi-disulfide-containing proteins from their fully reduced states, is a key step in the folding process, as the burial greatly accelerates the oxidative folding rate of the protein by sequestering the native disulfide bonds from thiol-disulfide exchange reactions. Nevertheless, several proteins retain solvent-exposed disulfide bonds in their native structures. Here, we have examined the impact of an easily reducible native disulfide bond on the oxidative folding rate of a protein. Our studies reveal that the susceptibility of the (40-95) disulfide bond of Y92G bovine pancreatic ribonuclease A (RNase A) to reduction results in a reduced rate of oxidative regeneration, compared with wild-type RNase A. In the native state of RNase A, Tyr 92 lies atop its (40-95) disulfide bond, effectively shielding this bond from the reducing agent, thereby promoting protein oxidative regeneration. Our work sheds light on the unique contribution of a local structural element in promoting the oxidative folding of a multi-disulfide-containing protein.     

Mutational and structural-based analyses of the osmolyte effect on protein stability

It is known that several naturally occurring substances known as osmolytes increase the conformational stability of proteins. Bolen and co-worker proposed the osmophobic theory, which asserts the osmolyte effect occurs because of an unfavorable interaction of osmolytes mainly with the protein backbone, based on the results on the transfer Gibbs energy of amino acids (Deltag) [Bolen and Baskakov (2001) J. Mol. Biol. 310, 955-963]. In this paper, we report the effect of sarcosine on the conformational stability (DeltaG) of RNase Sa (96 residues and one disulfide bond) and four mutant proteins. The thermal denaturation curves for RNase Sa in sarcosine fitted a two-state model on nonlinear least-squares analysis. All the RNase Sa proteins were stabilized by sarcosine. For example, the increase in stability of the wild-type protein in 4 M sarcosine due to the osmolyte effect (Delta(o)DeltaG) is 3.2 kcal/mol. Mutational analysis of the osmolyte effect indicated that the changed Delta(o)DeltaG values upon mutation (Delta(m)Delta(o)DeltaG), as estimated from the Deltag values, are similar to the experimental values. Structural-based analysis of the osmolyte effect was also performed using model denatured structures: (a) a fully extended model (single chain) with no disulfide bond, (b) two-part, unfolded models (two chains) with a disulfide bond constructed through molecular dynamic (MD) simulation, and (c) a two-part, folded model (two chains). The two-part, unfolded models were expected to be more suitable as denatured structures. The Delta(o)DeltaG values calculated using the two-part, unfolded models were more consistent with experimental values than those calculated using the fully extended and two-part, folded models. This suggests that MD simulation is useful for testing denatured structures. These results indicate that the osmophobic theory can explain the osmolyte effect on protein stability.     

Ribonuclease Sa conformational stability studied by NMR-monitored hydrogen exchange

The conformational stability of ribonuclease Sa (RNase Sa) has been measured at the per-residue level by NMR-monitored hydrogen exchange at pH* 5.5 and 30 degrees C. In these conditions, the exchange mechanism was found to be EXII. The conformational stability calculated from the slowest exchanging amide groups was found to be 8.8 kcal/mol, in close agreement with values determined by spectroscopic methods. RNase Sa is curiously rich in acidic residues (pI = 3.5) with most basic residues being concentrated in the active-site cleft. The effects of dissolved salts on the stability of RNase Sa was studied by thermal denaturation experiments in NaCl and GdmCl and by comparing hydrogen exchange rates in 0.25 M NaCl to water. The protein was found to be stabilized by salt, with the magnitude of the stabilization being influenced by the solvent exposure and local charge environment at individual amide groups. Amide hydrogen exchange was also measured in 0.25, 0.50, 0.75, and 1.00 M GdmCl to characterize the unfolding events that permit exchange. In contrast to other microbial ribonucleases studied to date, the most protected, globally exchanging amides in RNase Sa lie not chiefly in the central beta strands but in the 3/10 helix and an exterior beta strand. These structural elements are near the Cys7-Cys96 disulfide bond.     

The contribution of polar group burial to protein stability is strongly context-dependent

We previously suggested that proteins gain more stability from the burial and hydrogen bonding of polar groups than from the burial of nonpolar groups (Pace, C. N. (2001) Biochemistry 40, 310-313). To study this further, we prepared eight Thr-to-Val mutants of RNase Sa, four in which the Thr side chain is hydrogen-bonded and four in which it is not. We measured the stability of these mutants by analyzing their thermal denaturation curves. The four hydrogen-bonded Thr side chains contribute 1.3 +/- 0.9 kcal/mol to the stability; those that are not still contribute 0.4 +/- 0.9 kcal/mol to the stability. For 40 Thr-to-Val mutants of 11 proteins, the average decrease in stability is 1.0 +/- 1.0 kcal/mol when the Thr side chain is hydrogen-bonded and 0.0 +/- 0.5 kcal/mol when it is not. This is clear evidence that hydrogen bonds contribute favorably to protein stability. In addition, we prepared four Val-to-Thr mutants of RNase Sa, measured their stability, and determined their crystal structures. In all cases, the mutants are less stable than the wild-type protein, with the decreases in stability ranging from 0.5 to 4.4 kcal/mol. For 41 Val-to-Thr mutants of 11 proteins, the average decrease in stability is 1.8 +/- 1.3 kcal/mol and is unfavorable for 40 of 41 mutants. This shows that placing an [bond]OH group at a site designed for a [bond]CH3 group is very unfavorable. So, [bond]OH groups can contribute favorably to protein stability, even if they are not hydrogen-bonded, if the site was selected for an [bond]OH group, but they will make an unfavorable contribution to stability, even if they are hydrogen-bonded, when they are placed at a site selected for a [bond]CH3 group. The contribution that polar groups make to protein stability depends strongly on their environment.     

Contribution of single tryptophan residues to the fluorescence and stability of ribonuclease Sa

Ribonuclease Sa (RNase Sa) contains no tryptophan (Trp) residues. We have added single Trp residues to RNase Sa at sites where Trp is found in four other microbial ribonucleases, yielding the following variants of RNase Sa: Y52W, Y55W, T76W, and Y81W. We have determined crystal structures of T76W and Y81W at 1.1 and 1.0 A resolution, respectively. We have studied the fluorescence properties and stabilities of the four variants and compared them to wild-type RNase Sa and the other ribonucleases on which they were based. Our results should help others in selecting sites for adding Trp residues to proteins. The most interesting findings are: 1), Y52W is 2.9 kcal/mol less stable than RNase Sa and the fluorescence intensity emission maximum is blue-shifted to 309 nm. Only a Trp in azurin is blue-shifted to a greater extent (308 nm). This blue shift is considerably greater than observed for Trp71 in barnase, the Trp on which Y52W is based. 2), Y55W is 2.1 kcal/mol less stable than RNase Sa and the tryptophan fluorescence is almost completely quenched. In contrast, Trp59 in RNase T1, on which Y55W is based, has a 10-fold greater fluorescence emission intensity. 3), T76W is 0.7 kcal/mol more stable than RNase Sa, indicating that the Trp side chain has more favorable interactions with the protein than the threonine side chain. The fluorescence properties of folded Y76W are similar to those of the unfolded protein, showing that the tryptophan side chain in the folded protein is largely exposed to solvent. This is confirmed by the crystal structure of the T76W which shows that the side chain of the Trp is only approximately 7% buried. 4), Y81W is 0.4 kcal/mol less stable than RNase Sa. Based on the crystal structure of Y81W, the side chain of the Trp is 87% buried. Although all of the Trp side chains in the variants contribute to the unusual positive circular dichroism band observed near 235 nm for RNase Sa, the contribution is greatest for Y81W.     

Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI=3.5) and a basic variant (pI=10.2)

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 degrees C in 0.1M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge-dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds.A primary objective in this study was to investigate the influence of charge-charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the alpha-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge-dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge-charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pK(Sa)-pK(5K)) calculated using a simple model based on Coulomb's Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulomb's Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson-Boltzmann equation, which considers desolvation and charge-dipole interactions in addition to charge-charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge-charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.     

Elucidating protein secondary structures using alpha-carbon recurrence quantifications.

Secondary structures of proteins were studied by recurrence quantification analysis (RQA). High-resolution, 3-dimensional coordinates of alpha-carbon atoms comprising a set of 68 proteins were downloaded from the Protein Data Bank. By fine-tuning four recurrence parameters (radius, line, residue, separation), it was possible to establish excellent agreement between percent contribution of alpha-helix and beta-sheet structures determined independently by RQA and that of the DSSP algorithm (Define Secondary Structure of Proteins). These results indicate that there is an equivalency between these two techniques, which are based upon totally different pattern recognition strategies. RQA enhances qualitative contact maps by quantifying the arrangements of recurrent points of alpha carbons close in 3-dimensional space. For example, the radius was systematically increased, moving the analysis beyond local alpha-carbon neighborhoods in order to capture super-secondary and tertiary structures. However, differences between proteins could only be detected within distances up to about 6-11 A, but not higher. This result underscores the complexity of alpha-carbon spacing when super-secondary structures appear at larger distances. Finally, RQA-defined secondary structures were found to be robust against random displacement of alpha carbons upwards of 1 A. This finding has potential import for the dynamic functions of proteins in motion.     

Adenine recognition: a motif present in ATP-, CoA-, NAD-, NADP-, and FAD-dependent proteins.

Adenosine triphosphate (ATP) plays an essential role in energy transfer within the cell. In the form of NAD, adenine participates in multiple redox reactions. Phosphorylation and ATP-hydrolysis reactions have key roles in signal transduction and regulation of many proteins, especially enzymes. In each cell, proteins with many different functions use adenine and its derivatives as ligands; adenine, of course, is present in DNA and RNA. We show that an adenine binding motif, which differs according to the backbone chain direction of a loop that binds adenine (and in one variant by the participation of an aspartate side-chain), is common to many proteins; it was found from an analysis of all adenylate-containing protein structures from the Protein Data Bank. Indeed, 224 protein-ligand complexes (86 different proteins) from a total of 645 protein structure files bind ATP, CoA, NAD, NADP, FAD, or other adenine-containing ligands, and use the same structural elements to recognize adenine, regardless of whether the ligand is a coenzyme, cofactor, substrate, or an allosteric effector. The common adenine-binding motif shown in this study is simple to construct. It uses only (1) backbone polar interactions that are not dependent on the protein sequence or particular properties of amino acid side-chains, and (2) nonspecific hydrophobic interactions. This is probably why so many different proteins with different functions use this motif to bind an adenylate-containing ligand. The adenylate-binding motif reported is present in "ancient proteins" common to all living organisms, suggesting that adenine-containing ligands and the common motif for binding them were exploited very early in evolution. The geometry of adenine binding by this motif mimics almost exactly the geometry of adenine base-pairing seen in DNA and RNA.     

Crystal structure of protein Ph1481p in complex with protein Ph1877p of archaeal RNase P from Pyrococcus horikoshii OT3: implication of dimer formation of the holoenzyme

Ribonuclease P (RNase P) in the hyperthermophilic archaeon Pyrococcus horikoshii OT3 consists of a catalytic RNA and five protein subunits. We previously determined crystal structures of four protein subunits. Ph1481p, an archaeal homologue for human hPop5, is the protein component of the P.horikoshii RNase P for which no structural information is available. Here we report the crystal structure of Ph1481p in complex with another protein subunit, Ph1877p, determined at 2.0 A resolution. Ph1481p consists of a five-stranded antiparallel beta-sheet and five helices, which fold in a way that is topologically similar to the ribonucleoprotein (RNP) domain. Ph1481p is, however, distinct from the typical RNP domain in that it has additional helices at the C terminus, which pack against one face of the beta-sheet. The presence of two complexes in the asymmetric unit, together with gel filtration chromatography indicates that the heterotetramer is stable in solution and represents a fundamental building block in the crystals. In the heterotetrameric structure (Ph1877p-(Ph1481p)(2)-Ph1877p), a homodimer of Ph1481p sits between two Ph1877p monomers. Ph1481p dimerizes through hydrogen bonding interaction from the loop between alpha1 and alpha2 helices, and each Ph1481p interacts with two Ph1877p molecules, where alpha2 and alpha3 in Ph1481p interact with alpha7 in one Ph1877p and alpha8 in the other Ph1877p molecule, respectively. Deletion of the alpha1-alpha2 loop in Ph1481p caused heterodimerization with Ph1877p, and abolished ability to homodimerize itself and heterotetramerize with Ph1877p. Furthermore, the reconstituted particle containing the deletion mutant Ph1481p (mPh1481p) exhibited significantly reduced nuclease activity. These results suggest the presence of the heterotetramer of Ph1481p and Ph1877p in P.horikoshii RNase P.     

Crystal structure of type 1 ribonuclease H from hyperthermophilic archaeon Sulfolobus tokodaii: role of arginine 118 and C-terminal anchoring

The crystal structure of ribonuclease HI from the hyperthermophilic archaeon Sulfolobus tokodaii (Sto-RNase HI) was determined at 1.6 A resolution. Sto-RNase HI exhibits not only RNase H activity but also double-stranded RNA-dependent ribonuclease (dsRNase) activity. The main-chain fold and steric configurations of the four acidic active-site residues of Sto-RNase HI are very similar to those of other type 1 RNases H. However, Arg118 of Sto-RNase HI is located at the position in which His124 of E. coli RNase HI, His539 of HIV-1 RNase H, and Glu188 of Bacillus halodurans RNase H are located. The mutation of this residue to Ala considerably reduced both the RNase H and dsRNase activities without seriously affecting substrate binding, suggesting that Arg118 is involved in catalytic function. This residue may promote product release by perturbing the coordination of the metal ion A as proposed for Glu188 of B. halodurans RNase H. In addition, the extreme C-terminal region of Sto-RNase HI is anchored to its core region by one disulfide bond and several hydrogen bonds. Differential scanning calorimetry measurements indicated that Sto-RNase HI is a hyperstable protein with a melting temperature of 102 degrees C. The mutations of the cysteine residues forming disulfide bond or elimination of the extreme C-terminal region greatly destabilized the protein, indicating that anchoring of the C-terminal tail is responsible for hyperstabilization of Sto-RNase HI.     

Probing the dark state tertiary structure in the cytoplasmic domain of rhodopsin: proximities between amino acids deduced from spontaneous disulfide bond formation between cysteine pairs engineered in cytoplasmic loops 1, 3, and 4

To probe proximities between amino acids in the cytoplasmic domain by using mutants containing engineered cysteine pairs, three sets of rhodopsin mutants have been prepared. In the first two sets, a cysteine was placed, one at a time, at positions 311-314 in helix VIII, while the second cysteine was fixed at position 246 (set I) and at position 250 (set II) at the cytoplasmic end of helix VI. In the third set, one cysteine was fixed at position 65 while the second cysteine was varied between amino acid positions 306 and 321 located at the cytoplasmic end of helix VII and throughout in helix VIII. Rapid disulfide bond formation in the dark was found between the cysteine pairs in mutants A246C/Q312C,A246C/K311C and in mutants H65C/C316, H65C/315C and H65C/312C. Disulfide bond formation at much lower rates was found in mutants A246C/F313C, V250C/Q312C, H65C/N310C, H65C/K311C, H65C/F313C, and H65C/R314C; the remaining mutants showed no significant disulfide bond formation. Comparisons of the results from disulfide bond formation in solution with the distances observed in the rhodopsin crystal structure showed that the rates of disulfide bond formation in most cases were consistent with the amino acid proximities as revealed in crystal structure. However, deviations were also found, in particular, in the set containing fixed cysteine at position Cys246 and cysteines at positions 311-314. The results implicate significant effects of structural dynamics on disulfide bond formation in solution.     

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.     

Comparison of local and global stability of an analogue of a disulfide-folding intermediate with those of the wild-type protein in bovine pancreatic ribonuclease A: identification of specific regions of stable structure along the oxidative folding pathway

We have identified specific regions of the polypeptide chain of bovine pancreatic ribonuclease A (RNase A) that are critical for stabilizing the oxidative folding intermediate des-[40-95] (with three native disulfide bonds but lacking the fourth native Cys40-Cys95 disulfide bond) in an ensemble of largely disordered three-disulfide precursors (3S if des-[40-95]). A stable analogue of des-[40-95], viz., [C40A, C95A] RNase A, which contains three out of four native disulfide pairings, was previously found to have a three-dimensional structure very similar to that of the wild-type protein. However, it is determined here from GdnHCl denaturation experiments to have significantly reduced global stability, i.e., = 4.5 kcal /mol at 20 degrees C and pH 4.6. The local stability of [C40A, C95A] RNase A was also examined using site-specific amide (2)H/(1)H exchange measurements at pD 5.0 to determine the individual unfolding free energy of specific residues under both strongly native (12 degrees C) and more destabilizing (20 degrees C) conditions. Comparison of the relative stabilities at specific amide sites of [C40A, C95A] RNase A at both temperatures with the corresponding values for the wild-type protein at 35 degrees C corroborates previous experimental evidence that unidentified intramolecular contacts in the vicinity of the preferentially formed native one-disulfide (C65-C72) loop are crucial for stabilizing early folding intermediates, leading to des-[40-95]. Moreover, values of for residues at or near the third alpha-helix, and in part of the second beta-sheet of [C40A, C95A] RNase A, indicate that these two regions of regular backbone structure contribute to stabilizing the global chain fold of the des-[40-95] disulfide-folding intermediate in the wild-type protein. More significantly, we have identified numerous specific residues in the first alpha-helix and the first beta-sheet of the protein that are stabilized in the final step of the major oxidative regeneration pathway of RNase A (des-[40-95] --> N).     

Interactions that favor the native over the non-native disulfide bond among residues 58-72 in the oxidative folding of bovine pancreatic ribonuclease A

In the initial stages of the oxidative folding of both bovine pancreatic ribonuclease A (RNase A) and a 58-72 fragment thereof from the fully reduced, denatured state, the 65-72 correctly paired disulfide bond forms in preponderance over the incorrectly paired 58-65 disulfide bond. Since both disulfide-bonded loops contain the same number of amino acid residues, the question arises as to whether the native pairing results from interactions within the 58-72 segment that lead to a nativelike structure even in its fully reduced form. To answer this question, the chain buildup procedure, based on ECEPP, including a solvation treatment, was used to generate the low-energy structures for the 58-72 RNase segment, beginning with residue 72 and building back to residue 58; in this fragment, all three Cys residues (at positions 58, 65, and 72) initially exist in the reduced (CysH) state. After the open-chain energy minima of the 65-72 peptide were generated, these conformations were allowed to form the 65-72 disulfide bond, and the energies of the resulting oxidized conformations were reminimized and rehydrated. The global minimum of the loop-closed 65-72 structure and many of the low-lying loop-closed minima could be superimposed on the energy-minimized X-ray structure for residues 65-72. The low-energy structures for the full open chain 58-72 peptide were then computed and were allowed to form disulfide bonds either between residues 65 and 72 (native) or between residues 58 and 65 (non-native), and their energies were reminimized and rehydrated in the loop-closed state. Although the overall fold of the 65-72 loop-closed global minimum was the same as for the energy-minimized X-ray structure of these residues, the overall rms deviation was 3.9 A because of local deviations among residues 58-64. In contrast, the 65-72 segment of the global minimum of the 58-72 fragment could be superimposed on the corresponding residues of the energy-minimized X-ray structure. The lowest-energy structure for the 58-65 non-native paired 58-72 sequence was 6 kcal/mol higher in energy than that for the 58-72 peptide with the 65-72 disulfide bond formed. These results suggest that the native pairing of the 65-72 peptide arises from energetic determinants (adoption of left-handed single-residue conformations by Gly 68, and side chain interactions involving Gln 69) contained within this peptide sequence.