Intrinsically disordered electronegative clusters improve stability and binding specificity of RNA-binding proteins

RNA-binding proteins play crucial roles in various cellular functions and contain abundant disordered protein regions. The disordered regions in RNA-binding proteins are rich in repetitive sequences, such as poly-K/R, poly-N/Q, poly-A, and poly-G residues. Our bioinformatic analysis identified a largely neglected repetitive sequence family we define as electronegative clusters (ENCs) that contain acidic residues and/or phosphorylation sites. The abundance and length of ENCs exceed other known repetitive sequences. Despite their abundance, the functions of ENCs in RNA-binding proteins are still elusive. To investigate the impacts of ENCs on protein stability, RNA-binding affinity, and specificity, we selected one RNA-binding protein, the ribosomal biogenesis factor 15 (Nop15), as a model. We found that the Nop15 ENC increases protein stability and inhibits nonspecific RNA binding, but minimally interferes with specific RNA binding. To investigate the effect of ENCs on sequence specificity of RNA binding, we grafted an ENC to another RNA-binding protein, Ser/Arg-rich splicing factor 3. Using RNA Bind-n-Seq, we found that the engineered ENC inhibits disparate RNA motifs differently, instead of weakening all RNA motifs to the same extent. The motif site directly involved in electrostatic interaction is more susceptible to the ENC inhibition. These results suggest that one of functions of ENCs is to regulate RNA binding via electrostatic interaction. This is consistent with our finding that ENCs are also overrepresented in DNA-binding proteins, whereas underrepresented in halophiles, in which nonspecific nucleic acid binding is inhibited by high concentrations of salts.     

Intrinsically disordered inhibitor of glutamine synthetase is a functional protein with random‐coil‐like pKa values

The sequential action of glutamine synthetase (GS) and glutamate synthase (GOGAT) in cyanobacteria allows the incorporation of ammonium into carbon skeletons. In the cyanobacterium Synechocystis sp. PCC 6803, the activity of GS is modulated by the interaction with proteins, which include a 65‐residue‐long intrinsically disordered protein (IDP), the inactivating factor IF7. This interaction is regulated by the presence of charged residues in both IF7 and GS. To understand how charged amino acids can affect the binding of an IDP with its target and to provide clues on electrostatic interactions in disordered states of proteins, we measured the pK_a values of all IF7 acidic groups (Glu32, Glu36, Glu38, Asp40, Asp58, and Ser65, the backbone C‐terminus) at 100 mM NaCl concentration, by using NMR spectroscopy. We also obtained solution structures of IF7 through molecular dynamics simulation, validated them on the basis of previous experiments, and used them to obtain theoretical estimates of the pK_a values. Titration values for the two Asp and three Glu residues of IF7 were similar to those reported for random‐coil models, suggesting the lack of electrostatic interactions around these residues. Furthermore, our results suggest the presence of helical structure at the N‐terminus of the protein and of conformational changes at acidic pH values. The overall experimental and in silico findings suggest that local interactions and conformational equilibria do not play a role in determining the electrostatic features of the acidic residues of IF7.     

An N-terminal, 830 residues intrinsically disordered region of the cytoskeleton-regulatory protein supervillin contains Myosin II- and F-actin-binding sites

Supervillin, the largest member of the villin/gelsolin family, is a cytoskeleton regulating, peripheral membrane protein. Supervillin increases cell motility and promotes invasive activity in tumors. Major cytoskeletal interactors, including filamentous actin and myosin II, bind within the unique supervillin amino terminus, amino acids 1–830. The structural features of this key region of the supervillin polypeptide are unknown. Here, we utilize circular dichroism and bioinformatics sequence analysis to demonstrate that the N-terminal part of supervillin forms an extended intrinsically disordered region (IDR). Our combined data indicate that the N-terminus of human and bovine supervillin sequences (positions 1–830) represents an IDR, which is the largest IDR known to date in the villin/gelsolin family. Moreover, this result suggests a potentially novel mechanism of regulation of myosin II and F-actin via the intrinsically disordered N-terminal region of hub protein supervillin.     

Backbone-driven collapse in unfolded protein chains

Collapse of unfolded protein chains is an early event in folding. It affects structural properties of intrinsically disordered proteins, which take a considerable fraction of the human proteome. Collapse is generally believed to be driven by hydrophobic forces imposed by the presence of nonpolar amino acid side chains. Contributions from backbone hydrogen bonds to protein folding and stability, however, are controversial. To date, the experimental dissection of side-chain and backbone contributions has not yet been achieved because both types of interactions are integral parts of protein structure. Here, we realized this goal by applying mutagenesis and chemical modification on a set of disordered peptides and proteins. We measured the protein dimensions and kinetics of intra-chain diffusion of modified polypeptides at the level of individual molecules using fluorescence correlation spectroscopy, thereby avoiding artifacts commonly caused by aggregation of unfolded protein material in bulk. We found no contributions from side chains to collapse but, instead, identified backbone interactions as a source sufficient to form globules of native-like dimensions. The presence of backbone hydrogen bonds decreased polypeptide water solubility dramatically and accelerated the nanosecond kinetics of loop closure, in agreement with recent predictions from computer simulation. The presence of side chains, instead, slowed loop closure and modulated the dimensions of intrinsically disordered domains. It appeared that the transient formation of backbone interactions facilitates the diffusive search for productive conformations at the early stage of folding and within intrinsically disordered proteins.     

The pairwise energy content estimated from amino acid composition discriminates between folded and intrinsically unstructured proteins

The structural stability of a protein requires a large number of interresidue interactions. The energetic contribution of these can be approximated by low-resolution force fields extracted from known structures, based on observed amino acid pairing frequencies. The summation of such energies, however, cannot be carried out for proteins whose structure is not known or for intrinsically unstructured proteins. To overcome these limitations, we present a novel method for estimating the total pairwise interaction energy, based on a quadratic form in the amino acid composition of the protein. This approach is validated by the good correlation of the estimated and actual energies of proteins of known structure and by a clear separation of folded and disordered proteins in the energy space it defines. As the novel algorithm has not been trained on unstructured proteins, it substantiates the concept of protein disorder, i.e. that the inability to form a well-defined 3D structure is an intrinsic property of many proteins and protein domains. This property is encoded in their sequence, because their biased amino acid composition does not allow sufficient stabilizing interactions to form. By limiting the calculation to a predefined sequential neighborhood, the algorithm was turned into a position-specific scoring scheme that characterizes the tendency of a given amino acid to fall into an ordered or disordered region. This application we term IUPred and compare its performance with three generally accepted predictors, PONDR VL3H, DISOPRED2 and GlobPlot on a database of disordered proteins.     

Extended Structures in RNA Folding Intermediates Are Due to Nonnative Interactions Rather than Electrostatic Repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg²⁺ concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.     

Ensemble models of proteins and protein domains based on distance distribution restraints

Conformational ensembles of intrinsically disordered peptide chains are not fully determined by experimental observations. Uncertainty due to lack of experimental restraints and due to intrinsic disorder can be distinguished if distance distributions restraints are available. Such restraints can be obtained from pulsed dipolar electron paramagnetic resonance (EPR) spectroscopy applied to pairs of spin labels. Here, we introduce a Monte Carlo approach for generating conformational ensembles that are consistent with a set of distance distribution restraints, backbone dihedral angle statistics in known protein structures, and optionally, secondary structure propensities or membrane immersion depths. The approach is tested with simulated restraints for a terminal and an internal loop and for a protein with 69 residues by using sets of sparse restraints for underlying well‐defined conformations and for published ensembles of a premolten globule‐like and a coil‐like intrinsically disordered protein. Proteins 2016; 84:544–560. © 2016 Wiley Periodicals, Inc.     

Does Intrinsic Disorder in Proteins Favor Their Interaction with Lipids?

Intrinsically disordered proteins (IDPs) are implicated in a range of human diseases, some of which are associated with the ability to bind to lipids. Although the presence of solvent‐exposed hydrophobic regions in IDPs should favor their interactions with low‐molecular‐weight hydrophobic/amphiphilic compounds, this hypothesis has not been systematically explored as of yet. In this study, the analysis of the DisProt database with regard to the presence of lipid‐binding IDPs (LBIDPs) reveals that they comprise, at least, 15% of DisProt entries. LBIDPs are classified into four groups by ligand type, functional categories, domain structure, and conformational state. 57% of LBIDPs are classified as ordered according to the CH‐CDF analysis, and 70% of LBIDPs possess lengths of disordered regions below 50%. To investigate the lipid‐binding properties of IDPs for which lipid binding is not reported, three proteins from different conformational groups are rationally selected. They all are shown to bind linoleic (LA) and oleic (OA) acids with capacities ranging from 9 to 34 LA/OA molecules per protein molecule. The association with LA/OA causes the formation of high‐molecular‐weight lipid–protein complexes. These findings suggest that lipid binding is common among IDPs, which can favor their involvement in lipid metabolism.     

Structured proteins and proteins with intrinsic disorder

Until recently, the point of view that the unique tertiary structure is necessary for protein function has prevailed. However, recent data have demonstrated that many cell proteins do not possess such structure in isolation, although displaying a distinct function under physiological conditions. These proteins were named the naturally, or intrinsically, disordered proteins. The fraction of intrinsically disordered regions in such proteins may vary from several amino acid residues to a completely unordered sequence of several tens or even several hundreds of residues. The main distinction of these proteins from structured (globular) proteins is that they have no unique tertiary structure in isolation and acquire it only upon interaction with their partners. The conformation of these proteins in a complex is determined not only by their own amino acid sequence (as is typical of structured, or globular, proteins) but also by the interacting partner. This review discusses the structure-function relationships in structured and intrinsically disordered proteins. The intricateness of this problem and the possible ways to solve it are illustrated by the example of the EF1A elongation factor family.     

固有无序蛋白质及逆境胁迫下的分子功能

固有无序蛋白质是一类在生理条件下缺乏稳定三维结构而具有正常功能,参与信号转导、转录调控、胁迫应答等多种生物学过程的蛋白质.植物中许多逆境响应蛋白是固有无序蛋白质,通过其结构无序或部分无序区域在蛋白质 蛋白质、蛋白质 膜脂、蛋白质 核酸的互作中发挥重要作用.本文主要对固有无序蛋白质的类别、氨基酸组成和结构特点以及在逆境胁迫下其稳定细胞膜、保护核酸和蛋白质、调控基因表达等分子功能进行综述,以拓展对逆境胁迫下蛋白质作用分子机制的认识.     

AlphaFold and Implications for Intrinsically Disordered Proteins

Accurate predictions of the three-dimensional structures of proteins from their amino acid sequences have come of age. AlphaFold, a deep learning-based approach to protein structure prediction, shows remarkable success in independent assessments of prediction accuracy. A significant epoch in structural bioinformatics was the structural annotation of over 98% of protein sequences in the human proteome. Interestingly, many predictions feature regions of very low confidence, and these regions largely overlap with intrinsically disordered regions (IDRs). That over 30% of regions within the proteome are disordered is congruent with estimates that have been made over the past two decades, as intense efforts have been undertaken to generalize the structure–function paradigm to include the importance of conformational heterogeneity and dynamics. With structural annotations from AlphaFold in hand, there is the temptation to draw inferences regarding the “structures” of IDRs and their interactomes. Here, we offer a cautionary note regarding the misinterpretations that might ensue and highlight efforts that provide concrete understanding of sequence-ensemble-function relationships of IDRs. This perspective is intended to emphasize the importance of IDRs in sequence-function relationships (SERs) and to highlight how one might go about extracting quantitative SERs to make sense of how IDRs function.     

Cancer-Associated Mutations Perturb the Disordered Ensemble and Interactions of the Intrinsically Disordered p53 Transactivation Domain

Intrinsically disordered proteins (IDPs) are key components of regulatory networks that control crucial aspects of cell decision making. The intrinsically disordered transactivation domain (TAD) of tumor suppressor p53 mediates its interactions with multiple regulatory pathways to control the p53 homeostasis during the cellular response to genotoxic stress. Many cancer-associated mutations have been discovered in p53-TAD, but their structural and functional consequences are poorly understood. Here, by combining atomistic simulations, NMR spectroscopy, and binding assays, we demonstrate that cancer-associated mutations can significantly perturb the balance of p53 interactions with key activation and degradation regulators. Importantly, the four mutations studied in this work do not all directly disrupt the known interaction interfaces. Instead, at least three of these mutations likely modulate the disordered state of p53-TAD to perturb its interactions with regulators. Specifically, NMR and simulation analysis together suggest that these mutations can modulate the level of conformational expansion as well as rigidity of the disordered state. Our work suggests that the disordered conformational ensemble of p53-TAD can serve as a central conduit in regulating the response to various cellular stimuli at the protein–protein interaction level. Understanding how the disordered state of IDPs may be modulated by regulatory signals and/or disease associated perturbations will be essential in the studies on the role of IDPs in biology and diseases.     

Structural transitions and thermodynamics of a glycine-dependent riboswitch from Vibrio cholerae

Riboswitches are complex folded RNA domains found in noncoding regions of mRNA that regulate gene expression upon small molecule binding. Recently, Breaker and coworkers reported a tandem aptamer riboswitch (VCI-II) that binds glycine cooperatively. Here, we use hydroxyl radical footprinting and small-angle X-ray scattering (SAXS) to study the conformations of this tandem aptamer as a function of Mg(2+) and glycine concentration. We fit a simple three-state thermodynamic model that describes the energetic coupling between magnesium-induced folding and glycine binding. Furthermore, we characterize the structural conformations of each of the three states: In low salt with no magnesium present, the VCI-II construct has an extended overall conformation, presumably representing unfolded structures. Addition of millimolar concentrations of Mg(2+) in the absence of glycine leads to a significant compaction and partial folding as judged by hydroxyl radical protections. In the presence of millimolar Mg(2+) concentrations, the tandem aptamer binds glycine cooperatively. The glycine binding transition involves a further compaction, additional tertiary packing interactions and further uptake of magnesium ions relative to the state in high Mg(2+) but no glycine. Employing density reconstruction algorithms, we obtain low resolution 3-D structures for all three states from the SAXS measurements. These data provide a first glimpse into the structural conformations of the VCI-II aptamer, establish rigorous constraints for further modeling, and provide a framework for future mechanistic studies.     

The Ligand-Free State of the TPP Riboswitch: A Partially Folded RNA Structure

Riboswitches are elements of mRNA that regulate gene expression by undergoing structural changes upon binding of small ligands. Although the structures of several riboswitches have been solved with their ligands bound, the ligand-free states of only a few riboswitches have been characterized. The ligand-free state is as important for the functionality of the riboswitch as the ligand-bound form, but the ligand-free state is often a partially folded structure of the RNA, with conformational heterogeneity that makes it particularly challenging to study. Here, we present models of the ligand-free state of a thiamine pyrophosphate riboswitch that are derived from a combination of complementary experimental and computational modeling approaches. We obtain a global picture of the molecule using small-angle X-ray scattering data and use an RNA structure modeling software, MC-Sym, to fit local structural details to these data on an atomic scale. We have used two different approaches to obtaining these models. Our first approach develops a model of the RNA from the structures of its constituent junction fragments in isolation. The second approach treats the RNA as a single entity, without bias from the structure of its individual constituents. We find that both approaches give similar models for the ligand-free form, but the ligand-bound models differ for the two approaches, and only the models from the second approach agree with the ligand-bound structure known previously from X-ray crystallography. Our models provide a picture of the conformational changes that may occur in the riboswitch upon binding of its ligand. Our results also demonstrate the power of combining experimental small-angle X-ray scattering data with theoretical structure prediction tools in the determination of RNA structures beyond riboswitches.     

The carboxy terminus of yeast Atg13 binds phospholipid membrane via motifs that overlap with the Vac8-interacting domain

ABSTRACT

Macroautophagy/autophagy is a conserved catabolic recycling pathway involving the sequestration of cytoplasmic components within double-membrane vesicles termed autophagosomes. The autophagy-related (Atg) protein Atg13 is a key member of the autophagy initiation complex. The Atg13 C terminus is an intrinsically disordered region (IDR) harboring a binding site for the vacuolar membrane protein Vac8. Recent reports suggest Atg13 acts as a hub to assemble the initiation complex, and also participates in membrane recognition. Here we show that the Atg13 C terminus directly binds to lipid membranes via electrostatic interactions between positively charged residues in Atg13 and negatively charged phospholipids as well as a hydrophobic insertion of a Phe residue. We identified 2 sets of residues in the Atg13 IDR that affect its phospholipid-binding properties; these residues overlap with the Vac8-binding domain of Atg13. Our data indicate that Atg13 binding to phospholipids and Vac8 is mutually exclusive, and both are required for efficient autophagy.