Exploring the conformational properties of the sequence space between two proteins with different folds: an experimental study

A comparison of the folding characteristics of free and ribosome-tethered polypeptide chains using limited proteolysis and mass spectrometry

The kinetics and thermodynamics of protein folding are commonly studied in vitro by denaturing/renaturing intact protein sequences. How these folding mechanisms relate to de novo folding that occurs as the nascent polypeptide emerges from the ribosome is much less well understood. Here, we have employed limited proteolysis followed by mass spectrometry analyses to compare directly free and ribosome-tethered polypeptide chains of the Src-homology 3 (SH3) domain and its unfolded variant, SH3-m10. The disordered variant was found to undergo faster proteolysis than SH3. Furthermore, the trypsin cleavage patterns observed show minor, but significant, differences for the free and ribosome-bound nascent chains, with significantly fewer tryptic peptides detected in the presence of ribosome. The results highlight the utility of limited proteolysis coupled with mass spectrometry for the structural analysis of these complex systems, and pave the way for detailed future analyses by combining this technique with chemical labeling methods (for example, hydrogen-deuterium exchange, photochemical oxidation) to analyze protein folding in real time, including in the presence of additional ribosome-associated factors.

Studying protein fold evolution with hybrids of differently folded homologs

To study the sequence determinants governing protein fold evolution, we generated hybrid sequences from two homologous proteins with 40% identity but different folds: Pfl 6 Cro, which has a mixed α + β structure, and Xfaso 1 Cro, which has an all α-helical structure. First, we first examined eight chimeric hybrids in which the more structurally conserved N-terminal half of one protein was fused to the more structurally divergent C-terminal half of the other. None of these chimeras folded, as judged by circular dichroism spectra and thermal melts, suggesting that both halves have strong intrinsic preferences for the native global fold pattern, and/or that the interfaces between the halves are not readily interchangeable. Second, we examined 10 hybrids in which blocks of the structurally divergent C-terminal region were exchanged. These hybrids showed varying levels of thermal stability and suggested that the key residues in the Xfaso 1 C terminus specifying the all-α fold were concentrated near the end of helix 4 in Xfaso 1, which aligns to the end of strand 2 in Pfl 6. Finally, we generated hybrid substitutions for each individual residue in this critical region and measured thermal stabilities. The results suggested that R47 and V48 were the strongest factors that excluded formation of the α + β fold in the C-terminal region of Xfaso 1. In support of this idea, we found that the folding stability of one of the original eight chimeras could be rescued by back-substituting these two residues. Overall, the results show not only that the key factors for Cro fold specificity and evolution are global and multifarious, but also that some all-α Cro proteins have a C-terminal subdomain sequence within a few substitutions of switching to the α + β fold.

Effect of single-site mutations on hydrophobic-polar lattice proteins

We developed a heuristic method for determining the ground-state degeneracy of hydrophobic-polar (HP) lattice proteins, based on Wang-Landau and multicanonical sampling. It is applied during comprehensive studies of single-site mutations in specific HP proteins with different sequences. The effects in which we are interested include structural changes in ground states, changes of ground-state energy, degeneracy, and thermodynamic properties of the system. With respect to mutations, both extremely sensitive and insensitive positions in the HP sequence have been found. That is, ground-state energies and degeneracies, as well as other thermodynamic and structural quantities, may be either largely unaffected or may change significantly due to mutation.

Protein folding on the ribosome studied using NMR spectroscopy

NMR spectroscopy is a powerful tool for the investigation of protein folding and misfolding, providing a characterization of molecular structure, dynamics and exchange processes, across a very wide range of timescales and with near atomic resolution. In recent years NMR methods have also been developed to study protein folding as it might occur within the cell, in a de novo manner, by observing the folding of nascent polypeptides in the process of emerging from the ribosome during synthesis. Despite the 2.3 MDa molecular weight of the bacterial 70S ribosome, many nascent polypeptides, and some ribosomal proteins, have sufficient local flexibility that sharp resonances may be observed in solution-state NMR spectra. In providing information on dynamic regions of the structure, NMR spectroscopy is therefore highly complementary to alternative methods such as X-ray crystallography and cryo-electron microscopy, which have successfully characterized the rigid core of the ribosome particle. However, the low working concentrations and limited sample stability associated with ribosome-nascent chain complexes means that such studies still present significant technical challenges to the NMR spectroscopist. This review will discuss the progress that has been made in this area, surveying all NMR studies that have been published to date, and with a particular focus on strategies for improving experimental sensitivity.

Concerted action of the ribosome and the associated chaperone trigger factor confines nascent polypeptide folding

How nascent polypeptides emerging from ribosomes fold into functional structures is poorly understood. Here, we monitor disulfide bond formation, protease resistance, and enzymatic activity in nascent polypeptides to show that in close proximity to the ribosome, conformational space and kinetics of folding are restricted. Folding constraints decrease incrementally with distance from the ribosome surface. Upon ribosome binding, the chaperone Trigger Factor counters folding also of longer nascent chains, to extents varying between different chain segments. Trigger Factor even binds and unfolds pre-existing folded structures, the unfolding activity being limited by the thermodynamic stability of nascent chains. Folding retardation and unfolding activities are not shared by the DnaK chaperone assisting later folding steps. These ribosome- and Trigger Factor-specific activities together constitute an efficient mechanism to prevent or even revert premature folding, effectively limiting misfolded intermediates during protein synthesis.

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.

Cotranslational structure acquisition of nascent polypeptides monitored by NMR spectroscopy

The folding of proteins in living cells may start during their synthesis when the polypeptides emerge gradually at the ribosomal exit tunnel. However, our current understanding of cotranslational folding processes at the atomic level is limited. We employed NMR spectroscopy to monitor the conformation of the SH3 domain from alpha-spectrin at sequential stages of elongation via in vivo ribosome-arrested (15)N,(13)C-labeled nascent polypeptides. These nascent chains exposed either the entire SH3 domain or C-terminally truncated segments thereof, thus providing snapshots of the translation process. We show that nascent SH3 polypeptides remain unstructured during elongation but fold into a compact, native-like beta-sheet assembly when the entire sequence information is available. Moreover, the ribosome neither imposes major conformational constraints nor significantly interacts with exposed unfolded nascent SH3 domain moieties. Our data provide evidence for a domainwise folding of the SH3 domain on ribosomes without significant population of folding intermediates. The domain follows a thermodynamically favorable pathway in which sequential folding units are stabilized, thus avoiding kinetic traps during the process of cotranslational folding.

The antibiotic ADEP reprogrammes ClpP, switching it from a regulated to an uncontrolled protease

A novel class of antibiotic acyldepsipeptides (designated ADEPs) exerts its unique antibacterial activity by targeting the peptidase caseinolytic protease P (ClpP). ClpP forms proteolytic complexes with heat shock proteins (Hsp100) that select and process substrate proteins for ClpP-mediated degradation. Here, we analyse the molecular mechanism of ADEP action and demonstrate that ADEPs abrogate ClpP interaction with cooperating Hsp100 adenosine triphosphatases (ATPases). Consequently, ADEP treated bacteria are affected in ClpP-dependent general and regulatory proteolysis. At the same time, ADEPs also activate ClpP by converting it from a tightly regulated peptidase, which can only degrade short peptides, into a proteolytic machinery that recognizes and degrades unfolded polypeptides. In vivo nascent polypeptide chains represent the putative primary target of ADEP-activated ClpP, providing a rationale for the antibacterial activity of the ADEPs. Thus, ADEPs cause a complete functional reprogramming of the Clp–protease complex.

Insights into protein aggregation by NMR characterization of insoluble SH3 mutants solubilized in salt-free water

Protein aggregation in vivo has been extensively associated with a large spectrum of human diseases. On the other hand, mechanistic insights into protein aggregation in vitro were incomplete due to the inability in solubilizing insoluble proteins for high-resolution biophysical investigations. However, a new avenue may be opened up by our recent discovery that previously-thought insoluble proteins can in fact be solubilized in salt-free water. Here we use this approach to study the NMR structural and dynamic properties of an insoluble SH3 mutant with a naturally-occurring insertion of Val22 at the tip of the diverging turn. The obtained results reveal: 1) regardless of whether the residue is Val, Ala, Asp or Arg, the insertion will render the first hNck2 SH3 domain to be insoluble in buffers. Nevertheless, all four mutants could be solubilized in salt-free water and appear to be largely unfolded as evident from their CD and NMR HSQC spectra. 2) Comparison of the chemical shift deviations reveals that while in V22-SH3 the second helical region is similarly populated as in the wild-type SH3 at pH 2.0, the first helical region is largely unformed. 3) In V22-SH3, many non-native medium-range NOEs manifest to define non-native helical conformations. In the meanwhile a small group of native-like long-range NOEs still persists, indicating the existence of a rudimentary native-like tertiary topology. 4) Although overall, V22-SH3 has significantly increased backbone motions on the ps-ns time scale, some regions still own restricted backbone motions as revealed by analyzing ¹⁵N relaxation data. Our study not only leads to the establishment of the first high-resolution structural and dynamic picture for an insoluble protein, but also shed more light on the molecular events for the nonhierarchical folding mechanism. Furthermore, a general mechanism is also proposed for in vivo protein aggregation triggered by the genetic mutation and posttranslational modification.

Is protein classification necessary? Toward alternative approaches to function annotation

The current nonredundant protein sequence database contains over seven million entries and the number of individual functional domains is significantly larger than this value. The vast quantity of data associated with these proteins poses enormous challenges to any attempt at function annotation. Classification of proteins into sequence and structural groups has been widely used as an approach to simplifying the problem. In this article we question such strategies. We describe how the multifunctionality and structural diversity of even closely related proteins confounds efforts to assign function on the basis of overall sequence or structural similarity. Rather, we suggest that strategies that avoid classification may offer a more robust approach to protein function annotation.

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.

Stylus: a system for evolutionary experimentation based on a protein/proteome model with non-arbitrary functional constraints

The study of protein evolution is complicated by the vast size of protein sequence space, the huge number of possible protein folds, and the extraordinary complexity of the causal relationships between protein sequence, structure, and function. Much simpler model constructs may therefore provide an attractive complement to experimental studies in this area. Lattice models, which have long been useful in studies of protein folding, have found increasing use here. However, while these models incorporate actual sequences and structures (albeit non-biological ones), they incorporate no actual functions—relying instead on largely arbitrary structural criteria as a proxy for function. In view of the central importance of function to evolution, and the impossibility of incorporating real functional constraints without real function, it is important that protein-like models be developed around real structure–function relationships. Here we describe such a model and introduce open-source software that implements it. The model is based on the structure–function relationship in written language, where structures are two-dimensional ink paths and functions are the meanings that result when these paths form legible characters. To capture something like the hierarchical complexity of protein structure, we use the traditional characters of Chinese origin. Twenty coplanar vectors, encoded by base triplets, act like amino acids in building the character forms. This vector-world model captures many aspects of real proteins, including life-size sequences, a life-size structural repertoire, a realistic genetic code, secondary, tertiary, and quaternary structure, structural domains and motifs, operon-like genetic structures, and layered functional complexity up to a level resembling bacterial genomes and proteomes. Stylus is a full-featured implementation of the vector world for Unix systems. To demonstrate the utility of Stylus, we generated a sample set of homologous vector proteins by evolving successive lines from a single starting gene. These homologues show sequence and structure divergence resembling those of natural homologues in many respects, suggesting that the system may be sufficiently life-like for informative comparison to biology.

Molecular mechanism and structure of Trigger Factor bound to the translating ribosome

Ribosome-associated chaperone Trigger Factor (TF) initiates folding of newly synthesized proteins in bacteria. Here, we pinpoint by site-specific crosslinking the sequence of molecular interactions of Escherichia coli TF and nascent chains during translation. Furthermore, we provide the first full-length structure of TF associated with ribosome-nascent chain complexes by using cryo-electron microscopy. In its active state, TF arches over the ribosomal exit tunnel accepting nascent chains in a protective void. The growing nascent chain initially follows a predefined path through the entire interior of TF in an unfolded conformation, and even after folding into a domain it remains accommodated inside the protective cavity of ribosome-bound TF. The adaptability to accept nascent chains of different length and folding states may explain how TF is able to assist co-translational folding of all kinds of nascent polypeptides during ongoing synthesis. Moreover, we suggest a model of how TF's chaperoning function can be coordinated with the co-translational processing and membrane targeting of nascent polypeptides by other ribosome-associated factors.

Transitive homology-guided structural studies lead to discovery of Cro proteins with 40% sequence identity but different folds

Proteins that share common ancestry may differ in structure and function because of divergent evolution of their amino acid sequences. For a typical diverse protein superfamily, the properties of a few scattered members are known from experiment. A satisfying picture of functional and structural evolution in relation to sequence changes, however, may require characterization of a larger, well chosen subset. Here, we employ a "stepping-stone" method, based on transitive homology, to target sequences intermediate between two related proteins with known divergent properties. We apply the approach to the question of how new protein folds can evolve from preexisting folds and, in particular, to an evolutionary change in secondary structure and oligomeric state in the Cro family of bacteriophage transcription factors, initially identified by sequence-structure comparison of distant homologs from phages P22 and lambda. We report crystal structures of two Cro proteins, Xfaso 1 and Pfl 6, with sequences intermediate between those of P22 and lambda. The domains show 40% sequence identity but differ by switching of alpha-helix to beta-sheet in a C-terminal region spanning approximately 25 residues. Sedimentation analysis also suggests a correlation between helix-to-sheet conversion and strengthened dimerization.

Dynamics of trigger factor interaction with translating ribosomes

In all organisms ribosome-associated chaperones assist early steps of protein folding. To elucidate the mechanism of their action, we determined the kinetics of individual steps of the ribosome binding/release cycle of bacterial trigger factor (TF), using fluorescently labeled chaperone and ribosome-nascent chain complexes. Both the association and dissociation rates of TF-ribosome complexes are modulated by nascent chains, whereby their length, sequence, and folding status are influencing parameters. However, the effect of the folding status is modest, indicating that TF can bind small globular domains and accommodate them within its substrate binding cavity. In general, the presence of a nascent chain causes an up to 9-fold increase in the rate of TF association, which provides a kinetic explanation for the observed ability of TF to efficiently compete with other cytosolic chaperones for binding to nascent chains. Furthermore, a subset of longer nascent polypeptides promotes the stabilization of TF-ribosome complexes, which increases the half-life of these complexes from 15 to 50 s. Nascent chains thus regulate their folding environment generated by ribosome-associated chaperones.

The design and characterization of two proteins with 88% sequence identity but different structure and function

To identify a simplified code for conformational switching, we have redesigned two natural proteins to have 88% sequence identity but different tertiary structures: a 3-alpha helix fold and an alpha/beta fold. We describe the design of these homologous heteromorphic proteins, their structural properties as determined by NMR, their conformational stabilities, and their affinities for their respective ligands: IgG and serum albumin. Each of these proteins is completely folded at 25 degrees C, is monomeric, and retains the native binding activity. The complete binding epitope for both ligands is encoded within each of the proteins. The IgG-binding epitope is functional only in the alpha/beta fold, and the albumin-binding epitope is functional only in the 3-alpha fold. These results demonstrate that two monomeric folds and two different functions can be encoded with only 12% of the amino acids in a protein (7 of 56). The fact that 49 aa in these proteins are compatible with both folds shows that the essential information determining a fold can be highly concentrated in a few amino acids and that a very limited subset of interactions in the protein can tip the balance from one monomer fold to another. This delicate balance helps explain why protein structure prediction is so challenging. Furthermore, because a few mutations can result in both new conformation and new function, the evolution of new folds driven by natural selection for alternative functions may be much more probable than previously recognized.

Protein engineering: opportunities and challenges

The extraordinary properties of natural proteins demonstrate that life-like protein engineering is both achievable and valuable. Rapid progress and impressive results have been made towards this goal using rational design and random techniques or a combination of both. However, we still do not have a general theory on how to specify a structure that is suited to a target function nor can we specify a sequence that folds to a target structure. There is also overreliance on the Darwinian blind search to obtain practical results. In the long run, random methods cannot replace insight in constructing life-like proteins. For the near future, however, in enzyme development, we need to rely on a combination of both.

Solution NMR structures of IgG binding domains with artificially evolved high levels of sequence identity but different folds

We describe here the solution NMR structures of two IgG binding domains with highly homologous sequences but different three-dimensional structures. The proteins, G311 and A219, are derived from the IgG binding domains of their wild-type counterparts, protein G and protein A, respectively. Through a series of site-directed mutations and phage display selections, the sequences of G311 and A219 were designed to converge to a point of high-level sequence identity while keeping their respective wild-type tertiary folds. Structures of both artificially evolved sequences were determined by NMR spectroscopy. The main chain fold of G311 can be superimposed on the wild-type alpha/beta protein G structure with a backbone rmsd of 1.4 A, and the A219 structure can be overlaid on the wild-type three-alpha-helix protein A fold also with a backbone rmsd of 1.4 A. The structure of G311, in particular, accommodates a large number of mutational changes without undergoing a change in the overall fold of the main chain. The structural differences are maintained despite a high level (59%) of sequence identity. These proteins serve as starting points for further experiments that will probe basic concepts of protein folding and conformational switching.

Trigger factor forms a protective shield for nascent polypeptides at the ribosome

In prokaryotes, the ribosome-associated Trigger Factor is the first chaperone newly synthesized polypeptides encounter when they emerge from the ribosomal exit tunnel. The effects that Trigger Factor exerts on nascent polypeptides, however, remain unclear. Here we analyzed the potential of the Trigger Factor to shield nascent polypeptides at the ribosome. A set of arrested nascent polypeptides differing in origin, size, and folding status were synthesized in an Escherichia coli-based in vitro transcription/translation system and tested for sensitivity to degradation by the unspecific protease proteinase K. In the absence of Trigger Factor, nascent polypeptides exposed outside the ribosomal exit tunnel were rapidly degraded unless they were folded into a compact domain. The presence of Trigger Factor, as well as a Trigger Factor fragment lacking its peptidyl-prolyl isomerase domain, counteracted degradation of all unfolded nascent polypeptides tested. This protective function was specific for ribosome-tethered Trigger Factor, since neither non-ribosomal Trigger Factor nor the DnaK system, which cooperates with Trigger Factor in the folding process in vivo, revealed a comparable efficiency in protection. Furthermore, shielding by Trigger Factor was not restricted to short stretches of nascent chains but was evident for large, non-native nascent polypeptides exposing up to 41 kDa outside the ribosome. We suggest that Trigger Factor supports productive de novo folding by shielding nascent polypeptides on the ribosome thereby preventing untimely degradation or aggregation processes. This protected environment provided by Trigger Factor might be particularly important for large multidomain proteins to fold productively into their native states.

Rapid evolution in conformational space: a study of loop regions in a ubiquitous GTP binding domain

The rapidly evolving subsets of a protein are often evident in multiple sequence alignments as poorly defined, gap-containing regions. We investigated the 3D context of these regions observed in 28 protein structures containing a GTP-binding domain assumed to be homologous to the transforming factor p21-RAS. The phylogenetic depth of this data set is such that it is possible to observe lineages sharing a common protein core that diverged early in the eukaryotic cell history. The sequence variability among these homolog proteins is directly linked to the structural variability of surface loops. We demonstrate that these regions are self-contained and thus mostly free of the evolutionary constraints imposed by the conserved core of the domain. These intraloop interactions have the property to create stem-like structures. Interestingly, these stem-like structures can be observed in loops of varying size, up to the size of small protein domains. We propose a model under which the diversity of protein topologies observed in these loops can be the product of a stochastic sampling of sequence and conformational space in a near-neutral fashion, while the proximity of the functional features of the domain core allows novel beneficial traits to be fixed. Our comparative observations, limited here to the proteins containing the RAS-like GTP-binding domain, suggest that a stochastic process of insertion/deletion analogous to "budding" of loops is a likely mechanism of structural innovation. Such a framework could be experimentally exploited to investigate the folding of increasingly complex model inserts.

Estimating the Prevalence of Protein Sequences Adopting Functional Enzyme Folds

Proteins employ a wide variety of folds to perform their biological functions. How are these folds first acquired? An important step toward answering this is to obtain an estimate of the overall prevalence of sequences adopting functional folds. Since tertiary structure is needed for a typical enzyme active site to form, one way to obtain this estimate is to measure the prevalence of sequences supporting a working active site. Although the immense number of sequence combinations makes wholly random sampling unfeasible, two key simplifications may provide a solution. First, given the importance of hydrophobic interactions to protein folding, it seems likely that the sample space can be restricted to sequences carrying the hydropathic signature of a known fold. Second, because folds are stabilized by the cooperative action of many local interactions distributed throughout the structure, the overall problem of fold stabilization may be viewed reasonably as a collection of coupled local problems. This enables the difficulty of the whole problem to be assessed by assessing the difficulty of several smaller problems. Using these simplifications, the difficulty of specifying a working beta-lactamase domain is assessed here. An alignment of homologous domain sequences is used to deduce the pattern of hydropathic constraints along chains that form the domain fold. Starting with a weakly functional sequence carrying this signature, clusters of ten side-chains within the fold are replaced randomly, within the boundaries of the signature, and tested for function. The prevalence of low-level function in four such experiments indicates that roughly one in 10(64) signature-consistent sequences forms a working domain. Combined with the estimated prevalence of plausible hydropathic patterns (for any fold) and of relevant folds for particular functions, this implies the overall prevalence of sequences performing a specific function by any domain-sized fold may be as low as 1 in 10(77), adding to the body of evidence that functional folds require highly extraordinary sequences.

Searching for folded proteins in vitro and in silico

Understanding the sequence determinants of protein structure, stability and folding is critical for understanding how natural proteins have evolved and how proteins can be engineered to perform novel functions. The complexity of the protein folding problem requires the ability to search large volumes of sequence space for proteins with specific structural or functional characteristics. Here we describe our efforts to identify novel proteins using a phage-display selection strategy from a 'mini-exon' shuffling library generated from the yeast genome and from completely random sequence libraries, and compare the results to recent successes in generating novel proteins using in silico protein design.

An in silico exploration of the neutral network in protein sequence space

Designating amino-acid sequences that fold into a common main-chain structure as "neutral sequences" for the structure, regardless of their function or stability, we investigated the distribution of neutral sequences in protein sequence space. For four distinct target structures (alpha, beta,alpha/beta and alpha+beta types) with the same chain length of 108, we generated the respective neutral sequences by using the inverse folding technique with a knowledge-based potential function. We assumed that neutral sequences for a protein structure have Z scores higher than or equal to fixed thresholds, where thresholds are defined as the Z score for the corresponding native sequence (case 1) or much greater Z score (case 2). An exploring walk simulation suggested that the neutral sequences mapped into the sequence space were connected with each other through straight neutral paths and formed an inherent neutral network over the sequence space. Through another exploring walk simulation, we investigated contiguous regions between or among the neutral networks for the distinct protein structures and obtained the following results. The closest approach distance between the two neutral networks ranged from 5 to 29 on the Hamming distance scale, showing a linear increase against the threshold values. The sequences located at the "interchange" regions between the two neutral networks have intermediate sequence-profile-scores for both corresponding structures. Introducing a "ball" in the sequence space that contains at least one neutral sequence for each of the four structures, we found that the minimal radius of the ball that is centered at an arbitrary position ranged from 35 to 50, while the minimal radius of the ball that is centered at a certain special position ranged from 20 to 30, in the Hamming distance scale. The relatively small Hamming distances (5-30) may support an evolution mechanism by transferring from a network for a structure to another network for a more beneficial structure via the interchange regions.

Stabilization of human pancreatic ribonuclease through mutation at its N-terminal edge

Enzyme stability can be an important parameter in the design of recombinant toxins because unstable proteins are often degraded before they can reach their cellular target. There is great interest in the design of human pancreatic ribonuclease variants that could be cytotoxic against tumoral cells. To this end, some residues in the protein need to be substituted, but this may result in a loss of stability. Previous papers have reported the production of N- and C-terminal human pancreatic ribonuclease variants with increased thermal stability. Here, we investigated the contribution of the different amino acid changes at the N-terminus of the protein to its thermostability increase. We show that this increase correlates with the helical propensity of the first alpha-helix of the protein. On the other hand, deletion of the four last residues of the protein does not affect its thermal stability. These results set the basis for the design of a human pancreatic ribonuclease template on which amino acid substitutions can be made that could render the enzyme cytotoxic, without an important loss in its stability.

Fold change in evolution of protein structures

Typically, protein spatial structures are more conserved in evolution than amino acid sequences. However, the recent explosion of sequence and structure information accompanied by the development of powerful computational methods led to the accumulation of examples of homologous proteins with globally distinct structures. Significant sequence conservation, local structural resemblance, and functional similarity strongly indicate evolutionary relationships between these proteins despite pronounced structural differences at the fold level. Several mechanisms such as insertions/deletions/substitutions, circular permutations, and rearrangements in beta-sheet topologies account for the majority of detected structural irregularities. The existence of evolutionarily related proteins that possess different folds brings new challenges to the homology modeling techniques and the structure classification strategies and offers new opportunities for protein design in experimental studies.

Mapping protein sequence spaces by recurrence quantification analysis: a case study on chimeric structures

Recurrence quantification analysis (RQA) was used to characterize the folding properties of 22 chimeric sequences derived from two parent proteins of similar length but different three-dimensional arrangement. A non-linear relation between sequence data and their RQA representation was revealed, which points to new information carried by this method as compared with classical best-alignment methods. This new information is significantly correlated with the folding properties of the hybrid polypeptide chains, as substantiated by careful statistical analysis of the recurrence plots' numerical descriptors, thus encouraging their systematic use to complement sequence data in both proteomics and protein engineering tasks. Even the direct visual screening of the qualitative graphical features of recurrence plots is shown to provide useful hints to discriminate between different recurrence structures of protein sequences.

A tale of two secondary structure elements: when a beta-hairpin becomes an alpha-helix

In this work, we have analyzed the relative importance of secondary versus tertiary interactions in stabilizing and guiding protein folding. For this purpose, we have designed four different mutants to replace the alpha-helix of the GB1 domain by a sequence with strong beta-hairpin propensity in isolation. In particular, we have chosen the sequence of the second beta-hairpin of the GB1 domain, which populates the native conformation in aqueous solution to a significant extent. The resulting protein has roughly 30 % of its sequence duplicated and maintains the 3D-structure of the wild-type protein, but with lower stability (up to -5 kcal/mol). The loss of intrinsic helix stability accounts for about 80 % of the decrease in free energy, illustrating the importance of local interactions in protein stability. Interestingly enough, all the mutant proteins, included the one with the duplicated beta-hairpin sequence, fold with similar rates as the GB1 domain. Essentially, it is the nature of the rate-limiting step in the folding reaction that determines whether a particular interaction will speed up, or not, the folding rates. While local contacts are important in determining protein stability, residues involved in tertiary contacts in combination with the topology of the native fold, seem to be responsible for the specificity of protein structures. Proteins with non-native secondary structure tendencies can adopt stable folds and be as efficient in folding as those proteins with native-like propensities.