raf RBD and Ubiquitin Proteins Share Similar Folds, Folding Rates and Mechanisms Despite Having Unrelated Amino Acid Sequences

Specific inhibition of oncogenic RAS using cell-permeable RAS-binding domains

Oncogenic RAS proteins, common oncogenic drivers in many human cancers, have been refractory to conventional small-molecule and macromolecule inhibitors due to their intracellular localization and the lack of druggable pockets. Here, we present a feasible strategy for designing RAS inhibitors that involves intracellular delivery of RAS-binding domain (RBD), a nanomolar-affinity specific ligand of RAS. Screening of 51 different combinations of RBD and cell-permeable peptides has identified Pen-cRaf-v1 as a cell-permeable pan-RAS inhibitor capable of targeting both G12C and non-G12C RAS mutants. Pen-cRaf-v1 crosses the cell membrane via endocytosis, competitively inhibits RAS-effector interaction, and thereby exerts anticancer activity against several KRAS-mutant cancer cell lines. Moreover, Pen-cRaf-v1 exhibits excellent activity comparable with a leading pan-RAS inhibitor (BI-2852), as well as high target specificity in transcriptome analysis and alanine mutation analysis. These findings demonstrate that specific inhibition of oncogenic RAS, and possibly treatment of RAS-mutant cancer, is feasible by intracellular delivery of RBD.

PFDB: A standardized protein folding database with temperature correction

We constructed a standardized protein folding kinetics database (PFDB) in which the logarithmic rate constants of all listed proteins are calculated at the standard temperature (25 °C). A temperature correction based on the Eyring–Kramers equation was introduced for proteins whose folding kinetics were originally measured at temperatures other than 25 °C. We verified the temperature correction by comparing the logarithmic rate constants predicted and experimentally observed at 25 °C for 14 different proteins, and the results demonstrated improvement of the quality of the database. PFDB consists of 141 (89 two-state and 52 non-two-state) single-domain globular proteins, which has the largest number among the currently available databases of protein folding kinetics. PFDB is thus intended to be used as a standard for developing and testing future predictive and theoretical studies of protein folding. PFDB can be accessed from the following link: http://lee.kias.re.kr/~bala/PFDB.

Visualizing transient protein-folding intermediates by tryptophan-scanning mutagenesis

To understand how proteins fold, assemble and function, it is necessary to characterize the structure and dynamics of each state they adopt during their lifetime. Experimental characterization of the transient states of proteins remains a major challenge because high-resolution structural techniques, including NMR and X-ray crystallography, cannot be directly applied to study short-lived protein states. To circumvent this limitation, we show that transient states during protein folding can be characterized by measuring the fluorescence of tryptophan residues, introduced at many solvent-exposed positions to determine whether each position is native-like, denatured-like or non-native-like in the intermediate state. We use this approach to characterize a late-folding-intermediate state of the small globular mammalian protein ubiquitin, and we show the presence of productive non-native interactions that suggest a 'flycatcher' mechanism of concerted binding and folding.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

Helix mutations stabilize a late productive intermediate on the folding pathway of ubiquitin

We have investigated the relative placement of rate-limiting energy barriers and the role of productive or obstructive intermediates on the folding pathway of yeast wild-type ubiquitin ( wt-Ub) containing the F45W mutation. To manipulate the folding barriers, we have designed a family of mutants in which stabilizing substitutions have been introduced incrementally on the solvent-exposed surface of the main alpha-helix (residues 23-34), which has a low intrinsic helical propensity in the native sequence. Although the U --> I and I --> N transitions are not clearly delineated in the kinetics of wt-Ub, we show that an intermediate becomes highly populated and more clearly resolved as the predicted stability of the helix increases. The observed acceleration in the rate of folding correlates with helix stability and is consistent with the I-state representing a productive rather than misfolded state. A Leffler analysis of the effects on kinetics of changes in stability within the family of helix mutants results in a biphasic correlation in both the refolding and unfolding rates that suggest a shift from a nucleation-condensation mechanism (weakly stabilized helix) toward a diffusion-collision model (highly stabilized helix). Through the introduction of helix-stabilizing mutations, we are able to engineer a well-resolved I-state on the folding pathway of ubiquitin which is likely to be structurally distinct from that which is only weakly populated on the folding pathway of wild-type ubiquitin.

Multiple tryptophan probes reveal that ubiquitin folds via a late misfolded intermediate

Much of our understanding of protein folding mechanisms is derived from experiments using intrinsic fluorescence of natural or genetically inserted tryptophan (Trp) residues to monitor protein refolding and site-directed mutagenesis to determine the energetic role of amino acids in the native (N), intermediate (I) or transition (T) states. However, this strategy has limited use to study complex folding reactions because a single fluorescence probe may not detect all low-energy folding intermediates. To overcome this limitation, we suggest that protein refolding should be monitored with different solvent-exposed Trp probes. Here, we demonstrate the utility of this approach by investigating the controversial folding mechanism of ubiquitin (Ub) using Trp probes located at residue positions 1, 28, 45, 57, and 66. We first show that these Trp are structurally sensitive and minimally perturbing fluorescent probes for monitoring folding/unfolding of the protein. Using a conventional stopped-flow instrument, we show that ANS and Trp fluorescence detect two distinct transitions during the refolding of all five Trp mutants at low concentrations of denaturant: T(1), a denaturant-dependent transition and T(2), a slower transition, largely denaturant-independent. Surprisingly, some Trp mutants (Ub(M1W), Ub(S57W)) display Trp fluorescence changes during T(1) that are distinct from the expected U-->N transition suggesting that the denaturant-dependent refolding transition of Ub is not a U-->N transition but represents the formation of a structurally distinct I-state (U-->I). Alternatively, this U-->I transition could be also clearly distinguished by using a combination of two Trp mutations Ub(F45W-T66W) for which the two Trp probes that display fluorescence changes of opposite sign during T(1) and T(2) (Ub(F45W-T66W)). Global fitting of the folding/unfolding kinetic parameters and additional folding-unfolding double-jump experiments performed on Ub(M1W), a mutant with enhanced fluorescence in the I-state, demonstrate that the I-state is stable, compact, misfolded, and on-pathway. These results illustrate how transient low-energy I-states can be characterized efficiently in complex refolding reactions using multiple Trp probes.

Structural and mechanistic exploration of acid resistance: kinetic stability facilitates evolution of extremophilic behavior

Kinetically stable proteins are unique in that their stability is determined solely by kinetic barriers rather than by thermodynamic equilibria. To better understand how kinetic stability promotes protein survival under extreme environmental conditions, we analyzed the unfolding behavior and determined the structure of Nocardiopsis alba Protease A (NAPase), an acid-resistant, kinetically stable protease, and compared these results with a neutrophilic homolog, alpha-lytic protease (alphaLP). Although NAPase and alphaLP have the same number of acid-titratable residues, kinetic studies revealed that the height of the unfolding free energy barrier for NAPase is less sensitive to acid than that of alphaLP, thereby accounting for NAPase's improved tolerance of low pH. A comparison of the alphaLP and NAPase structures identified multiple salt-bridges in the domain interface of alphaLP that were relocated to outer regions of NAPase, suggesting a novel mechanism of acid stability in which acid-sensitive electrostatic interactions are rearranged to similarly affect the energetics of both the native state and the unfolding transition state. An acid-stable variant of alphaLP in which a single interdomain salt-bridge is replaced with a corresponding intradomain NAPase salt-bridge shows a dramatic >15-fold increase in acid resistance, providing further evidence for this hypothesis. These observations also led to a general model of the unfolding transition state structure for alphaLP protease family members in which the two domains separate from each other while remaining relatively intact themselves. These results illustrate the remarkable utility of kinetic stability as an evolutionary tool for developing longevity over a broad range of harsh conditions.

The transition state of the ras binding domain of Raf is structurally polarized based on Phi-values but is energetically diffuse

The ras binding domain (RBD) of the Ser/Thr kinase c-Raf/Raf-1 spans 78 residues and adopts a structure characteristic of the beta-grasp ubiquitin-like topology. Recently, the primary sequence of Raf RBD has been nearly exhaustively mutated experimentally by insertion of stretches of degenerate codons, which revealed sequence conservation and hydrophobic core organization similar to that found in an alignment of beta-grasp ubiquitin-like proteins. These results now allow us to examine the relationship between sequence conservation and the folding process, particularly viewed through the analysis of transition state (TS) structure. Specifically, we present herein a protein engineering study combining classic truncation (Ala/Gly) and atypical mutants to predict folding TS ensemble properties. Based on classical Phi-value analysis, Raf RBD TS structure is particularly polarized around the N-terminal beta-hairpin. However, all residues constituting the inner layer of the hydrophobic core are involved in TS stabilization, although they are clearly found in a less native-like environment. The TS structure can also be probed by a direct measure of its destabilization upon mutation, DeltaDeltaG(U-++). Viewed through this analysis, Raf RBD TS is a more diffuse structure, in which all residues of the hydrophobic core including beta-strands 1, 2, 3 and 5 and the major alpha-helix play similar roles in TS stabilization. In addition, Phi-values and DeltaDeltaG(U-++) reveal striking similarities in the TS of Raf RBD and ubiquitin, a structural analogue displaying insignificant sequence identity (<12%). However, ubiquitin TS appears more denatured-like and polarized around the N-terminal beta-hairpin. We suggest that analysis of Phi-values should also consider the direct impact of mutations on differences in free energy between the unfolded and TS (DeltaDeltaG(U-++)) to ensure that the description of TS properties is accurate. Finally, the impact of these findings on the modeling of protein folding is discussed.

Folding of S6 structures with divergent amino acid composition: pathway flexibility within partly overlapping foldons

Studies of circular permutants have demonstrated that the folding reaction of S6 from Thermus thermophilus (S6(T)) is malleable and responds in an ordered manner to changes of the sequence separation between interacting residues: the S6(T) permutants retain a common nucleation pattern in the form of a two-strand-helix motif that can be recruited from different parts of the structure. To further test the robustness of the two-strand-helix nucleus we have here determined the crystal structure and folding reaction of an evolutionary divergent S6 protein from the hyperthermophilic bacterium Aquifex aeolicus (S6(A)). Although the overall topology of S6(A) is very similar to that of S6(T) the architecture of the hydrophobic core is radically different by containing a large proportion of stacked Phe side-chains. Despite this disparate core composition, the folding rate constant and the kinetic m values of S6(A) are identical to those of S6(T). The folding nucleus of S6(A) is also found to retain the characteristic two-strand-helix motif of the S6(T) permutants, but with a new structural emphasis. The results suggest that the protein folding reaction is linked to topology only in the sense that the native-state topology determines the repertoire of accessible nucleation motifs. If the native structure allows several equivalent ways of recruiting a productive nucleus the folding reaction is free to redistribute within these topological constraints.

Massive sequence perturbation of the Raf ras binding domain reveals relationships between sequence conservation, secondary structure propensity, hydrophobic core organization and stability

The contributions of specific residues to the delicate balance between function, stability and folding rates could be determined, in part by [corrected] comparing the sequences of structures having identical folds, but insignificant sequence homology. Recently, we have devised an experimental strategy to thoroughly explore residue substitutions consistent with a specific class of structure. Using this approach, the amino acids tolerated at virtually all residues of the c-Raf/Raf1 ras binding domain (Raf RBD), an exemplar of the common beta-grasp ubiquitin-like topology, were obtained and used to define the sequence determinants of this fold. Herein, we present analyses suggesting that more subtle sequence selection pressure, including propensity for secondary structure, the hydrophobic core organization and charge distribution are imposed on the Raf RBD sequence. Secondly, using the Gibbs free energies (DeltaG(F-U)) obtained for 51 mutants of Raf RBD, we demonstrate a strong correlation between amino acid conservation and the destabilization induced by truncating mutants. In addition, four mutants are shown to significantly stabilize Raf RBD native structure. Two of these mutations, including the well-studied R89L, are known to severely compromise binding affinity for ras. Another stabilized mutant consisted of a deletion of amino acid residues E104-K106. This deletion naturally occurs in the homologues a-Raf and b-Raf and could indicate functional divergence. Finally, the combination of mutations affecting five of 78 residues of Raf RBD results in stabilization of the structure by approximately 12 kJ mol(-1) (DeltaG(F-U) is -22 and -34 kJ mol(-1) for wt and mutant, respectively). The sequence perturbation approach combined with sequence/structure analysis of the ubiquitin-like fold provide a basis for the identification of sequence-specific requirements for function, stability and folding rate of the Raf RBD and structural analogues, highlighting the utility of conservation profiles as predictive tools of structural organization.

Population of On-pathway Intermediates in the Folding of Ubiquitin

The role that intermediate states play in protein folding is the subject of intense investigation and in the case of ubiquitin has been controversial. We present fluorescence-detected kinetic data derived from single and double mixing stopped-flow experiments to show that the F45W mutant of ubiquitin (WT*), a well-studied single-domain protein and most recently regarded as a simple two-state system, folds via on-pathway intermediates. To account for the discrepancy we observe between equilibrium and kinetic stabilities and m-values, we show that the polypeptide chain undergoes rapid collapse to an intermediate whose presence we infer from a fast lag phase in interrupted refolding experiments. Double-jump kinetic experiments identify two direct folding phases that are not associated with slow isomerisation reactions in the unfolded state. These two phases are explained by kinetic partitioning which allows molecules to reach the native state from the collapsed state via two possible competing routes, which we further examine using two destabilised ubiquitin mutants. Interrupted refolding experiments allow us to observe the formation and decay of an intermediate along one of these pathways. A plausible model for the folding pathway of ubiquitin is presented that demonstrates that obligatory intermediates and/or chain collapse are important events in restricting the conformational search for the native state of ubiquitin.

Folding is coupled to dimerization of Tctex-1 dynein light chain

Equilibrium analyses have been performed to elucidate the role of dimerization in folding and stability of dynein light chain Tctex-1. The equilibrium unfolding transition was monitored by intrinsic fluorescence intensity, fluorescence anisotropy, and circular dichroism and was modeled as a two-state mechanism where a folded dimer dissociates to two unfolded monomers without populating thermodynamically stable monomeric or dimeric intermediates. Sedimentation equilibrium and chemical cross-linking experiments performed at increasing concentrations of denaturants show no change in the association state before the unfolding transition and are consistent with the two-state model of dissociation coupled to unfolding. A linear dependence on denaturant concentration is observed by fluorescence intensity and anisotropy before unfolding in the 0-2 M GdnCl, and 0-4 M urea concentration range. This change is not protein concentration-dependent and possibly reflects relief of quenching associated with premelting conformational disorder in the vicinity of Trp 83. The data clearly show that the dissociation-coupled unfolding mechanism of Tctex-1 is different from the three-state denaturation mechanism of its structural homologue light chain LC8. The absence of a stable monomer in Tctex-1 offers insight into its functional differences from LC8.

Ubiquitin: a small protein folding paradigm

For the past twenty years, the small, 76-residue protein ubiquitin has been used as a model system to study protein structure, stability, folding and dynamics. In this time, ubiquitin has become a paradigm for both the experimental and computational folding communities. The folding energy landscape is now uniquely characterised with a plethora of information available on not only the native and denatured states, but partially structured states, alternatively folded states and locally unfolded states, in addition to the transition state ensemble. This Perspective focuses on the experimental characterisation of ubiquitin using a comprehensive range of biophysical techniques.

Massive sequence perturbation of a small protein

Most protein topologies rarely occur in nature, thus limiting our ability to extract sequence information that could be used to predict structure, function, and evolutionary constraints on protein folds. In principle, the sequence diversity explored by a given protein topology could be expanded by introducing sequence perturbations and selecting variant proteins that fold correctly. However, our capacity to explore sequence space is intrinsically limited by the enormous number of sequences generated from the 20 amino acids and the limited number of variants likely to fold. Here we sought to test whether the sequence space for naturally existing proteins can be explored by simple, sequential degeneration of a complete set of short sequence segments of a model protein, without long-range covariation. Using the Raf ras binding domain as a model of a small protein capable of autonomous folding, we degenerated 72 of 76 positions of the primary structure for the 20 amino acids in segments of four to seven residues defined by secondary structure and selected the folded species for interaction with h-ras by using an in vivo survival-selection assay. The methodology presented allowed for rigorous statistical analysis and comparison of sequence diversity. The ensemble of sequence variants of Raf ras binding domain obtained have recaptured the diversity observed for the ubiquitin-roll topology. A signature sequence for this fold and the implication of this strategy to protein design and structure prediction are discussed.

Protein folding: defining a "standard" set of experimental conditions and a preliminary kinetic data set of two-state proteins

Recent years have seen the publication of both empirical and theoretical relationships predicting the rates with which proteins fold. Our ability to test and refine these relationships has been limited, however, by a variety of difficulties associated with the comparison of folding and unfolding rates, thermodynamics, and structure across diverse sets of proteins. These difficulties include the wide, potentially confounding range of experimental conditions and methods employed to date and the difficulty of obtaining correct and complete sequence and structural details for the characterized constructs. The lack of a single approach to data analysis and error estimation, or even of a common set of units and reporting standards, further hinders comparative studies of folding. In an effort to overcome these problems, we define here a "consensus" set of experimental conditions (25 degrees C at pH 7.0, 50 mM buffer), data analysis methods, and data reporting standards that we hope will provide a benchmark for experimental studies. We take the first step in this initiative by describing the folding kinetics of 30 apparently two-state proteins or protein domains under the consensus conditions. The goal of our efforts is to set uniform standards for the experimental community and to initiate an accumulating, self-consistent data set that will aid ongoing efforts to understand the folding process.

Two different T cell receptors use different thermodynamic strategies to recognize the same peptide/MHC ligand

A6 and B7 are two alphabeta T cell receptors (TCRs) that recognize the Tax peptide presented by the class I major histocompatibility molecule HLA-A2 (Tax/HLA-A2). Despite the fact that the two TCRs have different CDR loops and use different amino acid residues to contact their ligand, both receptors bind ligand with similar diagonal orientations. Here we show that they also bind with very similar binding affinities and kinetics (the DeltaDeltaG degrees for binding is approximately 0.3kcal/mol at 25 degrees C). The two receptors respond similarly to alterations in the MHC molecule, yet differ dramatically in their responses to ionic strength and temperature. The different responses to temperature indicate markedly different binding thermodynamics, which are not predictable from the surface area buried in the interfaces. A6 and B7 thus represent two TCRs that are both compatible with Tax/HLA-A2, although compatibility has been achieved through the use of different thermodynamic strategies. Finally, neither A6 nor B7 are predicted to undergo large conformational adaptations upon binding, distinguishing them from a number of other TCRs whose structure, thermodynamics, and kinetics have been characterized.

Molecular characterization of the 30-AA N-terminal mineral interaction domain of the biomineralization protein AP7

The AP7 protein is one of several mollusk shell proteins which are responsible for aragonite polymorph formation and stabilization within the nacre layer of the Pacific red abalone, H. rufescens. Previously, we demonstrated that the 30-AA N-terminal domain of AP7, denoted as AP7-1, exists as an unfolded sequence and possesses the capability of inhibiting calcium carbonate crystal growth in vitro via growth step frustration or interruption. However, very little is known with regard to the interactive capabilities of this sequence with Ca(II) and with calcium carbonates. Using multidisciplinary techniques, we determine that the AP7-1 polypeptide interacts with Ca(II) ions at the -DD- sequence clusters, yet retains its unfolded, conformationally labile structure in the presence of Ca(II) ions. Further, NMR experiments reveal that the extended structured sequence blocks, -GNGM-, -SVRTQG-, and -ISYL, exhibit motional, chemical exchange, and/or backbone geometry perturbations in response to Ca(II) interactions with AP7-1. Solid-state NMR magic angle spinning studies verify that during the course of in vitro calcium carbonate crystal growth, AP7-1 becomes bound to calcite fragments and cannot be entirely displaced from the mineral fragments using competitive Ca(II) washing. Finally, using a scrambled sequence version of the AP7-1 polypeptide, we observe that sequence scrambling does not adversely affect the crystal growth inhibitory activity of AP7-1, suggesting that the amino acid composition of AP7-1 may be more critical to growth step inhibition than the linear ordering of amino acids.