Next Article in Journal
Structure–Activity Relationship Study Reveals the Molecular Basis for Specific Sensing of Hydrophobic Amino Acids by the Campylobacter jejuni Chemoreceptor Tlp3
Next Article in Special Issue
Sequence-Dependent Correlated Segments in the Intrinsically Disordered Region of ChiZ
Previous Article in Journal
Blueberry-Derived Exosome-Like Nanoparticles Counter the Response to TNF-α-Induced Change on Gene Expression in EA.hy926 Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Targeting Intrinsically Disordered Proteins through Dynamic Interactions

1
Department of Hematology, Taizhou Central Hospital (Taizhou University Hospital), Taizhou 318000, Zhejiang, China
2
Department of Chemistry, University of Massachusetts, Amherst, MA 01003, USA
3
Department of Biochemistry and Molecular Biology, University of Massachusetts, Amherst, MA 01003, USA
*
Author to whom correspondence should be addressed.
Biomolecules 2020, 10(5), 743; https://doi.org/10.3390/biom10050743
Submission received: 13 April 2020 / Revised: 4 May 2020 / Accepted: 9 May 2020 / Published: 11 May 2020
(This article belongs to the Special Issue Computational Perspectives on Intrinsic Disorder-Based Functionality)

Abstract

:
Intrinsically disordered proteins (IDPs) are over-represented in major disease pathways and have attracted significant interest in understanding if and how they may be targeted using small molecules for therapeutic purposes. While most existing studies have focused on extending the traditional structure-centric drug design strategies and emphasized exploring pre-existing structure features of IDPs for specific binding, several examples have also emerged to suggest that small molecules could achieve specificity in binding IDPs and affect their function through dynamic and transient interactions. These dynamic interactions can modulate the disordered conformational ensemble and often lead to modest compaction to shield functionally important interaction sites. Much work remains to be done on further elucidation of the molecular basis of the dynamic small molecule–IDP interaction and determining how it can be exploited for targeting IDPs in practice. These efforts will rely critically on an integrated experimental and computational framework for disordered protein ensemble characterization. In particular, exciting advances have been made in recent years in enhanced sampling techniques, Graphic Processing Unit (GPU)-computing, and protein force field optimization, which have now allowed rigorous physics-based atomistic simulations to generate reliable structure ensembles for nontrivial IDPs of modest sizes. Such de novo atomistic simulations will play crucial roles in exploring the exciting opportunity of targeting IDPs through dynamic interactions.

1. Introduction

Proteins are central components of regulatory networks that dictate virtually all aspects of cellular decision-making [1]. Demand for more sophisticated signaling in complex multicellular organisms has been met with increasing utilization of proteins that are highly flexible [2,3,4]. In particular, so-called intrinsically disordered proteins (IDPs) account for ~50% of signaling-associated proteins in eukaryotes [5]. These proteins have lower sequence complexity compared to folded proteins, lacking large hydrophobic residues and enriched with charged and polar ones [6]. They do not have stable tertiary structures in the unbound state under physiological conditions, even though they frequently undergo folding transitions upon binding to specific targets [7]. The inherent thermodynamic instability of the structural features of this class of proteins allows their conformational properties to respond sensitively to numerous stimuli, including the binding of various small and large molecules, changes in cellular environments (e.g., pH), and post-translational modifications [8,9,10,11,12,13]. Multiple signals could also be naturally integrated through cooperative responses of the dynamic structure ensemble (such as coupled binding and folding) [14]. These properties make IDPs uniquely suitable for fulfilling the complex signaling need of higher organisms. At the same time, deregulation of IDPs has been associated with many human diseases, including cancers, neurodegenerative diseases, heart disease, and diabetes [5,15,16,17,18,19,20]. For example, over two-thirds of cancer-associated proteins have been predicted to contain extensive regions of intrinsic disorder [5], and predicted disordered regions have been estimated to house almost one quarter of disease-associated missense mutations [21]. There is thus tremendous interest in determining if and how IDPs may be targeted for therapeutic purposes.
The dynamic and heterogeneous nature of unbound IDPs presents substantial challenges for characterization and this has proven to be a major bottleneck for establishing a reliable sequence–structure–function–disease relationship of IDPs [14,22,23,24,25,26]. The lack of a clear understanding of the molecular basis of IDP function and deregulation in diseases has created significant ambiguity on the druggability of most IDPs, including transcription factors [16]. Most existing case studies of targeting IDPs have focused on extending the traditional structure-based screening and drug design strategies and emphasize exploiting residual structures and pre-existing potential binding pockets of the unbound state [27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43]. Nonetheless, it is clear that the disordered nature of IDPs would require novel strategies for targeting as well as new conceptual frameworks for thinking about how small molecule binding could modulate IDP structure and function. In particular, it has been recognized that it may be more useful to consider the problem of targeting IDPs in the context of structural ensemble modulation [44], even though it is generally believed that one still needs to achieve specific interactions, such as by exploiting pre-existing structural features [45]. Many outstanding reviews have already been dedicated towards existing examples along these lines and they also provide extensive discussion of the successes, opportunities, and challenges of targeting IDPs via specific interactions of small molecules in neurodegenerative diseases, cancers, and other diseases [18,45,46,47,48,49,50,51,52,53,54,55].
In this review, we will first summarize important recent advances in physics-based de novo simulations of disordered protein ensembles, including Graphic Processing Unit (GPU) computing, enhanced sampling, and re-balanced protein force fields, and then focus on emerging examples that suggest the exciting possibility of targeting IDPs by directly modulating the disordered ensembles through dynamic and transient interactions. We will discuss the promise of such a broader view of how IDPs may be targeted as well as key challenges and required methodological developments to support targeting IDPs via dynamic interactions.

2. Characterization of Disordered Protein Ensembles: A Crucial Role for Atomistic Simulations

A principal challenge in understanding the druggability and best targeting strategy of IDPs resides in the difficulty of detailed characterization of disordered protein states [14,23,24,56]. These states need to be represented using heterogeneous structure ensembles and are not amenable to traditional high-resolution structure determination methods. Here, we briefly discuss the current status and challenges of disordered protein ensemble determination, which has a direct impact on the ability to devise effective strategies of designing IDP binders and optimizing leads identified from traditional screening efforts.

2.1. Fundamental Challenges in Experimental Determination of Disordered Protein Ensembles

Experimentally, a wide range of biophysical methods can be applied to characterize disordered protein states, including NMR, circular dichroism (CD), small-angle X-ray scattering (SAXS), Förster resonance energy transfer (FRET), hydrogen/deuterium (H/D) exchange, mass spectrometry, and others [57,58]. These methods can provide complementary information on the local, intermediate, and long-range structural organizations of IDPs. NMR in particular is arguably the most powerful technique for structural studies of IDP. Many NMR observables can be measured at residue and atomic levels to infer secondary and tertiary structural properties. SAXS and FRET are highly complementary to NMR and provide information on the long-range global organization of the disordered ensemble. Yet, a fundamental limitation is that these experimental measurements generally reflects the average properties, which alone are not sufficient to uniquely define the underlying heterogeneous ensemble due to the severely underdetermined nature of the structural calculation problem [22,23,25,59,60,61]. At present, the most robust methods generally involve first generating a large number of candidate random structures and then using experimental structural restraints to select and construct optimal sub-ensembles according to various statistical criteria [62,63,64,65,66,67,68,69,70,71]. Nonetheless, these methods rely critically on the ability to generate initial candidate structures that are not only diverse enough to cover the range of accessible states of the protein but also specific enough to contain any nontrivial local (and long-range) structure features associated with a particular protein state. These two requirements are difficult to satisfy simultaneously for most proteins of moderate sizes (e.g., 50–100 residues or longer) and complexity. As a consequence, disordered ensembles constructed using these approaches depend critically on the underlying protein model and/or coil library. Furthermore, these ensembles are generally not proper thermodynamic ensembles. They should be considered just as ensemble models and cannot be used to reliably quantify statistical properties and extract thermodynamic parameters. As such, these experimental restraint-based ensemble construction approaches are likely inadequate for capturing potentially subtle effects of ligand binding on the disordered ensemble.

2.2. Recent Advances in de Novo Simulations of Disordered Protein Ensembles

Given the fundamental challenges of disordered ensemble modeling based on experimental restraints alone, physics-based atomistic simulations have a crucial role to play in helping elucidate the conformational properties of IDPs and establishing a reliable molecular basis of their function and regulation [72,73,74,75,76]. A particularly attractive approach is to first generate the atomistic ensemble using a transferable physics-based force field in absence of any experimental restraints and then use the experimental data for independent evaluation of the quality of the simulated ensemble. Such a de novo simulation approach can effectively overcome the under-determined nature of disordered ensemble calculation, by leveraging the laws of physics that govern the nature of conformational fluctuation of the protein. Successful simulations of disordered protein ensembles have proven to be very challenging, requiring both accurate description of the conformational dependence of energy and sufficient sampling of relevant conformational space of the protein. Early IDP simulations suffered from both systematic biases in the general-purpose protein force fields and severe limitations in conformational sampling [14]. Nonetheless, these physics-based simulation approaches promise to provide rigorous thermodynamic ensembles required for reliable description of IDP–ligand interactions. The accuracy and capability of de novo IDP simulations can be expected to improve continuously over time, benefiting from robust advances in molecular simulation methodologies and high-performance computing hardware. Indeed, important breakthroughs have been made in the last few years in force field accuracy and sampling capability, which arguably have now allowed reliable de novo simulations of at least moderate-sized IDPs in general.

2.2.1. Overcoming Sampling Bottleneck using Enhanced Sampling and GPU Computing

One of the most important recent advances in molecular dynamics (MD) simulations is the widespread availability of efficient GPU-enabled algorithms in virtually all major molecular simulation packages [77,78,79,80,81,82]. Modern GPUs can process thousands of threads in parallel to accelerate explicit solvent MD simulations by up to 100× compared to traditional CPU computing, significantly boosting the sampling capability. For example, the most efficient GPU-enabled MD codes can yield 100–200 ns per day for systems of ~1,000,000 atoms on a single NVIDIA RTX 2080Ti GPU card that costs only ~$1000. The ability to efficiently sample the protein conformational space has further benefitted from the emergence of enhanced sampling techniques [83,84,85,86,87,88,89,90], particularly various replica exchange (REX)-based methods. Among the various REX methods, replica exchange with solute tempering methods (REST) [86,91] is particularly suitable for atomistic simulations of IDPs in explicit solvent. In REST, a selected region of the system (e.g., the protein solute or a flexible segment of the protein) is subjected to tempering (i.e., random walk in the temperature space) while the rest of the system is maintained at a constant temperature (e.g., room temperature). This is achieved by scaling the interactions within the selected region and between it and the rest of the system. Because the number of required replicas of replica exchange simulations scales as the square root of the number of atoms, REST can significantly reduce the number of replicas required for covering the needed temperature space (by ~3-fold), thus overcoming a critical drawback in traditional temperature REX simulation in explicit solvent. The role of enhanced sampling in the recent successes of simulating dynamic IDP–ligand interactions will become evident in the later sections of this review.
The importance of enhanced sampling for the generation of converged ensembles cannot be over emphasized. For example, Figure 1A shows the evolution of the per-residue β-sheet structure during a 30-μs conventional MD simulation of an intrinsically disordered Aβ40 peptide in explicit solvent at 300 K, performed with Anton specialized hardware using the well optimized a99SB-disp force field [92]. Note that this trajectory is one to two orders of magnitude longer than typical MD simulations performed using general purpose CPU- or GPU-based high-performance computing platforms. Yet, very few reversible transitions are observed. For example, a transient β-hairpin spanning residues 15 to 36 persists for about 5 μs from ~10 to 15 μs and never appears again for the rest of the simulation. As a result, the final average residue helix and β-sheet probability profiles calculated from the first and second halves of the trajectory differ greatly, reflecting a very limited level of convergence. There is thus danger in relying on standard MD simulations in deriving quantitative characterizations of disordered protein ensembles. It should also be emphasized that achieving a sufficient level of convergence required for resolving potentially subtle changes in the disordered ensemble, such as upon ligand binding, can be extremely challenging even with enhanced sampling. It is critical to carefully analyze and establish the level of convergence for proper interpretation of simulated ensembles. Ideally, one should perform two or more independent simulations using distinct initial conformations and compare the resulting ensembles. The simulated ensemble from a single continuous run may appear to stop changing with respect to simulation time due to trapping in numerous local energy minima, giving rise to a misleading impression of convergence.

2.2.2. Balanced Explicit Protein Force Fields for Describing Disordered Protein Ensembles

The dramatically improved sampling capability has facilitated extensive efforts to reparametrize general-purpose protein force fields to achieve greater balance of describing protein conformational equilibria. Studies of disordered protein states have been a key driver of these developments and several well-characterized IDPs have been widely used as training systems and/or benchmarks for force field optimization [92,93,94,95,96,97,98,99]. Many force field variants have been developed in recent years, including Amber ff98SB [100], ff99SB*-ILDN [101] and variants [95,102], ff03ws [103], ff14SB [104], ff99SBnmr [105], CHARMM22* [106], CHARMM36m (and C36mw) [93], a99SB-disp [92], and others. A key focus of these optimization efforts has been to rebalance the protein–protein, protein–water, and water–water interactions. Earlier versions of general-purpose protein force fields consistently over-stabilize nonspecific protein–protein interactions and lead to overly compact conformational ensembles for disordered protein states [107,108]. It was demonstrated that such bias could be effectively compensated by directly increasing the strengths of protein–water dispersion interactions [95,103,109], even though other components of the force field should also be reparametrized for self-consistency. The latest CHARMM36m and a99SB-disp force fields, in particular, have been systematically reparametrized based on extensive simulations of tens of globular and disordered proteins and achieve impressive levels of accuracy for describing both structured and unstructured proteins.
In a recent benchmark study, six of the latest protein force fields were evaluated using the 61-residue N-terminal transactivation domain (TAD) of tumor suppressor p53, which is a very challenging system due to its size and complex conformational features [99]. It has been extensively characterized by NMR, SAXS, and single-molecule FRET and shown to contain a range of nontrivial local and long-range residual structures [110]. The disordered ensemble of p53-TAD was calculated using REST2-enhanced sampling using GPU-accelerated GROMACS 5.1.4 [80,111] patched with PLUMED 2.3.0 [112,113,114]. Each REST2 simulation lasted 1.0 μs per replica, representing one of the most extensive atomistic simulations of IDPs of similar sizes. The results show that the ensembles generated using the force field a99SB-disp yield the best agreement with the experimental data at both secondary and tertiary structure levels. For example, the back-calculated NMR paramagnetic relaxation enhancement (PRE) effects are highly consistent with the experimental results for all four available labelling sites (Figure 2). This suggests that the simulated ensembles not only have the proper overall chain dimension but also recapitulate much of the transient long-range ordering within the unbound state of p53-TAD. The latter is an extremely challenging task. The fact that this could be achieved by a99SB-disp represents an exciting breakthrough, suggesting that de novo atomistic simulations are now ready to provide a reliable approach for detailed characterization of the disordered ensembles of at least moderately sized IDPs.

2.2.3. Implicit Solvent as a Promising Alternative for Simulating Disordered Ensembles

We note that implicit solvent protein force fields have also been developed and deployed for atomistic simulations of IDPs with various levels of success [75,96,98,115]. Implicit treatment of solvent reduces the simulation system size ~10-fold by direct estimation of the solvation free energy. It could provide important advantages for satisfying the simultaneous requirements of adequate sampling and sufficient force field accuracy for simulating disordered protein states. The ABSINTH model, in particular, has demonstrated significant successes in mapping the sequence–conformational space relationship of IDPs [25]. An improved version named ABSINTH-C was recently developed by including the backbone torsion cross-terms optimized based on experimentally derived statistics [98]. Independently, the generalized Born with the molecular volume 2 (GMMV2) model, which is considered one of the most accurate implicit solvent models, was recently implemented on GPU within the CHARMM/OpenMM interface [116]. GPU-accelerated GBMV2 is about 60-fold faster and provides a competitive alternative to explicit solvent simulations for studying the IDP structure and interaction. This model was previously optimized based on enhanced sampling of model peptides and shown to be capable of accurately describing the conformational properties of both folded and unfolded peptides [96]. The development of the GPU-accelerated version thus removed a key bottleneck to broader application of GBMV2 for atomistic simulations of the IDP structure and interactions. Nonetheless, whether these implicit solvent models could provide a viable alternative to explicit solvent simulations in studies of IDP–ligand interactions is yet to be demonstrated.

3. Modulating Disordered Protein Ensembles via Dynamic Interactions

To date, small molecular targeting of IDPs has mostly focused on proteins involved in neurodegenerative diseases, such as amyloid β (Aβ) peptides, α-synuclein, and tau protein, and disordered regions of cancer-associated transcription factors, such as p53, c-Myc, EWS-FLI1, KLF5, and others [16,117]. Advances in experimental and computational methods for studying disordered protein ensembles have allowed examination in greater details the interactions between IDPs and small molecules. There has been an emergence of examples showing that small molecules could modulate the disordered ensemble itself through nonspecific and dynamic interactions and achieve specific functional effects, in complete contrast to the traditional paradigm of drug binding that emphasizes strong specific interactions. Such a dynamic mode of IDP–small molecule interactions is reminiscent of “fuzzy complexes” in protein–protein interactions involving IDPs [118,119,120]. The observation that small molecules could induce substantial effects through dynamic interactions is fascinating and suggests a broader and more effective strategy for targeting IDPs in general.

3.1. Dynamic Interactions of c-Myc Inhibitors

One of the earliest examples of dynamic interactions between small molecules and IDPs concerns c-Myc inhibitors [27,28]. NMR, CD, and fluorescence studies initially suggested that these inhibitors bound specifically to multiple independent sites in the monomeric and disordered c-Myc, which induced Max-binding-incompatible conformations to disrupt the c-Myc–Max interaction [27]. Subsequent explicit solvent simulations by Michel and Cuchillo in 2012 [121] revealed that the interaction between c-Myc and one of the inhibitors, 10058-F4, was actually dynamic and involved many short-lived contacts. Such a dynamic and nonspecific nature of the interaction could explain the observation that small modifications to the ligand had limited effects on the binding affinity [122]. A similar mode of interaction was proposed by Jin et al. in 2013 [123] for another c-Myc inhibitor, 10074-A4, where atomistic simulations in explicit and implicit solvent force fields revealed the ligand to form a “ligand cloud” and interact dynamically with the disordered “protein cloud”. Importantly, the dynamic nature of the interaction between 10058-F4 and c-Myc402-412 can confer sequence specificity, even though the interaction involves many transient contacts instead of well-defined specific ones [124]. It has been further suggested that such dynamic interactions may be driven by entropic expansion of the IDP conformational space [51]. Indeed, thermodynamics analysis showed that binding of 10058-F4 to c-Myc402-412 was dominated by the entropic contribution (−20.7± −4.2 kJ/mol out of the total binding free energy of −27.6 ± −8.5 kJ/mol) [124]. The caveat, however, is that binding of hydrophobic ligands is generally associated with large entropic contributions due to the release of restricted water molecules near the hydrophobic surface. It is thus not clear if c-Myc indeed undergoes entropic expansion upon 10058-F4 binding. An in-depth NMR study of the interaction of a small molecule with intrinsically disordered p27Kip1 found little evidence of ligand-induced conformational space expansion [31]. Instead, ligand binding was shown to mainly shift the populations of pre-existing states.

3.2. Inhibition of Aggregation by Induced Compaction

Dynamic interactions were also found to underlie the mechanism of α-synuclein aggregation inhibition by analogs of cyclized nordihydroguaiaretic acid (cNDGA) [125]. The structural basis of cNDGA inhibition was characterized using an array of biochemical and biophysical methods, including NMR and fluorescence correlation spectroscopy. The results revealed that cNDGA induced modest compaction of the conformational ensemble of monomeric α-synuclein, apparently mediated by dynamic and transient interactions with the protein and without hindering membrane association. cNDGA-treated α-synuclein is resistant to aggregation even when seeded with α-synuclein aggregates. Importantly, cNDGA was further shown to be effective in reducing α-synuclein-driven neurodegeneration in C. elegans. The observation that dynamic interactions between a small molecule and IDP could be functionally effective both in vitro and in vivo is very encouraging and supports the promise of targeting IDPs using dynamic interactions for therapeutics.
Induced compaction of the disordered conformational ensemble has also been predicted to underlie the inhibition mechanism of two drugs under clinical trials for treating Alzheimer’s diseases [126]. The disordered ensembles of Aβ42 with and without tramiprosate (homotaurine; HT) and scyllo-inositol (SI) were calculated using REST2 simulations in the CHARMM36m force field that lasted 10 μs per replica, making them the most extensive atomistic simulations of Aβ42 to date. The resulting ensembles are well converged and appear consistent with the NMR chemical shifts. Comparing the ensembles with and without the ligand showed that both HT and SI mainly reduced the β propensity in the C-terminal region with minimal secondary structure perturbation in the rest of the peptide. Intriguingly, both HT and SI were found to induce modest compaction of the conformational ensemble, particularly in the C-terminal segment that is known to be important for amyloid fibril formation (Figure 3a–c). Detailed analysis further revealed that the effects of both HT and SI binding were achieved via dynamic and nonspecific interactions with various backbone and sidechain moieties of the peptide. It is noteworthy that the conformational modulation effects of both HT and SI can be very difficult to detect at the ensemble level using bulk measurements, highlighting a critical role for reliable atomistic simulations that leverage recent advances in both protein force field quality and sampling capability. Nonetheless, additional validation is needed to support the predicted conformational shifts induced by drugs and to establish the roles of such conformational changes in the mechanisms of drug action.

3.3. Modulating Regulatory IDPs via Dynamic Interactions

De novo atomistic simulations have also been integrated with NMR and biophysical experiments to examine how an anticancer drug, epigallocatechin gallate (EGCG), modulates the disordered unbound state of p53-TAD [127]. EGCG is a major active ingredient of green tea and has been reported to have anticancer effects in both animal studies and clinical trials [128,129,130]. The results suggested that EGCG also interacted dynamically with p53-TAD through numerous transient and nonspecific interactions, which appeared consistent with NMR chemical shift titration results. Multiple hydrophobic and particularly aromatic sidechains contribute significantly to EGCG binding. The dynamic interaction with EGCG was predicted to induce significant conformational compaction of p53-TAD in the N-terminal region (Figure 3d,e), which appears consistent with the SAXS measurements. The compaction could shield the p53-TAD site required for interacting with MDM2 and thus inhibit p53 degradation to promote its anticancer activities. It is noteworthy that EGCG has also been shown to be active in inhibiting the aggregation of multiple proteins, including Aβ peptides and α-synuclein [131,132]. The implication is that dynamic interactions of EGCG could provide a molecular basis for promiscuous selectivity, which is a fascinating property to investigate further using a combination of experiments and simulation. Recently, Neira et al. screened a set of FDA-approved drugs and identified 15 compounds that could bind NUPR1, a multi-functional IDP involved in pancreatic ductal adenocarcinoma [133,134]. NMR chemical shift analysis suggested that NUPR1 remained disordered in complex with all compounds, which is consistent with an inhibition mechanism involving transient and dynamic IDP–ligand interactions. Importantly, these compounds showed efficacy in cell-based assays, the most effective of which was found to completely arrest tumor growth in a mouse model.

4. Concluding Discussions

IDPs have remained an extremely challenging class of proteins to target using small molecules. Albeit limited, successful inhibitors have been discovered and designed for several IDPs involved in cancers and neurodegenerative diseases, suggesting that IDPs are not undruggable. Nonetheless, the unstructured and dynamic nature of IDPs is distinct from typical protein targets with well-defined binding pockets. It requires new conceptional frameworks to guide the development of novel strategies for discovering and designing small molecules that can modulate IDP functions. Traditional structure-based screening and lead optimization strategies are clearly inadequate, even though some success has been demonstrated in deploying existing tools to identify possible binding pockets and target pre-existing structural elements. Such structural elements are lightly populated and generally too small to harbor significant pockets for small molecular binding that relies on specific interactions to achieve high affinity. In fact, there is a great uncertainty on whether high-affinity binding to IDPs is feasible with small molecules (e.g., to meet the typical industrial standard of dissociation constants of nM or lower). It is encouraging that examples are emerging that small molecules could modulate the IDP ensembles entirely through dynamic nonspecific interactions. Importantly, there is evidence that high-affinity binding may not be necessary to induce functional responses in vitro and in vivo. This may reflect a fundamental nature of how IDPs mediate function in biology, in that the disordered ensemble of an IDP is poised to respond sensitively to a wide array of cellular signals to support signal transduction and cellular regulation [10]. Therefore, there is a great potential and promise for targeting IDPs through dynamic interactions with small molecules.
It is noteworthy that IDPs have been found to play central roles in mediating liquid–liquid phase separation (LLPS) that underlies a range of cellular processes [135,136,137]. How small molecules may modulate the equilibrium and properties of these biological condensates is essentially unknown at this point. However, it is conceivable that LLPS may also be ideally targeted using dynamic interactions with small molecules that modulate the conformational flexibility and preference of disordered regions, which in turn modify the multivalent interaction profiles and entropic contribution that affect the condensation process as well as the properties of the condensate itself.
Elucidating the molecular details of dynamic interactions between IDPs and small molecules will require further development and integration of new experimental and computational methodologies. The capability for reliable disordered ensemble characterization with and without small molecules will almost certainly be required for any future design and optimization strategy to discover drugs that target IDPs through dynamic interactions. This unfortunately remains a formidable task. Bulk experimental measurements on average properties alone are not sufficient to uniquely define the disordered ensemble. The dynamic and transient nature of molecular contacts can be very difficult to detect and resolve experimentally [138,139]. Leveraging significant recent advances in the protein force field quality, sampling techniques, and GPU computing, de novo atomistic simulations are now poised to help meet these challenges and play a pivotal role in establishing the molecular basis of dynamic IDP–small molecule interactions.

Author Contributions

data analysis, X.L.; writing—draft and edit, J.C. (Jianhan Chen); writing—review and edit, J.C. (Jianlin Chen) and X.L.; All authors have read and agreed to the published version of the manuscript.

Funding

This work is partially supported by National Institutes of Health (GM114300 to JHC).

Acknowledgments

We thank D. E. Shaw Research for giving us access to the MD trajectories of Aβ40. The authors thank Chungwen Liang for helpful discussions. We would also link to thank anonymous reviewers for their careful reading and insightful suggestions that have greatly improved both the content and writing of this review.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bushweller, J.H. Targeting transcription factors in cancer—From undruggable to reality. Nat. Rev. Cancer 2019, 19, 611–624. [Google Scholar] [CrossRef]
  2. Wright, P.E.; Dyson, H.J. Intrinsically unstructured proteins: Re-assessing the protein structure-function paradigm. J. Mol. Biol. 1999, 293, 321–331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Tompa, P. Intrinsically unstructured proteins. Trends Biochem. Sci. 2002, 27, 527–533. [Google Scholar] [CrossRef]
  4. Dunker, A.K.; Cortese, M.S.; Romero, P.; Iakoucheva, L.M.; Uversky, V.N. Flexible nets—The roles of intrinsic disorder in protein interaction networks. FEBS J. 2005, 272, 5129–5148. [Google Scholar] [CrossRef] [PubMed]
  5. Iakoucheva, L.M.; Brown, C.J.; Lawson, J.D.; Obradovic, Z.; Dunker, A.K. Intrinsic disorder in cell-signaling and cancer-associated proteins. J. Mol. Biol. 2002, 323, 573–584. [Google Scholar] [CrossRef] [Green Version]
  6. Romero, P.; Obradovic, Z.; Li, X.H.; Garner, E.C.; Brown, C.J.; Dunker, A.K. Sequence complexity of disordered protein. Proteins-Struct. Funct. Genet. 2001, 42, 38–48. [Google Scholar] [CrossRef]
  7. Dyson, H.J.; Wright, P.E. Coupling of folding and binding for unstructured proteins. Curr. Opin. Struct. Biol. 2002, 12, 54–60. [Google Scholar] [CrossRef]
  8. Uversky, V.N.; Oldfield, C.J.; Dunker, A.K. Showing your ID: Intrinsic disorder as an ID for recognition, regulation and cell signaling. J. Mol. Recognit. 2005, 18, 343–384. [Google Scholar] [CrossRef]
  9. Dyson, H.J.; Wright, P.E. Intrinsically unstructured proteins and their functions. Nat. Rev. Mol. Cell Biol. 2005, 6, 197–208. [Google Scholar] [CrossRef]
  10. Smock, R.G.; Gierasch, L.M. Sending signals dynamically. Science 2009, 324, 198–203. [Google Scholar] [CrossRef] [Green Version]
  11. Hilser, V.J.; Thompson, E.B. Intrinsic disorder as a mechanism to optimize allosteric coupling in proteins. Proc. Natl. Acad. Sci. USA 2007, 104, 8311–8315. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. White, J.T.; Li, J.; Grasso, E.; Wrabl, J.O.; Hilser, V.J. Ensemble allosteric model: Energetic frustration within the intrinsically disordered glucocorticoid receptor. Philos. Trans. R. Soc. B-Biol. Sci. 2018, 373. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, L.H.; Li, M.D.; Liu, Z.R. A comprehensive ensemble model for comparing the allosteric effect of ordered and disordered proteins. PLoS Comput. Biol. 2018, 14, e1006393. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Chen, J. Towards the physical basis of how intrinsic disorder mediates protein function. Arch. Biochem. Biophys. 2012, 524, 123–131. [Google Scholar] [CrossRef] [PubMed]
  15. Uversky, V.N.; Oldfield, C.J.; Dunker, A.K. Intrinsically disordered proteins in human diseases: Introducing the D-2 concept. Annu. Rev. Biophys. 2008, 37, 215–246. [Google Scholar] [CrossRef] [PubMed]
  16. Tsafou, K.; Tiwari, P.B.; Forman-Kay, J.D.; Metallo, S.J.; Toretsky, J.A. Targeting Intrinsically Disordered Transcription Factors: Changing the Paradigm. J. Mol. Biol. 2018, 430, 2321–2341. [Google Scholar] [CrossRef] [PubMed]
  17. Giri, R.; Kumar, D.; Sharma, N.; Uversky, V.N. Intrinsically Disordered Side of the Zika Virus Proteome. Front. Cell. Infect. Microbiol. 2016, 6, 144. [Google Scholar] [CrossRef]
  18. Hu, G.; Wu, Z.H.; Wang, K.; Uversky, V.N.; Kurgan, L. Untapped Potential of Disordered Proteins in Current Druggable Human Proteome. Curr. Drug Targets 2016, 17, 1198–1205. [Google Scholar] [CrossRef]
  19. Sluchanko, N.N.; Bustos, D.M. Intrinsic disorder associated with 14-3-3 proteins and their partners. Prog. Mol. Biol. Trans. Sci. 2019, 166, 19–61. [Google Scholar]
  20. Uversky, V.N. Intrinsic disorder here, there, and everywhere, and nowhere to escape from it. Cell. Mol. Life Sci. 2017, 74, 3065–3067. [Google Scholar] [CrossRef]
  21. Vacic, V.; Iakoucheva, L.M. Disease mutations in disordered regions--exception to the rule? Mol. BioSyst. 2012, 8, 27–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Mittag, T.; Forman-Kay, J.D. Atomic-level characterization of disordered protein ensembles. Curr. Opin. Struct. Biol. 2007, 17, 3–14. [Google Scholar] [CrossRef] [PubMed]
  23. Fisher, C.K.; Stultz, C.M. Constructing ensembles for intrinsically disordered proteins. Curr. Opin. Struct. Biol. 2011, 21, 426–431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Click, T.H.; Ganguly, D.; Chen, J. Intrinsically Disordered Proteins in a Physics-Based World. Int. J. Mol. Sci. 2010, 11, 5292–5309. [Google Scholar] [CrossRef] [Green Version]
  25. Das, R.K.; Ruff, K.M.; Pappu, R.V. Relating sequence encoded information to form and function of intrinsically disordered proteins. Curr. Opin. Struct. Biol. 2015, 32, 102–112. [Google Scholar] [CrossRef] [PubMed]
  26. Yan, J.; Dunker, A.K.; Uversky, V.N.; Kurgan, L. Molecular recognition features (MoRFs) in three domains of life. Mol. BioSyst. 2016, 12, 697–710. [Google Scholar] [CrossRef] [Green Version]
  27. Hammoudeh, D.I.; Follis, A.V.; Prochownik, E.V.; Metallo, S.J. Multiple independent binding sites for small-molecule inhibitors on the oncoprotein c-Myc. J. Am. Chem. Soc. 2009, 131, 7390–7401. [Google Scholar] [CrossRef]
  28. Mustata, G.; Follis, A.V.; Hammoudeh, D.I.; Metallo, S.J.; Wang, H.; Prochownik, E.V.; Lazo, J.S.; Bahar, I. Discovery of novel Myc-Max heterodimer disruptors with a three-dimensional pharmacophore model. J. Med. Chem. 2009, 52, 1247–1250. [Google Scholar] [CrossRef] [Green Version]
  29. Krishnan, N.; Koveal, D.; Miller, D.H.; Xue, B.; Akshinthala, S.D.; Kragelj, J.; Jensen, M.R.; Gauss, C.M.; Page, R.; Blackledge, M.; et al. Targeting the disordered C terminus of PTP1B with an allosteric inhibitor. Nat. Chem. Biol. 2014, 10, 558–566. [Google Scholar] [CrossRef] [Green Version]
  30. Iconaru, L.I.; Ban, D.; Bharatham, K.; Ramanathan, A.; Zhang, W.X.; Shelat, A.A.; Zuo, J.; Kriwacki, R.W. Discovery of Small Molecules that Inhibit the Disordered Protein, p27(Kip1). Sci Rep-Uk 2015, 5, 15686. [Google Scholar] [CrossRef] [Green Version]
  31. Ban, D.; Iconaru, L.I.; Ramanathan, A.; Zuo, J.; Kriwacki, R.W. A Small Molecule Causes a Population Shift in the Conformational Landscape of an Intrinsically Disordered Protein. J. Am. Chem. Soc. 2017, 139, 13692–13700. [Google Scholar] [CrossRef] [PubMed]
  32. Zhang, Y.G.; Cao, H.Q.; Liu, Z.R. Binding cavities and druggability of intrinsically disordered proteins. Protein Sci. 2015, 24, 688–705. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Kim, H.; Han, K.H. PreSMo Target-Binding Signatures in Intrinsically Disordered Proteins. Mol. Cells 2018, 41, 889–899. [Google Scholar] [PubMed]
  34. Zhu, M.; De Simone, A.; Schenk, D.; Toth, G.; Dobson, C.M.; Vendruscolo, M. Identification of small-molecule binding pockets in the soluble monomeric form of the A beta 42 peptide. J. Chem. Phys. 2013, 139, 07B609_1. [Google Scholar] [CrossRef] [Green Version]
  35. Pickhardt, M.; Neumann, T.; Schwizer, D.; Callaway, K.; Vendruscolo, M.; Schenk, D.; St George-Hyslop, P.; Mandelkow, E.M.; Dobson, C.M.; McConlogue, L.; et al. Identification of Small Molecule Inhibitors of Tau Aggregation by Targeting Monomeric Tau As a Potential Therapeutic Approach for Tauopathies. Curr. Alzheimer Res. 2015, 12, 814–828. [Google Scholar] [CrossRef] [PubMed]
  36. Joshi, P.; Chia, S.; Habchi, J.; Knowles, T.P.J.; Dobson, C.M.; Vendruscolo, M. A Fragment-Based Method of Creating Small-Molecule Libraries to Target the Aggregation of Intrinsically Disordered Proteins. Acs Comb. Sci. 2016, 18, 144–153. [Google Scholar] [CrossRef]
  37. Yu, C.; Niu, X.G.; Jin, F.; Liu, Z.R.; Jin, C.W.; Lai, L.H. Structure-based Inhibitor Design for the Intrinsically Disordered Protein c-Myc. Sci. Rep. 2016, 6, 1–11. [Google Scholar] [CrossRef]
  38. Bayliss, R.; Burgess, S.G.; Leen, E.; Richards, M.W. A moving target: Structure and disorder in pursuit of Myc inhibitors. Biochem. Soc. Trans. 2017, 45, 709–717. [Google Scholar] [CrossRef]
  39. Chong, B.; Li, M.D.; Li, T.; Yu, M.; Zhang, Y.G.; Liu, Z.R. Conservation of Potentially Druggable Cavities in Intrinsically Disordered Proteins. Acs Omega 2018, 3, 15643–15652. [Google Scholar] [CrossRef]
  40. De Mol, E.; Fenwick, R.B.; Phang, C.T.W.; Buzon, V.; Szulc, E.; de la Fuente, A.; Escobedo, A.; Garcia, J.; Bertoncini, C.W.; Estebanez-Perpina, E.; et al. EPI-001, A Compound Active against Castration-Resistant Prostate Cancer, Targets Transactivation Unit 5 of the Androgen Receptor. ACS Chem. Biol. 2016, 11, 2499–2505. [Google Scholar] [CrossRef]
  41. Toth, G.; Gardai, S.J.; Zago, W.; Bertoncini, C.W.; Cremades, N.; Roy, S.L.; Tambe, M.A.; Rochet, J.C.; Galvagnion, C.; Skibinski, G.; et al. Targeting the Intrinsically Disordered Structural Ensemble of alpha-Synuclein by Small Molecules as a Potential Therapeutic Strategy for Parkinson’s Disease. PLoS ONE 2014, 9, e8713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Nygren, P.J.; Mehta, S.; Schweppe, D.K.; Langeberg, L.K.; Whiting, J.L.; Weisbrod, C.R.; Bruce, J.E.; Zhang, J.; Veesler, D.; Scott, J.D. Intrinsic disorder within AKAP79 fine-tunes anchored phosphatase activity toward substrates and drug sensitivity. Elife 2017, 6. [Google Scholar] [CrossRef]
  43. Basu, S.; Soderquist, F.; Wallner, B. Proteus: A random forest classifier to predict disorder-to-order transitioning binding regions in intrinsically disordered proteins. J. Comput.-Aided Mol. Des. 2017, 31, 453–466. [Google Scholar] [CrossRef]
  44. Heller, G.T.; Bonomi, M.; Vendruscolo, M. Structural Ensemble Modulation upon Small-Molecule Binding to Disordered Proteins. J. Mol. Biol. 2018, 430, 2288–2292. [Google Scholar] [CrossRef] [PubMed]
  45. Ruan, H.; Sun, Q.; Zhang, W.L.; Liu, Y.; Lai, L.H. Targeting intrinsically disordered proteins at the edge of chaos. Drug Discov. Today 2019, 24, 217–227. [Google Scholar] [CrossRef]
  46. Rezaei-Ghaleh, N.; Blackledge, M.; Zweckstetter, M. Intrinsically Disordered Proteins: From Sequence and Conformational Properties toward Drug Discovery. ChemBioChem 2012, 13, 930–950. [Google Scholar] [CrossRef] [PubMed]
  47. Dunker, A.K.; Uversky, V.N. Drugs for ’protein clouds’: Targeting intrinsically disordered transcription factors. Curr. Opin. Pharmacol. 2010, 10, 782–788. [Google Scholar] [CrossRef] [PubMed]
  48. Metallo, S.J. Intrinsically disordered proteins are potential drug targets. Curr. Opin. Chem. Biol. 2010, 14, 481–488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Uversky, V.N. Intrinsically disordered proteins and novel strategies for drug discovery. Expert Opin. Drug Discov. 2012, 7, 475–488. [Google Scholar] [CrossRef]
  50. Chen, C.Y.C.; Tou, W.L. How to design a drug for the disordered proteins. Drug Discov. Today 2013, 18, 910–915. [Google Scholar] [CrossRef]
  51. Heller, G.T.; Sormanni, P.; Vendruscolo, M. Targeting disordered proteins with small molecules using entropy. Trends Biochem. Sci. 2015, 40, 491–496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Joshi, P.; Vendruscolo, M. Druggability of Intrinsically Disordered Proteins. In Intrinsically Disordered Proteins Studied by Nmr Spectroscopy; Felli, I.C., Pierattelli, R., Eds.; Springer: Cham, Switzerland, 2015; Volume 870, pp. 383–400. [Google Scholar]
  53. Sammak, S.; Zinzalla, G. Targeting protein-protein interactions (PPIs) of transcription factors: Challenges of intrinsically disordered proteins (IDPs) and regions (IDRs). Prog. Biophys. Mol. Biol. 2015, 119, 41–46. [Google Scholar] [CrossRef] [PubMed]
  54. Ambadipudi, S.; Zweckstetter, M. Targeting intrinsically disordered proteins in rational drug discovery. Expert Opin. Drug Discov. 2016, 11, 65–77. [Google Scholar] [CrossRef] [PubMed]
  55. Chang, C.K.; Lo, S.C.; Wang, Y.S.; Hou, M.H. Recent insights into the development of therapeutics against coronavirus diseases by targeting N protein. Drug Discov. Today 2016, 21, 562–572. [Google Scholar] [CrossRef]
  56. Varadi, M.; Kosol, S.; Lebrun, P.; Valentini, E.; Blackledge, M.; Dunker, A.K.; Felli, I.C.; Forman-Kay, J.D.; Kriwacki, R.W.; Pierattelli, R.; et al. pE-DB: A database of structural ensembles of intrinsically disordered and of unfolded proteins. Nucleic Acids Res. 2014, 42, D326–D335. [Google Scholar] [CrossRef] [Green Version]
  57. Wright, P.E.; Dyson, H.J. Intrinsically disordered proteins in cellular signalling and regulation. Nat. Rev. Mol. Cell Biol. 2015, 16, 18–29. [Google Scholar] [CrossRef]
  58. Felli, I.C.; Pierattelli, R. Intrinsically Disordered Proteins Studied by NMR Spectroscopy, 1st ed.; Springer International Publishing: Cham, Switzerland, 2015; Volume 870. [Google Scholar]
  59. Ganguly, D.; Chen, J. Structural interpretation of paramagnetic relaxation enhancement-derived distances for disordered protein states. J. Mol. Biol. 2009, 390, 467–477. [Google Scholar] [CrossRef]
  60. Bonomi, M.; Heller, G.T.; Camilloni, C.; Vendruscolo, M. Principles of protein structural ensemble determination. Curr. Opin. Struct. Biol. 2017, 42, 106–116. [Google Scholar] [CrossRef] [Green Version]
  61. Ravera, E.; Sgheri, L.; Parigi, G.; Luchinat, C. A critical assessment of methods to recover information from averaged data. Phys. Chem. Chem. Phys. 2016, 18, 5686–5701. [Google Scholar] [CrossRef]
  62. Marsh, J.A.; Neale, C.; Jack, F.E.; Choy, W.Y.; Lee, A.Y.; Crowhurst, K.A.; Forman-Kay, J.D. Improved structural characterizations of the drkN SH3 domain unfolded state suggest a compact ensemble with native-like and non-native structure. J. Mol. Biol. 2007, 367, 1494–1510. [Google Scholar] [CrossRef]
  63. Salmon, L.; Nodet, G.; Ozenne, V.; Yin, G.W.; Jensen, M.R.; Zweckstetter, M.; Blackledge, M. NMR Characterization of Long-Range Order in Intrinsically Disordered Proteins. J. Am. Chem. Soc. 2010, 132, 8407–8418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Wells, M.; Tidow, H.; Rutherford, T.J.; Markwick, P.; Jensen, M.R.; Mylonas, E.; Svergun, D.I.; Blackledge, M.; Fersht, A.R. Structure of tumor suppressor p53 and its intrinsically disordered N-terminal transactivation domain. Proc. Natl. Acad. Sci. USA 2008, 105, 5762–5767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Berlin, K.; Castaneda, C.A.; Schneidman-Duhovny, D.; Sali, A.; Nava-Tudela, A.; Fushman, D. Recovering a Representative Conformational Ensemble from Underdetermined Macromolecular Structural Data. J. Am. Chem. Soc. 2013, 135, 16595–16609. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Krueger, S.; Shin, J.H.; Raghunandan, S.; Curtis, J.E.; Kelman, Z. Atomistic Ensemble Modeling and Small-Angle Neutron Scattering of Intrinsically Disordered Protein Complexes: Applied to Minichromosome Maintenance Protein. Biophys. J. 2011, 101, 2999–3007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Perkins, S.J.; Wright, D.W.; Zhang, H.L.; Brookes, E.H.; Chen, J.H.; Irving, T.C.; Krueger, S.; Barlow, D.J.; Edler, K.J.; Scott, D.J.; et al. Atomistic modelling of scattering data in the Collaborative Computational Project for Small Angle Scattering (CCP-SAS). J. Appl. Crystallogr. 2016, 49, 1861–1875. [Google Scholar] [CrossRef] [PubMed]
  68. Zhang, W.H.; Howell, S.C.; Wright, D.W.; Heindel, A.; Qiu, X.Y.; Chen, J.H.; Curtis, J.E. Combined Monte Carlo/torsion-angle molecular dynamics for ensemble modeling of proteins, nucleic acids and carbohydrates. J. Mol. Graphics Modell. 2017, 73, 179–190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Kragelj, J.; Ozenne, V.; Blackledge, M.; Jensen, M.R. Conformational propensities of intrinsically disordered proteins from NMR chemical shifts. ChemPhysChem 2013, 14, 3034–3045. [Google Scholar] [CrossRef]
  70. Fisher, C.K.; Huang, A.; Stultz, C.M. Modeling intrinsically disordered proteins with bayesian statistics. J. Am. Chem. Soc. 2010, 132, 14919–14927. [Google Scholar] [CrossRef]
  71. Bertini, I.; Giachetti, A.; Luchinat, C.; Parigi, G.; Petoukhov, M.V.; Pierattelli, R.; Ravera, E.; Svergun, D.I. Conformational space of flexible biological macromolecules from average data. J. Am. Chem. Soc. 2010, 132, 13553–13558. [Google Scholar] [CrossRef]
  72. Levine, Z.A.; Shea, J.E. Simulations of disordered proteins and systems with conformational heterogeneity. Curr. Opin. Struct. Biol. 2017, 43, 95–103. [Google Scholar] [CrossRef] [Green Version]
  73. Knott, M.; Best, R.B. A preformed binding interface in the unbound ensemble of an intrinsically disordered protein: Evidence from molecular simulations. PLoS Comput. Biol. 2012, 8, e1002605. [Google Scholar] [CrossRef]
  74. Mao, A.H.; Crick, S.L.; Vitalis, A.; Chicoine, C.L.; Pappu, R.V. Net charge per residue modulates conformational ensembles of intrinsically disordered proteins. Proc. Natl. Acad. Sci. USA 2010, 107, 8183–8188. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Ganguly, D.; Chen, J. Atomistic details of the disordered states of KID and pKID. implications in coupled binding and folding. J. Am. Chem. Soc. 2009, 131, 5214–5223. [Google Scholar] [CrossRef] [PubMed]
  76. Zhang, W.; Ganguly, D.; Chen, J. Residual structures, conformational fluctuations, and electrostatic interactions in the synergistic folding of two intrinsically disordered proteins. PLoS Comput. Biol. 2012, 8, e1002353. [Google Scholar] [CrossRef] [Green Version]
  77. Brooks, B.R.; Brooks, C.L.; Mackerell, A.D.; Nilsson, L.; Petrella, R.J.; Roux, B.; Won, Y.; Archontis, G.; Bartels, C.; Boresch, S.; et al. CHARMM: The Biomolecular Simulation Program. J. Comput. Chem. 2009, 30, 1545–1614. [Google Scholar] [CrossRef]
  78. Case, D.A.; Cerutti, D.S.; Cheatham III, T.E.; Darden, T.A.; Duke, R.E.; Giese, T.J.; Gohlke, H.; Goetz, A.W.; Greene, D.; Homeyer, N.; et al. Kollman AMBER 2017; University of California: San Francisco, CA, USA, 2017. [Google Scholar]
  79. Eastman, P.; Friedrichs, M.S.; Chodera, J.D.; Radmer, R.J.; Bruns, C.M.; Ku, J.P.; Beauchamp, K.A.; Lane, T.J.; Wang, L.-P.; Shukla, D.; et al. OpenMM 4: A Reusable, Extensible, Hardware Independent Library for High Performance Molecular Simulation. J. Chem. Theory Comput. 2012, 9, 461–469. [Google Scholar] [CrossRef]
  80. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef] [Green Version]
  81. Phillips, J.C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kal, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Gotz, A.W.; Williamson, M.J.; Xu, D.; Poole, D.; Le Grand, S.; Walker, R.C. Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born. J. Chem. Theory Comput. 2012, 8, 1542–1555. [Google Scholar] [CrossRef] [PubMed]
  83. Zhang, W.H.; Chen, J.H. Accelerate Sampling in Atomistic Energy Landscapes Using Topology-Based Coarse-Grained Models. J. Chem. Theory Comput. 2014, 10, 918–923. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Moritsugu, K.; Terada, T.; Kidera, A. Scalable free energy calculation of proteins via multiscale essential sampling. J. Chem. Phys. 2010, 133, 12B603. [Google Scholar] [CrossRef] [PubMed]
  85. Sugita, Y.; Okamoto, Y. Replica-exchange molecular dynamics method for protein folding. Chem. Phys. Lett. 1999, 314, 141–151. [Google Scholar] [CrossRef]
  86. Liu, P.; Kim, B.; Friesner, R.A.; Berne, B.J. Replica exchange with solute tempering: A method for sampling biological systems in explicit water. Proc. Natl. Acad. Sci. USA 2005, 102, 13749–13754. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Mittal, A.; Lyle, N.; Harmon, T.S.; Pappu, R.V. Hamiltonian Switch Metropolis Monte Carlo Simulations for Improved Conformational Sampling of Intrinsically Disordered Regions Tethered to Ordered Domains of Proteins. J. Chem. Theory Comput. 2014, 10, 3550–3562. [Google Scholar] [CrossRef] [PubMed]
  88. Peter, E.K.; Shea, J.E. A hybrid MD-kMC algorithm for folding proteins in explicit solvent. Phys. Chem. Chem. Phys. 2014, 16, 6430–6440. [Google Scholar] [CrossRef]
  89. Zhang, C.; Ma, J. Enhanced sampling and applications in protein folding in explicit solvent. J. Chem. Phys. 2010, 132, 244101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Zheng, L.Q.; Yang, W. Practically Efficient and Robust Free Energy Calculations: Double-Integration Orthogonal Space Tempering. J. Chem. Theory Comput. 2012, 8, 810–823. [Google Scholar] [CrossRef]
  91. Wang, L.; Friesner, R.A.; Berne, B.J. Replica exchange with solute scaling: A more efficient version of replica exchange with solute tempering (REST2). J. Phys. Chem. B 2011, 115, 9431–9438. [Google Scholar] [CrossRef] [Green Version]
  92. Robustelli, P.; Piana, S.; Shaw, D.E. Developing a molecular dynamics force field for both folded and disordered protein states. Proc. Natl. Acad. Sci. USA 2018, 115, E4758–E4766. [Google Scholar] [CrossRef] [Green Version]
  93. Huang, J.; Rauscher, S.; Nawrocki, G.; Ran, T.; Feig, M.; de Groot, B.L.; Grubmuller, H.; MacKerell, A.D., Jr. CHARMM36m: An improved force field for folded and intrinsically disordered proteins. Nat. Methods 2017, 14, 71–73. [Google Scholar] [CrossRef] [Green Version]
  94. Robertson, M.J.; Tirado-Rives, J.; Jorgensen, W.L. Improved Peptide and Protein Torsional Energetics with the OPLS-AA Force Field. J. Chem. Theory Comput. 2015, 11, 3499–3509. [Google Scholar] [CrossRef] [PubMed]
  95. Piana, S.; Donchev, A.G.; Robustelli, P.; Shaw, D.E. Water dispersion interactions strongly influence simulated structural properties of disordered protein states. J. Phys. Chem. B 2015, 119, 5113–5123. [Google Scholar] [CrossRef] [PubMed]
  96. Lee, K.H.; Chen, J.H. Optimization of the GBMV2 implicit solvent force field for accurate simulation of protein conformational equilibria. J. Comput. Chem. 2017, 38, 1332–1341. [Google Scholar] [CrossRef] [PubMed]
  97. Nerenberg, P.S.; Jo, B.; So, C.; Tripathy, A.; Head-Gordon, T. Optimizing Solute-Water van der Waals Interactions To Reproduce Solvation Free Energies. J. Phys. Chem. B 2012, 116, 4524–4534. [Google Scholar] [CrossRef] [Green Version]
  98. Choi, J.M.; Pappu, R.V. Improvements to the ABSINTH Force Field for Proteins Based on Experimentally Derived Amino Acid Specific Backbone Conformational Statistics. J. Chem. Theory Comput. 2019, 15, 1367–1382. [Google Scholar] [CrossRef]
  99. Liu, X.; Chen, J. Residual Structures and Transient Long-Range Interactions of p53 Transactivation Domain: Assessment of Explicit Solvent Protein Force Fields. J. Chem. Theory Comput. 2019, 15, 4708–4720. [Google Scholar] [CrossRef]
  100. Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins 2006, 65, 712–725. [Google Scholar] [CrossRef] [Green Version]
  101. Lindorff-Larsen, K.; Piana, S.; Palmo, K.; Maragakis, P.; Klepeis, J.L.; Dror, R.O.; Shaw, D.E. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins 2010, 78, 1950–1958. [Google Scholar] [CrossRef] [Green Version]
  102. Nerenberg, P.S.; Head-Gordon, T. Optimizing Protein-Solvent Force Fields to Reproduce Intrinsic Conformational Preferences of Model Peptides. J. Chem Theory. Comput. 2011, 7, 1220–1230. [Google Scholar] [CrossRef] [Green Version]
  103. Best, R.B.; Zheng, W.; Mittal, J. Balanced Protein-Water Interactions Improve Properties of Disordered Proteins and Non-Specific Protein Association. J. Chem. Theory Comput. 2014, 10, 5113–5124. [Google Scholar] [CrossRef] [Green Version]
  104. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Yu, L.; Li, D.-W.; Brüschweiler, R. Balanced Amino-Acid-Specific Molecular Dynamics Force Field for the Realistic Simulation of Both Folded and Disordered Proteins. J. Chem. Theory Comput. 2020, 16, 1311–1318. [Google Scholar] [CrossRef] [PubMed]
  106. Piana, S.; Lindorff-Larsen, K.; Shaw, D.E. How Robust Are Protein Folding Simulations with Respect to Force Field Parameterization? Biophys. J. 2011, 100, L47–L49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Best, R.B.; Zhu, X.; Shim, J.; Lopes, P.E.M.; Mittal, J.; Feig, M.; MacKerell, A.D. Optimization of the Additive CHARMM All-Atom Protein Force Field Targeting Improved Sampling of the Backbone ϕ, ψ and Side-Chain χ1 and χ2 Dihedral Angles. J. Chem. Theory Comput. 2012, 8, 3257–3273. [Google Scholar] [CrossRef] [Green Version]
  108. Mackerell, A.D.; Feig, M.; Brooks, C.L. Extending the treatment of backbone energetics in protein force fields: Limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations. J. Comput. Chem. 2004, 25, 1400–1415. [Google Scholar] [CrossRef]
  109. Best, R.B. Computational and theoretical advances in studies of intrinsically disordered proteins. Curr. Opin. Struct. Biol. 2017, 42, 147–154. [Google Scholar] [CrossRef]
  110. Ganguly, D.; Chen, J. Modulation of the disordered conformational ensembles of the p53 transactivation domain by cancer-associated mutations. PLoS Comput. Biol. 2015, 11, e1004247. [Google Scholar] [CrossRef] [Green Version]
  111. Pall, S.; Abraham, M.J.; Kutzner, C.; Hess, B.; Lindahl, E. Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS. Lect. Notes Comput. Sci. 2015, 8759, 3–27. [Google Scholar]
  112. Tribello, G.A.; Bonomi, M.; Branduardi, D.; Camilloni, C.; Bussi, G. PLUMED 2: New feathers for an old bird. Comput. Phys. Commun. 2014, 185, 604–613. [Google Scholar] [CrossRef] [Green Version]
  113. Terakawa, T.; Kameda, T.; Takada, S. On Easy Implementation of a Variant of the Replica Exchange with Solute Tempering in GROMACS. J. Comput. Chem. 2011, 32, 1228–1234. [Google Scholar] [CrossRef]
  114. Bussi, G. Hamiltonian replica exchange in GROMACS: A flexible implementation. Mol. Phys. 2014, 112, 379–384. [Google Scholar] [CrossRef] [Green Version]
  115. Chen, J.; Im, W.; Brooks, C.L. Balancing solvation and intramolecular interactions: Toward a consistent generalized born force field. J. Am. Chem. Soc. 2006, 128, 3728–3736. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  116. Gong, X.; Chiricotto, M.; Liu, X.; Nordquist, E.; Feig, M.; Brooks, C.L., 3rd; Chen, J. Accelerating the Generalized Born with Molecular Volume and Solvent Accessible Surface Area Implicit Solvent Model Using Graphics Processing Units. J. Comput. Chem. 2020, 41, 830–838. [Google Scholar] [CrossRef] [PubMed]
  117. Wojcik, S.; Birol, M.; Rhoades, E.; Miranker, A.D.; Levine, Z.A. Targeting the Intrinsically Disordered Proteome Using Small-Molecule Ligands. In Intrinsically Disordered Proteins; Rhoades, E., Ed.; Elsevier Inc.: London, UK, 2018; Volume 611, pp. 703–734. [Google Scholar]
  118. Tompa, P.; Fuxreiter, M. Fuzzy complexes: Polymorphism and structural disorder in protein-protein interactions. Trends Biochem. Sci. 2008, 33, 2–8. [Google Scholar] [CrossRef]
  119. McDowell, C.; Chen, J.; Chen, J. Potential Conformational Heterogeneity of p53 Bound to S100B(betabeta). J. Mol. Biol. 2013, 425, 999–1010. [Google Scholar] [CrossRef] [Green Version]
  120. Weng, J.; Wang, W. Dynamic multivalent interactions of intrinsically disordered proteins. Curr. Opin. Struct. Biol. 2019, 62, 9–13. [Google Scholar] [CrossRef]
  121. Michel, J.; Cuchillo, R. The impact of small molecule binding on the energy landscape of the intrinsically disordered protein C-myc. PLoS ONE 2012, 7, e41070. [Google Scholar] [CrossRef] [Green Version]
  122. Wang, H.; Hammoudeh, D.I.; Follis, A.V.; Reese, B.E.; Lazo, J.S.; Metallo, S.J.; Prochownik, E.V. Improved low molecular weight Myc-Max inhibitors. Mol. Cancer Ther. 2007, 6, 2399–2408. [Google Scholar] [CrossRef] [Green Version]
  123. Jin, F.; Yu, C.; Lai, L.; Liu, Z. Ligand clouds around protein clouds: A scenario of ligand binding with intrinsically disordered proteins. PLoS Comput. Biol. 2013, 9, e1003249. [Google Scholar] [CrossRef] [Green Version]
  124. Heller, G.T.; Aprile, F.A.; Bonomi, M.; Camilloni, C.; De Simone, A.; Vendruscolo, M. Sequence Specificity in the Entropy-Driven Binding of a Small Molecule and a Disordered Peptide. J. Mol. Biol. 2017, 429, 2772–2779. [Google Scholar] [CrossRef]
  125. Daniels, M.J.; Nourse, J.B.; Kim, H.; Sainati, V.; Schiavina, M.; Murrali, M.G.; Pan, B.Y.; Ferrie, J.J.; Haney, C.M.; Moons, R.; et al. Cyclized NDGA modifies dynamic alpha-synuclein monomers preventing aggregation and toxicity. Sci. Rep. 2019, 9, 1–17. [Google Scholar]
  126. Liang, C.; Savinov, S.N.; Fejzo, J.; Eyles, S.J.; Chen, J. Modulation of Amyloid-beta42 Conformation by Small Molecules Through Nonspecific Binding. J. Chem. Theory Comput. 2019, 15, 5169–5174. [Google Scholar] [CrossRef]
  127. Liu, X.; Chen, J. Modulation of p53 Transactivation Domain Conformations by Ligand Binding and Cancer-Associated Mutations. Pac. Symp. Biocomput. 2020, 25, 195–206. [Google Scholar]
  128. Du, G.J.; Zhang, Z.; Wen, X.D.; Yu, C.; Calway, T.; Yuan, C.S.; Wang, C.Z. Epigallocatechin Gallate (EGCG) is the most effective cancer chemopreventive polyphenol in green tea. Nutrients 2012, 4, 1679–1691. [Google Scholar] [CrossRef]
  129. Fujiki, H.; Sueoka, E.; Watanabe, T.; Suganuma, M. Primary cancer prevention by green tea, and tertiary cancer prevention by the combination of green tea catechins and anticancer compounds. J. Cancer Prev. 2015, 20, 1–4. [Google Scholar] [CrossRef] [PubMed]
  130. Shin, C.M.; Lee, D.H.; Seo, A.Y.; Lee, H.J.; Kim, S.B.; Son, W.C.; Kim, Y.K.; Lee, S.J.; Park, S.H.; Kim, N.; et al. Green tea extracts for the prevention of metachronous colorectal polyps among patients who underwent endoscopic removal of colorectal adenomas: A randomized clinical trial. Clin. Nutr. 2018, 37, 452–458. [Google Scholar] [CrossRef] [PubMed]
  131. Bieschke, J.; Russ, J.; Friedrich, R.P.; Ehrnhoefer, D.E.; Wobst, H.; Neugebauer, K.; Wanker, E.E. EGCG remodels mature alpha-synuclein and amyloid-beta fibrils and reduces cellular toxicity. Proc. Natl. Acad. Sci. USA 2010, 107, 7710–7715. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  132. Rezai-Zadeh, K.; Arendash, G.W.; Hou, H.; Fernandez, F.; Jensen, M.; Runfeldt, M.; Shytle, R.D.; Tan, J. Green tea epigallocatechin-3-gallate (EGCG) reduces beta-amyloid mediated cognitive impairment and modulates tau pathology in Alzheimer transgenic mice. Brain Res. 2008, 1214, 177–187. [Google Scholar] [CrossRef]
  133. Santofimia-Castano, P.; Rizzuti, B.; Xia, Y.; Abian, O.; Peng, L.; Velazquez-Campoy, A.; Iovanna, J.L.; Neira, J.L. Designing and repurposing drugs to target intrinsically disordered proteins for cancer treatment: Using NUPR1 as a paradigm. Mol. Cell. Oncol. 2019, 6, e1612678. [Google Scholar] [CrossRef]
  134. Neira, J.L.; Bintz, J.; Arruebo, M.; Rizzuti, B.; Bonacci, T.; Vega, S.; Lanas, A.; Velazquez-Campoy, A.; Iovanna, J.L.; Abian, O. Identification of a Drug Targeting an Intrinsically Disordered Protein Involved in Pancreatic Adenocarcinoma. Sci. Rep. 2017, 7, 39732. [Google Scholar] [CrossRef]
  135. Alberti, S.; Gladfelter, A.; Mittag, T. Considerations and Challenges in Studying Liquid-Liquid Phase Separation and Biomolecular Condensates. Cell 2019, 176, 419–434. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  136. Holehouse, A.S.; Pappu, R.V. Functional Implications of Intracellular Phase Transitions. Biochemistry 2018, 57, 2415–2423. [Google Scholar] [CrossRef] [PubMed]
  137. Brangwynne, C.P.; Tompa, P.; Pappu, R.V. Polymer physics of intracellular phase transitions. Nat. Phys. 2015, 11, 899–904. [Google Scholar] [CrossRef]
  138. Heller, G.T.; Aprile, F.A.; Vendruscolo, M. Methods of probing the interactions between small molecules and disordered proteins. Cell. Mol. Life Sci. 2017, 74, 3225–3243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  139. Dobrev, V.S.; Fred, L.M.; Gerhart, K.P.; Metallo, S.J. Characterization of the Binding of Small Molecules to Intrinsically Disordered Proteins. In Intrinsically Disordered Proteins; Rhoades, E., Ed.; Elsevier Inc.: London, UK, 2018; Volume 611, pp. 677–702. [Google Scholar]
Figure 1. (a). Evolution of the per-residue β-sheet structure during a 30-μs Anton MD simulation of intrinsically disordered Aβ40 in explicit solvent at 300 K. Residues assigned to be in the β-sheet conformation are colored in red. (b) Average residue helix and (c) β-sheet probability profiles derived from the first and second halves of the trajectory. Note that both pairs of profiles differ significantly, reflecting a lack of convergence in the simulated disordered ensemble. The original MD trajectory was generously provided by D. E. Shaw Research [92].
Figure 1. (a). Evolution of the per-residue β-sheet structure during a 30-μs Anton MD simulation of intrinsically disordered Aβ40 in explicit solvent at 300 K. Residues assigned to be in the β-sheet conformation are colored in red. (b) Average residue helix and (c) β-sheet probability profiles derived from the first and second halves of the trajectory. Note that both pairs of profiles differ significantly, reflecting a lack of convergence in the simulated disordered ensemble. The original MD trajectory was generously provided by D. E. Shaw Research [92].
Biomolecules 10 00743 g001
Figure 2. Calculated (lines) and experimental (grey bars) NMR Paramagnetic Relaxation Enhancement (PRE) effects induced by paramagnetic spin labelling at residues D7, E28, A39, and D61 of p53-TAD. Red and green traces were calculated from an independent control and folding REST2 simulations of p53-TAD in a99SB-disp, respectively, to evaluate the level of convergence. Control and folding simulations were initiated from helical and fully unstructured structures, respectively, and the length of REST2 simulations were 1 μs per replica. This figure was adapted from [99]. See [99] for details on the simulation and analysis.
Figure 2. Calculated (lines) and experimental (grey bars) NMR Paramagnetic Relaxation Enhancement (PRE) effects induced by paramagnetic spin labelling at residues D7, E28, A39, and D61 of p53-TAD. Red and green traces were calculated from an independent control and folding REST2 simulations of p53-TAD in a99SB-disp, respectively, to evaluate the level of convergence. Control and folding simulations were initiated from helical and fully unstructured structures, respectively, and the length of REST2 simulations were 1 μs per replica. This figure was adapted from [99]. See [99] for details on the simulation and analysis.
Biomolecules 10 00743 g002
Figure 3. Conformational ensembles of Aβ42 with and without the ligands (ac) and p53-TAD with and without ligands (d,e). The conformational space of Aβ42 is projected onto the number of backbone hydrogen bonds and end-to-end distance and that of p53-TAD is projected onto the first two principal components. The conformational ensembles were calculated using long timescale REST2 simulations in explicit solvent (10 and 1 μs per replica for Aβ42 and p53-TAD, respectively). Representative conformations are shown in backbone traces. This figure was adapted from [126,127]. See [126,127] for details on the simulation and analysis.
Figure 3. Conformational ensembles of Aβ42 with and without the ligands (ac) and p53-TAD with and without ligands (d,e). The conformational space of Aβ42 is projected onto the number of backbone hydrogen bonds and end-to-end distance and that of p53-TAD is projected onto the first two principal components. The conformational ensembles were calculated using long timescale REST2 simulations in explicit solvent (10 and 1 μs per replica for Aβ42 and p53-TAD, respectively). Representative conformations are shown in backbone traces. This figure was adapted from [126,127]. See [126,127] for details on the simulation and analysis.
Biomolecules 10 00743 g003

Share and Cite

MDPI and ACS Style

Chen, J.; Liu, X.; Chen, J. Targeting Intrinsically Disordered Proteins through Dynamic Interactions. Biomolecules 2020, 10, 743. https://doi.org/10.3390/biom10050743

AMA Style

Chen J, Liu X, Chen J. Targeting Intrinsically Disordered Proteins through Dynamic Interactions. Biomolecules. 2020; 10(5):743. https://doi.org/10.3390/biom10050743

Chicago/Turabian Style

Chen, Jianlin, Xiaorong Liu, and Jianhan Chen. 2020. "Targeting Intrinsically Disordered Proteins through Dynamic Interactions" Biomolecules 10, no. 5: 743. https://doi.org/10.3390/biom10050743

APA Style

Chen, J., Liu, X., & Chen, J. (2020). Targeting Intrinsically Disordered Proteins through Dynamic Interactions. Biomolecules, 10(5), 743. https://doi.org/10.3390/biom10050743

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop