The expanding view of protein-protein interactions: complexes involving intrinsically disordered proteins

AAontology: An Ontology of Amino Acid Scales for Interpretable Machine Learning

Amino acid scales are crucial for protein prediction tasks, many of them being curated in the AAindex database. Despite various clustering attempts to organize them and to better understand their relationships, these approaches lack the fine-grained classification necessary for satisfactory interpretability in many protein prediction problems. To address this issue, we developed AAontology-a two-level classification for 586 amino acid scales (mainly from AAindex) together with an in-depth analysis of their relations-using bag-of-word-based classification, clustering, and manual refinement over multiple iterations. AAontology organizes physicochemical scales into 8 categories and 67 subcategories, enhancing the interpretability of scale-based machine learning methods in protein bioinformatics. Thereby it enables researchers to gain a deeper biological insight. We anticipate that AAontology will be a building block to link amino acid properties with protein function and dysfunctions as well as aid informed decision-making in mutation analysis or protein drug design.

Chaotic aging: intrinsically disordered proteins in aging-related processes

The development of aging is associated with the disruption of key cellular processes manifested as well-established hallmarks of aging. Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) have no stable tertiary structure that provide them a power to be configurable hubs in signaling cascades and regulate many processes, potentially including those related to aging. There is a need to clarify the roles of IDPs/IDRs in aging. The dataset of 1702 aging-related proteins was collected from established aging databases and experimental studies. There is a noticeable presence of IDPs/IDRs, accounting for about 36% of the aging-related dataset, which is however less than the disorder content of the whole human proteome (about 40%). A Gene Ontology analysis of the used here aging proteome reveals an abundance of IDPs/IDRs in one-third of aging-associated processes, especially in genome regulation. Signaling pathways associated with aging also contain IDPs/IDRs on different hierarchical levels, revealing the importance of "structure-function continuum" in aging. Protein-protein interaction network analysis showed that IDPs present in different clusters associated with different aging hallmarks. Protein cluster with IDPs enrichment has simultaneously high liquid-liquid phase separation (LLPS) probability, "nuclear" localization and DNA-associated functions, related to aging hallmarks: genomic instability, telomere attrition, epigenetic alterations, and stem cells exhaustion. Intrinsic disorder, LLPS, and aggregation propensity should be considered as features that could be markers of pathogenic proteins. Overall, our analyses indicate that IDPs/IDRs play significant roles in aging-associated processes, particularly in the regulation of DNA functioning. IDP aggregation, which can lead to loss of function and toxicity, could be critically harmful to the cell. A structure-based analysis of aging and the identification of proteins that are particularly susceptible to disturbances can enhance our understanding of the molecular mechanisms of aging and open up new avenues for slowing it down.

Current research on viral proteins that interact with fibrillarin

The nucleolus is a multifunctional nuclear domain primarily dedicated to ribosome biogenesis. Certain viruses developed strategies to manipulate host nucleolar proteins to facilitate their replication by modulating ribosomal RNA (rRNA) processing. This association interferes with nucleolar functions resulting in overactivation or arrest of ribosome biogenesis, induction or inhibition of apoptosis, and affecting stress response. The nucleolar protein fibrillarin (FBL) is an important target of some plant and animal viruses. FBL is an essential and highly conserved S-adenosyl methionine (SAM) dependent methyltransferase, capable of rRNA degradation by its intrinsically disordered region (IDR), the glycine/arginine-rich (GAR) domain. It forms a ribonucleoprotein complex that directs 2'-O-methylations in more than 100 sites of pre-rRNAs. It is involved in multiple cellular processes, including initiation of transcription, oncogenesis, and apoptosis, among others. The interaction with animal viruses, including human viruses, triggered its redistribution to the nucleoplasm and cytoplasm, interfering with its role in pre-rRNA processing. Viral-encoded proteins with IDRs as nucleocapsids, matrix, Tat protein, and even a viral snoRNA, can associate with FBL, forcing the nucleolar protein to undergo atypical functions. Here we review the molecular mechanisms employed by animal and human viruses to usurp FBL functions and the effect on cellular processes, particularly in ribosome biogenesis.

Calpain cleavage of Junctophilin-2 generates a spectrum of calcium-dependent cleavage products and DNA-rich NT1-fragment domains in cardiomyocytes

Calpains are calcium-activated neutral proteases involved in the regulation of key signaling pathways. Junctophilin-2 (JP2) is a Calpain-specific proteolytic target and essential structural protein inside Ca²⁺ release units required for excitation-contraction coupling in cardiomyocytes. While downregulation of JP2 by Calpain cleavage in heart failure has been reported, the precise molecular identity of the Calpain cleavage sites and the (patho-)physiological roles of the JP2 proteolytic products remain controversial. We systematically analyzed the JP2 cleavage fragments as function of Calpain-1 versus Calpain-2 proteolytic activities, revealing that both Calpain isoforms preferentially cleave mouse JP2 at R565, but subsequently at three additional secondary Calpain cleavage sites. Moreover, we identified the Calpain-specific primary cleavage products for the first time in human iPSC-derived cardiomyocytes. Knockout of RyR2 in hiPSC-cardiomyocytes destabilized JP2 resulting in an increase of the Calpain-specific cleavage fragments. The primary N-terminal cleavage product NT₁ accumulated in the nucleus of mouse and human cardiomyocytes in a Ca²⁺-dependent manner, closely associated with euchromatic chromosomal regions, where NT₁ is proposed to function as a cardio-protective transcriptional regulator in heart failure. Taken together, our data suggest that stabilizing NT₁ by preventing secondary cleavage events by Calpain and other proteases could be an important therapeutic target for future studies.

Proteins, Small Peptides and Other Signaling Molecules Identified as Inconspicuous but Possibly Important Players in Microspores Reprogramming Toward Embryogenesis

In this review, we describe and integrate the latest knowledge on the signaling role of proteins and peptides in the stress-induced microspore embryogenesis (ME) in some crop plants with agricultural importance (i.e., oilseed rape, tobacco, barley, wheat, rice, triticale, rye). Based on the results received from the most advanced omix analyses, we have selected some inconspicuous but possibly important players in microspores reprogramming toward embryogenic development. We provide an overview of the roles and downstream effect of stress-related proteins (e.g., β-1,3-glucanases, chitinases) and small signaling peptides, especially cysteine—(e.g., glutathione, γ-thionins, rapid alkalinization factor, lipid transfer, phytosulfokine) and glycine-rich peptides and other proteins (e.g., fasciclin-like arabinogalactan protein) on acclimation ability of microspores and the cell wall reconstruction in a context of ME induction and haploids/doubled haploids (DHs) production. Application of these molecules, stimulating the induction and proper development of embryo-like structures and green plant regeneration, brings significant improvement of the effectiveness of DHs procedures and could result in its wider incorporation on a commercial scale. Recent advances in the design and construction of synthetic peptides–mainly cysteine-rich peptides and their derivatives–have accelerated the development of new DNA-free genome-editing techniques. These new systems are evolving incredibly fast and soon will find application in many areas of plant science and breeding.

Exploring Manually Curated Annotations of Intrinsically Disordered Proteins with DisProt

DisProt is the major repository of manually curated data for intrinsically disordered proteins collected from the literature. Although lacking a stable tertiary structure under physiological conditions, intrinsically disordered proteins carry out a plethora of biological functions, some of them directly arising from their flexible nature. A growing number of scientific studies have been published during the last few decades in an effort to shed light on their unstructured state, their binding modes, and their functions. DisProt makes use of a team of expert biocurators to provide up-to-date annotations of intrinsically disordered proteins from the literature, making them available to the scientific community. Here we present a comprehensive description on how to use DisProt in different contexts and provide a detailed explanation of how to explore and interpret manually curated annotations of intrinsically disordered proteins. We describe how to search DisProt annotations, using both the web interface and the API for programmatic access. Finally, we explain how to visualize and interpret a DisProt entry, p53, a widely studied protein characterized by the presence of unstructured N-terminal and C-terminal regions. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Performing a search in DisProt Support Protocol 1: Downloading options Support Protocol 2: Programmatic access with DisProt REST API Basic Protocol 2: Visualizing and interpreting DisProt entries: the p53 use case Basic Protocol 3: Providing feedback and submitting new intrinsic disorder-related data.

The Molecular Tweezer CLR01 Stabilizes a Disordered Protein-Protein Interface

Protein regions that are involved in protein-protein interactions (PPIs) very often display a high degree of intrinsic disorder, which is reduced during the recognition process. A prime example is binding of the rigid 14-3-3 adapter proteins to their numerous partner proteins, whose recognition motifs undergo an extensive disorder-to-order transition. In this context, it is highly desirable to control this entropy-costly process using tailored stabilizing agents. This study reveals how the molecular tweezer CLR01 tunes the 14-3-3/Cdc25CpS216 protein-protein interaction. Protein crystallography, biophysical affinity determination and biomolecular simulations unanimously deliver a remarkable finding: a supramolecular "Janus" ligand can bind simultaneously to a flexible peptidic PPI recognition motif and to a well-structured adapter protein. This binding fills a gap in the protein-protein interface, "freezes" one of the conformational states of the intrinsically disordered Cdc25C protein partner and enhances the apparent affinity of the interaction. This is the first structural and functional proof of a supramolecular ligand targeting a PPI interface and stabilizing the binding of an intrinsically disordered recognition motif to a rigid partner protein.

Under-folded proteins: Conformational ensembles and their roles in protein folding, function, and pathogenesis

For decades, protein function was intimately linked to the presence of a unique, aperiodic crystal‐like structure in a functional protein. The two only places for conformational ensembles of under‐folded (or partially folded) protein forms in this picture were either the end points of the protein denaturation processes or transiently populated folding intermediates. Recent years witnessed dramatic change in this perception and conformational ensembles, which the under‐folded proteins are, have moved from the shadow. Accumulated to date data suggest that a protein can exist in at least three global forms–functional and folded, functional and intrinsically disordered (nonfolded), and nonfunctional and misfolded/aggregated. Under‐folded protein states are crucial for each of these forms, serving as important folding intermediates of ordered proteins, or as functional states of intrinsically disordered proteins (IDPs) and IDP regions (IDPRs), or as pathology triggers of misfolded proteins. Based on these observations, conformational ensembles of under‐folded proteins can be classified as transient (folding and misfolding intermediates) and permanent (IDPs and stable misfolded proteins). Permanently under‐folded proteins can further be split into intentionally designed (IDPs and IDPRs) and unintentionally designed (misfolded proteins). Although intrinsic flexibility, dynamics, and pliability are crucial for all under‐folded proteins, the different categories of under‐foldedness are differently encoded in protein amino acid sequences. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 870–887, 2013.

The Proteasome Subunit Rpn8 Interacts with the Small Nucleolar RNA Protein (snoRNP) Assembly Protein Pih1 and Mediates Its Ubiquitin-independent Degradation in Saccharomyces cerevisiae

Pih1 is a scaffold protein of the Rvb1-Rvb2-Tah1-Pih1 (R2TP) protein complex, which is conserved in fungi and animals. The chaperone-like activity of the R2TP complex has been implicated in the assembly of multiple protein complexes, such as the small nucleolar RNA protein complex. However, the mechanism of the R2TP complex activity in vivo and the assembly of the complex itself are still largely unknown. Pih1 is an unstable protein and tends to aggregate when expressed alone. The C-terminal fragment of Pih1 contains multiple destabilization factors and acts as a degron when fused to other proteins. In this study, we investigated Pih1 interactors and identified a specific interaction between Pih1 and the proteasome subunit Rpn8 in yeast Saccharomyces cerevisiae when HSP90 co-chaperone Tah1 is depleted. By analyzing truncation mutants, we identified that the C-terminal 30 amino acids of Rpn8 are sufficient for the binding to Pih1 C terminus. With in vitro and in vivo degradation assays, we showed that the Pih1 C-terminal fragment Pih1(282-344) is able to induce a ubiquitin-independent degradation of GFP. Additionally, we demonstrated that truncation of the Rpn8 C-terminal disordered region does not affect proteasome assembly but specifically inhibits the degradation of the GFP-Pih1(282-344) fusion protein in vivo and Pih1 in vitro We propose that Pih1 is a ubiquitin-independent proteasome substrate, and the direct interaction with Rpn8 C terminus mediates its proteasomal degradation.

Recognition of Poly-Ubiquitins by the Proteasome through Protein Refolding Guided by Electrostatic and Hydrophobic Interactions

Specificity of protein degradation by cellular proteasomes comes from tetra-ubiquitin recognition. We carry out molecular dynamics simulations to characterize how the ubiquitin receptor Rpn10 recognizes in the 26S proteasome K48-linked tetra-ubiquitin. In the binding pose, ubiquitin and Rpn10 interact primarily through hydrophobic patches. However, K48-linked tetra-ubiquitin mostly assumes a closed form in solution prior to binding, in which its hydrophobic patches are not exposed to solvent. Likewise, the hydrophobic ubiquitin interacting motifs (UIMs) of Rpn10 are mostly protected prior to binding. As a result, ubiquitin recognition in the proteasome requires refolding of both K48-linked tetra-ubiquitin and Rpn10. Simulations suggest that conserved complementary electrostatic patterns of Rpn10 and ubiquitins guide protein association (stage 1 in the recognition process), which induces refolding (stage 2), and then facilitates formation of hydrophobic contacts (stage 3). The simulations also explain why Rpn10 has a higher affinity for K48-linked tetra-ubiquitin than for mono-ubiquitin and K48-linked di- and tri-ubiquitins. Simulation results expand on the current view that the flexible arm of Rpn10 acts as an extended fragment of α-helices and flexible coils in the recognition process.

High-resolution structural characterization of Noxa, an intrinsically disordered protein, by microsecond molecular dynamics simulations

High-resolution characterization of the structure and dynamics of intrinsically disordered proteins (IDPs) remains a challenging task. Consequently, a detailed understanding of the structural and functional features of IDPs remains limited, as very few full-length disordered proteins have been structurally characterized. We have performed microsecond-long molecular dynamics (MD) simulations of Noxa, the smallest member of the large Bcl-2 family of apoptosis regulating proteins, to characterize in atomic-level detail the structural features of a disordered protein. A 2.5 μs MD simulation starting from an unfolded state of the protein revealed the formation of a central antiparallel β-sheet structure flanked by two disordered segments at the N- and C-terminal ends. This topology is in reasonable agreement with protein disorder predictions and available experimental data. We show that this fold plays an essential role in the intracellular function and regulation of Noxa. We demonstrate that unbiased MD simulations in combination with a modern force field reveal structural and functional features of disordered proteins at atomic-level resolution.

Structure-based Inhibitor Design for the Intrinsically Disordered Protein c-Myc

Intrinsically disordered proteins (IDPs) are associated with various diseases and have been proposed as promising drug targets. However, conventional structure-based approaches cannot be applied directly to IDPs, due to their lack of ordered structures. Here, we describe a novel computational approach to virtually screen for compounds that can simultaneously bind to different IDP conformations. The test system used c-Myc, an oncoprotein containing a disordered basic helix-loop-helix-leucine zipper (bHLH-LZ) domain that adopts a helical conformation upon binding to Myc-associated factor X (Max). For the virtual screen, we used three binding pockets in representative conformations of c-Myc_370–409, which is part of the disordered bHLH-LZ domain. Seven compounds were found to directly bind c-Myc_370–409in vitro, and four inhibited the growth of the c-Myc-overexpressing cells by affecting cell cycle progression. Our approach of IDP conformation sampling, binding site identification, and virtual screening for compounds that can bind to multiple conformations provides a useful strategy for structure-based drug discovery targeting IDPs.

Application of linker technique to trap transiently interacting protein complexes for structural studies

Protein-protein interactions are key events controlling several biological processes. We have developed and employed a method to trap transiently interacting protein complexes for structural studies using glycine-rich linkers to fuse interacting partners, one of which is unstructured. Initial steps involve isothermal titration calorimetry to identify the minimum binding region of the unstructured protein in its interaction with its stable binding partner. This is followed by computational analysis to identify the approximate site of the interaction and to design an appropriate linker length. Subsequently, fused constructs are generated and characterized using size exclusion chromatography and dynamic light scattering experiments. The structure of the chimeric protein is then solved by crystallization, and validated both in vitro and in vivo by substituting key interacting residues of the full length, unlinked proteins with alanine. This protocol offers the opportunity to study crucial and currently unattainable transient protein interactions involved in various biological processes.

Rational design of antibodies targeting specific epitopes within intrinsically disordered proteins

Antibodies are powerful tools in life sciences research, as well as in diagnostic and therapeutic applications, because of their ability to bind given molecules with high affinity and specificity. Using current methods, however, it is laborious and sometimes difficult to generate antibodies to target specific epitopes within a protein, in particular if these epitopes are not effective antigens. Here we present a method to rationally design antibodies to enable them to bind virtually any chosen disordered epitope in a protein. The procedure consists in the sequence-based design of one or more complementary peptides targeting a selected disordered epitope and the subsequent grafting of such peptides on an antibody scaffold. We illustrate the method by designing six single-domain antibodies to bind different epitopes within three disease-related intrinsically disordered proteins and peptides (α-synuclein, Aβ42, and IAPP). Our results show that all these designed antibodies bind their targets with good affinity and specificity. As an example of an application, we show that one of these antibodies inhibits the aggregation of α-synuclein at substoichiometric concentrations and that binding occurs at the selected epitope. Taken together, these results indicate that the design strategy that we propose makes it possible to obtain antibodies targeting given epitopes in disordered proteins or protein regions.

Gene Model Annotations for Drosophila melanogaster: The Rule-Benders

In the context of the FlyBase annotated gene models in Drosophila melanogaster, we describe the many exceptional cases we have curated from the literature or identified in the course of FlyBase analysis. These range from atypical but common examples such as dicistronic and polycistronic transcripts, noncanonical splices, trans-spliced transcripts, noncanonical translation starts, and stop-codon readthroughs, to single exceptional cases such as ribosomal frameshifting and HAC1-type intron processing. In FlyBase, exceptional genes and transcripts are flagged with Sequence Ontology terms and/or standardized comments. Because some of the rule-benders create problems for handlers of high-throughput data, we discuss plans for flagging these cases in bulk data downloads.

Protein flexibility in the light of structural alphabets

Protein structures are valuable tools to understand protein function. Nonetheless, proteins are often considered as rigid macromolecules while their structures exhibit specific flexibility, which is essential to complete their functions. Analyses of protein structures and dynamics are often performed with a simplified three-state description, i.e., the classical secondary structures. More precise and complete description of protein backbone conformation can be obtained using libraries of small protein fragments that are able to approximate every part of protein structures. These libraries, called structural alphabets (SAs), have been widely used in structure analysis field, from definition of ligand binding sites to superimposition of protein structures. SAs are also well suited to analyze the dynamics of protein structures. Here, we review innovative approaches that investigate protein flexibility based on SAs description. Coupled to various sources of experimental data (e.g., B-factor) and computational methodology (e.g., Molecular Dynamic simulation), SAs turn out to be powerful tools to analyze protein dynamics, e.g., to examine allosteric mechanisms in large set of structures in complexes, to identify order/disorder transition. SAs were also shown to be quite efficient to predict protein flexibility from amino-acid sequence. Finally, in this review, we exemplify the interest of SAs for studying flexibility with different cases of proteins implicated in pathologies and diseases.

The intrinsic disorder alphabet. III. Dual personality of serine

Proteins are natural polypeptides consisting of 20 major amino acid residues, content and order of which in a given amino acid sequence defines the ability of a related protein to fold into unique functional state or to stay intrinsically disordered. Amino acid sequences code for both foldable (ordered) proteins/domains and for intrinsically disordered proteins (IDPs) and IDP regions (IDPRs), but these sequence codes are dramatically different. This difference starts with a very general property of the corresponding amino acid sequences, namely, their compositions. IDPs/IDPRs are enriched in specific disorder-promoting residues, whereas amino acid sequences of ordered proteins/domains typically contain more order-promoting residues. Therefore, the relative abundances of various amino acids in ordered and disordered proteins can be used to scale amino acids according to their disorder promoting potentials. This review continues a series of publications on the roles of different amino acids in defining the phenomenon of protein intrinsic disorder and represents serine, which is the third most disorder-promoting residue. Similar to previous publications, this review represents some physico-chemical properties of serine and the roles of this residue in structures and functions of ordered proteins, describes major posttranslational modifications tailored to serine, and finally gives an overview of roles of serine in structure and functions of intrinsically disordered proteins.

The N-Terminal Domain of SIRT1 Is a Positive Regulator of Endogenous SIRT1-Dependent Deacetylation and Transcriptional Outputs

The NAD⁺-dependent protein deacetylase SIRT1 regulates energy metabolism, responses to stress, and aging by deacetylating many different proteins, including histones and transcription factors. The mechanisms controlling SIRT1 enzymatic activity are complex and incompletely characterized, yet essential for understanding how to develop therapeutics that target SIRT1. Here, we demonstrate that the N-terminal domain of SIRT1 (NTERM) can trans-activate deacetylation activity by physically interacting with endogenous SIRT1 and promoting its association with the deacetylation substrate NF-κB p65. Two motifs within the NTERM domain contribute to activation of SIRT1-dependent activities, and expression of one of these motifs in mice is sufficient to lower fasting glucose levels and improve glucose tolerance in a manner similar to overexpression of SIRT1. Our results provide insights into the regulation of SIRT1 activity and a rationale for pharmacological control of SIRT1-dependent activities.

Determining and visualizing flexibility in protein structures

How to compare the structures of an ensemble of protein conformations is a fundamental problem in structural biology. As has been previously observed, the widely used RMSD measure due to Kabsch, in which a rigid-body superposition minimizing the least-squares positional deviations is performed, has its drawbacks when comparing and visualizing a set of flexible protein structures. Here, we develop a method, fleximatch, of protein structure comparison that takes flexibility into account. Based on a distance matrix measure of flexibility, a weighted superposition of distance matrices rather than of atomic coordinates is performed. Subsequently, this allows a consistent determination of (a) a superposition of structures for visualization, (b) a partitioning of the protein structure into rigid molecular components (core atoms), and (c) an atomic mobility measure. The method is suitable for highlighting both particularly flexible and rigid parts of a protein from structures derived from NMR, X-ray diffraction or molecular simulation.

Specificity and affinity quantification of flexible recognition from underlying energy landscape topography

Flexibility in biomolecular recognition is essential and critical for many cellular activities. Flexible recognition often leads to moderate affinity but high specificity, in contradiction with the conventional wisdom that high affinity and high specificity are coupled. Furthermore, quantitative understanding of the role of flexibility in biomolecular recognition is still challenging. Here, we meet the challenge by quantifying the intrinsic biomolecular recognition energy landscapes with and without flexibility through the underlying density of states. We quantified the thermodynamic intrinsic specificity by the topography of the intrinsic binding energy landscape and the kinetic specificity by association rate. We found that the thermodynamic and kinetic specificity are strongly correlated. Furthermore, we found that flexibility decreases binding affinity on one hand, but increases binding specificity on the other hand, and the decreasing or increasing proportion of affinity and specificity are strongly correlated with the degree of flexibility. This shows more (less) flexibility leads to weaker (stronger) coupling between affinity and specificity. Our work provides a theoretical foundation and quantitative explanation of the previous qualitative studies on the relationship among flexibility, affinity and specificity. In addition, we found that the folding energy landscapes are more funneled with binding, indicating that binding helps folding during the recognition. Finally, we demonstrated that the whole binding-folding energy landscapes can be integrated by the rigid binding and isolated folding energy landscapes under weak flexibility. Our results provide a novel way to quantify the affinity and specificity in flexible biomolecular recognition.

Flexibility in biomolecular recognition is crucial for the function. Flexibility often leads to moderate binding affinity but high binding specificity, challenging the conventional wisdom that high specificity is guaranteed by high affinity. Currently, understanding of the relationship between affinity and specificity in flexible biomolecular recognition is still obscure, even in a qualitative way. By exploring the intrinsic biomolecular recognition energy landscapes, we provided a novel way to quantify the thermodynamic intrinsic specificity by energy landscape topography and kinetic specificity by association rate. We show quantitatively that flexibility decreases binding affinity while increases binding specificity, and the relative changes in affinity and specificity are strongly correlated with the degree of flexibility. Our results show that more (less) flexibility leads to weaker (stronger) coupling between affinity and specificity. Importantly, we demonstrated that flexibility modulates affinity and specificity through the underlying energy landscape. Our study establishes the quantitative relationship among flexibility, affinity and specificity, bridging the gap between theory and experiments.

Discrete molecular dynamics can predict helical prestructured motifs in disordered proteins

Intrinsically disordered proteins (IDPs) lack a stable tertiary structure, but their short binding regions termed Pre-Structured Motifs (PreSMo) can form transient secondary structure elements in solution. Although disordered proteins are crucial in many biological processes and designing strategies to modulate their function is highly important, both experimental and computational tools to describe their conformational ensembles and the initial steps of folding are sparse. Here we report that discrete molecular dynamics (DMD) simulations combined with replica exchange (RX) method efficiently samples the conformational space and detects regions populating α-helical conformational states in disordered protein regions. While the available computational methods predict secondary structural propensities in IDPs based on the observation of protein-protein interactions, our ab initio method rests on physical principles of protein folding and dynamics. We show that RX-DMD predicts α-PreSMos with high confidence confirmed by comparison to experimental NMR data. Moreover, the method also can dissect α-PreSMos in close vicinity to each other and indicate helix stability. Importantly, simulations with disordered regions forming helices in X-ray structures of complexes indicate that a preformed helix is frequently the binding element itself, while in other cases it may have a role in initiating the binding process. Our results indicate that RX-DMD provides a breakthrough in the structural and dynamical characterization of disordered proteins by generating the structural ensembles of IDPs even when experimental data are not available.

From sequence and forces to structure, function, and evolution of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs), which lack persistent structure, are a challenge to structural biology due to the inapplicability of standard methods for characterization of folded proteins as well as their deviation from the dominant structure/function paradigm. Their widespread presence and involvement in biological function, however, has spurred the growing acceptance of the importance of IDPs and the development of new tools for studying their structure, dynamics, and function. The interplay of folded and disordered domains or regions for function and the existence of a continuum of protein states with respect to conformational energetics, motional timescales, and compactness are shaping a unified understanding of structure-dynamics-disorder/function relationships. In the 20(th) anniversary of Structure, we provide a historical perspective on the investigation of IDPs and summarize the sequence features and physical forces that underlie their unique structural, functional, and evolutionary properties.

Remarkably fast coupled folding and binding of the intrinsically disordered transactivation domain of cMyb to CBP KIX

Association rates for interactions between folded proteins have been investigated extensively, allowing the development of computational and theoretical prediction methods. Less is known about association rates for complexes where one or more partner is initially disordered, despite much speculation about how they may compare to those for folded proteins. We have attached a fluorophore to the N-terminus of the 25 amino acid cMyb peptide used previously in NMR and equilibrium studies (termed FITC-cMyb), and used this to monitor the kinetics of its interaction with the KIX protein. We have investigated the ionic strength and temperature dependence of the kinetics, and conclude that the association process is extremely fast, apparently exceeding the rates predicted by formulations applicable to interactions between pairs of folded proteins. This is despite the fact that not all collisions result in complex formation (there is an observable activation energy for the association process). We propose that this is partially a result of the disordered nature of the FITC-cMyb peptide itself.

Effects of phosphorylation on the structure and backbone dynamics of the intrinsically disordered connexin43 C-terminal domain

Phosphorylation of the connexin43 C-terminal (Cx43CT) domain regulates gap junction intercellular communication. However, an understanding of the mechanisms by which phosphorylation exerts its effects is lacking. Here, we test the hypothesis that phosphorylation regulates Cx43 gap junction intercellular communication by mediating structural changes in the C-terminal domain. Circular dichroism and nuclear magnetic resonance were used to characterize the effects of phosphorylation on the secondary structure and backbone dynamics of soluble and membrane-tethered Cx43CT domains. Cx43CT phospho-mimetic isoforms, which have Asp substitutions at specific Ser/Tyr sites, revealed phosphorylation alters the α-helical content of the Cx43CT domain only when attached to the membrane. The changes in secondary structure are due to variations in the conformational preference and backbone flexibility of residues adjacent and distal to the site(s) of modification. In addition to the known direct effects of phosphorylation on molecular partner interactions, the data presented here suggest phosphorylation may also indirectly regulate binding affinity by altering the conformational preference of the Cx43CT domain.

The unstructured domain of colicin N kills Escherichia coli

Bacteria often produce toxins which kill competing bacteria. Colicins, produced by and toxic to Escherichia coli bacteria are three-domain proteins so efficient that one molecule can kill a cell. The C-terminal domain carries the lethal activity and the central domain is required for surface receptor binding. The N-terminal domain, required for translocation across the outer membrane, is always intrinsically unstructured. It has always been assumed therefore that the C-terminal cytotoxic domain is required for the bactericidal activity. Here we report the unexpected finding that in isolation, the 90-residue unstructured N-terminal domain of colicin N is cytotoxic. Furthermore it causes ion leakage from cells but, unlike known antimicrobial peptides (AMPs) with this property, shows no membrane binding behaviour. Finally, its activity remains strictly dependent upon the same receptor proteins (OmpF and TolA) used by full-length colicin N. This mechanism of rapid membrane disruption, via receptor mediated binding of a soluble peptide, may reveal a new target for the development of highly specific antibacterials.

BeAtMuSiC: Prediction of changes in protein-protein binding affinity on mutations

The ability of proteins to establish highly selective interactions with a variety of (macro)molecular partners is a crucial prerequisite to the realization of their biological functions. The availability of computational tools to evaluate the impact of mutations on protein–protein binding can therefore be valuable in a wide range of industrial and biomedical applications, and help rationalize the consequences of non-synonymous single-nucleotide polymorphisms. BeAtMuSiC (http://babylone.ulb.ac.be/beatmusic) is a coarse-grained predictor of the changes in binding free energy induced by point mutations. It relies on a set of statistical potentials derived from known protein structures, and combines the effect of the mutation on the strength of the interactions at the interface, and on the overall stability of the complex. The BeAtMuSiC server requires as input the structure of the protein–protein complex, and gives the possibility to assess rapidly all possible mutations in a protein chain or at the interface, with predictive performances that are in line with the best current methodologies.

The alphabet of intrinsic disorder: II. Various roles of glutamic acid in ordered and intrinsically disordered proteins

The ability of a protein to fold into unique functional state or to stay intrinsically disordered is encoded in its amino acid sequence. Both ordered and intrinsically disordered proteins (IDPs) are natural polypeptides that use the same arsenal of 20 proteinogenic amino acid residues as their major building blocks. The exceptional structural plasticity of IDPs, their capability to exist as heterogeneous structural ensembles and their wide array of important disorder-based biological functions that complements functional repertoire of ordered proteins are all rooted within the peculiar differential usage of these building blocks by ordered proteins and IDPs. In fact, some residues (so-called disorder-promoting residues) are noticeably more common in IDPs than in sequences of ordered proteins, which, in their turn, are enriched in several order-promoting residues. Furthermore, residues can be arranged according to their “disorder promoting potencies,” which are evaluated based on the relative abundances of various amino acids in ordered and disordered proteins. This review continues a series of publications on the roles of different amino acids in defining the phenomenon of protein intrinsic disorder and concerns glutamic acid, which is the second most disorder-promoting residue.

The roles of conditional disorder in redox proteins

Cells are constantly exposed to various oxidants, either generated endogenously due to metabolic activity or exogenously. One way that cells respond to oxidants is through the action of redox-regulated proteins. These proteins also play important roles in oxidant signaling and protein biogenesis events. The key sensors built into redox-regulated proteins are cysteines, which undergo reversible thiol oxidation in response to changes in the oxidation status of the cellular environment. In this review, we discuss three examples of redox-regulated proteins found in bacteria, mitochondria, and chloroplasts. These proteins use oxidation of their redox-sensitive cysteines to reversibly convert large structural domains into more disordered regions or vice versa. These massive structural rearrangements are directly implicated in the functions of these proteins.

Native-based simulations of the binding interaction between RAP74 and the disordered FCP1 peptide

By dephosphorylating the C-terminal domain (CTD) of RNA polymerase II (Pol II), the Transcription Factor IIF (TFIIF)-associating CTD phosphatase (FCP1) performs an essential function in recycling Pol II for subsequent rounds of transcription. The interaction between FCP1 and TFIIF is mediated by the disordered C-terminal tail of FCP1, which folds to form an α-helix upon binding the RAP74 subunit of TFIIF. The present work reports a structure-based simulation study of this interaction between the folded winged-helix domain of RAP74 and the disordered C-terminal tail of FCP1. The comparison of measured and simulated chemical shifts suggests that the FCP1 peptide samples 40-60% of its native helical structure in the unbound disordered ensemble. Free energy calculations suggest that productive binding begins when RAP74 makes hydrophobic contacts with the C-terminal region of the FCP1 peptide. The FCP1 peptide then folds into an amphipathic helix by zipping up the binding interface. The relative plasticity of FCP1 results in a more cooperative binding mechanism, allows for a greater diversity of pathways leading to the bound complex, and may also eliminate the need for "backtracking" from contacts that form out of sequence.

Interaction of serum amyloid A with human cystatin C--identification of binding sites

Serum amyloid A (SAA) is a multifunctional acute-phase protein whose natural role seems to be participation in many physiologic and pathological processes. Prolonged increased SAA level in a number of chronic inflammatory and neoplastic diseases gives rise to reactive systemic amyloid A amyloidosis, where the N-terminal 76-amino acid residue-long segment of SAA is deposited as amyloid fibrils. Recently, a specific interaction between SAA and the ubiquitous inhibitor of cysteine proteases--human cystatin C (hCC)--has been described. Here, we report further evidence corroborating this interaction, and the identification of the SAA and hCC binding sites in the SAA-hCC complex, using a combination of selective proteolytic excision and high-resolution mass spectrometry. The shortest binding site in the SAA sequence was determined as SAA(86-104), whereas the binding site in hCC sequence was identified as hCC(96-102). Binding specificities of both interacting sequences were ascertained by affinity experiments (ELISA) and by registration of mass spectrum of SAA-hCC complex.

On the binding affinity of macromolecular interactions: daring to ask why proteins interact

Interactions between proteins are orchestrated in a precise and time-dependent manner, underlying cellular function. The binding affinity, defined as the strength of these interactions, is translated into physico-chemical terms in the dissociation constant (K_d), the latter being an experimental measure that determines whether an interaction will be formed in solution or not. Predicting binding affinity from structural models has been a matter of active research for more than 40 years because of its fundamental role in drug development. However, all available approaches are incapable of predicting the binding affinity of protein–protein complexes from coordinates alone. Here, we examine both theoretical and experimental limitations that complicate the derivation of structure–affinity relationships. Most work so far has concentrated on binary interactions. Systems of increased complexity are far from being understood. The main physico-chemical measure that relates to binding affinity is the buried surface area, but it does not hold for flexible complexes. For the latter, there must be a significant entropic contribution that will have to be approximated in the future. We foresee that any theoretical modelling of these interactions will have to follow an integrative approach considering the biology, chemistry and physics that underlie protein–protein recognition.

Intramolecular interactions stabilizing compact conformations of the intrinsically disordered kinase-inhibitor domain of Sic1: a molecular dynamics investigation

Cyclin-dependent kinase inhibitors (CKIs) are key regulatory proteins of the eukaryotic cell cycle, which modulate cyclin-dependent kinase (Cdk) activity. CKIs perform their inhibitory effect by the formation of ternary complexes with a target kinase and its cognate cyclin. These regulators generally belong to the class of intrinsically disordered proteins (IDPs), which lack a well-defined and organized three-dimensional (3D) structure in their free state, undergoing folding upon binding to specific partners. Unbound IDPs are not merely random-coil structures, but can present intrinsically folded structural units (IFSUs) and collapsed conformations. These structural features can be relevant to protein function in vivo. The yeast CKI Sic1 is a 284-amino acid IDP that binds to Cdk1 in complex with the Clb5,6 cyclins, preventing phosphorylation of G1 substrates and, therefore, entrance to the S phase. Sic1 degradation, triggered by multiple phosphorylation events, promotes cell-cycle progression. Previous experimental studies pointed out a propensity of Sic1 and its isolated domains to populate both extended and compact conformations. The present contribution provides models for compact conformations of the Sic1 kinase-inhibitory domain (KID) by all-atom molecular dynamics (MD) simulations in explicit solvent and in the absence of interactors. The results are integrated by spectroscopic and spectrometric data. Helical IFSUs are identified, along with networks of intramolecular interactions. The results identify a group of putative hub residues and networks of electrostatic interactions, which are likely to be involved in the stabilization of the globular states.

The stability of the small nucleolar ribonucleoprotein (snoRNP) assembly protein Pih1 in Saccharomyces cerevisiae is modulated by its C terminus

Pih1 is an unstable protein and a subunit of the R2TP complex that, in yeast Saccharomyces cerevisiae, also contains the helicases Rvb1, Rvb2, and the Hsp90 cofactor Tah1. Pih1 and the R2TP complex are required for the box C/D small nucleolar ribonucleoprotein (snoRNP) assembly and ribosomal RNA processing. Purified Pih1 tends to aggregate in vitro. Molecular chaperone Hsp90 and its cochaperone Tah1 are required for the stability of Pih1 in vivo. We had shown earlier that the C terminus of Pih1 destabilizes the protein and that the C terminus of Tah1 binds to the Pih1 C terminus to form a stable complex. Here, we analyzed the secondary structure of the Pih1 C terminus and identified two intrinsically disordered regions and five hydrophobic clusters. Site-directed mutagenesis indicated that one predicted intrinsically disordered region IDR2 is involved in Tah1 binding, and that the C terminus of Pih1 contains multiple destabilization or degron elements. Additionally, the Pih1 N-terminal domain, Pih1(1-230), was found to be able to complement the physiological role of full-length Pih1 at 37 °C. Pih1(1-230) as well as a shorter Pih1 N-terminal fragment Pih1(1-195) is able to bind Rvb1/Rvb2 heterocomplex. However, the sequence between the two disordered regions in Pih1 significantly enhances the Pih1 N-terminal domain binding to Rvb1/Rvb2. Based on these data, a model of protein-protein interactions within the R2TP complex is proposed.

Analysis of the relationships between evolvability, thermodynamics, and the functions of intrinsically disordered proteins/regions

The evolvability of proteins is not only restricted by functional and structural importance, but also by other factors such as gene duplication, protein stability, and an organism's robustness. Recently, intrinsically disordered proteins (IDPs)/regions (IDRs) have been suggested to play a role in facilitating protein evolution. However, the mechanisms by which this occurs remain largely unknown. To address this, we have systematically analyzed the relationship between the evolvability, stability, and function of IDPs/IDRs. Evolutionary analysis shows that more recently emerged IDRs have higher evolutionary rates with more functional constraints relaxed (or experiencing more positive selection), and that this may have caused accelerated evolution in the flanking regions and in the whole protein. A systematic analysis of observed stability changes due to single amino acid mutations in IDRs and ordered regions shows that while most mutations induce a destabilizing effect in proteins, mutations in IDRs cause smaller stability changes than in ordered regions. The weaker impact of mutations in IDRs on protein stability may have advantages for protein evolvability in the gain of new functions. Interestingly, however, an analysis of functional motifs in the PROSITE and ELM databases showed that motifs in IDRs are more conserved, characterized by smaller entropy and lower evolutionary rate, than in ordered regions. This apparently opposing evolutionary effect may be partly due to the flexible nature of motifs in IDRs, which require some key amino acid residues to engage in tighter interactions with other molecules. Our study suggests that the unique conformational and thermodynamic characteristics of IDPs/IDRs play an important role in the evolvability of proteins to gain new functions.

In vivo protein complex topologies: sights through a cross-linking lens

Proteins are a remarkable class of molecules that exhibit wide diversity of shapes or topological features that underpin protein interactions and give rise to biological function. In addition to quantitation of abundance levels of proteins in biological systems under a variety of conditions, the field of proteome research has as a primary mission the assignment of function for proteins and if possible, illumination of factors that enable function. For many years, chemical cross-linking methods have been used to provide structural data on single purified proteins and purified protein complexes. However, these methods also offer the alluring possibility to extend capabilities to complex biological samples such as cell lysates or intact living cells where proteins may exhibit native topological features that do not exist in purified form. Recent efforts are beginning to provide glimpses of protein complexes and topologies in cells that suggest continued development will yield novel capabilities to view functional topological features of many proteins and complexes as they exist in cells, tissues, or other complex samples. This review will describe rationale, challenges, and a few success stories along the path of development of cross-linking technologies for measurement of in vivo protein interaction topologies.

Virtual interactomics of proteins from biochemical standpoint

Virtual interactomics represents a rapidly developing scientific area on the boundary line of bioinformatics and interactomics. Protein-related virtual interactomics then comprises instrumental tools for prediction, simulation, and networking of the majority of interactions important for structural and individual reproduction, differentiation, recognition, signaling, regulation, and metabolic pathways of cells and organisms. Here, we describe the main areas of virtual protein interactomics, that is, structurally based comparative analysis and prediction of functionally important interacting sites, mimotope-assisted and combined epitope prediction, molecular (protein) docking studies, and investigation of protein interaction networks. Detailed information about some interesting methodological approaches and online accessible programs or databases is displayed in our tables. Considerable part of the text deals with the searches for common conserved or functionally convergent protein regions and subgraphs of conserved interaction networks, new outstanding trends and clinically interesting results. In agreement with the presented data and relationships, virtual interactomic tools improve our scientific knowledge, help us to formulate working hypotheses, and they frequently also mediate variously important in silico simulations.

Carbon-Detected (15)N NMR Spin Relaxation of an Intrinsically Disordered Protein: FCP1 Dynamics Unbound and in Complex with RAP74

Intrinsically disordered proteins (IDPs) lack unique 3D structures under native conditions and as such exist as highly dynamic ensembles in solution. We present two (13)C-direct detection experiments for the measurement of (15)N NMR spin relaxation called the CON(T1)-IPAP and CON(T2)-IPAP that quantify backbone dynamics on a per-residue basis for IDPs in solution. These experiments have been applied to the intrinsically disordered C-terminal of FCP1, both free in solution and while bound to the RAP74 winged-helix domain. The results provide evidence that most of FCP1 remains highly dynamic in both states, while the 20 residues forming direct contact with RAP74 become more ordered in the complex. Parallel analysis of RAP74 backbone (15)N NMR spin relaxation reveals only very limited ordering of RAP74 upon FCP1 binding. Taken together, these data show that folding-upon-binding is highly local in this system, with disorder prevailing even in the complex.

Backbone conformational preferences of an intrinsically disordered protein in solution

We have performed a 4-μs molecular dynamics simulation to investigate the native conformational preferences of the intrinsically disordered kinase-inducible domain (KID) of the transcription factor CREB in solution. There is solid experimental evidence showing that KID does not possess a bound-like structure in solution; however, it has been proposed that coil-to-helix transitions upon binding to its binding partner (CBP) are template-driven. While these studies indicate that IDPs possess a bias towards the bound structure, they do not provide direct evidence on the time-dependent conformational preferences of IDPs in atomic detail. Our simulation captured intrinsic conformational characteristics of KID that are in good agreement with experimental data such as a very small percentage of helical structure in its segment α(B) and structural disorder in solution. We used dihedral principal component analysis dPCA to map the conformations of KID in the microsecond timescale. By using principal components as reaction coordinates, we further constructed dPCA-based free energy landscapes of KID. Analysis of the free energy landscapes showed that KID is best characterized as a conformational ensemble of rapidly interconverting conformations. Interestingly, we found that despite the conformational heterogeneity of the backbone and the absence of substantial secondary structure, KID does not randomly sample the conformational space in solution: analysis of the (Φ, Ψ) dihedral angles showed that several individual residues of KID possess a strong bias toward the helical region of the Ramachandran plot. We suggest that the intrinsic conformational preferences of KID provide a bias toward the folded state without having to populate bound-like conformations before binding. Furthermore, we argue that these conformational preferences do not represent actual structural constraints which drive binding through a single pathway, which allows for specific interactions with multiple binding partners. Based on this evidence, we propose that the backbone conformational preferences of KID provide a thermodynamic advantage for folding and binding without negatively affecting the kinetics of binding. We further discuss the relation of our results to previous studies to rationalize the functional implications of the conformational preferences of IDPs, such as the optimization of structural disorder in protein-protein interactions. This study illustrates the importance in obtaining atomistic information of intrinsically disordered proteins in real time to reveal functional features arising from their complex conformational space.

Inherent structural disorder and dimerisation of murine norovirus NS1-2 protein

Human noroviruses are highly infectious viruses that cause the majority of acute, non-bacterial epidemic gastroenteritis cases worldwide. The first open reading frame of the norovirus RNA genome encodes for a polyprotein that is cleaved by the viral protease into six non-structural proteins. The first non-structural protein, NS1-2, lacks any significant sequence similarity to other viral or cellular proteins and limited information is available about the function and biophysical characteristics of this protein. Bioinformatic analyses identified an inherently disordered region (residues 1–142) in the highly divergent N-terminal region of the norovirus NS1-2 protein. Expression and purification of the NS1-2 protein of Murine norovirus confirmed these predictions by identifying several features typical of an inherently disordered protein. These were a biased amino acid composition with enrichment in the disorder promoting residues serine and proline, a lack of predicted secondary structure, a hydrophilic nature, an aberrant electrophoretic migration, an increased Stokes radius similar to that predicted for a protein from the pre-molten globule family, a high sensitivity to thermolysin proteolysis and a circular dichroism spectrum typical of an inherently disordered protein. The purification of the NS1-2 protein also identified the presence of an NS1-2 dimer in Escherichia coli and transfected HEK293T cells. Inherent disorder provides significant advantages including structural flexibility and the ability to bind to numerous targets allowing a single protein to have multiple functions. These advantages combined with the potential functional advantages of multimerisation suggest a multi-functional role for the NS1-2 protein.

Intrinsic protein disorder in human pathways

We analyze human-specific KEGG pathways trying to understand the functional role of intrinsic disorder in proteins. Pathways provide a comprehensive picture of biological processes and allow better understanding of a protein's function within the specific context of its surroundings. Our study pinpoints a few specific pathways significantly enriched in disorder-containing proteins and identifies the role of these proteins within the framework of pathway relationships. Three major categories of relations are shown to be significantly enriched in disordered proteins: gene expression, protein binding and to a lesser degree, protein phosphorylation. Finally we find that relations involving protein activation and to some extent inhibition are characterized by low disorder content.

AH-DB: collecting protein structure pairs before and after binding

This work presents the Apo-Holo DataBase (AH-DB, http://ahdb.ee.ncku.edu.tw/ and http://ahdb.csbb.ntu.edu.tw/), which provides corresponding pairs of protein structures before and after binding. Conformational transitions are commonly observed in various protein interactions that are involved in important biological functions. For example, copper-zinc superoxide dismutase (SOD1), which destroys free superoxide radicals in the body, undergoes a large conformational transition from an 'open' state (apo structure) to a 'closed' state (holo structure). Many studies have utilized collections of apo-holo structure pairs to investigate the conformational transitions and critical residues. However, the collection process is usually complicated, varies from study to study and produces a small-scale data set. AH-DB is designed to provide an easy and unified way to prepare such data, which is generated by identifying/mapping molecules in different Protein Data Bank (PDB) entries. Conformational transitions are identified based on a refined alignment scheme to overcome the challenge that many structures in the PDB database are only protein fragments and not complete proteins. There are 746,314 apo-holo pairs in AH-DB, which is about 30 times those in the second largest collection of similar data. AH-DB provides sophisticated interfaces for searching apo-holo structure pairs and exploring conformational transitions from apo structures to the corresponding holo structures.

Is there a biological cost of protein disorder? Analysis of cancer-associated mutations

As many diseases can be traced back to altered protein function, studying the effect of genetic variations at the level of proteins can provide a clue to understand how changes at the DNA level lead to various diseases. Cellular processes rely not only on proteins with well-defined structure but can also involve intrinsically disordered proteins (IDPs) that exist as highly flexible ensembles of conformations. Disordered proteins are mostly involved in signaling and regulatory processes, and their functional repertoire largely complements that of globular proteins. However, it was also suggested that protein disorder entails an increased biological cost. This notion was supported by a set of individual IDPs involved in various diseases, especially in cancer, and the increased amount of disorder observed among disease-associated proteins. In this work, we tested if there is any biological risk associated with protein disorder at the level of single nucleotide mutations. Specifically, we analyzed the distribution of mutations within ordered and disordered segments. Our results demonstrated that while neutral polymorphisms were more likely to occur within disordered segments, cancer-associated mutations had a preference for ordered regions. Additionally, we proposed an alternative explanation for the association of protein disorder and the involvement in cancer with the consideration of functional annotations. Individual examples also suggested that although disordered segments are fundamental functional elements, their presence is not necessarily accompanied with an increased mutation rate in cancer. The presented study can help to understand how the different structural properties of proteins influence the consequences of genetic mutations.