Biological Roles of Protein Kinetic Stability

Proteolytic stability and aggregation in a key metabolic enzyme of bacteria

Proteins that are kinetically stable are thought to be less prone to both aggregation and proteolysis. We demonstrate that the classical lac system of Escherichia coli can be leveraged as a model system to study this relation. β-galactosidase (LacZ) plays a critical role in lactose metabolism and is an extremely stable protein that can persist in growing cells for multiple generations after expression has stopped. By attaching degradation tags to the LacZ protein, we find that LacZ can be transiently degraded during lac operon expression but once expression has stopped functional LacZ is protected from degradation. We reversibly destabilize its tetrameric assembly using α-complementation, and show that unassembled LacZ monomers and dimers can either be degraded or lead to formation of aggregates within cells, while the tetrameric state protects against proteolysis and aggregation. We show that the presence of aggregates is associated with cell death, and that these proteotoxic stress phenotypes can be alleviated by attaching an ssrA tag to LacZ monomers which leads to their degradation. We unify our findings using a biophysical model that enables the interplay of protein assembly, degradation, and aggregation to be studied quantitatively in vivo. This work may yield approaches to reversing and preventing protein-misfolding disease states, while elucidating the functions of proteolytic stability in constant and fluctuating environments.

Protein engineering using circular permutation - structure, function, stability, and applications

Protein engineering is important for creating novel variants from natural proteins, enabling a wide range of applications. Approaches such as rational design and directed evolution are routinely used to make new protein variants. Computational tools like de novo design can introduce new protein folds. Expanding the amino acid repertoire to include unnatural amino acids with non-canonical side chains in vitro by native chemical ligation and in vivo via codon expansion methods broadens sequence and structural possibilities. Circular permutation (CP) is an invaluable approach to redesigning a protein by rearranging the amino acid sequence, where the connectivity of the secondary structural elements is altered without changing the overall structure of the protein. Artificial CP proteins (CPs) are employed in various applications such as biocatalysis, sensing of small molecules by fluorescence, genome editing, ligand-binding protein switches, and optogenetic engineering. Many studies have shown that CP can lead to either reduced or enhanced stability or catalytic efficiency. The effects of CP on a protein's energy landscape cannot be predicted a priori. Thus, it is important to understand how CP can affect the thermodynamic and kinetic stability of a protein. In this review, we discuss the discovery and advancement of techniques to create protein CP, and existing reviews on CP. We delve into the plethora of biological applications for designed CP proteins. We subsequently discuss the experimental and computational reports on the effects of CP on the thermodynamic and kinetic stabilities of proteins of various topologies. An understanding of the various aspects of CP will allow the reader to design robust CP proteins for their specific purposes.

From Milliseconds to Minutes: Melittin Self-Assembly from Concerted Non-Equilibrium Pressure-Jump and Equilibrium Relaxation Nuclear Magnetic Resonance

Non-equilibrium kinetics techniques like pressure-jump nuclear magnetic resonance (NMR) are powerful in tracking changes in oligomeric populations and are not limited by relaxation rates for the time scales of exchange that can be probed. However, these techniques are less sensitive to minor, transient populations than are Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments. We integrated non-equilibrium pressure-jump and equilibrium CPMG relaxation dispersion data to fully map the kinetic landscape of melittin tetramerization. While monomeric peptides weakly form dimers (Kd,D/M ≈ 26 mM) whose population never exceeds 1.6% at 288 K, dimers associate tightly to form stable tetrameric species (Kd,T/D ≈ 740 nM). Exchange between the monomer and dimer, along with exchange between the dimer and tetramer, occurs on the millisecond time scale. The NMR approach developed herein can be readily applied to studying the folding and misfolding of a wide range of oligomeric assemblies.

Correlation of In Vitro Kinetic Stability to Preclinical In Vivo Pharmacokinetics for a Panel of Anti-PD-1 Monoclonal Antibody Interleukin 21 Mutein Immunocytokines

The development of therapeutic fusion protein drugs is often impeded by the unintended consequences that occur from fusing together domains from independent naturally occurring proteins, consequences such as altered biodistribution, tissue uptake, or rapid clearance and potential immunogenicity. For therapeutic fusion proteins containing globular domains, we hypothesized that aberrant in vivo behavior could be related to low kinetic stability of these domains leading to local unfolding and susceptibility to partial proteolysis and/or salvage and uptake. Herein we describe an assay to measure kinetic stability of therapeutic fusion proteins by way of their sensitivity to the protease thermolysin. The results indicate that in vivo pharmacokinetics of a panel of anti-programmed cell death protein 1 monocolonal antibody:interleukin 21 immunocytokines in both mice and nonhuman primates are highly correlated with their in vitro susceptibility to thermolysin-mediated proteolysis. This assay can be used as a tool to quickly identify in vivo liabilities of globular domains of therapeutic proteins, thus aiding in the optimization and development of new multispecific drug candidates. SIGNIFICANCE STATEMENT: This work describes a novel assay utilizing protein kinetic stability to identify preclinical in vivo pharmacokinetic liabilities of multispecific therapeutic fusion proteins. This provides an efficient, inexpensive method to ascertain inherent protein stability in vitro before conducting in vivo studies, which can rapidly increase the speed of preclinical drug development.

Competition between inside-out unfolding and pathogenic aggregation in an amyloid-forming β-propeller

Studies of folded-to-misfolded transitions using model protein systems reveal a range of unfolding needed for exposure of amyloid-prone regions for subsequent fibrillization. Here, we probe the relationship between unfolding and aggregation for glaucoma-associated myocilin. Mutations within the olfactomedin domain of myocilin (OLF) cause a gain-of-function, namely cytotoxic intracellular aggregation, which hastens disease progression. Aggregation by wild-type OLF (OLFWT) competes with its chemical unfolding, but only below the threshold where OLF loses tertiary structure. Representative moderate (OLFD380A) and severe (OLFI499F) disease variants aggregate differently, with rates comparable to OLFWT in initial stages of unfolding, and variants adopt distinct partially folded structures seen along the OLFWT urea-unfolding pathway. Whether initiated with mutation or chemical perturbation, unfolding propagates outward to the propeller surface. In sum, for this large protein prone to amyloid formation, the requirement for a conformational change to promote amyloid fibrillization leads to direct competition between unfolding and aggregation.

Circular permutation at azurin's active site slows down its folding

Circular permutation (CP) is a technique by which the primary sequence of a protein is rearranged to create new termini. The connectivity of the protein is altered but the overall protein structure generally remains unperturbed. Understanding the effect of CP can help design robust proteins for numerous applications such as in genetic engineering, optoelectronics, and improving catalytic activity. Studies on different protein topologies showed that CP usually affects protein stability as well as unfolding rates. Though a significant number of proteins contain metals or other cofactors, reports of metalloprotein CPs are rare. Thus, we chose a bacterial metalloprotein, azurin, and its CP within the metal-binding site (cpF114). We studied the stabilities, folding, and unfolding rates of apo- and Zn2+-bound CP azurin using fluorescence and circular dichroism. The introduced CP had destabilizing effects on the protein. Also, the folding of the Zn2+-CP protein was much slower than that of the Zn2+-WT or apo-protein. We compared this study to our previously reported azurin-cpN42, where we had observed an equilibrium and kinetic intermediate. cpF114 exhibits an apparent two-state equilibrium unfolding but has an off-pathway kinetic intermediate. Our study hinted at CP as a method to modify the energy landscape of proteins to alter their folding pathways. WT azurin, being a faster folder, may have evolved to optimize the folding rate of metal-bound protein compared to its CPs, albeit all of them have the same structure and function. Our study underscores that protein sequence and protein termini positions are crucial for metalloproteins. TOC Figure. (Top) Zn2+-azurin WT structure (PDB code: 1E67) and 2-D topology diagram of Zn2+-cpF114 azurin. (Bottom) Cartoon diagram representing folding (red arrows) and unfolding (blue arrows) of apo- and Zn2+- WT and cpF114 azurins. The width of the arrows represents the rate of the corresponding processes.

Molecular mechanisms of phase separation and amyloidosis of ALS/FTD-linked FUS and TDP-43

FUS and TDP-43, two RNA-binding proteins from the heterogeneous nuclear ribonucleoprotein family, have gained significant attention in the field of neurodegenerative diseases due to their association with amyotrophic lateral sclerosis (ALS) and frontotemporal degeneration (FTD). They possess folded domains for binding ATP and various nucleic acids including DNA and RNA, as well as substantial intrinsically disordered regions (IDRs) including prion-like domains (PLDs) and RG-/RGG-rich regions. They play vital roles in various cellular processes, including transcription, splicing, microRNA maturation, RNA stability and transport and DNA repair. In particular, they are key components for forming ribonucleoprotein granules and stress granules (SGs) through homotypic or heterotypic liquid-liquid phase separation (LLPS). Strikingly, liquid-like droplets formed by FUS and TDP-43 may undergo aging to transform into less dynamic assemblies such as hydrogels, inclusions, and amyloid fibrils, which are the pathological hallmarks of ALS and FTD. This review aims to synthesize and consolidate the biophysical knowledge of the sequences, structures, stability, dynamics, and inter-domain interactions of FUS and TDP-43 domains, so as to shed light on the molecular mechanisms underlying their liquid-liquid phase separation (LLPS) and amyloidosis. The review further delves into the mechanisms through which ALS-causing mutants of the well-folded hPFN1 disrupt the dynamics of LLPS of FUS prion-like domain, providing key insights into a potential mechanism for misfolding/aggregation-prone proteins to cause neurodegenerative diseases and aging by gain of functions. With better understanding of different biophysical aspects of FUS and TDP-43, the ultimate goal is to develop drugs targeting LLPS and amyloidosis, which could mediate protein homeostasis within cells and lead to new treatments for currently intractable diseases, particularly neurodegenerative diseases such as ALS, FTD and aging. However, the study of membrane-less organelles and condensates is still in its infancy and therefore the review also highlights key questions that require future investigation.

Phylogenetic variations in a novel family of hyperstable apple snail egg proteins: insights into structural stability and functional trends

The relationship between protein stability and function evolution has not been explored in proteins from natural sources. Here, we investigate the phylogenetic differences of Perivitellin-1 (PV1) a novel family of hyperstable egg carotenoproteins crucial to the reproductive success of Pomacea snails, as they have evolved clade-specific protective functions. We studied P. patula PV1 (PpaPV1) from Flagellata clade eggs, the most basal of Pomacea and compared it with PV1s orthologs from derived clades. PpaPV1 stands as the most stable, with longer unfolding half-life, resistance to detergent unfolding, and therefore higher kinetic stability than PV1s from derived clades. In fact, PpaPV1 is among the most hyperstable proteins described in nature. In addition, its spectral characteristics providing a pale egg coloration, mild lectin activity and glycan specificity are narrower than derived clades. Our results provide evidence indicating large structural and functional changes along the evolution of the genus. Notably, the lectin binding of PpaPV1 is less pronounced, and its glycan specificity is narrower compared to PV1s in the sister Bridgesii clade. Our findings underscore the phylogenetic disparities in terms of structural and kinetic stability, as well as defensive traits like a potent lectin activity affecting the gut morphology of potential predators within the Bridgesii clade or a conspicuous, likely warning coloration, within the Canaliculata clade. This work provides a comprehensive comparison of the structural attributes, stability profiles, and functional roles of apple snail egg PV1s from multiple species within a phylogenetic context. Furthermore, it proposes an evolutionary hypothesis suggesting a trade-off between structural stability and the functional aspects of apple snail's major egg defense protein.

In-cell kinetic stability is an essential trait in metallo-β-lactamase evolution

Protein stability is an essential property for biological function. In contrast to the vast knowledge on protein stability in vitro, little is known about the factors governing in-cell stability. Here we show that the metallo-β-lactamase (MBL) New Delhi MBL-1 (NDM-1) is a kinetically unstable protein on metal restriction that has evolved by acquiring different biochemical traits that optimize its in-cell stability. The nonmetalated (apo) NDM-1 is degraded by the periplasmic protease Prc that recognizes its partially unstructured C-terminal domain. Zn(II) binding renders the protein refractory to degradation by quenching the flexibility of this region. Membrane anchoring makes apo-NDM-1 less accessible to Prc and protects it from DegP, a cellular protease degrading misfolded, nonmetalated NDM-1 precursors. NDM variants accumulate substitutions at the C terminus that quench its flexibility, enhancing their kinetic stability and bypassing proteolysis. These observations link MBL-mediated resistance with the essential periplasmic metabolism, highlighting the importance of the cellular protein homeostasis.

Metal-binding and circular permutation-dependent thermodynamic and kinetic stability of azurin

Native topology is known to determine the folding kinetics and the energy landscape of proteins. Furthermore, the circular permutation (CP) of proteins alters the order of the secondary structure connectivity while retaining the three-dimensional structure, making it an elegant and powerful approach to altering native topology. Previous studies elucidated the influence of CP in proteins with different folds such as Greek key β-barrel, β-sandwich, β-α-β, and all α-Greek key. CP mainly affects the protein stability and unfolding kinetics, while folding kinetics remains mostly unaltered. However, the effect of CP on metalloproteins is yet to be elaborately studied. The active site of metalloproteins poses an additional complexity in studying protein folding. Here, we investigate a CP variant (cpN42) of azurin-in both metal-free and metal-bound (holo) forms. As observed earlier in other proteins, apo-forms of wild-type (WT) and cpN42 fold with similar rates. In contrast, zinc-binding accelerates the folding of WT but decelerates the folding of cpN42. On zinc-binding, the spontaneous folding rate of WT increases by >250 times that of cpN42, which is unprecedented and the highest for any CP to date. On the other hand, zinc-binding reduces the spontaneous unfolding rate of cpN42 by ~100 times, making the WT and CP azurins unfold at similar rates. Our study demonstrates metal binding as a novel way to modulate the unfolding and folding rates of CPs compared to their WT counterparts. We hope our study increases the understanding of the effect of CP on the folding mechanism and energy landscape of metalloproteins.

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.

Investigating the effects of cyclic topology on the performance of a plastic degrading enzyme for polyethylene terephthalate degradation

Agitation is a commonly encountered stress for enzymes during all stages of production and application, but investigations that aim to improve their tolerance using topological engineering have yet to be reported. Here, the plastic-degrading enzyme IsPETase was cyclized in a range of topologies including a cyclic monomer, cyclic dimer and catenane using SpyTag/SpyCatcher technologies, and their tolerance towards different stresses including mechanical agitation was investigated. The cyclic dimer and catenane topologies were less susceptible to agitation-induced inactivation resulting in enhancement of polyethylene terephthalate (PET) degradation. While contrary to conventional belief, cyclic topologies did not improve tolerance of IsPETase towards heat or proteolytic treatment, the close proximity of active sites in the dimeric and catenane variants was found to enhance PET conversion into small soluble products. Together, these findings illustrate that it is worthwhile to explore the topology engineering of enzymes used in heterogeneous catalysis as it improves factors that are often overlooked in homogeneous catalysis studies.

Structural and functional adaptation in extremophilic microbial α-amylases

Maintaining stable native conformation of a protein under a given ecological condition is the prerequisite for survival of organisms. Extremophilic bacteria and archaea have evolved to adapt under extreme conditions of temperature, pH, salt, and pressure. Molecular adaptations of proteins under these conditions are essential for their survival. These organisms have the capability to maintain stable, native conformations of proteins under extreme conditions. The enzymes produced by the extremophiles are also known as extremozyme, which are used in several industries. Stability and functionality of extremozymes under varying temperature, pH, and solvent conditions are the most desirable requirement of industry. α-Amylase is one of the most important enzymes used in food, pharmaceutical, textile, and detergent industries. This enzyme is produced by diverse microorganisms including various extremophiles. Therefore, understanding its stability is important from fundamental as well as an applied point of view. Each class of extremophiles has a distinctive set of dominant non-covalent interactions which are important for their stability. Static information obtained by comparative analysis of amino acid sequence and atomic resolution structure provides information on the prevalence of particular amino acids or a group of non-covalent interactions. Protein folding studies give the information about thermodynamic and kinetic stability in order to understand dynamic aspect of molecular adaptations. In this review, we have summarized information on amino acid sequence, structure, stability, and adaptability of α-amylases from different classes of extremophiles.

Confinements of Thermoresponsive Dendronized Polymers to Proteins

Crowding environment created by host polymers plays crucial roles in manipulating interactions with proteins and modulating their bioactivity. Here, we report our investigation on interactions between polymers and proteins in...

Expanding the toolbox for predictive parameters describing antibody stability considering thermodynamic and kinetic determinants

PURPOSE

Introduction of the activation energy (E_a) as a kinetic parameter to describe and discriminate monoclonal antibody (mAb) stability.

METHODS

E_a is derived from intrinsic fluorescence (IF) unfolding thermograms. An apparent irreversible three-state fit model based on the Arrhenius integral is developed to determine E_a of respective unfolding transitions. These activation energies are compared to the thermodynamic parameter of van´t Hoff enthalpies (∆H_vh). Using a set of 34 mAbs formulated in four different formulations, both the apparent thermodynamic and kinetic parameters together with apparent melting temperatures are correlated collectively with each other to storage stabilities to evaluate its predictive power with respect to long-term effects potentially reflected in shelf-life.

RESULTS

E_a allows for the discrimination of (i) different parent mAbs, (ii) different variants that originate from parent mAbs, and (iii) different formulations. Interestingly, we observed that the E_a of the CH₂ unfolding transition shows strongest correlations with monomer and aggregate content after storage at accelerated and stress conditions when collectively compared to ∆H_vh and T_m of the CH₂ transition. Moreover, the predictive parameters determined for the CH₂ domain show generally stronger correlations with monomer and aggregate content than those derived for the Fab. Qualitative assessment by ranking E_a of the Fab domain showed good agreement with monomer content in storage stabilities of individual mAb sub-sets.

CONCLUSION

E_a from IF unfolding transitions can be used in addition to other commonly used thermodynamic predictive parameters to discriminate and characterize thermal stability of different mAbs in different formulations. Hence, it shows great potential for antibody engineering and formulation scientists.

Transient disorder along pathways to amyloid

High-resolution structures of amyloid fibrils formed from normally-folded proteins have revealed non-native conformations of the polypeptide chains. Attaining these conformations apparently requires transition from the native state via a highly disordered conformation, in contrast to earlier models that posited a role for assembly of partially folded proteins. Modifications or interactions that extend the lifetime or constrain the conformations of these disordered states could act to enhance or suppress amyloid formation. Understanding how the properties of both the folded and transiently disordered structural ensembles influence the process of amyloid formation is a substantial challenge, but research into the properties of intrinsically disordered proteins will deliver important insights.

Protein analysis and stability: Overcoming trial-and-error by grouping according to physicochemical properties

Today proteins are possibly the most important class of substances. Yet new tasks for proteins are still often solved by trial-and-error approaches. However, in some areas these euphemistically called "screening approaches" are not suitable. E.g. stability tests just take too long and therefore require a more strategic, target-orientated concept. This concept is available by grouping proteins according to their physicochemical properties and then pulling out the right drawer for new tasks. These properties include size, then charge and hydrophobicity as well as their patchinesses, and the degree of order. In addition, solubility, the content of (free) enthalpy, aromatic-amino-acid- and α/β-frequency as well as helix capping, and corresponding patchiness, the number of specific motifs and domains as well as the typical concentration range can be helpful to discriminate between different groups of proteins. Analyzing correlations will reduce the necessary amount of parameters and additional ones, which may be still undiscovered at the present time, can be identified looking at protein subgroups with similar physicochemical properties which still behave heterogeneously. Step-by-step the methodology will be improved. Possibly protein stability will be the driver of this process, but all other areas such as production, purification and analytics including sample pre-treatment and the choice of appropriate separation conditions for e.g. chromatography and electrophoresis will profit from a rational strategy.

Exploring the Evolutionary History of Kinetic Stability in the α-Lytic Protease Family

In addition to encoding the tertiary fold and stability, the primary sequence of a protein encodes the folding trajectory and kinetic barriers that determine the speed of folding. How these kinetic barriers are encoded is not well understood. Here, we use evolutionary sequence variation in the α-lytic protease (αLP) protein family to probe the relationship between sequence and energy landscape. αLP has an unusual energy landscape: the native state of αLP is not the most thermodynamically favored conformation and, instead, remains folded due to a large kinetic barrier preventing unfolding. To fold, αLP utilizes an N-terminal pro region similar in size to the protease itself that functions as a folding catalyst. Once folded, the pro region is removed, and the native state does not unfold on a biologically relevant time scale. Without the pro region, αLP folds on the order of millennia. A phylogenetic search uncovers αLP homologs with a wide range of pro region sizes, including some with no pro region at all. In the resulting phylogenetic tree, these homologs cluster by pro region size. By studying homologs naturally lacking a pro region, we demonstrate they can be thermodynamically stable, fold much faster than αLP, yet retain the same fold as αLP. Key amino acids thought to contribute to αLP's extreme kinetic stability are lost in these homologs, supporting their role in kinetic stability. This study highlights how the entire energy landscape plays an important role in determining the evolutionary pressures on the protein sequence.

Gray code representation of the universal genetic code: Generation of never born protein sequences using Toeplitz matrix approach

In this paper, we identify all possible Gray Code and Partitioned Gray Code representations of the Universal Genetic Code for n = 2-bit and 3-bit binary numbers. We analyse the Hamming Distance matrices of all these Gray code and Partitioned Gray Code possibilities for which we obtain the Toeplitz and Partitioned Toeplitz Matrices, respectively. We use this Gray Code and Partitioned Gray Code representations of the Universal Genetic Code combined with the novel Toeplitz matrix approach to generate many Never Born Protein (NBP) Sequences, which exhibit intrinsic structural stability. In general, Never Born Protein sequences may have many potential applications in synthetic biology and opens a new vista in understanding this new subset of proteins for better applications in drug discovery, synthesis of fine chemicals, etc.

Non-conservation of folding rates in the thioredoxin family reveals degradation of ancestral unassisted-folding

Evolution involves not only adaptation, but also the degradation of superfluous features. Many examples of degradation at the morphological level are known (vestigial organs, for instance). However, the impact of degradation on molecular evolution has been rarely addressed. Thioredoxins serve as general oxidoreductases in all cells. Here, we report extensive mutational analyses on the folding of modern and resurrected ancestral bacterial thioredoxins. Contrary to claims from recent literature, in vitro folding rates in the thioredoxin family are not evolutionarily conserved, but span at least a ∼100-fold range. Furthermore, modern thioredoxin folding is often substantially slower than ancestral thioredoxin folding. Unassisted folding, as probed in vitro, thus emerges as an ancestral vestigial feature that underwent degradation, plausibly upon the evolutionary emergence of efficient cellular folding assistance. More generally, our results provide evidence that degradation of ancestral features shapes, not only morphological evolution, but also the evolution of individual proteins.

Comparative Structural Analysis of 20S Proteasome Ortholog Protein Complexes by Native Mass Spectrometry

Ortholog protein complexes are responsible for equivalent functions in different organisms. However, during evolution, each organism adapts to meet its physiological needs and the environmental challenges imposed by its niche. This selection pressure leads to structural diversity in protein complexes, which are often difficult to specify, especially in the absence of high-resolution structures. Here, we describe a multilevel experimental approach based on native mass spectrometry (MS) tools for elucidating the structural preservation and variations among highly related protein complexes. The 20S proteasome, an essential protein degradation machinery, served as our model system, wherein we examined five complexes isolated from different organisms. We show that throughout evolution, from the Thermoplasma acidophilum archaeal prokaryotic complex to the eukaryotic 20S proteasomes in yeast (Saccharomyces cerevisiae) and mammals (rat - Rattus norvegicus, rabbit - Oryctolagus cuniculus and human - HEK293 cells), the proteasome increased both in size and stability. Native MS structural signatures of the rat and rabbit 20S proteasomes, which heretofore lacked high-resolution, three-dimensional structures, highly resembled that of the human complex. Using cryoelectron microscopy single-particle analysis, we were able to obtain a high-resolution structure of the rat 20S proteasome, allowing us to validate the MS-based results. Our study also revealed that the yeast complex, and not those in mammals, was the largest in size and displayed the greatest degree of kinetic stability. Moreover, we also identified a new proteoform of the PSMA7 subunit that resides within the rat and rabbit complexes, which to our knowledge have not been previously described. Altogether, our strategy enables elucidation of the unique structural properties of protein complexes that are highly similar to one another, a framework that is valid not only to ortholog protein complexes, but also for other highly related protein assemblies.

Evaluating Protein Engineering Thermostability Prediction Tools Using an Independently Generated Dataset

Engineering proteins to enhance thermal stability is a widely utilized approach for creating industrially relevant biocatalysts. The development of new experimental datasets and computational tools to guide these engineering efforts remains an active area of research. Thus, to complement the previously reported measures of T 50 and kinetic constants, we are reporting an expansion of our previously published dataset of mutants for β-glucosidase to include both measures of T M and ΔΔG. For a set of 51 mutants, we found that T 50 and T M are moderately correlated, with a Pearson correlation coefficient and Spearman's rank coefficient of 0.58 and 0.47, respectively, indicating that the two methods capture different physical features. The performance of predicted stability using nine computational tools was also evaluated on the dataset of 51 mutants, none of which are found to be strong predictors of the observed changes in T 50, T M, or ΔΔG. Furthermore, the ability of the nine algorithms to predict the production of isolatable soluble protein was examined, which revealed that Rosetta ΔΔG, FoldX, DeepDDG, PoPMuSiC, and SDM were capable of predicting if a mutant could be produced and isolated as a soluble protein. These results further highlight the need for new algorithms for predicting modest, yet important, changes in thermal stability as well as a new utility for current algorithms for prescreening designs for the production of mutants that maintain fold and soluble production properties.

Stable calcium-free myocilin olfactomedin domain variants reveal challenges in differentiating between benign and glaucoma-causing mutations

Nonsynonymous gene mutations can be beneficial, neutral, or detrimental to the stability, structure, and biological function of the encoded protein, but the effects of these mutations are often not readily predictable. For example, the β-propeller olfactomedin domain of myocilin (mOLF) exhibits a complex interrelationship among structure(s), stability, and aggregation. Numerous mutations within mOLF are linked to glaucoma; the resulting variants are less stable, aggregation-prone, and sequestered intracellularly, causing cytotoxicity. Here, we report the first stable mOLF variants carrying substitutions in the calcium-binding site that exhibit solution characteristics indistinguishable from those of glaucoma variants. Crystal structures of these stable variants at 1.8-2.0-Å resolution revealed features that we could not predict by molecular dynamics simulations, including loss of loop structure, helix unwinding, and a blade shift. Double mutants that combined a stabilizing substitution and a selected glaucoma-causing single-point mutant rescued in vitro folding and stability defects. In the context of full-length myocilin, secretion of stable single variants was indistinguishable from that of the WT protein, and the double mutants were secreted to varying extents. In summary, our finding that mOLF can tolerate particular substitutions that render the protein stable despite a conformational switch emphasizes the complexities in differentiating between benign and glaucoma-causing variants and provides new insight into the possible biological function of myocilin.

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.

A strategy based on thermal flexibility to design triosephosphate isomerase proteins with increased or decreased kinetic stability

Kinetic stability of proteins determines their susceptibility to irreversibly unfold in a time-dependent process, and therefore its half-life. A residue displacement analysis of temperature-induced unfolding molecular dynamics simulations was recently employed to define the thermal flexibility of proteins. This property was found to be correlated with the activation energy barrier (E_act) separating the native from the transition state in the denaturation process. The E_act was determined from the application of a two-state irreversible model to temperature unfolding experiments using differential scanning calorimetry (DSC). The contribution of each residue to the thermal flexibility of proteins is used here to propose multiple mutations in triosephosphate isomerase (TIM) from Trypanosoma brucei (TbTIM) and Trypanosoma cruzi (TcTIM), two parasites closely related by evolution. These two enzymes, taken as model systems, have practically identical structure but large differences in their kinetic stability. We constructed two functional TIM variants with more than twice and less than half the activation energy of their respective wild-type reference structures. The results show that the proposed strategy is able to identify the crucial residues for the kinetic stability in these enzymes. As it occurs with other protein properties reflecting their complex behavior, kinetic stability appears to be the consequence of an extensive network of inter-residue interactions, acting in a concerted manner. The proposed strategy to design variants can be used with other proteins, to increase or decrease their functional half-life.

Evolutionary Divergent Suppressor Mutations in Conformational Diseases

Neutral and adaptive mutations are key players in the evolutionary dynamics of proteins at molecular, cellular and organismal levels. Conversely, largely destabilizing mutations are rarely tolerated by evolution, although their occurrence in diverse human populations has important roles in the pathogenesis of conformational diseases. We have recently proposed that divergence at certain sites from the consensus (amino acid) state during mammalian evolution may have rendered some human proteins more vulnerable towards disease-associated mutations, primarily by decreasing their conformational stability. We herein extend and refine this hypothesis discussing results from phylogenetic and structural analyses, structure-based energy calculations and structure-function studies at molecular and cellular levels. As proof-of-principle, we focus on different mammalian orthologues of the NQO1 (NAD(P)H:quinone oxidoreductase 1) and AGT (alanine:glyoxylate aminotransferase) proteins. We discuss the different loss-of-function pathogenic mechanisms associated with diseases involving the two enzymes, including enzyme inactivation, accelerated degradation, intracellular mistargeting, and aggregation. Last, we take into account the potentially higher robustness of mammalian orthologues containing certain consensus amino acids as suppressors of human disease, and their relation with different intracellular post-translational modifications and protein quality control capacities, to be discussed as sources of phenotypic variability between human and mammalian models of disease and as tools for improving current therapeutic approaches.