Adenosine triphosphate energy‐independently controls protein homeostasis with unique structure and diverse mechanisms

Proteins function in the crowded cellular environments with high salt concentrations, thus facing tremendous challenges of misfolding/aggregation which represents a pathological hallmark of aging and an increasing spectrum of human diseases. Recently, intrinsically disordered regions (IDRs) were recognized to drive liquid–liquid phase separation (LLPS), a common principle for organizing cellular membraneless organelles (MLOs). ATP, the universal energy currency for all living cells, mysteriously has concentrations of 2–12 mM, much higher than required for its previously‐known functions. Only recently, ATP was decoded to behave as a biological hydrotrope to inhibit protein LLPS and aggregation at mM. We further revealed that ATP also acts as a bivalent binder, which not only biphasically modulates LLPS driven by IDRs of human and viral proteins, but also bind to the conserved nucleic‐acid‐binding surfaces of the folded proteins. Most unexpectedly, ATP appears to act as a hydration mediator to antagonize the crowding‐induced destabilization as well as to enhance folding of proteins without significant binding. Here, this review focuses on summarizing the results of these biophysical studies and discussing their implications in an evolutionary context. By linking triphosphate with unique hydration property to adenosine, ATP appears to couple the ability for establishing hydrophobic, π‐π, π‐cation and electrostatic interactions to the capacity in mediating hydration of proteins, which is at the heart of folding, dynamics, stability, phase separation and aggregation. Consequently, ATP acquired a category of functions at ~mM to energy‐independently control protein homeostasis with diverse mechanisms, thus implying a link between cellular ATP concentrations and protein‐aggregation diseases.     

CTD of SARS‐CoV‐2 N protein is a cryptic domain for binding ATP and nucleic acid that interplay in modulating phase separation

SARS‐CoV‐2 nucleocapsid (N) protein plays essential roles in many steps of the viral life cycle, thus representing a key drug target. N protein contains the folded N‐/C‐terminal domains (NTD/CTD) and three intrinsically disordered regions, while its functions including liquid–liquid phase separation (LLPS) depend on the capacity in binding various viral/host‐cell RNA/DNA of diverse sequences. Previously NTD was established to bind various RNA/DNA while CTD to dimerize/oligomerize for forming high‐order structures. By NMR, here for the first time we decrypt that CTD is not only capable of binding S2m, a specific probe derived from SARS‐CoV‐2 gRNA but with the affinity even higher than that of NTD. Very unexpectedly, ATP, the universal energy currency for all living cells with high cellular concentrations (2–16 mM), specifically binds CTD with Kd of 1.49 ± 0.28 mM. Strikingly, the ATP‐binding residues of NTD/CTD are identical in the SARS‐CoV‐2 variants while ATP and S2m interplay in binding NTD/CTD, as well as in modulating LLPS critical for the viral life cycle. Results together not only define CTD as a novel binding domain for ATP and nucleic acid, but enforce our previous proposal that ATP has been evolutionarily exploited by SARS‐CoV‐2 to complete its life cycle in the host cell. Most importantly, the unique ATP‐binding pockets on NTD/CTD may offer promising targets for design of specific anti‐SARS‐CoV‐2 molecules to fight the pandemic. Fundamentally, ATP emerges to act at mM as a cellular factor to control the interface between the host cell and virus lacking the ability to generate ATP.     

RGG/RG Motif Regions in RNA Binding and Phase Separation

RGG/RG motifs are RNA binding segments found in many proteins that can partition into membraneless organelles. They occur in the context of low-complexity disordered regions and often in multiple copies. Although short RGG/RG-containing regions can sometimes form high-affinity interactions with RNA structures, multiple RGG/RG repeats are generally required for high-affinity binding, suggestive of the dynamic, multivalent interactions that are thought to underlie phase separation in formation of cellular membraneless organelles. Arginine can interact with nucleotide bases via hydrogen bonding and π-stacking; thus, nucleotide conformers that provide access to the bases provide enhanced opportunities for RGG interactions. Methylation of RGG/RG regions, which is accomplished by protein arginine methyltransferase enzymes, occurs to different degrees in different cell types and may regulate the behavior of proteins containing these regions.     

SARS‐CoV‐2 nucleocapsid protein phase‐separates with RNA and with human hnRNPs

Tightly packed complexes of nucleocapsid protein and genomic RNA form the core of viruses and assemble within viral factories, dynamic compartments formed within the host cells associated with human stress granules. Here, we test the possibility that the multivalent RNA‐binding nucleocapsid protein (N) from severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) condenses with RNA via liquid–liquid phase separation (LLPS) and that N protein can be recruited in phase‐separated forms of human RNA‐binding proteins associated with SG formation. Robust LLPS with RNA requires two intrinsically disordered regions (IDRs), the N‐terminal IDR and central‐linker IDR, as well as the folded C‐terminal oligomerization domain, while the folded N‐terminal domain and the C‐terminal IDR are not required. N protein phase separation is induced by addition of non‐specific RNA. In addition, N partitions in vitro into phase‐separated forms of full‐length human hnRNPs (TDP‐43, FUS, hnRNPA2) and their low‐complexity domains (LCs). These results provide a potential mechanism for the role of N in SARS‐CoV‐2 viral genome packing and in host‐protein co‐opting necessary for viral replication and infectivity.     

Structural mechanism for modulation of functional amyloid and biofilm formation by Staphylococcal Bap protein switch

The Staphylococcal Bap proteins sense environmental signals (such as pH, [Ca²⁺]) to build amyloid scaffold biofilm matrices via unknown mechanisms. We here report the crystal structure of the aggregation‐prone region of Staphylococcus aureus Bap which adopts a dumbbell‐shaped fold. The middle module (MM) connecting the N‐terminal and C‐terminal lobes consists of a tandem of novel double‐Ca²⁺‐binding motifs involved in cooperative interaction networks, which undergoes Ca²⁺‐dependent order–disorder conformational switches. The N‐terminal lobe is sufficient to mediate amyloid aggregation through liquid–liquid phase separation and maturation, and subsequent biofilm formation under acidic conditions. Such processes are promoted by disordered MM at low [Ca²⁺] but inhibited by ordered MM stabilized by Ca²⁺ binding, with inhibition efficiency depending on structural integrity of the interaction networks. These studies illustrate a novel protein switch in pathogenic bacteria and provide insights into the mechanistic understanding of Bap proteins in modulation of functional amyloid and biofilm formation, which could be implemented in the anti‐biofilm drug design.     

Frontotemporal dementia‐linked P112H mutation of TDP‐43 induces protein structural change and impairs its RNA binding function

TDP‐43 forms the primary constituents of the cytoplasmic inclusions contributing to various neurodegenerative diseases, including amyotrophic lateral sclerosis and frontotemporal dementia (FTD). Over 60 TDP‐43 mutations have been identified in patients suffering from these two diseases, but most variations are located in the protein''s disordered C‐terminal glycine‐rich region. P112H mutation of TDP‐43 has been uniquely linked to FTD, and is located in the first RNA recognition motif (RRM1). This mutation is thought to be pathogenic, but its impact on TDP‐43 at the protein level remains unclear. Here, we compare the biochemical and biophysical properties of TDP‐43 truncated proteins with or without P112H mutation. We show that P112H‐mutated TDP‐43 proteins exhibit higher thermal stability, impaired RNA‐binding activity, and a reduced tendency to aggregate relative to wild‐type proteins. Near‐UV CD, 2D‐nuclear‐magnetic resonance, and intrinsic fluorescence spectrometry further reveal that the P112H mutation in RRM1 generates local conformational changes surrounding the mutational site that disrupt the stacking interactions of the W113 side chain with nucleic acids. Together, these results support the notion that P112H mutation of TDP‐43 contributes to FTD through functional impairment of RNA metabolism and/or structural changes that curtail protein clearance.     

“Protein aggregates” contain RNA and DNA,entrapped by misfolded proteins but largely rescued by slowing translational elongation

All neurodegenerative diseases feature aggregates, which usually contain disease‐specific diagnostic proteins; non‐protein constituents, however, have rarely been explored. Aggregates from SY5Y‐APP_Sw neuroblastoma, a cell model of familial Alzheimer''s disease, were crosslinked and sequences of linked peptides identified. We constructed a normalized “contactome” comprising 11 subnetworks, centered on 24 high‐connectivity hubs. Remarkably, all 24 are nucleic acid‐binding proteins. This led us to isolate and sequence RNA and DNA from Alzheimer''s and control aggregates. RNA fragments were mapped to the human genome by RNA‐seq and DNA by ChIP‐seq. Nearly all aggregate RNA sequences mapped to specific genes, whereas DNA fragments were predominantly intergenic. These nucleic acid mappings are all significantly nonrandom, making an artifactual origin extremely unlikely. RNA (mostly cytoplasmic) exceeded DNA (chiefly nuclear) by twofold to fivefold. RNA fragments recovered from AD tissue were ~1.5‐to 2.5‐fold more abundant than those recovered from control tissue, similar to the increase in protein. Aggregate abundances of specific RNA sequences were strikingly differential between cultured SY5Y‐APP_Sw glioblastoma cells expressing APOE3 vs. APOE4, consistent with APOE4 competition for E‐box/CLEAR motifs. We identified many G‐quadruplex and viral sequences within RNA and DNA of aggregates, suggesting that sequestration of viral genomes may have driven the evolution of disordered nucleic acid‐binding proteins. After RNA‐interference knockdown of the translational‐procession factor EEF2 to suppress translation in SY5Y‐APP_Sw cells, the RNA content of aggregates declined by >90%, while reducing protein content by only 30% and altering DNA content by ≤10%. This implies that cotranslational misfolding of nascent proteins may ensnare polysomes into aggregates, accounting for most of their RNA content.     

DEAD‐box polypeptide 43 facilitates piRNA amplification by actively liberating RNA from Ago3‐piRISC

The piRNA amplification pathway in Bombyx is operated by Ago3 and Siwi in their piRISC form. The DEAD‐box protein, Vasa, facilitates Ago3‐piRISC production by liberating cleaved RNAs from Siwi‐piRISC in an ATP hydrolysis‐dependent manner. However, the Vasa‐like factor facilitating Siwi‐piRISC production along this pathway remains unknown. Here, we identify DEAD‐box polypeptide 43 (DDX43) as the Vasa‐like protein functioning in Siwi‐piRISC production. DDX43 belongs to the helicase superfamily II along with Vasa, and it contains a similar helicase core. DDX43 also contains a K‐homology (KH) domain, a prevalent RNA‐binding domain, within its N‐terminal region. Biochemical analyses show that the helicase core is responsible for Ago3‐piRISC interaction and ATP hydrolysis, while the KH domain enhances the ATPase activity of the helicase core. This enhancement is independent of the RNA‐binding activity of the KH domain. For maximal DDX43 RNA‐binding activity, both the KH domain and helicase core are required. This study not only provides new insight into the piRNA amplification mechanism but also reveals unique collaborations between the two domains supporting DDX43 function within the pathway.     

Electrostatic modulation of hnRNPA1 low‐complexity domain liquid–liquid phase separation and aggregation

Membrane‐less organelles and RNP granules are enriched in RNA and RNA‐binding proteins containing disordered regions. Heterogeneous nuclear ribonucleoprotein A1 (hnRNPA1), a key regulating protein in RNA metabolism, localizes to cytoplasmic RNP granules including stress granules. Dysfunctional nuclear‐cytoplasmic transport and dynamic phase separation of hnRNPA1 leads to abnormal amyloid aggregation and neurodegeneration. The intrinsically disordered C‐terminal domain (CTD) of hnRNPA1 mediates both dynamic liquid–liquid phase separation (LLPS) and aggregation. While cellular phase separation drives the formation of membrane‐less organelles, aggregation within phase‐separated compartments has been linked to neurodegenerative diseases. To understand some of the underlying mechanisms behind protein phase separation and LLPS‐mediated aggregation, we studied LLPS of hnRNPA1 CTD in conditions that probe protein electrostatics, modulated specifically by varying pH conditions, and protein, salt and RNA concentrations. In the conditions investigated, we observed LLPS to be favored in acidic conditions, and by high protein, salt and RNA concentrations. We also observed that conditions that favor LLPS also enhance protein aggregation and fibrillation, which suggests an aggregation pathway that is LLPS‐mediated. The results reported here also suggest that LLPS can play a direct role in facilitating protein aggregation, and that changes in cellular environment that affect protein electrostatics can contribute to the pathological aggregation exhibited in neurodegeneration.     

Systematic discovery of linear binding motifs targeting an ancient protein interaction surface on MAP kinases

Mitogen‐activated protein kinases (MAPK) are broadly used regulators of cellular signaling. However, how these enzymes can be involved in such a broad spectrum of physiological functions is not understood. Systematic discovery of MAPK networks both experimentally and in silico has been hindered because MAPKs bind to other proteins with low affinity and mostly in less‐characterized disordered regions. We used a structurally consistent model on kinase‐docking motif interactions to facilitate the discovery of short functional sites in the structurally flexible and functionally under‐explored part of the human proteome and applied experimental tools specifically tailored to detect low‐affinity protein–protein interactions for their validation in vitro and in cell‐based assays. The combined computational and experimental approach enabled the identification of many novel MAPK‐docking motifs that were elusive for other large‐scale protein–protein interaction screens. The analysis produced an extensive list of independently evolved linear binding motifs from a functionally diverse set of proteins. These all target, with characteristic binding specificity, an ancient protein interaction surface on evolutionarily related but physiologically clearly distinct three MAPKs (JNK, ERK, and p38). This inventory of human protein kinase binding sites was compared with that of other organisms to examine how kinase‐mediated partnerships evolved over time. The analysis suggests that most human MAPK‐binding motifs are surprisingly new evolutionarily inventions and newly found links highlight (previously hidden) roles of MAPKs. We propose that short MAPK‐binding stretches are created in disordered protein segments through a variety of ways and they represent a major resource for ancient signaling enzymes to acquire new regulatory roles.     

The free energy folding penalty accompanying binding of intrinsically disordered α‐helical motifs

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes and preform critical roles in many cellular processes, most often through the association with globular proteins. Despite lacking a stable three‐dimensional structure by themselves, they may acquire a defined conformation upon binding globular targets. The most common type of secondary structure acquired by these binding motifs entails formation of an α‐helix. It has been hypothesized that such disorder‐to‐order transitions are associated with a significant free energy penalty due to IDP folding, which reduces the overall IDP‐target affinity. However, the exact magnitude of IDP folding penalty in α‐helical binding motifs has not been systematically estimated. Here, we report the folding penalty contributions for 30 IDPs undergoing folding‐upon‐binding and find that the average IDP folding penalty is +2.0 kcal/mol and ranges from 0.7 to 3.5 kcal/mol. We observe that the folding penalty scales approximately linearly with the change in IDP helicity upon binding, which provides a simple empirical way to estimate folding penalty. We analyze to what extent do pre‐structuring and target‐bound IDP dynamics (fuzziness) reduce the folding penalty and find that these effects combined, on average, reduce the folding cost by around half. Taken together, the presented analysis provides a quantitative basis for understanding the role of folding penalty in IDP‐target interactions and introduces a method estimate this quantity. Estimation and reduction of IDP folding penalty may prove useful in the rational design of helix‐stabilized inhibitors of IDP‐target interactions.StatementThe α‐helical binding motifs are ubiquitous among the intrinsically disordered proteins (IDPs). Upon binding their targets, they undergo a disorder‐to‐order transition, which is accompanied by a significant folding penalty whose magnitude is generally not known. Here, we use recently developed statistical‐thermodynamic model to estimate the folding penalties for 30 IDPs and clarify the roles of IDP pre‐folding and bound‐state dynamics in reducing the folding penalty.     

NONO enhances mRNA processing of super‐enhancer‐associated GATA2 and HAND2 genes in neuroblastoma

High‐risk neuroblastoma patients have poor survival rates and require better therapeutic options. High expression of a multifunctional DNA and RNA‐binding protein, NONO, in neuroblastoma is associated with poor patient outcome; however, there is little understanding of the mechanism of NONO‐dependent oncogenic gene regulatory activity in neuroblastoma. Here, we used cell imaging, biochemical and genome‐wide molecular analysis to reveal complex NONO‐dependent regulation of gene expression. NONO forms RNA‐ and DNA‐tethered condensates throughout the nucleus and undergoes phase separation in vitro, modulated by nucleic acid binding. CLIP analyses show that NONO mainly binds to the 5′ end of pre‐mRNAs and modulates pre‐mRNA processing, dependent on its RNA‐binding activity. NONO regulates super‐enhancer‐associated genes, including HAND2 and GATA2. Abrogating NONO RNA binding, or phase separation activity, results in decreased expression of HAND2 and GATA2. Thus, future development of agents that target RNA‐binding activity of NONO may have therapeutic potential in this cancer context.     

Disease‐linked TDP‐43 hyperphosphorylation suppresses TDP‐43 condensation and aggregation

Post‐translational modifications (PTMs) have emerged as key modulators of protein phase separation and have been linked to protein aggregation in neurodegenerative disorders. The major aggregating protein in amyotrophic lateral sclerosis and frontotemporal dementia, the RNA‐binding protein TAR DNA‐binding protein (TDP‐43), is hyperphosphorylated in disease on several C‐terminal serine residues, a process generally believed to promote TDP‐43 aggregation. Here, we however find that Casein kinase 1δ‐mediated TDP‐43 hyperphosphorylation or C‐terminal phosphomimetic mutations reduce TDP‐43 phase separation and aggregation, and instead render TDP‐43 condensates more liquid‐like and dynamic. Multi‐scale molecular dynamics simulations reveal reduced homotypic interactions of TDP‐43 low‐complexity domains through enhanced solvation of phosphomimetic residues. Cellular experiments show that phosphomimetic substitutions do not affect nuclear import or RNA regulatory functions of TDP‐43, but suppress accumulation of TDP‐43 in membrane‐less organelles and promote its solubility in neurons. We speculate that TDP‐43 hyperphosphorylation may be a protective cellular response to counteract TDP‐43 aggregation.     

Thermodynamic coupling between neighboring binding sites in homo‐oligomeric ligand sensing proteins from mass resolved ligand‐dependent population distributions

Homo‐oligomeric ligand‐activated proteins are ubiquitous in biology. The functions of such molecules are commonly regulated by allosteric coupling between ligand‐binding sites. Understanding the basis for this regulation requires both quantifying the free energy ΔG transduced between sites, and the structural basis by which it is transduced. We consider allostery in three variants of the model ring‐shaped homo‐oligomeric trp RNA‐binding attenuation protein (TRAP). First, we developed a nearest‐neighbor statistical thermodynamic binding model comprising microscopic free energies for ligand binding to isolated sites ΔG ₀, and for coupling between adjacent sites, ΔG _α . Using the resulting partition function (PF) we explored the effects of these parameters on simulated population distributions for the 2 ^N possible liganded states. We then experimentally monitored ligand‐dependent population shifts using conventional spectroscopic and calorimetric methods and using native mass spectrometry (MS). By resolving species with differing numbers of bound ligands by their mass, native MS revealed striking differences in their ligand‐dependent population shifts. Fitting the populations to a binding polynomial derived from the PF yielded coupling free energy terms corresponding to orders of magnitude differences in cooperativity. Uniquely, this approach predicts which of the possible 2 ^N liganded states are populated at different ligand concentrations, providing necessary insights into regulation. The combination of statistical thermodynamic modeling with native MS may provide the thermodynamic foundation for a meaningful understanding of the structure–thermodynamic linkage that drives cooperativity.