Molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins are biomolecules that do not have a definite 3D structure; therefore, their dynamical simulation cannot start from a known list of atomistic positions, such as a Protein Data Bank file. We describe a method to start a computer simulation of these proteins. The first step of the procedure is the creation of a multi-rod configuration of the molecule, derived from its primary sequence. This structure is dynamically evolved in vacuo until its gyration radius reaches the experimental average value; at this point solvent molecules, in explicit or implicit implementation, are added to the protein and a regular molecular dynamics simulation follows. We have applied this procedure to the simulation of tau, one of the largest totally disordered proteins.     

Binding cavities and druggability of intrinsically disordered proteins

To assess the potential of intrinsically disordered proteins (IDPs) as drug design targets, we have analyzed the ligand-binding cavities of two datasets of IDPs (containing 37 and 16 entries, respectively) and compared their properties with those of conventional ordered (folded) proteins. IDPs were predicted to possess more binding cavity than ordered proteins at similar length, supporting the proposed advantage of IDPs economizing genome and protein resources. The cavity number has a wide distribution within each conformation ensemble for IDPs. The geometries of the cavities of IDPs differ from the cavities of ordered proteins, for example, the cavities of IDPs have larger surface areas and volumes, and are more likely to be composed of a single segment. The druggability of the cavities was examined, and the average druggable probability is estimated to be 9% for IDPs, which is almost twice that for ordered proteins (5%). Some IDPs with druggable cavities that are associated with diseases are listed. The optimism versus obstacles for drug design for IDPs is also briefly discussed.     

Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding

The sequence–structure–function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well‐defined three‐dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.     

Context-dependent resistance to proteolysis of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs), also known as intrinsically unstructured proteins (IUPs), lack a well-defined 3D structure in vitro and, in some cases, also in vivo. Here, we discuss the question of proteolytic sensitivity of IDPs, with a view to better explaining their in vivo characteristics. After an initial assessment of the status of IDPs in vivo, we briefly survey the intracellular proteolytic systems. Subsequently, we discuss the evidence for IDPs being inherently sensitive to proteolysis. Such sensitivity would not, however, result in enhanced degradation if the protease-sensitive sites were sequestered. Accordingly, IDP access to and degradation by the proteasome, the major proteolytic complex within eukaryotic cells, are discussed in detail. The emerging picture appears to be that IDPs are inherently sensitive to proteasomal degradation along the lines of the "degradation by default" model. However, available data sets of intracellular protein half-lives suggest that intrinsic disorder does not imply a significantly shorter half-life. We assess the power of available systemic half-life measurements, but also discuss possible mechanisms that could protect IDPs from intracellular degradation. Finally, we discuss the relevance of the proteolytic sensitivity of IDPs to their function and evolution.     

Resilience of death: intrinsic disorder in proteins involved in the programmed cell death

It is recognized now that intrinsically disordered proteins (IDPs), which do not have unique 3D structures as a whole or in noticeable parts, constitute a significant fraction of any given proteome. IDPs are characterized by an astonishing structural and functional diversity that defines their ability to be universal regulators of various cellular pathways. Programmed cell death (PCD) is one of the most intricate cellular processes where the cell uses specialized cellular machinery and intracellular programs to kill itself. This cell-suicide mechanism enables metazoans to control cell numbers and to eliminate cells that threaten the animal''s survival. PCD includes several specific modules, such as apoptosis, autophagy, and programmed necrosis (necroptosis). These modules are not only tightly regulated but also intimately interconnected and are jointly controlled via a complex set of protein–protein interactions. To understand the role of the intrinsic disorder in controlling and regulating the PCD, several large sets of PCD-related proteins across 28 species were analyzed using a wide array of modern bioinformatics tools. This study indicates that the intrinsic disorder phenomenon has to be taken into consideration to generate a complete picture of the interconnected processes, pathways, and modules that determine the essence of the PCD. We demonstrate that proteins involved in regulation and execution of PCD possess substantial amount of intrinsic disorder. We annotate functional roles of disorder across and within apoptosis, autophagy, and necroptosis processes. Disordered regions are shown to be implemented in a number of crucial functions, such as protein–protein interactions, interactions with other partners including nucleic acids and other ligands, are enriched in post-translational modification sites, and are characterized by specific evolutionary patterns. We mapped the disorder into an integrated network of PCD pathways and into the interactomes of selected proteins that are involved in the p53-mediated apoptotic signaling pathway.     

Temporary secondary structures in tau,an intrinsically disordered protein

The tau protein belongs to the category of intrinsically disordered proteins, which in their native state do not have an average stable structure and fluctuate between many conformations. In its physiological state, tau helps nucleating and stabilising the microtubules in the axons of the neurons. On the other hand, the same tau is involved in the development of Alzheimer disease, when it aggregates in paired helical filaments forming fibrils, which form insoluble tangles. The beginning of the pathological aggregation of tau has been attributed to a local transition of protein portions from random coil to a β-sheet. These structures would very likely be transient; therefore, we performed a molecular dynamics simulation of tau to gather information on the existence of segments of tau endowed with a secondary structure. We combined the results of our simulation with small-angle X-ray scattering experimental data to extract from the dynamics a set of most probable conformations of tau. The analysis of these conformations highlights the presence of transient secondary structures such as turns, β-bridges, β-sheets and α-helices. It also shows that a large segment of the N-terminal region is found near the repeats domain in a globular-like shape.     

Structural divergence is more extensive than sequence divergence for a family of intrinsically disordered proteins

The p53 transactivation domain (p53TAD) is an intrinsically disordered protein (IDP) domain that undergoes coupled folding and binding when interacting with partner proteins like the E3 ligase, MDM2, and the 70 kDa subunit of replication protein A, RPA70. The secondary structure and dynamics of six closely related mammalian homologues of p53TAD were investigated using nuclear magnetic resonance (NMR) spectroscopy. Differences in both transient secondary structure and backbone dynamics were observed for the homologues. Many of these differences were localized to the binding sites for MDM2 and RPA70. The amount of transient helical secondary structure observed for the MDM2 binding site was lower for the dog and mouse homologues, compared with human, and the amount of transient helical secondary structure observed for the RPA70 binding site was higher for guinea pig and rabbit, compared with human. Differences in the amount of transient helical secondary structure observed for the MDM2 binding site were directly related to amino acid substitutions occurring on the solvent exposed side of the amphipathic helix that forms during the p53TAD/MDM2 interaction. Differences in the amount of transient helical secondary structure were not as easily explained for the RPA70 binding site because of its extensive sequence divergence. Clustering analysis shows that the divergence in the transient secondary structure of the p53TAD homologues exceeds the amino acid sequence divergence. In contrast, strong correlations were observed between the backbone dynamics of the homologues and the sequence identity matrix, suggesting that the dynamic behavior of IDPs is a conserved evolutionary feature. Proteins 2013; 81:1686–1698. © 2013 Wiley Periodicals, Inc.     

Exploring the binding diversity of intrinsically disordered proteins involved in one‐to‐many binding

Molecular recognition features (MoRFs) are intrinsically disordered protein regions that bind to partners via disorder‐to‐order transitions. In one‐to‐many binding, a single MoRF binds to two or more different partners individually. MoRF‐based one‐to‐many protein–protein interaction (PPI) examples were collected from the Protein Data Bank, yielding 23 MoRFs bound to 2–9 partners, with all pairs of same‐MoRF partners having less than 25% sequence identity. Of these, 8 MoRFs were bound to 2–9 partners having completely different folds, whereas 15 MoRFs were bound to 2–5 partners having the same folds but with low sequence identities. For both types of partner variation, backbone and side chain torsion angle rotations were used to bring about the conformational changes needed to enable close fits between a single MoRF and distinct partners. Alternative splicing events (ASEs) and posttranslational modifications (PTMs) were also found to contribute to distinct partner binding. Because ASEs and PTMs both commonly occur in disordered regions, and because both ASEs and PTMs are often tissue‐specific, these data suggest that MoRFs, ASEs, and PTMs may collaborate to alter PPI networks in different cell types. These data enlarge the set of carefully studied MoRFs that use inherent flexibility and that also use ASE‐based and/or PTM‐based surface modifications to enable the same disordered segment to selectively associate with two or more partners. The small number of residues involved in MoRFs and in their modifications by ASEs or PTMs may simplify the evolvability of signaling network diversity.     

Molecular principles of the interactions of disordered proteins

Thorough knowledge of the molecular principles of protein-protein recognition is essential to our understanding of protein function at the cellular level. Whereas interactions of ordered proteins have been analyzed in great detail, complexes of intrinsically unstructured/disordered proteins (IUPs) have hardly been addressed so far. Here, we have collected a database of 39 complexes of experimentally verified IUPs, and compared their interfaces with those of 72 complexes of ordered, globular proteins. The characteristic differences found between the two types of complexes suggest that IUPs represent a distinct molecular implementation of the principles of protein-protein recognition. The interfaces do not differ in size, but those of IUPs cover a much larger part of the surface of the protein than for their ordered counterparts. Moreover, IUP interfaces are significantly more hydrophobic relative to their overall amino acid composition, but also in absolute terms. They rely more on hydrophobic-hydrophobic than on polar-polar interactions. Their amino acids in the interface realize more intermolecular contacts, which suggests a better fit with the partner due to induced folding upon binding that results in a better adaptation to the partner. The two modes of interaction also differ in that IUPs usually use only a single continuous segment for partner binding, whereas the binding sites of ordered proteins are more segmented. Probably, all these features contribute to the increased evolutionary conservation of IUP interface residues. These noted molecular differences are also manifested in the interaction energies of IUPs. Our approximation of these by low-resolution force-fields shows that IUPs gain much more stabilization energy from intermolecular contacts, than from folding, i.e. they use their binding energy for folding. Overall, our findings provide a structural rationale to the prior suggestions that many IUPs are specialized for functions realized by protein-protein interactions.     

CDF it all: Consensus prediction of intrinsically disordered proteins based on various cumulative distribution functions

Many biologically active proteins are intrinsically disordered. A reasonable understanding of the disorder status of these proteins may be beneficial for better understanding of their structures and functions. The disorder contents of disordered proteins vary dramatically, with two extremes being fully ordered and fully disordered proteins. Often, it is necessary to perform a binary classification and classify a whole protein as ordered or disordered. Here, an improved error estimation technique was applied to develop the cumulative distribution function (CDF) algorithms for several established disorder predictors. A consensus binary predictor, based on the artificial neural networks, NN-CDF, was developed by using output of the individual CDFs. The consensus method outperforms the individual predictors by 4-5% in the averaged accuracy.     

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.     

Systematic analysis of tropomodulin/tropomyosin interactions uncovers fine-tuned binding specificity of intrinsically disordered proteins

An intriguing regulatory mechanism is the ability of some proteins to recognize their binding partners in an isoform‐specific manner. In this study we undertook a systematic analysis of the specificity of the tropomodulin (Tmod) interaction with tropomyosin (TM) to show that affinities of different Tmod isoforms to TM are isoform‐dependent. Intrinsic disorder predictions, alignment of sequences, and circular dichroism were utilized to establish a structural basis for these isoform‐specific interactions. The affinity of model peptides derived from the N‐terminus of different TM isoforms to protein fragments that correspond to the two TM‐binding sites of different Tmod isoforms were analyzed. Several residues were determined to be responsible for the isoform‐dependent differences in affinity. We suggest that changing a set of residues rather than a single residue is needed to alter the binding affinity of one isoform to mimic the affinity of another isoform. The general intrinsic disorder predictor, PONDR® VLXT, was shown to be a useful tool for analyzing regions involved in isoform‐specific binding and for predicting the residues important for isoform differences in binding. Knowing the residues responsible for isoform‐specific affinity creates a tool suitable for studying the influence of Tmod/TM interactions on sarcomere assembly in muscle cells or actin dynamics in non‐muscle cells. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrogen exchange of disordered proteins in Escherichia coli

A truly disordered protein lacks a stable fold and its backbone amide protons exchange with solvent at rates predicted from studies of unstructured peptides. We have measured the exchange rates of two model disordered proteins, FlgM and α-synuclein, in buffer and in Escherichia coli using the NMR experiment, SOLEXSY. The rates are similar in buffer and cells and are close to the rates predicted from data on small, unstructured peptides. This result indicates that true disorder can persist inside the crowded cellular interior and that weak interactions between proteins and macromolecules in cells do not necessarily affect intrinsic rates of exchange.     

Beta turn propensity and a model polymer scaling exponent identify intrinsically disordered phase-separating proteins

The complex cellular milieu can spontaneously demix, or phase separate, in a process controlled in part by intrinsically disordered (ID) proteins. A protein''s propensity to phase separate is thought to be driven by a preference for protein–protein over protein–solvent interactions. The hydrodynamic size of monomeric proteins, as quantified by the polymer scaling exponent (v), is driven by a similar balance. We hypothesized that mean v, as predicted by protein sequence, would be smaller for proteins with a strong propensity to phase separate. To test this hypothesis, we analyzed protein databases containing subsets of proteins that are folded, disordered, or disordered and known to spontaneously phase separate. We find that the phase-separating disordered proteins, on average, had lower calculated values of v compared with their non-phase-separating counterparts. Moreover, these proteins had a higher sequence-predicted propensity for β-turns. Using a simple, surface area-based model, we propose a physical mechanism for this difference: transient β-turn structures reduce the desolvation penalty of forming a protein-rich phase and increase exposure of atoms involved in π/sp² valence electron interactions. By this mechanism, β-turns could act as energetically favored nucleation points, which may explain the increased propensity for turns in ID regions (IDRs) utilized biologically for phase separation. Phase-separating IDRs, non-phase-separating IDRs, and folded regions could be distinguished by combining v and β-turn propensity. Finally, we propose a new algorithm, ParSe (partition sequence), for predicting phase-separating protein regions, and which is able to accurately identify folded, disordered, and phase-separating protein regions based on the primary sequence.     

Sequence patterns associated with disordered regions in proteins

The relationship between amino acid sequence and intrinsic disorder in proteins is investigated. Two databases, one of disordered proteins and the other of globular proteins, are analyzed and compared in order to extract simple sequence patterns of a few amino acids or amino acid properties that characterize disordered segments. It is found that a number of reliable, nonrandom associations exists. In particular, two types of patterns appear to be recurrent: a proline-rich pattern and a (positively or negatively) charged pattern. These results indicate that local sequence information can determine disordered regions in proteins. The derived patterns provide some insights into the physical reasons for disordered structures. They should also be helpful in improving currently available prediction methods.     

Transient tertiary structures in tau,an intrinsically disordered protein

An intrinsically disordered protein (IDP) does not have a definite 3D structure, and because of its highly flexible nature it evolves dynamically in very large and diverse regions of the phase space. A standard molecular dynamics run can sample only a limited region of the latter; even though this kind of simulation may be effective in sampling local temporary secondary structures, it is not sufficient to highlight properties that require a larger sampling of the phase space to be detected, like transient tertiary structures. But if the structure of an IDP is dynamically evolved using metadynamics (an algorithm that keeps track of the regions of the phase space already sampled), the system can be forced to wander in a much larger region of the phase space. We have applied this procedure to the simulation of tau, one of the largest totally disordered proteins. Combining the results of the simulation with small-angle X-ray scattering yields a significant improvement in the sampling of the phase space in comparison with standard molecular dynamics, and provides evidence of extended hairpin- and paperclip-like transient tertiary structures of the molecule. The more persistent tertiary pattern is a hairpin folding encompassing part of the N-terminal, the proline-rich domain, the former repeat and a functionally relevant part of the second repeat.