Comparative analysis of essential collective dynamics and NMR-derived flexibility profiles in evolutionarily diverse prion proteins

Essential collective dynamics analysis reveals nonlocal interactions of alpha-synuclein38-95 monomers with fibrillar seeds

Alpha-synuclein (α-syn) is a cytoplasmic protein commonly found in the nervous system. In solution, α-syn adopts disordered unfolded conformations, although it can also form α-helices upon binding to lipid membranes. Under conditions that are not yet fully understood, α-syn can misfold and aggregate, giving rise to β-sheet rich amyloid fibrils, which then tend to accumulate in degenerating neurons. This leads to Parkinson's disease (PD) and several other conditions collectively termed synucleinopathies. Development of disease-modifying treatments requires detailed understanding of structure and dynamics of α-syn's misfolded aggregates. We have employed 1000 ns long, all-atom molecular dynamics simulations to investigate the interaction of monomeric α-syn_38-95 fragments, which contain the most important amyloidogenic regions, with preformed fibrillar seeds composed of staggered, β-sheet rich α-syn chains of matching length. The simulations indicate that α-syn_38-95 monomers tend to form aggregates with the fibrillar seeds, although we have not observed alignment of the monomeric chains with β-strands of the fibril. To analyze the stability of these aggregates, we have employed the essential collective dynamics method, which allows making accurate assessment of dynamical coupling across individual atoms in macromolecules and supramolecular complexes. The analysis revealed extensive dynamical coupling across initially monomeric α-syn chains and the fibrillar seeds including distal regions thereof that did not contact the monomer directly. We have discussed structural origins of these long-range interactions, their impacts for the stability of α-syn aggregates, and potential implications for the development of anti-PD treatments.

Aggregation of Aβ40/42 chains in the presence of cyclic neuropeptides investigated by molecular dynamics simulations

Alzheimer’s disease is associated with the formation of toxic aggregates of amyloid beta (Aβ) peptides. Despite tremendous efforts, our understanding of the molecular mechanisms of aggregation, as well as cofactors that might influence it, remains incomplete. The small cyclic neuropeptide somatostatin-14 (SST₁₄) was recently found to be the most selectively enriched protein in human frontal lobe extracts that binds Aβ₄₂ aggregates. Furthermore, SST₁₄’s presence was also found to promote the formation of toxic Aβ₄₂ oligomers in vitro. In order to elucidate how SST₁₄ influences the onset of Aβ oligomerization, we performed all-atom molecular dynamics simulations of model mixtures of Aβ₄₂ or Aβ₄₀ peptides with SST₁₄ molecules and analyzed the structure and dynamics of early-stage aggregates. For comparison we also analyzed the aggregation of Aβ₄₂ in the presence of arginine vasopressin (AVP), a different cyclic neuropeptide. We observed the formation of self-assembled aggregates containing the Aβ chains and small cyclic peptides in all mixtures of Aβ₄₂–SST₁₄, Aβ₄₂–AVP, and Aβ₄₀–SST₁₄. The Aβ₄₂–SST₁₄ mixtures were found to develop compact, dynamically stable, but small aggregates with the highest exposure of hydrophobic residues to the solvent. Differences in the morphology and dynamics of aggregates that comprise SST₁₄ or AVP appear to reflect distinct (1) regions of the Aβ chains they interact with; (2) propensities to engage in hydrogen bonds with Aβ peptides; and (3) solvent exposures of hydrophilic and hydrophobic groups. The presence of SST₁₄ was found to impede aggregation in the Aβ₄₂–SST₁₄ system despite a high hydrophobicity, producing a stronger “sticky surface” effect in the aggregates at the onset of Aβ₄₂–SST₁₄ oligomerization.

Improper folding of proteins causes disorders known as protein misfolding diseases. Under normal conditions most proteins adopt particular folds, which allow them functioning properly. However, for reasons that are not yet fully understood, proteins may misfold and aggregate, forming deposits known as amyloid fibrils, which accumulate in the brain or other tissues. This process affects functioning of the nervous system, gradually causing loss of cognitive abilities. Alzheimer’s disease is one of the most common diseases from this group. A better understanding of the aggregation of peptides implicated in Alzheimer’s disease, known as amyloid beta (Aβ) peptides, may facilitate the development of treatments that ameliorate or prevent the disease. We use detailed molecular dynamics simulations to investigate the influence of somatostatin-14 (SST₁₄), a cyclic neuropeptide that might be involved in the amyloidogenic aggregation of Aβ, on molecular processes occurring at the onset of Aβ aggregation. Results of these simulations explain how the presence of SST₁₄ might alter pathways of aggregation of Aβ, shedding light upon the possible role of extrinsic factors in the aggregation at a molecular level.

Combining molecular dynamics simulations and experimental analyses in protein misfolding

The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Detecting early stage structural changes in wild type, pathogenic and non-pathogenic prion variants using Markov state model

The conversion of prion protein from normal to scrapie followed by the aggregation and deposition of this scrapie form leads to various neurodegenerative diseases. A few studies carried out by researchers suggest that E219K prion mutant (glutamate to lysine mutation at residue position 219) is more stable than wild type protein. However a similar point mutation E200K (glutamate to lysine mutation at residue position 200) is pathogenic. In this study we have carried out detailed atomistic simulation of the wild type, pathogenic mutant E200K and E219K mutant which provides more stability. The aim of the study was to detect the early structural changes present in all the three variants which might be responsible for the stability or for their conversion from PrPC to PrPSc. MSM based analyses have been carried out to find out the differences between WT, E200K and E219K systems. Markov state model (MSM) analysis was able to predict the intermediate states which helped to understand the effect of same mutation at two different locations. The MSM analysis was able to show that the extra stability of E219K mutant may be a result of the increase in number of native contacts, strong salt bridges and less random motions. While pathogenicity of E200K mutant can be attributed to loss of some crucial salt-bridge interactions, increased random motions between helix 2 and helix 3.

Markov state model to find out the differences between WT, E200K and E219K systems.

A novel Gerstmann-Sträussler-Scheinker disease mutation defines a precursor for amyloidogenic 8 kDa PrP fragments and reveals N-terminal structural changes shared by other GSS alleles

To explore pathogenesis in a young Gerstmann-Sträussler-Scheinker Disease (GSS) patient, the corresponding mutation, an eight-residue duplication in the hydrophobic region (HR), was inserted into the wild type mouse PrP gene. Transgenic (Tg) mouse lines expressing this mutation (Tg.HRdup) developed spontaneous neurologic syndromes and brain extracts hastened disease in low-expressor Tg.HRdup mice, suggesting de novo formation of prions. While Tg.HRdup mice exhibited spongiform change, PrP aggregates and the anticipated GSS hallmark of a proteinase K (PK)-resistant 8 kDa fragment deriving from the center of PrP, the LGGLGGYV insertion also imparted alterations in PrP's unstructured N-terminus, resulting in a 16 kDa species following thermolysin exposure. This species comprises a plausible precursor to the 8 kDa PK-resistant fragment and its detection in adolescent Tg.HRdup mice suggests that an early start to accumulation could account for early disease of the index case. A 16 kDa thermolysin-resistant signature was also found in GSS patients with P102L, A117V, H187R and F198S alleles and has coordinates similar to GSS stop codon mutations. Our data suggest a novel shared pathway of GSS pathogenesis that is fundamentally distinct from that producing structural alterations in the C-terminus of PrP, as observed in other prion diseases such as Creutzfeldt-Jakob Disease and scrapie.

Prion diseases can be sporadic, infectious or genetic. The central event of all prion diseases is the structural conversion of the cellular prion protein (PrP^C) to its disease associated conformer, PrP^Sc. Gerstmann-Sträussler-Scheinker Disease (GSS) is a genetic prion disease presenting as a multi-systemic neurological syndrome. A novel mutation, an eight amino acid insertion, was discovered in a young GSS patient. We created transgenic mice expressing this mutation and found that they recapitulate key features of the disease; namely PrP deposition in the brain and a low molecular weight proteinase K (PK) resistant internal PrP fragment. While structural investigations did not reveal a gross alteration in the conformation of this mutant PrP, the insertion lying at the boundary of the globular domain causes alterations in the unstructured amino terminal portion of the protein such that it becomes resistant to digestion by the enzyme thermolysin. We demonstrate by kinetic analysis and sequential digestion that this novel thermolysin resistant species is a precursor to the pathognomonic PK resistant fragment. Analysis of samples from other GSS patients revealed this same signature, suggesting a common molecular pathway.

Probing oligomerization of amyloid beta peptide in silico

Aggregation of amyloid β (Aβ) peptide is implicated in fatal Alzheimer's disease, for which no cure is available. Understanding the mechanisms responsible for this aggregation is required in order for therapies to be developed. In an effort to better understand the molecular mechanisms involved in spontaneous aggregation of Aβ peptide, extensive molecular dynamics simulations are reported, and the results are analyzed through a combination of structural biology tools and a novel essential collective dynamics method. Several model systems composed of ten or twelve Aβ_17-42 chains in water are investigated, and the influence of metal ions is probed. The results suggest that Aβ monomers tend to aggregate into stable globular-like oligomers with 13-23% of β-sheet content. Two stages of oligomer formation have been identified: quick collapse within the first 40 ns of the simulation, characterized by a decrease in inter-chain separation and build-up of β-sheets, and the subsequent slow relaxation of the oligomer structure. The resulting oligomers comprise a stable, coherently moving sub-aggregate of 6-9 strongly inter-correlated chains. Cu²⁺ and Fe²⁺ ions have been found to develop coordination bonds with carboxylate groups of E22, D23 and A42, which remain stable during 200 ns simulations. The presence of Fe²⁺, and particularly Cu²⁺ ions, in negatively charged cavities has been found to cause significant changes in the structure and dynamics of the oligomers. The results indicate, in particular, that formation of non-fibrillar oligomers might be involved in early template-free aggregation of Aβ_17-42 monomers, with charged species such as Cu²⁺ or Fe²⁺ ions playing an important role.

Understanding the dynamics of monomeric, dimeric, and tetrameric α-synuclein structures in water

Human α‐synuclein (αS) is an intrinsically disordered protein associated with Parkinson's disease. Molecular mechanisms of corruptive misfolding and aggregation of αS resulting in the disease, as well as the structure and other properties of the corresponding oligomers are not entirely understood yet, preventing the development of efficient therapies. In this study, we investigate the folding dynamics of initially unfolded hypothetical αS constructs in water using all‐atom molecular dynamics simulations. We also employ the novel essential collective dynamics method to analyze the results obtained from the simulations. Our comparative analysis of monomeric, dimeric, and tetrameric αS models reveals pronounced differences in their structure and stability, emphasizing the importance of small oligomers, particularly dimers, in the process of misfolding.

Understanding interactions of functionalized nanoparticles with proteins: a case study on lactate dehydrogenase

Nanomaterials in biological solutions are known to interact with proteins and have been documented to affect protein function, such as enzyme activity. Understanding the interactions of nanoparticles with biological components at the molecular level will allow for rational designs of nanomaterials for use in medical technologies. Here we present the first detailed molecular mechanics model of functionalized gold nanoparticle (NP) interacting with an enzyme (L-lactate dehydrogenase (LDH) enzyme). Molecular dynamics (MD) simulations of the response of LDH to the NP binding demonstrate that although atomic motions (dynamics) of the main chain exhibit only a minor response to the binding, the dynamics of side chains are significantly constrained in all four active sites that predict alteration in kinetic properties of the enzyme. It is also demonstrated that the 5 nm gold NPs cause a decrease in the maximal velocity of the enzyme reaction (V(max)) and a trend towards a reduced affinity (increased K(m)) for the β-NAD binding site, while pyruvate enzyme kinetics (K(m) and V(max)) are not significantly altered in the presence of the gold NPs. These results demonstrate that modeling of NP:protein interactions can be used to understand alterations in protein function.

Evolutionary implications of metal binding features in different species' prion protein: an inorganic point of view

Prion disorders are a group of fatal neurodegenerative conditions of mammals. The key molecular event in the pathogenesis of such diseases is the conformational conversion of prion protein, PrPC, into a misfolded form rich in β-sheet structure, PrPSc, but the detailed mechanistic aspects of prion protein conversion remain enigmatic. There is uncertainty on the precise physiological function of PrPC in healthy individuals. Several evidences support the notion of its role in copper homeostasis. PrPC binds Cu2+ mainly through a domain composed by four to five repeats of eight amino acids. In addition to mammals, PrP homologues have also been identified in birds, reptiles, amphibians and fish. The globular domain of protein is retained in the different species, suggesting that the protein carries out an essential common function. However, the comparison of amino acid sequences indicates that prion protein has evolved differently in each vertebrate class. The primary sequences are strongly conserved in each group, but these exhibit a low similarity with those of mammals. The N-terminal domain of different prions shows tandem amino acid repeats with an increasing amount of histidine residues going from amphibians to mammals. The difference in the sequence affects the number of copper binding sites, the affinity and the coordination environment of metal ions, suggesting that the involvement of prion in metal homeostasis may be a specific characteristic of mammalian prion protein. In this review, we describe the similarities and the differences in the metal binding of different species' prion protein, as revealed by studies carried out on the entire protein and related peptide fragments.

Molecular mechanisms in the selective basal activation of pyrabactin receptor 1: Comparative analysis of mutants

Pyrabactin receptors (PYR) play a central role in abscisic acid (ABA) signal transduction; they are ABA receptors that inhibit type 2C protein phosphatases (PP2C). Molecular aspects contributing to increased basal activity of PYR against PP2C are studied by molecular dynamics (MD) simulations. An extensive series of MD simulations of the apo-form of mutagenized PYR1 as a homodimer and in complex with homology to ABA-insensitive 1 (HAB1) phosphatase are reported. In order to investigate the detailed molecular mechanisms mediating PYR1 activity, the MD data was analyzed by essential collective dynamics (ECD), a novel approach that allows the identification, with atomic resolution, of persistent dynamic correlations based on relatively short MD trajectories. Employing the ECD method, the effects of select mutations on the structure and dynamics of the PYR1 complexes were investigated and considered in the context of experimentally determined constitutive activities against HAB1. Approaches to rationally design constitutively active PYR1 constructs to increase PP2C inhibition are discussed.

Molecular mechanisms in the activation of abscisic acid receptor PYR1

The pyrabactin resistance 1 (PYR1)/PYR1-like (PYL)/regulatory component of abscisic acid (ABA) response (RCAR) proteins comprise a well characterized family of ABA receptors. Recent investigations have revealed two subsets of these receptors that, in the absence of ABA, either form inactive homodimers (PYR1 and PYLs 1–3) or mediate basal inhibition of downstream target type 2C protein phosphatases (PP2Cs; PYLs 4–10) respectively in vitro. Addition of ABA has been shown to release the apo-homodimers yielding ABA-bound monomeric holo-receptors that can interact with PP2Cs; highlighting a competitive-interaction process. Interaction selectivity has been shown to be mediated by subtle structural variations of primary sequence and ligand binding effects. Now, the dynamical contributions of ligand binding on interaction selectivity are investigated through extensive molecular dynamics (MD) simulations of apo and holo-PYR1 in monomeric and dimeric form as well as in complex with a PP2C, homology to ABA insensitive 1 (HAB1). Robust comparative interpretations were enabled by a novel essential collective dynamics approach. In agreement with recent experimental findings, our analysis indicates that ABA-bound PYR1 should efficiently bind to HAB1. However, both ABA-bound and ABA-extracted PYR1-HAB1 constructs have demonstrated notable similarities in their dynamics, suggesting that apo-PYR1 should also be able to make a substantial interaction with PP2Cs, albeit likely with slower complex formation kinetics. Further analysis indicates that both ABA-bound and ABA-free PYR1 in complex with HAB1 exhibit a higher intra-molecular structural stability and stronger inter-molecular dynamic correlations, in comparison with either holo- or apo-PYR1 dimers, supporting a model that includes apo-PYR1 in complex with HAB1. This possibility of a conditional functional apo-PYR1-PP2C complex was validated in vitro. These findings are generally consistent with the competitive-interaction model for PYR1 but highlight dynamical contributions of the PYR1 structure in mediating interaction selectivity suggesting added degrees of complexity in the regulation of the competitive-inhibition.

Protein pyrabactin resistance 1 (PYR1) belongs to a group of PYR1-like (PYL) proteins that regulate plant development and responses to conditions of drought and salinity. Recent studies have reported characterization of their molecular structures as well as elucidation of important aspects of their function; highlighting their roles as receptors for the stress responsive phytohormone, abscisic acid (ABA). However details of the molecular mechanisms regulating their receptor signalling remain enigmatic. In this work, we use molecular dynamics simulations complemented by a sophisticated statistical-mechanical analysis to investigate structural and dynamical properties of PYR1 protein and how its interaction with ABA modifies receptor-protein complex formation. Our results provide detailed insight into how the PYR1-mediated inactivation of its downstream phosphatase target is regulated by homodimer formation and yield new hypotheses, supported by in vitro experiments, for further investigation. Ultimately, this knowledge provides insight into how plants respond to stress, with potential applications in the development of crops with improved growth characteristics and higher stress tolerance.

Helices 2 and 3 are the initiation sites in the PrP(C) → PrP(SC) transition

It is established that prion protein is the sole causative agent in a number of diseases in humans and animals. However, the nature of conformational changes that the normal cellular form, PrP(C), undergoes in its conversion to a self-replicating state is still not fully understood. The ordered C-terminus of PrP(C) proteins has three helices (H1-H3). Here, we use statistical coupling analysis (SCA) to infer covariations at various locations using a family of evolutionarily related sequences and the response of mouse and human PrP(C)s to mechanical force to decipher the initiation sites for the transition from PrP(C) to an aggregation-prone PrP* state. Sequence-based SCA predicts that the clustered residues in nonmammals are localized in the stable core (near H1) of PrP(C), whereas in mammalian PrP(C), they are localized in frustrated helices H2 and H3 where most of the pathogenic mutations are found. Force-extension curves and free energy profiles as a function of extension of mouse and human PrP(C) in the absence of a disulfide (SS) bond between residues Cys179 and Cys214, generated by applying mechanical force to the ends of the molecule, show a sequence of unfolding events starting first with rupture of H2 and H3. This is followed by disruption of structure in two strands. Helix H1, stabilized by three salt bridges, resists substantial force before unfolding. Force extension profiles and the dynamics of rupture of tertiary contacts also show that even in the presence of an SS bond the instabilities in most of H3 and parts of H2 still determine the propensity to form the PrP* state. In mouse PrP(C) with an SS bond, there are ∼10 residues that retain their order even at high forces. Both SCA and single-molecule force simulations show that in the conversion from PrP(C) to PrP(SC) major conformational changes occur (at least initially) in H2 and H3, which because of their sequence compositions are frustrated in the helical state. Implications of our findings for the structural model for the scrapie form of PrP(C) are discussed.

The landscape of the prion protein's structural response to mutation revealed by principal component analysis of multiple NMR ensembles

Prion Proteins (PrP) are among a small number of proteins for which large numbers of NMR ensembles have been resolved for sequence mutants and diverse species. Here, we perform a comprehensive principle components analysis (PCA) on the tertiary structures of PrP globular proteins to discern PrP subdomains that exhibit conformational change in response to point mutations and clade-specific evolutionary sequence mutation trends. This is to our knowledge the first such large-scale analysis of multiple NMR ensembles of protein structures, and the first study of its kind for PrPs. We conducted PCA on human (n = 11), mouse (n = 14), and wildtype (n = 21) sets of PrP globular structures, from which we identified five conformationally variable subdomains within PrP. PCA shows that different non-local patterns and rankings of variable subdomains arise for different pathogenic mutants. These subdomains may thus be key areas for initiating PrP conversion during disease. Furthermore, we have observed the conformational clustering of divergent TSE-non-susceptible species pairs; these non-phylogenetic clusterings indicate structural solutions towards TSE resistance that do not necessarily coincide with evolutionary divergence. We discuss the novelty of our approach and the importance of PrP subdomains in structural conversion during disease.

Prion Proteins (PrP) cause a variety of incurable TSE diseases, and are among a small number of proteins for which large numbers of NMR ensembles have been resolved for sequence mutants and diverse species. Here, we perform a comprehensive PCA study to assess conformational variation and discern the landscape of the PrP structural response to sequence mutation. This is to our knowledge the first large-scale analysis of multiple NMR ensembles for a specific protein, and the first study to perform a multivariate PCA on the native globular structures of PrP. We conducted exhaustive PCA on three PrP subsets: human and mouse subsets that include structures of sequence mutants, and the set of wild-type PrP (16 PrP species). PCA shows that different non-local patterns of variable subdomains arise for different pathogenic mutants. These subdomains may thus be key areas for initiating PrP conversion during disease. Furthermore, we observed that some evolutionarily divergent species that are non-susceptible to TSEs have surprising structural similarities in their PrPs. We discuss the novelty of our approach with respect to prions, and the advantage of this analysis as a fast, reliable starting point to identify interesting domains that may warrant further experimental and computational analysis.

Exploring the essential collective dynamics of interacting proteins: application to prion protein dimers

Essential collective dynamics is a promising and robust approach for analysing the slow motions of macromolecules from short molecular dynamics trajectories. In this study, an extension of the method to treat a collection of interacting protein molecules is presented. The extension is applied to investigate the effects of dimerization on the collective dynamics of ovine prion protein molecules in two different arrangements. Examination of the structural plasticity shows that aggregation has a restricting effect on the local mobility of the prion protein molecules in the interfacial regions. Domain motions of the two dimeric ovine prion protein conformations are distinctly different and can be related to interatomic correlations at the interface. Correlated motions are among the slow collective modes extensively analysed by considering both main-chain and side-chain atoms. Correlation maps reveal the existence of a vast network of dynamically correlated side groups, which extends beyond individual subunits via interfacial interconnections. The network is formed by a core of hydrophobic side chains surrounded by hydrophilic groups at the periphery. The relevance of these findings are discussed in the context of mutations associated with prion diseases. The binding free energy of the dimeric conformations is evaluated to probe their thermodynamic stability. The descriptions afforded by the analysis of the essential collective dynamics of the prion dimers are consistent with their binding free energies. The agreement validates the extension of the methodology and provides a means of interpreting the collective dynamics in terms of the thermodynamic stability of ovine prion proteins.

Effects of temperature on the p53-DNA binding interactions and their dynamical behavior: comparing the wild type to the R248Q mutant

Background

The protein p53 plays an active role in the regulation of cell cycle. In about half of human cancers, the protein is inactivated by mutations located primarily in its DNA-binding domain. Interestingly, a number of these mutations possess temperature-induced DNA-binding characteristics. A striking example is the mutation of Arg248 into glutamine or tryptophan. These mutants are defective for binding to DNA at 310 K although they have been shown to bind specifically to several p53 response elements at sub-physiological temperatures (298–306 K).

Methodology/Principal Findings

This important experimental finding motivated us to examine the effects of temperature on the structure and configuration of R248Q mutant and compare it to the wild type protein. Our aim is to determine how and where structural changes of mutant variants take place due to temperature changes. To answer these questions, we compared the mutant to the wild-type proteins from two different aspects. First, we investigated the systems at the atomistic level through their DNA-binding affinity, hydrogen bond networks and spatial distribution of water molecules. Next, we assessed changes in their long-lived conformational motions at the coarse-grained level through the collective dynamics of their side-chain and backbone atoms separately.

Conclusions

The experimentally observed effect of temperature on the DNA-binding properties of p53 is reproduced. Analysis of atomistic and coarse-grained data reveal that changes in binding are determined by a few key residues and provide a rationale for the mutant-loss of binding at physiological temperatures. The findings can potentially enable a rescue strategy for the mutant structure.