An Atomistic View on the Mechanism of Diatom Peptide-Guided Biomimetic Silica Formation.

Deciphering nature's remarkable way of encoding functions in its biominerals holds the potential to enable the rational development of nature-inspired materials with tailored properties. However, the complex processes that convert solution-state precursors into solid biomaterials remain largely unknown. In this study, an unconventional approach is presented to characterize these precursors for the diatom-derived peptides R5 and synthetic Silaffin-1A1 (synSil-1A1). These molecules can form defined supramolecular assemblies in solution, which act as templates for solid silica structures. Using a tailored structural biology toolbox, the structure-function relationships of these self-assemblies are unveiled. NMR-derived constraints are employed to enable a recently developed fractal-cluster formalism and then reveal the architecture of the peptide assemblies in atomistic detail. Finally, by monitoring the self-assembly activities during silica formation at simultaneous high temporal and residue resolution using real-time spectroscopy, the mechanism is elucidated underlying template-driven silica formation. Thus, it is demonstrated how to exercise morphology control over bioinorganic solids by manipulating the template architectures. It is found that the morphology of the templates is translated into the shape of bioinorganic particles via a mechanism that includes silica nucleation on the solution-state complexes' surfaces followed by complete surface coating and particle precipitation.

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.

On a mechanistic impact of transmembrane tetramerization in the pathological activation of RTKs.

Constitutive activation of receptor tyrosine kinases (RTKs) via different mutations has a strong impact on the development of severe human disorders, including cancer. Here we propose a putative activation scenario of RTKs, whereby transmembrane (TM) mutations can also promote higher-order oligomerization of the receptors that leads to the subsequent ligand-free activation. We illustrate this scenario using a computational modelling framework comprising sequence-based structure prediction and all-atom 1 µs molecular dynamics (MD) simulations in a lipid membrane for a previously characterised oncogenic TM mutation V536E in platelet-derived growth factor receptor alpha (PDGFRA). We show that in the course of MD simulations the mutant TM tetramer retains stable and compact configuration strengthened by tight protein-protein interactions, while the wild type TM tetramer demonstrates looser packing and a tendency to dissociate. Moreover, the mutation affects the characteristic motions of mutated TM helical segments by introducing additional non-covalent crosslinks in the middle of the TM tetramer, which operate as mechanical hinges. This leads to dynamic decoupling of the C-termini from the rigidified N-terminal parts and facilitates more pronounced possible displacement between the C-termini of the mutant TM helical regions that can provide more freedom for mutual rearrangement of the kinase domains located downstream. Our results for the V536E mutation in the context of PDGFRA TM tetramer allow for the possibility that the effect of oncogenic TM mutations can go beyond alternating the structure and dynamics of TM dimeric states and might also promote the formation of higher-order oligomers directly contributing to ligand-independent signalling effectuated by PDGFRA and other RTKs.

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.

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.

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.

A Uniquely Stable Trimeric Model of SARS-CoV-2 Spike Transmembrane Domain.

Understanding fusion mechanisms employed by SARS-CoV-2 spike protein entails realistic transmembrane domain (TMD) models, while no reliable approaches towards predicting the 3D structure of transmembrane (TM) trimers exist. Here, we propose a comprehensive computational framework to model the spike TMD only based on its primary structure. We performed amino acid sequence pattern matching and compared the molecular hydrophobicity potential (MHP) distribution on the helix surface against TM homotrimers with known 3D structures and selected an appropriate template for homology modeling. We then iteratively built a model of spike TMD, adjusting "dynamic MHP portraits" and residue variability motifs. The stability of this model, with and without palmitoyl modifications downstream of the TMD, and several alternative configurations (including a recent NMR structure), was tested in all-atom molecular dynamics simulations in a POPC bilayer mimicking the viral envelope. Our model demonstrated unique stability under the conditions applied and conforms to known basic principles of TM helix packing. The original computational framework looks promising and could potentially be employed in the construction of 3D models of TM trimers for a wide range of membrane proteins.

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.

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.

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.

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.

A novel non-canonical PIP-box mediates PARG interaction with PCNA.

Poly(ADP-ribose) glycohydrolase (PARG) regulates cellular poly(ADP-ribose) (PAR) levels by rapidly cleaving glycosidic bonds between ADP-ribose units. PARG interacts with proliferating cell nuclear antigen (PCNA) and is strongly recruited to DNA damage sites in a PAR- and PCNA-dependent fashion. Here we identified PARG acetylation site K409 that is essential for its interaction with PCNA, its localization within replication foci and its recruitment to DNA damage sites. We found K409 to be part of a non-canonical PIP-box within the PARG disordered regulatory region. The previously identified putative N-terminal PIP-box does not bind PCNA directly but contributes to PARG localization within replication foci. X-ray structure and MD simulations reveal that the PARG non-canonical PIP-box binds PCNA in a manner similar to other canonical PIP-boxes and may represent a new type of PIP-box. While the binding of previously described PIP-boxes is based on hydrophobic interactions, PARG PIP-box binds PCNA via both stabilizing hydrophobic and fine-tuning electrostatic interactions. Our data explain the mechanism of PARG-PCNA interaction through a new PARG PIP-box that exhibits non-canonical sequence properties but a canonical mode of PCNA binding.

Dependence of Binding Free Energies between RNA Nucleobases and Protein Side Chains on Local Dielectric Properties.

In order to fully understand the microscopic origins of binding specificity between nucleic acids and proteins, it is imperative to study the dependence of the binding preferences between nucleobases and protein side chains on the properties of the environment. Here, we employ molecular dynamics simulations and umbrella sampling to derive the potentials of mean force and the associated absolute binding free energies between the four standard RNA nucleobases and the side chains of aspartic acid and tryptophan in water/methanol mixtures exhibiting a wide range of dielectric constants. In addition to their opposing character when it comes to hydrophobicity, aspartate and tryptophan side chains were chosen because they exhibit the greatest change in binding free energies with nucleobases between pure water and methanol environments. We exploit a strong linear dependence of the derived ΔG values on the mole fraction of methanol to estimate the binding free energies of all possible combinations of different standard RNA nucleobases and side chains at multiple values of dielectric constants. Finally, we critically assess the recently proposed complementarity hypothesis concerning direct, coaligned binding between mRNAs and their cognate proteins in light of the present results.

Fuento: functional enrichment for bioinformatics.

SUMMARY: The currently available functional enrichment software focuses mostly on gene expression analysis, whereby server- and graphical-user-interface-based tools with specific scope dominate the field. Here we present an efficient, user-friendly, multifunctional command-line-based functional enrichment tool (fu-en-to), tailored for the bioinformatics researcher. AVAILABILITY AND IMPLEMENTATION: Source code and binaries freely available for download at github.com/DavidWeichselbaum/fuento, implemented in C ++ and supported on Linux and OS X. CONTACT: newant@gmail.com or bojan.zagrovic@univie.ac.at.

The Conformation of the Epidermal Growth Factor Receptor Transmembrane Domain Dimer Dynamically Adapts to the Local Membrane Environment.

The epidermal growth factor receptor (EGFR) family is an important class of receptor tyrosine kinases, mediating a variety of cellular responses in normal biological processes and in pathological states of multicellular organisms. Different modes of dimerization of the human EGFR transmembrane domain (TMD) in different membrane mimetics recently prompted us to propose a novel signal transduction mechanism based on protein-lipid interaction. However, the experimental evidence for it was originally obtained with slightly different TMD fragments used in the two different mimetics, compromising the validity of the comparison. To eliminate ambiguity, we determined the nuclear magnetic resonance (NMR) structure of the bicelle-incorporated dimer of the EGFR TMD fragment identical to the one previously used in micelles. The NMR results augmented by molecular dynamics simulations confirm the mutual influence of the TMD and lipid environment, as is required for the proposed lipid-mediated activation mechanism. They also reveal the possible functional relevance of a subtle interplay between the concurrent processes in the lipid and protein during signal transduction.

Structural mechanism for the recognition and ubiquitination of a single nucleosome residue by Rad6-Bre1.

Cotranscriptional ubiquitination of histone H2B is key to gene regulation. The yeast E3 ubiquitin ligase Bre1 (human RNF20/40) pairs with the E2 ubiquitin conjugating enzyme Rad6 to monoubiquitinate H2B at Lys123. How this single lysine residue on the nucleosome core particle (NCP) is targeted by the Rad6-Bre1 machinery is unknown. Using chemical cross-linking and mass spectrometry, we identified the functional interfaces of Rad6, Bre1, and NCPs in a defined in vitro system. The Bre1 RING domain cross-links exclusively with distinct regions of histone H2B and H2A, indicating a spatial alignment of Bre1 with the NCP acidic patch. By docking onto the NCP surface in this distinct orientation, Bre1 positions the Rad6 active site directly over H2B Lys123. The Spt-Ada-Gcn5 acetyltransferase (SAGA) H2B deubiquitinase module competes with Bre1 for binding to the NCP acidic patch, indicating regulatory control. Our study reveals a mechanism that ensures site-specific NCP ubiquitination and fine-tuning of opposing enzymatic activities.