A fine balance of hydrophobic-electrostatic communication pathways in a pH-switching protein.

Allostery is the phenomenon of coupling between distal binding sites in a protein. Such coupling is at the crux of protein function and regulation in a myriad of scenarios, yet determining the molecular mechanisms of coupling networks in proteins remains a major challenge. Here, we report mechanisms governing pH-dependent myristoyl switching in monomeric hisactophilin, whereby the myristoyl moves between a sequestered state, i.e., buried within the core of the protein, to an accessible state, in which the myristoyl has increased accessibility for membrane binding. Measurements of the pH and temperature dependence of amide chemical shifts reveal protein local structural stability and conformational heterogeneity that accompany switching. An analysis of these measurements using a thermodynamic cycle framework shows that myristoyl-proton coupling at the single-residue level exists in a fine balance and extends throughout the protein. Strikingly, small changes in the stereochemistry or size of core and surface hydrophobic residues by point mutations readily break, restore, or tune myristoyl switch energetics. Synthesizing the experimental results with those of molecular dynamics simulations illuminates atomistic details of coupling throughout the protein, featuring a large network of hydrophobic interactions that work in concert with key electrostatic interactions. The simulations were critical for discerning which of the many ionizable residues in hisactophilin are important for switching and identifying the contributions of nonnative interactions in switching. The strategy of using temperature-dependent NMR presented here offers a powerful, widely applicable way to elucidate the molecular mechanisms of allostery in proteins at high resolution.

Methodological advances and strategies for high resolution structure determination of cellular protein aggregates.

Aggregation of proteins is at the nexus of molecular processes crucial to aging, disease, and employing proteins for biotechnology and medical applications. There has been much recent progress in determining the structural features of protein aggregates that form in cells; yet, owing to prevalent heterogeneity in aggregation, many aspects remain obscure and often experimentally intractable to define. Here, we review recent results of structural studies for cell-derived aggregates of normally globular proteins, with a focus on high-resolution methods for their analysis and prediction. Complementary results obtained by solid-state NMR spectroscopy, FTIR spectroscopy and microspectroscopy, cryo-EM, and amide hydrogen/deuterium exchange measured by NMR and mass spectrometry, applied to bacterial inclusion bodies and disease inclusions, are uncovering novel information on in-cell aggregation patterns as well as great diversity in the structural features of useful and aberrant protein aggregates. Using these advances as a guide, this review aims to advise the reader on which combination of approaches may be the most appropriate to apply to their unique system.

Quenched hydrogen-deuterium amide exchange optimization for high-resolution structural analysis of cellular protein aggregates.

Inclusion bodies (IBs) are large, insoluble aggregates that often form during the overexpression of proteins in bacteria. These aggregates are of broad fundamental and practical significance, for recombinant protein preparation and due to their relevance to aggregation-related medical conditions and their recent emergence as promising functional nanomaterials. Despite their significance, high resolution knowledge of IB structure remains very limited. Such knowledge will advance understanding and control of IB formation and properties in myriad practical applications. Here, we report a detailed quenched hydrogen-deuterium amide exchange (qHDX) method with NMR readout to define the structure of IBs at the level of individual residues throughout the protein. Applying proper control of experimental conditions, such as sample pH, water content, temperature, and intrinsic rate of amide exchange, yields in depth results for these cellular protein aggregates. qHDX results illustrated for Cu, Zn superoxide dismutase 1 (SOD1) and Adnectins show their IBs include native-like structure and some but not all mutations alter IB structure.

High-Resolution NMR H/D Exchange of Human Superoxide Dismutase Inclusion Bodies Reveals Significant Native Features Despite Structural Heterogeneity.

Protein aggregation is central to aging, disease and biotechnology. While there has been recent progress in defining structural features of cellular protein aggregates, many aspects remain unclear due to heterogeneity of aggregates presenting obstacles to characterization. Here we report high-resolution analysis of cellular inclusion bodies (IBs) of immature human superoxide dismutase (SOD1) mutants using NMR quenched amide hydrogen/deuterium exchange (qHDX), FTIR and Congo red binding. The extent of aggregation is correlated with mutant global stability and, notably, the free energy of native dimer dissociation, indicating contributions of native-like monomer associations to IB formation. This is further manifested by a common pattern of extensive protection against H/D exchange throughout nine mutant SOD1s despite their diverse characteristics. These results reveal multiple aggregation-prone regions in SOD1 and illuminate how aggregation may occur via an ensemble of pathways.

TNF receptor-associated factor 6 interacts with ALS-linked misfolded superoxide dismutase 1 and promotes aggregation.

Amyotrophic lateral sclerosis (ALS) is a fatal disease, characterized by the selective loss of motor neurons leading to paralysis. Mutations in the gene encoding superoxide dismutase 1 (SOD1) are the second most common cause of familial ALS, and considerable evidence suggests that these mutations result in an increase in toxicity due to protein misfolding. We previously demonstrated in the SOD1G93A rat model that misfolded SOD1 exists as distinct conformers and forms deposits on mitochondrial subpopulations. Here, using SOD1G93A rats and conformation-restricted antibodies specific for misfolded SOD1 (B8H10 and AMF7-63), we identified the interactomes of the mitochondrial pools of misfolded SOD1. This strategy identified binding proteins that uniquely interacted with either AMF7-63 or B8H10-reactive SOD1 conformers as well as a high proportion of interactors common to both conformers. Of this latter set, we identified the E3 ubiquitin ligase TNF receptor-associated factor 6 (TRAF6) as a SOD1 interactor, and we determined that exposure of the SOD1 functional loops facilitates this interaction. Of note, this conformational change was not universally fulfilled by all SOD1 variants and differentiated TRAF6 interacting from TRAF6 noninteracting SOD1 variants. Functionally, TRAF6 stimulated polyubiquitination and aggregation of the interacting SOD1 variants. TRAF6 E3 ubiquitin ligase activity was required for the former but was dispensable for the latter, indicating that TRAF6-mediated polyubiquitination and aggregation of the SOD1 variants are independent events. We propose that the interaction between misfolded SOD1 and TRAF6 may be relevant to the etiology of ALS.

Wild-type and mutant SOD1 localizes to RNA-rich structures in cells and mice but does not bind RNA.

Many of the genes whose mutation causes Amyotrophic Lateral Sclerosis (ALS) are RNA-binding proteins which localize to stress granules, while others impact the assembly, stability, and elimination of stress granules. This has led to the hypothesis that alterations in the dynamics of stress granules and RNA biology cause ALS. Genetic mutations in Superoxide Dismutase 1 (SOD1) also cause ALS. Evidence demonstrates that SOD1 harboring ALS-linked mutations is recruited to stress granules, induces changes in alternative splicing, and could be an RNA-binding protein. Whether SOD1 inclusions contain RNA in disease models and whether SOD1 directly binds RNA remains uncertain. We applied methods including cross-linking immunoprecipitation and in vitro gel shift assays to detect binding of SOD1 to RNA in vitro, in cells with and without stress granules, and in mice expressing human SOD1 G93A. We find that SOD1 localizes to RNA-rich structures including stress granules, and SOD1 inclusions in mice contain mRNA. However, we find no evidence that SOD1 directly binds RNA. This suggests that SOD1 may impact stress granules, alternative splicing and RNA biology without binding directly to RNA.

Effects of maturation on the conformational free-energy landscape of SOD1.

Amyotrophic lateral sclerosis (ALS) is a devastating fatal syndrome characterized by very rapid degeneration of motor neurons. A leading hypothesis is that ALS is caused by toxic protein misfolding and aggregation, as also occurs in many other neurodegenerative disorders, such as prion, Alzheimer's, Parkinson's, and Huntington's diseases. A prominent cause of familial ALS is mutations in the protein superoxide dismutase (SOD1), which promote the formation of misfolded SOD1 conformers that are prone to aberrant interactions both with each other and with other cellular components. We have shown previously that immature SOD1, lacking bound Cu and Zn metal ions and the intrasubunit disulfide bond (apoSOD12SH), has a rugged free-energy surface (FES) and exchanges with four other conformations (excited states) that have millisecond lifetimes and sparse populations on the order of a few percent. Here, we examine further states of SOD1 along its maturation pathway, as well as those off-pathway resulting from metal loss that have been observed in proteinaceous inclusions. Metallation and disulfide bond formation lead to structural transformations including local ordering of the electrostatic loop and native dimerization that are observed in rare conformers of apoSOD12SH; thus, SOD1 maturation may occur via a population-switch mechanism whereby posttranslational modifications select for preexisting structures on the FES. Metallation and oxidation of SOD1 stabilize the native, mature conformation and decrease the number of detected excited conformational states, suggesting that it is the immature forms of the protein that contribute to misfolded conformations in vivo rather than the highly stable enzymatically active dimer.

Spectrophotometric method for simultaneous measurement of zinc and copper in metalloproteins using 4-(2-pyridylazo)resorcinol.

Bound metals are observed in a great many natural proteins, where they perform diverse roles in determining protein folding, stability and function. Due to the broad impact of bound metals on biophysical and biochemical properties of proteins, it is valuable to have accurate and facile methods for determining the metal content of proteins. Here we describe an optimized methodology using 4-(2-pyridylazo)resorcinol (PAR) to simultaneously quantify two metal ions in solution. The assay is demonstrated for quantification of Cu2+ and Zn2+ ions in human Cu, Zn superoxide dismutases (SOD1s); however, the method is general and can be applied to various combinations of metal ions. Advantages of the assay are that it is rapid and inexpensive, requires little sample and preparation, and has simple data analysis. We show that spectral decomposition software can accurately resolve the absorption bands of Cu2+ and Zn2+ with high accuracy and precision. Using the PAR assay, we determined that metal binding is altered in disease-associated mutants of SOD1, with comparable results to those determined by ICP-AES. In addition, we highlight key issues for using spectrophotometric chelators such as PAR for metal analysis of proteins.

Probing the free energy landscapes of ALS disease mutants of SOD1 by NMR spectroscopy.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease that, in some cases, has been linked with mutations to the antioxidant metalloenzyme superoxide dismutase (SOD1). Although the mature form of this enzyme is highly stable and resistant to aggregation, the most immature form, lacking metal and a stabilizing intrasubunit disulfide bond, apoSOD12SH, is dynamic and hypothesized to be a major cause of toxicity in vivo. Previous solution NMR studies of wild-type apoSOD12SH have shown that the ground state interconverts with a series of sparsely populated and transiently formed conformers, some of which have aberrant nonnative structures. Here, we study seven disease mutants of apoSOD12SH and characterize their free energy landscapes as a first step in understanding the initial stages of disease progression and, more generally, to evaluate the plasticity of low-lying protein conformational states. The mutations lead to little change in the structures and dynamics of the ground states of the mutant proteins. By contrast, the numbers of low-lying excited states that are accessible to each of the disease mutants can vary significantly, with additional conformers accessed in some cases. Our study suggests that the diversity of these structures can provide alternate interaction motifs for different mutants, establishing additional pathways for new and often aberrant intra- and intermolecular contacts. Further, it emphasizes the potential importance of conformationally excited states in directing both folding and misfolding processes.

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.

Designed protein reveals structural determinants of extreme kinetic stability.

The design of stable, functional proteins is difficult. Improved design requires a deeper knowledge of the molecular basis for design outcomes and properties. We previously used a bioinformatics and energy function method to design a symmetric superfold protein composed of repeating structural elements with multivalent carbohydrate-binding function, called ThreeFoil. This and similar methods have produced a notably high yield of stable proteins. Using a battery of experimental and computational analyses we show that despite its small size and lack of disulfide bonds, ThreeFoil has remarkably high kinetic stability and its folding is specifically chaperoned by carbohydrate binding. It is also extremely stable against thermal and chemical denaturation and proteolytic degradation. We demonstrate that the kinetic stability can be predicted and modeled using absolute contact order (ACO) and long-range order (LRO), as well as coarse-grained simulations; the stability arises from a topology that includes many long-range contacts which create a large and highly cooperative energy barrier for unfolding and folding. Extensive data from proteomic screens and other experiments reveal that a high ACO/LRO is a general feature of proteins with strong resistances to denaturation and degradation. These results provide tractable approaches for predicting resistance and designing proteins with sufficient topological complexity and long-range interactions to accommodate destabilizing functional features as well as withstand chemical and proteolytic challenge.

Asymmetric Anchoring Is Required for Efficient Ω-Loop Opening and Closing in Cytosolic Phosphoenolpyruvate Carboxykinase.

Mobile Ω-loops play essential roles in the function of many enzymes. Here we investigated the importance of a residue lying outside of the mobile Ω-loop element in the catalytic function of an H477R variant of cytosolic phosphoenolpyruvate carboxykinase using crystallographic, kinetic, and computational analysis. The crystallographic data suggest that the efficient transition of the Ω-loop to the closed conformation requires stabilization of the N-terminus of the loop through contacts between R461 and E588. In contrast, the C-terminal end of the Ω-loop undergoes changing interactions with the enzyme body through contacts between H477 at the C-terminus of the loop and E591 located on the enzyme body. Potential of mean force calculations demonstrated that altering the anchoring of the C-terminus of the Ω-loop via the H477R substitution results in the destabilization of the closed state of the Ω-loop by 3.4 kcal mol-1. The kinetic parameters for the enzyme were altered in an asymmetric fashion with the predominant effect being observed in the direction of oxaloacetate synthesis. This is exemplified by a reduction in kcat for the H477R mutant by an order of magnitude in the direction of OAA synthesis, while in the direction of PEP synthesis, it decreased by a factor of only 2. The data are consistent with a mechanism for loop conformational exchange between open and closed states in which a balance between fixed anchoring of the N-terminus of the Ω-loop and a flexible, unattached C-terminus drives the transition between a disordered (open) state and an ordered (closed) state.

Concurrent Increases and Decreases in Local Stability and Conformational Heterogeneity in Cu, Zn Superoxide Dismutase Variants Revealed by Temperature-Dependence of Amide Chemical Shifts.

The chemical shifts of backbone amide protons in proteins are sensitive reporters of local structural stability and conformational heterogeneity, which can be determined from their readily measured linear and nonlinear temperature-dependences, respectively. Here we report analyses of amide proton temperature-dependences for native dimeric Cu, Zn superoxide dismutase (holo pWT SOD1) and structurally diverse mutant SOD1s associated with amyotrophic lateral sclerosis (ALS). Holo pWT SOD1 loses structure with temperature first at its periphery and, while having extremely high global stability, nevertheless exhibits extensive conformational heterogeneity, with ∼1 in 5 residues showing evidence for population of low energy alternative states. The holo G93A and E100G ALS mutants have moderately decreased global stability, whereas V148I is slightly stabilized. Comparison of the holo mutants as well as the marginally stable immature monomeric unmetalated and disulfide-reduced (apo(2SH)) pWT with holo pWT shows that changes in the local structural stability of individual amides vary greatly, with average changes corresponding to differences in global protein stability measured by differential scanning calorimetry. Mutants also exhibit altered conformational heterogeneity compared to pWT. Strikingly, substantial increases as well as decreases in local stability and conformational heterogeneity occur, in particular upon maturation and for G93A. Thus, the temperature-dependence of amide shifts for SOD1 variants is a rich source of information on the location and extent of perturbation of structure upon covalent changes and ligand binding. The implications for potential mechanisms of toxic misfolding of SOD1 in disease and for general aspects of protein energetics, including entropy-enthalpy compensation, are discussed.

Combined Isothermal Titration and Differential Scanning Calorimetry Define Three-State Thermodynamics of fALS-Associated Mutant Apo SOD1 Dimers and an Increased Population of Folded Monomer.

Many proteins are naturally homooligomers, homodimers most frequently. The overall stability of oligomeric proteins may be described in terms of the stability of the constituent monomers and the stability of their association; together, these stabilities determine the populations of different monomer and associated species, which generally have different roles in the function or dysfunction of the protein. Here we show how a new combined calorimetry approach, using isothermal titration calorimetry to define monomer association energetics together with differential scanning calorimetry to measure total energetics of oligomer unfolding, can be used to analyze homodimeric unmetalated (apo) superoxide dismutase (SOD1) and determine the effects on the stability of structurally diverse mutations associated with amyotrophic lateral sclerosis (ALS). Despite being located throughout the protein, all mutations studied weaken the dimer interface, while concomitantly either decreasing or increasing the marginal stability of the monomer. Analysis of the populations of dimer, monomer, and unfolded monomer under physiological conditions of temperature, pH, and protein concentration shows that all mutations promote the formation of folded monomers. These findings may help rationalize the key roles proposed for monomer forms of SOD1 in neurotoxic aggregation in ALS, as well as roles for other forms of SOD1. Thus, the results obtained here provide a valuable approach for the quantitative analysis of homooligomeric protein stabilities, which can be used to elucidate the natural and aberrant roles of different forms of these proteins and to improve methods for predicting protein stabilities.

Evolution of magnetization due to asymmetric dimerization: theoretical considerations and application to aberrant oligomers formed by apoSOD1(2SH).

A set of coupled differential equations is presented describing the evolution of magnetization due to an exchange reaction whereby a pair of identical monomers form an asymmetric dimer. In their most general form the equations describe a three-site exchange process that reduces to two-site exchange under certain limiting conditions that are discussed. An application to the study of sparsely populated, transiently formed sets of aberrant dimers, symmetric and asymmetric, of superoxide dismutase is presented. Fits of concentration dependent CPMG relaxation dispersion profiles provide measures of the dimer dissociation constants and both on- and off-rates. Dissociation constants on the order of 70 mM are extracted from fits of the data, with dimeric populations of ∼2% and lifetimes of ∼6 and ∼2 ms for the symmetric and asymmetric complexes, respectively. This work emphasizes the important role that NMR relaxation experiments can play in characterizing very weak molecular complexes that remain invisible to most biophysical approaches.

Nonnative interactions regulate folding and switching of myristoylated protein.

We present an integrated experimental and computational study of the molecular mechanisms by which myristoylation affects protein folding and function, which has been little characterized to date. Myristoylation, the covalent linkage of a hydrophobic C14 fatty acyl chain to the N-terminal glycine in a protein, is a common modification that plays a critical role in vital regulated cellular processes by undergoing reversible energetic and conformational switching. Coarse-grained folding simulations for the model pH-dependent actin- and membrane-binding protein hisactophilin reveal that nonnative hydrophobic interactions of the myristoyl with the protein as well as nonnative electrostatic interactions have a pronounced effect on folding rates and thermodynamic stability. Folding measurements for hydrophobic residue mutations of hisactophilin and atomistic simulations indicate that the nonnative interactions of the myristoyl group in the folding transition state are nonspecific and robust, and so smooth the energy landscape for folding. In contrast, myristoyl interactions in the native state are highly specific and tuned for sensitive control of switching functionality. Simulations and amide hydrogen exchange measurements provide evidence for increases as well as decreases in stability localized on one side of the myristoyl binding pocket in the protein, implicating strain and altered dynamics in switching. The effects of folding and function arising from myristoylation are profoundly different from the effects of other post-translational modifications.

Decreased stability and increased formation of soluble aggregates by immature superoxide dismutase do not account for disease severity in ALS.

Protein aggregation is a hallmark of many diseases, including amyotrophic lateral sclerosis (ALS), where aggregation of Cu/Zn superoxide dismutase (SOD1) is implicated in causing neurodegeneration. Recent studies have suggested that destabilization and aggregation of the most immature form of SOD1, the disulfide-reduced, unmetallated (apo) protein is particularly important in causing ALS. We report herein in depth analyses of the effects of chemically and structurally diverse ALS-associated mutations on the stability and aggregation of reduced apo SOD1. In contrast with previous studies, we find that various reduced apo SOD1 mutants undergo highly reversible thermal denaturation with little aggregation, enabling quantitative thermodynamic stability analyses. In the absence of ALS-associated mutations, reduced apo SOD1 is marginally stable but predominantly folded. Mutations generally result in slight decreases to substantial increases in the fraction of unfolded protein. Calorimetry, ultracentrifugation, and light scattering show that all mutations enhance aggregation propensity, with the effects varying widely, from subtle increases in most cases, to pronounced formation of 40-100 nm soluble aggregates by A4V, a mutation that is associated with particularly short disease duration. Interestingly, although there is a correlation between observed aggregation and stability, there is minimal to no correlation between observed aggregation, predicted aggregation propensity, and disease characteristics. These findings suggest that reduced apo SOD1 does not play a dominant role in modulating disease. Rather, additional and/or multiple forms of SOD1 and additional biophysical and biological factors are needed to account for the toxicity of mutant SOD1 in ALS.

Energetics of oligomeric protein folding and association.

In nature, proteins most often exist as complexes, with many of these consisting of identical subunits. Understanding of the energetics governing the folding and misfolding of such homooligomeric proteins is central to understanding their function and misfunction, in disease or biotechnology. Much progress has been made in defining the mechanisms and thermodynamics of homooligomeric protein folding. In this review, we outline models as well as calorimetric and spectroscopic methods for characterizing oligomer folding, and describe extensive results obtained for diverse proteins, ranging from dimers to octamers and higher order aggregates. To our knowledge, this area has not been reviewed comprehensively in years, and the collective progress is impressive. The results provide evolutionary insights into the development of subunit interfaces, mechanisms of oligomer folding, and contributions of oligomerization to protein stability, function and regulation. Thermodynamic analyses have also proven valuable for understanding protein misfolding and aggregation mechanisms, suggesting new therapeutic avenues. Successful recent designs of novel, functional proteins demonstrate increased understanding of oligomer folding. Further rigorous analyses using multiple experimental and computational approaches are still required, however, to achieve consistent and accurate prediction of oligomer folding energetics. Modeling the energetics remains challenging but is a promising avenue for future advances.

Energetics and mechanisms of folding and flipping the myristoyl switch in the {beta}-trefoil protein, hisactophilin.

Myristoylation, the covalent linkage of a saturated, C(14) fatty acyl chain to the N-terminal glycine in a protein, plays a vital role in reversible membrane binding and signaling by the modified proteins. Currently, little is known about the effects of myristoylation on protein folding and stability, or about the energetics and molecular mechanisms of switching involving states with sequestered versus accessible myristoyl group. Our analysis of these effects in hisactophilin, a histidine-rich protein that binds cell membranes and actin in a pH-dependent manner, shows that myristoylation significantly increases hisactophilin stability, while also markedly increasing global protein folding and unfolding rates. The switching between sequestered and accessible states is pH dependent, with an apparent pK(switch) of 6.95, and an apparent free energy change of 2.0 kcal·mol(-1). The myristoyl switch is linked to the reversible uptake of ∼1.5 protons, likely by histidine residues. This pH dependence of switching appears to be the physical basis of the sensitive, pH-dependent regulation of membrane binding observed in vivo. We conclude that an increase in protein stability upon modification and burial of the attached group is likely to occur in numerous proteins modified with fatty acyl or other hydrophobic groups, and that the biophysical effects of such modification are likely to play an important role in their functional switches. In addition, the increased global dynamics caused by myristoylation of hisactophilin reveals a general mechanism whereby hydrophobic moieties can make nonnative interactions or relieve strain in transition states, thereby increasing the rates of interconversion between different states.

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.