N6-methyladenosine in 5' UTR does not promote translation initiation.

The most abundant N6-methyladenosine (m6A) modification on mRNAs is installed non-stoichiometrically across transcripts, with 5' untranslated regions (5' UTRs) being the least conductive. 5' UTRs are essential for translation initiation, yet the molecular mechanisms orchestrated by m6A remain poorly understood. Here, we combined structural, biochemical, and single-molecule approaches and show that at the most common position, a single m6A does not affect translation yields, the kinetics of translation initiation complex assembly, or start codon recognition both under permissive growth and following exposure to oxidative stress. Cryoelectron microscopy (cryo-EM) structures of the late preinitiation complex reveal that m6A purine ring established stacking interactions with an arginine side chain of the initiation factor eIF2α, although with only a marginal energy contribution, as estimated computationally. These findings provide molecular insights into m6A interactions with the initiation complex and suggest that the subtle stabilization is unlikely to affect the translation dynamics under homeostatic conditions or stress.

The SPOC proteins DIDO3 and PHF3 co-regulate gene expression and neuronal differentiation.

Transcription is regulated by a multitude of activators and repressors, which bind to the RNA polymerase II (Pol II) machinery and modulate its progression. Death-inducer obliterator 3 (DIDO3) and PHD finger protein 3 (PHF3) are paralogue proteins that regulate transcription elongation by docking onto phosphorylated serine-2 in the C-terminal domain (CTD) of Pol II through their SPOC domains. Here, we show that DIDO3 and PHF3 form a complex that bridges the Pol II elongation machinery with chromatin and RNA processing factors and tethers Pol II in a phase-separated microenvironment. Their SPOC domains and C-terminal intrinsically disordered regions are critical for transcription regulation. PHF3 and DIDO exert cooperative and antagonistic effects on the expression of neuronal genes and are both essential for neuronal differentiation. In the absence of PHF3, DIDO3 is upregulated as a compensatory mechanism. In addition to shared gene targets, DIDO specifically regulates genes required for lipid metabolism. Collectively, our work reveals multiple layers of gene expression regulation by the DIDO3 and PHF3 paralogues, which have specific, co-regulatory and redundant functions in transcription.

Driving forces behind phase separation of the carboxy-terminal domain of RNA polymerase II.

Eukaryotic gene regulation and pre-mRNA transcription depend on the carboxy-terminal domain (CTD) of RNA polymerase (Pol) II. Due to its highly repetitive, intrinsically disordered sequence, the CTD enables clustering and phase separation of Pol II. The molecular interactions that drive CTD phase separation and Pol II clustering are unclear. Here, we show that multivalent interactions involving tyrosine impart temperature- and concentration-dependent self-coacervation of the CTD. NMR spectroscopy, molecular ensemble calculations and all-atom molecular dynamics simulations demonstrate the presence of diverse tyrosine-engaging interactions, including tyrosine-proline contacts, in condensed states of human CTD and other low-complexity proteins. We further show that the network of multivalent interactions involving tyrosine is responsible for the co-recruitment of the human Mediator complex and CTD during phase separation. Our work advances the understanding of the driving forces of CTD phase separation and thus provides the basis to better understand CTD-mediated Pol II clustering in eukaryotic gene transcription.

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

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

The SPOC domain is a phosphoserine binding module that bridges transcription machinery with co- and post-transcriptional regulators.

The heptad repeats of the C-terminal domain (CTD) of RNA polymerase II (Pol II) are extensively modified throughout the transcription cycle. The CTD coordinates RNA synthesis and processing by recruiting transcription regulators as well as RNA capping, splicing and 3'end processing factors. The SPOC domain of PHF3 was recently identified as a CTD reader domain specifically binding to phosphorylated serine-2 residues in adjacent CTD repeats. Here, we establish the SPOC domains of the human proteins DIDO, SHARP (also known as SPEN) and RBM15 as phosphoserine binding modules that can act as CTD readers but also recognize other phosphorylated binding partners. We report the crystal structure of SHARP SPOC in complex with CTD and identify the molecular determinants for its specific binding to phosphorylated serine-5. PHF3 and DIDO SPOC domains preferentially interact with the Pol II elongation complex, while RBM15 and SHARP SPOC domains engage with writers and readers of m6A, the most abundant RNA modification. RBM15 positively regulates m6A levels and mRNA stability in a SPOC-dependent manner, while SHARP SPOC is essential for its localization to inactive X-chromosomes. Our findings suggest that the SPOC domain is a major interface between the transcription machinery and regulators of transcription and co-transcriptional processes.

Coding From Binding? Molecular Interactions at the Heart of Translation.

The mechanism and the evolution of DNA replication and transcription, the key elements of the central dogma of biology, are fundamentally well explained by the physicochemical complementarity between strands of nucleic acids. However, the determinants that have shaped the third part of the dogma-the process of biological translation and the universal genetic code-remain unclear. We review and seek parallels between different proposals that view the evolution of translation through the prism of weak, noncovalent interactions between biological macromolecules. In particular, we focus on a recent proposal that there exists a hitherto unrecognized complementarity at the heart of biology, that between messenger RNA coding regions and the proteins that they encode, especially if the two are unstructured. Reflecting the idea that the genetic code evolved from intrinsic binding propensities between nucleotides and amino acids, this proposal promises to forge a link between the distant past and the present of biological systems.

Interaction preferences between protein side chains and key epigenetic modifications 5-methylcytosine, 5-hydroxymethycytosine and N6-methyladenine.

Covalent modifications of standard DNA/RNA nucleobases affect epigenetic regulation of gene expression by modulating interactions between nucleic acids and protein readers. We derive here the absolute binding free energies and analyze the binding modalities between key modified nucleobases 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) and N6-methyladenine (m6A) and all non-prolyl/non-glycyl protein side chains using molecular dynamics simulations and umbrella sampling in both water and methanol, the latter mimicking the low dielectric environment at the dehydrated nucleic-acid/protein interfaces. We verify the derived affinities by comparing against a comprehensive set of high-resolution structures of nucleic-protein complexes involving 5mC. Our analysis identifies protein side chains that are highly tuned for detecting cytosine methylation as a function of the environment and can thus serve as microscopic readers of epigenetic marks. Conversely, we show that the relative ordering of sidechain affinities for 5hmC and m6A does not differ significantly from those for their precursor bases, cytosine and adenine, respectively, especially in the low dielectric environment. For those two modified bases, the effect is more nuanced and manifests itself primarily at the level of absolute changes in the binding free energy. Our results contribute towards establishing a quantitative foundation for understanding, predicting and modulating the interactions between modified nucleic acids and proteins at the atomistic level.

Widespread autogenous mRNA-protein interactions detected by CLIP-seq.

Autogenous interactions between mRNAs and the proteins they encode are implicated in cellular feedback-loop regulation, but their extent and mechanistic foundation are unclear. It was recently hypothesized that such interactions may be common, reflecting the role of intrinsic nucleobase-amino acid affinities in shaping the genetic code's structure. Here we analyze a comprehensive set of CLIP-seq experiments involving multiple protocols and report on widespread autogenous interactions across different organisms. Specifically, 230 of 341 (67%) studied RNA-binding proteins (RBPs) interact with their own mRNAs, with a heavy enrichment among high-confidence hits and a preference for coding sequence binding. We account for different confounding variables, including physical (overexpression and proximity during translation), methodological (difference in CLIP protocols, peak callers and cell types) and statistical (treatment of null backgrounds). In particular, we demonstrate a high statistical significance of autogenous interactions by sampling null distributions of fixed-margin interaction matrices. Furthermore, we study the dependence of autogenous binding on the presence of RNA-binding motifs and structured domains in RBPs. Finally, we show that intrinsic nucleobase-amino acid affinities favor co-aligned binding between mRNA coding regions and the proteins they encode. Our results suggest a central role for autogenous interactions in RBP regulation and support the possibility of a fundamental connection between coding and binding.

FIBCD1 is an endocytic GAG receptor associated with a novel neurodevelopmental disorder.

Whole-exome sequencing of two patients with idiopathic complex neurodevelopmental disorder (NDD) identified biallelic variants of unknown significance within FIBCD1, encoding an endocytic acetyl group-binding transmembrane receptor with no known function in the central nervous system. We found that FIBCD1 preferentially binds and endocytoses glycosaminoglycan (GAG) chondroitin sulphate-4S (CS-4S) and regulates GAG content of the brain extracellular matrix (ECM). In silico molecular simulation studies and GAG binding analyses of patient variants determined that such variants are loss-of-function by disrupting FIBCD1-CS-4S association. Gene knockdown in flies resulted in morphological disruption of the neuromuscular junction and motor-related behavioural deficits. In humans and mice, FIBCD1 is expressed in discrete brain regions, including the hippocampus. Fibcd1 KO mice exhibited normal hippocampal neuronal morphology but impaired hippocampal-dependent learning. Further, hippocampal synaptic remodelling in acute slices from Fibcd1 KO mice was deficient but restored upon enzymatically modulating the ECM. Together, we identified FIBCD1 as an endocytic receptor for GAGs in the brain ECM and a novel gene associated with an NDD, revealing a critical role in nervous system structure, function and plasticity.

Compositional complementarity between genomic RNA and coat proteins in positive-sense single-stranded RNA viruses.

During packaging in positive-sense single-stranded RNA (+ssRNA) viruses, coat proteins (CPs) interact directly with multiple regions in genomic RNA (gRNA), but the underlying physicochemical principles remain unclear. Here we analyze the high-resolution cryo-EM structure of bacteriophage MS2 and show that the gRNA/CP binding sites, including the known packaging signal, overlap significantly with regions where gRNA nucleobase-density profiles match the corresponding CP nucleobase-affinity profiles. Moreover, we show that the MS2 packaging signal corresponds to the global minimum in gRNA/CP interaction energy in the unstructured state as derived using a linearly additive model and knowledge-based nucleobase/amino-acid affinities. Motivated by this, we predict gRNA/CP interaction sites for a comprehensive set of 1082 +ssRNA viruses. We validate our predictions by comparing them with site-resolved information on gRNA/CP interactions derived in SELEX and CLIP experiments for 10 different viruses. Finally, we show that in experimentally studied systems CPs frequently interact with autologous coding regions in gRNA, in agreement with both predicted interaction energies and a recent proposal that proteins in general tend to interact with own mRNAs, if unstructured. Our results define a self-consistent framework for understanding packaging in +ssRNA viruses and implicate interactions between unstructured gRNA and CPs in the process.

POSTAR3: an updated platform for exploring post-transcriptional regulation coordinated by RNA-binding proteins.

RNA-binding proteins (RBPs) play key roles in post-transcriptional regulation. Accurate identification of RBP binding sites in multiple cell lines and tissue types from diverse species is a fundamental endeavor towards understanding the regulatory mechanisms of RBPs under both physiological and pathological conditions. Our POSTAR annotation processes make use of publicly available large-scale CLIP-seq datasets and external functional genomic annotations to generate a comprehensive map of RBP binding sites and their association with other regulatory events as well as functional variants. Here, we present POSTAR3, an updated database with improvements in data collection, annotation infrastructure, and analysis that support the annotation of post-transcriptional regulation in multiple species including: we made a comprehensive update on the CLIP-seq and Ribo-seq datasets which cover more biological conditions, technologies, and species; we added RNA secondary structure profiling for RBP binding sites; we provided miRNA-mediated degradation events validated by degradome-seq; we included RBP binding sites at circRNA junction regions; we expanded the annotation of RBP binding sites, particularly using updated genomic variants and mutations associated with diseases. POSTAR3 is freely available at http://postar.ncrnalab.org.

Order from disorder in the sarcomere: FATZ forms a fuzzy but tight complex and phase-separated condensates with α-actinin.

In sarcomeres, α-actinin cross-links actin filaments and anchors them to the Z-disk. FATZ (filamin-, α-actinin-, and telethonin-binding protein of the Z-disk) proteins interact with α-actinin and other core Z-disk proteins, contributing to myofibril assembly and maintenance. Here, we report the first structure and its cellular validation of α-actinin-2 in complex with a Z-disk partner, FATZ-1, which is best described as a conformational ensemble. We show that FATZ-1 forms a tight fuzzy complex with α-actinin-2 and propose an interaction mechanism via main molecular recognition elements and secondary binding sites. The obtained integrative model reveals a polar architecture of the complex which, in combination with FATZ-1 multivalent scaffold function, might organize interaction partners and stabilize α-actinin-2 preferential orientation in Z-disk. Last, we uncover FATZ-1 ability to phase-separate and form biomolecular condensates with α-actinin-2, raising the question whether FATZ proteins can create an interaction hub for Z-disk proteins through membraneless compartmentalization during myofibrillogenesis.

RNA polymerase II-binding aptamers in human ACRO1 satellites disrupt transcription in cis.

Transcription elongation is a highly regulated process affected by many proteins, RNAs and the underlying DNA. Here we show that the nascent RNA can interfere with transcription in human cells, extending our previous findings from bacteria and yeast. We identified a variety of Pol II-binding aptamers (RAPs), prominent in repeat elements such as ACRO1 satellites, LINE1 retrotransposons and CA simple repeats, and also in several protein-coding genes. ACRO1 repeat, when translated in silico, exhibits ~50% identity with the Pol II CTD sequence. Taken together with a recent proposal that proteins in general tend to interact with RNAs similar to their cognate mRNAs, this suggests a mechanism for RAP binding. Using a reporter construct, we show that ACRO1 potently inhibits Pol II elongation in cis. We propose a novel mode of transcriptional regulation in humans, in which the nascent RNA binds Pol II to silence its own expression.

Frameshifting preserves key physicochemical properties of proteins.

Frameshifts in protein coding sequences are widely perceived as resulting in either nonfunctional or even deleterious protein products. Indeed, frameshifts typically lead to markedly altered protein sequences and premature stop codons. By analyzing complete proteomes from all three domains of life, we demonstrate that, in contrast, several key physicochemical properties of protein sequences exhibit significant robustness against +1 and -1 frameshifts. In particular, we show that hydrophobicity profiles of many protein sequences remain largely invariant upon frameshifting. For example, over 2,900 human proteins exhibit a Pearson's correlation coefficient R between the hydrophobicity profiles of the original and the +1-frameshifted variants greater than 0.7, despite an average sequence identity between the two of only 6.5% in this group. We observe a similar effect for protein sequence profiles of affinity for certain nucleobases as well as protein sequence profiles of intrinsic disorder. Finally, analysis of significance and optimality demonstrates that frameshift stability is embedded in the structure of the universal genetic code and may have contributed to shaping it. Our results suggest that frameshifting may be a powerful evolutionary mechanism for creating new proteins with vastly different sequences, yet similar physicochemical properties to the proteins from which they originate.

Direct interplay between stereochemistry and conformational preferences in aminoacylated oligoribonucleotides.

To address the structural and dynamical consequences of amino-acid attachment at 2'- or 3'-hydroxyls of the terminal ribose in oligoribonucleotides, we have performed an extensive set of molecular dynamics simulations of model aminoacylated RNA trinucleotides. Our simulations suggest that 3'-modified trinucleotides exhibit higher solvent exposure of the aminoacylester bond and may be more susceptible to hydrolysis than their 2' counterparts. Moreover, we observe an invariant adoption of well-defined collapsed and extended conformations for both stereoisomers. We show that the average conformational preferences of aminoacylated trinucleotides are determined by their nucleotide composition and are fine-tuned by amino-acid attachment. Conversely, solvent exposure of the aminoacylester bond depends on the attachment site, the nature of attached amino acid and the strength of its interactions with the bases. Importantly, aminoacylated CCA trinucleotides display a systematically higher solvent exposure of the aminoacylester bond and a weaker dependence of such exposure on sidechain interactions than other trinucleotides. These features could facilitate hydrolytic release of the amino acid, especially for 3' attachment, and may have contributed to CCA becoming the universal acceptor triplet in tRNAs. Our results provide novel atomistic details about fundamental aspects of biological translation and furnish clues about its primordial origins.

Configurational Entropy Components and Their Contribution to Biomolecular Complex Formation.

Configurational entropy change is a central constituent of the free energy change in noncovalent interactions between biomolecules. Due to both experimental and computational limitations, however, the impact of individual contributions to configurational entropy change remains underexplored. Here, we develop a novel, fully analytical framework to dissect the configurational entropy change of binding into contributions coming from molecular internal and external degrees of freedom. Importantly, this framework accounts for all coupled and uncoupled contributions in the absence of an external field. We employ our parallel implementation of the maximum information spanning tree algorithm to provide a comprehensive numerical analysis of the importance of the individual contributions to configurational entropy change on an extensive set of molecular dynamics simulations of protein binding processes. Contrary to commonly accepted assumptions, we show that different coupling terms contribute significantly to the overall configurational entropy change. Finally, while the magnitude of individual terms may be largely unpredictable a priori, the total configurational entropy change can be well approximated by rescaling the sum of uncoupled contributions from internal degrees of freedom only, providing support for NMR-based approaches for configurational entropy change estimation.

VOLPES: an interactive web-based tool for visualizing and comparing physicochemical properties of biological sequences.

The structure, dynamics and, ultimately, biological function of proteins and nucleic acids are determined by the physicochemical properties of their primary sequences. Such properties are frequently captured via one-dimensional profile plots depicting a given physicochemical variable as a function of sequence position. Hydrophobicity, charge or structural disorder in proteins or nucleobase-density in nucleic acids are routinely visualized in this manner to analyze sequences at a glance. Such visualizations, however, are typically created case-by-case in a purely static manner, employ fixed visualization parameters only and do not enable a quantitative comparison between different sequences. Here, we present VOLPES (volpes.univie.ac.at), a user-friendly web server and the corresponding JavaScript library that enable a fully interactive, multifunctional visualization, analysis and comparison of the physicochemical properties of protein and nucleic-acid sequences, allowing unprecedented insight into biological sequence data and creating a starting point for further in-depth exploration.

Atomistic mechanism of the constitutive activation of PDGFRA via its transmembrane domain.

Single-point mutations in the transmembrane (TM) region of receptor tyrosine kinases (RTKs) can lead to abnormal ligand-independent activation. We use a combination of computational modeling, NMR spectroscopy and cell experiments to analyze in detail the mechanism of how TM domains contribute to the activation of wild-type (WT) PDGFRA and its oncogenic V536E mutant. Using a computational framework, we scan all positions in PDGFRA TM helix for identification of potential functional mutations for the WT and the mutant and reveal the relationship between the receptor activity and TM dimerization via different interfaces. This strategy also allows us design a novel activating mutation in the WT (I537D) and a compensatory mutation in the V536E background eliminating its constitutive activity (S541G). We show both computationally and experimentally that single-point mutations in the TM region reshape the TM dimer ensemble and delineate the structural and dynamic determinants of spontaneous activation of PDGFRA via its TM domain. Our atomistic picture of the coupling between TM dimerization and PDGFRA activation corroborates the data obtained for other RTKs and provides a foundation for developing novel modulators of the pathological activity of PDGFRA.

RNA-protein interactions in an unstructured context.

Despite their importance, our understanding of noncovalent RNA-protein interactions is incomplete. This especially concerns the binding between RNA and unstructured protein regions, a widespread class of such interactions. Here, we review the recent experimental and computational work on RNA-protein interactions in an unstructured context with a particular focus on how such interactions may be shaped by the intrinsic interaction affinities between individual nucleobases and protein side chains. Specifically, we articulate the claim that the universal genetic code reflects the binding specificity between nucleobases and protein side chains and that, in turn, the code may be seen as the Rosetta stone for understanding RNA-protein interactions in general.

Self-Consistent Framework Connecting Experimental Proxies of Protein Dynamics with Configurational Entropy.

The recently developed NMR techniques enable estimation of protein configurational entropy change from the change in the average methyl order parameters. This experimental observable, however, does not directly measure the contribution of intramolecular couplings, protein main-chain motions, or angular dynamics. Here, we carry out a self-consistent computational analysis of the impact of these missing contributions on an extensive set of molecular dynamics simulations of different proteins undergoing binding. Specifically, we compare the configurational entropy change in protein complex formation as obtained by the maximum information spanning tree approximation (MIST), which treats the above entropy contributions directly, and the change in the average NMR methyl and NH order parameters. Our parallel implementation of MIST allows us to treat hard angular degrees of freedom as well as couplings up to full pairwise order explicitly, while still involving a high degree of sampling and tackling molecules of biologically relevant sizes. First, we demonstrate a remarkably strong linear relationship between the total configurational entropy change and the average change in both methyl and backbone-NH order parameters. Second, in contrast to canonical assumptions, we show that the main-chain and angular terms contribute significantly to the overall configurational entropy change and also scale linearly with it. Consequently, linear models starting from the average methyl order parameters are able to capture the contribution of main-chain and angular terms well. After applying the quantum-mechanical harmonic oscillator entropy formalism, we establish a similarly strong linear relationship for X-ray crystallographic B-factors. Finally, we demonstrate that the observed linear relationships remain robust against drastic undersampling and argue that they reflect an intrinsic property of compact proteins. Despite their remarkable strength, however, the above linear relationships yield estimates of configurational entropy change whose accuracy appears to be sufficient for qualitative applications only.