General Characterization of Properties of Ordered and Disordered Proteins by Wide-Line 1H NMR.

Wide-line 1H NMR is an efficient spectroscopic method to determine the disorder tendency of a protein. It directly measures the properties of the hydration shell of proteins, delivering exact and measurable values of their disorder/order content. A comparison is performed between several globular and disordered proteins. The common properties of the subzero mobile hydration water of these two groups were investigated. The amount of the mobile hydration water and the shape of the melting diagram at subzero temperatures together provide a possibility to distinguish globular proteins from disordered proteins. The shape of the melting diagram also gives information about the presence of secondary structural elements. The disordered and globular protein regions' fundamentally different structures are reflected in their melting diagrams, allowing one to directly determine the level of disorder in a specific protein structure. Intrinsically disordered proteins bind water more strongly than globular proteins, which is shown by the somewhat higher temperature values where mobile hydration water first appears but with a significantly higher heterogeneity in the energy distributions of protein-water interactions.

Morphological Changes Induced by TKS4 Deficiency Can Be Reversed by EZH2 Inhibition in Colorectal Carcinoma Cells.

BACKGROUND: The scaffold protein tyrosine kinase substrate 4 (TKS4) undergoes tyrosine phosphorylation by the epidermal growth factor receptor (EGFR) pathway via Src kinase. The TKS4 deficiency in humans is responsible for the manifestation of a genetic disorder known as Frank-Ter Haar syndrome (FTHS). Based on our earlier investigation, the absence of TKS4 triggers migration, invasion, and epithelial-mesenchymal transition (EMT)-like phenomena while concurrently suppressing cell proliferation in HCT116 colorectal carcinoma cells. This indicates that TKS4 may play a unique role in the progression of cancer. In this study, we demonstrated that the enhancer of zeste homolog 2 (EZH2) and the histone methyltransferase of polycomb repressive complex 2 (PRC2) are involved in the migration, invasion, and EMT-like changes in TKS4-deficient cells (KO). EZH2 is responsible for the maintenance of the trimethylated lysine 27 on histone H3 (H3K27me3). METHODS: We performed transcriptome sequencing, chromatin immunoprecipitation, protein and RNA quantitative studies, cell mobility, invasion, and proliferation studies combined with/without the EZH2 activity inhibitor 3-deazanoplanocine (DZNep). RESULTS: We detected an elevation of global H3K27me3 levels in the TKS4 KO cells, which could be reduced with treatment with DZNep, an EZH2 inhibitor. Inhibition of EZH2 activity reversed the phenotypic effects of the knockout of TKS4, reducing the migration speed and wound healing capacity of the cells as well as decreasing the invasion capacity, while the decrease in cell proliferation became stronger. In addition, inhibition of EZH2 activity also reversed most epithelial and mesenchymal markers. We investigated the wider impact of TKS4 deletion on the gene expression profile of colorectal cancer cells using transcriptome sequencing of wild-type and TKS4 knockout cells, particularly before and after treatment with DZNep. Additionally, we observed changes in the expression of several protein-coding genes and long non-coding RNAs that showed a recovery in expression levels following EZH2 inhibition. CONCLUSIONS: Our results indicate that the removal of TKS4 causes a notable disruption in the gene expression pattern, leading to the disruption of several signal transduction pathways. Inhibiting the activity of EZH2 can restore most of these transcriptomics and phenotypic effects in colorectal carcinoma cells.

A cyclin D1 intrinsically disordered domain accesses modified histone motifs to govern gene transcription.

The essential G1-cyclin, CCND1, is frequently overexpressed in cancer, contributing to tumorigenesis by driving cell-cycle progression. D-type cyclins are rate-limiting regulators of G1-S progression in mammalian cells via their ability to bind and activate CDK4 and CDK6. In addition, cyclin D1 conveys kinase-independent transcriptional functions of cyclin D1. Here we report that cyclin D1 associates with H2BS14 via an intrinsically disordered domain (IDD). The same region of cyclin D1 was necessary for the induction of aneuploidy, induction of the DNA damage response, cyclin D1-mediated recruitment into chromatin, and CIN gene transcription. In response to DNA damage H2BS14 phosphorylation occurs, resulting in co-localization with γH2AX in DNA damage foci. Cyclin D1 ChIP seq and γH2AX ChIP seq revealed ~14% overlap. As the cyclin D1 IDD functioned independently of the CDK activity to drive CIN, the IDD domain may provide a rationale new target to complement CDK-extinction strategies.

KMT2D preferentially binds mRNAs of the genes it regulates, suggesting a role in RNA processing.

Histone lysine methyltransferases (HKMTs) perform vital roles in cellular life by controlling gene expression programs through the posttranslational modification of histone tails. Since many of them are intimately involved in the development of different diseases, including several cancers, understanding the molecular mechanisms that control their target recognition and activity is vital for the treatment and prevention of such conditions. RNA binding has been shown to be an important regulatory factor in the function of several HKMTs, such as the yeast Set1 and the human Ezh2. Moreover, many HKMTs are capable of RNA binding in the absence of a canonical RNA binding domain. Here, we explored the RNA binding capacity of KMT2D, one of the major H3K4 monomethyl transferases in enhancers, using RNA immunoprecipitation followed by sequencing. We identified a broad range of coding and non-coding RNAs associated with KMT2D and confirmed their binding through RNA immunoprecipitation and quantitative PCR. We also showed that a separated RNA binding region within KMT2D is capable of binding a similar RNA pool, but differences in the binding specificity indicate the existence of other regulatory elements in the sequence of KMT2D. Analysis of the bound mRNAs revealed that KMT2D preferentially binds co-transcriptionally to the mRNAs of the genes under its control, while also interacting with super enhancer- and splicing-related non-coding RNAs. These observations, together with the nuclear colocalization of KMT2D with differentially phosphorylated forms of RNA Polymerase II suggest a so far unexplored role of KMT2D in the RNA processing of the nascent transcripts.

In Vivo and In Vitro Characterization of the RNA Binding Capacity of SETD1A (KMT2F).

For several histone lysine methyltransferases (HKMTs), RNA binding has been already shown to be a functionally relevant feature, but detailed information on the RNA interactome of these proteins is not always known. Of the six human KMT2 proteins responsible for the methylation of the H3K4 residue, two-SETD1A and SETD1B-contain RNA recognition domains (RRMs). Here we investigated the RNA binding capacity of SETD1A and identified a broad range of interacting RNAs within HEK293T cells. Our analysis revealed that similar to yeast Set1, SETD1A is also capable of binding several coding and non-coding RNAs, including RNA species related to RNA processing. We also show direct RNA binding activity of the individual RRM domain in vitro, which is in contrast with the RRM domain found in yeast Set1. Structural modeling revealed important details on the possible RNA recognition mode of SETD1A and highlighted some fundamental differences between SETD1A and Set1, explaining the differences in the RNA binding capacity of their respective RRMs.

Absence of Scaffold Protein Tks4 Disrupts Several Signaling Pathways in Colon Cancer Cells.

Tks4 is a large scaffold protein in the EGFR signal transduction pathway that is involved in several cellular processes, such as cellular motility, reactive oxygen species-dependent processes, and embryonic development. It is also implicated in a rare developmental disorder, Frank-ter Haar syndrome. Loss of Tks4 resulted in the induction of an EMT-like process, with increased motility and overexpression of EMT markers in colorectal carcinoma cells. In this work, we explored the broader effects of deletion of Tks4 on the gene expression pattern of HCT116 colorectal carcinoma cells by transcriptome sequencing of wild-type and Tks4 knockout (KO) cells. We identified several protein coding genes with altered mRNA levels in the Tks4 KO cell line, as well as a set of long non-coding RNAs, and confirmed these changes with quantitative PCR on a selected set of genes. Our results show a significant perturbation of gene expression upon the deletion of Tks4, suggesting the involvement of different signal transduction pathways over the well-known EGFR signaling.

Intrinsic protein disorder uncouples affinity from binding specificity.

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins often function by molecular recognition, in which they undergo induced folding. Based on prior generalizations, the idea prevails in the IDP field that due to the entropic penalty of induced folding, the major functional advantage associated with this binding mode is "uncoupling" specificity from binding strength. Nevertheless, both weaker binding and high specificity of IDPs/IDRs rest on limited experimental observations, making these assumptions more speculations than evidence-supported facts. The issue is also complicated by the rather vague concept of specificity that lacks an exact measure, such as the Kd for binding strength. We addressed these issues by creating and analyzing a comprehensive dataset of well-characterized ID/globular protein complexes, for which both the atomic structure of the complex and free energy (ΔG, Kd ) of interaction is known. Through this analysis, we provide evidence that the affinity distributions of IDP/globular and globular/globular complexes show different trends, whereas specificity does not connote to weaker binding strength of IDPs/IDRs. Furthermore, protein disorder extends the spectrum in the direction of very weak interactions, which may have important regulatory consequences and suggest that, in a biological sense, strict correlation of specificity and binding strength are uncoupled by structural disorder.

Functional Characterization of the N-Terminal Disordered Region of the piggyBac Transposase.

The piggyBac DNA transposon is an active element initially isolated from the cabbage looper moth, but members of this superfamily are also present in most eukaryotic evolutionary lineages. The functionally important regions of the transposase are well described. There is an RNase H-like fold containing the DDD motif responsible for the catalytic DNA cleavage and joining reactions and a C-terminal cysteine-rich domain important for interaction with the transposon DNA. However, the protein also contains a ~100 amino acid long N-terminal disordered region (NTDR) whose function is currently unknown. Here we show that deletion of the NTDR significantly impairs piggyBac transposition, although the extent of decrease is strongly cell-type specific. Moreover, replacing the NTDR with scrambled but similarly disordered sequences did not rescue transposase activity, indicating the importance of sequence conservation. Cell-based transposon excision and integration assays reveal that the excision step is more severely affected by NTDR deletion. Finally, bioinformatic analyses indicated that the NTDR is specific for the piggyBac superfamily and is also present in domesticated, transposase-derived proteins incapable of catalyzing transposition. Our results indicate an essential role of the NTDR in the "fine-tuning" of transposition and its significance in the functions of piggyBac-originated co-opted genes.

Wide-Line NMR Melting Diagrams, Their Thermodynamic Interpretation, and Secondary Structure Predictions for A30P and E46K α-Synuclein.

Parkinson's disease is thought to be caused by aggregation of the intrinsically disordered protein, α-synuclein. Two amyloidogenic variants, A30P, and E46K familial mutants were investigated by wide-line 1H NMR spectrometry as a completion of our earlier work on wild-type and A53T α-synuclein (Bokor M. et al. WT and A53T α-synuclein systems: melting diagram and its new interpretation. Int. J. Mol. Sci.2020, 21, 3997.). A monolayer of mobile water molecules hydrates A30P α-synuclein at the lowest potential barriers (temperatures), while E46K α-synuclein has here third as much mobile hydration, insufficient for functionality. According to wide-line 1H NMR results and secondary structure predictions, E46K α-synuclein is more compact than the A30P variant and they are more compact than the wild type (WT) and A53T variants. Linear hydration vs potential barrier sections of A30P α-synuclein shows one and E46K shows two slopes. The different slopes of the latter between potential barriers E a,1 and E a,2 reflect a change in water-protein interactions. The 31-32% homogeneous potential barrier distribution of the protein-water bonds refers to a non-negligible amount of secondary structures in all four α-synuclein variants. The secondary structures detected by wide-line 1H NMR are solvent-exposed α-helices, which are predicted by secondary structure models. ß-sheets are only minor components of the protein structures as three- and eight-state predicted secondary structures are dominated by α-helices and coils.

The Disordered EZH2 Loop: Atomic Level Characterization by 1HN- and 1Hα-Detected NMR Approaches, Interaction with the Long Noncoding HOTAIR RNA.

The 96-residue-long loop of EZH2 is proposed to play a role in the interaction with long non-coding RNAs (lncRNAs) and to contribute to EZH2 recruitment to the chromatin. However, molecular details of RNA recognition have not been described so far. Cellular studies have suggested that phosphorylation of the Thr345 residue localized in this loop influences RNA binding; however, no mechanistic explanation has been offered. To address these issues, a systematic NMR study was performed. As the 1HN-detected NMR approach presents many challenges under physiological conditions, our earlier developed, as well as improved, 1Hα-detected experiments were used. As a result of the successful resonance assignment, the obtained chemical shift values indicate the highly disordered nature of the EZH2 loop, with some nascent helical tendency in the Ser407-Ser412 region. Further investigations conducted on the phosphomimetic mutant EZH2T345D showed that the mutation has only a local effect, and that the loop remains disordered. On the other hand, the mutation influences the cis/trans Pro346 equilibrium. Interactions of both the wild-type and the phosphomimetic mutant with the lncRNA HOTAIR140 (1-140 nt) highlight that the Thr367-Ser375 region is affected. This segment does not resemble any of the previously reported RNA-binding motifs, therefore the identified binding region is unique. As no structural changes occur in the EZH2 loop upon RNA binding, we can consider the protein-RNA interaction as a "fuzzy" complex.

BeStSel: webserver for secondary structure and fold prediction for protein CD spectroscopy.

Circular dichroism (CD) spectroscopy is widely used to characterize the secondary structure composition of proteins. To derive accurate and detailed structural information from the CD spectra, we have developed the Beta Structure Selection (BeStSel) method (PNAS, 112, E3095), which can handle the spectral diversity of ß-structured proteins. The BeStSel webserver provides this method with useful accessories to the community with the main goal to analyze single or multiple protein CD spectra. Uniquely, BeStSel provides information on eight secondary structure components including parallel ß-structure and antiparallel ß-sheets with three different groups of twist. It overperforms any available method in accuracy and information content, moreover, it is capable of predicting the protein fold down to the topology/homology level of the CATH classification. A new module of the webserver helps to distinguish intrinsically disordered proteins by their CD spectrum. Secondary structure calculation for uploaded PDB files will help the experimental verification of protein MD and in silico modelling using CD spectroscopy. The server also calculates extinction coefficients from the primary sequence for CD users to determine the accurate protein concentrations which is a prerequisite for reliable secondary structure determination. The BeStSel server can be freely accessed at https://bestsel.elte.hu.

Disordered-Ordered Protein Binary Classification by Circular Dichroism Spectroscopy.

Intrinsically disordered proteins lack a stable tertiary structure and form dynamic conformational ensembles due to their characteristic physicochemical properties and amino acid composition. They are abundant in nature and responsible for a large variety of cellular functions. While numerous bioinformatics tools have been developed for in silico disorder prediction in the last decades, there is a need for experimental methods to verify the disordered state. CD spectroscopy is widely used for protein secondary structure analysis. It is usable in a wide concentration range under various buffer conditions. Even without providing high-resolution information, it is especially useful when NMR, X-ray, or other techniques are problematic or one simply needs a fast technique to verify the structure of proteins. Here, we propose an automatized binary disorder-order classification method by analyzing far-UV CD spectroscopy data. The method needs CD data at only three wavelength points, making high-throughput data collection possible. The mathematical analysis applies the k-nearest neighbor algorithm with cosine distance function, which is independent of the spectral amplitude and thus free of concentration determination errors. Moreover, the method can be used even for strong absorbing samples, such as the case of crowded environmental conditions, if the spectrum can be recorded down to the wavelength of 212 nm. We believe the classification method will be useful in identifying disorder and will also facilitate the growth of experimental data in IDP databases. The method is implemented on a webserver and freely available for academic users.

Identification of Intrinsically Disordered Proteins and Regions in a Non-Model Insect Species Ostrinia nubilalis (Hbn.).

Research in previous decades has shown that intrinsically disordered proteins (IDPs) and regions in proteins (IDRs) are as ubiquitous as highly ordered proteins. Despite this, research on IDPs and IDRs still has many gaps left to fill. Here, we present an approach that combines wet lab methods with bioinformatics tools to identify and analyze intrinsically disordered proteins in a non-model insect species that is cold-hardy. Due to their known resilience to the effects of extreme temperatures, these proteins likely play important roles in this insect's adaptive mechanisms to sub-zero temperatures. The approach involves IDP enrichment by sample heating and double-digestion of proteins, followed by peptide and protein identification. Next, proteins are bioinformatically analyzed for disorder content, presence of long disordered regions, amino acid composition, and processes they are involved in. Finally, IDP detection is validated with an in-house 2D PAGE. In total, 608 unique proteins were identified, with 39 being mostly disordered, 100 partially disordered, 95 nearly ordered, and 374 ordered. One-third contain at least one long disordered segment. Functional information was available for only 90 proteins with intrinsic disorders out of 312 characterized proteins. Around half of the 90 proteins are cytoskeletal elements or involved in translational processes.

Deep structural insights into RNA-binding disordered protein regions.

Recent efforts to identify RNA binding proteins in various organisms and cellular contexts have yielded a large collection of proteins that are capable of RNA binding in the absence of conventional RNA recognition domains. Many of the recently identified RNA interaction motifs fall into intrinsically disordered protein regions (IDRs). While the recognition mode and specificity of globular RNA binding elements have been thoroughly investigated and described, much less is known about the way IDRs can recognize their RNA partners. Our aim was to summarize the current state of structural knowledge on the RNA binding modes of disordered protein regions and to propose a classification system based on their sequential and structural properties. Through a detailed structural analysis of the complexes that contain disordered protein regions binding to RNA, we found two major binding modes that represent different recognition strategies and, most likely, functions. We compared these examples with DNA binding disordered proteins and found key differences stemming from the nucleic acids as well as similar binding strategies, implying a broader substrate acceptance by these proteins. Due to the very limited number of known structures, we integrated molecular dynamics simulations in our study, whose results support the proposed structural preferences of specific RNA-binding IDRs. To broaden the scope of our review, we included a brief analysis of RNA-binding small molecules and compared their structural characteristics and RNA recognition strategies to the RNA-binding IDRs. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions.

Cellular Chaperone Function of Intrinsically Disordered Dehydrin ERD14.

Disordered plant chaperones play key roles in helping plants survive in harsh conditions, and they are indispensable for seeds to remain viable. Aside from well-known and thoroughly characterized globular chaperone proteins, there are a number of intrinsically disordered proteins (IDPs) that can also serve as highly effective protecting agents in the cells. One of the largest groups of disordered chaperones is the group of dehydrins, proteins that are expressed at high levels under different abiotic stress conditions, such as drought, high temperature, or osmotic stress. Dehydrins are characterized by the presence of different conserved sequence motifs that also serve as the basis for their categorization. Despite their accepted importance, the exact role and relevance of the conserved regions have not yet been formally addressed. Here, we explored the involvement of each conserved segment in the protective function of the intrinsically disordered stress protein (IDSP) A. thaliana's Early Response to Dehydration (ERD14). We show that segments that are directly involved in partner binding, and others that are not, are equally necessary for proper function and that cellular protection emerges from the balanced interplay of different regions of ERD14.

Protein-Protein Connections-Oligomer, Amyloid and Protein Complex-By Wide Line 1H NMR.

The amount of bonds between constituting parts of a protein aggregate were determined in wild type (WT) and A53T α-synuclein (αS) oligomers, amyloids and in the complex of thymosin-ß4-cytoplasmic domain of stabilin-2 (Tß4-stabilin CTD). A53T αS aggregates have more extensive ßsheet contents reflected by constant regions at low potential barriers in difference (to monomers) melting diagrams (MDs). Energies of the intermolecular interactions and of secondary structures bonds, formed during polymerization, fall into the 5.41 kJ mol-1 ≤ Ea ≤ 5.77 kJ mol-1 range for αS aggregates. Monomers lose more mobile hydration water while forming amyloids than oligomers. Part of the strong mobile hydration water-protein bonds break off and these bonding sites of the protein form intermolecular bonds in the aggregates. The new bonds connect the constituting proteins into aggregates. Amyloid-oligomer difference MD showed an overall more homogeneous solvent accessible surface of A53T αS amyloids. From the comparison of the nominal sum of the MDs of the constituting proteins to the measured MD of the Tß4-stabilin CTD complex, the number of intermolecular bonds connecting constituent proteins into complex is 20(1) H2O/complex. The energies of these bonds are in the 5.40(3) kJ mol-1 ≤ Ea ≤ 5.70(5) kJ mol-1 range.

Secondary Structures of Proteins: A Comparison of Models and Experimental Results.

Secondary structure predictions of proteins were compared to experimental results by wide-line 1H NMR. IUPred2A was used to generate predictions of disordered protein or binding regions. Thymosin-ß4 and the stabilin-2 cytoplasmic domain were found to be mainly disordered, in agreement with the experimental results. α-Synuclein variants were predicted to be disordered, as in the experiments, but the A53T mutant showed less predicted disorder, in contrast with the wide-line 1H NMR result. A disordered binding site was found for thymosin-ß4, whereas the stabilin-2 cytoplasmic domain was indicated as such in its entire length. The last third of the α-synuclein variant's sequence was a disordered binding site. Thymosin-ß4 and the stabilin-2 cytoplasmic domain contained only coils and helices according to five secondary structure prediction methods (SPIDER3-SPOT-1D, PSRSM, MUFold-SSW, Porter 5, and RaptorX). ß-Sheets are present in α-synucleins, and they extend to more amino acid residues in the A53T mutant according to the predictions. The latter is verified by experiments. The comparison of the predictions with the experiments suggests that helical parts are buried.

A generic approach to study the kinetics of liquid-liquid phase separation under near-native conditions.

Understanding the kinetics, thermodynamics, and molecular mechanisms of liquid-liquid phase separation (LLPS) is of paramount importance in cell biology, requiring reproducible methods for studying often severely aggregation-prone proteins. Frequently applied approaches for inducing LLPS, such as dilution of the protein from an urea-containing solution or cleavage of its fused solubility tag, often lead to very different kinetic behaviors. Here we demonstrate that at carefully selected pH values proteins such as the low-complexity domain of hnRNPA2, TDP-43, and NUP98, or the stress protein ERD14, can be kept in solution and their LLPS can then be induced by a jump to native pH. This approach represents a generic method for studying the full kinetic trajectory of LLPS under near native conditions that can be easily controlled, providing a platform for the characterization of physiologically relevant phase-separation behavior of diverse proteins.