The structure and regulation of human muscle α-actinin.

The spectrin superfamily of proteins plays key roles in assembling the actin cytoskeleton in various cell types, crosslinks actin filaments, and acts as scaffolds for the assembly of large protein complexes involved in structural integrity and mechanosensation, as well as cell signaling. α-actinins in particular are the major actin crosslinkers in muscle Z-disks, focal adhesions, and actin stress fibers. We report a complete high-resolution structure of the 200 kDa α-actinin-2 dimer from striated muscle and explore its functional implications on the biochemical and cellular level. The structure provides insight into the phosphoinositide-based mechanism controlling its interaction with sarcomeric proteins such as titin, lays a foundation for studying the impact of pathogenic mutations at molecular resolution, and is likely to be broadly relevant for the regulation of spectrin-like 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.

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.

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.

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.

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.

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.

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. Here we present for the first time the absolute binding free energies and the associated interaction modes between hypoxanthine and all standard, non-glycyl/non-prolyl amino acid side chain analogs as derived from molecular dynamics simulations and umbrella sampling in high- and low-dielectric environments. We illustrate the biological relevance of the derived affinities by providing a quantitative explanation for the specificity of hypoxanthine-guanine phosphoribosyltransferase, a key enzyme in the purine salvage pathway. Our results demonstrate that in its affinities for protein side chains, hypoxanthine closely matches guanine, much more so than its precursor adenine.

PREDDIMER: a web server for prediction of transmembrane helical dimers.

SUMMARY: Here we present PREDDIMER, a web tool for prediction of dimer structure of transmembrane (TM) helices. PREDDIMER allows (i) reconstruction of a number of dimer structures for given sequence(s) of TM protein fragments, (ii) ranking and filtering of predicted structures according to respective values of a scoring function, (iii) visualization of predicted 3D dimer structures and (iv) visualization of surface hydrophobicity of TM helices and their contacting (interface) regions represented as 2D maps. RESULTS: We implemented online the original PREDDIMER algorithm and benchmarked the server on 11 TM sequences, whose 3D dimer conformations were obtained previously by nuclear magnetic resonance spectroscopy. In the most of tested cases backbone root-mean-square deviations of closest predicted conformations from the experimental reference are below 3 Å. A randomization test displays good anticorrelation (-0.82) between values of the scoring function and statistical significance of the prediction 'by chance'. Going beyond a single dimer conformation, our web tool predicts an ensemble of possible conformations, which may be useful for explanation of a functioning of bitopic membrane proteins, e.g. receptor tyrosine kinases. AVAILABILITY AND IMPLEMENTATION: PREDDIMER can be accessed for free on the web at http://model.nmr.ru/preddimer/ CONTACT: newant@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evidence of direct complementary interactions between messenger RNAs and their cognate proteins.

Recently, the ability to interact with messenger RNA (mRNA) has been reported for a number of known RNA-binding proteins, but surprisingly also for different proteins without recognizable RNA binding domains including several transcription factors and metabolic enzymes. Moreover, direct binding to cognate mRNAs has been detected for multiple proteins, thus creating a strong impetus to search for functional significance and basic physico-chemical principles behind such interactions. Here, we derive interaction preferences between amino acids and RNA bases by analyzing binding interfaces in the known 3D structures of protein-RNA complexes. By applying this tool to human proteome, we reveal statistically significant matching between the composition of mRNA sequences and base-binding preferences of protein sequences they code for. For example, purine density profiles of mRNA sequences mirror guanine affinity profiles of cognate protein sequences with quantitative accuracy (median Pearson correlation coefficient R = -0.80 across the entire human proteome). Notably, statistically significant anti-matching is seen only in the case of adenine. Our results provide strong evidence for the stereo-chemical foundation of the genetic code and suggest that mRNAs and cognate proteins may in general be directly complementary to each other and associate, especially if unstructured.

A Lys-Trp cation-π interaction mediates the dimerization and function of the chloride intracellular channel protein 1 transmembrane domain.

Chloride intracellular channel protein 1 (CLIC1) is a dual-state protein that can exist either as a soluble monomer or in an integral membrane form. The oligomerization of the transmembrane domain (TMD) remains speculative despite it being implicated in pore formation. The extent to which electrostatic and van der Waals interactions drive folding and association of the dimorphic TMD is unknown and is complicated by the requirement of interactions favorable in both aqueous and membrane environments. Here we report a putative Lys37-Trp35 cation-π interaction and show that it stabilizes the dimeric form of the CLIC1 TMD in membranes. A synthetic 30-mer peptide comprising a K37M TMD mutant was examined in 2,2,2-trifluoroethanol, sodium dodecyl sulfate micelles, and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine liposomes using far-ultraviolet (UV) circular dichroism, fluorescence, and UV absorbance spectroscopy. Our data suggest that Lys37 is not implicated in the folding, stability, or membrane insertion of the TMD peptide. However, removal of this residue impairs the formation of dimers and higher-order oligomers. This is accompanied by a 30-fold loss of chloride influx activity, suggesting that dimerization modulates the rate of chloride conductance. We propose that, within membranes, individual TMD helices associate via a Lys37-mediated cation-π interaction to form active dimers. The latter findings are also supported by results of modeling a putative TMD dimer conformation in which Lys37 and Trp35 form cation-π pairs at the dimer interface. Dimeric helix bundles may then associate to form fully active ion channels. Thus, within a membrane-like environment, aromatic interactions involving a polar lysine side chain provide a thermodynamic driving force for helix-helix association.

Modular organization of α-toxins from scorpion venom mirrors domain structure of their targets, sodium channels.

To gain success in the evolutionary "arms race," venomous animals such as scorpions produce diverse neurotoxins selected to hit targets in the nervous system of prey. Scorpion α-toxins affect insect and/or mammalian voltage-gated sodium channels (Na(v)s) and thereby modify the excitability of muscle and nerve cells. Although more than 100 α-toxins are known and a number of them have been studied into detail, the molecular mechanism of their interaction with Na(v)s is still poorly understood. Here, we employ extensive molecular dynamics simulations and spatial mapping of hydrophobic/hydrophilic properties distributed over the molecular surface of α-toxins. It is revealed that despite the small size and relatively rigid structure, these toxins possess modular organization from structural, functional, and evolutionary perspectives. The more conserved and rigid "core module" is supplemented with the "specificity module" (SM) that is comparatively flexible and variable and determines the taxon (mammal versus insect) specificity of α-toxin activity. We further show that SMs in mammal toxins are more flexible and hydrophilic than in insect toxins. Concomitant sequence-based analysis of the extracellular loops of Na(v)s suggests that α-toxins recognize the channels using both modules. We propose that the core module binds to the voltage-sensing domain IV, whereas the more versatile SM interacts with the pore domain in repeat I of Na(v)s. These findings corroborate and expand the hypothesis on different functional epitopes of toxins that has been reported previously. In effect, we propose that the modular structure in toxins evolved to match the domain architecture of Na(v)s.

Sequence signatures of direct complementarity between mRNAs and cognate proteins on multiple levels.

A potential connection between physico-chemical properties of mRNAs and cognate proteins, with implications concerning both the origin of the genetic code and mRNA-protein interactions, is unexplored. We compare pyrimidine content of naturally occurring mRNA coding sequences with the propensity of cognate protein sequences to interact with pyrimidines. The latter is captured by polar requirement, a measure of solubility of amino acids in aqueous solutions of pyridines, heterocycles closely related to pyrimidines. We find that the higher the pyrimidine content of an mRNA, the stronger the average propensity of its cognate protein's amino acids to interact with pyridines. Moreover, window-averaged pyrimidine profiles of individual mRNAs strongly mirror polar-requirement profiles of cognate protein sequences. For example, 4953 human proteins exhibit a correlation between the two with |R| > 0.8. In other words, pyrimidine-rich mRNA regions quantitatively correspond to regions in cognate proteins containing residues soluble in pyrimidine mimetics and vice versa. Finally, by studying randomized genetic code variants we show that the universal genetic code is highly optimized to preserve these correlations. Overall, our findings redefine the stereo-chemical hypothesis concerning code's origin and provide evidence of direct complementary interactions between mRNAs and cognate proteins before development of ribosomal decoding, but also presently, especially if both are unstructured.

Hydrophobic matching controls the tilt and stability of the dimeric platelet-derived growth factor receptor (PDGFR) ß transmembrane segment.

The platelet-derived growth factor receptor ß is a member of the cell surface receptor tyrosine kinase family and dimerizes upon activation. We determined the structure of the transmembrane segment in dodecylphosphocholine micelles by liquid-state NMR and found that it forms a stable left-handed helical dimer. Solid-state NMR and oriented circular dichroism were used to measure the tilt angle of the helical segments in macroscopically aligned model membranes with different acyl chain lengths. Both methods showed that decreasing bilayer thickness (DEPC-POPC-DMPC) led to an increase in the helix tilt angle from 10° to 30° with respect to the bilayer normal. At the same time, reconstitution of the comparatively long hydrophobic segment became less effective, eventually resulting in complete protein aggregation in the short-chain lipid DLPC. Unrestrained molecular dynamics simulations of the dimer were carried out in explicit lipid bilayers (DEPC, POPC, DMPC, sphingomyelin), confirming the observed dependence of the helix tilt angle on bilayer thickness. Notably, molecular dynamics revealed that the left-handed dimer gets tilted en bloc, whereas conformational transitions to alternative (e.g. right-handed dimeric) states were not supported. The experimental data along with the simulation results demonstrate a pronounced interplay between the platelet-directed growth factor receptor ß transmembrane segment and the bilayer thickness. The effect of hydrophobic mismatch might play a key role in the redistribution and activation of the receptor within different lipid microdomains of the plasma membrane in vivo.

Role of dimerization efficiency of transmembrane domains in activation of fibroblast growth factor receptor 3.

Mutations in transmembrane (TM) domains of receptor tyrosine kinases are shown to cause a number of inherited diseases and cancer development. Here, we use a combined molecular modeling approach to understand molecular mechanism of effect of G380R and A391E mutations on dimerization of TM domains of human fibroblast growth factor receptor 3 (FGFR3). According to results of Monte Carlo conformational search in the implicit membrane and further molecular dynamics simulations, TM dimer of this receptor is able to form a number of various conformations, which differ significantly by the free energy of association in a full-atom model bilayer. The aforementioned mutations affect dimerization efficiency of TM segments and lead to repopulation of conformational ensemble for the dimer. Particularly, both mutations do not change the dimerization free energy of the predominant (putative "non-active") symmetric conformation of TM dimer, while affect dimerization efficiency of its asymmetric ("intermediate") and alternative symmetric (putative "active") models. Results of our simulations provide novel atomistic prospective of the role of G380 and A391E mutations in dimerization of TM domains of FGFR3 and their consecutive contributions to the activation pathway of the receptor.

Proteome-wide analysis reveals clues of complementary interactions between mRNAs and their cognate proteins as the physicochemical foundation of the genetic code.

Despite more than 50 years of effort, the origin of the genetic code remains enigmatic. Among different theories, the stereochemical hypothesis suggests that the code evolved as a consequence of direct interactions between amino acids and appropriate bases. If indeed true, such physicochemical foundation of the mRNA/protein relationship could also potentially lead to novel principles of protein-mRNA interactions in general. Inspired by this promise, we have recently explored the connection between the physicochemical properties of mRNAs and their cognate proteins at the proteome level. Using experimentally and computationally derived measures of solubility of amino acids in aqueous solutions of pyrimidine analogs together with knowledge-based interaction preferences of amino acids for different nucleobases, we have revealed a statistically significant matching between the composition of mRNA coding sequences and the base-binding preferences of their cognate protein sequences. Our findings provide strong support for the stereochemical hypothesis of genetic code's origin and suggest the possibility of direct complementary interactions between mRNAs and cognate proteins even in present-day cells.

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.

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.

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.