Anchoring intrinsically disordered proteins to multiple targets: lessons from N-terminus of the p53 protein

Structural investigation of CDCA3-Cdh1 protein-protein interactions using in vitro studies and molecular dynamics simulation

The anaphase-promoting complex/cyclosome (APC/C) ubiquitin ligase and its cofactor, Cdh1, regulate the expression of several cell-cycle proteins and their functions during mitosis. Levels of the protein cell division cycle-associated protein 3 (CDCA3), which is functionally required for mitotic entry, are regulated by APC/CCdh1 . CDCA3 is an intrinsically disordered protein and contains both C-terminal KEN box and D-box recognition motifs, enabling binding to Cdh1. Our previous findings demonstrate that CDCA3 has a phosphorylation-dependent non-canonical ABBA-like motif within the linker region bridging these two recognition motifs and is required for efficient binding to Cdh1. Here, we sought to identify and further characterize additional residues that participate within this ABBA-like motif using detailed in vitro experiments and in silico modeling studies. We identified the role of H-bonds, hydrophobic and ionic interactions across the CDCA3 ABBA-like motif in the linker region between KEN and D-box motifs. This linker region adopts a well-defined structure when bound to Cdh1 in the presence of phosphorylation. Upon alanine mutation, the structure of this region is lost, leading to higher flexibility, and alteration in affinities due to binding to alternate sites on Cdh1. Our findings identify roles for the anchoring residues in the non-canonical ABBA-like motif to promote binding to the APC/CCdh1 and regulation of CDCA3 protein levels.

Anchor Residues Govern Binding and Folding of an Intrinsically Disordered Domain

Minimal protein mimics have yielded novel classes of protein-protein interaction inhibitors; however, this success has not been extended to targeting intrinsically disordered proteins, which represent a significant proportion of important therapeutic targets. We sought to determine the requirements for binding an intrinsically disordered region (IDR) by its native binding partner as a prelude to developing minimal protein mimics that regulate IDR interactions. Our analysis reinforces the hypothesis that IDRs reside on a fulcrum between unfolded and folded states and that a handful of key binding residues on partner protein surfaces dictate their folding. Our studies also suggest that minimal mimics of protein surfaces may not offer specific ligands for IDRs and that it would be more judicious to target the globular protein partners of IDRs.

Understanding disorder-to-order transitions in protein-RNA complexes using molecular dynamics simulations

Intrinsically disordered regions (IDRs) in proteins are characterized by their flexibilities and low complexity regions, which lack unique 3 D structures in solution. IDRs play a significant role in signaling, regulation, and binding multiple partners, including DNA, RNA, and proteins. Although various experiments have shown the role of disordered regions in binding with RNA, a detailed computational analysis is required to understand their binding and recognition mechanism. In this work, we performed molecular dynamics simulations of 10 protein-RNA complexes to understand the binding governed by intrinsically disordered regions. The simulation results show that most of the disordered regions are important for RNA-binding and have a transition from disordered-to-ordered conformation upon binding, which often contribute significantly towards the binding affinity. Interestingly, most of the disordered residues are present at the interface or located as a linker between two regions having similar movements. The DOT regions are overlaped or flanked with experimentally reported functionally important residues in the recognition of protein-RNA complexes. This study provides additional insights for understanding the role and recognition mechanism of disordered regions in protein-RNA complexes.Communicated by Ramaswamy H. Sarma.

Modulation of Disordered Proteins with a Focus on Neurodegenerative Diseases and Other Pathologies

Intrinsically disordered proteins (IDPs) do not have rigid 3D structures, showing changes in their folding depending on the environment or ligands. Intrinsically disordered proteins are widely spread in eukaryotic genomes, and these proteins participate in many cell regulatory metabolism processes. Some IDPs, when aberrantly folded, can be the cause of some diseases such as Alzheimer′s, Parkinson′s, and prionic, among others. In these diseases, there are modifications in parts of the protein or in its entirety. A common conformational variation of these IDPs is misfolding and aggregation, forming, for instance, neurotoxic amyloid plaques. In this review, we discuss some IDPs that are involved in neurodegenerative diseases (such as beta amyloid, alpha synuclein, tau, and the “IDP-like” PrP), cancer (p53, c-Myc), and diabetes (amylin), focusing on the structural changes of these IDPs that are linked to such pathologies. We also present the IDP modulation mechanisms that can be explored in new strategies for drug design. Lastly, we show some candidate drugs that can be used in the future for the treatment of diseases caused by misfolded IDPs, considering that cancer therapy has more advanced research in comparison to other diseases, while also discussing recent and future developments in this area of research. Therefore, we aim to provide support to the study of IDPs and their modulation mechanisms as promising approaches to combat such severe diseases.

Folding and binding energy of a calmodulin-binding cell antiproliferative peptide

We carry out a computational study of a calmodulin-binding peptide shown to be effective in reducing cell proliferation. We find several folded states for two short variants of different length of the peptide and determine the location of the binding site on calmodulin, the binding free energy for the different conformers and structural details that play a role in optimal binding. Binding to a hydrophobic pocket in calmodulin occurs via an anchoring phenylalanine residue of the natively disordered peptide, and is enhanced when a neighbouring hydrophobic residue acts as a co-anchor. The shorter sequence possesses better binding to calmodulin, which is encouraging in terms of the development of non-peptide analogues as therapeutic agents.

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.

Advantages of proteins being disordered

The past decade has witnessed great advances in our understanding of protein structure-function relationships in terms of the ubiquitous existence of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs). The structural disorder of IDPs/IDRs enables them to play essential functions that are complementary to those of ordered proteins. In addition, IDPs/IDRs are persistent in evolution. Therefore, they are expected to possess some advantages over ordered proteins. In this review, we summarize and survey nine possible advantages of IDPs/IDRs: economizing genome/protein resources, overcoming steric restrictions in binding, achieving high specificity with low affinity, increasing binding rate, facilitating posttranslational modifications, enabling flexible linkers, preventing aggregation, providing resistance to non-native conditions, and allowing compatibility with more available sequences. Some potential advantages of IDPs/IDRs are not well understood and require both experimental and theoretical approaches to decipher. The connection with protein design is also briefly discussed.

Ligand clouds around protein clouds: a scenario of ligand binding with intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) were found to be widely associated with human diseases and may serve as potential drug design targets. However, drug design targeting IDPs is still in the very early stages. Progress in drug design is usually achieved using experimental screening; however, the structural disorder of IDPs makes it difficult to characterize their interaction with ligands using experiments alone. To better understand the structure of IDPs and their interactions with small molecule ligands, we performed extensive simulations on the c-Myc₃₇₀₋₄₀₉ peptide and its binding to a reported small molecule inhibitor, ligand 10074-A4. We found that the conformational space of the apo c-Myc₃₇₀₋₄₀₉ peptide was rather dispersed and that the conformations of the peptide were stabilized mainly by charge interactions and hydrogen bonds. Under the binding of the ligand, c-Myc₃₇₀₋₄₀₉ remained disordered. The ligand was found to bind to c-Myc₃₇₀₋₄₀₉ at different sites along the chain and behaved like a 'ligand cloud'. In contrast to ligand binding to more rigid target proteins that usually results in a dominant bound structure, ligand binding to IDPs may better be described as ligand clouds around protein clouds. Nevertheless, the binding of the ligand and a non-ligand to the c-Myc₃₇₀₋₄₀₉ target could be clearly distinguished. The present study provides insights that will help improve rational drug design that targets IDPs.

Molecular dynamics of the full-length p53 monomer

The p53 protein is frequently mutated in a very large proportion of human tumors, where it seems to acquire gain-of-function activity that facilitates tumor onset and progression. A possible mechanism is the ability of mutant p53 proteins to physically interact with other proteins, including members of the same family, namely p63 and p73, inactivating their function. Assuming that this interaction might occurs at the level of the monomer, to investigate the molecular basis for this interaction, here, we sample the structural flexibility of the wild-type p53 monomeric protein. The results show a strong stability up to 850 ns in the DNA binding domain, with major flexibility in the N-terminal transactivations domains (TAD1 and TAD2) as well as in the C-terminal region (tetramerization domain). Several stable hydrogen bonds have been detected between N-terminal or C-terminal and DNA binding domain, and also between N-terminal and C-terminal. Essential dynamics analysis highlights strongly correlated movements involving TAD1 and the proline-rich region in the N-terminal domain, the tetramerization region in the C-terminal domain; Lys120 in the DNA binding region. The herein presented model is a starting point for further investigation of the whole protein tetramer as well as of its mutants.

Inherent relationships among different biophysical prediction methods for intrinsically disordered proteins

Intrinsically disordered proteins do not have stable secondary and/or tertiary structures but still function. More than 50 prediction methods have been developed and inherent relationships may be expected to exist among them. To investigate this, we conducted molecular simulations and algorithmic analyses on a minimal coarse-grained polypeptide model and discovered a common basis for the charge-hydropathy plot and packing-density algorithms that was verified by correlation analysis. The correlation analysis approach was applied to realistic datasets, which revealed correlations among some physical-chemical properties (charge-hydropathy plot, packing density, pairwise energy). The correlations indicated that these biophysical methods find a projected direction to discriminate ordered and disordered proteins. The optimized projection was determined and the ultimate accuracy limit of the existing algorithms is discussed.

The most important thing is the tail: multitudinous functionalities of intrinsically disordered protein termini

Many functional proteins do not have well-folded structures in their substantial parts, representing hybrids that possess both ordered and disordered regions. Disorder is unevenly distributed within these hybrid proteins and is typically more common at protein termini. Disordered tails are engaged in a wide range of functions, some of which are unique for termini and cannot be found in other disordered parts of a protein. This review covers some of the key functions of disordered protein termini and emphasizes that these tails are not simple flexible protrusions but are evolved to serve.

Binding of two intrinsically disordered peptides to a multi-specific protein: a combined Monte Carlo and molecular dynamics study

The unique ability of intrinsically disordered proteins (IDPs) to fold upon binding to partner molecules makes them functionally well-suited for cellular communication networks. For example, the folding-binding of different IDP sequences onto the same surface of an ordered protein provides a mechanism for signaling in a many-to-one manner. Here, we study the molecular details of this signaling mechanism by applying both Molecular Dynamics and Monte Carlo methods to S100B, a calcium-modulated homodimeric protein, and two of its IDP targets, p53 and TRTK-12. Despite adopting somewhat different conformations in complex with S100B and showing no apparent sequence similarity, the two IDP targets associate in virtually the same manner. As free chains, both target sequences remain flexible and sample their respective bound, natively -helical states to a small extent. Association occurs through an intermediate state in the periphery of the S100B binding pocket, stabilized by nonnative interactions which are either hydrophobic or electrostatic in nature. Our results highlight the importance of overall physical properties of IDP segments, such as net charge or presence of strongly hydrophobic amino acids, for molecular recognition via coupled folding-binding.

A substantial fraction of our proteins are believed to be partly or completely disordered, meaning that they contain regions that lack a stable folded structure under typical physiological conditions. This is a feature which plays a key role in their functions. For example, it allows them to have many structurally different binding partners which in turn permits the construction of the intricate signaling and regulatory networks necessary to sustain complex biological organisms such as ourselves. Whereas measuring the binding strengths of associations involving disordered proteins is routine, the binding process itself is today still not fully understood. We use two different computational models to study the interactions of a folded protein, S100B, which can bind various disordered peptides. In particular, we compare two peptides whose structures are known when in complex with S100B. Our results suggest that, although the peptides assume different structures in the bound state, there are similarities in how they associate with S100B. The possibility to computationally model the interplay between proteins is an important complement to experiments, by identifying crucial steps in the binding process. This is essential to understand, e.g., how single mutations sometimes lead to serious diseases.

A preformed binding interface in the unbound ensemble of an intrinsically disordered protein: evidence from molecular simulations

Intrinsically disordered proteins play an important role in cellular signalling, mediated by their interactions with other biomolecules. A key question concerns the nature of their binding mechanism, and whether the bound structure is induced only by proximity to the binding partner. This is difficult to answer through experiment alone because of the very heterogeneous nature of the unbound ensemble, and the probable rapid interconversion of the various unbound structures. Here we report the most extensive set of simulations on NCBD to date: we use large-scale replica exchange molecular dynamics to explore the unbound state. An important feature of the study is the use of an atomistic force field that has been parametrised against experimental data for weakly structured peptides, together with an accurate explicit water model. Neither the force field nor the starting conformations are biased towards a particular structure. The regions of NCBD that have high helical propensity in the simulations correspond closely to helices in the 'core' unbound conformation determined by NMR, although no single member of the simulated unbound ensemble closely resembles the core conformation, or either of the two known bound conformations. We have validated the results against NMR spectroscopy and SAXS measurements, obtaining reasonable agreement. The two helices which most stabilise the binding of NCBD with ACTR are formed readily; the third helix, which is less important for binding but is involved in most of the intraprotein contacts of NCBD in the bound conformation, is formed more rarely, and tends not to coexist with the other helices. These results support a mechanism by which NCBD gains the advantages of disorder, while forming binding-competent structures in the unbound state. We obtain support for this mechanism from coarse-grained simulations of NCBD with, and without, its binding partner.

Binding free energy landscape of domain-peptide interactions

Peptide recognition domains (PRDs) are ubiquitous protein domains which mediate large numbers of protein interactions in the cell. How these PRDs are able to recognize peptide sequences in a rapid and specific manner is incompletely understood. We explore the peptide binding process of PDZ domains, a large PRD family, from an equilibrium perspective using an all-atom Monte Carlo (MC) approach. Our focus is two different PDZ domains representing two major PDZ classes, I and II. For both domains, a binding free energy surface with a strong bias toward the native bound state is found. Moreover, both domains exhibit a binding process in which the peptides are mostly either bound at the PDZ binding pocket or else interact little with the domain surface. Consistent with this, various binding observables show a temperature dependence well described by a simple two-state model. We also find important differences in the details between the two domains. While both domains exhibit well-defined binding free energy barriers, the class I barrier is significantly weaker than the one for class II. To probe this issue further, we apply our method to a PDZ domain with dual specificity for class I and II peptides, and find an analogous difference in their binding free energy barriers. Lastly, we perform a large number of fixed-temperature MC kinetics trajectories under binding conditions. These trajectories reveal significantly slower binding dynamics for the class II domain relative to class I. Our combined results are consistent with a binding mechanism in which the peptide C terminal residue binds in an initial, rate-limiting step.

The complex biological processes occurring in living organisms are enabled by numerous networks of interacting proteins. It is therefore of great interest to understand the physical interplay between proteins and, in particular, how this process gives rise to highly specific network connectivities. For a long time, the dominant molecular view of protein-protein interactions was the docking of more or less static folded structures, with specificity obtained from a complementarity in shape and charge distributions. Lately it has been realized that many of the links in protein networks are mediated by interactions between folded domains, on the one hand, and disordered polypeptide segments, on the other. We use an all-atom Monte Carlo based approach which attempts to capture this domain-peptide binding process in full and apply it to representative members of a common domain family. This allows us to examine and compare detailed aspects of the binding free energy landscapes which underlie specificity and affinity. Being able to model domain-peptide binding in a physically sound, yet computationally tractable way is essential for identifying molecular binding mechanisms and opens up possibilities for modifying interaction networks in a controlled way.