Pressure- and temperature-induced unfolding studies: thermodynamics of core hydrophobicity and packing of ribonuclease A

Comprehensive mutagenesis to identify amino acid residues contributing to the difference in thermostability between two originally thermostable ancestral proteins

Further improvement of the thermostability of inherently thermostable proteins is an attractive challenge because more thermostable proteins are industrially more useful and serve as better scaffolds for protein engineering. To establish guidelines that can be applied for the rational design of hyperthermostable proteins, we compared the amino acid sequences of two ancestral nucleoside diphosphate kinases, Arc1 and Bac1, reconstructed in our previous study. Although Bac1 is a thermostable protein whose unfolding temperature is around 100°C, Arc1 is much more thermostable with an unfolding temperature of 114°C. However, only 12 out of 139 amino acids are different between the two sequences. In this study, one or a combination of amino acid(s) in Bac1 was/were substituted by a residue(s) found in Arc1 at the same position(s). The best mutant, which contained three amino acid substitutions (S108D, G116A and L120P substitutions), showed an unfolding temperature more than 10°C higher than that of Bac1. Furthermore, a combination of the other nine amino acid substitutions also led to improved thermostability of Bac1, although the effects of individual substitutions were small. Therefore, not only the sum of the contributions of individual amino acids, but also the synergistic effects of multiple amino acids are deeply involved in the stability of a hyperthermostable protein. Such insights will be helpful for future rational design of hyperthermostable proteins.

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.

Generation of new cytotoxic human ribonuclease variants directed to the nucleus

Ribonucleases are promising agents for use in anticancer therapy. Engineering a nuclear localization signal into the sequence of the human pancreatic ribonuclease has been revealed as a new strategy to endow this enzyme with cytotoxic activity against tumor cells. We previously described a cytotoxic human pancreatic ribonuclease variant, named PE5, which is able to cleave nuclear RNA, inducing the apoptosis of cancer cells and reducing the amount of P-glycoprotein in different multidrug-resistant cell lines. These results open the opportunity to use this ribonuclease in combination with other chemotherapeutics. In this work, we have investigated how to improve the properties of PE5 as an antitumor drug candidate. When attempting to develop a recombinant protein as a drug, two of the main desirable attributes are minimum immunogenicity and maximum potency. The improvements of PE5 have been designed in both senses. First, in order to reduce the potential immunogenicity of the protein, we have studied which residues mutated on PE5 can be reverted to those of the wild-type human pancreatic ribonuclease sequence without affecting its cytotoxicity. Second, we have investigated the effect of introducing an additional nuclear localization signal at different sites of PE5 in an effort to obtain a more cytotoxic enzyme. We show that the nuclear localization signal location is critical for the cytotoxicity. One of these variants, named NLSPE5, presents about a 10-fold increase in cytotoxicity respective to PE5. This variant induces apoptosis and kills the cells using the same mechanism as PE5.

The nuclear transport capacity of a human-pancreatic ribonuclease variant is critical for its cytotoxicity

We have previously described a human pancreatic-ribonuclease variant, named PE5, which carries a non-contiguous extended bipartite nuclear localization signal. This signal comprises residues from at least three regions of the protein. We postulated that the introduction of this signal in the ribonuclease provides it with cytotoxic activity because although the variant poorly evades the ribonuclease inhibitor in vitro, it is routed to the nucleus, which is devoid of the inhibitor. In this work, we have investigated the relationship between the cytotoxicity produced by PE5 and its ability to reach the nucleus. First, we show that this enzyme, when incubated with HeLa cells, specifically cleaves nuclear RNA while it leaves cytoplasmic RNA unaffected. On the other hand, we have created new variants in which the residues of the nuclear localization signal that are important for the nuclear transport have been replaced. As expected, the individual changes produce a significant decrease in the cytotoxicity of the resulting variants. We conclude that the nuclear transport of PE5 is critical for its cytotoxicity. Therefore, routing a ribonuclease to the nucleus is an alternative strategy to endow it with cytotoxic activity.

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.

Mapping the stability clusters in bovine pancreatic ribonuclease A

In the present work, we have thermodynamically characterized the thermally induced unfolding of 20 variants of bovine pancreatic ribonuclease A (RNase A) to experimentally describe the residues and the regions that are critical for the stability of the enzyme. The achieved results, complemented with previous studies by our group, allowed us to define the significance of the two hydrophobic nuclei present in the RNase A structure, as well as the contribution of the participating residues within each nucleus, to the global enzyme stability. We propose a structural model for the major and the minor hydrophobic nuclei of RNase A. The major nucleus is composite and located in the cavity delimited by alpha-helices 1 and 3, and the beta-sheet that is formed by strands 2, 3, 5, and 6. It consists of a central tight packed part constituted by residues Phe8, Met13, Val54, Val57, Ile106, Val108, and Phe120. This central part is surrounded by a layer formed by residues Val63, Tyr73, Met79, Ile107, Val116, and Val118. The minor nucleus, although less complex, is also constituted by a tight packing that involves the side chains of residues Tyr25, Met29, Met30, Leu35, Phe46, and Tyr97, which fill the cavity that originates the beta-sheet formed by beta-strands 1, 4, and 5 together with alpha-helix2.

Contribution of the C30/C75 disulfide bond to the biological properties of onconase

Onconase, a member of the pancreatic type ribonuclease family, is currently used as a chemotherapeutic agent for the treatment of different types of cancer. It is widely accepted that one of the properties that renders this enzyme cytotoxic is its ability to evade the cytosolic ribonuclease inhibitor (RI). In the present work, we produced and characterized an onconase variant that lacks the disulfide bond C30/C75. This variant mimics the stable unfolding intermediate des(30–75) produced in the reductive unfolding pathway of onconase. We found that the reduction of the C30/C75 disulfide bond does not significantly alter the cytotoxic properties of onconase, although the variant possesses a notably reduced conformational stability. Interestingly, both its catalytic activity and its ability to evade RI are comparable to wild-type onconase under mild reductive conditions in which the three disulfide containing intermediate des(30–75) is present. These results suggest that the C30/C75 disulfide bond could easily be reduced under physiological redox conditions.

Destabilizing mutations alter the hydrogen exchange mechanism in ribonuclease A

The effect of strongly destabilizing mutations, I106A and V108G of Ribonuclease A (RNase A), on its structure and stability has been determined by NMR. The solution structures of these variants are essentially equivalent to RNase A. The exchange rates of the most protected amide protons in RNase A (35 degrees C), the I106A variant (35 degrees C), and the V108G variant (10 degrees C) yield stability values of 9.9, 6.0, and 6.8 kcal/mol, respectively, when analyzed assuming an EX2 exchange mechanism. Thus, the destabilization induced by these mutations is propagated throughout the protein. Simulation of RNase A hydrogen exchange indicates that the most protected protons in RNase A and the V108G variant exchange via the EX2 regime, whereas those of I106A exchange through a mixed EX1 + EX2 process. It is striking that a single point mutation can alter the overall exchange mechanism. Thus, destabilizing mutations joins high temperatures, high pH and the presence of denaturating agents as a factor that induces EX1 exchange in proteins. The calculations also indicate a shift from the EX2 to the EX1 mechanism for less protected groups within the same protein. This should be borne in mind when interpreting exchange data as a measure of local stability in less protected regions.

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 contribution of the residues from the main hydrophobic core of ribonuclease A to its pressure-folding transition state

The role of hydrophobic interactions established by the residues that belong to the main hydrophobic core of ribonuclease A in its pressure-folding transition state was investigated using the Phi-value method. The folding kinetics was studied using pressure-jump techniques both in the pressurization and depressurization directions. The ratio between the folding activation volume and the reaction volume (beta p-value), which is an index of the compactness or degree of solvation of the transition state, was calculated. All the positions analyzed presented fractional Phi f-values, and the lowest were those corresponding to the most critical positions for the ribonuclease A stability. The structure of the transition state of the hydrophobic core of ribonuclease A, from the point of view of formed interactions, is a relatively, uniformly expanded form of the folded structure with a mean Phi f-value of 0.43. This places it halfway between the folded and unfolded states. On the other hand, for the variants, the average of beta p-values is 0.4, suggesting a transition state that is 40% native-like. Altogether the results suggest that the pressure-folding transition state of ribonuclease A looks like a collapsed globule with some secondary structure and a weakened hydrophobic core. A good correlation was found between the Phi f-values and the Deltabeta p-values. Although the nature of the transition state inferred from pressure-induced folding studies and the results of the protein engineering method have been reported to be consistent for other proteins, to the best of our knowledge this is the first direct comparison using a set of mutants.

Pressure as a tool to study protein-unfolding/refolding processes: The case of ribonuclease A

This paper gives an overview of the application of high-pressure to study the folding/unfolding processes of proteins using Ribonuclease A as a model protein. A particular focus is the study of pressure-equilibrium unfolding and folding kinetics using variants and the information obtained by comparing these with the wild-type enzyme.