Progress in studying intrinsically disordered proteins with atomistic simulations

Map conformational landscapes of intrinsically disordered proteins with polymer physics quantities

Disordered proteins are conformationally flexible proteins that are biologically important and have been implicated in devastating diseases such as Alzheimer's disease and cancer. Unlike stably folded structured proteins, disordered proteins sample a range of different conformations that needs to be accounted for. Here, we treat disordered proteins as polymer chains, and compute a dimensionless quantity called instantaneous shape ratio (Rs), as Rs = Ree2/Rg2, where Ree is end-to-end distance and Rg is radius of gyration. Extended protein conformations tend to have high Ree compared with Rg, and thus have high Rs values, whereas compact conformations have smaller Rs values. We use a scatter plot of Rs (representing shape) against Rg (representing size) as a simple map of conformational landscapes. We first examine the conformational landscape of simple polymer models such as Random Walk, Self-Avoiding Walk, and Gaussian Walk (GW), and we notice that all protein/polymer maps lie within the boundaries of the GW map. We thus use the GW map as a reference and, to assess conformational diversity, we compute the fraction of the GW conformations (fC) covered by each protein/polymer. Disordered proteins all have high fC scores, consistent with their disordered nature. Each disordered protein accesses a different region of the reference map, revealing differences in their conformational ensembles. We additionally examine the conformational maps of the nonviral gene delivery vector polyethyleneimine at various protonation states, and find that they resemble disordered proteins, with coverage of the reference map decreasing with increasing protonation state, indicating decreasing conformational diversity. We propose that our method of combining Rs and Rg in a scatter plot generates a simple, meaningful map of the conformational landscape of a disordered protein, which in turn can be used to assess conformational diversity of disordered proteins.

Exploring the Role of Globular Domain Locations on an Intrinsically Disordered Region of p53: A Molecular Dynamics Investigation

The pre-tetramerization loop (PTL) of the human tumor suppressor protein p53 is an intrinsically disordered region (IDR) necessary for the tetramerization process, and its flexibility contributes to the essential conformational changes needed. Although the IDR can be accurately simulated in the traditional manner of molecular dynamics (MD) with the end-to-end distance (EEdist) unhindered, we sought to explore the effects of restraining the EEdist to the values predicted by electron microscopy (EM) and other distances. Simulating the PTL trajectory with a restrained EEdist , we found an increased agreement of nuclear magnetic resonance (NMR) chemical shifts with experiments. Additionally, we observed a plethora of secondary structures and contacts that only appear when the trajectory is restrained. Our findings expand the understanding of the tetramerization of p53 and provide insight into how mutations could make the protein impotent. In particular, our findings demonstrate the importance of restraining the EEdist in studying IDRs and how their conformations change under different conditions. Our results provide a better understanding of the PTL and the conformational dynamics of IDRs in general, which are useful for further studies regarding mutations and their effects on the activity of p53.

Interaction Mechanism between α-Lactalbumin and Caffeic Acid: A Multispectroscopic and Molecular Docking Study

Caffeic acid (CA) is a phenolic acid found in a variety of foods. In this study, the interaction mechanism between α-lactalbumin (ALA) and CA was explored with the use of spectroscopic and computational techniques. The Stern-Volmer quenching constant data suggest a static mode of quenching between CA and ALA, depicting a gradual decrease in quenching constants with temperature rise. The binding constant, Gibbs free energy, enthalpy, and entropy values at 288, 298, and 310 K were calculated, and the obtained values suggest that the reaction is spontaneous and exothermic. Both in vitro and in silico studies show that hydrogen bonding is the dominant force in the CA-ALA interaction. Ser112 and Lys108 of ALA are predicted to form three hydrogen bonds with CA. The UV-visible spectroscopy measurements demonstrated that the absorbance peak A_280nm increased after addition of CA due to conformational change. The secondary structure of ALA was also slightly modified due to CA interaction. The circular dichroism (CD) studies showed that ALA gains more α-helical structure in response to increasing concentration of CA. The surface hydrophobicity of ALA is not changed in the presence of ethanol and CA. The present findings shown herein are helpful in understanding the binding mechanism of CA with whey proteins for the dairy processing industry and food nutrition security.

Exploring Structures and Dynamics of Protamine Molecules through Molecular Dynamics Simulations

Protamines are arginine-rich proteins that condense DNA in sperm. Despite their importance in reproduction, information on protamine structure is scarce. We, therefore, used molecular dynamics to examine the structures of salmon, bull P1, and human P1 protamines. The sizes and shapes of each protamine varied widely, indicating that they were disordered with structures covering a broad conformational landscape, from hairpin loop structures to extended coils. Despite their general disorder, the protamines did form secondary structures, including helices and hairpin loops. In eutherians, hairpins may promote disulfide bonding that facilitates protamine-DNA condensation, but the specifics of this bonding is not well established. We examined inter-residue distances in the simulations to predict residue pairs likely to form intramolecular bonds, leading to the identification of bonding pairs consistent with previous results in bull and human. These results support a model for eutherian protamine structures where a highly charged center is surrounded by disulfide-bond-stabilized loops.

Exploring the conformational dynamics and flexibility of intrinsically disordered HIV-1 Nef protein using molecular dynamic network approaches

Intrinsically disordered proteins represent a class of proteins that lack fixed and well-defined three-dimensional structures in solution. HIV-1 Nef is an intrinsically disordered peripheral membrane protein involved in the replication and pathogenesis of HIV-1 infection. Nef controls expression levels of cell surface CD4 molecules that are essential for adaptive immunity. Despite the lack of fixed and stable structures, Nef physically interacts with the host cellular proteins (AP-1/MHC-I) and modulates intracellular trafficking pathways. Therefore, it is essential to understand how this dynamic conformational flexibility affects Nef structures and function. In this study, we combined all-atom molecular dynamics (MD) simulations and dynamic network approaches to better understand the structure and dynamics of Nef in two different forms, the free unbound and the bound state. Using the MD simulation approach, we show that the intrinsically disordered Nef exhibit a large dynamic field with more atomic fluctuations and lesser thermodynamic stability in the unbound conditions. The conformations of Nef change over time, and this protein remains more compact, folded, and stable in the bound form. The dynamic network analysis revealed regions of the protein capable of modulating the conformational behavior of the disordered Nef. The average betweenness centrality (BC) unveiled residues that are critical for mediating protein-protein interactions. The average shortest path length (L) and the perturbation response scanning exposed residues that are likely to be important in steering protein conformational changes. Overall, the study demonstrates how all-atom MD simulations combined with the dynamic network approach can be used to gain further insights into the structure and dynamics-function relationship of intrinsically disordered HIV-1 Nef.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02698-8.

Recent advances in atomic molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) play important roles in cellular functions. The inherent structural heterogeneity of IDPs makes the high-resolution experimental characterization of IDPs extremely difficult. Molecular dynamics (MD) simulation could provide the atomic-level description of the structural and dynamic properties of IDPs. This perspective reviews the recent progress in atomic MD simulation studies of IDPs, including the development of force fields and sampling methods, as well as applications in IDP-involved protein-protein interactions. The employment of large-scale simulations and advanced sampling techniques allows more accurate estimation of the thermodynamics and kinetics of IDP-mediated protein interactions, and the holistic landscape of the binding process of IDPs is emerging.

Intrinsic Disorder in the T Cell Receptor Creates Cooperativity and Controls ZAP70 Binding

Many immunoreceptors have cytoplasmic domains that are intrinsically disordered (i.e., have high configurational entropy), have multiple sites of posttranslational modification (e.g., tyrosine phosphorylation), and participate in nonlinear signaling pathways (e.g., exhibiting switch-like behavior). Several hypotheses to explain the origin of these nonlinearities fall under the broad hypothesis that modification at one site changes the immunoreceptor’s entropy, which in turn changes further modification dynamics. Here, we use coarse-grain simulation to study three scenarios, all related to the chains that constitute the T cell receptor (TCR). We find that first, if phosphorylation induces local changes in the flexibility of the TCR ζ-chain, this naturally leads to rate enhancements and cooperativity. Second, we find that TCR CD3ɛ can provide a switch by modulating its residence in the plasma membrane. By constraining our model to be consistent with the previous observation that both basic residues and phosphorylation control membrane residence, we find that there is only a moderate rate enhancement of 10% between first and subsequent phosphorylation events. Third, we find that volume constraints do not limit the number of ZAP70s that can bind the TCR but that entropic penalties lead to a 200-fold decrease in binding rate by the seventh ZAP70, potentially explaining the observation that each TCR has around six ZAP70 molecules bound after receptor triggering. In all three scenarios, our results demonstrate that phenomena that change an immunoreceptor chain’s entropy (stiffening, confinement to a membrane, and multiple simultaneous binding) can lead to nonlinearities (rate enhancement, switching, and negative cooperativity) in how the receptor participates in signaling. These polymer-entropy-driven nonlinearities may augment the nonlinearities that arise from, e.g., kinetic proofreading and cluster formation. They also suggest different design strategies for engineered receptors, e.g., whether or not to put signaling modules on one chain or multiple clustered chains.

Significant compaction of H4 histone tail upon charge neutralization by acetylation and its mimics, possible effects on chromatin structure

The intrinsically disordered, positively charged H4 histone tail is important for chromatin structure and function. We have explored conformational ensembles of human H4 tail in solution, with varying levels of charge neutralization via acetylation or amino-acid substitutions such as K→Q. We have employed an explicit water model shown recently to be well suited for simulations of intrinsically disordered proteins. Upon progressive neutralization of the H4, its radius of gyration decreases linearly with the tail charge q, the trend is explained using a simple polymer model. While the wild type state (q=+8) is essentially a random coil, hyper-acetylated H4 (q=+3) is virtually as compact and stable as a globular protein of the same number of amino-acids. Conformational ensembles of acetylated H4 match the corresponding K→X substitutions only approximately: based on the ensemble similarity, we propose K→M as a possible alternative to the commonly used K→Q. Possible effects of the H4 tail compaction on chromatin structure are discussed within a qualitative model in which the chromatin is highly heterogeneous, easily inter-converting between various structural forms. We predict that upon progressive charge neutralization of the H4 tail, the least compact sub-states of chromatin de-condense first, followed by de-condensation of more compact structures, e.g. those that harbor a high fraction of stacked di-nucleosomes. The predicted hierarchy of DNA accessibility increase upon progressive acetylation of H4 might be utilized by the cell for selective DNA accessibility control.

Sequence specificity despite intrinsic disorder: How a disease-associated Val/Met polymorphism rearranges tertiary interactions in a long disordered protein

The role of electrostatic interactions and mutations that change charge states in intrinsically disordered proteins (IDPs) is well-established, but many disease-associated mutations in IDPs are charge-neutral. The Val66Met single nucleotide polymorphism (SNP) in precursor brain-derived neurotrophic factor (BDNF) is one of the earliest SNPs to be associated with neuropsychiatric disorders, and the underlying molecular mechanism is unknown. Here we report on over 250 μs of fully-atomistic, explicit solvent, temperature replica-exchange molecular dynamics (MD) simulations of the 91 residue BDNF prodomain, for both the V66 and M66 sequence. The simulations were able to correctly reproduce the location of both local and non-local secondary structure changes due to the Val66Met mutation, when compared with NMR spectroscopy. We find that the change in local structure is mediated via entropic and sequence specific effects. We developed a hierarchical sequence-based framework for analysis and conceptualization, which first identifies “blobs” of 4-15 residues representing local globular regions or linkers. We use this framework within a novel test for enrichment of higher-order (tertiary) structure in disordered proteins; the size and shape of each blob is extracted from MD simulation of the real protein (RP), and used to parameterize a self-avoiding heterogenous polymer (SAHP). The SAHP version of the BDNF prodomain suggested a protein segmented into three regions, with a central long, highly disordered polyampholyte linker separating two globular regions. This effective segmentation was also observed in full simulations of the RP, but the Val66Met substitution significantly increased interactions across the linker, as well as the number of participating residues. The Val66Met substitution replaces β-bridging between V66 and V94 (on either side of the linker) with specific side-chain interactions between M66 and M95. The protein backbone in the vicinity of M95 is then free to form β-bridges with residues 31-41 near the N-terminus, which condenses the protein. A significant role for Met/Met interactions is consistent with previously-observed non-local effects of the Val66Met SNP, as well as established interactions between the Met66 sequence and a Met-rich receptor that initiates neuronal growth cone retraction.

Intrinsically disordered proteins are proteins that have no well-defined structure in at least one functional form. Mutations in one amino acid may still affect their function significantly, especially in subtle ways with cumulative adverse effects on health. Here we report on molecular dynamics simulations of a protein that is critical for neuronal health throughout adulthood (brain-derived neurotrophic factor). We investigate the effects of a mutation carried by 30% of human population, which has been widely studied for its association with aging-related and stress-related disorders, reduced volume of the hippocampus, and variations in episodic memory. We identify a molecular mechanism in which the mutation may change the global conformations of the protein and its ability to bind to receptors.

An intrinsically disordered proteins community for ELIXIR

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) are now recognised as major determinants in cellular regulation. This white paper presents a roadmap for future e-infrastructure developments in the field of IDP research within the ELIXIR framework. The goal of these developments is to drive the creation of high-quality tools and resources to support the identification, analysis and functional characterisation of IDPs. The roadmap is the result of a workshop titled “An intrinsically disordered protein user community proposal for ELIXIR” held at the University of Padua. The workshop, and further consultation with the members of the wider IDP community, identified the key priority areas for the roadmap including the development of standards for data annotation, storage and dissemination; integration of IDP data into the ELIXIR Core Data Resources; and the creation of benchmarking criteria for IDP-related software. Here, we discuss these areas of priority, how they can be implemented in cooperation with the ELIXIR platforms, and their connections to existing ELIXIR Communities and international consortia. The article provides a preliminary blueprint for an IDP Community in ELIXIR and is an appeal to identify and involve new stakeholders.

Novel phosphorylation-dependent regulation in an unstructured protein

This work explores how phosphorylation of an unstructured protein region in inhibitor-2 (I2) regulates protein phosphatase-1 (PP1) enzyme activity using molecular dynamics (MD). Free I2 is largely unstructured; however, when bound to PP1, three segments adopt a stable structure. In particular, an I2 helix (i-helix) blocks the PP1 active site and inhibits phosphatase activity. I2 phosphorylation in the PP1-I2 complex activates phosphatase activity without I2 dissociation. The I2 Thr74 regulatory phosphorylation site is in an unstructured domain in PP1-I2. PP1-I2 MD demonstrated that I2 phosphorylation promotes early steps of PP1-I2 activation in explicit solvent models. Moreover, phosphorylation-dependent activation occurred in PP1-I2 complexes derived from I2 orthologs with diverse sequences from human, yeast, worm, and protozoa. This system allowed exploration of features of the 73-residue unstructured human I2 domain critical for phosphorylation-dependent activation. These studies revealed that components of I2 unstructured domain are strategically positioned for phosphorylation responsiveness including a transient α-helix. There was no evidence that electrostatic interactions of I2 phosphothreonine74 influenced PP1-I2 activation. Instead, phosphorylation altered the conformation of residues around Thr74. Phosphorylation uncurled the distance between I2 residues Glu71 to Tyr76 to promote PP1-I2 activation, whereas reduced distances reduced activation. This I2 residue Glu71 to Tyr76 distance distribution, independently from Thr74 phosphorylation, controls I2 i-helix displacement from the PP1 active site leading to PP1-I2 activation.

Studies on Molecular Dynamics of Intrinsically Disordered Proteins and Their Fuzzy Complexes: A Mini-Review

The molecular dynamics (MD) method is a promising approach toward elucidating the molecular mechanisms of intrinsically disordered regions (IDRs) of proteins and their fuzzy complexes. This mini-review introduces recent studies that apply MD simulations to investigate the molecular recognition of IDRs. Firstly, methodological issues by which MD simulations treat IDRs, such as developing force fields, treating periodic boundary conditions, and enhanced sampling approaches, are discussed. Then, several examples of the applications of MD to investigate molecular interactions of IDRs in terms of the two kinds of complex formations; coupled-folding and binding and fuzzy complex. MD simulations provide insight into the molecular mechanisms of these binding processes by sampling conformational ensembles of flexible IDRs. In particular, we focused on all-atom explicit-solvent MD simulations except for studies of higher-order assembly of IDRs. Recent advances in MD methods, and computational power make it possible to dissect the molecular details of realistic molecular systems involving the dynamic behavior of IDRs.

General Purpose Water Model Can Improve Atomistic Simulations of Intrinsically Disordered Proteins

Unconstrained atomistic simulations of intrinsically disordered proteins and peptides (IDP) remain a challenge: widely used, "general purpose" water models tend to favor overly compact structures relative to experiment. Here we have performed a total of 93 μs of unrestrained MD simulations to explore, in the context of IDPs, a recently developed "general-purpose" 4-point rigid water model OPC, which describes liquid state of water close to experiment. We demonstrate that OPC, together with a popular AMBER force field ff99SB, offers a noticeable improvement over TIP3P in producing more realistic structural ensembles of three common IDPs benchmarks: 55-residue apo N-terminal zinc-binding domain of HIV-1 integrase ("protein IN"), amyloid β-peptide (Aβ₄₂) (residues 1-42), and 26-reside H4 histone tail. As a negative control, computed folding profile of a regular globular miniprotein (CLN025) in OPC water is in appreciably better agreement with experiment than that obtained in TIP3P, which tends to overstabilize the compact native state relative to the extended conformations. We employed Aβ₄₂ peptide to investigate the possible influence of the solvent box size on simulation outcomes. We advocate a cautious approach for simulations of IDPs: we suggest that the solvent box size should be at least four times the radius of gyration of the random coil corresponding to the IDP. The computed free energy landscape of protein IN in OPC resembles a shallow "tub" - conformations with substantially different degrees of compactness that are within 2 k_B T of each other. Conformations with very different secondary structure content coexist within 1 k_B T of the global free energy minimum. States with higher free energy tend to have less secondary structure. Computed low helical content of the protein has virtually no correlation with its degree of compactness, which calls into question the possibility of using the helicity as a metric for assessing performance of water models for IDPs, when the helicity is low. Predicted radius of gyration ( R _g) of H4 histone tail in OPC water falls in-between that of a typical globular protein and a fully denatured protein of the same size; the predicted R _g is consistent with two independent predictions. In contrast, H4 tail in TIP3P water is as compact as the corresponding globular protein. The computed free energy landscape of H4 tail in OPC is relatively flat over a significant range of compactness, which, we argue, is consistent with its biological function as facilitator of internucleosome interactions.

Extreme Fuzziness: Direct Interactions between Two IDPs

Protein interactions involving intrinsically disordered proteins (IDPs) greatly extend the range of binding mechanisms available to proteins. In interactions employing coupled folding and binding, IDPs undergo disorder-to-order transitions to form a complex with a well-defined structure. In many other cases, IDPs retain structural plasticity in the final complexes, which have been defined as the fuzzy complexes. While a large number of fuzzy complexes have been characterized with variety of fuzzy patterns, many of the interactions are between an IDP and a structured protein. Thus, whether two IDPs can interact directly to form a fuzzy complex without disorder-to-order transition remains an open question. Recently, two studies of interactions between IDPs (4.1G-CTD/NuMA and H1/ProTα) have found a definite answer to this question. Detailed characterizations combined with nuclear magnetic resonance (NMR), single-molecule Förster resonance energy transfer (smFRET) and molecular dynamics (MD) simulation demonstrate that direct interactions between these two pairs of IDPs do form fuzzy complexes while retaining the conformational dynamics of the isolated proteins, which we name as the extremely fuzzy complexes. Extreme fuzziness completes the full spectrum of protein-protein interaction modes, suggesting that a more generalized model beyond existing binding mechanisms is required. Previous models of protein interaction could be applicable to some aspects of the extremely fuzzy interactions, but in more general sense, the distinction between native and nonnative contacts, which was used to understand protein folding and binding, becomes obscure. Exploring the phenomenon of extreme fuzziness may shed new light on molecular recognition and drug design.

Validating Molecular Dynamics Simulations against Experimental Observables in Light of Underlying Conformational Ensembles

Far from the static, idealized conformations deposited into structural databases, proteins are highly dynamic molecules that undergo conformational changes on temporal and spatial scales that may span several orders of magnitude. These conformational changes, often intimately connected to the functional roles that proteins play, may be obscured by traditional biophysical techniques. Over the past 40 years, molecular dynamics (MD) simulations have complemented these techniques by providing the "hidden" atomistic details that underlie protein dynamics. However, there are limitations of the degree to which molecular simulations accurately and quantitatively describe protein motions. Here we show that although four molecular dynamics simulation packages (AMBER, GROMACS, NAMD, and ilmm) reproduced a variety of experimental observables for two different proteins (engrailed homeodomain and RNase H) equally well overall at room temperature, there were subtle differences in the underlying conformational distributions and the extent of conformational sampling obtained. This leads to ambiguity about which results are correct, as experiment cannot always provide the necessary detailed information to distinguish between the underlying conformational ensembles. However, the results with different packages diverged more when considering larger amplitude motion, for example, the thermal unfolding process and conformational states sampled, with some packages failing to allow the protein to unfold at high temperature or providing results at odds with experiment. While most differences between MD simulations performed with different packages are attributed to the force fields themselves, there are many other factors that influence the outcome, including the water model, algorithms that constrain motion, how atomic interactions are handled, and the simulation ensemble employed. Here four different MD packages were tested each using best practices as established by the developers, utilizing three different protein force fields and three different water models. Differences between the simulated protein behavior using two different packages but the same force field, as well as two different packages with different force fields but the same water models and approaches to restraining motion, show how other factors can influence the behavior, and it is incorrect to place all the blame for deviations and errors on force fields or to expect improvements in force fields alone to solve such problems.

Looking at the Disordered Proteins through the Computational Microscope

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.

Computational and theoretical advances in studies of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are increasingly realized to play diverse biological roles, ranging from molecular signaling to the formation of membraneless organelles. Their high degree of disorder makes them more challenging to study using the techniques of conventional structural biology, because any observable will be averaged over a heterogeneous ensemble of structures. Molecular simulations and theory are therefore a natural complement to experiment for studying the structure, dynamics and function of IDPs. The diverse time and length scales relevant to the roles played by IDPs require flexibility in the techniques applied. Here, I summarize some of the developments in simulation and theory in recent years, which have been driven by the desire to better capture IDP properties at different time- and length-scales. I also provide an outlook for how methods can be improved in the future and emerging problems which may be addressed by theory and simulation.

Computer Simulations of Intrinsically Disordered Proteins

The investigation of intrinsically disordered proteins (IDPs) is a new frontier in structural and molecular biology that requires a new paradigm to connect structural disorder to function. Molecular dynamics simulations and statistical thermodynamics potentially offer ideal tools for atomic-level characterizations and thermodynamic descriptions of this fascinating class of proteins that will complement experimental studies. However, IDPs display sensitivity to inaccuracies in the underlying molecular mechanics force fields. Thus, achieving an accurate structural characterization of IDPs via simulations is a challenge. It is also daunting to perform a configuration-space integration over heterogeneous structural ensembles sampled by IDPs to extract, in particular, protein configurational entropy. In this review, we summarize recent efforts devoted to the development of force fields and the critical evaluations of their performance when applied to IDPs. We also survey recent advances in computational methods for protein configurational entropy that aim to provide a thermodynamic link between structural disorder and protein activity.

Comparison of force fields for Alzheimer's A β42: A case study for intrinsically disordered proteins

Intrinsically disordered proteins are essential for biological processes such as cell signalling, but are also associated to devastating diseases including Alzheimer's disease, Parkinson's disease or type II diabetes. Because of their lack of a stable three-dimensional structure, molecular dynamics simulations are often used to obtain atomistic details that cannot be observed experimentally. The applicability of molecular dynamics simulations depends on the accuracy of the force field chosen to represent the underlying free energy surface of the system. Here, we use replica exchange molecular dynamics simulations to test five modern force fields, OPLS, AMBER99SB, AMBER99SB*ILDN, AMBER99SBILDN-NMR and CHARMM22*, in their ability to model Aβ42 , an intrinsically disordered peptide associated with Alzheimer's disease, and compare our results to nuclear magnetic resonance (NMR) experimental data. We observe that all force fields except AMBER99SBILDN-NMR successfully reproduce local NMR observables, with CHARMM22* being slightly better than the other force fields.

Molecular Dynamics Simulation of Tau Peptides for the Investigation of Conformational Changes Induced by Specific Phosphorylation Patterns

The Tau protein plays an important role due to its biomolecular interactions in neurodegenerative diseases. The lack of stable structure and various posttranslational modifications such as phosphorylation at various sites in the Tau protein pose a challenge for many experimental methods that are traditionally used to study protein folding and aggregation. Atomistic molecular dynamics (MD) simulations can help around deciphering relationship between phosphorylation and various intermediate and stable conformations of the Tau protein which occur on longer timescales. This chapter outlines protocols for the preparation, execution, and analysis of all-atom MD simulations of a 21-amino acid-long phosphorylated Tau peptide with the aim of generating biologically relevant structural and dynamic information. The simulations are done in explicit solvent and starting from nearly extended configurations of the peptide. The scaled MD method implemented in AMBER14 was chosen to achieve enhanced conformational sampling in addition to a conventional MD approach, thereby allowing the characterization of folding for such an intrinsically disordered peptide at 293 K. Emphasis is placed on the analysis of the simulation trajectories to establish correlations with NMR data (i.e., chemical shifts and NOEs). Finally, in-depth discussions are provided for commonly encountered problems.

Comprehensive structural and dynamical view of an unfolded protein from the combination of single-molecule FRET, NMR, and SAXS

The properties of unfolded proteins are essential both for the mechanisms of protein folding and for the function of the large group of intrinsically disordered proteins. However, the detailed structural and dynamical characterization of these highly dynamic and conformationally heterogeneous ensembles has remained challenging. Here we combine and compare three of the leading techniques for the investigation of unfolded proteins, NMR spectroscopy (NMR), small-angle X-ray scattering (SAXS), and single-molecule Förster resonance energy transfer (FRET), with the goal of quantitatively testing their consistency and complementarity and for obtaining a comprehensive view of the unfolded-state ensemble. Using unfolded ubiquitin as a test case, we find that its average dimensions derived from FRET and from structural ensembles calculated using the program X-PLOR-NIH based on NMR and SAXS restraints agree remarkably well; even the shapes of the underlying intramolecular distance distributions are in good agreement, attesting to the reliability of the approaches. The NMR-based results provide a highly sensitive way of quantifying residual structure in the unfolded state. FRET-based nanosecond fluorescence correlation spectroscopy allows long-range distances and chain dynamics to be probed in a time range inaccessible by NMR. The combined techniques thus provide a way of optimally using the complementarity of the available methods for a quantitative structural and dynamical description of unfolded proteins both at the global and the local level.

Binding of Disordered Peptides to Kelch: Insights from Enhanced Sampling Simulations

Keap1 protein plays an essential role in regulating cellular oxidative stress response and is a crucial binding hub for multiple proteins, several of which are intrinsically disordered proteins (IDP). Among Kelch's IDP binding partners, NRF2 and PTMA are the two most interesting cases. They share a highly similar binding motif; however, NRF2 binds to Kelch with a binding affinity of approximately 100-fold higher than that of PTMA. In this study, we perform an exhaustive sampling composed of 6 μs well-tempered metadynamics and 2 μs unbiased molecular dynamics (MD) simulations aiming at characterizing the binding mechanisms and structural properties of these two peptides. Our results agree with previous experimental observations that PTMA is remarkably more disordered than NRF2 in both the free and bound states. This explains PTMA's lower binding affinity. Our extensive sampling also provides valuable insights into the vast conformational ensembles of both NRF2 and PTMA, supports the hypothesis of coupled folding-binding, and confirms the essential role of linear motifs in IDP binding.