Designed protein reveals structural determinants of extreme kinetic stability

Integrated strategies for efficient production of Streptomyces mobaraensis transglutaminase in Komagataella phaffii

Transglutaminase (TGase) from Streptomyces mobaraensis commonly used to improve protein-based foods due to its unique enzymatic reactions, which imply considerable attention in its production. Recently, TGase exhibit broad market potential in non-food industries. However, achieving efficient synthesis of TGase remains a significant challenge. Herein, we achieved a substantial amount of a fully functional and kinetically stable TGase produced by Komagataella phaffii (Pichia pastoris) using multiple strategies including Geneticin (G418) screening, combinatorial mutations, promoter optimization, and co-expression. The active TGase expression reached a maximum of 10.1 U mL-1 in shake flask upon 96 h of induction, which was 3.8-fold of the wild type. Also, the engineered strain exhibited a 6.4-fold increase in half-life and a 2-fold increase in specific activity, reaching 172.67 min at 60 °C (t1/2(60 °C)) and 65.3 U mg-1, respectively. Moreover, the high-cell density cultivation in 5-L fermenter was also applied to test the productivity at large scale. Following optimization at a fermenter, the secretory yield of TGase reached 47.96 U mL-1 in the culture supernatant. Given the complexity inherent in protein expression and secretion, our research is of great significance and offers a comprehensive guide for improving the production of a wide range of heterologous proteins.

Engineering the kinetic stability of a β-trefoil protein by tuning its topological complexity

Kinetic stability, defined as the rate of protein unfolding, is central to determining the functional lifetime of proteins, both in nature and in wide-ranging medical and biotechnological applications. Further, high kinetic stability is generally correlated with high resistance against chemical and thermal denaturation, as well as proteolytic degradation. Despite its significance, specific mechanisms governing kinetic stability remain largely unknown, and few studies address the rational design of kinetic stability. Here, we describe a method for designing protein kinetic stability that uses protein long-range order, absolute contact order, and simulated free energy barriers of unfolding to quantitatively analyze and predict unfolding kinetics. We analyze two β-trefoil proteins: hisactophilin, a quasi-three-fold symmetric natural protein with moderate stability, and ThreeFoil, a designed three-fold symmetric protein with extremely high kinetic stability. The quantitative analysis identifies marked differences in long-range interactions across the protein hydrophobic cores that partially account for the differences in kinetic stability. Swapping the core interactions of ThreeFoil into hisactophilin increases kinetic stability with close agreement between predicted and experimentally measured unfolding rates. These results demonstrate the predictive power of readily applied measures of protein topology for altering kinetic stability and recommend core engineering as a tractable target for rationally designing kinetic stability that may be widely applicable.

Potential Self-Peptide Inhibitors of the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) plays an essential role in viral replication, cleaving viral polyproteins into functional proteins. This makes Mpro an important drug target. Mpro consists of an N-terminal catalytic domain and a C-terminal α-helical domain (MproC). Previous studies have shown that peptides derived from a given protein sequence (self-peptides) can affect the folding and, in turn, the function of that protein. Since the SARS-CoV-1 MproC is known to stabilize its Mpro and regulate its function, we hypothesized that SARS-CoV-2 MproC-derived self-peptides may modulate the folding and the function of SARS-CoV-2 Mpro. To test this, we studied the folding of MproC in the presence of various self-peptides using coarse-grained structure-based models and molecular dynamics simulations. In these simulations of MproC and one self-peptide, we found that two self-peptides, the α1-helix and the loop between α4 and α5 (loop4), could replace the equivalent native sequences in the MproC structure. Replacement of either sequence in full-length Mpro should, in principle, be able to perturb Mpro function albeit through different mechanisms. Some general principles for the rational design of self-peptide inhibitors emerge: The simulations show that prefolded self-peptides are more likely to replace native sequences than those which do not possess structure. Additionally, the α1-helix self-peptide is kinetically stable and once inserted rarely exchanges with the native α1-helix, while the loop4 self-peptide is easily replaced by the native loop4, making it less useful for modulating function. In summary, a prefolded α1-derived peptide should be able to inhibit SARS-CoV-2 Mpro function.

Effects of Surface Tethering on the Thermodynamics and Kinetics of Frustrated Protein Folding

The interaction between the protein and surface plays an important role in biology and biotechnology. To understand how surface tethering influences the folding behavior of frustrated proteins, in this work, we systematically study the thermodynamics and folding kinetics of the bacterial immunity protein Im7 and Fyn SH3 domain tethered to a surface using Langevin dynamics simulations. Upon surface tethering, the stabilization often results from the entropic effect, whereas the destabilization is usually caused by either an energetic or entropic effect. For the Fyn SH3 domain with a two-state folding manner, the influence of nonnative interactions on thermodynamic stability is not significant, while nonnative interactions can weaken the effect of surface tethering on the change in the folding rate. By contrast, for the frustrated protein Im7, depending on where the protein is tethered, the surface tethering can promote or suppress misfolding by modulating specific nonnative contacts, thereby altering the folding rate and folding mechanism. Because surface tethering can change the intrachain diffusivity of unfolding, the kinetic stability cannot be well captured by the thermodynamic stability at some tether points. This study should be helpful in general to understand how surface tethering affects the folding energy landscape of frustrated proteins.

Spontaneous Orthogonal Protein Crosslinking via a Genetically Encoded 2-Carboxy-4-Aryl-1,2,3-Triazole

Here we report the design of N2 -carboxy-4-aryl-1,2,3-triazole-lysines (CATKs) and their site-specific incorporation into proteins via genetic code expansion. When introduced into the protein dimer interface, CATKs permitted spontaneous, proximity-driven, site-selective crosslinking to generate covalent protein dimers in living cells, with phenyl-bearing CATK-1 exhibiting high reactivity toward the proximal Lys and Tyr. Furthermore, when introduced into the N-terminal β-strand of either a single-chain VHH antibody or a supercharged monobody, CATK-1 enabled site-specific, inter-strand, orthogonal crosslinking with a proximal Tyr located on the opposing β-strand. Compared with a non-crosslinked monobody, the orthogonally crosslinked monobody displayed improved cellular uptake and enhanced proteolytic stability against an endosomal enzyme. The robust crosslinking reactivity of CATKs should facilitate the design of novel protein topologies with improved physicochemical properties.

Variable and Conserved Regions of Secondary Structure in the β-Trefoil Fold: Structure Versus Function

β-trefoil proteins exhibit an approximate C₃ rotational symmetry. An analysis of the secondary structure for members of this diverse superfamily of proteins indicates that it is comprised of remarkably conserved β-strands and highly-divergent turn regions. A fundamental “minimal” architecture can be identified that is devoid of heterogenous and extended turn regions, and is conserved among all family members. Conversely, the different functional families of β-trefoils can potentially be identified by their unique turn patterns (or turn “signature”). Such analyses provide clues as to the evolution of the β-trefoil family, suggesting a folding/stability role for the β-strands and a functional role for turn regions. This viewpoint can also guide de novo protein design of β-trefoil proteins having novel functionality.

Proteome-wide landscape of solubility limits in a bacterial cell

Proteins are prone to aggregate when expressed above their solubility limits. Aggregation may occur rapidly, potentially as early as proteins emerge from the ribosome, or slowly, following synthesis. However, in vivo data on aggregation rates are scarce. Here, we classified the Escherichia coli proteome into rapidly and slowly aggregating proteins using an in vivo image-based screen coupled with machine learning. We find that the majority (70%) of cytosolic proteins that become insoluble upon overexpression have relatively low rates of aggregation and are unlikely to aggregate co-translationally. Remarkably, such proteins exhibit higher folding rates compared to rapidly aggregating proteins, potentially implying that they aggregate after reaching their folded states. Furthermore, we find that a substantial fraction (~ 35%) of the proteome remain soluble at concentrations much higher than those found naturally, indicating a large margin of safety to tolerate gene expression changes. We show that high disorder content and low surface stickiness are major determinants of high solubility and are favored in abundant bacterial proteins. Overall, our study provides a global view of aggregation rates and hence solubility limits of proteins in a bacterial cell.

Combining Ancestral Reconstruction with Folding-Landscape Simulations to Engineer Heterologous Protein Expression

Obligate symbionts typically exhibit high evolutionary rates. Consequently, their proteins may differ considerably from their modern and ancestral homologs in terms of both sequence and properties, thus providing excellent models to study protein evolution. Also, obligate symbionts are challenging to culture in the lab and proteins from uncultured organisms must be produced in heterologous hosts using recombinant DNA technology. Obligate symbionts thus replicate a fundamental scenario of metagenomics studies aimed at the functional characterization and biotechnological exploitation of proteins from the bacteria in soil. Here, we use the thioredoxin from Candidatus Photodesmus katoptron, an uncultured symbiont of flashlight fish, to explore evolutionary and engineering aspects of protein folding in heterologous hosts. The symbiont protein is a standard thioredoxin in terms of 3D-structure, stability and redox activity. However, its folding outside the original host is severely impaired, as shown by a very slow refolding in vitro and an inefficient expression in E. coli that leads mostly to insoluble protein. By contrast, resurrected Precambrian thioredoxins express efficiently in E. coli, plausibly reflecting an ancient adaptation to unassisted folding. We have used a statistical-mechanical model of the folding landscape to guide back-to-ancestor engineering of the symbiont protein. Remarkably, we find that the efficiency of heterologous expression correlates with the in vitro (i.e., unassisted) folding rate and that the ancestral expression efficiency can be achieved with only 1-2 back-to-ancestor replacements. These results demonstrate a minimal-perturbation, sequence-engineering approach to rescue inefficient heterologous expression which may potentially be useful in metagenomics efforts targeting recent adaptations.

Understanding the Folding Mediated Assembly of the Bacteriophage MS2 Coat Protein Dimers

The capsids of RNA viruses such as MS2 are great models for studying protein self-assembly because they are made almost entirely of multiple copies of a single coat protein (CP). Although CP is the minimal repeating unit of the capsid, previous studies have shown that CP exists as a homodimer (CP2) even in an acid-disassembled system, indicating that CP2 is an obligate dimer. Here, we investigate the molecular basis of this obligate dimerization using coarse-grained structure-based models and molecular dynamics simulations. We find that, unlike monomeric proteins of similar size, CP populates a single partially folded ensemble whose "foldedness" is sensitive to denaturing conditions. In contrast, CP2 folds similarly to single-domain proteins populating only the folded and the unfolded ensembles, separated by a prominent folding free energy barrier. Several intramonomer contacts form early, but the CP2 folding barrier is crossed only when the intermonomer contacts are made. A dissection of the structure of CP2 through mutant folding simulations shows that the folding barrier arises both from the topology of CP and the interface contacts of CP2. Together, our results show that CP2 is an obligate dimer because of kinetic stability, that is, dimerization induces a folding barrier and that makes it difficult for proteins in the dimer minimum to partially unfold and access the monomeric state without completely unfolding. We discuss the advantages of this obligate dimerization in the context of dimer design and virus stability.

Cooperative hydrophobic core interactions in the β-trefoil architecture

Symmetric protein architectures have a compelling aesthetic that suggests a plausible evolutionary process (i.e., gene duplication/fusion) yielding complex architecture from a simpler structural motif. Furthermore, symmetry inspires a practical approach to computational protein design that substantially reduces the combinatorial explosion problem, and may provide practical solutions for structure optimization. Despite such broad relevance, the role of structural symmetry in the key area of hydrophobic core-packing cooperativity has not been adequately studied. In the present report, the threefold rotational symmetry intrinsic to the β-trefoil architecture is shown to form a geometric basis for highly-cooperative core-packing interactions that both stabilize the local repeating motif and promote oligomerization/long-range contacts in the folding process. Symmetry in the β-trefoil structure also permits tolerance towards mutational drift that involves a structural quasi-equivalence at several key core positions.

Conserved buried water molecules enable the β-trefoil architecture

Available high-resolution crystal structures for the family of β-trefoil proteins in the structural databank were queried for buried waters. Such waters were classified as either: (a) unique to a particular domain, family, or superfamily or (b) conserved among all β-trefoil folds. Three buried waters conserved among all β-trefoil folds were identified. These waters are related by the threefold rotational pseudosymmetry characteristic of this protein architecture (representing three instances of an identical structural environment within each repeating trefoil-fold motif). The structural properties of this buried water are remarkable and include: residing in a cavity space no larger than a single water molecule, exhibiting a positional uncertainty (i.e., normalized B-factor) substantially lower than the average Cα atom, providing essentially ideal H-bonding geometry with three solvent-inaccessible main chain groups, simultaneously serving as a bridging H-bond for three different β-strands at a point of secondary structure divergence, and orienting conserved hydrophobic side chains to form a nascent core-packing group. Other published work supports an interpretation that these interactions are key to the formation of an efficient folding nucleus and folded thermostability. The fundamental threefold symmetric structural element of the β-trefoil fold is therefore, surprisingly, a buried water molecule.

Molecular Switch between Structural Compaction and Thermodynamic Stability by the Xxx-Pro Interface in Transmembrane β-Barrels

Transmembrane β-barrel scaffolds found in outer membrane proteins are formed and stabilized by a defined pattern of interstrand intraprotein H-bonds, in hydrophobic lipid bilayers. Introducing the conformationally constrained proline in β-barrels can cause significant destabilization of these structural regions that require H-bonding, with proline additionally acting as a secondary structure breaker. Membrane protein β-barrels are therefore expected to show poor tolerance to the presence of a transmembrane proline. Here, we assign a thermodynamic measure for the extent to which a single proline can be tolerated at the C-terminal interface of the model transmembrane β-barrel PagP. We find that proline drastically destabilizes PagP by 7.0 kcal mol^-1 with respect to wild-type PagP (F¹⁶¹ → P¹⁶¹). Interestingly, strategic modulation of the preceding residue can modify the measured energetics. Placing a hydrophobic or bulky side chain as the preceding residue increases the thermodynamic stability by ≤8.0 kcal mol^-1. While polar substituents at the preceding residue decrease the PagP stability, these residues demonstrate stronger tertiary packing interactions in the barrel and retain the catalytic activity of native PagP. This biophysical interplay between enhanced thermodynamic stability and attaining a structurally compact functional β-barrel scaffold highlights the detrimental effect caused by proline incorporation. Our findings also reveal alternative mechanisms that protein sequences can employ to salvage the structural integrity of transmembrane protein structures.

Structure and engineering of tandem repeat lectins

Through their ability to bind complex glycoconjugates, lectins have unique specificity and potential for biomedical and biotechnological applications. In particular, lectins with short repeated peptides forming carbohydrate-binding domains are not only of high interest for understanding protein evolution but can also be used as scaffold for engineering novel receptors. Synthetic glycobiology now provides the tools for engineering the specificity of lectins as well as their structure, multivalency and topologies. This review focuses on the structure and diversity of two families of tandem-repeat lectins, that is, β-trefoils and β-propellers, demonstrated as the most promising scaffold for engineering novel lectins.

An extensively glycosylated archaeal pilus survives extreme conditions

Pili on the surface of Sulfolobus islandicus are used for many functions, and serve as receptors for certain archaeal viruses. The cells grow optimally at pH 3 and ~80 °C, exposing these extracellular appendages to a very harsh environment. The pili, when removed from cells, resist digestion by trypsin or pepsin, and survive boiling in sodium dodecyl sulfate or 5 M guanidine hydrochloride. We used electron cryo-microscopy to determine the structure of these filaments at 4.1 Å resolution. An atomic model was built by combining the electron density map with bioinformatics without previous knowledge of the pilin sequence-an approach that should prove useful for assemblies where all of the components are not known. The atomic structure of the pilus was unusual, with almost one-third of the residues being either threonine or serine, and with many hydrophobic surface residues. While the map showed extra density consistent with glycosylation for only three residues, mass measurements suggested extensive glycosylation. We propose that this extensive glycosylation renders these filaments soluble and provides the remarkable structural stability. We also show that the overall fold of the archaeal pilin is remarkably similar to that of archaeal flagellin, establishing common evolutionary origins.

Sarkosyl: A milder detergent than SDS for identifying proteins with moderately high hyperstability using gel electrophoresis

Sodium dodecyl sulfate (SDS) is a detergent used as a strong denaturant of proteins in gel electrophoresis. It has previously been shown that certain hyperstable, also known as kinetically stable, proteins are resistant to SDS and thus require heating for their denaturation in the presence of SDS. Because of its high denaturing strength, relatively few proteins are resistant to SDS thereby limiting the current use of SDS-PAGE for identifying hyperstable degradation-resistant proteins. In this study, we show that sarkosyl, a milder detergent than SDS, is able to identify proteins with moderately high kinetic stability that lack SDS-resistance. Our assay involves running and subsequently comparing boiled and unheated protein samples containing sarkosyl, instead of SDS, on PAGE gels and identifying subsequent differences in protein migration. Our results also show that sarkosyl and SDS may be combined in PAGE experiments at varying relative percentages to obtain semi-quantitative information about a protein's kinetic stability in a range inaccessible by probing through native- or SDS-PAGE. Using protein extracts from various legumes as model systems, we detected proteins with a range of protein stability from nearly SDS-resistant to barely sarkosyl resistant.

Computational design of symmetrical eight-bladed β-propeller proteins

β-Propeller proteins form one of the largest families of protein structures, with a pseudo-symmetrical fold made up of subdomains called blades. They are not only abundant but are also involved in a wide variety of cellular processes, often by acting as a platform for the assembly of protein complexes. WD40 proteins are a subfamily of propeller proteins with no intrinsic enzymatic activity, but their stable, modular architecture and versatile surface have allowed evolution to adapt them to many vital roles. By computationally reverse-engineering the duplication, fusion and diversification events in the evolutionary history of a WD40 protein, a perfectly symmetrical homologue called Tako8 was made. If two or four blades of Tako8 are expressed as single polypeptides, they do not self-assemble to complete the eight-bladed architecture, which may be owing to the closely spaced negative charges inside the ring. A different computational approach was employed to redesign Tako8 to create Ika8, a fourfold-symmetrical protein in which neighbouring blades carry compensating charges. Ika2 and Ika4, carrying two or four blades per subunit, respectively, were found to assemble spontaneously into a complete eight-bladed ring in solution. These artificial eight-bladed rings may find applications in bionanotechnology and as models to study the folding and evolution of WD40 proteins.

Design of tryptophan-containing mutants of the symmetrical Pizza protein for biophysical studies

β-propeller proteins are highly symmetrical, being composed of a repeated motif with four anti-parallel β-sheets arranged around a central axis. Recently we designed the first completely symmetrical β-propeller protein, Pizza6, consisting of six identical tandem repeats. Pizza6 is expected to prove a useful building block for bionanotechnology, and also a tool to investigate the folding and evolution of β-propeller proteins. Folding studies are made difficult by the high stability and the lack of buried Trp residues to act as monitor fluorophores, so we have designed and characterized several Trp-containing Pizza6 derivatives. In total four proteins were designed, of which three could be purified and characterized. Crystal structures confirm these mutant proteins maintain the expected structure, and a clear redshift of Trp fluorescence emission could be observed upon denaturation. Among the derivative proteins, Pizza6-AYW appears to be the most suitable model protein for future folding/unfolding kinetics studies as it has a comparable stability as natural β-propeller proteins.

Intrinsic Disorder in a Well-Folded Globular Protein

The folded structure of the heterodimeric sweet protein monellin mimics single-chain proteins with topology β1-α1-β2-β3-β4-β5 (chain A: β3-β4-β5; chain B: β1-α1-β2). Furthermore, like naturally occurring single-chain proteins of a similar size, monellin folds cooperatively with no detectable intermediates. However, the two monellin chains, A and B, are marginally structured in isolation and fold only upon binding to each other. Thus, monellin presents a unique opportunity to understand the design of intrinsically disordered proteins that fold upon binding. Here, we study the folding of a single-chain variant of monellin (scMn) using simulations of an all heavy-atom structure-based model. These simulations can explain mechanistic details derived from scMn experiments performed using several different structural probes. scMn folds cooperatively in our structure-based simulations, as is also seen in experiments. We find that structure formation near the transition-state ensemble of scMn is not uniformly distributed but is localized to a hairpin-like structure which contains one strand from each chain (β2, β3). Thus, the sequence and the underlying energetics of heterodimeric monellin promote the early formation of the interchain interface (β2-β3). By studying computational scMn mutants whose "interchain" interactions are deleted, we infer that this energy distribution allows the two protein chains to remain largely disordered when this interface is not folded. From these results, we suggest that cutting the protein backbone of a globular protein between residues which lie within its folding nucleus may be one way to construct two disordered fragments which fold upon binding.

Kinetic stability of membrane proteins

Although membrane proteins constitute an important class of biomolecules involved in key cellular processes, study of the thermodynamic and kinetic stability of their structures is far behind that of soluble proteins. It is known that many membrane proteins become unstable when removed by detergent extraction from the lipid environment. In addition, most of them undergo irreversible denaturation, even under mild experimental conditions. This process was found to be associated with partial unfolding of the polypeptide chain exposing hydrophobic regions to water, and it was proposed that the formation of kinetically trapped conformations could be involved. In this review, we will describe some of the efforts toward understanding the irreversible inactivation of membrane proteins. Furthermore, its modulation by phospholipids, ligands, and temperature will be herein discussed.

Origin of the Reflectin Gene and Hierarchical Assembly of Its Protein

Cephalopods, the group of animals including octopus, squid, and cuttlefish, have remarkable ability to instantly modulate body coloration and patterns so as to blend into surrounding environments [1, 2] or send warning signals to other animals [3]. Reflectin is expressed exclusively in cephalopods, filling the lamellae of intracellular Bragg reflectors that exhibit dynamic iridescence and structural color change [4]. Here, we trace the possible origin of the reflectin gene back to a transposon from the symbiotic bioluminescent bacterium Vibrio fischeri and report the hierarchical structural architecture of reflectin protein. Intrinsic self-assembly, and higher-order assembly tightly modulated by aromatic compounds, provide insights into the formation of multilayer reflectors in iridophores and spherical microparticles in leucophores and may form the basis of structural color change in cephalopods. Self-assembly and higher-order assembly in reflectin originated from a core repeating octapeptide (here named protopeptide), which may be from the same symbiotic bacteria. The origin of the reflectin gene and assembly features of reflectin protein are of considerable biological interest. The hierarchical structural architecture of reflectin and its domain and protopeptide not only provide insights for bioinspired photonic materials but also serve as unique "assembly tags" and feasible molecular platforms in biotechnology.

Computational tools help improve protein stability but with a solubility tradeoff

Accurately predicting changes in protein stability upon amino acid substitution is a much sought after goal. Destabilizing mutations are often implicated in disease, whereas stabilizing mutations are of great value for industrial and therapeutic biotechnology. Increasing protein stability is an especially challenging task, with random substitution yielding stabilizing mutations in only ∼2% of cases. To overcome this bottleneck, computational tools that aim to predict the effect of mutations have been developed; however, achieving accuracy and consistency remains challenging. Here, we combined 11 freely available tools into a meta-predictor (meieringlab.uwaterloo.ca/stabilitypredict/). Validation against ∼600 experimental mutations indicated that our meta-predictor has improved performance over any of the individual tools. The meta-predictor was then used to recommend 10 mutations in a previously designed protein of moderate thermodynamic stability, ThreeFoil. Experimental characterization showed that four mutations increased protein stability and could be amplified through ThreeFoil's structural symmetry to yield several multiple mutants with >2-kcal/mol stabilization. By avoiding residues within functional ties, we could maintain ThreeFoil's glycan-binding capacity. Despite successfully achieving substantial stabilization, however, almost all mutations decreased protein solubility, the most common cause of protein design failure. Examination of the 600-mutation data set revealed that stabilizing mutations on the protein surface tend to increase hydrophobicity and that the individual tools favor this approach to gain stability. Thus, whereas currently available tools can increase protein stability and combining them into a meta-predictor yields enhanced reliability, improvements to the potentials/force fields underlying these tools are needed to avoid gaining protein stability at the cost of solubility.

Evidence for the principle of minimal frustration in the evolution of protein folding landscapes

Theoretical and experimental studies have firmly established that protein folding can be described by a funneled energy landscape. This funneled energy landscape is the result of foldable protein sequences evolving following the principle of minimal frustration, which allows proteins to rapidly fold to their native biologically functional conformations. For a protein family with a given functional fold, the principle of minimal frustration suggests that, independent of sequence, all proteins within this family should fold with similar rates. However, depending on the optimal living temperature of the organism, proteins also need to modulate their thermodynamic stability. Consequently, the difference in thermodynamic stability should be primarily caused by differences in the unfolding rates. To test this hypothesis experimentally, we performed comprehensive thermodynamic and kinetic analyses of 15 different proteins from the thioredoxin family. Eight of these thioredoxins were extant proteins from psychrophilic, mesophilic, or thermophilic organisms. The other seven protein sequences were obtained using ancestral sequence reconstruction and can be dated back over 4 billion years. We found that all studied proteins fold with very similar rates but unfold with rates that differ up to three orders of magnitude. The unfolding rates correlate well with the thermodynamic stability of the proteins. Moreover, proteins that unfold slower are more resistant to proteolysis. These results provide direct experimental support to the principle of minimal frustration hypothesis.

Kinetic Stability of Proteins in Beans and Peas: Implications for Protein Digestibility, Seed Germination, and Plant Adaptation

Kinetically stable proteins (KSPs) are resistant to the denaturing detergent sodium dodecyl sulfate (SDS). Such resilience makes KSPs resistant to proteolytic degradation and may have arisen in nature as a mechanism for organismal adaptation and survival against harsh conditions. Legumes are well-known for possessing degradation-resistant proteins that often diminish their nutritional value. Here we applied diagonal two-dimensional (D2D) SDS-polyacrylamide gel electrophoresis (PAGE), a method that allows for the proteomics-level identification of KSPs, to a group of 12 legumes (mostly beans and peas) of agricultural and nutritional importance. Our proteomics results show beans that are more difficult to digest, such as soybean, lima beans, and various common beans, have high contents of KSPs. In contrast, mung bean, red lentil, and various peas that are highly digestible contain low amounts of KSPs. Identified proteins with high kinetic stability are associated with warm-season beans, which germinate at higher temperatures. In contrast, peas and red lentil, which are cool-season legumes, contain low levels of KSPs. Thus, our results show protein kinetic stability is an important factor in the digestibility of legume proteins and may relate to nutrition efficiency, timing of seed germination, and legume resistance to biotic stressors. Furthermore, we show D2D SDS-PAGE is a powerful method that could be applied for determining the abundance and identity of KSPs in engineered and wild legumes and for advancing basic research and associated applications.

Using natural sequences and modularity to design common and novel protein topologies

Protein design is still a challenging undertaking, often requiring multiple attempts or iterations for success. Typically, the source of failure is unclear, and scoring metrics appear similar between successful and failed cases. Nevertheless, the use of sequence statistics, modularity and symmetry from natural proteins, combined with computational design both at the coarse-grained and atomistic levels is propelling a new wave of design efforts to success. Here we highlight recent examples of design, showing how the wealth of natural protein sequence and topology data may be leveraged to reduce the search space and increase the likelihood of achieving desired outcomes.

Quality over quantity: optimizing co-translational protein folding with non-'optimal' synonymous codons

Protein folding occurs on a time scale similar to peptide bond formation by the ribosome, which has long sparked speculation that altering translation rate could alter the folding mechanism or even the final folded structure of a protein in vivo. Recent results have provided strong support for this model: synonymous substitutions to codons with different usage frequency, which are often translated at different rates, have been shown to significantly alter the co-translational folding mechanism of some proteins, leading to altered cell function. Here we review recent progress towards understanding the connections between synonymous codon usage, translation rate and co-translational protein folding mechanisms.

How to Build a Complex, Functional Propeller Protein, From Parts

By combining ancestral sequence reconstruction and in vitro evolution, Smock et al. identified single motifs that assemble into a functional five-bladed β-propeller, and a likely route for conversion into the more complex, extant single chain fusion. Interestingly, although sequence diversification destabilized five-motif fusions, it also destabilized aggregation-prone intermediates, increasing the level of functional protein in vivo.