IUPred2A: context-dependent prediction of protein disorder as a function of redox state and protein binding

MMPatho: Leveraging Multilevel Consensus and Evolutionary Information for Enhanced Missense Mutation Pathogenic Prediction

Understanding the pathogenicity of missense mutation (MM) is essential for shed light on genetic diseases, gene functions, and individual variations. In this study, we propose a novel computational approach, called MMPatho, for enhancing missense mutation pathogenic prediction. First, we established a large-scale nonredundant MM benchmark data set based on the entire Ensembl database, complemented by a focused blind test set specifically for pathogenic GOF/LOF MM. Based on this data set, for each mutation, we utilized Ensembl VEP v104 and dbNSFP v4.1a to extract variant-level, amino acid-level, individuals' outputs, and genome-level features. Additionally, protein sequences were generated using ENSP identifiers with the Ensembl API, and then encoded. The mutant sites' ESM-1b and ProtTrans-T5 embeddings were subsequently extracted. Then, our model group (MMPatho) was developed by leveraging upon these efforts, which comprised ConsMM and EvoIndMM. To be specific, ConsMM employs individuals' outputs and XGBoost with SHAP explanation analysis, while EvoIndMM investigates the potential enhancement of predictive capability by incorporating evolutionary information from ESM-1b and ProtT5-XL-U50, large protein language embeddings. Through rigorous comparative experiments, both ConsMM and EvoIndMM were capable of achieving remarkable AUROC (0.9836 and 0.9854) and AUPR (0.9852 and 0.9902) values on the blind test set devoid of overlapping variations and proteins from the training data, thus highlighting the superiority of our computational approach in the prediction of MM pathogenicity. Our Web server, available at http://csbio.njust.edu.cn/bioinf/mmpatho/, allows researchers to predict the pathogenicity (alongside the reliability index score) of MMs using the ConsMM and EvoIndMM models and provides extensive annotations for user input. Additionally, the newly constructed benchmark data set and blind test set can be accessed via the data page of our web server.

SARS-CoV-2 Orphan Gene ORF10 Contributes to More Severe COVID-19 Disease

The orphan gene of SARS-CoV-2, ORF10, is the least studied gene in the virus responsible for the COVID-19 pandemic. Recent experimentation indicated ORF10 expression moderates innate immunity in vitro. However, whether ORF10 affects COVID-19 in humans remained unknown. We determine that the ORF10 sequence is identical to the Wuhan-Hu-1 ancestral haplotype in 95% of genomes across five variants of concern (VOC). Four ORF10 variants are associated with less virulent clinical outcomes in the human host: three of these affect ORF10 protein structure, one affects ORF10 RNA structural dynamics. RNA-Seq data from 2070 samples from diverse human cells and tissues reveals ORF10 accumulation is conditionally discordant from that of other SARS-CoV-2 transcripts. Expression of ORF10 in A549 and HEK293 cells perturbs immune-related gene expression networks, alters expression of the majority of mitochondrially-encoded genes of oxidative respiration, and leads to large shifts in levels of 14 newly-identified transcripts. We conclude ORF10 contributes to more severe COVID-19 clinical outcomes in the human host.

Time-resolved proteomic profiling reveals compositional and functional transitions across the stress granule life cycle

Stress granules (SGs) are dynamic, membrane-less organelles. With their formation and disassembly processes characterized, it remains elusive how compositional transitions are coordinated during prolonged stress to meet changing functional needs. Here, using time-resolved proteomic profiling of the acute to prolonged heat-shock SG life cycle, we identify dynamic SG proteins, further segregated into early and late proteins. Comparison of different groups of SG proteins suggests that their biochemical properties help coordinate SG compositional and functional transitions. In particular, early proteins, with high phase-separation-propensity, drive the rapid formation of the initial SG platform, while late proteins are subsequently recruited as discrete modules to further functionalize SGs. This model, supported by immunoblotting and immunofluorescence imaging, provides a conceptual framework for the compositional transitions throughout the acute to prolonged SG life cycle. Additionally, an early SG constituent, non-muscle myosin II, is shown to promote SG formation by increasing SG fusion, underscoring the strength of this dataset in revealing the complexity of SG regulation.

Augmenting L3MBTL2-induced condensates suppresses tumor growth in osteosarcoma

Osteosarcoma is a highly aggressive cancer and lacks effective therapeutic targets. We found that L3MBTL2 acts as a tumor suppressor by transcriptionally repressing IFIT2 in osteosarcoma. L3MBTL2 recruits the components of Polycomb repressive complex 1.6 to form condensates via both Pho-binding pockets and polybasic regions within carboxyl-terminal intrinsically disordered regions; the L3MBTL2-induced condensates are required for its tumor suppression. Multi-monoubiquitination of L3MBTL2 by UBE2O results in its proteasomal degradation, and the UBE2O/L3MBTL2 axis was crucial for osteosarcoma growth. There is a reverse correlation between L3MBTL2 and UBE2O in osteosarcoma tissues, and higher UBE2O and lower L3MBTL2 are associated with poorer prognosis in osteosarcoma. Pharmacological blockage of UBE2O by arsenic trioxide can enhance L3MBTL2-induced condensates and consequently suppress osteosarcoma growth. Our findings unveil a crucial biological function of L3MBTL2-induced condensates in mediating tumor suppression, proposing the UBE2O-L3MBTL2 axis as a potential cancer therapeutic target in osteosarcoma.

Molecular insights into the VIRESCENS amino acid sequence and its implication in anthocyanin production in red- and yellow-fruited cultivars of date palm

This study examined the amino acid sequence of the VIRESCENS gene (VIR), which regulates the production of anthocyanin in 12 cultivars of the date palm (Phoenix dactylifera L.), grown in Al-Madinah Al-Munawarah of the Kingdom of Saudi Arabia. The gene products were amplified via polymerase chain reactions, amplifying both exons and introns. The products were sequenced for the reconstruction of a phylogenetic tree, which used the associated amino acid sequences. The ripening stages of Khalal, Rutab, and Tamar varied among the cultivars. Regarding VIR genotype, the red date had the wild-type gene (VIR+), while the yellow date carried a dominant mutation (VIRIM), i.e., long terminal repeat retrotransposons (LTR-RTs). The DNA sequence of VIRIM revealed that the insertion length of the LTR-RTs ranged between 386 and 476 bp. The R2 and R3 motifs in both VIR+ and VIRIM were conserved. The C-terminus motifs S6A, S6B, and S6C were found in the VIR+ protein sequence. However, the amino acids at positions 123, 161, 166, and 168 differed between VIR+ and VIRIM, and were not included in the C-terminus motifs. Within the VIR+ allele, the lysine at position 187 in the C-terminus was located immediately after S6B, with a protein binding score of 0.3, which was unique to the dark, red-fruited cultivars Ajwah, Anbarah, and Safawi. In the lighter, red-fruited cultivars, the presence of glutamic acid at the same position suggested that the anthocyanin regulation of date palm might be outside the R2 and R3 domains in the N-terminus.

Prediction of folding patterns for intrinsic disordered protein

The conformation flexibility of natural protein causes both complexity and difficulty to understand the relationship between structure and function. The prediction of intrinsically disordered protein primarily is focusing on to disclose the regions with structural flexibility involving relevant biological functions and various diseases. The order of amino acids in protein sequence determines possible conformations, folding flexibility and biological function. Although many methods provided the information of intrinsically disordered protein (IDP), but the results are mainly limited to determine the locations of regions without knowledge of possible folding conformations. Here, the developed protein folding fingerprint adopted the protein folding variation matrix (PFVM) to reveal all possible folding patterns for the intrinsically disordered protein along its sequence. The PFVM integrally exhibited the intrinsically disordered protein with disordering regions, degree of disorder as well as folding pattern. The advantage of PFVM will not only provide rich information for IDP, but also may promote the study of protein folding problem.

Functional Consequences of Shifting Transcript Boundaries in Glucose Starvation

Glucose is a major source of carbon and essential for the survival of many organisms, ranging from yeast to human. A sudden 60-fold reduction of glucose in exponentially growing fission yeast induces transcriptome-wide changes in gene expression. This regulation is multilayered, and the boundaries of transcripts are known to vary, with functional consequences at the protein level. By combining direct RNA sequencing with 5'-CAGE and short-read sequencing, we accurately defined the 5'- and 3'-ends of transcripts that are both poly(A) tailed and 5'-capped in glucose starvation, followed by proteome analysis. Our results confirm previous experimentally validated loci with alternative isoforms and reveal several transcriptome-wide patterns. First, we show that sense-antisense gene pairs are more strongly anticorrelated when a time lag is taken into account. Second, we show that the glucose starvation response initially elicits a shortening of 3'-UTRs and poly(A) tails, followed by a shortening of the 5'-UTRs at later time points. These result in domain gains and losses in proteins involved in the stress response. Finally, the relatively poor overlap both between differentially expressed genes (DEGs), differential transcript usage events (DTUs), and differentially detected proteins (DDPs) highlight the need for further study on post-transcriptional regulation mechanisms in glucose starvation.

Intra- and inter-molecular regulation by intrinsically-disordered regions governs PUF protein RNA binding

PUF proteins are characterized by globular RNA-binding domains. They also interact with partner proteins that modulate their RNA-binding activities. Caenorhabditis elegans PUF protein fem-3 binding factor-2 (FBF-2) partners with intrinsically disordered Lateral Signaling Target-1 (LST-1) to regulate target mRNAs in germline stem cells. Here, we report that an intrinsically disordered region (IDR) at the C-terminus of FBF-2 autoinhibits its RNA-binding affinity by increasing the off rate for RNA binding. Moreover, the FBF-2 C-terminal region interacts with its globular RNA-binding domain at the same site where LST-1 binds. This intramolecular interaction restrains an electronegative cluster of amino acid residues near the 5' end of the bound RNA to inhibit RNA binding. LST-1 binding in place of the FBF-2 C-terminus therefore releases autoinhibition and increases RNA-binding affinity. This regulatory mechanism, driven by IDRs, provides a biochemical and biophysical explanation for the interdependence of FBF-2 and LST-1 in germline stem cell self-renewal.

JL. Additional principles that govern the release of pre-ribosomes from the nucleolus into the nucleoplasm in yeast

During eukaryotic ribosome biogenesis, pre-ribosomes travel from the nucleolus, where assembly is initiated, to the nucleoplasm and then are exported to the cytoplasm, where assembly concludes. Although nuclear export of pre-ribosomes has been extensively investigated, the release of pre-ribosomes from the nucleolus is an understudied phenomenon. Initial data indicate that unfolded rRNA interacts in trans with nucleolar components and that, when rRNA folds due to ribosomal protein (RP) binding, the number of trans interactions drops below the threshold necessary for nucleolar retention. To validate and expand on this idea, we performed a bioinformatic analysis of the protein components of the Saccharomyces cerevisiae ribosome assembly pathway. We found that ribosome biogenesis factors (RiBi factors) contain significantly more predicted trans interacting regions than RPs. We also analyzed cryo-EM structures of ribosome assembly intermediates to determine how nucleolar pre-ribosomes differ from post-nucleolar pre-ribosomes, specifically the capacity of RPs, RiBi factors, and rRNA components to interact in trans. We observed a significant decrease in the theoretical trans-interacting capability of pre-ribosomes between nucleolar and post-nucleolar stages of assembly due to the release of RiBi factors from particles and the folding of rRNA. Here, we provide a mechanism for the release of pre-ribosomes from the nucleolus.

Structural impairment of p53 C-terminal due to the effect of phosphorylation and acetylation: a study on the interdependence of PTM

The C-terminal of tumor suppressor protein p53 is intrinsically disordered while unbound. This particular segment often shows structural plasticity when bound to other binding partners. The disordered component undergoes a disordered to ordered transition upon recognition. Post-translational modifications (PTMs), namely phosphorylation and acetylation, significantly alter the structural motifs of the segment. Among the various types of PTMs, phosphorylation, and acetylation of p53 at both N- and C- terminals lead to stabilization and activation. It has been noted experimentally that phosphorylation often regulates (enhances or reduces) the acetylation at specific sites. The phosphorylation of Thr377 and Ser378 reduces the acetylation of Lys373 and Lys382. Mutations of Thr377 and Ser378 to neutral Ala enhance and phospho mimic Asp reduce the acetylation of Lys373 and Lys382. Simulations of several single-point and pair-wise mutated systems have been generated to compare how the presence or absence of phosphorylation favors or disfavors the acetylation by thermodynamic and conformational analysis. We are using implicit solvent replica exchange molecular dynamics simulations to get 200 ns well-converged conformational ensembles of each system. Different sets of systems having both single and double PTMs are simulated. The results admit the appreciable change in the secondary structural level upon specific PTM. Also, the residual structure of the unbound p53 with single-point PTM varies significantly with pair-wise modifications. These observations further shed light on the relationship between the interdependencies of the specific PTM sites and the secondary structural levels.Communicated by Ramaswamy H. Sarma.

SGF29 nuclear condensates reinforce cellular aging

Phase separation, a biophysical segregation of subcellular milieus referred as condensates, is known to regulate transcription, but its impacts on physiological processes are less clear. Here, we demonstrate the formation of liquid-like nuclear condensates by SGF29, a component of the SAGA transcriptional coactivator complex, during cellular senescence in human mesenchymal progenitor cells (hMPCs) and fibroblasts. The Arg 207 within the intrinsically disordered region is identified as the key amino acid residue for SGF29 to form phase separation. Through epigenomic and transcriptomic analysis, our data indicated that both condensate formation and H3K4me3 binding of SGF29 are essential for establishing its precise chromatin location, recruiting transcriptional factors and co-activators to target specific genomic loci, and initiating the expression of genes associated with senescence, such as CDKN1A. The formation of SGF29 condensates alone, however, may not be sufficient to drive H3K4me3 binding or achieve transactivation functions. Our study establishes a link between phase separation and aging regulation, highlighting nuclear condensates as a functional unit that facilitate shaping transcriptional landscapes in aging.

Liquid-liquid phase separation of α-synuclein is highly sensitive to sequence complexity

The Parkinson's-associated protein α-synuclein (α-syn) can undergo liquid-liquid phase separation (LLPS), which typically leads to the formation of amyloid fibrils. The coincidence of LLPS and amyloid formation has complicated the identification of the molecular determinants unique to LLPS of α-syn. Moreover, the lack of strategies to selectively perturb LLPS makes it difficult to dissect the biological roles specific to α-syn LLPS, independent of fibrillation. Herein, using a combination of subtle missense mutations, we show that LLPS of α-syn is highly sensitive to its sequence complexity. In fact, we find that even a highly conservative mutation (V16I) that increases sequence complexity without perturbing physicochemical and structural properties, is sufficient to reduce LLPS by 75%; this effect can be reversed by an adjacent V-to-I mutation (V15I) that restores the original sequence complexity. A18T, a complexity-enhancing PD-associated mutation, was likewise found to reduce LLPS, implicating sequence complexity in α-syn pathogenicity. Furthermore, leveraging the differences in LLPS propensities among different α-syn variants, we demonstrate that fibrillation of α-syn does not necessarily correlate with its LLPS. In fact, we identify mutations that selectively perturb LLPS or fibrillation of α-syn, unlike previously studied mutations. The variants and design principles reported herein should therefore empower future studies to disentangle these two phenomena and distinguish their (patho)biological roles.

CoMemMoRFPred: Sequence-based Prediction of MemMoRFs by Combining Predictors of Intrinsic Disorder, MoRFs and Disordered Lipid-binding Regions

Molecular recognition features (MoRFs) are a commonly occurring type of intrinsically disordered regions (IDRs) that undergo disorder-to-order transition upon binding to partner molecules. We focus on recently characterized and functionally important membrane-binding MoRFs (MemMoRFs). Motivated by the lack of computational tools that predict MemMoRFs, we use a dataset of experimentally annotated MemMoRFs to conceptualize, design, evaluate and release an accurate sequence-based predictor. We rely on state-of-the-art tools that predict residues that possess key characteristics of MemMoRFs, such as intrinsic disorder, disorder-to-order transition and lipid-binding. We identify and combine results from three tools that include flDPnn for the disorder prediction, DisoLipPred for the prediction of disordered lipid-binding regions, and MoRFCHiBiLight for the prediction of disorder-to-order transitioning protein binding regions. Our empirical analysis demonstrates that combining results produced by these three methods generates accurate predictions of MemMoRFs. We also show that use of a smoothing operator produces predictions that closely mimic the number and sizes of the native MemMoRF regions. The resulting CoMemMoRFPred method is available as an easy-to-use webserver at http://biomine.cs.vcu.edu/servers/CoMemMoRFPred. This tool will aid future studies of MemMoRFs in the context of exploring their abundance, cellular functions, and roles in pathologic phenomena.

Tutorial: a guide for the selection of fast and accurate computational tools for the prediction of intrinsic disorder in proteins

Intrinsic disorder is instrumental for a wide range of protein functions, and its analysis, using computational predictions from primary structures, complements secondary and tertiary structure-based approaches. In this Tutorial, we provide an overview and comparison of 23 publicly available computational tools with complementary parameters useful for intrinsic disorder prediction, partly relying on results from the Critical Assessment of protein Intrinsic Disorder prediction experiment. We consider factors such as accuracy, runtime, availability and the need for functional insights. The selected tools are available as web servers and downloadable programs, offer state-of-the-art predictions and can be used in a high-throughput manner. We provide examples and instructions for the selected tools to illustrate practical aspects related to the submission, collection and interpretation of predictions, as well as the timing and their limitations. We highlight two predictors for intrinsically disordered proteins, flDPnn as accurate and fast and IUPred as very fast and moderately accurate, while suggesting ANCHOR2 and MoRFchibi as two of the best-performing predictors for intrinsically disordered region binding. We link these tools to additional resources, including databases of predictions and web servers that integrate multiple predictive methods. Altogether, this Tutorial provides a hands-on guide to comparatively evaluating multiple predictors, submitting and collecting their own predictions, and reading and interpreting results. It is suitable for experimentalists and computational biologists interested in accurately and conveniently identifying intrinsic disorder, facilitating the functional characterization of the rapidly growing collections of protein sequences.

AAK1-like: A putative pseudokinase with potential roles in cargo uptake in bloodstream form Trypanosoma brucei parasites

Selection and internalization of cargo via clathrin-mediated endocytosis requires adaptor protein complexes. One complex, AP-2, acts during cargo selection at the plasma membrane. African trypanosomes lack all components of the AP-2 complex, except for a recently identified orthologue of the AP-2-associated protein kinase 1, AAK1. In characterized eukaryotes, AAK1 phosphorylates the μ2 subunit of the AP-2 complex to enhance cargo recognition and uptake into clathrin-coated vesicles. Here, we show that kinetoplastids encode not one, but two AAK1 orthologues: one (AAK1L2) is absent from salivarian trypanosomes, while the other (AAK1L1) lacks important kinase-specific residues in a range of trypanosomes. These AAK1L1 and AAK1L2 novelties reinforce suggestions of functional divergence in endocytic uptake within salivarian trypanosomes. Despite this, we show that AAK1L1 null mutant Trypanosoma brucei, while viable, display slowed proliferation, morphological abnormalities including swelling of the flagellar pocket, and altered cargo uptake. In summary, our data suggest an unconventional role for a putative pseudokinase during endocytosis and/or vesicular trafficking in T. brucei, independent of AP-2.

Coronaviruses spike glycoprotein endodomains: The sequence and structure-based comprehensive study

Any protein's flexibility or region makes it available to interact with many biomolecules in the cell. Specifically, such interactions in viruses help them to perform more functions despite having a smaller genome. Therefore, these flexible regions can be exciting and essential targets to be explored for their role in pathogenicity and therapeutic developments as they achieve essential interactions. In the continuation with our previous study on disordered analysis of SARS-CoV-2 spike protein's cytoplasmic tail (CTR), or endodomain, here we have explored the endodomain's disordered potential of six other coronaviruses using multiple bioinformatics approaches and molecular dynamics simulations. Based on the comprehensive analysis of its sequence and structural composition, we report the varying disorder propensity in endodomains of spike proteins of coronaviruses. The observations of this study may help to understand the importance of spike glycoprotein endodomain and creating therapeutic interventions against them.

Prediction of the interaction between Calloselasma rhodostoma venom-derived peptides and cancer-associated hub proteins: A computational study

The use of peptide drugs to treat cancer is gaining popularity because of their efficacy, fewer side effects, and several advantages over other properties. Identifying the peptides that interact with cancer proteins is crucial in drug discovery. Several approaches related to predicting peptide-protein interactions have been conducted. However, problems arise due to the high costs of resources and time and the smaller number of studies. This study predicts peptide-protein interactions using Random Forest, XGBoost, and SAE-DNN. Feature extraction is also performed on proteins and peptides using intrinsic disorder, amino acid sequences, physicochemical properties, position-specific assessment matrices, amino acid composition, and dipeptide composition. Results show that all algorithms perform equally well in predicting interactions between peptides derived from venoms and target proteins associated with cancer. However, XGBoost produces the best results with accuracy, precision, and area under the receiver operating characteristic curve of 0.859, 0.663, and 0.697, respectively. The enrichment analysis revealed that peptides from the Calloselasma rhodostoma venom targeted several proteins (ESR1, GOPC, and BRD4) related to cancer.

The Ustilago maydis AA10 LPMO is active on fungal cell wall chitin

Lytic polysaccharide monooxygenases (LPMOs) can perform oxidative cleavage of glycosidic bonds in carbohydrate polymers (e.g., cellulose, chitin), making them more accessible to hydrolytic enzymes. While most studies have so far mainly explored the role of LPMOs in a (plant) biomass conversion context, alternative roles and paradigms begin to emerge. The AA10 LPMOs are active on chitin and/or cellulose and mostly found in bacteria and in some viruses and archaea. Interestingly, AA10-encoding genes are also encountered in some pathogenic fungi of the Ustilaginomycetes class, such as Ustilago maydis, responsible for corn smut disease. Transcriptomic studies have shown the overexpression of the AA10 gene during the infectious cycle of U. maydis. In fact, U. maydis has a unique AA10 gene that codes for a catalytic domain appended with a C-terminal disordered region. To date, there is no public report on fungal AA10 LPMOs. In this study, we successfully produced the catalytic domain of this LPMO (UmAA10_cd) in Pichia pastoris and carried out its biochemical characterization. Our results show that UmAA10_cd oxidatively cleaves α- and β-chitin with C1 regioselectivity and boosts chitin hydrolysis by a GH18 chitinase from U. maydis (UmGH18A). Using a biologically relevant substrate, we show that UmAA10_cd exhibits enzymatic activity on U. maydis fungal cell wall chitin and promotes its hydrolysis by UmGH18A. These results represent an important step toward the understanding of the role of LPMOs in the fungal cell wall remodeling process during the fungal life cycle.IMPORTANCELytic polysaccharide monooxygenases (LPMOs) have been mainly studied in a biotechnological context for the efficient degradation of recalcitrant polysaccharides. Only recently, alternative roles and paradigms begin to emerge. In this study, we provide evidence that the AA10 LPMO from the phytopathogen Ustilago maydis is active against fungal cell wall chitin. Given that chitin-active LPMOs are commonly found in microbes, it is important to consider fungal cell wall as a potential target for this enigmatic class of enzymes.

Systematic identification of conditionally folded intrinsically disordered regions by AlphaFold2

The AlphaFold Protein Structure Database contains predicted structures for millions of proteins. For the majority of human proteins that contain intrinsically disordered regions (IDRs), which do not adopt a stable structure, it is generally assumed that these regions have low AlphaFold2 confidence scores that reflect low-confidence structural predictions. Here, we show that AlphaFold2 assigns confident structures to nearly 15% of human IDRs. By comparison to experimental NMR data for a subset of IDRs that are known to conditionally fold (i.e., upon binding or under other specific conditions), we find that AlphaFold2 often predicts the structure of the conditionally folded state. Based on databases of IDRs that are known to conditionally fold, we estimate that AlphaFold2 can identify conditionally folding IDRs at a precision as high as 88% at a 10% false positive rate, which is remarkable considering that conditionally folded IDR structures were minimally represented in its training data. We find that human disease mutations are nearly fivefold enriched in conditionally folded IDRs over IDRs in general and that up to 80% of IDRs in prokaryotes are predicted to conditionally fold, compared to less than 20% of eukaryotic IDRs. These results indicate that a large majority of IDRs in the proteomes of human and other eukaryotes function in the absence of conditional folding, but the regions that do acquire folds are more sensitive to mutations. We emphasize that the AlphaFold2 predictions do not reveal functionally relevant structural plasticity within IDRs and cannot offer realistic ensemble representations of conditionally folded IDRs.

KIF1C, an RNA transporting kinesin-3, undergoes liquid-liquid phase separation through its C-terminal disordered domain

The spatial distribution of mRNA is critical for local control of protein production. Recent studies have identified the kinesin-3 family member KIF1C as an RNA transporter. However, it is not clear how KIF1C interacts with RNA molecules. Here, we show that KIF1C's C-terminal tail domain is an intrinsically disordered region (IDR) containing a prion-like domain (PLD) that is unique compared to the C-terminal tails of other kinesin family members. In cells, KIF1C constructs undergo reversible formation of dynamic puncta that display physical properties of liquid condensates and incorporate RNA molecules in a sequence-selective manner. The IDR is necessary and sufficient for driving liquid-liquid phase separation (LLPS) but the condensate properties can be modulated by adjacent coiled-coil segments. The purified KIF1C IDR domain undergoes LLPS in vitro at near-endogenous nM concentrations in a salt-dependent manner. Deletion of the IDR abolished the ability of KIF1C to undergo LLPS and disrupted the distribution of mRNA cargoes to the cell periphery. Our work thus uncovers an intrinsic correlation between the LLPS activity of KIF1C and its role as an RNA transporter. In addition, as the first kinesin motor reported to undergo LLPS, our work reveals a previously uncharacterized mode of motor-cargo interaction that extends our understanding of the behavior of cytoskeletal motor proteins.

Evolution of Intrinsic Disorder in Protein Loops

Intrinsic disorder accounts for the flexibility of protein loops, molecular building blocks that are largely responsible for the processes and molecular functions of the living world. While loops likely represent early structural forms that served as intermediates in the emergence of protein structural domains, their origin and evolution remain poorly understood. Here, we conduct a phylogenomic survey of disorder in loop prototypes sourced from the ArchDB classification. Tracing prototypes associated with protein fold families along an evolutionary chronology revealed that ancient prototypes tended to be more disordered than their derived counterparts, with ordered prototypes developing later in evolution. This highlights the central evolutionary role of disorder and flexibility. While mean disorder increased with time, a minority of ordered prototypes exist that emerged early in evolutionary history, possibly driven by the need to preserve specific molecular functions. We also revealed the percolation of evolutionary constraints from higher to lower levels of organization. Percolation resulted in trade-offs between flexibility and rigidity that impacted prototype structure and geometry. Our findings provide a deep evolutionary view of the link between structure, disorder, flexibility, and function, as well as insights into the evolutionary role of intrinsic disorder in loops and their contribution to protein structure and function.

BRCA1 frameshift variants leading to extended incorrect protein C termini

Carriers of BRCA1 germline pathogenic variants are at substantially higher risk of developing breast and ovarian cancer than the general population. Accurate identification of at-risk individuals is crucial for risk stratification and the implementation of targeted preventive and therapeutic interventions. Despite significant progress in variant classification efforts, a sizable portion of reported BRCA1 variants remain as variants of uncertain clinical significance (VUSs). Variants leading to premature protein termination and loss of essential functional domains are typically classified as pathogenic. However, the impact of frameshift variants that result in an extended incorrect terminus is not clear. Using validated functional assays, we conducted a systematic functional assessment of 17 previously reported BRCA1 extended incorrect terminus variants (EITs) and concluded that 16 constitute loss-of-function variants. This suggests that most EITs are likely to be pathogenic. However, one variant, c.5578dup, displayed a protein expression level, affinity to known binding partners, and activity in transcription and homologous recombination assays comparable to the wild-type BRCA1 protein. Twenty-three additional carriers of c.5578dup were identified at a US clinical diagnostic lab and assessed using a family history likelihood model providing, in combination with the functional data, a likely benign interpretation. These results, consistent with family history data in the current study and available data from ClinVar, indicate that most, but not all, BRCA1 variants leading to an extended incorrect terminus constitute loss-of-function variants and underscore the need for comprehensive assessment of individual variants.

The Intrinsically Disordered N Terminus in Atg12 from Yeast Is Necessary for the Functional Structure of the Protein

The Atg12 protein in yeast is an indispensable polypeptide in the highly conserved ubiquitin-like conjugation system operating in the macroautophagy/autophagy pathway. Atg12 is covalently conjugated to Atg5 through the action of Atg7 and Atg10; the Atg12-Atg5 conjugate binds Atg16 to form an E3 ligase that functions in a separate conjugation pathway involving Atg8. Atg12 is comprised of a ubiquitin-like (UBL) domain preceded at the N terminus by an intrinsically disordered protein region (IDPR), a domain that comprises a major portion of the protein but remains elusive in its conformation and function. Here, we show that the IDPR in unconjugated Atg12 is positioned in proximity to the UBL domain, a configuration that is important for the functional structure of the protein. A major deletion in the IDPR disrupts intactness of the UBL domain at the unconjugated C terminus, and a mutation in the predicted α0 helix in the IDPR prevents Atg12 from binding to Atg7 and Atg10, which ultimately affects the protein function in the ubiquitin-like conjugation cascade. These findings provide evidence that the IDPR is an indispensable part of the Atg12 protein from yeast.

A cyclin-dependent kinase-mediated phosphorylation switch of disordered protein condensation

Cell cycle transitions result from global changes in protein phosphorylation states triggered by cyclin-dependent kinases (CDKs). To understand how this complexity produces an ordered and rapid cellular reorganisation, we generated a high-resolution map of changing phosphosites throughout unperturbed early cell cycles in single Xenopus embryos, derived the emergent principles through systems biology analysis, and tested them by biophysical modelling and biochemical experiments. We found that most dynamic phosphosites share two key characteristics: they occur on highly disordered proteins that localise to membraneless organelles, and are CDK targets. Furthermore, CDK-mediated multisite phosphorylation can switch homotypic interactions of such proteins between favourable and inhibitory modes for biomolecular condensate formation. These results provide insight into the molecular mechanisms and kinetics of mitotic cellular reorganisation.

Ohnologs and SSD Paralogs Differ in Genomic and Expression Features Related to Dosage Constraints

Gene duplication is recognized as a critical process in genome evolution; however, many questions about this process remain unanswered. Although gene duplicability has been observed to differ by duplication mechanism and evolutionary rate, there is so far no broad characterization of its determinants. Many features correlate with this difference in duplicability; however, our ability to exploit these observations to advance our understanding of the role of duplication in evolution is hampered by limitations within existing work. In particular, the existence of methodological differences across studies impedes meaningful comparison. Here, we use consistent definitions of duplicability in the human lineage to explore these associations, allow resolution of the impact of confounding factors, and define the overall relevance of individual features. Using a classifier approach and controlling for the confounding effect of duplicate longevity, we find a subset of gene features important in differentiating genes duplicable by small-scale duplication from those duplicable by whole-genome duplication, revealing critical roles for gene dosage and expression costs in duplicability. We further delve into patterns of functional enrichment and find a lack of constraint on duplicate retention in any context for genes duplicable by small-scale duplication.

Intermolecular interactions underlie protein/peptide phase separation irrespective of sequence and structure at crowded milieu

Liquid-liquid phase separation (LLPS) has emerged as a crucial biological phenomenon underlying the sequestration of macromolecules (such as proteins and nucleic acids) into membraneless organelles in cells. Unstructured and intrinsically disordered domains are known to facilitate multivalent interactions driving protein LLPS. We hypothesized that LLPS could be an intrinsic property of proteins/polypeptides but with distinct phase regimes irrespective of their sequence and structure. To examine this, we studied many (a total of 23) proteins/polypeptides with different structures and sequences for LLPS study in the presence and absence of molecular crowder, polyethylene glycol (PEG-8000). We showed that all proteins and even highly charged polypeptides (under study) can undergo liquid condensate formation, however with different phase regimes and intermolecular interactions. We further demonstrated that electrostatic, hydrophobic, and H-bonding or a combination of such intermolecular interactions plays a crucial role in individual protein/peptide LLPS.

Promoter R-Loops Recruit U2AF1 to Modulate Its Phase Separation and RNA Splicing

R-loops and guanine quadruplexes (G4s) are secondary structures of nucleic acids that are ubiquitously present in cells and are enriched in promoter regions of genes. By employing a bioinformatic approach based on overlap analysis of transcription factor chromatin immunoprecipitation sequencing (ChIP-seq) data sets, we found that many splicing factors, including U2AF1 whose recognition of the 3' splicing site is crucial for pre-mRNA splicing, exhibit pronounced enrichment at endogenous R-loop- and DNA G4-structure loci in promoter regions of human genes. We also revealed that U2AF1 binds directly to R-loops and DNA G4 structures at a low-nM binding affinity. Additionally, we showed the ability of U2AF1 to undergo phase separation, which could be stimulated by binding with R-loops, but not duplex DNA, RNA/DNA hybrid, DNA G4, or single-stranded RNA. We also demonstrated that U2AF1 binds to promoter R-loops in human cells, and this binding competes with U2AF1's interaction with 3' splicing site and leads to augmented distribution of RNA polymerase II (RNAPII) to promoters over gene bodies, thereby modulating cotranscriptional pre-mRNA splicing. Together, we uncovered a group of candidate proteins that can bind to both R-loops and DNA G4s, revealed the direct and strong interactions of U2AF1 with these nucleic acid structures, and established a biochemical rationale for U2AF1's occupancy in gene promoters. We also unveiled that interaction with R-loops promotes U2AF1's phase separation, and our work suggests that U2AF1 modulates pre-mRNA splicing by regulating RNAPII's partition in transcription initiation versus elongation.

Tejas functions as a core component in nuage assembly and precursor processing in Drosophila piRNA biogenesis

PIWI-interacting RNAs (piRNAs), which protect genome from the attack by transposons, are produced and amplified in membraneless granules called nuage. In Drosophila, PIWI family proteins, Tudor-domain-containing (Tdrd) proteins, and RNA helicases are assembled and form nuage to ensure piRNA production. However, the molecular functions of the Tdrd protein Tejas (Tej) in piRNA biogenesis remain unknown. Here, we conduct a detailed analysis of the subcellular localization of fluorescently tagged nuage proteins and behavior of piRNA precursors. Our results demonstrate that Tej functions as a core component that recruits Vasa (Vas) and Spindle-E (Spn-E) into nuage granules through distinct motifs, thereby assembling nuage and engaging precursors for further processing. Our study also reveals that the low-complexity region of Tej regulates the mobility of Vas. Based on these results, we propose that Tej plays a pivotal role in piRNA precursor processing by assembling Vas and Spn-E into nuage and modulating the mobility of nuage components.

Intrinsic disorder may drive the interaction of PROS1 and MERTK in uveal melanoma

BACKGROUND

Class 2 uveal melanomas are associated with the inactivation of the BRCA1 ((breast cancer type 1 susceptibility protein)-associated protein 1 (BAP1)) gene. Inactivation of BAP1 promotes the upregulation of vitamin K-dependent protein S (PROS1), which interacts with the tyrosine-protein kinase Mer (MERTK) receptor on M2 macrophages to induce an immunosuppressive environment.

METHODS

We simulated the interaction of PROS1 with MERTK with ColabFold. We evaluated PROS1 and MERTK for the presence of intrinsically disordered protein regions (IDPRs) and disorder-to-order (DOT) regions to understand their protein-protein interaction (PPI). We first evaluated the structure of each protein with AlphaFold. We then analyzed specific sequence-based features of the PROS1 and MERTK with a suite of bioinformatics tools.

RESULTS

With high-resolution, moderate confidence, we successfully modeled the interaction between PROS1 and MERTK (predicted local distance difference test score (pDLLT) = 70.68). Our structural analysis qualitatively demonstrated IDPRs (i.e., spaghetti-like entities) in PROS1 and MERK. PROS1 was 23.37 % disordered, and MERTK was 23.09 % disordered, classifying them as moderately disordered and flexible proteins. PROS1 was significantly enriched in cysteine, the most order-promoting residue (p-value <0.05). Our IUPred analysis demonstrated that there are two disorder-to-order transition (DOT) regions in PROS1. MERTK was significantly enriched in proline, the most disorder-promoting residue (p-value <0.05), but did not contain DOT regions. Our STRING analysis demonstrated that the PPI network between PROS1 and MERTK is more complex than their assumed one-to-one binding (p-value <2.0 × 10-6).

CONCLUSION

Our findings present a novel prediction for the interaction between PROS1 and MERTK. Our findings show that PROS1 and MERTK contain elements of intrinsic disorder. PROS1 has two DOT regions that are attractive immunotherapy targets. We recommend that IDPRs and DOT regions found in PROS1 and MERTK should be considered when developing immunotherapies targeting this PPI.

Computational analysis of regulatory regions in human protein kinases

Eukaryotic proteins often feature modular domain structures comprising globular domains that are connected by linker regions and intrinsically disordered regions that may contain important functional motifs. The intramolecular interactions of globular domains and nonglobular regions can play critical roles in different aspects of protein function. However, studying these interactions and their regulatory roles can be challenging due to the flexibility of nonglobular regions, the long insertions separating interacting modules, and the transient nature of some interactions. Obtaining the experimental structures of multiple domains and functional regions is more difficult than determining the structures of individual globular domains. High-quality structural models generated by AlphaFold offer a unique opportunity to study intramolecular interactions in eukaryotic proteins. In this study, we systematically explored intramolecular interactions between human protein kinase domains (KDs) and potential regulatory regions, including globular domains, N- and C-terminal tails, long insertions, and distal nonglobular regions. Our analysis identified intramolecular interactions between human KDs and 35 different types of globular domains, exhibiting a variety of interaction modes that could contribute to orthosteric or allosteric regulation of kinase activity. We also identified prevalent interactions between human KDs and their flanking regions (N- and C-terminal tails). These interactions exhibit group-specific characteristics and can vary within each specific kinase group. Although long-range interactions between KDs and nonglobular regions are relatively rare, structural details of these interactions offer new insights into the regulation mechanisms of several kinases, such as HASPIN, MAPK7, MAPK15, and SIK1B.

Disordered regions endow structural flexibility to shell proteins and function towards shell-enzyme interactions in 1,2-propanediol utilization microcompartment

Intrinsically disordered regions in proteins have been functionally linked to the protein-protein interactions and genesis of several membraneless organelles. Depending on their residual makeup, hydrophobicity or charge distribution they may remain in extended form or may assume certain conformations upon biding to a partner protein or peptide. The present work investigates the distribution and potential roles of disordered regions in the integral proteins of 1,2-propanediol utilization microcompartments. We use bioinformatics tools to identify the probable disordered regions in the shell proteins and enzyme of the 1,2-propanediol utilization microcompartment. Using a combination of computational modelling and biochemical techniques we elucidate the role of disordered terminal regions of a major shell protein and enzyme. Our findings throw light on the importance of disordered regions in the self-assembly, providing flexibility to shell protein and mediating its interaction with a native enzyme.Communicated by Ramaswamy H. Sarma.

SAFB associates with nascent RNAs and can promote gene expression in mouse embryonic stem cells

Scaffold attachment factor B (SAFB) is a conserved RNA-binding protein that is essential for early mammalian development. However, the functions of SAFB in mouse embryonic stem cells (ESCs) have not been characterized. Using RNA immunoprecipitation followed by RNA-seq (RIP-seq), we examined the RNAs associated with SAFB in wild-type and SAFB/SAFB2 double-knockout ESCs. SAFB predominantly associated with introns of protein-coding genes through purine-rich motifs. The transcript most enriched in SAFB association was the lncRNA Malat1, which also contains a purine-rich region in its 5' end. Knockout of SAFB/SAFB2 led to differential expression of approximately 1000 genes associated with multiple biological processes, including apoptosis, cell division, and cell migration. Knockout of SAFB/SAFB2 also led to splicing changes in a set of genes that were largely distinct from those that exhibited changes in expression level. The spliced and nascent transcripts of many genes whose expression levels were positively regulated by SAFB also associated with high levels of SAFB, implying that SAFB binding promotes their expression. Reintroduction of SAFB into double-knockout cells restored gene expression toward wild-type levels, an effect again observable at the level of spliced and nascent transcripts. Proteomics analysis revealed a significant enrichment of nuclear speckle-associated and RS domain-containing proteins among SAFB interactors. Neither Xist nor Polycomb functions were dramatically altered in SAFB/2 knockout ESCs. Our findings suggest that among other potential functions in ESCs, SAFB promotes the expression of certain genes through its ability to bind nascent RNA.

Defining the condensate landscape of fusion oncoproteins

Fusion oncoproteins (FOs) arise from chromosomal translocations in ~17% of cancers and are often oncogenic drivers. Although some FOs can promote oncogenesis by undergoing liquid-liquid phase separation (LLPS) to form aberrant biomolecular condensates, the generality of this phenomenon is unknown. We explored this question by testing 166 FOs in HeLa cells and found that 58% formed condensates. The condensate-forming FOs displayed physicochemical features distinct from those of condensate-negative FOs and segregated into distinct feature-based groups that aligned with their sub-cellular localization and biological function. Using Machine Learning, we developed a predictor of FO condensation behavior, and discovered that 67% of ~3000 additional FOs likely form condensates, with 35% of those predicted to function by altering gene expression. 47% of the predicted condensate-negative FOs were associated with cell signaling functions, suggesting a functional dichotomy between condensate-positive and -negative FOs. Our Datasets and reagents are rich resources to interrogate FO condensation in the future.

The Origin of Discrepancies between Predictions and Annotations in Intrinsically Disordered Proteins

Disorder prediction methods that can discriminate between ordered and disordered regions have contributed fundamentally to our understanding of the properties and prevalence of intrinsically disordered proteins (IDPs) in proteomes as well as their functional roles. However, a recent large-scale assessment of the performance of these methods indicated that there is still room for further improvements, necessitating novel approaches to understand the strengths and weaknesses of individual methods. In this study, we compared two methods, IUPred and disorder prediction, based on the pLDDT scores derived from AlphaFold2 (AF2) models. We evaluated these methods using a dataset from the DisProt database, consisting of experimentally characterized disordered regions and subsets associated with diverse experimental methods and functions. IUPred and AF2 provided consistent predictions in 79% of cases for long disordered regions; however, for 15% of these cases, they both suggested order in disagreement with annotations. These discrepancies arose primarily due to weak experimental support, the presence of intermediate states, or context-dependent behavior, such as binding-induced transitions. Furthermore, AF2 tended to predict helical regions with high pLDDT scores within disordered segments, while IUPred had limitations in identifying linker regions. These results provide valuable insights into the inherent limitations and potential biases of disorder prediction methods.

Systematic identification of inter-chromosomal interaction networks supports the existence of RNA factories

Most studies of genome organization have focused on intra-chromosomal (cis) contacts because they harbor key features such as DNA loops and topologically associating domains. Inter-chromosomal (trans) contacts have received much less attention, and tools for interrogating potential biologically relevant trans structures are lacking. Here, we develop a computational framework to identify sets of loci that jointly interact in trans from Hi-C data. This method, trans-C, initiates probabilistic random walks with restarts from a set of seed loci to traverse an input Hi-C contact network, thereby identifying sets of trans-contacting loci. We validate trans-C in three increasingly complex models of established trans contacts: the Plasmodium falciparum var genes, the mouse olfactory receptor "Greek islands", and the human RBM20 cardiac splicing factory. We then apply trans-C to systematically test the hypothesis that genes co-regulated by the same trans-acting element (i.e., a transcription or splicing factor) co-localize in three dimensions to form "RNA factories" that maximize the efficiency and accuracy of RNA biogenesis. We find that many loci with multiple binding sites of the same transcription factor interact with one another in trans, especially those bound by transcription factors with intrinsically disordered domains. Similarly, clustered binding of a subset of RNA binding proteins correlates with trans interaction of the encoding loci. These findings support the existence of trans interacting chromatin domains (TIDs) driven by RNA biogenesis. Trans-C provides an efficient computational framework for studying these and other types of trans interactions, empowering studies of a poorly understood aspect of genome architecture.

The architecture of binding cooperativity between densely bound transcription factors

The binding of transcription factors (TFs) along genomes is restricted to a subset of sites containing their preferred motifs. TF-binding specificity is often attributed to the co-binding of interacting TFs; however, apart from specific examples, this model remains untested. Here, we define dependencies among budding yeast TFs that localize to overlapping promoters by profiling the genome-wide consequences of co-depleting multiple TFs. We describe unidirectional interactions, revealing Msn2 as a central factor allowing TF binding at its target promoters. By contrast, no case of mutual cooperation was observed. Particularly, Msn2 retained binding at its preferred promoters upon co-depletion of fourteen similarly bound TFs. Overall, the consequences of TF co-depletions were moderate, limited to a subset of promoters, and failed to explain the role of regions outside the DNA-binding domain in directing TF-binding preferences. Our results call for re-evaluating the role of cooperative interactions in directing TF-binding preferences.

The RNA-binding domain of hnRNP U extends beyond the RGG/RG motifs

Heterogeneous nuclear ribonucleoprotein U (hnRNP U) is a ubiquitously expressed protein that regulates chromatin architecture through its interactions with numerous DNA, protein, and RNA partners. The RNA-binding domain (RBD) of hnRNP U was previously mapped to an RGG/RG element within its disordered C-terminal region, but little is understood about its binding mode and potential for selective RNA recognition. Analysis of publicly available hnRNP U enhanced UV crosslinking and immunoprecipitation (eCLIP) data identified high-confidence binding sites within human RNAs. We synthesized a set of diverse RNAs encompassing eleven of these identified crosslink sites for biochemical characterization using a combination of fluorescence anisotropy and electrophoretic mobility shift assays. These in vitro binding experiments with a rationally designed set of RNAs and hnRNP U domains revealed that the RGG/RG element is a small part of a more expansive RBD that encompasses most of the disordered C-terminal region. This RBD contains a second, previously experimentally uncharacterized RGG/RG element with RNA-binding properties comparable to the canonical RGG/RG element. These RGG/RG elements serve redundant functions, with neither serving as the primary RBD. While in isolation each RGG/RG element has modest affinity for RNA, together they significantly enhance the association of hnRNP U with RNA, enabling binding of most of the designed RNA set with low to mid-nanomolar binding affinities. Identification and characterization of the complete hnRNP U RBD highlights the perils of a reductionist approach to defining biochemical activities in this system and paves the way for a detailed investigation of its RNA-binding specificity.

Flanking regions, amyloid cores, and polymorphism: the potential interplay underlying structural diversity

The β-sheet-rich amyloid core is the defining feature of protein aggregates associated with neurodegenerative disorders. Recent investigations have revealed that there exist multiple examples of the same protein, with the same sequence, forming a variety of amyloid cores with distinct structural characteristics. These structural variants, termed as polymorphs, are hypothesized to influence the pathological profile and the progression of different neurodegenerative diseases, giving rise to unique phenotypic differences. Thus, identifying the origin and properties of these structural variants remain a focus of studies, as a preliminary step in the development of therapeutic strategies. Here, we review the potential role of the flanking regions of amyloid cores in inducing polymorphism. These regions, adjacent to the amyloid cores, show a preponderance for being structurally disordered, imbuing them with functional promiscuity. The dynamic nature of the flanking regions can then manifest in the form of conformational polymorphism of the aggregates. We take a closer look at the sequences flanking the amyloid cores, followed by a review of the polymorphic aggregates of the well-characterized proteins amyloid-β, α-synuclein, Tau, and TDP-43. We also consider different factors that can potentially influence aggregate structure and how these regions can be viewed as novel targets for therapeutic strategies by utilizing their unique structural properties.

Impact of Asp/Glu-ADP-ribosylation on protein-protein interaction and protein function

PARylation plays critical role in regulating multiple cellular processes such as DNA damage response and repair, transcription, RNA processing, and stress response. More than 300 human proteins have been found to be modified by PARylation on acidic residues, that is, Asp (D) and Glu (E). We used the deep-learning tool AlphaFold to predict protein-protein interactions (PPIs) and their interfaces for these proteins based on coevolution signals from joint multiple sequence alignments (MSAs). AlphaFold predicted 260 confident PPIs involving PARylated proteins, and about one quarter of these PPIs have D/E-PARylation sites in their predicted PPI interfaces. AlphaFold predictions offer novel insights into the mechanisms of PARylation regulations by providing structural details of the PPI interfaces. D/E-PARylation sites have a preference to occur in coil regions and disordered regions, and PPI interfaces containing D/E-PARylation sites tend to occur between short linear sequence motifs in disordered regions and globular domains. The hub protein PCNA is predicted to interact with more than 20 proteins via the common PIP box motif and the structurally variable flanking regions. D/E-PARylation sites were found in the interfaces of key components of the RNA transcription and export complex, the SF3a spliceosome complex, and H/ACA and C/D small nucleolar ribonucleoprotein complexes, suggesting that systematic PARylation have a profound effect in regulating multiple RNA-related processes such as RNA nuclear export, splicing, and modification. Finally, PARylation of SUMO2 could modulate its interaction with CHAF1A, thereby representing a potential mechanism for the cross-talk between PARylation and SUMOylation in regulation of chromatin remodeling.

Analysis and Prediction of Pathogen Nucleic Acid Specificity for Toll-like Receptors in Vertebrates

Identification of key sequence, expression and function related features of nucleic acid-sensing host proteins is of fundamental importance to understand the dynamics of pathogen-specific host responses. To meet this objective, we considered toll-like receptors (TLRs), a representative class of membrane-bound sensor proteins, from 17 vertebrate species covering mammals, birds, reptiles, amphibians, and fishes in this comparative study. We identified the molecular signatures of host TLRs that are responsible for sensing pathogen nucleic acids or other pathogen-associated molecular patterns (PAMPs), and potentially play important roles in host defence mechanism. Interestingly, our findings reveal that such host-specific features are directly related to the strand (single or double) specificity of nucleic acid from pathogens. However, during host-pathogen interactions, such features were unable to explain the pathogenic PAMP (i.e., DNA, RNA or other) selectivity, suggesting a more complex mechanism. Using these features, we developed a number of machine learning models, of which Random Forest achieved a high performance (94.57% accuracy) to predict strand specificity of TLRs from protein-derived features. We applied the trained model to propose strand specificity of some previously uncharacterized distinct fish-specific novel TLRs (TLR18, TLR23, TLR24, TLR25, TLR27).

Biology of the hepatitis B virus (HBV) core and capsid assembly modulators (CAMs) for chronic hepatitis B (CHB) cure

Hepatitis B virus (HBV) is a hepadnavirus, a small DNA virus that infects liver tissue, with some unusual replication steps that share similarities to retroviruses. HBV infection can lead to chronic hepatitis B (CHB), a life-long infection associated with significant risks of liver disease, especially if untreated. HBV is a significant global health problem, with hundreds of millions currently living with CHB. Currently approved strategies to prevent or inhibit HBV are highly effective, however, a cure for CHB has remained elusive. To achieve a cure, elimination of the functionally integrated HBV covalently closed chromosomal DNA (cccDNA) genome is required. The capsid core is an essential component of HBV replication, serving roles when establishing infection and in creating new virions. Over the last two and a half decades, significant efforts have been made to find and characterize antivirals that target the capsid, specifically the HBV core protein (Cp). The antivirals that interfere with the kinetics and morphology of the capsid, termed capsid assembly modulators (CAMs), are extremely potent, and clinical investigations indicate they are well tolerated and highly effective. Several CAMs offer the potential to cure CHB by decreasing the cccDNA pools. Here, we review the biology of the HBV capsid, focused on Cp, and the development of inhibitors that target it.

Different mechanisms of synapsin-induced vesicle clustering at inhibitory and excitatory synapses

Synapsins cluster synaptic vesicles (SVs) to provide a reserve pool (RP) of SVs that maintains synaptic transmission during sustained activity. However, it is unclear how synapsins cluster SVs. Here we show that either liquid-liquid phase separation (LLPS) or tetramerization-dependent cross-linking can cluster SVs, depending on whether a synapse is excitatory or inhibitory. Cell-free reconstitution reveals that both mechanisms can cluster SVs, with tetramerization being more effective. At inhibitory synapses, perturbing synapsin-dependent LLPS impairs SV clustering and synchronization of gamma-aminobutyric acid (GABA) release, while preventing synapsin tetramerization does not. At glutamatergic synapses, the opposite is true: synapsin tetramerization enhances clustering of glutamatergic SVs and mobilization of these SVs from the RP, while synapsin LLPS does not. Comparison of inhibitory and excitatory transmission during prolonged synaptic activity reveals that synapsin LLPS serves as a brake to limit GABA release, while synapsin tetramerization enables rapid mobilization of SVs from the RP to sustain glutamate release.

Intrinsic disorder in PRAME and its role in uveal melanoma

INTRODUCTION

The PReferentially expressed Antigen in MElanoma (PRAME) protein has been shown to be an independent biomarker for increased risk of metastasis in Class 1 uveal melanomas (UM). Intrinsically disordered proteins and regions of proteins (IDPs/IDPRs) are proteins that do not have a well-defined three-dimensional structure and have been linked to neoplastic development. Our study aimed to evaluate the presence of intrinsic disorder in PRAME and the role these structureless regions have in PRAME( +) Class 1 UM.

METHODS

A bioinformatics study to characterize PRAME's propensity for the intrinsic disorder. We first used the AlphaFold tool to qualitatively assess the protein structure of PRAME. Then we used the Compositional Profiler and a set of per-residue intrinsic disorder predictors to quantify the intrinsic disorder. The Database of Disordered Protein Prediction (D2P2) platform, IUPred, FuzDrop, fIDPnn, AUCpred, SPOT-Disorder2, and metapredict V2 allowed us to evaluate the potential functional disorder of PRAME. Additionally, we used the Search Tool for the Retrieval of Interacting Genes (STRING) to analyze PRAME's potential interactions with other proteins.

RESULTS

Our structural analysis showed that PRAME contains intrinsically disordered protein regions (IDPRs), which are structureless and flexible. We found that PRAME is significantly enriched with serine (p-value < 0.05), a disorder-promoting amino acid. PRAME was found to have an average disorder score of 16.49% (i.e., moderately disordered) across six per-residue intrinsic disorder predictors. Our IUPred analysis revealed the presence of disorder-to-order transition (DOT) regions in PRAME near the C-terminus of the protein (residues 475-509). The D2P2 platform predicted a region from approximately 140 and 175 to be highly concentrated with post-translational modifications (PTMs). FuzDrop predicted the PTM hot spot of PRAME to be a droplet-promoting region and an aggregation hotspot. Finally, our analysis using the STRING tool revealed that PRAME has significantly more interactions with other proteins than expected for randomly selected proteins of the same size, with the ability to interact with 84 different partners (STRING analysis result: p-value < 1.0 × 10-16; model confidence: 0.400).

CONCLUSION

Our study revealed that PRAME has IDPRs that are possibly linked to its functionality in the context of Class 1 UM. The regions of functionality (i.e., DOT regions, PTM sites, droplet-promoting regions, and aggregation hotspots) are localized to regions of high levels of disorder. PRAME has a complex protein-protein interaction (PPI) network that may be secondary to the structureless features of the polypeptide. Our findings contribute to our understanding of UM and suggest that IDPRs and DOT regions in PRAME may be targeted in developing new therapies for this aggressive cancer. Video Abstract.

TRGdb: a universal resource for the exploration of taxonomically restricted genes in bacteria

The TRGdb database is a resource dedicated to taxonomically restricted genes (TRGs) in bacteria. It provides a comprehensive collection of genes that are specific to different genera and species, according to the latest release of bacterial taxonomy. The user interface allows for easy browsing and searching as well as sequence similarity exploration. The website also provides information on each TRG protein sequence, including its level of disorder, complexity and tendency to aggregate. TRGdb is a valuable resource for gaining a deeper understanding of the TRG-associated, unique features, and characteristics of bacterial organisms. Database URL www.combio.pl/trgdb.

DDX24 Mutation Alters NPM1 Phase Behavior and Disrupts Nucleolar Homeostasis in Vascular Malformations

Point mutations in the DEAD-box helicase DDX24 are associated with vascular malformations such as multi-organ venous and lymphatic defect (MOVLD) syndrome and Budd-Chiari syndrome, with the pathogenesis largely uncharacterized. DDX24 is mainly located in the nucleolus, where nucleophosmin (NPM1) regulates nucleolar homeostasis via liquid-liquid phase separation (LLPS). However, the connection between DDX24 and NPM1 in vascular malformation remains elusive. Here we demonstrated that DDX24 formed biomolecular condensates in vitro and the mutated DDX24 protein, DDX24E271K, partitioned less into the nucleoli in tissues from patients with MOVLD syndrome and cultured endothelial cells (ECs), altering nucleolar morphology. Furthermore, DDX24 was directly associated with NPM1 to regulate its phase behavior as a client in the nucleolar granular component (GC). Functionally, we showed that DDX24 was essential in maintaining nucleolar homeostasis of ECs and that either mutation or knockdown of DDX24 led to the dysfunction of ribosome biogenesis and the elevated capability of cell migration and tube formation. Our findings illustrate how DDX24 mutation affects nucleolar structure and function by regulating the phase behavior of NPM1 in the setting of vascular malformation.

Over-expression of Dyrk1A affects bleeding by modulating plasma fibronectin and fibrinogen level in mice

Down syndrome is the most common chromosomal abnormality in humans. Patients with Down syndrome have hematologic disorders, including mild to moderate thrombocytopenia. In case of Down syndrome, thrombocytopenia is not associated with bleeding, and it remains poorly characterized regarding molecular mechanisms. We investigated the effects of overexpression of Dyrk1A, an important factor contributing to some major Down syndrome phenotypes, on platelet number and bleeding in mice. Mice overexpressing Dyrk1A have a decrease in platelet number by 20%. However, bleeding time was found to be reduced by 50%. The thrombocytopenia and the decreased bleeding time observed were not associated to an abnormal platelet receptors expression, to a defect of platelet activation by ADP, thrombin or convulxin, to the presence of activated platelets in the circulation or to an abnormal half-life of the platelets. To propose molecular mechanisms explaining this discrepancy, we performed a network analysis of Dyrk1A interactome and demonstrated that Dyrk1A, fibronectin and fibrinogen interact indirectly through two distinct clusters of proteins. Moreover, in mice overexpressing Dyrk1A, increased plasma fibronectin and fibrinogen levels were found, linked to an increase of the hepatic fibrinogen production. Our results indicate that overexpression of Dyrk1A in mice induces decreased bleeding consistent with increased plasma fibronectin and fibrinogen levels, revealing a new role of Dyrk1A depending on its indirect interaction with these two proteins.

Mechanosensitive ion channels MSL8, MSL9, and MSL10 have environmentally sensitive intrinsically disordered regions with distinct biophysical characteristics in vitro

Intrinsically disordered protein regions (IDRs) are highly dynamic sequences that rapidly sample a collection of conformations over time. In the past several decades, IDRs have emerged as a major component of many proteomes, comprising ~30% of all eukaryotic protein sequences. Proteins with IDRs function in a wide range of biological pathways and are notably enriched in signaling cascades that respond to environmental stresses. Here, we identify and characterize intrinsic disorder in the soluble cytoplasmic N-terminal domains of MSL8, MSL9, and MSL10, three members of the MscS-like (MSL) family of mechanosensitive ion channels. In plants, MSL channels are proposed to mediate cell and organelle osmotic homeostasis. Bioinformatic tools unanimously predicted that the cytosolic N-termini of MSL channels are intrinsically disordered. We examined the N-terminus of MSL10 (MSL10N) as an exemplar of these IDRs and circular dichroism spectroscopy confirms its disorder. MSL10N adopted a predominately helical structure when exposed to the helix-inducing compound trifluoroethanol (TFE). Furthermore, in the presence of molecular crowding agents, MSL10N underwent structural changes and exhibited alterations to its homotypic interaction favorability. Lastly, interrogations of collective behavior via in vitro imaging of condensates indicated that MSL8N, MSL9N, and MSL10N have sharply differing propensities for self-assembly into condensates, both inherently and in response to salt, temperature, and molecular crowding. Taken together, these data establish the N-termini of MSL channels as intrinsically disordered regions with distinct biophysical properties and the potential to respond uniquely to changes in their physiochemical environment.

Nondomain biopolymers: Flexible molecular strategies to acquire biological functions

A long-standing assumption in molecular biology posits that the conservation of protein and nucleic acid sequences emphasizes the functional significance of biomolecules. These conserved sequences fold into distinct secondary and tertiary structures, enable highly specific molecular interactions, and regulate complex yet organized molecular processes within living cells. However, recent evidence suggests that biomolecules can also function through primary sequence regions that lack conservation across species or gene families. These regions typically do not form rigid structures, and their inherent flexibility is critical for their functional roles. This review examines the emerging roles and molecular mechanisms of "nondomain biomolecules," whose functions are not easily predicted due to the absence of conserved functional domains. We propose the hypothesis that both domain- and nondomain-type molecules work together to enable flexible and efficient molecular processes within the highly crowded intracellular environment.

Solution structure and behaviour of the Arabidopsis thaliana HYL1 protein

In plants, microRNA biogenesis involves the complex assembly of molecular processes that are mostly governed by three proteins: RNase III protein DCL1 and two RNA binding proteins, SERRATE and HYL1. HYL1 protein is a double stranded RNA binding protein that is needed for the precise excision of miRNA/miRNA* duplex from the stem-loop containing primary miRNA gene transcripts. Moreover, HYL1 protein partners with HSP90 and CARP9 proteins to load the miRNA molecules onto the AGO1 endonuclease. HYL1 protein as a crucial player in the biogenesis pathway is regulated by its phosphorylation status to fine tune the levels of miRNA in various physiological conditions. HYL1 protein consists of two dsRNA binding domains (dsRBD) that are involved in RNA binding and dimerization and a C-terminal disordered tail of unknown function. Although the spatial structures of the individual dsRBDs have been determined there is a lack of information about the behaviour and structure of the full length protein. Using small the angle X-ray scattering (SAXS) technique we investigated the structure and dynamic of the HYL1 protein from Arabidopsis thaliana in solution. We show that the C-terminal domain is disordered and dynamic in solution and that HYL1 protein dimerization is dependent on the concentration. HYL1 protein lacking a C-terminal tail and a nuclear localisation signal (NLS) fragment is almost exclusively monomeric and similarly to full-length protein has a dynamic nature in solution. Our results point for the first time to the role of the C-terminal fragment in stabilisation of HYL1 dimer formation.

Regulation of Biomolecular Condensates by Poly(ADP-ribose)

Biomolecular condensates are reversible compartments that form through a process called phase separation. Post-translational modifications like ADP-ribosylation can nucleate the formation of these condensates by accelerating the self-association of proteins. Poly(ADP-ribose) (PAR) chains are remarkably transient modifications with turnover rates on the order of minutes, yet they can be required for the formation of granules in response to oxidative stress, DNA damage, and other stimuli. Moreover, accumulation of PAR is linked with adverse phase transitions in neurodegenerative diseases, including Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. In this review, we provide a primer on how PAR is synthesized and regulated, the diverse structures and chemistries of ADP-ribosylation modifications, and protein-PAR interactions. We review substantial progress in recent efforts to determine the molecular mechanism of PAR-mediated phase separation, and we further delineate how inhibitors of PAR polymerases may be effective treatments for neurodegenerative pathologies. Finally, we highlight the need for rigorous biochemical interrogation of ADP-ribosylation in vivo and in vitro to clarify the exact pathway from PARylation to condensate formation.