Looking at the Disordered Proteins through the Computational

Mar 22, 2018 - Intrinsically disordered proteins (IDPs) have attracted wide interest over the past decade due to their surprising prevalence in the pr...
1 downloads 4 Views 2MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Outlook Cite This: ACS Cent. Sci. XXXX, XXX, XXX−XXX

Looking at the Disordered Proteins through the Computational Microscope Payel Das,*,† Silvina Matysiak,‡ and Jeetain Mittal§ †

IBM Thomas J. Watson Research Center, Yorktown Heights, New York 10598, United States Fischell Department of Bioengineering, University of Maryland, College Park, Maryland 20742, United States § Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States ‡

ABSTRACT: Intrinsically disordered proteins (IDPs) have attracted wide interest over the past decade due to their surprising prevalence in the proteome and versatile roles in cell physiology and pathology. A large selection of IDPs has been identified as potential targets for therapeutic intervention. Characterizing the structure− function relationship of disordered proteins is therefore an essential but daunting task, as these proteins can adapt transient structure, necessitating a new paradigm for connecting structural disorder to function. Molecular simulation has emerged as a natural complement to experiments for atomic-level characterizations and mechanistic investigations of this intriguing class of proteins. The diverse range of length and time scales involved in IDP function requires performing simulations at multiple levels of resolution. In this Outlook, we focus on summarizing available simulation methods, along with a few interesting example applications. We also provide an outlook on how these simulation methods can be further improved in order to provide a more accurate description of IDP structure, binding, and assembly.



INTRODUCTION Intrinsically disordered proteins (IDPs) are a class of proteins that are an exception to the widely accepted Anfinsen’s sequence-structure−function dogma, in which the amino acid sequence of a protein defines its unique three-dimensional fold that is inherently linked to the function. Instead, IDPs exist as a diverse ensemble of rapidly interconverting structures.1 The structural disorder, unstructuredness, or fuzziness is abundant in eukaryotic proteome. About 25% of the eukaryotic proteins are predicted to be IDPs, and as high as 55% of eukaryotic proteins are predicted to consist of an intrinsically disordered region (IDR) longer than 30 residues.2 IDPs, typically composed of 30 to 1000 residues, exhibit significant chain compaction, contrary to the typical random-coil-like behavior of chemically denatured proteins.3 The function of an IDP often involves interaction with partners resulting in “fuzzy complexes”,4 which encompasses molecular recognition, membraneless organelle formation, amyloidogenesis, cellular signaling, and regulation.5 In addition, a number of IDPs are pathological targets in the age-associated amyloid disorders.6 The known examples include Amyloid beta peptide (Aβ) and tau peptide in Alzheimer’s disease (AD), α-synuclein in Parkinson’s disease (PD) and synucleopathies, huntingtin protein in Huntington’s disease (HD), and islet amyloid polypeptide (hIAPP) in type-2 diabetes. Accumulation of the intracellular and extracellular protein aggregates, often with amyloid-like characteristics, is a pathological hallmark of those disorders.7,8 The vast literature supports higher toxicity of the prefibrillar, early aggregate species compared to the mature fibril.9−14 In healthy individuals, formation of those potentially © XXXX American Chemical Society

toxic aggregates is prevented via cellular protection mechanism(s), such as molecular chaperones and protein degradation machinery.15 The detailed structural characterization of disordered proteins using conventional experimental techniques remains extremely challenging, due to their structural heterogeneity and dynamic nature. IDPs undergo conformational transitions at a time scale spanning several orders of magnitude (ps-ms), which is often difficult to address using a single experimental approach. Smallangle X-ray scattering (SAXS)16 and single molecule Forster resonance energy transfer (FRET)17 experiments have provided valuable information on the chain dimension, while NMR18 has been useful in extracting secondary structure. However, the interpretation of those ensemble-averaged experimental observations has been hindered by the lack of a molecular model. Additionally, identification of transient structural epitope(s) within an IDP ensemble appears to be crucial, as those might have functional implication or disease relevance. Molecular simulations play an essential complementary role by allowing generation of sufficiently accurate structural ensembles that can directly interpret and/or predict experimental results. Simulations give mechanistic and dynamic insights into the IDP assembly as well as binding with partners. Additionally, simulations provide an in silico platform for designing and testing therapeutic strategies. In this Outlook, we will summarize the existing computational approaches for IDP structure characterization along with some Received: December 31, 2017

A

DOI: 10.1021/acscentsci.7b00626 ACS Cent. Sci. XXXX, XXX, XXX−XXX

ACS Central Science

Outlook

The detailed structural characterization of disordered proteins using conventional experimental techniques remains extremely challenging, due to their structural heterogeneity and dynamic nature.

In this Outlook, we will summarize the existing computational approaches for IDP structure characterization along with some example applications, accentuate their strengths and limitations, and comment on future directions for improving those approaches.

example applications, accentuate their strengths and limitations, and comment on future directions for improving those approaches. Computational approaches for generating IDP conformational ensembles can be essentially classified into two primary categories: (1) knowledge-based methods that use experimental data, mostly NMR, and (2) first principle or de novo approaches that combine empirical force fields with Molecular Dynamics or Monte Carlo simulations. Before getting into the details, it must be mentioned that these methods have been used for studying folded globular proteins for years. However, an extensive amount of research shows the need for modifications and improvements of existing tools, as well as calls for development of new methods.

oligomeric form (for a review of recent work on this topic see ref 33 and references cited hereafter). The central idea behind conducting such simulations is to identify the intrinsic propensities of the IDP sequence and changes thereof due to disease mutations, environmental factors (solvent quality, proximity to lipid membranes, post-translational modification, etc.), which can affect function, aggregation, and toxicity. Thermodynamics and kinetics of the higher-order assembly of shorter amyloidogenic peptide fragments have been also investigated.34−36 Simulations have been extensively used to investigate IDP interactions with various partners, for example, with small molecule,37,38 peptide inhibitor,39 membrane,40,41 and chaperone. 42 Use of a higher than experimental concentration of the peptide and/or enhanced sampling protocols is becoming routine to ensure adequate sampling of an IDP’s conformational landscape. While standard molecular dynamics simulations can provide both structural and dynamical information on IDPs, sufficient convergence of the equilibrated structural ensembles can be more efficiently accomplished using enhanced sampling techniques, such as replica exchange molecular dynamics,43,44 metadynamics,45 and related methods, e.g. bias-exchange metadynamics,46 solute tempering,47−49 etc. Recently, parallel tempering in the welltempered ensemble (PT-WTE) method has become available, which allows for temperature based sampling enhancement with reduced computational cost, as compared to standard REMD.50,51 In the following, we will discuss examples of how simulations can help in characterizing IDP ensembles, both native and modified. Transient Structure Identification. Transient local structure formation within monomeric IDP has been postulated to impact aggregation, loss-of-function, and toxicity, which can be efficiently investigated using de novo simulations. For example, atomistic simulations of the monomeric TDP-43 protein revealed local helix stabilization at the C-terminus in the presence of ALS mutations, consistent with NMR experiments.52 Another simulation study showed that the known aggregation propensity of four amylin sequences, implicated in type-2 diabetes, can potentially be explained by their local helical propensity differences.53 A computationally well-studied IDP is Aβ peptide that is attractive due to its length (40−43 residues), disease-relevance (Alzheimer’s), and wealth of available experimental data on the wild-type (WT) and mutated sequences (Figure 1). A number of familial mutations reported in Aβ sequence result in early onset of AD. Simulations from several groups have shed light onto the structural alterations resulting from those point mutations in the full-length and fragment Aβ, particularly within the monomeric and dimeric landscape. In the following, we will review in detail the differential effect of two recently



KNOWLEDGE-BASED APPROACHES In the knowledge-based approaches, random coil ensembles are first generated using methods such as TraDES19 and FlexibleMeccano.20 The generated ensemble is then reweighted to match experimentally derived conformational constraints. Examples include the Monte Carlo optimization approach (e.g., ENSEMBLE21 and ASTEROIDS22), Bayesian weighting approach,23 and maximum entropy approach.24 Albeit useful, lowering the generation bias toward statistical coil ensembles, improvements in back-calculation algorithms, and simultaneous optimization with respect to local, global, and dynamic experimental observables are some of the improvements needed for accurate identification of existing structures in an IDP ensemble.25



DE NOVO SIMULATIONS Molecular simulation grounded in statistical mechanical principles is a compelling alternative to the knowledge-based approaches, which generates a Boltzmann-weighted ensemble of conformational subpopulations irrespective of experimental data. It is a common practice to compare the simulationgenerated ensemble with experimental data by back-calculating the structural (e.g., chain dimension) and dynamical (e.g., relaxation) properties. Once validated, simulations predict a complete description of the IDP ensemble and allow further analysis of important transient structures. • Explicit Solvent Simulations. The most accurate way to model a system is by including both protein and water at their atomistic detail. Different combinations of popular protein force fields (GROMOS, 26 OPLS/AA, 27 CHARMM, 28 AMBER29) and water models (TIP3P,30 SPC,31 SPC/E32) are commonly used for this purpose. However, the large spatiotemporal scales associated with IDP dynamics and aggregation often make the conventional atomistic explicit solvent simulations computationally unaffordable. Until recently, full-blown atomistic simulations in explicit water have been successfully used to study the structural landscape of a number of IDP sequences in their monomeric and small B

DOI: 10.1021/acscentsci.7b00626 ACS Cent. Sci. XXXX, XXX, XXX−XXX

ACS Central Science

Outlook

studied using atomistic simulations.64−68 While aggregationprone WT/WT and A2V/A2V homodimers are found to populate highly stabilized interpeptide structures with a high content of interpeptide β-sheets or β-hairpin conformations, the AD-protective WT/A2V and WT/A2T Aβ42 heterotypical interactions favored formation of weak dimers that are highly disordered and collapsed in nature.67−69 These results suggest that A2V and A2T mutations differentially alter the misfolding into the transient monomeric β-hairpin structures that are thought to play a crucial role in aggregation, as suggested by earlier simulations70,71 and experiments.72 The overly stabilized double β-hairpin structures might explain the longer lag phase as well as the atypical spherical oligomer population during A2V Aβ42 aggregation,73 which might be linked to its diseasecausing effect in homozygous individuals. β-Hairpin lowering in the A2T Aβ42 ensemble is interesting, given its relatively lower long-term potentiation (LTP) inhibition ability observed in rat CA1 hippocampal CA1 neurons73 and the reported toxicity relevance of a monomeric β-hairpin structure.72 The weak dimer formation observed in protective WT/A2V and WT/ A2T mixture simulations agrees with the oligomeric properties reported in FT-IR and ANS fluorescence experiments.74 Additionally, the enhanced disorder in the heterodimeric population might explain the absence of toxic dodecamer in the protective WT/A2T(A2V) mixture, as revealed in IM-MS experiments.75 Taken together, existing simulations imply that the mutated Aβ N-terminus can modulate an aggregationprone, toxic β-hairpin population in a sequence-dependent manner. Toward Structure-Based Therapeutic Design. Designing Aβ-targeting small molecules,77−79 peptide fragments,80,39,81 or antibodies82 is an active area of AD research due to the widely accepted amyloid cascade hypothesis. It has been suggested that a better understanding of structure-toxicity relationships, much earlier detection of AD, and drugs that consider multiple targets (Aβ, tau) and can pass blood-brain barrier are needed for effective drug design.83 Atomistic simulations are constantly being employed in designing mechanism-based antiamyloid molecules against monomeric, oligomeric, and fibrillated forms of Aβ. The atypical interactions triggered by an altered Nterminus, as seen within monomeric and dimeric Aβ, call for further simulations toward Aβ ΝΤR-based therapeutic design.

Figure 1. Sequence (top) and simulated ensemble (bottom) of the wild-type (WT), causative A2V, and protective A2T Aβ42 variant monomer.54

reported N-terminal mutations, one AD-causative and one protective, on the monomeric and dimeric Aβ landscapes (Figure 1). The A2V mutation is associated with the early onset of AD in homozygous carriers, whereas it offers some protection in the heterozygous state.55 The A2T mutation reportedly has an AD-protective effect.56 Correspondingly, a number of atomistic REMD studies using different combination of protein force fields and water models have reported differences in the monomeric configuration between the wildtype, A2V, and A2T variants of full-length Aβ peptide as well as of Aβ fragments.54,57,58 Although all three variant monomers were found to form collapsed, disordered ensembles in water (Figure 1, bottom), they were strikingly different in terms of transient secondary and tertiary structure formation (Figure 2).

Figure 2. Changes in the Aβ42 conformational landscape due to A2V and A2T mutation.54

Taken together, it is becoming more evident that synergistic simulation-experimental investigations are going to play a crucial role in rational design of mechanism-based therapeutics for IDPassociated diseases in the forthcoming days.

Simulations indicate the transient presence of a β-hairpin motif within the WT Aβ monomer, which is comprised of the central and C-terminal hydrophobic residues (Figure 2, magenta rectangle),54,59−61 similar to the β-hairpin structure found in the complex of Aβ40 with an affibody protein ZAβ3.62 The Nterminus forms limited interaction with the rest of the peptide in those short-lived β-hairpin structures. Interestingly, an enhanced population of a double β-hairpin motif involving the NTR has been observed in simulations of a A2V monomer,54,57 which contains a strong resemblance to a solid-state NMR model of neurotoxic Aβ42 oligomers63 (Figure 2, red rectangle). In contrast, the A2T mutation enhances disorder within the Aβ42 ensemble due to atypical long-range interactions between the NTR, CHC, and 22−28 turn region (Figure 2, green rectangle).54,58 The effect of homotypic and heterotypic interactions between WT, A2V, and A2T monomers on the Aβ40/42 dimeric landscape has been also

In fact, recent experiments have reported aggregation and/or toxicity inhibiting property of Aβ N-terminal fragments.84−86 Experiments also showed that linking NTR Aβ fragments with an arginine-rich TAT sequence can facilitate their delivery to brain and cell membrane translocation.84 The 1-6A2VAβ-TAT molecule exhibited robust antiamyloidogenic properties in vitro and protected transgenic C. elegans and human neuroblastoma cells from Aβ toxicity.84 Preclinical studies in AD mouse C

DOI: 10.1021/acscentsci.7b00626 ACS Cent. Sci. XXXX, XXX, XXX−XXX

ACS Central Science

Outlook

Figure 3. Aβ N-terminal hexapeptide binding stabilizes alternative Aβ42 structures.76 (a) Unstructured Aβ42 bound to 1-6A2V and (b) helical Aβ42 bound to 1-6A2T.

charge patterning in the “RAM” disordered region of the Notch receptor can modulate its transcriptional activity.94 • Coarse-Grained Models. Coarse-grained (CG) models further decrease computational requirements (usually by ∼100× compared to fully atomistic explicit solvent simulations) by eliminating few degrees of freedom of the peptide and thus allow simulations of larger length- or time scales. These models are sometimes simulated with the Discrete Molecular Dynamics (DMD) method that involves discontinuous potentials, in order to further gain computational efficiency.95 CG protein models can assume various levels of reduced representation, from a single monomer (peptide unit), to an intermediate level (multiple beads per residue), to a single residue (single bead per residue) resolution. The most coarse-grained are the lattice models in which a single residue occupies a single site of the lattice. While the lattice models allow efficient sampling of the complete phase diagram, due to their inability to capture realistic kinetics and structures they can only be used to answer very generic questions. Therefore, coarse-grained protein models are mostly off-lattice in nature. The interaction parameters in those models are derived either by matching with the atomistic simulations and/or with experimental data, as in the MARTINI model,96 or by using a knowledge-based approach, as in PRIME20.97 The lack of granularity in CG models leads to obvious limitations, such as failure to correctly estimate secondary structure. For example, one limitation of the widely used MARTINI model is the use of secondary and tertiary structure constraints, which prevents exploration of protein conformational fluctuations. It has been shown that inclusion of peptide backbone polarization effects within the MARTINI model framework can alleviate that limitation.98,99 A combination of physics-based and knowledge-based potentials has been also used, such as in OPEP.100 AWSEM (associated memory, water mediated, structure and energy model) developed by Wolynes, Papoian, and co-workers is another CG model that is comprised of both physics-based and knowledge-based potentials.101 Over the years, these CG models have provided valuable insights onto the mechanism of peptide assembly. For example, an AWSEM simulation study of Aβ40 aggregation has characterized the “amyloid funnel” that reveals a prefibrillar to fibrillar interconversion bottleneck at the pentamer level.102 Many coarse-grained models are built on adaptations of the so-called Go-like model that uses a nativecentric potential.103 Further reparametrization and inclusion of additional terms in the energy function were performed to better suit the task of IDP structure modeling, e.g. fine-tuning

models indicated that short-term treatment with 1-6A2VAβTAT impedes Aβ aggregation and cerebral amyloid deposition.84 Inspired by those experiments, the interaction of WT, 16A2V, and 1-6A2T Aβ hexapeptide with the full-length Aβ42 monomer and its effect on the monomeric Aβ42 landscape have been recently investigated using atomistic simulations76 (Figure 3). While all three hexapeptide bindings led to a robust lowering of transient WT hairpin population in simulations, interesting alternative structural features emerged. A2V-bound Aβ42 preferably visited a disordered population (Figure 3a), while a central helix was more frequently populated within the A2T-bound monomer (Figure 3b). Such redirection of an IDP to atypical disordered conformations induced by an inhibitor/ drug interaction may have implications in the lowering of IDP aggregation and toxicity in general, as found in experiments for α-synuclein, tau protein, and Aβ.79,77,87,88 Taken together, it is becoming more evident that synergistic simulation-experimental investigations are going to play a crucial role in rational design of mechanism-based therapeutics for IDP-associated diseases in the forthcoming days. • Implicit Solvent Models. A cheaper alternative way to simulate IDP conformational dynamics is by representing solvent in an implicit manner. In implicit solvent models, the solute−solvent interaction is presented using a mean field approach, in which the solvation free energy is represented as a function of the atomic coordinates of the solute. Approximating the solvation free energy relies on the solvent-accessible surface area (SASA) models or the solvent contact models. The accurate and efficient formulation of the solvent effect remains a nontrivial theoretical challenge, and detailed discussion of this topic is beyond the scope of this Outlook. Apparently, many implicit solvent models suffer from the same problems as the explicit solvent models do, such as overly collapsed disordered states.89 Implicit solvent models that appear promising for IDP simulations are based on the EEF1 model originally developed by Lazaridis and Karplus.90 The EEF1 model decouples and independently estimates the contributions of the charged and nonpolar groups toward the solvation free energy of the solute. A conceptually similar, but much superior, model is the ABSINTH force field of Vitalis and Pappu91 that additionally leverages the strengths of the Poisson framework and better represents polyelectrolyte and polyampholyte IDPs. The ABSINTH model has been extensively and successfully used in recent years to delineate the sequence-structure−function paradigm of IDPs.92,93 For example, ABSINTH simulations combined with experiments have recently uncovered how D

DOI: 10.1021/acscentsci.7b00626 ACS Cent. Sci. XXXX, XXX, XXX−XXX

ACS Central Science

Outlook

protein force field, mainly by using a refined backbone CMAP potential and improved description of salt-bridge interaction, in order to correct the overestimation of left-handed α-helix population. This study also highlights the role of protein−water dispersion interactions in correcting the general bias of physicsbased atomistic models toward overly compact IDP ensembles. Future force field development research should also focus on accurate modeling of the highly dynamic metal-IDP interactions, due to their pivotal roles in neurodegenerative diseases.116 For example, metal coordination is known to alter the aggregation kinetics and morphology of Aβ aggregates.117 Metal ions in unpolarized atomistic simulations are treated as simple van der Waals spheres, nonbonded models with atoms, or with bonded models by employing artificial bonds.118 Nevertheless, the polarization effects are crucial, as the polarization energy can range from very small values to values around 10−20% for the total interaction energy.119 Therefore, a number of polarizable modeling approaches are continuously being developed, such as fluctuating charge, Drude oscillator, and the induced dipole models.118 However, the usage of polarizable force fields increases the computational requirements compared to unpolarized approaches. The development of more accurate force fields will benefit from an extensive, open-source database comprising simulated ensembles and experimental data of various functional and toxic IDPs. To achieve such a platform, a crowd-sourcing approach is likely needed. This will also help tackling the “Big Data” problem, as with increasing computational power,120,121 emerging efficient sampling schemes,122 and wide interest in protein simulations, the cumulative amount of simulation data is becoming massive. How to merge, analyze, and build predictive models from this massive amount of data is becoming a research challenge.123 Application of machine learning and statistical methods will likely play a critical role in tackling this big data challenge. The open simulation data platform can also be helpful in developing novel coarse-grained models that are transferable (e.g., globular proteins to disordered proteins, in solution to on membrane) and more accurate. Such CG models can bias and guide all-atom simulations for further validation, establish complementary tools to experiments, and improve our current understanding of IDP behavior in cellular environments. For example, explicit representation of protein hydration water should be considered in CG models,98,124 as hydration water surrounding IDPs appears crucial for determining their structure and function.125,126 In addition, both atomistic and CG models should be calibrated as well as combined with various, often fragmentary, experimental data, for enabling deeper interpretation.

the local structural propensity and accounting for the nonnative and electrostatic interactions. Inclusion of those essential factors in the CG model has allowed dissecting the mechanism of coupled folding and binding of disordered proteins, which is highly valuable, as experimental characterization as well as atomistic simulation of IDP binding is nontrivial. One important observation is formation of a transient encounter complex at the initial stage of IDP binding with a partner, which is driven by non-native, intermolecular interactions, often electrostatic in nature.104−106 While CG models allow exploration of dimensions, time scales, and rare events that are inaccessible to atomistic simulations, they often lack transferability and fail to accurately reproduce experimental observables (structure, thermodynamics, kinetics, solvation effect), when compared to all-atom models. Further development, refinement, and optimization of CG models to make them more reliable and predictive in the context of IDP simulations is an active area of research. For example, a CG modeling framework at a residue-level resolution has been recently introduced to study the long time scale liquid−liquid phase separation process of IDPs.107



OUTLOOK The emerging synergy between simulations of various resolution and experiments has established simulations as a promising tool for delineating the sequence-structure-(mis)function relationship of disordered proteins. The structural transience and/or heterogeneity of soluble oligomers revealed in simulations further encourages use of computational modeling toward developing strategies for rational manipulation of IDP assembly. Nevertheless, the need for reassessment and reparametrization of existing all-atom protein force fields and water models is becoming increasingly evident, as those parameters were originally purposed for reproducing the properties of the folded, globular proteins. For example, two separate studies (refs 58 and 108) compared multiple modern force fields and reported comparable performance in characterizing Aβ conformations. All force fields, despite their fundamental differences, yielded similar tertiary structure sampling as well as adequate agreement with NMR data such as J-coupling. While convergence among different force fields at accurately representing the IDP conformational ensemble is encouraging, more work needs to be done to further fine-tune force fields or update solvent models109 for maintaining accurate secondary structure balance and chain dimension. The crucial role of hydration water and hydrodynamic interactions in IDP misfolding and aggregation is evident.110,111 One effort in this direction involves designing force field(s) optimized for right secondary structure balance (AMBER ff03*)112 to capturing temperature cooperativity in protein folding (AMBER03w)113 to reproducing unfolded protein dimensions consistent with experiment (AMBERws).114 The key difference between ff03* and ff03w is the use of the TIP4P/2005 water model in the latter, which is expected to provide a better representation of solvation effects such as temperature and pressure dependence. In ff0ws, protein−water interactions are strengthened by 10% from standard mixing rules to achieve the right balance of cross-interactions. The improvements in the resulting force field highlight the important role of capturing water-mediated interactions appropriately in a biomolecular force field. Recent development of the CHARMM36m115 force field is another notable effort in this direction, which involves refinement of the CHARMM36

The development of more accurate force fields will benefit from an extensive, open-source database comprising simulated ensembles and experimental data of various functional and toxic IDPs. In order to obtain the full mechanistic picture of IDP function, the peptides need to be simulated in a crowded, celllike environment. In addition, the short time scale structural dynamics should be combined with the assembly/binding E

DOI: 10.1021/acscentsci.7b00626 ACS Cent. Sci. XXXX, XXX, XXX−XXX

ACS Central Science

Outlook

(5) Wright, P. E.; Dyson, H. J. Intrinsically Disordered Proteins in Cellular Signaling and Regulation. Nat. Rev. Mol. Cell Biol. 2015, 16 (1), 18−29. (6) Dobson, C. M. Protein misfolding, evolution and disease. Trends Biochem. Sci. 1999, 24 (9), 329−332. (7) Knowles, T. P.; Vendruscolo, M.; Dobson, C. M. The amyloid state and its association with protein misfolding diseases. Nat. Rev. Mol. Cell Biol. 2014, 15 (6), 384−396. (8) Stefani, M.; Dobson, C. M. Protein aggregation and aggregate toxicity: new insights into protein folding, misfolding diseases and biological evolution. J. Mol. Med. (Heidelberg, Ger.) 2003, 81 (11), 678−699. (9) Haass, C.; Selkoe, D. J. Soluble protein oligomers in neurodegeneration: lessons from the Alzheimer’s amyloid betapeptide. Nat. Rev. Mol. Cell Biol. 2007, 8 (2), 101−112. (10) Campioni, S.; Mannini, B.; Zampagni, M.; Pensalfini, A.; Parrini, C.; Evangelisti, E.; Relini, A.; Stefani, M.; Dobson, C. M.; Cecchi, C.; Chiti, F. A causative link between the structure of aberrant protein oligomers and their toxicity. Nat. Chem. Biol. 2010, 6 (2), 140−147. (11) Winner, B.; Jappelli, R.; Maji, S. K.; Desplats, P. A.; Boyer, L.; Aigner, S.; Hetzer, C.; Loher, T.; Vilar, M.; Campioni, S.; Tzitzilonis, C.; Soragni, A.; Jessberger, S.; Mira, H.; Consiglio, A.; Pham, E.; Masliah, E.; Gage, F. H.; Riek, R. In vivo demonstration that alphasynuclein oligomers are toxic. Proc. Natl. Acad. Sci. U. S. A. 2011, 108 (10), 4194−4199. (12) Takahashi, T.; Kikuchi, S.; Katada, S.; Nagai, Y.; Nishizawa, M.; Onodera, O. Soluble polyglutamine oligomers formed prior to inclusion body formation are cytotoxic. Hum. Mol. Genet. 2008, 17 (3), 345−356. (13) Cheng, I. H.; Scearce-Levie, K.; Legleiter, J.; Palop, J. J.; Gerstein, H.; Bien-Ly, N.; Puolivali, J.; Lesne, S.; Ashe, K. H.; Muchowski, P. J.; Mucke, L. Accelerating amyloid-beta fibrillization reduces oligomer levels and functional deficits in Alzheimer disease mouse models. J. Biol. Chem. 2007, 282 (33), 23818−23828. (14) Larson, M. E.; Lesne, S. E. Soluble Abeta oligomer production and toxicity. J. Neurochem. 2012, 120, 125−139. (15) Chen, B.; Retzlaff, M.; Roos, T.; Frydman, J. Cellular strategies of protein quality control. Cold Spring Harbor Perspect. Biol. 2011, 3 (8), a004374. (16) 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 (39), 13692−13700. (17) Schuler, B.; Soranno, A.; Hofmann, H.; Nettels, D. Singlemolecule FRET spectroscopy and the polymer physics of unfolded and intrinsically disordered proteins. Annu. Rev. Biophys. 2016, 45, 207− 231. (18) Salvi, N.; Abyzov, A.; Blackledge, M. Atomic resolution conformational dynamics of intrinsically disordered proteins from NMR spin relaxation. Prog. Nucl. Magn. Reson. Spectrosc. 2017, 102103, 43−60. (19) Feldman, H. J.; Hogue, C. W. V. Probabilistic sampling of protein conformations: New hope for brute force? Proteins: Struct., Funct., Genet. 2002, 46 (1), 8−23. (20) Ozenne, V.; Bauer, F.; Salmon, L.; Huang, J.-r.; Jensen, M. R.; Segard, S.; Bernadó, P.; Charavay, C.; Blackledge, M. Flexiblemeccano: a tool for the generation of explicit ensemble descriptions of intrinsically disordered proteins and their associated experimental observables. Bioinformatics 2012, 28 (11), 1463−1470. (21) Krzeminski, M.; Marsh, J. A.; Neale, C.; Choy, W.-Y.; FormanKay, J. D. Characterization of disordered proteins with ENSEMBLE. Bioinformatics 2013, 29 (3), 398−399. (22) Salmon, L.; Nodet, G.; Ozenne, V.; Yin, G.; Jensen, M. R.; Zweckstetter, M.; Blackledge, M. NMR Characterization of LongRange Order in Intrinsically Disordered Proteins. J. Am. Chem. Soc. 2010, 132 (24), 8407−8418. (23) Fisher, C. K.; Stultz, C. M. Constructing Ensembles for Intrinsically Disordered Proteins. Curr. Opin. Struct. Biol. 2011, 21 (3), 426−431.

occurring at the larger spatiotemporal scale (Figure 4). As shown in Figure 4, such investigations may shed light onto the

Figure 4. A schematic diagram depicting potential states populated on the IDP assembly landscape.

relation between functional droplets and amyloid-related toxicity. While the computational efficiency of coarse-grained models makes them an attractive choice for this purpose, they need to be carefully combined with higher resolution models, highlighting the need for constructing novel multiscale modeling strategies. One way to achieve that goal is by allowing smooth, on-the-fly exchange of particles between different resolutions.127,128



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Payel Das: 0000-0002-7288-0516 Silvina Matysiak: 0000-0003-3824-9787 Jeetain Mittal: 0000-0002-9725-6402 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors are grateful to the many collaborators who are listed in relevant cited references and to the anonymous reviewers. The work mentioned in this Outlook was supported as part of the IBM Blue Gene Science Program (P.D.), NIGMS R01GM118530 (J.M.), and National Science Foundation under the Grant CHE-1454948 (S.M.). The simulations were performed using in-house IBM Blue Gene Supercomputers (P.D.), XSEDE facility (NSF Project No. TG-MCB120014), and TG-MCB120045 (S.M.). The authors thank Srirupa Chakraborty for assistance with making Figure 4.



REFERENCES

(1) Uversky, V. N.; Gillespie, J. R.; Fink, A. L. Why are ″natively unfolded″ proteins unstructured under physiologic conditions? Proteins: Struct., Funct., Genet. 2000, 41 (3), 415−427. (2) Gall, T. L.; Romero, P. R.; Cortese, M. S.; Uversky, V. N.; Dunker, A. K. Intrinsic disorder in the protein data bank. J. Biomol. Struct. Dyn. 2007, 24 (4), 325−341. (3) Marsh, J. A.; Forman-Kay, J. D. Sequence determinants of compaction in intrinsically disordered proteins. Biophys. J. 2010, 98 (10), 2383−2390. (4) Tompa, P.; Fuxreiter, M. Fuzzy complexes: polymorphism and structural disorder in protein-protein interactions. Trends Biochem. Sci. 2008, 33 (1), 2−8. F

DOI: 10.1021/acscentsci.7b00626 ACS Cent. Sci. XXXX, XXX, XXX−XXX

ACS Central Science

Outlook

systems: Application to protein structure prediction. J. Chem. Phys. 2002, 116 (2C0), 9058−9067. (44) Sugita, Y.; Okamoto, Y. Replica-exchange molecular dynamics method for protein folding. Chem. Phys. Lett. 1999, 314 (1-2), 141− 151. (45) Barducci, A.; Bonomi, M.; Parrinello, M. Metadynamics. Wiley Interdisciplinary Reviews: Computational Molecular Science 2011, 1 (5), 826−843. (46) Zerze, C. l. H.; Miller, C. M.; Granata, D.; Mittal, J. Free energy surface of an intrinsically disordered protein: comparison between temperature replica exchange molecular dynamics and bias-exchange metadynamics. J. Chem. Theory Comput. 2015, 11 (6), 2776−2782. (47) Liu, P.; Kim, B.; Friesner, R. A.; Berne, B. Replica exchange with solute tempering: A method for sampling biological systems in explicit water. Proc. Natl. Acad. Sci. U. S. A. 2005, 102 (39), 13749−13754. (48) Wang, L.; Friesner, R. A.; Berne, B. Replica exchange with solute scaling: a more efficient version of replica exchange with solute tempering (REST2). J. Phys. Chem. B 2011, 115 (30), 9431−9438. (49) Miller, C. M.; Brown, A. C.; Mittal, J. Disorder in cholesterolbinding functionality of CRAC peptides: A molecular dynamics study. J. Phys. Chem. B 2014, 118 (46), 13169−13174. (50) Bonomi, M.; Parrinello, M. Enhanced sampling in the welltempered ensemble. Phys. Rev. Lett. 2010, 104 (19), 190601. (51) Deighan, M.; Bonomi, M.; Pfaendtner, J. Efficient simulation of explicitly solvated proteins in the well-tempered ensemble. J. Chem. Theory Comput. 2012, 8 (7), 2189−2192. (52) Conicella, A. E.; Zerze, G. H.; Mittal, J.; Fawzi, N. L. ALS mutations disrupt phase separation mediated by -helical structure in the TDP-43 low complexity C-terminal domain. Structure 2016, 24 (9), 1537−1549. (53) Miller, C.; Zerze, G. l. H.; Mittal, J. Molecular simulations indicate marked differences in the structure of amylin mutants, correlated with known aggregation propensity. J. Phys. Chem. B 2013, 117 (50), 16066−16075. (54) Das, P.; Murray, B.; Belfort, G. Alzheimer’s protective A2T mutation changes the conformational landscape of the Abeta1−42 monomer differently than does the A2V mutation. Biophys. J. 2015, 108 (3), 738−747. (55) Di Fede, G.; Catania, M.; Morbin, M.; Rossi, G.; Suardi, S.; Mazzoleni, G.; Merlin, M.; Giovagnoli, A. R.; Prioni, S.; Erbetta, A.; Falcone, C.; Gobbi, M.; Colombo, L.; Bastone, A.; Beeg, M.; Manzoni, C.; Francescucci, B.; Spagnoli, A.; Cantu, L.; Del Favero, E.; Levy, E.; Salmona, M.; Tagliavini, F. A Recessive Mutation in the APP Gene with Dominant-Negative Effect on Amyloidogenesis. Science 2009, 323 (5920), 1473−1477. (56) Jonsson, T.; Atwal, J. K.; Steinberg, S.; Snaedal, J.; Jonsson, P. V.; Bjornsson, S.; Stefansson, H.; Sulem, P.; Gudbjartsson, D.; Maloney, J. A mutation in APP protects against Alzheimer/’s disease and age-related cognitive decline. Nature 2012, 488 (7409), 96−99. (57) Nguyen, P. H.; Tarus, B.; Derreumaux, P. Familial Alzheimer A2V mutation reduces the intrinsic disorder and completely changes the free energy landscape of the Aβ1−28 monomer. J. Phys. Chem. B 2014, 118 (2), 501−510. (58) Rosenman, D. J.; Wang, C.; García, A. E. Characterization of Aβ monomers through the convergence of ensemble properties among simulations with multiple force fields. J. Phys. Chem. B 2016, 120 (2), 259−277. (59) Sgourakis, N. G.; Merced-Serrano, M.; Boutsidis, C.; Drineas, P.; Du, Z.; Wang, C.; Garcia, A. E. Atomic-level characterization of the ensemble of the Aβ(1−42) monomer in water using unbiased Molecular Dynamics simulations and spectral algorithms. J. Mol. Biol. 2011, 405 (2), 570−583. (60) Rosenman, D. J.; Connors, C. R.; Chen, W.; Wang, C.; Garcia, A. E. Abeta Monomers Transiently Sample Oligomer and Fibril-Like Configurations: Ensemble Characterization Using a Combined MD/ NMR Approach. J. Mol. Biol. 2013, 425 (18), 3338−3359. (61) Ball, K. A.; Phillips, A. H.; Wemmer, D. E.; Head-Gordon, T. Differences in beta-strand Populations of Monomeric A beta 40 and A beta 42. Biophys. J. 2013, 104 (12), 2714−2724.

(24) Cavalli, A.; Camilloni, C.; Vendruscolo, M. Molecular dynamics simulations with replica-averaged structural restraints generate structural ensembles according to the maximum entropy principle. J. Chem. Phys. 2013, 138 (9), 094112. (25) Ball, K. A.; Wemmer, D. E.; Head-Gordon, T. Comparison of structure determination methods for intrinsically disordered amyloid-β peptides. J. Phys. Chem. B 2014, 118 (24), 6405−6416. (26) Oostenbrink, C.; Villa, A.; Mark, A. E.; Van Gunsteren, W. F. A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force-field parameter sets 53A5 and 53A6. J. Comput. Chem. 2004, 25 (13), 1656−1676. (27) Jorgensen, W. L.; Maxwell, D.; Tirado-Rives, J. Development and testing of the OPLS all-atom force field on confromational energetics and properties of organic liquids. J. Am. Chem. Soc. 1996, 118, 11225−11236. (28) MacKerell, A. D., Jr; Bashford, D.; Bellott, M.; Dunbrack, R. L., Jr; Evanseck, J. D.; Field, M. J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S. Allatom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B 1998, 102 (18), 3586−3616. (29) Cornell, W. D.; Cieplak, P.; Bayly, C. I.; Gould, I. R.; Merz, K. M.; Ferguson, D. M.; Spellmeyer, D. C.; Fox, T.; Caldwell, J. W.; Kollman, P. A. A second generation force field for the simulation of proteins, nucleic acids, and organic molecules. J. Am. Chem. Soc. 1995, 117 (19), 5179−5197. (30) Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926−935. (31) Berendsen, H. J.; Postma, J. P.; van Gunsteren, W. F.; Hermans, J. Interaction models for water in relation to protein hydration. In Intermolecular forces; Springer: 1981; pp 331−342, DOI: 10.1007/97894-015-7658-1_21. (32) Berendsen, H.; Grigera, J.; Straatsma, T. The missing term in effective pair potentials. J. Phys. Chem. 1987, 91 (24), 6269−6271. (33) Stanley, N.; Esteban-Martín, S.; De Fabritiis, G. Progress in studying intrinsically disordered proteins with atomistic simulations. Prog. Biophys. Mol. Biol. 2015, 119 (1), 47−52. (34) Nguyen, P. H.; Li, M. S.; Stock, G.; Straub, J. E.; Thirumalai, D. Monomer adds to preformed structured oligomers of Aβ-peptides by a two-stage dock−lock mechanism. Proc. Natl. Acad. Sci. U. S. A. 2007, 104 (1), 111−116. (35) Krone, M. G.; Hua, L.; Soto, P.; Zhou, R.; Berne, B.; Shea, J.-E. Role of water in mediating the assembly of Alzheimer amyloid-β Aβ16− 22 protofilaments. J. Am. Chem. Soc. 2008, 130 (33), 11066− 11072. (36) Zheng, J.; Ma, B.; Tsai, C.-J.; Nussinov, R. Structural stability and dynamics of an amyloid-forming peptide GNNQQNY from the yeast prion sup-35. Biophys. J. 2006, 91 (3), 824−833. (37) Mo, Y.; Lei, J.; Sun, Y.; Zhang, Q.; Wei, G. Conformational Ensemble of hIAPP Dimer: Insight into the Molecular Mechanism by which a Green Tea Extract inhibits hIAPP Aggregation. Sci. Rep. 2016, 6, 33076. (38) Zhang, T.; Zhang, J.; Derreumaux, P.; Mu, Y. Molecular Mechanism of the Inhibition of EGCG on the Alzheimer Aβ1−42 Dimer. J. Phys. Chem. B 2013, 117 (15), 3993−4002. (39) Zheng, X.; Wu, C.; Liu, D.; Li, H.; Bitan, G.; Shea, J. E.; Bowers, M. T. Mechanism of C-Terminal Fragments of Amyloid beta-Protein as Abeta Inhibitors: Do C-Terminal Interactions Play a Key Role in Their Inhibitory Activity? J. Phys. Chem. B 2016, 120 (8), 1615−1623. (40) Brown, A. M; Bevan, D. R. Molecular Dynamics Simulations of Amyloid β-Peptide (1−42): Tetramer Formation and Membrane Interactions. Biophys. J. 2016, 111 (5), 937−949. (41) Nagarajan, A.; Jawahery, S.; Matysiak, S. The effects of flanking sequences in the interaction of polyglutamine peptides with a membrane bilayer. J. Phys. Chem. B 2014, 118 (24), 6368−6379. (42) Das, P.; Kang, S.-g.; Temple, S.; Belfort, G. Interaction of Amyloid Inhibitor Proteins with Amyloid Beta Peptides: Insight from Molecular Dynamics Simulations. PLoS One 2014, 9 (11), e113041. (43) Fukunishi, H.; Watanabe, O.; Takada, S. On the Hamiltonian replica exchange method for efficient sampling of biomolecular G

DOI: 10.1021/acscentsci.7b00626 ACS Cent. Sci. XXXX, XXX, XXX−XXX

ACS Central Science

Outlook

(62) Hoyer, W.; Grönwall, C.; Jonsson, A.; Ståhl, S.; Härd, T. Stabilization of beta-hairpin in monomeric Alzheimer’s amyloid-beta peptide inhibits amyloid formation. Proc. Natl. Acad. Sci. U. S. A. 2008, 105 (13), 5099−5104. (63) Ahmed, M.; Davis, J.; Aucoin, D.; Sato, T.; Ahuja, S.; Aimoto, S.; Elliott, J. I.; Van Nostrand, W. E.; Smith, S. O. Structural conversion of neurotoxic amyloid-[beta]1−42 oligomers to fibrils. Nat. Struct. Mol. Biol. 2010, 17 (5), 561−567. (64) Man, V. H.; Nguyen, P. H.; Derreumaux, P. Conformational Ensembles of the Wild-Type and S8C Aβ1−42 Dimers. J. Phys. Chem. B 2017, 121 (11), 2434−2442. (65) Man, V. H.; Nguyen, P. H.; Derreumaux, P. High-Resolution Structures of the Amyloid-β 1−42 Dimers from the Comparison of Four Atomistic Force Fields. J. Phys. Chem. B 2017, 121 (24), 5977− 5987. (66) Tarus, B.; Tran, T. T.; Nasica-Labouze, J.; Sterpone, F.; Nguyen, P. H.; Derreumaux, P. Structures of the Alzheimer’s Wild-Type Aβ1− 40 Dimer from Atomistic Simulations. J. Phys. Chem. B 2015, 119 (33), 10478−10487. (67) Nguyen, P. H.; Sterpone, F.; Campanera, J. M.; Nasica-Labouze, J.; Derreumaux, P. Impact of the A2V mutation on the Heterozygous and Homozygous Aβ1−40 Dimer Structures from Atomistic Simulations. ACS Chem. Neurosci. 2016, 7 (6), 823−832. (68) Nguyen, P. H.; Sterpone, F.; Pouplana, R.; Derreumaux, P.; Campanera, J. M. Dimerization mechanism of Alzheimer Aβ40 peptides: the high content of intrapeptide-stabilized conformations in A2V and A2T heterozygous dimers retards amyloid fibril formation. J. Phys. Chem. B 2016, 120 (47), 12111−12126. (69) Das, P.; Chacko, A. R.; Belfort, G. Alzheimer’s Protective CrossInteraction between Wild-Type and A2T Variants Alters Aβ42 Dimer Structure. ACS Chem. Neurosci. 2017, 8 (6), 606−618. (70) Larini, L.; Shea, J.-E. Role of beta-Hairpin Formation in Aggregation: The Self-Assembly of the Amyloid-beta(25−35) Peptide. Biophys. J. 2012, 103 (3), 576−586. (71) Reddy, G.; Straub, J. E.; Thirumalai, D. Influence of Preformed Asp23−Lys28 Salt Bridge on the Conformational Fluctuations of Monomers and Dimers of Aβ Peptides with Implications for Rates of Fibril Formation. J. Phys. Chem. B 2009, 113 (4), 1162−1172. (72) Sandberg, A.; Luheshi, L. M.; Sollvander, S.; Pereira de Barros, T.; Macao, B.; Knowles, T. P. J.; Biverstal, H.; Lendel, C.; EkholmPetterson, F.; Dubnovitsky, A.; Lannfelt, L.; Dobson, C. M.; Hard, T. Stabilization of neurotoxic Alzheimer amyloid-beta oligomers by protein engineering. Proc. Natl. Acad. Sci. U. S. A. 2010, 107 (35), 15595−15600. (73) Murray, B.; Sorci, M.; Rosenthal, J.; Lippens, J.; Isaacson, D.; Das, P.; Fabris, D.; Li, S.; Belfort, G. A2T and A2V Abeta peptides exhibit different aggregation kinetics, primary nucleation, morphology, structure, and LTP inhibition. Proteins: Struct., Funct., Genet. 2016, 84 (4), 488−500. (74) Benilova, I.; Gallardo, R.; Ungureanu, A.-A.; Castillo Cano, V.; Snellinx, A.; Ramakers, M.; Bartic, C.; Rousseau, F.; Schymkowitz, J.; De Strooper, B. The Alzheimer Disease Protective Mutation Ala2Thr Modulates Kinetic and Thermodynamic Properties of Abeta Aggregation. J. Biol. Chem. 2014, 289 (45), 30977−30989. (75) Zheng, X.; Liu, D.; Roychaudhuri, R.; Teplow, D. B.; Bowers, M. T. Amyloid β-Protein Assembly: Differential Effects of the Protective A2T Mutation and Recessive A2V Familial Alzheimer’s Disease Mutation. ACS Chem. Neurosci. 2015, 6 (10), 1732−1740. (76) Chakraborty, S.; Das, P. Emergence of Alternative Structures in Amyloid Beta 1−42 Monomeric Landscape by N-terminal Hexapeptide Amyloid Inhibitors. Sci. Rep. 2017, 7 (1), 9941. (77) Ehrnhoefer, D. E.; Bieschke, J.; Boeddrich, A.; Herbst, M.; Masino, L.; Lurz, R.; Engemann, S.; Pastore, A.; Wanker, E. E. EGCG redirects amyloidogenic polypeptides into unstructured, off-pathway oligomers. Nat. Struct. Mol. Biol. 2008, 15 (6), 558−566. (78) Frid, P.; Anisimov, S. V.; Popovic, N. Congo red and protein aggregation in neurodegenerative diseases. Brain Res. Rev. 2007, 53 (1), 135−160.

(79) Ladiwala, A. R.; Lin, J. C.; Bale, S. S.; Marcelino-Cruz, A. M.; Bhattacharya, M.; Dordick, J. S.; Tessier, P. M. Resveratrol selectively remodels soluble oligomers and fibrils of amyloid Abeta into offpathway conformers. J. Biol. Chem. 2010, 285 (31), 24228−24237. (80) Tjernberg, L. O.; Naslund, J.; Lindqvist, F.; Johansson, J.; Karlstrom, A. R.; Thyberg, J.; Terenius, L.; Nordstedt, C. Arrest of beta-amyloid fibril formation by a pentapeptide ligand. J. Biol. Chem. 1996, 271 (15), 8545−8548. (81) Li, H.; Rahimi, F.; Bitan, G. Modulation of Amyloid betaProtein (Abeta) Assembly by Homologous C-Terminal Fragments as a Strategy for Inhibiting Abeta Toxicity. ACS Chem. Neurosci. 2016, 7 (7), 845−856. (82) Valera, E.; Spencer, B.; Masliah, E. Immunotherapeutic approaches targeting amyloid-β, α-synuclein, and tau for the treatment of neurodegenerative disorders. Neurotherapeutics 2016, 13 (1), 179− 189. (83) Doig, A. J.; del Castillo-Frias, M. P.; Berthoumieu, O.; Tarus, B.; Nasica-Labouze, J.; Sterpone, F.; Nguyen, P. H.; Hooper, N. M.; Faller, P.; Derreumaux, P. Why Is Research on Amyloid-β Failing to Give New Drugs for Alzheimer’s Disease? ACS Chem. Neurosci. 2017, 8, 1435−1437. (84) Di Fede, G.; Catania, M.; Maderna, E.; Morbin, M.; Moda, F.; Colombo, L.; Rossi, A.; Cagnotto, A.; Virgilio, T.; Palamara, L. Tackling amyloidogenesis in Alzheimer’s disease with A2V variants of Amyloid-β. Sci. Rep. 2016, 6, 20949. (85) Lin, T. W.; Chang, C. F.; Chang, Y. J.; Liao, Y. H.; Yu, H. M.; Chen, Y. R. Alzheimer’s amyloid-beta A2T variant and its N-terminal peptides inhibit amyloid-beta fibrillization and rescue the induced cytotoxicity. PLoS One 2017, 12 (3), e0174561. (86) Forest, K. H.; Alfulaij, N.; Arora, K.; Taketa, R.; Sherrin, T.; Todorovic, C.; Lawrence, J. L.; Yoshikawa, G. T.; Ng, H. L.; Hruby, V. J. Protection against β-amyloid neurotoxicity by a non-toxic endogenous N-terminal β-amyloid fragment and its active hexapeptide core sequence. J. Neurochem. 2018, 144 (2), 201−217. (87) Ait-Bouziad, N.; Lv, G.; Mahul-Mellier, A.-L.; Xiao, S.; Zorludemir, G.; Eliezer, D.; Walz, T.; Lashuel, H. A. Discovery and characterization of stable and toxic Tau/phospholipid oligomeric complexes. Nat. Commun. 2017, 8 (1), 1678. (88) Seidler, P.; Boyer, D.; Rodriguez, J.; Sawaya, M.; Cascio, D.; Murray, K.; Gonen, T.; Eisenberg, D. Structure-based inhibitors of tau aggregation. Nat. Chem. 2018, 10 (2), 170−176. (89) Bottaro, S.; Lindorff-Larsen, K.; Best, R. B. Variational optimization of an all-atom implicit solvent force field to match explicit solvent simulation data. J. Chem. Theory Comput. 2013, 9 (12), 5641−5652. (90) Lazaridis, T.; Karplus, M. Effective energy function for proteins in solution. Proteins: Struct., Funct., Genet. 1999, 35 (2), 133−152. (91) Vitalis, A.; Pappu, R. V. ABSINTH: A new continuum solvation model for simulations of polypeptides in aqueous solutions. J. Comput. Chem. 2009, 30 (5), 673−699. (92) Mao, A. H.; Lyle, N.; Pappu, R. V. Describing sequenceensemble relationships for intrinsically disordered proteins. Biochem. J. 2013, 449 (2), 307−318. (93) Das, R. K.; Pappu, R. V. Conformations of intrinsically disordered proteins are influenced by linear sequence distributions of oppositely charged residues. Proc. Natl. Acad. Sci. U. S. A. 2013, 110 (33), 13392−13397. (94) Sherry, K. P.; Das, R. K.; Pappu, R. V.; Barrick, D. Control of transcriptional activity by design of charge patterning in the intrinsically disordered RAM region of the Notch receptor. Proc. Natl. Acad. Sci. U. S. A. 2017, 114 (44), E9243−E9252. (95) Urbanc, B.; Cruz, L.; Ding, F.; Sammond, D.; Khare, S.; Buldyrev, S. V.; Stanley, H. E.; Dokholyan, N. V. Molecular Dynamics Simulation of Amyloid beta Dimer Formation. Biophys. J. 2004, 87 (4), 2310−2321. (96) Monticelli, L.; Kandasamy, S. K.; Periole, X.; Larson, R. G.; Tieleman, D. P.; Marrink, S. J. The MARTINI Coarse-Grained Force Field: Extension to Proteins. J. Chem. Theory Comput. 2008, 4 (5), 819−834. H

DOI: 10.1021/acscentsci.7b00626 ACS Cent. Sci. XXXX, XXX, XXX−XXX

ACS Central Science

Outlook

(97) Cheon, M.; Chang, I.; Hall, C. K. Extending the PRIME model for protein aggregation to all 20 amino acids. Proteins: Struct., Funct., Genet. 2010, 78 (14), 2950−2960. (98) Ganesan, S. J.; Matysiak, S. Interplay between the hydrophobic effect and dipole interactions in peptide aggregation at interfaces. Phys. Chem. Chem. Phys. 2016, 18 (4), 2449−2458. (99) Ganesan, S. J.; Xu, H.; Matysiak, S. Effect of lipid head group interactions on membrane properties and membrane-induced cationic β-hairpin folding. Phys. Chem. Chem. Phys. 2016, 18 (27), 17836− 17850. (100) Sterpone, F.; Melchionna, S.; Tuffery, P.; Pasquali, S.; Mousseau, N.; Cragnolini, T.; Chebaro, Y.; St-Pierre, J.-F.; Kalimeri, M.; Barducci, A. The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems. Chem. Soc. Rev. 2014, 43 (13), 4871−4893. (101) Davtyan, A.; Schafer, N. P.; Zheng, W.; Clementi, C.; Wolynes, P. G.; Papoian, G. A. AWSEM-MD: protein structure prediction using coarse-grained physical potentials and bioinformatically based local structure biasing. J. Phys. Chem. B 2012, 116 (29), 8494−8503. (102) Chen, M.; Tsai, M.; Zheng, W.; Wolynes, P. G. The Aggregation Free Energy Landscapes of Polyglutamine Repeats. J. Am. Chem. Soc. 2016, 138 (46), 15197−15203. (103) Clementi, C.; Nymeyer, H.; Onuchic, J. N. Topological and energetic factors: what determines the structural details of the transition state ensemble and “en-route” intermediates for protein folding? An investigation for small globular proteins. J. Mol. Biol. 2000, 298 (5), 937−953. (104) Chu, W.-T.; Clarke, J.; Shammas, S. L.; Wang, J. Role of nonnative electrostatic interactions in the coupled folding and binding of PUMA with Mcl-1. PLoS Comput. Biol. 2017, 13 (4), e1005468. (105) Ganguly, D.; Zhang, W.; Chen, J. Electrostatically Accelerated Encounter and Folding for Facile Recognition of Intrinsically Disordered Proteins. PLoS Comput. Biol. 2013, 9 (11), e1003363. (106) Turjanski, A. G.; Gutkind, J. S.; Best, R. B.; Hummer, G. Binding-induced folding of a natively unstructured transcription factor. PLoS Comput. Biol. 2008, 4 (4), e1000060. (107) Dignon, G. L.; Zheng, W.; Kim, Y. C.; Best, R. B.; Mittal, J. Sequence determinants of protein phase behavior from a coarsegrained model. PLoS Comput. Biol. 2018, 14 (1), e1005941. (108) Carballo-Pacheco, M.; Strodel, B. Comparison of force fields for Alzheimer’s A beta42: A case study for intrinsically disordered proteins. Protein Sci. 2017, 26 (2), 174−185. (109) 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 (16), 5113−5123. (110) Chong, S.-H.; Ham, S. Distinct role of hydration water in protein misfolding and aggregation revealed by fluctuating thermodynamics analysis. Acc. Chem. Res. 2015, 48 (4), 956−965. (111) Chiricotto, M.; Melchionna, S.; Derreumaux, P.; Sterpone, F. Hydrodynamic effects on β-amyloid (16−22) peptide aggregation. J. Chem. Phys. 2016, 145 (3), 035102. (112) Best, R. B.; Mittal, J. Protein simulations with an optimized water model: cooperative helix formation and temperature-induced unfolded state collapse. J. Phys. Chem. B 2010, 114 (46), 14916− 14923. (113) Best, R. B.; Hummer, G. Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides. J. Phys. Chem. B 2009, 113 (26), 9004−9015. (114) Best, R. B.; Zheng, W.; Mittal, J. Balanced protein−water interactions improve properties of disordered proteins and nonspecific protein association. J. Chem. Theory Comput. 2014, 10 (11), 5113−5124. (115) Huang, J.; Rauscher, S.; Nawrocki, G.; Ran, T.; Feig, M.; de Groot, B. L.; Grubmüller, H.; MacKerell, A. D., Jr CHARMM36m: an improved force field for folded and intrinsically disordered proteins. Nat. Methods 2017, 14 (1), 71−73.

(116) Faller, P.; Hureau, C.; La Penna, G. Metal ions and intrinsically disordered proteins and peptides: from Cu/Zn amyloid-β to general principles. Acc. Chem. Res. 2014, 47 (8), 2252−2259. (117) Tahmasebinia, F.; Emadi, S. Effect of metal chelators on the aggregation of beta-amyloid peptides in the presence of copper and iron. BioMetals 2017, 30 (2), 285−293. (118) Li, P.; Merz, K. M., Jr Metal ion modeling using classical mechanics. Chem. Rev. 2017, 117 (3), 1564−1686. (119) Jeziorski, B.; Moszynski, R.; Szalewicz, K. Perturbation theory approach to intermolecular potential energy surfaces of van der Waals complexes. Chem. Rev. 1994, 94 (7), 1887−1930. (120) Shaw, D. E.; Grossman, J.; Bank, J. A.; Batson, B.; Butts, J. A.; Chao, J. C.; Deneroff, M. M.; Dror, R. O.; Even, A.; Fenton, C. H. Anton 2: raising the bar for performance and programmability in a special-purpose molecular dynamics supercomputer. In Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis; IEEE Press: 2014; pp 41−53, DOI: 10.1109/SC.2014.9. (121) Stone, J. E.; Hardy, D. J.; Ufimtsev, I. S.; Schulten, K. GPUaccelerated molecular modeling coming of age. J. Mol. Graphics Modell. 2010, 29 (2), 116−125. (122) Pan, A. C.; Weinreich, T. M.; Piana, S.; Shaw, D. E. Demonstrating an order-of-magnitude sampling enhancement in molecular dynamics simulations of complex protein systems. J. Chem. Theory Comput. 2016, 12 (3), 1360−1367. (123) Lane, T. J.; Shukla, D.; Beauchamp, K. A.; Pande, V. S. To milliseconds and beyond: challenges in the simulation of protein folding. Curr. Opin. Struct. Biol. 2013, 23 (1), 58−65. (124) Custer, G. S.; Das, P.; Matysiak, S. Interplay between Conformational Heterogeneity and Hydration in the Folding Landscape of a Designed Three-Helix Bundle. J. Phys. Chem. B 2017, 121 (13), 2731−2738. (125) Fichou, Y.; Schirò, G.; Gallat, F.-X.; Laguri, C.; Moulin, M.; Combet, J.; Zamponi, M.; Härtlein, M.; Picart, C.; Mossou, E.; LortatJacob, H.; Colletier, J.-P.; Tobias, D. J.; Weik, M. Hydration water mobility is enhanced around tau amyloid fibers. Proc. Natl. Acad. Sci. U. S. A. 2015, 112 (20), 6365−6370. (126) Dalal, V.; Arya, S.; Mukhopadhyay, S. Confined Water in Amyloid-Competent Oligomers of the Prion Protein. ChemPhysChem 2016, 17 (18), 2804−2807. (127) Praprotnik, M.; Delle Site, L.; Kremer, K. Adaptive resolution molecular-dynamics simulation: Changing the degrees of freedom on the fly. J. Chem. Phys. 2005, 123 (22), 224106. (128) Fiorentini, R.; Kremer, K.; Potestio, R.; Fogarty, A. C. Using force-based adaptive resolution simulations to calculate solvation free energies of amino acid sidechain analogues. J. Chem. Phys. 2017, 146 (24), 244113.

I

DOI: 10.1021/acscentsci.7b00626 ACS Cent. Sci. XXXX, XXX, XXX−XXX