Shift of the insoluble content of the proteome in the aging mouse brain.

During aging, the cellular response to unfolded proteins is believed to decline, resulting in diminished proteostasis. In model organisms, such as Caenorhabditis elegans, proteostatic decline with age has been linked to proteome solubility shifts and the onset of protein aggregation. However, this correlation has not been extensively characterized in aging mammals. To uncover age-dependent changes in the insoluble portion of a mammalian proteome, we analyzed the detergent-insoluble fraction of mouse brain tissue by mass spectrometry. We identified a group of 171 proteins, including the small heat shock protein α-crystallin, that become enriched in the detergent-insoluble fraction obtained from old mice. To enhance our ability to detect features associated with proteins in that fraction, we complemented our data with a meta-analysis of studies reporting the detergent-insoluble proteins in various mouse models of aging and neurodegeneration. Strikingly, insoluble proteins from young and old mice are distinct in several features in our study and across the collected literature data. In younger mice, proteins are more likely to be disordered, part of membraneless organelles, and involved in RNA binding. These traits become less prominent with age, as an increased number of structured proteins enter the pellet fraction. This analysis suggests that age-related changes to proteome organization lead a group of proteins with specific features to become detergent-insoluble. Importantly, these features are not consistent with those associated with proteins driving membraneless organelle formation. We see no evidence in our system of a general increase of condensate proteins in the detergent-insoluble fraction with age.

NAD depletion mediates cytotoxicity in human neurons with autophagy deficiency.

Autophagy is a homeostatic process critical for cellular survival, and its malfunction is implicated in human diseases including neurodegeneration. Loss of autophagy contributes to cytotoxicity and tissue degeneration, but the mechanistic understanding of this phenomenon remains elusive. Here, we generated autophagy-deficient (ATG5-/-) human embryonic stem cells (hESCs), from which we established a human neuronal platform to investigate how loss of autophagy affects neuronal survival. ATG5-/- neurons exhibit basal cytotoxicity accompanied by metabolic defects. Depletion of nicotinamide adenine dinucleotide (NAD) due to hyperactivation of NAD-consuming enzymes is found to trigger cell death via mitochondrial depolarization in ATG5-/- neurons. Boosting intracellular NAD levels improves cell viability by restoring mitochondrial bioenergetics and proteostasis in ATG5-/- neurons. Our findings elucidate a mechanistic link between autophagy deficiency and neuronal cell death that can be targeted for therapeutic interventions in neurodegenerative and lysosomal storage diseases associated with autophagic defect.

Comparison of Biomolecular Condensate Localization and Protein Phase Separation Predictors.

Research in the field of biochemistry and cellular biology has entered a new phase due to the discovery of phase separation driving the formation of biomolecular condensates, or membraneless organelles, in cells. The implications of this novel principle of cellular organization are vast and can be applied at multiple scales, spawning exciting research questions in numerous directions. Of fundamental importance are the molecular mechanisms that underly biomolecular condensate formation within cells and whether insights gained into these mechanisms provide a gateway for accurate predictions of protein phase behavior. Within the last six years, a significant number of predictors for protein phase separation and condensate localization have emerged. Herein, we compare a collection of state-of-the-art predictors on different tasks related to protein phase behavior. We show that the tested methods achieve high AUCs in the identification of biomolecular condensate drivers and scaffolds, as well as in the identification of proteins able to phase separate in vitro. However, our benchmark tests reveal that their performance is poorer when used to predict protein segments that are involved in phase separation or to classify amino acid substitutions as phase-separation-promoting or -inhibiting mutations. Our results suggest that the phenomenological approach used by most predictors is insufficient to fully grasp the complexity of the phenomenon within biological contexts and make reliable predictions related to protein phase behavior at the residue level.

Pulse labeling reveals the tail end of protein folding by proteome profiling.

Accurate and efficient folding of nascent protein sequences into their native states requires support from the protein homeostasis network. Herein we probe which newly translated proteins are thermo-sensitive, making them susceptible to misfolding and aggregation under heat stress using pulse-SILAC mass spectrometry. We find a distinct group of proteins that is highly sensitive to this perturbation when newly synthesized but not once matured. These proteins are abundant and highly structured. Notably, they display a tendency to form ß sheet secondary structures, have more complex folding topology, and are enriched for chaperone-binding motifs, suggesting a higher demand for chaperone-assisted folding. These polypeptides are also more often components of stable protein complexes in comparison with other proteins. Combining these findings suggests the existence of a specific subset of proteins in the cell that is particularly vulnerable to misfolding and aggregation following synthesis before reaching the native state.

GraPES: The Granule Protein Enrichment Server for prediction of biological condensate constituents.

Phase separation-based condensate formation is a novel working paradigm in biology, helping to rationalize many important cellular phenomena including the assembly of membraneless organelles. Uncovering the functional impact of cellular condensates requires a better knowledge of these condensates' constituents. Herein, we introduce the webserver GraPES (Granule Protein Enrichment Server), a user-friendly online interface containing the MaGS and MaGSeq predictors, which provide propensity scores for proteins' localization into cellular condensates. Our webpage contains models trained on human (Homo sapiens) and yeast (Saccharomyces cerevisiae) stress granule proteins. MaGS utilizes experimentally-based protein features for prediction, whereas MaGSeq is an entirely protein sequence-based implementation. GraPES is implemented in HTML/CSS and Javascript and is freely available for public use at https://grapes.msl.ubc.ca/. Documentation for using the provided webtools, descriptions of their methodology, and implementation notes can be found on the webpage.

Protein feature analysis of heat shock induced ubiquitination sites reveals preferential modification site localization.

Protein aggregation is indicative of failing protein quality control systems. These systems are responsible for the refolding or degradation of aberrant and misfolded proteins. Heat stress can cause proteins to misfold, triggering cellular responses including a marked increase in the ubiquitination of proteins. This response has been characterized in yeast, however more studies are needed within mammalian cells. Herein, we examine proteins that become ubiquitinated during heat shock in human tissue culture cells using diGly enrichment coupled with mass spectrometry. A majority of these proteins are localized in the nucleus or cytosol. Proteins which are conjugated under stress display longer sequence lengths, more interaction partners, and more hydrophobic patches than controls but do not show lower melting temperatures. Furthermore, heat-induced conjugation sites occur less frequently in disordered regions and are closer to hydrophobic patches than other ubiquitination sites; perhaps providing novel insight into the molecular mechanism mediating this response. Nuclear and cytosolic pools of modified proteins appear to have different protein features. Using a pulse-SILAC approach, we found that both long-lived and newly-synthesized proteins are conjugated under stress. Modified long-lived proteins are predominately nuclear and were distinct from newly-synthesized proteins, indicating that different pathways may mediate the heat-induced increase of polyubiquitination. SIGNIFICANCE: The maintenance of protein homeostasis requires a balance of protein synthesis, folding, and degradation. Under stress conditions, the cell must rapidly adapt by increasing its folding capacity to eliminate aberrant proteins. A major pathway for proteolysis is mediated by the ubiquitin proteasome system. While increased ubiquitination after heat stress was observed over 30 years ago, it remains unclear which proteins are conjugated during heat shock in mammalian cells and by what means this conjugation occurs. In this study, we combined SILAC-based mass spectrometry with computational analyses to reveal features associated to proteins ubiquitinated while under heat shock. Interestingly, we found that conjugation sites induced by the stress are less often located within disordered regions and more often located near hydrophobic patches. Our study showcases how proteomics can reveal distinct feature associated to a cohort of proteins that are modified post translationally and how the ubiquitin conjugation sites are preferably selected in these conditions. Our work opens a new path for delineating the molecular mechanisms leading to the heat stress response and the regulation of protein homeostasis.

Protein-Protein Interactions Mediated by Intrinsically Disordered Protein Regions Are Enriched in Missense Mutations.

Because proteins are fundamental to most biological processes, many genetic diseases can be traced back to single nucleotide variants (SNVs) that cause changes in protein sequences. However, not all SNVs that result in amino acid substitutions cause disease as each residue is under different structural and functional constraints. Influential studies have shown that protein-protein interaction interfaces are enriched in disease-associated SNVs and depleted in SNVs that are common in the general population. These studies focus primarily on folded (globular) protein domains and overlook the prevalent class of protein interactions mediated by intrinsically disordered regions (IDRs). Therefore, we investigated the enrichment patterns of missense mutation-causing SNVs that are associated with disease and cancer, as well as those present in the healthy population, in structures of IDR-mediated interactions with comparisons to classical globular interactions. When comparing the different categories of interaction interfaces, division of the interface regions into solvent-exposed rim residues and buried core residues reveal distinctive enrichment patterns for the various types of missense mutations. Most notably, we demonstrate a strong enrichment at the interface core of interacting IDRs in disease mutations and its depletion in neutral ones, which supports the view that the disruption of IDR interactions is a mechanism underlying many diseases. Intriguingly, we also found an asymmetry across the IDR interaction interface in the enrichment of certain missense mutation types, which may hint at an increased variant tolerance and urges further investigations of IDR interactions.

Proteomic analysis reveals the direct recruitment of intrinsically disordered regions to stress granules in S. cerevisiae.

Stress granules (SGs) are stress-induced membraneless condensates that store non-translating mRNA and stalled translation initiation complexes. Although metazoan SGs are dynamic compartments where proteins can rapidly exchange with their surroundings, yeast SGs seem largely static. To gain a better understanding of yeast SGs, we identified proteins that sediment after heat shock using mass spectrometry. Proteins that sediment upon heat shock are biased toward a subset of abundant proteins that are significantly enriched in intrinsically disordered regions (IDRs). Heat-induced SG localization of over 80 proteins were confirmed using microscopy, including 32 proteins not previously known to localize to SGs. We found that several IDRs were sufficient to mediate SG recruitment. Moreover, the dynamic exchange of IDRs can be observed using fluorescence recovery after photobleaching, whereas other components remain immobile. Lastly, we showed that the IDR of the Ubp3 deubiquitinase was critical for yeast SG formation. This work shows that IDRs can be sufficient for SG incorporation, can remain dynamic in vitrified SGs, and can play an important role in cellular compartmentalization upon stress.This article has an associated First Person interview with the first author of the paper.

Distinct Features of Stress Granule Proteins Predict Localization in Membraneless Organelles.

Recently generated proteomic data provides unprecedented insight into stress granule composition and stands as fruitful ground for further analysis. Stress granules are stress-induced biological assemblies that are of keen interest due to being linked to both long-term cell viability and a variety of protein aggregation-based diseases. Herein, we compile recently published stress granule composition data, formed specifically through heat and oxidative stress, for both mammalian (Homo sapiens) and yeast (Saccharomyces cerevisiae) cells. Interrogation of the data reveals that stress granule proteins are enriched in features that favor protein liquid-liquid phase separation, being highly disordered, soluble, and abundant while maintaining a high level of protein-protein interactions under basal conditions. Furthermore, these "stress granuleomes" are shown to be enriched for multidomained, RNA-binding proteins with increased potential for post-translational modifications. Findings are consistent with the notion that stress granule formation is driven by protein liquid-liquid phase separation. Furthermore, stress granule proteins appear poised near solubility limits while possessing the ability to dynamically alter their phase behavior in response to external threat. Interestingly, several features, such as protein disorder, are more prominent among stress granule proteins that share homologs between yeast and mammalian systems also found within stress-induced foci. We culminate results from our stress granule analysis into novel predictors for granule incorporation and validate the mammalian predictor's performance against multiple types of membraneless condensates and by colocalization microscopy.

VR-SCOSMO: A smooth conductor-like screening model with charge-dependent radii for modeling chemical reactions.

To better represent the solvation effects observed along reaction pathways, and of ionic species in general, a charge-dependent variable-radii smooth conductor-like screening model (VR-SCOSMO) is developed. This model is implemented and parameterized with a third order density-functional tight binding quantum model, DFTB3/3OB-OPhyd, a quantum method which was developed for organic and biological compounds, utilizing a specific parameterization for phosphate hydrolysis reactions. Unlike most other applications with the DFTB3/3OB model, an auxiliary set of atomic multipoles is constructed from the underlying DFTB3 density matrix which is used to interact the solute with the solvent response surface. The resulting method is variational, produces smooth energies, and has analytic gradients. As a baseline, a conventional SCOSMO model with fixed radii is also parameterized. The SCOSMO and VR-SCOSMO models shown have comparable accuracy in reproducing neutral-molecule absolute solvation free energies; however, the VR-SCOSMO model is shown to reduce the mean unsigned errors (MUEs) of ionic compounds by half (about 2-3 kcal/mol). The VR-SCOSMO model presents similar accuracy as a charge-dependent Poisson-Boltzmann model introduced by Hou et al. [J. Chem. Theory Comput. 6, 2303 (2010)]. VR-SCOSMO is then used to examine the hydrolysis of trimethylphosphate and seven other phosphoryl transesterification reactions with different leaving groups. Two-dimensional energy landscapes are constructed for these reactions and calculated barriers are compared to those obtained from ab initio polarizable continuum calculations and experiment. Results of the VR-SCOSMO model are in good agreement in both cases, capturing the rate-limiting reaction barrier and the nature of the transition state.

Charge-dependent many-body exchange and dispersion interactions in combined QM/MM simulations.

Accurate modeling of the molecular environment is critical in condensed phase simulations of chemical reactions. Conventional quantum mechanical/molecular mechanical (QM/MM) simulations traditionally model non-electrostatic non-bonded interactions through an empirical Lennard-Jones (LJ) potential which, in violation of intuitive chemical principles, is bereft of any explicit coupling to an atom's local electronic structure. This oversight results in a model whereby short-ranged exchange-repulsion and long-ranged dispersion interactions are invariant to changes in the local atomic charge, leading to accuracy limitations for chemical reactions where significant atomic charge transfer can occur along the reaction coordinate. The present work presents a variational, charge-dependent exchange-repulsion and dispersion model, referred to as the charge-dependent exchange and dispersion (QXD) model, for hybrid QM/MM simulations. Analytic expressions for the energy and gradients are provided, as well as a description of the integration of the model into existing QM/MM frameworks, allowing QXD to replace traditional LJ interactions in simulations of reactive condensed phase systems. After initial validation against QM data, the method is demonstrated by capturing the solvation free energies of a series of small, chlorine-containing compounds that have varying charge on the chlorine atom. The model is further tested on the SN2 attack of a chloride anion on methylchloride. Results suggest that the QXD model, unlike the traditional LJ model, is able to simultaneously obtain accurate solvation free energies for a range of compounds while at the same time closely reproducing the experimental reaction free energy barrier. The QXD interaction model allows explicit coupling of atomic charge with many-body exchange and dispersion interactions that are related to atomic size and provides a more accurate and robust representation of non-electrostatic non-bonded QM/MM interactions.

Multiscale methods for computational RNA enzymology.

RNA catalysis is of fundamental importance to biology and yet remains ill-understood due to its complex nature. The multidimensional "problem space" of RNA catalysis includes both local and global conformational rearrangements, changes in the ion atmosphere around nucleic acids and metal ion binding, dependence on potentially correlated protonation states of key residues, and bond breaking/forming in the chemical steps of the reaction. The goal of this chapter is to summarize and apply multiscale modeling methods in an effort to target the different parts of the RNA catalysis problem space while also addressing the limitations and pitfalls of these methods. Classical molecular dynamics simulations, reference interaction site model calculations, constant pH molecular dynamics (CpHMD) simulations, Hamiltonian replica exchange molecular dynamics, and quantum mechanical/molecular mechanical simulations will be discussed in the context of the study of RNA backbone cleavage transesterification. This reaction is catalyzed by both RNA and protein enzymes, and here we examine the different mechanistic strategies taken by the hepatitis delta virus ribozyme and RNase A.

Quantum mechanical study of solvent effects in a prototype SN2 reaction in solution: Cl- attack on CH3Cl.

The nucleophilic attack of a chloride ion on methyl chloride is an important prototype SN2 reaction in organic chemistry that is known to be sensitive to the effects of the surrounding solvent. Herein, we develop a highly accurate Specific Reaction Parameter (SRP) model based on the Austin Model 1 Hamiltonian for chlorine to study the effects of solvation into an aqueous environment on the reaction mechanism. To accomplish this task, we apply high-level quantum mechanical calculations to study the reaction in the gas phase and combined quantum mechanical/molecular mechanical simulations with TIP3P and TIP4P-ew water models and the resulting free energy profiles are compared with those determined from simulations using other fast semi-empirical quantum models. Both gas phase and solution results with the SRP model agree very well with experiment and provide insight into the specific role of solvent on the reaction coordinate. Overall, the newly parameterized SRP Hamiltonian is able to reproduce both the gas phase and solution phase barriers, suggesting it is an accurate and robust model for simulations in the aqueous phase at greatly reduced computational cost relative to comparably accurate ab initio and density functional models.

A variational linear-scaling framework to build practical, efficient next-generation orbital-based quantum force fields.

We introduce a new hybrid molecular orbital/density-functional modified divide-and-conquer (mDC) approach that allows the linear-scaling calculation of very large quantum systems. The method provides a powerful framework from which linear-scaling force fields for molecular simulations can be developed. The method is variational in the energy, and has simple, analytic gradients and essentially no break-even point with respect to the corresponding full electronic structure calculation. Furthermore, the new approach allows intermolecular forces to be properly balanced such that non-bonded interactions can be treated, in some cases, to much higher accuracy than the full calculation. The approach is illustrated using the second-order self-consistent charge density-functional tight-binding model (DFTB2). Using this model as a base Hamiltonian, the new mDC approach is applied to a series of water systems, where results show that geometries and interaction energies between water molecules are greatly improved relative to full DFTB2. In order to achieve substantial improvement in the accuracy of intermolecular binding energies and hydrogen bonded cluster geometries, it was necessary to extend the DFTB2 model to higher-order atom-centered multipoles for the second-order self-consistent intermolecular electrostatic term. Using generalized, linear-scaling electrostatic methods, timings demonstrate that the method is able to calculate a water system of 3000 atoms in less than half of a second, and systems of up to one million atoms in only a few minutes using a conventional desktop workstation.