Communication between protein surfaces and their buried cores is central to protein structure and function, yet this phenomenon remains challenging to predict and control at high resolution. Changes in the protonation of surface ionizable residues communicate with the hydrophobic core, for example, in diverse pH-dependent protein functions. Hisactophilin, a histidine-rich actin- and membrane-binding protein, provides a general model for exploring such communication as it exhibits a finely tuned pH-regulated myristoyl-switching function. Upon reversible proton binding, the myristoyl group shifts between being sequestered in the hydrophobic core and more solvent accessible. In the current study we utilize experimental and computational approaches we uncover how binding of ~1.5 net protons alters electrostatic interactions involving ionizable residues distributed across much of the protein surface. These changes are transmitted to the hydrophobic core through dynamic communities of ionizable and hydrophobic residues which substantially rearrange upon switching. The effects of mutating individual ionizable residues are weaker than those of core hydrophobic residues, and only combined mutation of multiple ionizable residues caused substantial functional change. Together, these results reveal how communication between surface ionizable residues and the hydrophobic core is mediated by extensive interaction networks that reorganize in response to changes in protonation. These results may provide general insights for understanding protein cooperativity and the coupling of surface and core residues in protein function, disease, evolution, engineering, and design. Significance Statement:How changes on the protein surface, such as proton binding to ionizable amino acids, are communicated to the protein core to regulate protein stability and function remains ill-defined. Synthesis of experimental and computational analyses resolves the distributed networks of surface ionizable residue interactions coupled to the hydrophobic core that control pH-dependent myristoyl switching in hisactophilin. Small changes in protonation that create and alleviate local electrostatic repulsion give rise to protein-wide changes in fluctuating surface-core interactions. This distributed electrostatics-core coupling mechanism may help explain the often underrecognized and long-range impacts of ionizable residues in proteins and provide a framework for interpreting the effects of mutations in fundamental and applied protein science.

Abstract Protein aggregation is of broad importance in biotechnology and disease, yet the structural heterogeneity of cellular aggregates has confounded high-resolution structural analysis. Inclusion bodies (IBs) formed in Escherichia coli are an attractive, controllable system for unravelling the complexities of protein aggregation in a cellular context. Here, a multimodal analysis integrating residue-resolved quenched amide hydrogen–deuterium exchange (qHDX), proteolysis, FTIR, Congo red binding, and chemical denaturation is applied to IBs formed by proteins encompassing stable β- and α-globular folds, a partially structured protein fragment, and intrinsically disordered low complexity domains (LCDs). Remarkable conformational diversity is observed: IBs formed by well-folded proteins are extensively structured and include substantial local native-like features, whereas proteins with decreased access to stable native conformations form more heterogeneous and dynamic aggregates increasingly shaped by intrinsic sequence features. Strikingly, qHDX protection of TDP-43 LCD IBs strongly aligns with the core of cryo-EM structures of ex vivo pathological fibrils; however, peripheral regions that appear fully hydrogen-bonded in the cryo-EM structures exhibit little protection. The results reveal that individual protein IBs contain distinct mixtures of native-like, disordered, and amyloid-like conformers, informing the prediction and control of cellular aggregate structure and stability. Entry for the Table of Contents Inclusion body conformational heterogeneity. Left: aggregates retaining local native-like structure. Middle: assemblies populated by α-helical intermediates coexisting with antiparallel β-sheet oligomers. Right: initially disordered states undergo helix–helix interactions and progress towards β-sheet–rich conformations, illustrating distinct routes of cellular aggregation.

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.

We have recently reported that coadministration of sulpiride, an antipsychotic drug, and fluvoxamine, a selective serotonin (5-HT) reuptake inhibitor, selectively increases in vivo dopamine release in the prefrontal cortex. This study examined the effects of coadministration of these drugs on duration of immobility in the tail suspension test using mice. Neither sulpiride (3 or 10 mg/kg) nor fluvoxamine (10 or 20 mg/kg) alone affected immobility time, whereas coadministration significantly reduced immobility time. WAY 100635, a 5-HT(1A) receptor antagonist, did not affect the effects of sulpiride and fluvoxamine coadministration, but reduced immobility time in combination with fluvoxamine (20 mg/kg). A high dose of fluvoxamine alone (60 mg/kg) also reduced immobility time. These results suggest that the antidepressant-like effects of fluvoxamine in combination with sulpiride or WAY 100635 in the tail suspension test are mediated by the activation of dopamine or 5-HT systems, respectively.

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.

Protein misfolding and aggregation are hallmarks of many diseases, including amyotrophic lateral sclerosis (ALS). In familial ALS, aberrant self-association of mutant Cu,Zn-superoxide dismutase (SOD1) is implicated as a key contributor to disease. Mutations have the largest impacts on the stability of the most immature form of SOD1, the unmetallated, disulfide-reduced monomer (apoSH SOD1). Here we demonstrate that, despite the marginal stability of apoSH SOD1, aggregation is little correlated with the degree of protein unfolding, and multiple modes of aggregation occur, depending on the mutation and solution conditions. Light scattering and atomic force microscopy reveal two distinct mutant SOD1 behaviours: high aggregator mutants form abundant small assemblies, while low aggregator mutants form fewer, more fibre-like aggregates. Attenuated total reflectance-Fourier transform infrared spectroscopy and Thioflavin T binding show the aggregates maintain native-like anti-parallel beta structure. These results provide new evidence that ALS-associated mutations promote the aggregation of apoSH SOD1 through multiple pathways, with broad implications for understanding mechanisms of protein self-association in disease and biotechnology.

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.

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.

There has been a remarkable increase in the number of biologics, especially monoclonal antibodies, in the market over the last decade. In addition to attaining the desired binding to their targets, a crucial aspect is the 'developability' of these drugs, which includes several desirable properties such as high solubility, low viscosity and aggregation, physico-chemical stability, low immunogenicity and low poly-specificity. The lack of any of these desirable properties can lead to significant hurdles in advancing them to the clinic and are often discovered only during late stages of drug development. Hence, in silico methods for early detection of these properties, particularly the ones that affect aggregation and solubility in the earlier stages can be highly beneficial. We have developed a computational framework based on a large and diverse set of protein specific descriptors that is ideal for making liability predictions using a QSPR (quantitative structure-property relationship) approach. This set offers a high degree of feature diversity that may coarsely be classified based on (1) sequence (2) structure and (3) surface patches. We assess the sensitivity and applicability of these descriptors in four dedicated case studies that are believed to be representative of biophysical characterizations commonly employed during the development process of a biologics drug candidate. In addition to data sets obtained from public sources, we have validated the descriptors on novel experimental data sets in order to address antibody developability and to generate prospective predictions on Adnectins. The results show that the descriptors are well suited to assist in the improvement of protein properties of systems that exhibit poor solubility or aggregation.

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.

Inclusion bodies (IBs) – insoluble structures that can form upon overexpression of proteins- have wide relevance in research, industrial settings, and disease. In addition to their use in protein purification, IBs provide a tractable model to elucidate the molecular mechanisms of protein aggregation in cells, as also occurs in neurodegenerative disorders such as amyotrophic lateral sclerosis (ALS), Parkinson's, and Alzheimer's diseases. The aggregation of Cu, Zn superoxide dismutase (SOD1) variants is associated with ALS, yet the relationships between mutant characteristics and disease properties remain obscure. Here, we report a systematic investigation of SOD1 IBs using two powerful complementary methods to analyze the structures of these cellular aggregates and their changes upon mutation. Quenched amide hydrogen-deuterium exchange (HDX) measurements by NMR for individual residues throughout SOD1 reveal IB structural features and similarities, whereas quantitative conformation specific antibody binding assays highlight notable differences in surface features between mutant IBs. Taken together, these data provide a valuable new high-resolution view of cellular aggregation. These powerful HDX and quantitative antibody binding methods are widely applicable to other proteins to define the molecular determinants of their aggregation and solubility.

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 SOD1(G93A) rat model that misfolded SOD1 exists as distinct conformers and forms deposits on mitochondrial subpopulations. Here, using SOD1(G93A) 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.

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.

About 20% of all familial amyotrophic lateral sclerosis (ALS) cases are associated with mutations in superoxide dismutase (SOD1), a homodimeric protein. The disease has an autosomal-dominant inheritance pattern. It is, therefore, important to determine whether wild-type and mutant SOD1 subunits self-associate randomly or preferentially. A measure for the extent of bias in subunit association is the coupling constant determined in a double-mutant cycle type analysis. Here, cell lysates containing co-expressed wild-type and mutant SOD1 subunits were analyzed by native mass spectrometry to determine these coupling constants. Strikingly, we find a linear positive correlation between the coupling constant and the reported average duration of the disease. Our results indicate that inter-subunit communication and a preference for heterodimerization greatly increase the disease severity.

Kinetic stability is central to protein lifetime and survival under a wide range of challenging environmental conditions and is intricately related to protein function. Despite this centrality, the molecular determinants of kinetic stability are typically only crudely defined at a molecular level. We analyze the underpinnings of kinetic stability using proteins with a beta-trefoil superfold: the superstable designed multivalent carbohydrate binding protein, ThreeFoil, and the myristoylated pH-regulated histidine-rich actin-binding protein, Hisactophilin. Despite adopting the same fold, ThreeFoil has extreme kinetic stability and moderate thermodynamic stability whereas Hisactophlin is much more dynamic while having higher thermodynamic stability. The differences in stability are related to distinct requirements for function and to structural features that can be generalized to molecular mechanisms applying to a wide variety of proteins. These findings have broad implications for understanding natural proteins and for protein design.

Accurate modeling of the effects of mutations on protein stability is central to understanding and controlling proteins in myriad natural and applied contexts. Here, we reveal through rigorous quantitative analysis that stability prediction tools often favor mutations that increase stability at the expense of solubility. Moreover, while these tools may accurately identify strongly destabilizing mutations, the experimental effect of mutations predicted to stabilize is actually near neutral on average. The commonly used “classification accuracy” metric obscures this reality; accordingly, we recommend performance measures, such as the Matthews correlation coefficient (MCC). We demonstrate that an absurdly simple machine-learning algorithm—a neural network of just two neurons—unexpectedly achieves high classification accuracy, but its inadequacies are revealed by a low MCC. Despite the above limitations, making multiple mutations markedly improves the prospects for achieving a stabilization target, and modest improvements in the precision of future tools may yield disproportionate gains.