The contribution of the residues from the main hydrophobic core of ribonuclease A to its pressure-folding transition state

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

Volume and compressibility differences between protein conformations revealed by high-pressure NMR

Proteins often interconvert between different conformations in ways critical to their function. Although manipulating such equilibria for biophysical study is often challenging, the application of pressure is a potential route to achieve such control by favoring the population of lower volume states. Here, we use this feature to study the interconversion of ARNT PAS-B Y456T, which undergoes a dramatic +3 slip in the β-strand register as it switches between two stably folded conformations. Using high-pressure biomolecular NMR approaches, we obtained the first, to our knowledge, quantitative data testing two key hypotheses of this process: the slipped conformation is both smaller and less compressible than the wild-type equivalent, and the interconversion proceeds through a chiefly unfolded intermediate state. Data collected in steady-state pressure and time-resolved pressure-jump modes, including observed pressure-dependent changes in the populations of the two conformers and increased rate of interconversion between conformers, support both hypotheses. Our work exemplifies how these approaches, which can be generally applied to protein conformational switches, can provide unique information that is not easily accessible through other techniques.

Construction of Highly Stable Cytotoxic Nuclear-Directed Ribonucleases

Ribonucleases are proteins whose use is promising in anticancer therapy. We have previously constructed different human pancreatic ribonuclease variants that are selectively cytotoxic for tumor cells by introducing a nuclear localization signal into their sequence. However, these modifications produced an important decrease in their stability compromising their behavior in vivo. Here, we show that we can significantly increase the thermal stability of these cytotoxic proteins by introducing additional disulfide bonds by site-directed mutagenesis. One of these variants increases its thermal stability by around 17 °C, without affecting its catalytic activity while maintaining the cytotoxic activity against tumor cells. We also show that the most stable variant is significantly more resistant to proteolysis when incubated with proteinase K or with human sera, suggesting that its half-live could be increased in vivo once administered.

Insights into the mechanism of Apoptin's exquisitely selective anti-tumor action from atomic level characterization of its conformation and dynamics

Apoptin is a 121 residue protein which forms large, soluble aggregates and possesses an exceptionally selectively cytotoxic action on cancer cells. In the accompanying paper, we described the design, production and initial characterization of an Apoptin truncated variant called H₆-ApopΔProΔLeu. Whereas both the variant and wild type protein possess similar selective cytotoxicity against cancer cells following transfection, only the variant is cytotoxic when added externally. Remarkably, as observed by gel filtration chromatography and dynamic light scattering, H₆-ApopΔProΔLeu lacks the tendency of wild type Apoptin to form large aggregates, which greatly facilitated the study of its biological properties. Here, we characterize the conformation and dynamics of H₆-ApopΔProΔLeu. Using a battery of 2D, 3D and (4,2)D NMR spectra, the essentially complete ¹H, ¹³C and ¹⁵N resonance assignments of H₆-ApopΔProΔLeu were obtained. The analysis of these data shows that the variant is an intrinsically disordered protein, which lacks a preferred conformation. This conclusion is corroborated by a lack of protection against proteolytic cleavage and hydrogen/deuterium exchange. Moreover, the CD spectra are dominated by random coil contributions. Finally, ¹H-¹⁵N NOE ratios are low, which indicates flexibility on the ps-ns time scale. Interestingly, H₆-ApopΔProΔLeu's intrinsically disordered ensemble is not significantly altered by the redox state of its Cys residues or by Thr phosphorylation, which has been proposed to play a key role in Apoptin's selective cytotoxicity. These results serve to better comprehend Apoptin's remarkably selective anticancer action and provide a framework for the future design of improved Apoptin variants.

Inhibitor and substrate binding induced stability of HIV-1 protease against sequential dissociation and unfolding revealed by high pressure spectroscopy and kinetics

High-pressure methods have become an interesting tool of investigation of structural stability of proteins. They are used to study protein unfolding, but dissociation of oligomeric proteins can be addressed this way, too. HIV-1 protease, although an interesting object of biophysical experiments, has not been studied at high pressure yet. In this study HIV-1 protease is investigated by high pressure (up to 600 MPa) fluorescence spectroscopy of either the inherent tryptophan residues or external 8-anilino-1-naphtalenesulfonic acid at 25°C. A fast concentration-dependent structural transition is detected that corresponds to the dimer-monomer equilibrium. This transition is followed by a slow concentration independent transition that can be assigned to the monomer unfolding. In the presence of a tight-binding inhibitor none of these transitions are observed, which confirms the stabilizing effect of inhibitor. High-pressure enzyme kinetics (up to 350 MPa) also reveals the stabilizing effect of substrate. Unfolding of the protease can thus proceed only from the monomeric state after dimer dissociation and is unfavourable at atmospheric pressure. Dimer-destabilizing effect of high pressure is caused by negative volume change of dimer dissociation of -32.5 mL/mol. It helps us to determine the atmospheric pressure dimerization constant of 0.92 μM. High-pressure methods thus enable the investigation of structural phenomena that are difficult or impossible to measure at atmospheric pressure.

Investigation of an anomalously accelerating substitution in the folding of a prototypical two-state protein

The folding rates of two-state single-domain proteins are generally resistant to small-scale changes in amino acid sequence. For example, having surveyed here over 700 single-residue substitutions in 24 well-characterized two-state proteins, we find that the majority (55%) of these substitutions affect folding rates by less than a factor of 2, and that only 9% affect folding rates by more than a factor of 8. Among those substitutions that significantly affect folding rates, we find that accelerating substitutions are an order of magnitude less common than those that decelerate the process. One of the most extreme outliers in this data set, an arginine-to-phenylalanine substitution at position 48 (R48F) of chymotrypsin inhibitor 2 (CI2), accelerates the protein's folding rate by a factor of 36 relative to that of the wild-type protein and is the most accelerating substitution reported to date in a two-state protein. In order to better understand the origins of this anomalous behavior, we have characterized the kinetics of multiple additional substitutions at this position. We find that substitutions at position 48 in CI2 fall into two distinct classes. The first, comprising residues that ablate the charge of the wild-type arginine but retain the hydrophobicity of its alkane chain, accelerate folding by at least 10-fold. The second class, comprising all other residues, produces folding rates within a factor of two of the wild-type rate. A significant positive correlation between hydrophobicity and folding rate across all of the residues we have characterized at this position suggests that the hydrophobic methylene units of the wild-type arginine play a significant role in stabilizing the folding transition state. Likewise, studies of the pH dependence of the histidine substitution indicate a strong correlation between folding rate and charge state. Thus, mutations that ablate the arginine's positive charge while retaining the hydrophobic contacts of its methylene units tend to dramatically accelerate folding. Previous studies have suggested that arginine 48 plays an important functional role in CI2, which may explain why it is highly conserved despite the anomalously large deceleration it produces in the folding of this protein.

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.

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.

Predicting protein folding cores by empirical potential functions

Theoretical and in vitro experiments suggest that protein folding cores form early in the process of folding, and that proteins may have evolved to optimize both folding speed and native-state stability. In our previous work (Chen et al., Structure, 14 (2006) 1401), we developed a set of empirical potential functions and used them to analyze interaction energies among secondary-structure elements in two beta-sandwich proteins. Our work on this group of proteins demonstrated that the predicted folding core also harbors residues that form native-like interactions early in the folding reaction. In the current work, we have tested our empirical potential functions on structurally-different proteins for which the folding cores have been revealed by protein hydrogen-deuterium exchange experiments. Using a set of 29 unrelated proteins, which have been extensively studied in the literature, we demonstrate that the average prediction result from our method is significantly better than predictions based on other computational methods. Our study is an important step towards the ultimate goal of understanding the correlation between folding cores and native structures.

How general is the nucleation-condensation mechanism?

We investigate the structures of the major folding transition states of nine proteins by correlation of published Phi-values with inter-residue contact maps. Combined with previous studies on six proteins, the analysis suggests that at least 10 of the 15 small globular proteins fold via a nucleation-condensation mechanism with a concurrent build-up of secondary and tertiary structure contacts, but a structural consolidation that is clearly nonuniformly distributed over the molecule and most intense in a single structural region suggesting the occurrence of a single folding nucleus. However, on average helix- and sheet-forming residues show somewhat larger Phi-values in the major transition state, suggesting that secondary structure formation is one important driving force in the nucleation-condensation in many proteins and that secondary-structure forming residues tend to be more prominent in folding nuclei. We synthesize the combined information on these 10 of 15 proteins into a unified nucleation-condensation mechanism which also accounts for effects described by the framework, hydrophobic collapse, zipper, and funnel models.

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.