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-1A (synSil-1A).

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.

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.

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 Lge1 protein is required for histone H2B ubiquitination and its N-terminal intrinsically disordered fragment (Lge1) undergoes robust phase separation.

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).

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.

Interaction preferences between protein side chains and key epigenetic modifications 5-methylcytosine, 5-hydroxymethycytosine and N-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 N-methyladenine (mA) 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Inosine Nucleobase Acts as Guanine in Interactions with Protein Side Chains.

A central intermediate in purine catabolism, the inosine nucleobase hypoxanthine is also one of the most abundant modified nucleobases in RNA and plays key roles in the regulation of gene expression and determination of cell fate. It is known that hypoxanthine acts as guanine when interacting with other nucleobases and base pairs most favorably with cytosine. However, its preferences when it comes to interactions with amino acids remain unknown.

PARENT: A Parallel Software Suite for the Calculation of Configurational Entropy in Biomolecular Systems.

Accurate estimation of configurational entropy from the in silico-generated biomolecular ensembles, e.g., from molecular dynamics (MD) trajectories, is dependent strongly on exhaustive sampling for physical reasons.

Adaptable Lipid Matrix Promotes Protein-Protein Association in Membranes.

The cell membrane is "stuffed" with proteins, whose transmembrane (TM) helical domains spontaneously associate to form functionally active complexes. For a number of membrane receptors, a modulation of TM domains' oligomerization has been shown to contribute to the development of severe pathological states, thus calling for detailed studies of the atomistic aspects of the process. Despite considerable progress achieved so far, several crucial questions still remain: How do the helices recognize each other in the membrane? What is the driving force of their association? Here, we assess the dimerization free energy of TM helices along with a careful consideration of the interplay between the structure and dynamics of protein and lipids using atomistic molecular dynamics simulations in the hydrated lipid bilayer for three different model systems - TM fragments of glycophorin A, polyalanine and polyleucine peptides.

The evolutionarily conserved transcription factor PRDM12 controls sensory neuron development and pain perception.

PR homology domain-containing member 12 (PRDM12) belongs to a family of conserved transcription factors implicated in cell fate decisions. Here we show that PRDM12 is a key regulator of sensory neuronal specification in Xenopus. Modeling of human PRDM12 mutations that cause hereditary sensory and autonomic neuropathy (HSAN) revealed remarkable conservation of the mutated residues in evolution.

Modeling transmembrane domain dimers/trimers of plexin receptors: implications for mechanisms of signal transmission across the membrane.

Single-pass transmembrane (TM) receptors transmit signals across lipid bilayers by helix association or by configurational changes within preformed dimers. The structure determination for such TM regions is challenging and has mostly been accomplished by NMR spectroscopy. Recently, the computational prediction of TM dimer structures is becoming recognized for providing models, including alternate conformational states, which are important for receptor regulation.