Reversible switching between two common protein folds in a designed system using only temperature.

Naturally occurring metamorphic proteins have the ability to interconvert from one folded state to another through either a limited set of mutations or by way of a change in the local environment. Here, we show in a designed system that it is possible to switch reversibly between two of the most common monomeric folds employing only temperature changes. We demonstrate that a latent 3α state can be unmasked from an α/ß-plait topology with a single V90T amino acid substitution, populating both forms simultaneously. The equilibrium between these two states exhibits temperature dependence, such that the 3α state is predominant (>90%) at 5 °C, while the α/ß-plait fold is the major species (>90%) at 30 °C. We describe the structure and dynamics of these topologies, how mutational changes affect the temperature dependence, and the energetics and kinetics of interconversion. Additionally, we demonstrate how ligand-binding function can be tightly regulated by large amplitude changes in protein structure over a relatively narrow temperature range that is relevant to biology. The 3α/αß switch thus represents a potentially useful approach for designing proteins that alter their fold topologies in response to environmental triggers. It may also serve as a model for computational studies of temperature-dependent protein stability and fold switching.

Phosphorylation-induced Conformational Ensemble Switching in an Intrinsically Disordered Cancer/Testis Antigen.

Prostate-associated gene 4 (PAGE4) is an intrinsically disordered cancer/testis antigen that is up-regulated in the fetal and diseased human prostate. Knocking down PAGE4 expression results in cell death, whereas its overexpression leads to a growth advantage of prostate cancer cells (Zeng, Y., He, Y., Yang, F., Mooney, S. M., Getzenberg, R. H., Orban, J., and Kulkarni, P. (2011) The cancer/testis antigen prostate-associated gene 4 (PAGE4) is a highly intrinsically disordered protein. J. Biol. Chem. 286, 13985-13994). Phosphorylation of PAGE4 at Thr-51 is critical for potentiating c-Jun transactivation, an important factor in controlling cell growth, apoptosis, and stress response. Using NMR spectroscopy, we show that the PAGE4 polypeptide chain has local and long-range conformational preferences that are perturbed by site-specific phosphorylation at Thr-51. The population of transient turn-like structures increases upon phosphorylation in an ∼20-residue acidic region centered on Thr-51. This central region therefore becomes more compact and more negatively charged, with increasing intramolecular contacts to basic sequence motifs near the N and C termini. Although flexibility is decreased in the central region of phospho-PAGE4, the polypeptide chain remains highly dynamic overall. PAGE4 utilizes a transient helical structure adjacent to the central acidic region to bind c-Jun with low affinity in vitro. The binding interaction is attenuated by phosphorylation at Thr-51, most likely because of masking the effects of the more compact phosphorylated state. Therefore, phosphorylation of PAGE4 leads to conformational shifts in the dynamic ensemble, with large functional consequences. The changes in the structural ensemble induced by posttranslational modifications are similar conceptually to the conformational switching events seen in some marginally stable ("metamorphic") folded proteins in response to mutation or environmental triggers.

Subdomain interactions foster the design of two protein pairs with ∼80% sequence identity but different folds.

Metamorphic proteins, including proteins with high levels of sequence identity but different folds, are exceptions to the long-standing rule-of-thumb that proteins with as little as 30% sequence identity adopt the same fold. Which topologies can be bridged by these highly identical sequences remains an open question. Here we bridge two 3-α-helix bundle proteins with two radically different folds. Using a straightforward approach, we engineered the sequences of one subdomain within maltose binding protein (MBP, α/ß/α-sandwich) and another within outer surface protein A (OspA, ß-sheet) to have high sequence identity (80 and 77%, respectively) with engineered variants of protein G (GA, 3-α-helix bundle). Circular dichroism and nuclear magnetic resonance spectra of all engineered variants demonstrate that they maintain their native conformations despite substantial sequence modification. Furthermore, the MBP variant (80% identical to GA) remained active. Thermodynamic analysis of numerous GA and MBP variants suggests that the key to our approach involved stabilizing the modified MBP and OspA subdomains via external interactions with neighboring substructures, indicating that subdomain interactions can stabilize alternative folds over a broad range of sequence variation. These findings suggest that it is possible to bridge one fold with many other topologies, which has implications for protein folding, evolution, and misfolding diseases.

Ionizable lipid nanoparticles for RAS protease delivery to inhibit cancer cell proliferation.

Mutations in RAS, a family of proteins found in all human cells, drive a third of cancers, including many pancreatic, colorectal, and lung cancers. However, there is a lack of clinical therapies that can effectively prevent RAS from causing tumor growth. Recently, a protease was engineered that specifically degrades active RAS, offering a promising new tool for treating these cancers. However, like many other intracellularly acting protein-based therapies, this protease requires a delivery vector to reach its site of action within the cell. In this study, we explored the incorporation of cationic lipids into ionizable lipid nanoparticles (LNPs) to develop a RAS protease delivery platform capable of inhibiting cancer cell proliferation in vitro and in vivo. A library of 13 LNPs encapsulating RAS protease was designed, and each formulation was evaluated for in vitro delivery efficiency and toxicity. A subset of four top-performing LNP formulations was identified and further evaluated for their impact on cancer cell proliferation in human colorectal cancer cells with mutated KRAS in vitro and in vivo, as well as their in vivo biodistribution and toxicity. In vivo, both the concentration of cationic lipid and type of cargo influenced LNP and cargo distribution. All lead candidate LNPs showed RAS protease functionality in vitro, and the top-performing formulation achieved effective intracellular RAS protease delivery in vivo, decreasing cancer cell proliferation in an in vivo xenograft model and significantly reducing tumor growth and size. Overall, this work demonstrates the use of LNPs as an effective delivery platform for RAS proteases, which could potentially be utilized for cancer therapies.

Design and characterization of a protein fold switching network.

To better understand how amino acid sequence encodes protein structure, we engineered mutational pathways that connect three common folds (3α, ß-grasp, and α/ß-plait). The structures of proteins at high sequence-identity intersections in the pathways (nodes) were determined using NMR spectroscopy and analyzed for stability and function. To generate nodes, the amino acid sequence encoding a smaller fold is embedded in the structure of an ~50% larger fold and a new sequence compatible with two sets of native interactions is designed. This generates protein pairs with a 3α or ß-grasp fold in the smaller form but an α/ß-plait fold in the larger form. Further, embedding smaller antagonistic folds creates critical states in the larger folds such that single amino acid substitutions can switch both their fold and function. The results help explain the underlying ambiguity in the protein folding code and show that new protein structures can evolve via abrupt fold switching.

The denatured state dictates the topology of two proteins with almost identical sequence but different native structure and function.

The protein folding problem is often studied by comparing the mechanisms of proteins sharing the same structure but different sequence. The recent design of the two proteins G(A)88 and G(B)88, displaying different structures and functions while sharing 88% sequence identity (49 out of 56 amino acids), allows the unique opportunity for a complementary approach. At which stage of its folding pathway does a protein commit to a given topology? Which residues are crucial in directing folding mechanisms to a given structure? By using a combination of biophysical and computational techniques, we have characterized the folding of both G(A)88 and G(B)88. We show that, contrary to expectation, G(B)88, characterized by a native α+ß fold, displays in the denatured state a content of native-like helical structure greater than G(A)88, which is all-α in its native state. Both experiments and simulations indicate that such residual structure may be tuned by changing pH. Thus, despite the high sequence identity, the folding pathways for these two proteins appear to diverge as early as in the denatured state. Our results suggest a mechanism whereby protein topology is committed very early along the folding pathway, being imprinted in the residual structure of the denatured state.

Solution NMR structure of a sheddase inhibitor prodomain from the malarial parasite Plasmodium falciparum.

Plasmodium subtilisin 2 (Sub2) is a multidomain protein that plays an important role in malaria infection. Here, we describe the solution NMR structure of a conserved region of the inhibitory prodomain of Sub2 from Plasmodium falciparum, termed prosub2. Despite the absence of any detectable sequence homology, the protozoan prosub2 has structural similarity to bacterial and mammalian subtilisin-like prodomains. Comparison with the three-dimensional structures of these other prodomains suggests a likely binding interface with the catalytic domain of Sub2 and provides insights into the locations of primary and secondary processing sites in Plasmodium prodomains.

A minimal sequence code for switching protein structure and function.

We present here a structural and mechanistic description of how a protein changes its fold and function, mutation by mutation. Our approach was to create 2 proteins that (i) are stably folded into 2 different folds, (ii) have 2 different functions, and (iii) are very similar in sequence. In this simplified sequence space we explore the mutational path from one fold to another. We show that an IgG-binding, 4beta+alpha fold can be transformed into an albumin-binding, 3-alpha fold via a mutational pathway in which neither function nor native structure is completely lost. The stabilities of all mutants along the pathway are evaluated, key high-resolution structures are determined by NMR, and an explanation of the switching mechanism is provided. We show that the conformational switch from 4beta+alpha to 3-alpha structure can occur via a single amino acid substitution. On one side of the switch point, the 4beta+alpha fold is >90% populated (pH 7.2, 20 degrees C). A single mutation switches the conformation to the 3-alpha fold, which is >90% populated (pH 7.2, 20 degrees C). We further show that a bifunctional protein exists at the switch point with affinity for both IgG and albumin.

NMR structures of two designed proteins with high sequence identity but different fold and function.

How protein sequence codes for 3D structure remains a fundamental question in biology. One approach to understanding the folding code is to design a pair of proteins with maximal sequence identity but retaining different folds. Therefore, the nonidentities must be responsible for determining which fold topology prevails and constitute a fold-specific folding code. We recently designed two proteins, G(A)88 and G(B)88, with 88% sequence identity but different folds and functions [Alexander et al. (2007) Proc Natl Acad Sci USA 104:11963-11968]. Here, we describe the detailed 3D structures of these proteins determined in solution by NMR spectroscopy. Despite a large number of mutations taking the sequence identity level from 16 to 88%, G(A)88 and G(B)88 maintain their distinct wild-type 3-alpha and alpha/beta folds, respectively. To our knowledge, the 3D-structure determination of two monomeric proteins with such high sequence identity but different fold topology is unprecedented. The geometries of the seven nonidentical residues (of 56 total) provide insights into the structural basis for switching between 3-alpha and alpha/beta conformations. Further mutation of a subset of these nonidentities, guided by the G(A)88 and G(B)88 structures, leads to proteins with even higher levels of sequence identity (95%) and different folds. Thus, conformational switching to an alternative monomeric fold of comparable stability can be effected with just a handful of mutations in a small protein. This result has implications for understanding not only the folding code but also the evolution of new folds.

Engineering subtilisin proteases that specifically degrade active RAS.

We describe the design, kinetic properties, and structures of engineered subtilisin proteases that degrade the active form of RAS by cleaving a conserved sequence in switch 2. RAS is a signaling protein that, when mutated, drives a third of human cancers. To generate high specificity for the RAS target sequence, the active site was modified to be dependent on a cofactor (imidazole or nitrite) and protease sub-sites were engineered to create a linkage between substrate and cofactor binding. Selective proteolysis of active RAS arises from a 2-step process wherein sub-site interactions promote productive binding of the cofactor, enabling cleavage. Proteases engineered in this way specifically cleave active RAS in vitro, deplete the level of RAS in a bacterial reporter system, and also degrade RAS in human cell culture. Although these proteases target active RAS, the underlying design principles are fundamental and will be adaptable to many target proteins.

Structure of a switchable subtilisin complexed with a substrate and with the activator azide.

An engineered variant of the protease subtilisin from Bacillus amyloliquefaciens, in which the D32A mutation renders the enzyme's activity dependent on the presence of certain small anions such as fluoride or azide, has been produced. This modified enzyme has applications as an azide or fluoride-triggered expression-purification tool. We report activity measurements showing that the enzyme is activated more than 3000-fold by azide and describe the 1.8 A resolution structure of an inactive form (by replacing the catalytic nucleophile Ser 221 with alanine) of the protease, in complex with azide and with a substrate that spans the active site. Both enzyme and substrate have been engineered to increase their stability and the affinity of their interaction. The substrate is based on a stabilized subtilisin prodomain, extended across the active site by the addition of four residues at its C-terminus. In the crystal structure, the substrate is well-ordered across the active site, and the azide anion is observed bound adjacent to Ala 32. The structures of the substrate complex in three different crystals (anion-free, fluoride-soaked, and azide-soaked) are compared. These structures provide extensive information for understanding subtilisin's substrate binding and catalytic mechanism, and for the development of biotechnology tools based on anion-activated proteolysis. The mechanism of anion-dependent proteolysis appears to be a slight modification of the accepted charge-relay mechanism for serine proteases.

Engineering substrate preference in subtilisin: structural and kinetic analysis of a specificity mutant.

Bacillus subtilisin has been a popular model protein for engineering altered substrate specificity. Although some studies have succeeded in increasing the specificity of subtilisin, they also demonstrate that high specificity is difficult to achieve solely by engineering selective substrate binding. In this paper, we analyze the structure and transient state kinetic behavior of Sbt160, a subtilisin engineered to strongly prefer substrates with phenylalanine or tyrosine at the P4 position. As in previous studies, we measure improvements in substrate affinity and overall specificity. Structural analysis of an inactive version of Sbt160 in complex with its cognate substrate reveals improved interactions at the S4 subsite with a P4 tyrosine. Comparison of transient state kinetic behavior against an optimal sequence (DFKAM) and a similar, but suboptimal, sequence (DVRAF) reveals the kinetic and thermodynamic basis for increased specificity, as well as the limitations of this approach. While highly selective substrate binding is achieved in Sbt160, several factors cause sequence specificity to fall short of that observed with natural processing subtilisins. First, for substrate sequences which are nearly optimal, the acylation reaction becomes faster than substrate dissociation. As a result, the level of discrimination among these substrates diminishes due to the coupling between substrate binding and the first chemical step (acylation). Second, although Sbt160 has 24-fold higher substrate affinity for the optimal substrate DFKAM than for DVRAF, the increased substrate binding energy is not translated into improved transition state stabilization of the acylation reaction. Finally, as interactions at subsites become stronger, the rate-determining step in peptide hydrolysis changes from acylation to product release. Thus, the release of the product becomes sluggish and leads to a low k(cat) for the reaction. This also leads to strong product inhibition of substrate turnover as the reaction progresses. The structural and kinetic analysis reveals that differences in the binding modes at subsites for substrates, transition states, and products are subtle and difficult to manipulate via straightforward protein engineering. These findings suggest several new strategies for engineering highly sequence selective enzymes.

Structural metamorphism and polymorphism in proteins on the brink of thermodynamic stability.

The classical view of the structure-function paradigm advanced by Anfinsen in the 1960s is that a protein's function is inextricably linked to its three-dimensional structure and is encrypted in its amino acid sequence. However, it is now known that a significant fraction of the proteome consists of intrinsically disordered proteins (IDPs). These proteins populate a polymorphic ensemble of conformations rather than a unique structure but are still capable of performing biological functions. At the boundary, between well-ordered and inherently disordered states are proteins that are on the brink of stability, either weakly stable ordered systems or disordered but on the verge of being stable. In such marginal states, even relatively minor changes can significantly alter the energy landscape, leading to large-scale conformational remodeling. Some proteins on the edge of stability are metamorphic, with the capacity to switch from one fold topology to another in response to an environmental trigger (e.g., pH, temperature/salt, redox). Many IDPs, on the other hand, are marginally unstable such that small perturbations (e.g., phosphorylation, ligands) tip the balance over to a range of ordered, partially ordered, or even more disordered states. In general, the structural transitions described by metamorphic fold switches and polymorphic IDPs possess a number of common features including low or diminished stability, large-scale conformational changes, critical disordered regions, latent or attenuated binding sites, and expansion of function. We suggest that these transitions are, therefore, conceptually and mechanistically analogous, representing adjacent regions in the continuum of order/disorder transitions.

An artificially evolved albumin binding module facilitates chemical shift epitope mapping of GA domain interactions with phylogenetically diverse albumins.

Protein G-related albumin-binding (GA) modules occur on the surface of numerous Gram-positive bacterial pathogens and their presence may promote bacterial growth and virulence in mammalian hosts. We recently used phage display selection to evolve a GA domain, PSD-1 (phage selected domain-1), which tightly bound phylogenetically diverse albumins. With respect to PSD-1's broad albumin binding specificity, it remained unclear how the evolved binding epitope compared to those of naturally occurring GA domains and whether PSD-1's binding mode was the same for different albumins. We investigate these questions here using chemical shift perturbation measurements of PSD-1 with rabbit serum albumin (RSA) and human serum albumin (HSA) and put the results in the context of previous work on structure and dynamics of GA domains. Combined, these data provide insights into the requirements for broad binding specificity in GA-albumin interactions. Moreover, we note that using the phage-optimized PSD-1 protein significantly diminishes the effects of exchange broadening at the binding interface between GA modules and albumin, presumably through stabilization of a ligand-bound conformation. The employment of artificially evolved domains may be generally useful in NMR structural studies of other protein-protein complexes.

Offset recombinant PCR: a simple but effective method for shuffling compact heterologous domains.

DNA shuffling and other in vitro recombination strategies have proven highly effective at generating complex libraries for mutagenesis studies. While most recombination techniques employ DNA polymerases in part of a multi-step process, few seek to exploit the natural recombinogenic tendencies and exponential amplification rates of PCR. Here, we characterize a simple but effective method for using standard PCR to promote high recombination frequencies among compact heterologous domains by locating the domains near one end of the template. In a typical amplification reaction, Pfu polymerase generated chimeric crossover events in 13% of the population when markers were separated by only 70 nt. The fraction of recombinant sequences reached 42% after six consecutive rounds of PCR, a value close to the 50% expected from a fully shuffled population. When homology within the recombinant region was reduced to 82%, the recombination frequency dropped by nearly half for a single amplification reaction and crossover events were clustered toward one end of the domain. Surprisingly, recombination frequencies for template populations with high and low sequence homologies converged after just four rounds of PCR, suggesting that the exponential accumulation of chimeric molecules in the PCR mixture serves to promote recombination within heterologous domains.

G148-GA3: a streptococcal virulence module with atypical thermodynamics of folding optimally binds human serum albumin at physiological temperatures.

The third albumin binding domain of streptococcal protein G strain 148 (G148-GA3) belongs to a novel class of prokaryotic albumin binding modules that is thought to support virulence in several bacterial species. Here, we characterize G148-GA3 folding and albumin binding by using differential scanning calorimetry and isothermal titration calorimetry to obtain the most complete set of thermodynamic state functions for any member of this medically significant module. When buffered at pH 7.0 the 46-amino acid alpha-helical domain melts at 72 degrees C and exhibits marginal stability (15 kJ/mol) at 37 degrees C. G148-GA3 unfolding is characterized by small contributions to entropy from non-hydrophobic forces and a low DeltaCp (1.1 kJ/(deg mol)). Isothermal titration calorimetry reveals that the domain has evolved to optimally bind human serum albumin near 37 degrees C with a binding constant of 1.4 x 10 7 M(-1). Analysis of G148-GA3 thermodynamics suggests that the domain experiences atypically small per residue changes in structural dynamics and heat capacity while transiting between folded and unfolded states.

Solution structure of the pro-hormone convertase 1 pro-domain from Mus musculus.

The solution structure of the mouse pro-hormone convertase (PC) 1 pro-domain was determined using heteronuclear NMR spectroscopy and is the first structure to be obtained for any of the domains in the convertase family. The ensemble of NMR-derived structures shows a well-ordered core consisting of a four-stranded antiparallel beta-sheet with two alpha-helices packed against one side of this sheet. Sequence homology suggests that the other eukaryotic PC pro-domains will have the same overall fold and most of the residues forming the hydrophobic core of PC1 are highly conserved within the PC family. However, some of the core residues are predicted by homology to be replaced by polar amino acid residues in other PC pro-domains and this may help to explain their marginal stability. Interestingly, the folding topology observed here is also seen for the pro-domain of bacterial subtilisin despite little or no sequence homology. Both the prokaryotic and eukaryotic structures have hydrophobic residues clustered on the solvent-accessible surface of their beta-sheets although the individual residue types differ. In the bacterial case this region is buried at the binding interface with the catalytic domain and, in the eukaryotic PC family, these surface residues are conserved. We therefore propose that the hydrophobic patch in the PC1 pro-domain is involved in the binding interface with its cognate catalytic domain in a similar manner to that seen for the bacterial system. The PC1 pro-domain structure also reveals potential mechanisms for the acid-induced dissociation of the complex between pro- and catalytic domains.

Mutational tipping points for switching protein folds and functions.

While disordered to ordered rearrangements are relatively common, the ability of proteins to switch from one ordered fold to a completely different fold is generally regarded as rare, and few fold switches have been characterized. Here, in a designed system, we examine the mutational requirements for transitioning between folds and functions. We show that switching between monomeric 3α and 4ß+α folds can occur in multiple ways with successive single amino acid changes at diverse residue positions, raising the likelihood that such transitions occur in the evolution of new folds. Even mutations on the periphery of the core can tip the balance between alternatively folded states. Ligand-binding studies illustrate that a new immunoglobulin G-binding function can be gained well before the relevant 4ß+α fold is appreciably populated in the unbound protein. The results provide new insights into the evolution of fold and function.

Proteins that switch folds.

An increasing number of proteins demonstrate the ability to switch between very different fold topologies, expanding their functional utility through new binding interactions. Recent examples of fold switching from naturally occurring and designed systems have a number of common features: (i) The structural transitions require states with diminished stability; (ii) Switching involves flexible regions in one conformer or the other; (iii) A new binding surface is revealed in the alternate fold that can lead to both stabilization of the alternative state and expansion of biological function. Fold switching not only provides insight into how new folds evolve, but also indicates that an amino acid sequence has more information content than previously thought. A polypeptide chain can encode a stable fold while simultaneously hiding latent propensities for alternative states with novel functions.

De novo structure generation using chemical shifts for proteins with high-sequence identity but different folds.

Proteins with high-sequence identity but very different folds present a special challenge to sequence-based protein structure prediction methods. In particular, a 56-residue three-helical bundle protein (GA(95)) and an alpha/beta-fold protein (GB(95)), which share 95% sequence identity, were targets in the CASP-8 structure prediction contest. With only 12 out of 300 submitted server-CASP8 models for GA(95) exhibiting the correct fold, this protein proved particularly challenging despite its small size. Here, we demonstrate that the information contained in NMR chemical shifts can readily be exploited by the CS-Rosetta structure prediction program and yields adequate convergence, even when input chemical shifts are limited to just amide (1)H(N) and (15)N or (1)H(N) and (1)H(alpha) values.