Intrinsic Disorder in the Protein Data Bank

Abstract

The Protein Data Bank (PDB) is the preeminent source of protein structural information. PDB contains over 32,500 experimentally determined 3-D structures solved using X-ray crystallography or nuclear magnetic resonance spectroscopy. Intrinsically disordered regions fail to form a fixed 3-D structure under physiological conditions. In this study, we compare the amino-acid sequences of proteins whose structures are determined by X-ray crystallography with the corresponding sequences from the Swiss-Prot database. The analyzed dataset includes 16,370 structures, which represent 18,101 PDB chains and 5,434 different proteins from 910 different organisms (2,793 eukaryotic, 2,109 bacterial, 288 viral, and 244 archaeal). In this dataset, on average, each Swiss-Prot protein is represented by 7 PDB chains with 76% of the crystallized regions being represented by more than one structure. Intriguingly, the complete sequences of only ~7% of proteins are observed in the corresponding PDB structures, and only ~25% of the total dataset have >95% of their lengths observed in the corresponding PDB structures. This suggests that the vast majority of PDB proteins is shorter than their corresponding Swiss-Prot sequences and/or contain numerous residues, which are not observed in maps of electron density. To determine the prevalence of disordered regions in PDB, the residues in the Swiss-Prot sequences were grouped into four general categories, “Observed” (which correspond to structured regions), “Not observed” (regions with missing electron density, potentially disordered), “Uncharacterized,” and “Ambiguous,” depending on their appearance in the corresponding PDB entries. This non-redundant set of residues can be viewed as a ‘fragment’ or empirical domain database that contains a set of experimentally determined structured regions or domains and a set of experimentally verified disordered regions or domains. We studied the propensities and properties of residues in these four categories and analyzed their relations to the predictions of disorder using several algorithms. “Non-observed,” “Ambiguous,” and “Uncharacterized” regions were shown to possess the amino acid compositional biases typical of intrinsically disordered proteins. The application of four different disorder predictors (PONDR® VL-XT, VL3-BA, VSL1P, and IUPred) revealed that the vast majority of residues in the “Observed” dataset are ordered, and that the “Not observed” regions are mostly disordered. The “Uncharacterized” regions possess some tendency toward order, whereas the predictions for the short “Ambiguous” regions are really ambiguous. Long “Ambiguous” regions (>70 amino acid residues) are mostly predicted to be ordered, suggesting that they are likely to be “wobbly” domains.

Overall, we showed that completely ordered proteins are not highly abundant in PDB and many PDB sequences have disordered regions. In fact, in the analyzed dataset ~10% of the PDB proteins contain regions of consecutive missing or ambiguous residues longer than 30 amino-acids and ~40% of the proteins possess short regions (≥10 and <30 amino-acid long) of missing and ambiguous residues.     

Structure and Intrinsic Disorder in Protein Autoinhibition

1. Download : Download high-res image (184KB)
2. Download : Download full-size image

Highlights? Autoinhibitory regions in proteins are enriched in intrinsic disorder ? Modulation of intrinsic disorder contributes to the fine-tuning of autoinhibition ? Disordered autoinhibitory regions are more often phosphorylated and change structure ? Intrinsic disorder in inhibitory regions is exploited to fine-tune inhibition     

SPOT-Disorder2: Improved Protein Intrinsic Disorder Prediction by Ensembled Deep Learning

Intrinsically disordered or unstructured proteins (or regions in proteins) have been found to be important in a wide range of biological functions and implicated in many diseases. Due to the high cost and low efficiency of experimental determination of intrinsic disorder and the exponential increase of unannotated protein sequences, developing complementary computational prediction methods has been an active area of research for several decades. Here, we employed an ensemble of deep Squeeze-and-Excitation residual inception and long short-term memory (LSTM) networks for predicting protein intrinsic disorder with input from evolutionary information and predicted one-dimensional structural properties. The method, called SPOT-Disorder2, offers substantial and consistent improvement not only over our previous technique based on LSTM networks alone, but also over other state-of-the-art techniques in three independent tests with different ratios of disordered to ordered amino acid residues, and for sequences with either rich or limited evolutionary information. More importantly, semi-disordered regions predicted in SPOT-Disorder2 are more accurate in identifying molecular recognition features (MoRFs) than methods directly designed for MoRFs prediction. SPOT-Disorder2 is available as a web server and as a standalone program at https://sparks-lab.org/server/spot-disorder2/.     

Intrinsic Disorder in Pathogen Effectors: Protein Flexibility as an Evolutionary Hallmark in a Molecular Arms Race

Effector proteins represent a refined mechanism of bacterial pathogens to overcome plants’ innate immune systems. These modular proteins often manipulate host physiology by directly interfering with immune signaling of plant cells. Even if host cells have developed efficient strategies to perceive the presence of pathogenic microbes and to recognize intracellular effector activity, it remains an open question why only few effectors are recognized directly by plant resistance proteins. Based on in-silico genome-wide surveys and a reevaluation of published structural data, we estimated that bacterial effectors of phytopathogens are highly enriched in long-disordered regions (>50 residues). These structurally flexible segments have no secondary structure under physiological conditions but can fold in a stimulus-dependent manner (e.g., during protein–protein interactions). The high abundance of intrinsic disorder in effectors strongly suggests positive evolutionary selection of this structural feature and highlights the dynamic nature of these proteins. We postulate that such structural flexibility may be essential for (1) effector translocation, (2) evasion of the innate immune system, and (3) host function mimicry. The study of these dynamical regions will greatly complement current structural approaches to understand the molecular mechanisms of these proteins and may help in the prediction of new effectors.Plants and pathogens are entangled in a continual arms race. While host organisms have developed complex and dynamic immune systems able to recognize a wide range of pathogens and to discriminate them from beneficial microbes (Jones and Dangl, 2006; Medzhitov, 2007), bacterial pathogens have evolved refined adaptation strategies to overcome the plant’s innate immune system. Among these ingenious adaptations are effector proteins. Most of these proteins are secreted via the type III secretion system (TTSS) into the host cytoplasm, where they manipulate the immune signaling and the physiology of plant cells and thereby improve bacterial fitness within the host (Dean, 2011).Plant–pathogen interactions are highly dynamic processes, both from the evolutionary and the physiological point of view. Here, we postulate that they are equally dynamic at the protein-structure level. This is based on our finding that numerous effector proteins are predicted to be intrinsically disordered (ID) and that this feature may be essential for (1) effector translocation, (2) evasion of the innate immune system, and (3) host function mimicry. Intrinsic disorder has so far been postulated to preferentially occur in eukaryotic proteins. While on average ∼20% of the eukaryotic proteome harbors long (>50 residues) ID segments, these regions are only predicted at low abundance (8% on average) in bacterial proteomes (Dunker et al., 2000). The most likely reason for this discrepancy is the lack of efficient mechanisms to protect unfolded proteins from degradation (Ward et al., 2004). However, when surveying genomes of pathogenic bacteria with the widely used PONDR VL-XT program (Romero et al., 2001), we observed that not only the average percentage of sequence disorder, but most strikingly long (>50 residues) stretches of intrinsic disorder are highly overrepresented in secreted effectors, with especially high levels in phytopathogenic bacteria (Pseudomonas syringae, ∼39%; Ralstonia solanacearum, ∼70%; Xanthomonas spp, ∼77%) (Supplemental Table 1 online). This striking enrichment of unstructured regions strongly suggests positive evolutionary selection of intrinsic disorder in effector proteins and highlights their dynamic nature.

Table 1.

Predictions of Intrinsic Disorder in Effectors and Whole Proteomes of Different Bacterial Species

Intrinsic Disorder of the Neuronal SNARE Protein SNAP25a in its Pre-fusion Conformation

The neuronal SNARE protein SNAP25a (isoform 2) forms part of the SNARE complex eliciting synaptic vesicle fusion during neuronal exocytosis. While the post-fusion cis-SNARE complex has been studied extensively, little is known about the pre-fusion conformation of SNAP25a. Here we analyze monomeric SNAP25a by NMR spectroscopy, further supported by small-angle X-ray scattering (SAXS) experiments. SAXS data indicate that monomeric SNAP25 is more compact than a Gaussian chain but still a random coil. NMR shows that for monomeric SNAP25a, before SNAP25a interacts with its SNARE partners to drive membrane fusion, only the N-terminal part (region A5 to V36) of the first SNARE motif, SN1 (L11 - L81), is helical, comprising two α-helices (ranging from A5 to Q20 and S25 toV36). From E37 onwards, SNAP25a is mostly disordered and displays high internal flexibility, including the C-terminal part of SN1, almost the entire second SNARE motif (SN2, N144-A199), and the connecting loop region. Apart from the N-terminal helices, only the C-termini of both SN1 (E73 - K79) and SN2 (region T190 - A199), as well as two short regions in the connecting loop (D99 - K102 and E123 - M127) show a weak α-helical propensity (α-helical population < 25%). We speculate that the N-terminal helices (A5 to Q20 and S25 to V36) which constitute the N-terminus of SN1 act as a nucleation site for initiating SNARE zippering.     

Structural Analysis of Mitochondrial Mutations Reveals a Role for Bigenomic Protein Interactions in Human Disease

Mitochondria are the energy producing organelles of the cell, and mutations within their genome can cause numerous and often severe human diseases. At the heart of every mitochondrion is a set of five large multi-protein machines collectively known as the mitochondrial respiratory chain (MRC). This cellular machinery is central to several processes important for maintaining homeostasis within cells, including the production of ATP. The MRC is unique due to the bigenomic origin of its interacting proteins, which are encoded in the nucleus and mitochondria. It is this, in combination with the sheer number of protein-protein interactions that occur both within and between the MRC complexes, which makes the prediction of function and pathological outcome from primary sequence mutation data extremely challenging. Here we demonstrate how 3D structural analysis can be employed to predict the functional importance of mutations in mtDNA protein-coding genes. We mined the MITOMAP database and, utilizing the latest structural data, classified mutation sites based on their location within the MRC complexes III and IV. Using this approach, four structural classes of mutation were identified, including one underexplored class that interferes with nuclear-mitochondrial protein interactions. We demonstrate that this class currently eludes existing predictive approaches that do not take into account the quaternary structural organization inherent within and between the MRC complexes. The systematic and detailed structural analysis of disease-associated mutations in the mitochondrial Complex III and IV genes significantly enhances the predictive power of existing approaches and our understanding of how such mutations contribute to various pathologies. Given the general lack of any successful therapeutic approaches for disorders of the MRC, these findings may inform the development of new diagnostic and prognostic biomarkers, as well as new drugs and targets for gene therapy.     

Structural Insights into Alternate Aggregated Prion Protein Forms

The conversion of the cellular form of the prion protein (PrP^C) to an abnormal, alternatively folded isoform (PrP^Sc) is the central event in prion diseases or transmissible spongiform encephalopathies. Recent studies have demonstrated de novo generation of murine prions from recombinant prion protein (recPrP) after inoculation into transgenic and wild-type mice. These so-called synthetic prions lead to novel prion diseases with unique neuropathological and biochemical features. Moreover, the use of recPrP in an amyloid seeding assay can specifically detect and amplify various strains of prions. We employed this assay in our experiments and analyzed in detail the morphology of aggregate structures produced under defined chemical constraints. Our results suggest that changes in the concentration of guanidine hydrochloride can lead to different kinetic traces in a typical thioflavin T(ThT) assay. Morphological and structural analysis of these aggregates by atomic force microscopy indicates a variation in the structure of the PrP molecular assemblies.In particular, ThT positive PrP aggregates produced from rec mouse PrP residues 89 to 230 lead to mostly oligomeric structures at low concentrations of guanidine hydrochloride, while more amyloidal structures were observed at higher concentrations of the denaturant. These findings highlight the presence of numerous and complex pathways in deciphering prion constraints for infectivity and toxicity.     

Prediction of Spontaneous Protein Deamidation from Sequence-Derived Secondary Structure and Intrinsic Disorder

Asparagine residues in proteins undergo spontaneous deamidation, a post-translational modification that may act as a molecular clock for the regulation of protein function and turnover. Asparagine deamidation is modulated by protein local sequence, secondary structure and hydrogen bonding. We present NGOME, an algorithm able to predict non-enzymatic deamidation of internal asparagine residues in proteins in the absence of structural data, using sequence-based predictions of secondary structure and intrinsic disorder. Compared to previous algorithms, NGOME does not require three-dimensional structures yet yields better predictions than available sequence-only methods. Four case studies of specific proteins show how NGOME may help the user identify deamidation-prone asparagine residues, often related to protein gain of function, protein degradation or protein misfolding in pathological processes. A fifth case study applies NGOME at a proteomic scale and unveils a correlation between asparagine deamidation and protein degradation in yeast. NGOME is freely available as a webserver at the National EMBnet node Argentina, URL: http://www.embnet.qb.fcen.uba.ar/ in the subpage “Protein and nucleic acid structure and sequence analysis”.     

Structural Insights into the Interactions between Platelet Receptors and Fibrillar Collagen

Collagen peptides have been used to identify binding sites for several important collagen receptors, including integrin α₂β₁, glycoprotein VI, and von Willebrand factor. In parallel, the structures of these collagen receptors have been reported, and their interactions with collagen peptides have been studied. Recently, the three-dimensional structure of the intact type I collagen fiber from rat tail tendon has been resolved by fiber diffraction. It is now possible to map the binding sites of platelet collagen receptors onto the intact collagen fiber in three dimensions. This minireview will discuss these recent findings and their implications for platelet activation by collagen.     

Structural Insights into the Inhibition of Actin-Capping Protein by Interactions with Phosphatidic Acid and Phosphatidylinositol (4,5)-Bisphosphate

The actin cytoskeleton is a dynamic structure that coordinates numerous fundamental processes in eukaryotic cells. Dozens of actin-binding proteins are known to be involved in the regulation of actin filament organization or turnover and many of these are stimulus-response regulators of phospholipid signaling. One of these proteins is the heterodimeric actin-capping protein (CP) which binds the barbed end of actin filaments with high affinity and inhibits both addition and loss of actin monomers at this end. The ability of CP to bind filaments is regulated by signaling phospholipids, which inhibit the activity of CP; however, the exact mechanism of this regulation and the residues on CP responsible for lipid interactions is not fully resolved. Here, we focus on the interaction of CP with two signaling phospholipids, phosphatidic acid (PA) and phosphatidylinositol (4,5)-bisphosphate (PIP₂). Using different methods of computational biology such as homology modeling, molecular docking and coarse-grained molecular dynamics, we uncovered specific modes of high affinity interaction between membranes containing PA/phosphatidylcholine (PC) and plant CP, as well as between PIP₂/PC and animal CP. In particular, we identified differences in the binding of membrane lipids by animal and plant CP, explaining previously published experimental results. Furthermore, we pinpoint the critical importance of the C-terminal part of plant CPα subunit for CP–membrane interactions. We prepared a GST-fusion protein for the C-terminal domain of plant α subunit and verified this hypothesis with lipid-binding assays in vitro.     

Structural Insights into RbmA,a Biofilm Scaffolding Protein of V. Cholerae

V. cholerae can form sessile biofilms associated with abiotic surfaces, cyanobacteria, zoo-plankton, mollusks, or crustaceans. Along with the vibrio polysaccharide, secreted proteins of the rbm gene cluster are key to the biofilm ultrastructure. Here we provide a thorough structural characterization of RbmA, a protein involved in mediating cell-cell and cell-biofilm contacts. We correlate our structural findings with initial ligand specificity screening results, NMR protein-ligand interaction analysis, and complement our results with a full biocomputational study.     

Intrinsic Disorder in Human RNA-Binding Proteins

Although RNA-binding proteins (RBPs) are known to be enriched in intrinsic disorder, no previous analysis focused on RBPs interacting with specific RNA types. We fill this gap with a comprehensive analysis of the putative disorder in RBPs binding to six common RNA types: messenger RNA (mRNA), transfer RNA (tRNA), small nuclear RNA (snRNA), non-coding RNA (ncRNA), ribosomal RNA (rRNA), and internal ribosome RNA (irRNA). We also analyze the amount of putative intrinsic disorder in the RNA-binding domains (RBDs) and non-RNA-binding-domain regions (non-RBD regions). Consistent with previous studies, we show that in comparison with human proteome, RBPs are significantly enriched in disorder. However, closer examination finds significant enrichment in predicted disorder for the mRNA-, rRNA- and snRNA-binding proteins, while the proteins that interact with ncRNA and irRNA are not enriched in disorder, and the tRNA-binding proteins are significantly depleted in disorder. We show a consistent pattern of significant disorder enrichment in the non-RBD regions coupled with low levels of disorder in RBDs, which suggests that disorder is relatively rarely utilized in the RNA-binding regions. Our analysis of the non-RBD regions suggests that disorder harbors posttranslational modification sites and is involved in the putative interactions with DNA. Importantly, we utilize experimental data from DisProt and independent data from Pfam to validate the above observations that rely on the disorder predictions. This study provides new insights into the distribution of disorder across proteins that bind different RNA types and the functional role of disorder in the regions where it is enriched.     

Structural Analysis Using Chiroptical Spectroscopy: Insights and Cautions

Chiroptical spectroscopy has evolved into a promising tool for chiral molecular structural determination in the last four decades. Determination of the absolute configurations (ACs) of bromochlorofluoromethane and [²H₁,²H₂,²H₃]‐neopentane demonstrated the enviable advantages of chiroptical spectroscopy. Furthermore, uncovering the errors in the ACs reported in the literature established a glimpse of what can be accomplished with the modern chiroptical spectroscopic methods. Despite these triumphs, it is important to exercise caution in the practice of chiroptical spectroscopic methods, because certain widely practiced approaches can lead to erroneous conclusions. Selected major accomplishments and special precautions needed for future applications are emphasized. Chirality 28:445–452, 2016. © 2016 Wiley Periodicals, Inc.