Protein Abundance Biases the Amino Acid Composition of Disordered Regions to Minimize Non-functional Interactions

Domain architecture of the Mycobacterium tuberculosis MabR (Rv2242), a member of the PucR transcription factor family

MabR (Rv2242), a PucR-type transcription factor, plays a crucial role in regulating mycolic acid biosynthesis in Mycobacterium tuberculosis. To understand its regulatory mechanisms, we determined the crystal structures of its N-terminal and C-terminal domains. The N-terminal domain adopts a globin-like fold, while the C-terminal domain comprises an α/β GGDEF domain and an all-α effector domain with a helix-turn-helix DNA-binding motif. This unique domain combination is specific to Actinomycetes. Biochemical and computational studies suggest that full-length MabR forms both dimeric and tetrameric assemblies in solution. Structural analysis revealed two distinct dimerization interfaces within the N- and C-terminal domains, further supporting a tetrameric organization. These findings provide valuable insights into the domain architecture, oligomeric state, and potential regulatory mechanisms of MabR.

The kinase NEK6 positively regulates LSD1 activity and accumulation in local chromatin sub-compartments

LSD1 plays a crucial role in mammalian biology, regulated through interactions with coregulators and post-translational modifications. Here we show that the kinase NEK6 stimulates LSD1 activity in cells and observe a strong colocalization of NEK6 and LSD1 at distinct chromatin sub-compartments (CSCs). We demonstrate that LSD1 is a substrate for NEK6 phosphorylation at the N-terminal intrinsically disordered region (IDR) of LSD1, which shows phase separation behavior in vitro and in cells. The LSD1-IDR is important for LSD1 activity and functions to co-compartmentalize NEK6, histone peptides and DNA. The subsequent phosphorylation of LSD1 by NEK6 supports the concentration of LSD1 at these distinct CSCs, which is imperative for dynamic control of transcription. This suggest that phase separation is crucial for the regulatory function of LSD1 and our findings highlight the role of NEK6 in modulating LSD1 activity and phase separation, expanding our understanding of LSD1 regulation and its implications in cellular processes.

SEC-SAXS/MC Ensemble Structural Studies of the Microtubule Binding Protein Cdt1 Show Monomeric, Folded-Over Conformations

Cdt1 is a mixed folded protein critical for DNA replication licensing and it also has a "moonlighting" role at the kinetochore via direct binding to microtubules and the Ndc80 complex. However, it is unknown how the structure and conformations of Cdt1 could allow it to participate in these multiple, unique sets of protein complexes. While robust methods exist to study entirely folded or unfolded proteins, structure-function studies of combined, mixed folded/disordered proteins remain challenging. In this work, we employ orthogonal biophysical and computational techniques to provide structural characterization of mitosis-competent human Cdt1. Thermal stability analyses shows that both folded winged helix domains1 are unstable. CD and NMR show that the N-terminal and linker regions are intrinsically disordered. DLS shows that Cdt1 is monomeric and polydisperse, while SEC-MALS confirms that it is monomeric at high concentrations, but without any apparent inter-molecular self-association. SEC-SAXS enabled computational modeling of the protein structures. Using the program SASSIE, we performed rigid body Monte Carlo simulations to generate a conformational ensemble of structures. We observe that neither fully extended nor extremely compact Cdt1 conformations are consistent with SAXS. The best-fit models have the N-terminal and linker disordered regions extended into the solution and the two folded domains close to each other in apparent "folded over" conformations. We hypothesize the best-fit Cdt1 conformations could be consistent with a function as a scaffold protein that may be sterically blocked without binding partners. Our study also provides a template for combining experimental and computational techniques to study mixed-folded proteins.

The substrate quality of CK2 target sites has a determinant role on their function and evolution

Most biological processes are regulated by signaling modules that bind to short linear motifs. For protein kinases, substrates may have full or only partial matches to the kinase recognition motif, a property known as "substrate quality." However, it is not clear whether differences in substrate quality represent neutral variation or if they have functional consequences. We examine this question for the kinase CK2, which has many fundamental functions. We show that optimal CK2 sites are phosphorylated at maximal stoichiometries and found in many conditions, whereas minimal substrates are more weakly phosphorylated and have regulatory functions. Optimal CK2 sites tend to be more conserved, and substrate quality is often tuned by selection. For intermediate sites, increases or decreases in substrate quality may be deleterious, as we demonstrate for a CK2 substrate at the kinetochore. The results together suggest a strong role for substrate quality in phosphosite function and evolution. A record of this paper's transparent peer review process is included in the supplemental information.

SAXS/MC studies of the mixed-folded protein Cdt1 reveal monomeric, folded over conformations

Cdt1 is a protein critical for DNA replication licensing and is well-established to be a binding partner of the minichromosome maintenance (MCM) complex. Cdt1 has also been demonstrated to have an emerging, "moonlighting" role at the kinetochore via direct binding to microtubules and to the Ndc80 complex. However, it is not known how the structure and conformations of Cdt1 could allow for these multiple, completely unique sets of protein complexes. And while there exist multiple robust methods to study entirely folded or entirely unfolded proteins, structure-function studies of combined, mixed folded/disordered proteins remain challenging. It this work, we employ multiple orthogonal biophysical and computational techniques to provide a detailed structural characterization of human Cdt1 92-546. DSF and DSCD show both folded winged helix (WH) domains of Cdt1 are relatively unstable. CD and NMR show the N-terminal and the linker regions are intrinsically disordered. Using DLS and SEC-MALS, we show that Cdt1 is polydisperse, monomeric at high concentrations, and without any apparent inter-molecular self-association. SEC-SAXS of the monomer in solution enabled computational modeling of the protein in silico. Using the program SASSIE, we performed rigid body Monte Carlo simulations to generate a conformational ensemble. Using experimental SAXS data, we filtered for conformations which did and did not fit our data. We observe that neither fully extended nor extremely compact Cdt1 conformations are consistent with our SAXS data. The best fit models have the N-terminal and linker regions extended into solution and the two folded domains close to each other in apparent "folded over" conformations. The best fit Cdt1 conformations are consistent with a function as a scaffold protein which may be sterically blocked without the presence of binding partners. Our studies also provide a template for combining experimental and computational biophysical techniques to study mixed-folded proteins.

Photo-caged 2-butene-1,4-dial as an efficient, target-specific photo-crosslinker for covalent trapping of DNA-binding proteins

Covalent trapping of DNA-binding proteins via photo-crosslinking is an advantageous method for studying DNA-protein interactions. However, traditional photo-crosslinkers generate highly reactive intermediates that rapidly and non-selectively react with nearby functional groups, resulting in low target-capture yields and high non-target background capture. Herein, we report that photo-caged 2-butene-1,4-dial (PBDA) is an efficient photo-crosslinker for trapping DNA-binding proteins. Photo-irradiation (360 nm) of PBDA-modified DNA generates 2-butene-1,4-dial (BDA), a small, long-lived intermediate that reacts selectively with Lys residues of DNA-binding proteins, leading in minutes to stable DNA-protein crosslinks in up to 70% yield. In addition, BDA exhibits high specificity for target proteins, leading to low non-target background capture. The high photo-crosslinking yield and target specificity make PBDA a powerful tool for studying DNA-protein interactions.

"Order from disordered": Potential role of intrinsically disordered regions in phytopathogenic oomycete intracellular effector proteins

There is a continuous arms race between pathogens and their host plants. However, successful pathogens, such as phytopathogenic oomycetes, secrete effector proteins to manipulate host defense responses for disease development. Structural analyses of these effector proteins reveal the existence of regions that fail to fold into three-dimensional structures, intrinsically disordered regions (IDRs). Because of their flexibility, these regions are involved in important biological functions of effector proteins, such as effector-host protein interactions that perturb host immune responses. Despite their significance, the role of IDRs in phytopathogenic oomycete effector-host protein interactions is not clear. This review, therefore, searched the literature for functionally characterized oomycete intracellular effectors with known host interactors. We further classify regions that mediate effector-host protein interactions into globular or disordered binding sites in these proteins. To fully appreciate the potential role of IDRs, five effector proteins encoding potential disordered binding sites were used as case studies. We also propose a pipeline that can be used to identify, classify as well as characterize potential binding regions in effector proteins. Understanding the role of IDRs in these effector proteins can aid in the development of new disease-control strategies.

mt-tRNAs in the polymerase gamma mutant heart

Introduction

Mice harboring a D257A mutation in the proofreading domain of the mitochondrial DNA polymerase, Polymerase Gamma (POLG), experience severe metabolic dysfunction and display hallmarks of accelerated aging. We previously reported a mitochondrial unfolded protein response (UPTmt) - like (UPRmt-like) gene and protein expression pattern in the right ventricular tissue of POLG mutant mice.

Aim

We sought to determine if POLG mutation altered the expression of genes encoded by the mitochondria in a way that might also reduce proteotoxic stress.

Methods and Results

The expression of genes encoded by the mitochondrial DNA was interrogated via RNA-seq and northern blot analysis. A striking, location-dependent effect was seen in the expression of mitochondrial-encoded tRNAs in the POLG mutant as assayed by RNA-seq. These expression changes were negatively correlated with the tRNA partner amino acid's amyloidogenic potential. Direct measurement by northern blot was conducted on candidate mt-tRNAs identified from the RNA-seq. This analysis confirmed reduced expression of MT-TY in the POLG mutant but failed to show increased expression of MT-TP, which was dramatically increased in the RNA-seq data.

Conclusion

We conclude that reduced expression of amyloid-associated mt-tRNAs is another indication of adaptive response to severe mitochondrial dysfunction in the POLG mutant. Incongruence between RNA-seq and northern blot measurement of MT-TP expression points towards the existence of mt-tRNA post-transcriptional modification regulation in the POLG mutant that alters either polyA capture or cDNA synthesis in RNA-seq library generation. Together, these data suggest that 1) evolution has distributed mt-tRNAs across the circular mitochondrial genome to allow chromosomal location-dependent mt-tRNA regulation (either by expression or PTM) and 2) this regulation is cognizant of the tRNA partner amino acid's amyloidogenic properties.

The mRNA binding-mediated self-regulatory function of small heat shock protein IbpA in γ-proteobacteria is conferred by a conserved arginine

Bacterial small heat shock proteins, such as inclusion body-associated protein A (IbpA) and IbpB, coaggregate with denatured proteins and recruit other chaperones for the processing of aggregates thereby assisting in protein refolding. In addition, as a recently revealed uncommon feature, Escherichia coli IbpA self-represses its own translation through interaction with the 5'-untranslated region of the ibpA mRNA, enabling IbpA to act as a mediator of negative feedback regulation. Although IbpA also suppresses the expression of IbpB, IbpB does not have this self-repression activity despite the two Ibps being highly homologous. In this study, we demonstrate that the self-repression function of IbpA is conserved in other γ-proteobacterial IbpAs. Moreover, we show a cationic residue-rich region in the α-crystallin domain of IbpA, which is not conserved in IbpB, is critical for the self-suppression activity. Notably, we found arginine 93 (R93) located within the α-crystallin domain is an essential residue that cannot be replaced by any of the other 19 amino acids including lysine. We observed that IbpA-R93 mutants completely lost the interaction with the 5' untranslated region of the ibpA mRNA, but retained almost all chaperone activity and were able to sequester denatured proteins. Taken together, we propose the conserved Arg93-mediated translational control of IbpA through RNA binding would be beneficial for a rapid and massive supply of the chaperone on demand.

Proteins released into the plant apoplast by the obligate parasitic protist Albugo selectively repress phyllosphere-associated bacteria

Biotic and abiotic interactions shape natural microbial communities. The mechanisms behind microbe-microbe interactions, particularly those protein based, are not well understood. We hypothesize that released proteins with antimicrobial activity are a powerful and highly specific toolset to shape and defend plant niches. We have studied Albugo candida, an obligate plant parasite from the protist Oomycota phylum, for its potential to modulate the growth of bacteria through release of antimicrobial proteins into the apoplast. Amplicon sequencing and network analysis of Albugo-infected and uninfected wild Arabidopsis thaliana samples revealed an abundance of negative correlations between Albugo and other phyllosphere microbes. Analysis of the apoplastic proteome of Albugo-colonized leaves combined with machine learning predictors enabled the selection of antimicrobial candidates for heterologous expression and study of their inhibitory function. We found for three candidate proteins selective antimicrobial activity against Gram-positive bacteria isolated from A. thaliana and demonstrate that these inhibited bacteria are precisely important for the stability of the community structure. We could ascribe the antibacterial activity of the candidates to intrinsically disordered regions and positively correlate it with their net charge. This is the first report of protist proteins with antimicrobial activity under apoplastic conditions that therefore are potential biocontrol tools for targeted manipulation of the microbiome.

Protein compactness and interaction valency define the architecture of a biomolecular condensate across scales

Non-membrane-bound biomolecular condensates have been proposed to represent an important mode of subcellular organization in diverse biological settings. However, the fundamental principles governing the spatial organization and dynamics of condensates at the atomistic level remain unclear. The Saccharomyces cerevisiae Lge1 protein is required for histone H2B ubiquitination and its N-terminal intrinsically disordered fragment (Lge11-80) undergoes robust phase separation. This study connects single- and multi-chain all-atom molecular dynamics simulations of Lge11-80 with the in vitro behavior of Lge11-80 condensates. Analysis of modeled protein-protein interactions elucidates the key determinants of Lge11-80 condensate formation and links configurational entropy, valency, and compactness of proteins inside the condensates. A newly derived analytical formalism, related to colloid fractal cluster formation, describes condensate architecture across length scales as a function of protein valency and compactness. In particular, the formalism provides an atomistically resolved model of Lge11-80 condensates on the scale of hundreds of nanometers starting from individual protein conformers captured in simulations. The simulation-derived fractal dimensions of condensates of Lge11-80 and its mutants agree with their in vitro morphologies. The presented framework enables a multiscale description of biomolecular condensates and embeds their study in a wider context of colloid self-organization.

Buffering of genetic dominance by allele-specific protein complex assembly

Protein complex assembly often occurs while subunits are being translated, resulting in complexes whose subunits were translated from the same mRNA in an allele-specific manner. It has thus been hypothesized that such cotranslational assembly may counter the assembly-mediated dominant-negative effect, whereby co-assembly of mutant and wild-type subunits "poisons" complex activity. Here, we show that cotranslationally assembling subunits are much less likely to be associated with autosomal dominant relative to recessive disorders, and that subunits with dominant-negative disease mutations are significantly depleted in cotranslational assembly compared to those associated with loss-of-function mutations. We also find that complexes with known dominant-negative effects tend to expose their interfaces late during translation, lessening the likelihood of cotranslational assembly. Finally, by combining complex properties with other features, we trained a computational model for predicting proteins likely to be associated with non-loss-of-function disease mechanisms, which we believe will be of considerable utility for protein variant interpretation.

Prediction of Quality-control Degradation Signals in Yeast Proteins

Effective proteome homeostasis is key to cellular and organismal survival, and cells therefore contain efficient quality control systems to monitor and remove potentially toxic misfolded proteins. Such general protein quality control to a large extent relies on the efficient and robust delivery of misfolded or unfolded proteins to the ubiquitin-proteasome system. This is achieved via recognition of so-called degradation motifs-degrons-that are assumed to become exposed as a result of protein misfolding. Despite their importance, the nature and sequence properties of quality-control degrons remain elusive. Here, we have used data from a yeast-based screen of 23,600 17-residue peptides to build a predictor of quality-control degrons. The resulting model, QCDPred (Quality Control Degron Prediction), achieves good accuracy using only the sequence composition of the peptides as input. Our analysis reveals that strong degrons are enriched in hydrophobic amino acids and depleted in negatively charged amino acids, in line with the expectation that they are buried in natively folded proteins. We applied QCDPred to the yeast proteome, enabling us to analyse more widely the potential effects of degrons. As an example, we show a correlation between cellular abundance and degron potential in disordered regions of proteins. Together with recent results on membrane proteins, our work suggest that the recognition of exposed hydrophobic residues is a key and generic mechanism for proteome homeostasis. QCDPred is freely available as open source code and via a web interface.

HSP70-binding motifs function as protein quality control degrons

Protein quality control (PQC) degrons are short protein segments that target misfolded proteins for proteasomal degradation, and thus protect cells against the accumulation of potentially toxic non-native proteins. Studies have shown that PQC degrons are hydrophobic and rarely contain negatively charged residues, features which are shared with chaperone-binding regions. Here we explore the notion that chaperone-binding regions may function as PQC degrons. When directly tested, we found that a canonical Hsp70-binding motif (the APPY peptide) functioned as a dose-dependent PQC degron both in yeast and in human cells. In yeast, Hsp70, Hsp110, Fes1, and the E3 Ubr1 target the APPY degron. Screening revealed that the sequence space within the chaperone-binding region of APPY that is compatible with degron function is vast. We find that the number of exposed Hsp70-binding sites in the yeast proteome correlates with a reduced protein abundance and half-life. Our results suggest that when protein folding fails, chaperone-binding sites may operate as PQC degrons, and that the sequence properties leading to PQC-linked degradation therefore overlap with those of chaperone binding.

Nonspecific Binding of Adenosine Tripolyphosphate and Tripolyphosphate Modulates the Phase Behavior of Lysozyme

Adenosine tripolyphosphate (ATP) is a small polyvalent anion that has recently been shown to interact with proteins and have a major impact on assembly processes involved in biomolecular condensate formation and protein aggregation. However, the nature of non-specific protein-ATP interactions and their effects on protein solubility are largely unknown. Here, the binding of ATP to the globular model protein is characterized in detail using X-ray crystallography and nuclear magnetic resonance (NMR). Using NMR, we identified six ATP binding sites on the lysozyme surface, with one known high-affinity nucleic acid binding site and five non-specific previously unknown sites with millimolar affinities that also bind tripolyphosphate (TPP). ATP binding occurs primarily through the polyphosphate moiety, which was confirmed by the X-ray structure of the lysozyme-ATP complex. Importantly, ATP binds preferentially to arginine over lysine in non-specific binding sites. ATP and TPP have similar effects on solution-phase protein-protein interactions. At low salt concentrations, ion binding to lysozyme causes precipitation, while at higher salt concentrations, redissolution occurs. The addition of an equimolar concentration of magnesium to ATP does not alter ATP binding affinities but prevents lysozyme precipitation. These findings have important implications for both protein crystallization and cell biology. Crystallization occurs readily in ATP solutions outside the well-established crystallization window. In the context of cell biology, the findings suggest that ATP binds non-specifically to folded proteins in physiological conditions. Based on the nature of the binding sites identified by NMR, we propose several mechanisms for how ATP binding can prevent the aggregation of natively folded proteins.

Correlation between protein abundance and sequence conservation: what do recent experiments say?

Cells evolve in a space of parameter values set by physical and chemical forces. These constraints create associations among cellular properties. A particularly strong association is the negative correlation between the rate of evolution of proteins and their abundance in the cell. Highly expressed proteins evolve slower than lowly expressed ones. Multiple hypotheses have been put forward to explain this relationship, including, for instance, the requirement for higher mRNA stability, misfolding avoidance, and misinteraction avoidance for highly expressed proteins. Here, we review some of these hypotheses, their predictions, and how they are supported to finally discuss recent experiments that have been performed to test these predictions.

Prediction of protein-protein interaction sites in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) participate in many biological processes by interacting with other proteins, including the regulation of transcription, translation, and the cell cycle. With the increasing amount of disorder sequence data available, it is thus crucial to identify the IDP binding sites for functional annotation of these proteins. Over the decades, many computational approaches have been developed to predict protein-protein binding sites of IDP (IDP-PPIS) based on protein sequence information. Moreover, there are new IDP-PPIS predictors developed every year with the rapid development of artificial intelligence. It is thus necessary to provide an up-to-date overview of these methods in this field. In this paper, we collected 30 representative predictors published recently and summarized the databases, features and algorithms. We described the procedure how the features were generated based on public data and used for the prediction of IDP-PPIS, along with the methods to generate the feature representations. All the predictors were divided into three categories: scoring functions, machine learning-based prediction, and consensus approaches. For each category, we described the details of algorithms and their performances. Hopefully, our manuscript will not only provide a full picture of the status quo of IDP binding prediction, but also a guide for selecting different methods. More importantly, it will shed light on the inspirations for future development trends and principles.

Physics-driven coarse-grained model for biomolecular phase separation with near-quantitative accuracy

Various physics- and data-driven sequence-dependent protein coarse-grained models have been developed to study biomolecular phase separation and elucidate the dominant physicochemical driving forces. Here, we present Mpipi, a multiscale coarse-grained model that describes almost quantitatively the change in protein critical temperatures as a function of amino-acid sequence. The model is parameterised from both atomistic simulations and bioinformatics data and accounts for the dominant role of π-π and hybrid cation-π/π-π interactions and the much stronger attractive contacts established by arginines than lysines. We provide a comprehensive set of benchmarks for Mpipi and seven other residue-level coarse-grained models against experimental radii of gyration and quantitative in-vitro phase diagrams; Mpipi predictions agree well with experiment on both fronts. Moreover, it can account for protein-RNA interactions, correctly predicts the multiphase behaviour of a charge-matched poly-arginine/poly-lysine/RNA system, and recapitulates experimental LLPS trends for sequence mutations on FUS, DDX4 and LAF-1 proteins.

Altered Protein Abundance and Localization Inferred from Sites of Alternative Modification by Ubiquitin and SUMO

Protein modification by ubiquitin or SUMO can alter the function, stability or activity of target proteins. Previous studies have identified thousands of substrates that were modified by ubiquitin or SUMO on the same lysine residue. However, it remains unclear whether such overlap could result from a mere higher solvent accessibility, whether proteins containing those sites are associated with specific functional traits, and whether selectively perturbing their modification by ubiquitin or SUMO could result in different phenotypic outcomes. Here, we mapped reported lysine modification sites across the human proteome and found an enrichment of sites reported to be modified by both ubiquitin and SUMO. Our analysis uncovered thousands of proteins containing such sites, which we term Sites of Alternative Modification (SAMs). Among more than 36,000 sites reported to be modified by SUMO, 51.8% have also been reported to be modified by ubiquitin. SAM-containing proteins are associated with diverse biological functions including cell cycle, DNA damage, and transcriptional regulation. As such, our analysis highlights numerous proteins and pathways as putative targets for further elucidating the crosstalk between ubiquitin and SUMO. Comparing the biological and biochemical properties of SAMs versus other non-overlapping modification sites revealed that these sites were associated with altered cellular localization or abundance of their host proteins. Lastly, using S. cerevisiae as model, we show that mutating the SAM motif in a protein can influence its ubiquitination as well as its localization and abundance.

The Distinct Properties of the Consecutive Disordered Regions Inside or Outside Protein Domains and Their Functional Significance

The consecutive disordered regions (CDRs) are the basis for the formation of intrinsically disordered proteins, which contribute to various biological functions and increasing organism complexity. Previous studies have revealed that CDRs may be present inside or outside protein domains, but a comprehensive analysis of the property differences between these two types of CDRs and the proteins containing them is lacking. In this study, we investigated this issue from three viewpoints. Firstly, we found that in-domain CDRs are more hydrophilic and stable but have less stickiness and fewer post-translational modification sites compared with out-domain CDRs. Secondly, at the protein level, we found that proteins with only in-domain CDRs originated late, evolved rapidly, and had weak functional constraints, compared with the other two types of CDR-containing proteins. Proteins with only in-domain CDRs tend to be expressed spatiotemporal specifically, but they tend to have higher abundance and are more stable. Thirdly, we screened the CDR-containing protein domains that have a strong correlation with organism complexity. The CDR-containing domains tend to be evolutionarily young, or they changed from a domain without CDR to a CDR-containing domain during evolution. These results provide valuable new insights about the evolution and function of CDRs and protein domains.

On the specificity of protein-protein interactions in the context of disorder

With the increased focus on intrinsically disordered proteins (IDPs) and their large interactomes, the question about their specificity — or more so on their multispecificity — arise. Here we recapitulate how specificity and multispecificity are quantified and address through examples if IDPs in this respect differ from globular proteins. The conclusion is that quantitatively, globular proteins and IDPs are similar when it comes to specificity. However, compared with globular proteins, IDPs have larger interactome sizes, a phenomenon that is further enabled by their flexibility, repetitive binding motifs and propensity to adapt to different binding partners. For IDPs, this adaptability, interactome size and a higher degree of multivalency opens for new interaction mechanisms such as facilitated exchange through trimer formation and ultra-sensitivity via threshold effects and ensemble redistribution. IDPs and their interactions, thus, do not compromise the definition of specificity. Instead, it is the sheer size of their interactomes that complicates its calculation. More importantly, it is this size that challenges how we conceptually envision, interpret and speak about their specificity.

Abundance Imparts Evolutionary Constraints of Similar Magnitude on the Buried, Surface, and Disordered Regions of Proteins

An understanding of the forces shaping protein conservation is key, both for the fundamental knowledge it represents and to allow for optimal use of evolutionary information in practical applications. Sequence conservation is typically examined at one of two levels. The first is a residue-level, where intra-protein differences are analyzed and the second is a protein-level, where inter-protein differences are studied. At a residue level, we know that solvent-accessibility is a prime determinant of conservation. By inverting this logic, we inferred that disordered regions are slightly more solvent-accessible on average than the most exposed surface residues in domains. By integrating abundance information with evolutionary data within and across proteins, we confirmed a previously reported strong surface-core association in the evolution of structured regions, but we found a comparatively weak association between disordered and structured regions. The facts that disordered and structured regions experience different structural constraints and evolve independently provide a unique setup to examine an outstanding question: why is a protein’s abundance the main determinant of its sequence conservation? Indeed, any structural or biophysical property linked to the abundance-conservation relationship should increase the relative conservation of regions concerned with that property (e.g., disordered residues with mis-interactions, domain residues with misfolding). Surprisingly, however, we found the conservation of disordered and structured regions to increase in equal proportion with abundance. This observation implies that either abundance-related constraints are structure-independent, or multiple constraints apply to different regions and perfectly balance each other.

Integration of Data from Liquid-Liquid Phase Separation Databases Highlights Concentration and Dosage Sensitivity of LLPS Drivers

Liquid–liquid phase separation (LLPS) is a molecular process that leads to the formation of membraneless organelles, representing functionally specialized liquid-like cellular condensates formed by proteins and nucleic acids. Integrating the data on LLPS-associated proteins from dedicated databases revealed only modest agreement between them and yielded a high-confidence dataset of 89 human LLPS drivers. Analysis of the supporting evidence for our dataset uncovered a systematic and potentially concerning difference between protein concentrations used in a good fraction of the in vitro LLPS experiments, a key parameter that governs the phase behavior, and the proteomics-derived cellular abundance levels of the corresponding proteins. Closer scrutiny of the underlying experimental data enabled us to offer a sound rationale for this systematic difference, which draws on our current understanding of the cellular organization of the proteome and the LLPS process. In support of this rationale, we find that genes coding for our human LLPS drivers tend to be dosage-sensitive, suggesting that their cellular availability is tightly regulated to preserve their functional role in direct or indirect relation to condensate formation. Our analysis offers guideposts for increasing agreement between in vitro and in vivo studies, probing the roles of proteins in LLPS.

Universal and taxon-specific trends in protein sequences as a function of age

Extant protein-coding sequences span a huge range of ages, from those that emerged only recently to those present in the last universal common ancestor. Because evolution has had less time to act on young sequences, there might be 'phylostratigraphy' trends in any properties that evolve slowly with age. A long-term reduction in hydrophobicity and hydrophobic clustering was found in previous, taxonomically restricted studies. Here we perform integrated phylostratigraphy across 435 fully sequenced species, using sensitive HMM methods to detect protein domain homology. We find that the reduction in hydrophobic clustering is universal across lineages. However, only young animal domains have a tendency to have higher structural disorder. Among ancient domains, trends in amino acid composition reflect the order of recruitment into the genetic code, suggesting that the composition of the contemporary descendants of ancient sequences reflects amino acid availability during the earliest stages of life, when these sequences first emerged.