BIAPSS: A Comprehensive Physicochemical Analyzer of Proteins Undergoing Liquid–Liquid Phase Separation

The liquid–liquid phase separation (LLPS) of biomolecules is a phenomenon which is nowadays recognized as the driving force for the biogenesis of numerous functional membraneless organelles and cellular bodies. The interplay between the protein primary sequence and phase separation remains poorly understood, despite intensive research. To uncover the sequence-encoded signals of protein capable of undergoing LLPS, we developed a novel web platform named BIAPSS (Bioinformatics Analysis of LLPS Sequences). This web server provides on-the-fly analysis, visualization, and interpretation of the physicochemical and structural features for the superset of curated LLPS proteins.     

Liquid–Liquid Phase Separation in the Presence of Macromolecular Crowding and State-dependent Kinetics

Biomolecular condensates formed via liquid–liquid phase separation (LLPS) are increasingly being shown to play major roles in cellular self-organization dynamics in health and disease. It is well established that macromolecular crowding has a profound impact on protein interactions, particularly those that lead to LLPS. Although synthetic crowding agents are used during in vitro LLPS experiments, they are considerably different from the highly crowded nucleo-/cytoplasm and the effects of in vivo crowding remain poorly understood. In this work, we applied computational modeling to investigate the effects of macromolecular crowding on LLPS. To include biologically relevant LLPS dynamics, we extended the conventional Cahn–Hilliard model for phase separation by coupling it to experimentally derived macromolecular crowding dynamics and state-dependent reaction kinetics. Through extensive field-theoretic computer simulations, we show that the inclusion of macromolecular crowding results in late-stage coarsening and the stabilization of relatively smaller condensates. At a high crowding concentration, there is an accelerated growth and late-stage arrest of droplet formation, effectively resulting in anomalous labyrinthine morphologies akin to protein gelation observed in experiments. These results not only elucidate the crowder effects observed in experiments, but also highlight the importance of including state-dependent kinetics in LLPS models, and may help in designing further experiments to probe the intricate roles played by LLPS in self-organization dynamics of cells.     

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.     

Different Intermolecular Interactions Drive Nonpathogenic Liquid–Liquid Phase Separation and Potentially Pathogenic Fibril Formation by TDP-43

The liquid–liquid phase separation (LLPS) of proteins has been found ubiquitously in eukaryotic cells, and is critical in the control of many biological processes by forming a temporary condensed phase with different bimolecular components. TDP-43 is recruited to stress granules in cells and is the main component of TDP-43 granules and proteinaceous amyloid inclusions in patients with amyotrophic lateral sclerosis (ALS). TDP-43 low complexity domain (LCD) is able to de-mix in solution, forming the protein condensed droplets, and amyloid aggregates would form from the droplets after incubation. The molecular interactions regulating TDP-43 LCD LLPS were investigated at the protein fusion equilibrium stage, when the droplets stopped growing after incubation. We found the molecules in the droplet were still liquid-like, but with enhanced intermolecular helix–helix interactions. The protein would only start to aggregate after a lag time and aggregate slower than at the condition when the protein does not phase separately into the droplets, or the molecules have a reduced intermolecular helix–helix interaction. In the protein condensed droplets, a structural transition intermediate toward protein aggregation was discovered involving a decrease in the intermolecular helix–helix interaction and a reduction in the helicity. Our results therefore indicate that different intermolecular interactions drive LLPS and fibril formation. The discovery that TDP-43 LCD aggregation was faster through the pathway without the first protein phase separation supports that LLPS and the intermolecular helical interaction could help maintain the stability of TDP-43 LCD.     

Probing Protein Folding with Sequence-Reversed α-Helical Bundles

Recurrent protein folding motifs include various types of helical bundles formed by α-helices that supercoil around each other. While specific patterns of amino acid residues (heptad repeats) characterize the highly versatile folding motif of four-α-helical bundles, the significance of the polypeptide chain directionality is not sufficiently understood, although it determines sequence patterns, helical dipoles, and other parameters for the folding and oligomerization processes of bundles. To investigate directionality aspects in sequence-structure relationships, we reversed the amino acid sequences of two well-characterized, highly regular four-α-helical bundle proteins and studied the folding, oligomerization, and structural properties of the retro-proteins, using Circular Dichroism Spectroscopy (CD), Size Exclusion Chromatography combined with Multi-Angle Laser Light Scattering (SEC-MALS), and Small Angle X-ray Scattering (SAXS). The comparison of the parent proteins with their retro-counterparts reveals that while the α-helical character of the parents is affected to varying degrees by sequence reversal, the folding states, oligomerization propensities, structural stabilities, and shapes of the new molecules strongly depend on the characteristics of the heptad repeat patterns. The highest similarities between parent and retro-proteins are associated with the presence of uninterrupted heptad patterns in helical bundles sequences.     

Tau N-Terminal Inserts Regulate Tau Liquid-Liquid Phase Separation and Condensates Maturation in a Neuronal Cell Model

The microtubule-associated protein tau can undergo liquid–liquid phase separation (LLPS) to form membraneless condensates in neurons, yet the underlying molecular mechanisms and functions of tau LLPS and tau droplets remain to be elucidated. The human brain contains mainly 6 tau isoforms with different numbers of microtubule-binding repeats (3R, 4R) and N-terminal inserts (0N, 1N, 2N). However, little is known about the role of N-terminal inserts. Here we observed the dynamics of three tau isoforms with different N-terminal inserts in live neuronal cell line HT22. We validated tau LLPS in cytoplasm and found that 2N-tau forms liquid-like, hollow-shell droplets. Tau condensates became smaller in 1N-tau comparing with 2N-tau, while no obvious tau accumulated dots were shown in 0N-tau. The absence of N-terminal inserts significantly affected condensate colocalization of tau and p62. The results reveal insights into the tau LLPS assembly mechanism and functional effects of N-terminal inserts in tau.     

The Patterning and Proportion of Charged Residues in the Arginine-Rich Mixed-Charge Domain Determine the Membrane-Less Organelle Targeted by the Protein

Membrane-less organelles (MLOs) are formed by biomolecular liquid–liquid phase separation (LLPS). Proteins with charged low-complexity domains (LCDs) are prone to phase separation and localize to MLOs, but the mechanism underlying the distributions of such proteins to specific MLOs remains poorly understood. Recently, proteins with Arg-enriched mixed-charge domains (R-MCDs), primarily composed of R and Asp (D), were found to accumulate in nuclear speckles via LLPS. However, the process by which R-MCDs selectively incorporate into nuclear speckles is unknown. Here, we demonstrate that the patterning of charged amino acids and net charge determines the targeting of specific MLOs, including nuclear speckles and the nucleolus, by proteins. The redistribution of R and D residues from an alternately sequenced pattern to uneven blocky sequences caused a shift in R-MCD distribution from nuclear speckles to the nucleolus. In addition, the incorporation of basic residues in the R-MCDs promoted their localization to the MLOs and their apparent accumulation in the nucleolus. The R-MCD peptide with alternating amino acids did not undergo LLPS, whereas the blocky R-MCD peptide underwent LLPS with affinity to RNA, acidic poly-Glu, and the acidic nucleolar protein nucleophosmin, suggesting that the clustering of R residues helps avoid their neutralization by D residues and eventually induces R-MCD migration to the nucleolus. Therefore, the distribution of proteins to nuclear speckles requires the proximal positioning of D and R for the mutual neutralization of their charges.     

Pressure and Chemical Unfolding of an α-Helical Bundle Protein: The GH2 Domain of the Protein Adaptor GIPC1

When combined with NMR spectroscopy, high hydrostatic pressure is an alternative perturbation method used to destabilize globular proteins that has proven to be particularly well suited for exploring the unfolding energy landscape of small single-domain proteins. To date, investigations of the unfolding landscape of all-β or mixed-α/β protein scaffolds are well documented, whereas such data are lacking for all-α protein domains. Here we report the NMR study of the unfolding pathways of GIPC1-GH2, a small α-helical bundle domain made of four antiparallel α-helices. High-pressure perturbation was combined with NMR spectroscopy to unravel the unfolding landscape at three different temperatures. The results were compared to those obtained from classical chemical denaturation. Whatever the perturbation used, the loss of secondary and tertiary contacts within the protein scaffold is almost simultaneous. The unfolding transition appeared very cooperative when using high pressure at high temperature, as was the case for chemical denaturation, whereas it was found more progressive at low temperature, suggesting the existence of a complex folding pathway.     

Thermodynamics of π–π Interactions of Benzene and Phenol in Water

The π–π interaction is a major driving force that stabilizes protein assemblies during protein folding. Recent studies have additionally demonstrated its involvement in the liquid–liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs). As the participating residues in IDPs are exposed to water, π–π interactions for LLPS must be modeled in water, as opposed to the interactions that are often established at the hydrophobic domains of folded proteins. Thus, we investigated the association of free energies of benzene and phenol dimers in water by integrating van der Waals (vdW)-corrected density functional theory (DFT) and DFT in classical explicit solvents (DFT-CES). By comparing the vdW-corrected DFT and DFT-CES results with high-level wavefunction calculations and experimental solvation free energies, respectively, we established the quantitative credibility of these approaches, enabling a reliable prediction of the benzene and phenol dimer association free energies in water. We discovered that solvation influences dimer association free energies, but not significantly when no direct hydrogen-bond-type interaction exists between two monomeric units, which can be explained by the enthalpy–entropy compensation. Our comprehensive computational study of the solvation effect on π–π interactions in water could help us understand the molecular-level driving mechanism underlying the IDP phase behaviors.     

RNA-Binding Protein MEX3A Interacting with DVL3 Stabilizes Wnt/β-Catenin Signaling in Endometrial Carcinoma

Disease recurrence and metastasis lead to poor prognosis in patients with advanced endometrial carcinoma (EC). RNA-binding proteins (RBPs) are closely associated with tumor initiation and metastasis, but the function and molecular mechanisms of RBPs in EC are unclear. RBPs were screened and identified using the TCGA, GEO, and RBPTD databases. The effect of MEX3A on EC was verified by in vitro and in vivo experiments. Gene set enrichment analysis (GSEA), immunofluorescence (IF), and co-immunoprecipitation (Co-IP) were used to identify potential molecular mechanisms of action. We identified 148 differentially expressed RBPs in EC. MEX3A was upregulated and related to poor prognosis in patients with EC. In vitro and vivo experiments demonstrated that MEX3A promoted the growth, migration, and invasion capacities of EC cells. Mechanistically, DVL3, a positive regulator of the Wnt/β-catenin pathway, also increased the proliferation and metastasis of EC cells. MEX3A enhanced EMT and played a pro-carcinogenic role by interacting with DVL3 to stabilize β-catenin and upregulated the expression of its downstream target genes. MEX3A is upregulated in EC and promotes tumor progression by activating EMT and regulating the Wnt/β-catenin pathway via DVL3. MEX3A may therefore be a novel therapeutic target for EC.     

The Role of Natural Polymorphic Variants of DNA Polymerase β in DNA Repair

DNA polymerase β (Polβ) is considered the main repair DNA polymerase involved in the base excision repair (BER) pathway, which plays an important part in the repair of damaged DNA bases usually resulting from alkylation or oxidation. In general, BER involves consecutive actions of DNA glycosylases, AP endonucleases, DNA polymerases, and DNA ligases. It is known that protein–protein interactions of Polβ with enzymes from the BER pathway increase the efficiency of damaged base repair in DNA. However natural single-nucleotide polymorphisms can lead to a substitution of functionally significant amino acid residues and therefore affect the catalytic activity of the enzyme and the accuracy of Polβ action. Up-to-date databases contain information about more than 8000 SNPs in the gene of Polβ. This review summarizes data on the in silico prediction of the effects of Polβ SNPs on DNA repair efficacy; available data on cancers associated with SNPs of Polβ; and experimentally tested variants of Polβ. Analysis of the literature indicates that amino acid substitutions could be important for the maintenance of the native structure of Polβ and contacts with DNA; others affect the catalytic activity of the enzyme or play a part in the precise and correct attachment of the required nucleotide triphosphate. Moreover, the amino acid substitutions in Polβ can disturb interactions with enzymes involved in BER, while the enzymatic activity of the polymorphic variant may not differ significantly from that of the wild-type enzyme. Therefore, investigation regarding the effect of Polβ natural variants occurring in the human population on enzymatic activity and protein–protein interactions is an urgent scientific task.     

miR-29b Regulates TGF-β1-Induced Epithelial–Mesenchymal Transition by Inhibiting Heat Shock Protein 47 Expression in Airway Epithelial Cells

Tissue remodeling contributes to ongoing inflammation and refractoriness of chronic rhinosinusitis (CRS). During this process, epithelial-mesenchymal transition (EMT) plays an important role in dysregulated remodeling and both microRNA (miR)-29b and heat shock protein 47 (HSP47) may be engaged in the pathophysiology of CRS. This study aimed to determine the role of miR-29b and HSP47 in modulating transforming growth factor (TGF)-β1-induced EMT and migration in airway epithelial cells. Expression levels of miR-29b, HSP47, E-cadherin, α-smooth muscle actin (α-SMA), vimentin and fibronectin were assessed through real-time PCR, Western blotting, and immunofluorescence staining. Small interfering RNA (siRNA) targeted against miR-29b and HSP47 were transfected to regulate the expression of EMT-related markers. Cell migration was evaluated with wound scratch and transwell migration assay. miR-29b mimic significantly inhibited the expression of HSP47 and TGF-β1-induced EMT-related markers in A549 cells. However, the miR-29b inhibitor more greatly induced the expression of them. HSP47 knockout suppressed TGF-β1-induced EMT marker levels. Functional studies indicated that TGF-β1-induced EMT was regulated by miR-29b and HSP47 in A549 cells. These findings were further verified in primary nasal epithelial cells. miR-29b modulated TGF-β1-induced EMT-related markers and migration via HSP47 expression modulation in A549 and primary nasal epithelial cells. These results suggested the importance of miR-29b and HSP47 in pathologic tissue remodeling progression in CRS.     

Indomethacin Disrupts the Formation of β-Amyloid Plaques via an α2-Macroglobulin-Activating lrp1-Dependent Mechanism

Epidemiological studies have implied that the nonsteroidal anti-inflammatory drug (NSAID) indomethacin slows the development and progression of Alzheimer’s disease (AD). However, the underlying mechanisms are notably understudied. Using a chimeric mouse/human amyloid precursor protein (Mo/HuAPP695swe) and a mutant human presenilin 1 (PS1-dE9) (APP/PS1) expressing transgenic (Tg) mice and neuroblastoma (N) 2a cells as in vivo and in vitro models, we revealed the mechanisms of indomethacin in ameliorating the cognitive decline of AD. By screening AD-associated genes, we observed that a marked increase in the expression of α₂-macroglobulin (A2M) was markedly induced after treatment with indomethacin. Mechanistically, upregulation of A2M was caused by the inhibition of cyclooxygenase-2 (COX-2) and lipocalin-type prostaglandin D synthase (L-PGDS), which are responsible for the synthesis of prostaglandin (PG)H₂ and PGD₂, respectively. The reduction in PGD₂ levels induced by indomethacin alleviated the suppression of A2M expression through a PGD₂ receptor 2 (CRTH2)-dependent mechanism. Highly activated A2M not only disrupted the production and aggregation of β-amyloid protein (Aβ) but also induced Aβ efflux from the brain. More interestingly, indomethacin decreased the degradation of the A2M receptor, low-density lipoprotein receptor-related protein 1 (LRP1), which facilitated the brain efflux of Aβ. Through the aforementioned mechanisms, indomethacin ameliorated cognitive decline in APP/PS1 Tg mice by decreasing Aβ production and clearing Aβ from the brains of AD mice.     

Various Approaches to Studying the Phase Transition in an Octamethylcyclotetrasiloxane Crystal: From X-ray Structural Analysis to Metadynamics

The structure, thermodynamic parameters, and the character of thermal motion in octamethylcyclotetrasiloxane (D4) were investigated using the combination of experimental (single-crystal X-ray diffraction, thermochemistry) and theoretical (density functional theory calculations, ab initio molecular dynamics and metadynamics) methods. Single crystals of D4 were grown in a glass capillary in situ and the structures of high- (238–270 K) and low-temperature (100–230 K) phases were studied in detail. In the temperature interval 230–238 K, a phase transition with rather low enthalpy (−1.04(7) kcal/mol) was detected. It was found that phase transition is accompanied by change of conformation of cyclosiloxane moiety from boat-saddle (cradle) to chair. According to PBE0/6-311G(d,p) calculation of isolated D4, such conformation changes are characterized by a low barrier (0.07 kcal/mol). The character of molecular thermal motion and the path of phase transition were established with combination of periodic DFT calculations, including molecular dynamics and metadynamics. The effect of crystal field led to an increase in the calculated phase transition barrier (4.27 kcal/mol from low- to high-temperature phase and 3.20 kcal/mol in opposite direction).     

The Multivalent Polyampholyte Domain of Nst1, a P-Body-Associated Saccharomyces cerevisiae Protein,Provides a Platform for Interacting with P-Body Components

The condensation of nuclear promyelocytic leukemia bodies, cytoplasmic P-granules, P-bodies (PBs), and stress granules is reversible and dynamic via liquid–liquid phase separation. Although each condensate comprises hundreds of proteins with promiscuous interactions, a few key scaffold proteins are required. Essential scaffold domain sequence elements, such as poly-Q, low-complexity regions, oligomerizing domains, and RNA-binding domains, have been evaluated to understand their roles in biomolecular condensation processes. However, the underlying mechanisms remain unclear. We analyzed Nst1, a PB-associated protein that can intrinsically induce PB component condensations when overexpressed. Various Nst1 domain deletion mutants with unique sequence distributions, including intrinsically disordered regions (IDRs) and aggregation-prone regions, were constructed based on structural predictions. The overexpression of Nst1 deletion mutants lacking the aggregation-prone domain (APD) significantly inhibited self-condensation, implicating APD as an oligomerizing domain promoting self-condensation. Remarkably, cells overexpressing the Nst1 deletion mutant of the polyampholyte domain (PD) in the IDR region (Nst1_∆PD) rarely accumulate endogenous enhanced green fluorescent protein (EGFP)-tagged Dcp2. However, Nst1_∆PD formed self-condensates, suggesting that Nst1 requires PD to interact with Dcp2, regardless of its self-condensation. In Nst1_∆PD-overexpressing cells treated with cycloheximide (CHX), Dcp2, Xrn1, Dhh1, and Edc3 had significantly diminished condensation compared to those in CHX-treated Nst1-overexpressing cells. These observations suggest that the PD of the IDR in Nst1 functions as a hub domain interacting with other PB components.     

Seeing Keratinocyte Proteins through the Looking Glass of Intrinsic Disorder

Epidermal keratinocyte proteins include many with an eccentric amino acid content (compositional bias), atypical ultrastructural fate (built-in protease sensitivity), or assembly visible at the light microscope level (cytoplasmic granules). However, when considered through the looking glass of intrinsic disorder (ID), these apparent oddities seem quite expected. Keratinocyte proteins with highly repetitive motifs are of low complexity but high adaptation, providing polymers (e.g., profilaggrin) for proteolysis into bioactive derivatives, or monomers (e.g., loricrin) repeatedly cross-linked to self and other proteins to shield underlying tissue. Keratohyalin granules developing from liquid–liquid phase separation (LLPS) show that unique biomolecular condensates (BMC) and proteinaceous membraneless organelles (PMLO) occur in these highly customized cells. We conducted bioinformatic and in silico assessments of representative keratinocyte differentiation-dependent proteins. This was conducted in the context of them having demonstrated potential ID with the prospect of that characteristic driving formation of distinctive keratinocyte structures. Intriguingly, while ID is characteristic of many of these proteins, it does not appear to guarantee LLPS, nor is it required for incorporation into certain keratinocyte protein condensates. Further examination of keratinocyte-specific proteins will provide variations in the theme of PMLO, possibly recognizing new BMC for advancements in understanding intrinsically disordered proteins as reflected by keratinocyte biology.     

Highly Diverse Protein Library Based on the Ubiquitous (β/α)8 Enzyme Fold Yields Well‐Structured Proteins through in Vitro Folding Selection

Proper protein folding is a prerequisite for protein stability and enzymatic activity. Although directed evolution can be a powerful tool to investigate enzymatic function and to isolate novel activities, well‐designed libraries of folded proteins are essential. In vitro selection methods are particularly capable of searching for enzymatic activities in libraries of trillions of protein variants, yet high‐quality libraries of well‐folded enzymes with such high diversity are lacking. We describe the construction and detailed characterization of a folding‐enriched protein library based on the ubiquitous (β/α)₈ barrel fold, which is found in five of the six enzyme classes. We introduced seven randomized loops on the catalytic face of the monomeric, thermostable (β/α)₈ barrel of glycerophosphodiester phosphodiesterase (GDPD) from Thermotoga maritima. We employed in vitro folding selection based on protease digestion to enrich intermediate libraries containing three to four randomized loops for folded variants, and then combined them to assemble the final library (10¹⁴ DNA sequences). The resulting library was analyzed by using the in vitro protease assay and an in vivo GFP‐folding assay; it contains ～10¹² soluble monomeric protein variants. We isolated six library members and demonstrated that these proteins are soluble, monomeric and show (β/α)₈‐barrel fold‐like secondary and tertiary structure. The quality of the folding‐enriched library improved up to 50‐fold compared to a control library that was assembled without the folding selection. To the best of our knowledge, this work is the first example of combining the ultra‐high throughput mRNA display method with selection for folding. The resulting (β/α)₈ barrel libraries provide a valuable starting point to study the unique catalytic capabilities of the (β/α)₈ fold, and to isolate novel enzymes.     

Split Enzyme-Based Biosensors for Structural Characterization of Soluble and Insoluble β-Glucans

β-Glucan is widely distributed in various plants and microorganisms and is composed of β-1,3-linked d-glucose units. It may have a branched short or long side chain of glucose units with β-1,6- or β-1,4-linkage. Numerous studies have investigated different β-glucans and revealed their bioactivities. To understand the structure-function relationship of β-glucan, we constructed a split-luciferase complementation assay for the structural analysis of long-chain β-1,6-branched β-1,3-glucan. The N- and C-terminal fragments of luciferase from deep-sea shrimp were fused to insect-derived β-1,3-glucan recognition protein and fungal endo-β-1,6-glucanase (Neg1)-derived β-1,6-glucan recognition protein, respectively. In this approach, two β-glucan recognition proteins bound to β-glucan molecules come into close proximity, resulting in the assembly of the full-length reporter enzyme and induction of transient luciferase activity, indicative of the structure of β-glucan. To test the applicability of this assay, β-glucan and two β-glucan recognition proteins were mixed, resulting in an increase in the luminescence intensity in a β-1,3-glucan with a long polymer of β-1,6-glucan in a dose-dependent manner. This simple test also allows the monitoring of real-time changes in the side chain structure and serves as a convenient method to distinguish between β-1,3-glucan and long-chain β-1,6-branched β-1,3-glucan in various soluble and insoluble β-glucans.