Conformational recognition of an intrinsically disordered protein

Model-Dependent Solvation of the K-18 Domain of the Intrinsically Disordered Protein Tau

A known imbalance between intra-protein and protein-water interactions in many empirical force fields results in collapsed conformational ensembles of intrinsically disordered proteins in explicit solvent simulations that disagree with experiments. Multiple strategies have been introduced in the literature to modify protein-water interactions, which improve agreement between experiments and simulations. In this work, we combine simulations with standard and modified force fields with a spatially resolved analysis of solvation free energy contributions and compare the consequences of each strategy. We find that enhanced Lennard-Jones (LJ) interactions between protein atoms and water oxygens primarily improve the solvation of nonpolar functional groups of the protein. In contrast, modified electrostatics in the water model or strengthened LJ interactions between the protein and water hydrogens mainly affect the hydration of polar functional groups. Modified electrostatics further impact the average orientation of water molecules in the hydration shell. As a result, protein-water interactions with the first hydration layers are strengthened, while interactions with water molecules in higher hydration shells are weakened. Hence, distinct strategies to balance intra-protein and protein-water interactions in simulations have qualitatively different effects on protein solvation. These differences are not necessarily captured by comparisons to experiments that report on global parameters describing protein conformational ensembles, e.g., the radius of gyration, but will influence the tendency of a protein to form aggregates or phase-separated droplets.

The influence of random-coil chemical shifts on the assessment of structural propensities in folded proteins and IDPs

In studying secondary structural propensities of proteins by nuclear magnetic resonance (NMR) spectroscopy, secondary chemical shifts (SCSs) serve as the primary atomic scale observables. For SCS calculation, the selection of an appropriate random coil chemical shift (RCCS) dataset is a crucial step, especially when investigating intrinsically disordered proteins (IDPs). The scientific literature is abundant in such datasets, however, the effect of choosing one over all the others in a concrete application has not yet been studied thoroughly and systematically. Hereby, we review the available RCCS prediction methods and to compare them, we conduct statistical inference by means of the nonparametric sum of ranking differences and comparison of ranks to random numbers (SRD-CRRN) method. We try to find the RCCS predictors best representing the general consensus regarding secondary structural propensities. The existence and the magnitude of resulting differences on secondary structure determination under varying sample conditions (temperature, pH) are demonstrated and discussed for globular proteins and especially IDPs.

Evolutionary fine-tuning of residual helix structure in disordered proteins manifests in complex structure and lifetime

Transcription depends on complex networks, where folded hub proteins interact with intrinsically disordered transcription factors undergoing coupled folding and binding. For this, local residual structure, a prototypical feature of intrinsic disorder, is key. Here, we dissect the unexplored functional potential of residual structure by comparing structure, kinetics, and thermodynamics within the model system constituted of the DREB2A transcription factor interacting with the αα-hub RCD1-RST. To maintain biological relevance, we developed an orthogonal evolutionary approach for the design of variants with varying amounts of structure. Biophysical analysis revealed a correlation between the amount of residual helical structure and binding affinity, manifested in altered complex lifetime due to changed dissociation rate constants. It also showed a correlation between helical structure in free and bound DREB2A variants. Overall, this study demonstrated how evolution can balance and fine-tune residual structure to regulate complexes in coupled folding and binding, potentially affecting transcription factor competition.

Intrinsically disordered proteins play diverse roles in cell signaling

Abstract

Signaling pathways allow cells to detect and respond to a wide variety of chemical (e.g. Ca²⁺ or chemokine proteins) and physical stimuli (e.g., sheer stress, light). Together, these pathways form an extensive communication network that regulates basic cell activities and coordinates the function of multiple cells or tissues. The process of cell signaling imposes many demands on the proteins that comprise these pathways, including the abilities to form active and inactive states, and to engage in multiple protein interactions. Furthermore, successful signaling often requires amplifying the signal, regulating or tuning the response to the signal, combining information sourced from multiple pathways, all while ensuring fidelity of the process. This sensitivity, adaptability, and tunability are possible, in part, due to the inclusion of intrinsically disordered regions in many proteins involved in cell signaling. The goal of this collection is to highlight the many roles of intrinsic disorder in cell signaling. Following an overview of resources that can be used to study intrinsically disordered proteins, this review highlights the critical role of intrinsically disordered proteins for signaling in widely diverse organisms (animals, plants, bacteria, fungi), in every category of cell signaling pathway (autocrine, juxtacrine, intracrine, paracrine, and endocrine) and at each stage (ligand, receptor, transducer, effector, terminator) in the cell signaling process. Thus, a cell signaling pathway cannot be fully described without understanding how intrinsically disordered protein regions contribute to its function. The ubiquitous presence of intrinsic disorder in different stages of diverse cell signaling pathways suggest that more mechanisms by which disorder modulates intra- and inter-cell signals remain to be discovered.

Graphical abstract

Force Field Effects in Simulations of Flexible Peptides with Varying Polyproline II Propensity

Five peptides previously suggested to possess polyproline II (PPII) structure have here been investigated by using atomistic molecular dynamics simulations to compare how well four different force fields known for simulating intrinsically disordered proteins relatively well (Amber ff99SB-disp, Amber ff99SB-ILDN, CHARM36IDPSFF, and CHARMM36m) can capture this secondary structure element. The results revealed that all force fields sample PPII structures but to different extents and with different propensities toward other secondary structure elements, in particular, the β-sheet and "random coils". A cluster analysis of the simulations of histatin 5 also revealed that the conformational ensembles of the force fields are quite different. We compared the simulations to circular dichroism and nuclear magnetic resonance spectroscopy experiments and conclude that further experiments and methods for interpreting them are needed to assess the accuracy of force fields in determining PPII structure.

The diversity of molecular interactions involving intrinsically disordered proteins: A molecular modeling perspective

Intrinsically Disordered Proteins and Regions (IDPs/IDRs) are key components of a multitude of biological processes. Conformational malleability enables IDPs/IDRs to perform very specialized functions that cannot be accomplished by globular proteins. The functional role for most of these proteins is related to the recognition of other biomolecules to regulate biological processes or as a part of signaling pathways. Depending on the extent of disorder, the number of interacting sites and the type of partner, very different architectures for the resulting assemblies are possible. More recently, molecular condensates with liquid-like properties composed of multiple copies of IDPs and nucleic acids have been proven to regulate key processes in eukaryotic cells. The structural and kinetic details of disordered biomolecular complexes are difficult to unveil experimentally due to their inherent conformational heterogeneity. Computational approaches, alone or in combination with experimental data, have emerged as unavoidable tools to understand the functional mechanisms of this elusive type of assemblies. The level of description used, all-atom or coarse-grained, strongly depends on the size of the molecular systems and on the timescale of the investigated mechanism. In this mini-review, we describe the most relevant architectures found for molecular interactions involving IDPs/IDRs and the computational strategies applied for their investigation.

Tackling Covid-19 using disordered-to-order transition of residues in the spike protein upon angiotensin-converting enzyme 2 binding

The 2019-novel coronavirus also known as severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) is a common threat to animals and humans, and is responsible for the human SARS pandemic in 2019 to 2021. The infection of SARS-CoV-2 in humans involves a viral surface glycoprotein named as spike proteins, which bind to the human angiotensin-converting enzyme 2 (ACE2) proteins. Particularly, the receptor binding domains (RBDs) mediate the interaction and contain several disordered regions, which help in the binding. Investigations on the influence of disordered residues/regions in stability and binding of spike protein with ACE2 help to understand the disease pathogenesis, which has not yet been studied. In this study, we have used molecular-dynamics simulations to characterize the structural changes in disordered regions of the spike protein that result from ACE2 binding. We observed that the disordered regions undergo disorder-to-order transition (DOT) upon binding with ACE2, and the DOT residues are located at functionally important regions of RBD. Although the RBD is having rigid structure, DOT residues make conformational rearrangements for the spike protein to attach with ACE2. The binding is strengthened via hydrophilic and aromatic amino acids mainly present in the DOTs. The positively correlated motions of the DOT residues with its nearby residues also explain the binding profile of RBD with ACE2, and the residues are observed to be contributing more favorable binding energies for the spike-ACE2 complex formation. This study emphasizes that intrinsically disordered residues in the RBD of spike protein may provide insights into its etiology and be useful for drug and vaccine discovery.

A disordered encounter complex is central to the yeast Abp1p SH3 domain binding pathway

Protein-protein interactions are involved in a wide range of cellular processes. These interactions often involve intrinsically disordered proteins (IDPs) and protein binding domains. However, the details of IDP binding pathways are hard to characterize using experimental approaches, which can rarely capture intermediate states present at low populations. SH3 domains are common protein interaction domains that typically bind proline-rich disordered segments and are involved in cell signaling, regulation, and assembly. We hypothesized, given the flexibility of SH3 binding peptides, that their binding pathways include multiple steps important for function. Molecular dynamics simulations were used to characterize the steps of binding between the yeast Abp1p SH3 domain (AbpSH3) and a proline-rich IDP, ArkA. Before binding, the N-terminal segment 1 of ArkA is pre-structured and adopts a polyproline II helix, while segment 2 of ArkA (C-terminal) adopts a 310 helix, but is far less structured than segment 1. As segment 2 interacts with AbpSH3, it becomes more structured, but retains flexibility even in the fully engaged state. Binding simulations reveal that ArkA enters a flexible encounter complex before forming the fully engaged bound complex. In the encounter complex, transient nonspecific hydrophobic and long-range electrostatic contacts form between ArkA and the binding surface of SH3. The encounter complex ensemble includes conformations with segment 1 in both the forward and reverse orientation, suggesting that segment 2 may play a role in stabilizing the correct binding orientation. While the encounter complex forms quickly, the slow step of binding is the transition from the disordered encounter ensemble to the fully engaged state. In this transition, ArkA makes specific contacts with AbpSH3 and buries more hydrophobic surface. Simulating the binding between ApbSH3 and ArkA provides insight into the role of encounter complex intermediates and nonnative hydrophobic interactions for other SH3 domains and IDPs in general.

Entropy and Information within Intrinsically Disordered Protein Regions

Bioinformatics and biophysical studies of intrinsically disordered proteins and regions (IDRs) note the high entropy at individual sequence positions and in conformations sampled in solution. This prevents application of the canonical sequence-structure-function paradigm to IDRs and motivates the development of new methods to extract information from IDR sequences. We argue that the information in IDR sequences cannot be fully revealed through positional conservation, which largely measures stable structural contacts and interaction motifs. Instead, considerations of evolutionary conservation of molecular features can reveal the full extent of information in IDRs. Experimental quantification of the large conformational entropy of IDRs is challenging but can be approximated through the extent of conformational sampling measured by a combination of NMR spectroscopy and lower-resolution structural biology techniques, which can be further interpreted with simulations. Conformational entropy and other biophysical features can be modulated by post-translational modifications that provide functional advantages to IDRs by tuning their energy landscapes and enabling a variety of functional interactions and modes of regulation. The diverse mosaic of functional states of IDRs and their conformational features within complexes demands novel metrics of information, which will reflect the complicated sequence-conformational ensemble-function relationship of IDRs.

Mapping the allosteric network within a SH3 domain

SH3 domains are very abundant protein-protein interactions modules, involved in the regulation of several cellular processes. Whilst they have been associated to allosteric communication pathways between contiguous domains in multi-domain proteins, there is lack of information regarding the intra-domain allosteric cross-talk within the SH3 moiety. Here we scrutinize the presence of an allosteric network in the C-terminal SH3 domain of Grb2 protein, upon binding the Grb2-associated binding 2 protein. To explore allostery, we performed double mutant cycle analysis, a powerful quantitative approach based on mutagenesis in conjunction with kinetic experiments. Data reveal the presence of an unexpected allosteric sparse network that modulates the affinity between the SH3 domain and its physiological partner.

Intrinsically disordered proteins and phenotypic switching: Implications in cancer

It is now well established that intrinsically disordered proteins (IDPs) that constitute a large part of the proteome across the three kingdoms, play critical roles in several biological processes including phenotypic switching. However, dysregulated expression of IDPs that engage in promiscuous interactions can lead to pathological states. In this chapter, using cancer as a paradigm, we discuss how IDP conformational dynamics and the resultant conformational noise can modulate phenotypic switching. Thus, contrary to the prevailing wisdom that phenotypic switching is highly deterministic (has a genetic underpinning) in cancer, emerging evidence suggests that non-genetic mechanisms, at least in part due to the conformational noise, may also be a confounding factor in phenotypic switching.

Molecular Dynamics Simulations Combined with Nuclear Magnetic Resonance and/or Small-Angle X-ray Scattering Data for Characterizing Intrinsically Disordered Protein Conformational Ensembles

The concept of intrinsically disordered proteins (IDPs) has emerged relatively slowly, but over the past 20 years, it has become an intense research area in structural biology. Indeed, because of their considerable flexibility and structural heterogeneity, the determination of IDP conformational ensemble is particularly challenging and often requires a combination of experimental measurements and computational approaches. With the improved accuracy of all-atom force fields and the increasing computing performances, molecular dynamics (MD) simulations have become more and more reliable to generate realistic conformational ensembles. And the combination of MD simulations with experimental approaches, such as nuclear magnetic resonance (NMR) and/or small-angle X-ray scattering (SAXS) allows one to converge toward a more accurate and exhaustive description of IDP structures. In this Review, we discuss the state of the art of MD simulations of IDP conformational ensembles, with a special focus on studies that back-calculated and directly compared theoretical and experimental NMR or SAXS observables, such as chemical shifts (CS), ³J-couplings (³Jc), residual dipolar couplings (RDC), or SAXS intensities. We organize the review in three parts. In the first section, we discuss the studies which used NMR and/or SAXS data to test and validate the development of force fields or enhanced sampling techniques. In the second part, we explore different methods for the refinement of MD-derived structural ensembles, such as NMR or SAXS data-restrained MD simulations or ensemble reweighting to better fit experiments. Finally, we survey some recent studies combining MD simulations with NMR and/or SAXS measurements to investigate the relationship between IDP conformational ensemble and biological activity, as well as their implication in human diseases. From this review, we noticed that quite a few studies compared MD-generated conformational ensembles with both NMR and SAXS measurements to validate IDP structures at both local and global levels. Yet, beside the IDP propensity to form local secondary structures, their dynamic extension or compactness also appears important for their activity. Thus, we believe that a close synergy between MD simulations, NMR, and SAXS experiments would be greatly appropriate to address the challenges of characterizing the disordered structures of proteins and their complexes, relative to their biological functions.

The Singular NMR Fingerprint of a Polyproline II Helical Bundle

Polyproline II (PPII) helices play vital roles in biochemical recognition events and structures like collagen and form part of the conformational landscapes of intrinsically disordered proteins (IDPs). Nevertheless, this structure is generally hard to detect and quantify. Here, we report the first thorough NMR characterization of a PPII helical bundle protein, the Hypogastrura harveyi "snow flea" antifreeze protein (sfAFP). J-couplings and nuclear Overhauser enhancement spectroscopy confirm a natively folded structure consisting of six PPII helices. NMR spectral analyses reveal quite distinct Hα2 versus Hα3 chemical shifts for 28 Gly residues as well as ¹³Cα, ¹⁵N, and ¹HN conformational chemical shifts (Δδ) unique to PPII helical bundles. The ¹⁵N Δδ and ¹HN Δδ values and small negative ¹HN temperature coefficients evince hydrogen-bond formation. ¹H-¹⁵N relaxation measurements reveal that the backbone structure is generally highly rigid on ps-ns time scales. NMR relaxation parameters and biophysical characterization reveal that sfAFP is chiefly a dimer. For it, a structural model featuring the packing of long, flat hydrophobic faces at the dimer interface is advanced. The conformational stability, measured by amide H/D exchange to be 6.24 ± 0.2 kcal·mol^-1, is elevated. These are extraordinary findings considering the great entropic cost of fixing Gly residues and, together with the remarkable upfield chemical shifts of 28 Gly ¹Hα, evidence significant stabilizing contributions from CαHα ||| O═C hydrogen bonds. These stabilizing interactions are corroborated by density functional theory calculations and natural bonding orbital analysis. The singular conformational chemical shifts, J-couplings, high hNOE ratios, small negative temperature coefficients, and slowed H/D exchange constitute a unique set of fingerprints to identify PPII helical bundles, which may be formed by hundreds of Gly-rich motifs detected in sequence databases. These results should aid the quantification of PPII helices in IDPs, the development of improved antifreeze proteins, and the incorporation of PPII helices into novel designed proteins.

Affinity and specificity of motif-based protein-protein interactions

It is becoming increasingly clear that eukaryotic cell physiology is largely controlled by protein-protein interactions involving disordered protein regions, which usually interact with globular domains in a coupled binding and folding reaction. Several protein recognition domains are part of large families where members can interact with similar peptide ligands. Because of this, much research has been devoted to understanding how specificity can be achieved. A combination of interface complementarity, interactions outside of the core binding site, avidity from multidomain architecture and spatial and temporal regulation of expression resolves the conundrum. Here, we review recent advances in molecular aspects of affinity and specificity in such protein-protein interactions.

Molecular dynamics ensemble refinement of the heterogeneous native state of NCBD using chemical shifts and NOEs

Many proteins display complex dynamical properties that are often intimately linked to their biological functions. As the native state of a protein is best described as an ensemble of conformations, it is important to be able to generate models of native state ensembles with high accuracy. Due to limitations in sampling efficiency and force field accuracy it is, however, challenging to obtain accurate ensembles of protein conformations by the use of molecular simulations alone. Here we show that dynamic ensemble refinement, which combines an accurate atomistic force field with commonly available nuclear magnetic resonance (NMR) chemical shifts and NOEs, can provide a detailed and accurate description of the conformational ensemble of the native state of a highly dynamic protein. As both NOEs and chemical shifts are averaged on timescales up to milliseconds, the resulting ensembles reflect the structural heterogeneity that goes beyond that probed, e.g., by NMR relaxation order parameters. We selected the small protein domain NCBD as object of our study since this protein, which has been characterized experimentally in substantial detail, displays a rich and complex dynamical behaviour. In particular, the protein has been described as having a molten-globule like structure, but with a relatively rigid core. Our approach allowed us to describe the conformational dynamics of NCBD in solution, and to probe the structural heterogeneity resulting from both short- and long-timescale dynamics by the calculation of order parameters on different time scales. These results illustrate the usefulness of our approach since they show that NCBD is rather rigid on the nanosecond timescale, but interconverts within a broader ensemble on longer timescales, thus enabling the derivation of a coherent set of conclusions from various NMR experiments on this protein, which could otherwise appear in contradiction with each other.

Structural Basis for Natural Product Selection and Export by Bacterial ABC Transporters

Bacteria under stress produce ribosomally synthesized and post-translationally modified peptides (RiPPs) to target closely related species, such as the lasso peptide microcin J25 (MccJ25). These peptides are also toxic to the producing organisms that utilize dedicated ABC transporters to achieve self-immunity. MccJ25 is exported by the Escherichia coli ABC transporter McjD through a complex mechanism of recognition that has remained elusive. Here, we used biomolecular NMR to study this interaction and identified a region of the toxic peptide that is crucial to its recognition by the ABC transporter. Our study provides evidence that McjD is highly specific to MccJ25 and not to other RiPPs or antibiotics, unlike multidrug ABC transporters. Additionally, we show that MccJ25 is not exported by another natural product ABC transporter. Therefore, we propose that specific interactions between natural product ABC transporters and their substrate provides them with their high degree of specificity. Taken together, these findings suggest that ABC transporters might have acquired structural elements in their binding cavity to recognize and allow promiscuous export of a larger variety of compounds.

Characterizing the conformational landscape of MDM2-binding p53 peptides using Molecular Dynamics simulations

The conformational landscapes of p53 peptide variants and phage derived peptide (12/1) variants, all known to bind to MDM2, are studied using hamiltonian replica exchange molecular dynamics simulations. Complementing earlier observations, the current study suggests that the p53 peptides largely follow the ‘conformational selection’ paradigm in their recognition of and complexation by MDM2 while the 12/1 peptides likely undergo some element of conformational selection but are mostly driven by ‘binding induced folding’. This hypothesis is further supported by pulling simulations that pull the peptides away from their bound states with MDM2. This data extends the earlier mechanisms proposed to rationalize the entropically driven binding of the p53 set and the enthalpically driven binding of the 12/1 set. Using our hypothesis, we suggest mutations to the 12/1 peptide that increase its helicity in simulations and may, in turn, shift the binding towards conformational selection. In summary, understanding the conformational landscapes of the MDM2-binding peptides may suggest new peptide designs with bespoke binding mechanisms.

The PROSECCO server for chemical shift predictions in ordered and disordered proteins

The chemical shifts measured in solution-state and solid-state nuclear magnetic resonance (NMR) are powerful probes of the structure and dynamics of protein molecules. The exploitation of chemical shifts requires methods to correlate these data with the protein structures and sequences. We present here an approach to calculate accurate chemical shifts in both ordered and disordered proteins using exclusively the information contained in their sequences. Our sequence-based approach, protein sequences and chemical shift correlations (PROSECCO), achieves the accuracy of the most advanced structure-based methods in the characterization of chemical shifts of folded proteins and improves the state of the art in the study of disordered proteins. Our analyses revealed fundamental insights on the structural information carried by NMR chemical shifts of structured and unstructured protein states.

Understanding the mechanism of binding between Gab2 and the C terminal SH3 domain from Grb2

Gab2 is a large disordered protein that regulates several cellular signalling pathways and is overexpressed in different forms of cancer. Because of its disordered nature, a detailed characterization of the mechanisms of recognition between Gab2 and its physiological partners is particularly difficult. Here we provide a detailed kinetic characterization of the binding reaction between Gab2 and the C-terminal SH3 domain of the growth factor receptor-bound protein 2 (Grb2). We demonstrate that Gab2 folds upon binding following an induced fit type mechanism, whereby recognition is characterized by the formation of an intermediate, in which Gab2 is primarily disordered. In this scenario, folding of Gab2 into the bound conformation occurs only after binding. However, an alanine scanning of the proline residues of Gab2 suggests that the intermediate contains some degree of native-like structure, which might play a role for the recognition event to take place. The results, which represent a fundamental step forward in the understanding of this functional protein-protein interaction, are discussed on the light of previous structural works on these proteins.

Docking and molecular dynamics simulations of the Fyn-SH3 domain with free and phospholipid bilayer-associated 18.5-kDa myelin basic protein (MBP)-Insights into a noncanonical and fuzzy interaction

The molecular details of the association between the human Fyn-SH3 domain, and the fragment of 18.5-kDa myelin basic protein (MBP) spanning residues S38-S107 (denoted as xα2-peptide, murine sequence numbering), were studied in silico via docking and molecular dynamics over 50-ns trajectories. The results show that interaction between the two proteins is energetically favorable and heavily dependent on the MBP proline-rich region (P93-P98) in both aqueous and membrane environments. In aqueous conditions, the xα2-peptide/Fyn-SH3 complex adopts a "sandwich""-like structure. In the membrane context, the xα2-peptide interacts with the Fyn-SH3 domain via the proline-rich region and the β-sheets of Fyn-SH3, with the latter wrapping around the proline-rich region in a form of a clip. Moreover, the simulations corroborate prior experimental evidence of the importance of upstream segments beyond the canonical SH3-ligand. This study thus provides a more-detailed glimpse into the context-dependent interaction dynamics and importance of the β-sheets in Fyn-SH3 and proline-rich region of MBP. Proteins 2017; 85:1336-1350. © 2017 Wiley Periodicals, Inc.

Conserved Helix-Flanking Prolines Modulate Intrinsically Disordered Protein:Target Affinity by Altering the Lifetime of the Bound Complex

Appropriate integration of cellular signals requires a delicate balance of ligand–target binding affinities. Increasing the level of residual structure in intrinsically disordered proteins (IDPs), which are overrepresented in these cellular processes, has been shown previously to enhance binding affinities and alter cellular function. Conserved proline residues are commonly found flanking regions of IDPs that become helical upon interacting with a partner protein. Here, we mutate these helix-flanking prolines in p53 and MLL and find opposite effects on binding affinity upon an increase in free IDP helicity. In both cases, changes in affinity were due to alterations in dissociation, not association, rate constants, which is inconsistent with conformational selection mechanisms. We conclude that, contrary to previous suggestions, helix-flanking prolines do not regulate affinity by modulating the rate of complex formation. Instead, they influence binding affinities by controlling the lifetime of the bound complex.

Quarterly intrinsic disorder digest (April-May-June, 2014)

This is the 6^th issue of the Digested Disorder series that continues to use only 2 criteria for inclusion of a paper to this digest: The publication date (a paper should be published within the covered time frame) and the topic (a paper should be dedicated to any aspect of protein intrinsic disorder). The current digest issue covers papers published during the second quarter of 2014; i.e., during the period of April, May, and June of 2014. Similar to previous issues, the papers are grouped hierarchically by topics they cover, and for each of the included papers a short description is given on its major findings.

Bayesian inference of protein ensembles from SAXS data

The inherent flexibility of intrinsically disordered proteins (IDPs) and multi-domain proteins with intrinsically disordered regions (IDRs) presents challenges to structural analysis. These macromolecules need to be represented by an ensemble of conformations, rather than a single structure. Small-angle X-ray scattering (SAXS) experiments capture ensemble-averaged data for the set of conformations. We present a Bayesian approach to ensemble inference from SAXS data, called Bayesian ensemble SAXS (BE-SAXS). We address two issues with existing methods: the use of a finite ensemble of structures to represent the underlying distribution, and the selection of that ensemble as a subset of an initial pool of structures. This is achieved through the formulation of a Bayesian posterior of the conformational space. BE-SAXS modifies a structural prior distribution in accordance with the experimental data. It uses multi-step expectation maximization, with alternating rounds of Markov-chain Monte Carlo simulation and empirical Bayes optimization. We demonstrate the method by employing it to obtain a conformational ensemble of the antitoxin PaaA2 and comparing the results to a published ensemble.

Binding Mechanisms of Intrinsically Disordered Proteins: Theory, Simulation, and Experiment

In recent years, protein science has been revolutionized by the discovery of intrinsically disordered proteins (IDPs). In contrast to the classical paradigm that a given protein sequence corresponds to a defined structure and an associated function, we now know that proteins can be functional in the absence of a stable three-dimensional structure. In many cases, disordered proteins or protein regions become structured, at least locally, upon interacting with their physiological partners. Many, sometimes conflicting, hypotheses have been put forward regarding the interaction mechanisms of IDPs and the potential advantages of disorder for protein-protein interactions. Whether disorder may increase, as proposed, e.g., in the “fly-casting” hypothesis, or decrease binding rates, increase or decrease binding specificity, or what role pre-formed structure might play in interactions involving IDPs (conformational selection vs. induced fit), are subjects of intense debate. Experimentally, these questions remain difficult to address. Here, we review experimental studies of binding mechanisms of IDPs using NMR spectroscopy and transient kinetic techniques, as well as the underlying theoretical concepts and numerical methods that can be applied to describe these interactions at the atomic level. The available literature suggests that the kinetic and thermodynamic parameters characterizing interactions involving IDPs can vary widely and that there may be no single common mechanism that can explain the different binding modes observed experimentally. Rather, disordered proteins appear to make combined use of features such as pre-formed structure and flexibility, depending on the individual system and the functional context.

Structural Ensembles of Membrane-bound α-Synuclein Reveal the Molecular Determinants of Synaptic Vesicle Affinity

A detailed characterisation of the molecular determinants of membrane binding by α-synuclein (αS), a 140-residue protein whose aggregation is associated with Parkinson's disease, is of fundamental significance to clarify the manner in which the balance between functional and dysfunctional processes are regulated for this protein. Despite its biological relevance, the structural nature of the membrane-bound state αS remains elusive, in part because of the intrinsically dynamic nature of the protein and also because of the difficulties in studying this state in a physiologically relevant environment. In the present study we have used solid-state NMR and restrained MD simulations to refine structure and topology of the N-terminal region of αS bound to the surface of synaptic-like membranes. This region has fundamental importance in the binding mechanism of αS as it acts as to anchor the protein to lipid bilayers. The results enabled the identification of the key elements for the biological properties of αS in its membrane-bound state.

Insights into Coupled Folding and Binding Mechanisms from Kinetic Studies

Intrinsically disordered proteins (IDPs) are characterized by a lack of persistent structure. Since their identification more than a decade ago, many questions regarding their functional relevance and interaction mechanisms remain unanswered. Although most experiments have taken equilibrium and structural perspectives, fewer studies have investigated the kinetics of their interactions. Here we review and highlight the type of information that can be gained from kinetic studies. In particular, we show how kinetic studies of coupled folding and binding reactions, an important class of signaling event, are needed to determine mechanisms.

Sequence- and Temperature-Dependent Properties of Unfolded and Disordered Proteins from Atomistic Simulations

We use all-atom molecular simulation with explicit solvent to study the properties of selected intrinsically disordered proteins and unfolded states of foldable proteins, which include chain dimensions and shape, secondary structure propensity, solvent accessible surface area, and contact formation. We find that the qualitative scaling behavior of the chains matches expectations from theory under ambient conditions. In particular, unfolded globular proteins tend to be more collapsed under the same conditions than charged disordered sequences of the same length. However, inclusion of explicit solvent in addition naturally captures temperature-dependent solvation effects, which results in an initial collapse of the chains as temperature is increased, in qualitative agreement with experiment. There is a universal origin to the collapse, revealed in the change of hydration of individual residues as a function of temperature: namely, that the initial collapse is driven by unfavorable solvation free energy of individual residues, which in turn has a strong temperature dependence. We also observe that in unfolded globular proteins, increased temperature also initially favors formation of native-like (rather than non-native-like) structure. Our results help to establish how sequence encodes the degree of intrinsic disorder or order as well as its response to changes in environmental conditions.

Differential Recognition Preferences of the Three Src Homology 3 (SH3) Domains from the Adaptor CD2-associated Protein (CD2AP) and Direct Association with Ras and Rab Interactor 3 (RIN3)

CD2AP is an adaptor protein involved in membrane trafficking, with essential roles in maintaining podocyte function within the kidney glomerulus. CD2AP contains three Src homology 3 (SH3) domains that mediate multiple protein-protein interactions. However, a detailed comparison of the molecular binding preferences of each SH3 remained unexplored, as well as the discovery of novel interactors. Thus, we studied the binding properties of each SH3 domain to the known interactor Casitas B-lineage lymphoma protein (c-CBL), conducted a peptide array screen based on the recognition motif PxPxPR and identified 40 known or novel candidate binding proteins, such as RIN3, a RAB5-activating guanine nucleotide exchange factor. CD2AP SH3 domains 1 and 2 generally bound with similar characteristics and specificities, whereas the SH3-3 domain bound more weakly to most peptide ligands tested yet recognized an unusually extended sequence in ALG-2-interacting protein X (ALIX). RIN3 peptide scanning arrays revealed two CD2AP binding sites, recognized by all three SH3 domains, but SH3-3 appeared non-functional in precipitation experiments. RIN3 recruited CD2AP to RAB5a-positive early endosomes via these interaction sites. Permutation arrays and isothermal titration calorimetry data showed that the preferred binding motif is Px(P/A)xPR. Two high-resolution crystal structures (1.65 and 1.11 Å) of CD2AP SH3-1 and SH3-2 solved in complex with RIN3 epitopes 1 and 2, respectively, indicated that another extended motif is relevant in epitope 2. In conclusion, we have discovered novel interaction candidates for CD2AP and characterized subtle yet significant differences in the recognition preferences of its three SH3 domains for c-CBL, ALIX, and RIN3.

Computational approaches for inferring the functions of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are ubiquitously involved in cellular processes and often implicated in human pathological conditions. The critical biological roles of these proteins, despite not adopting a well-defined fold, encouraged structural biologists to revisit their views on the protein structure-function paradigm. Unfortunately, investigating the characteristics and describing the structural behavior of IDPs is far from trivial, and inferring the function(s) of a disordered protein region remains a major challenge. Computational methods have proven particularly relevant for studying IDPs: on the sequence level their dependence on distinct characteristics determined by the local amino acid context makes sequence-based prediction algorithms viable and reliable tools for large scale analyses, while on the structure level the in silico integration of fundamentally different experimental data types is essential to describe the behavior of a flexible protein chain. Here, we offer an overview of the latest developments and computational techniques that aim to uncover how protein function is connected to intrinsic disorder.

Globular and disordered-the non-identical twins in protein-protein interactions

In biology proteins from different structural classes interact across and within classes in ways that are optimized to achieve balanced functional outputs. The interactions between intrinsically disordered proteins (IDPs) and other proteins rely on changes in flexibility and this is seen as a strong determinant for their function. This has fostered the notion that IDP's bind with low affinity but high specificity. Here we have analyzed available detailed thermodynamic data for protein-protein interactions to put to the test if the thermodynamic profiles of IDP interactions differ from those of other protein-protein interactions. We find that ordered proteins and the disordered ones act as non-identical twins operating by similar principles but where the disordered proteins complexes are on average less stable by 2.5 kcal mol(-1).

Kirkwood-Buff Approach Rescues Overcollapse of a Disordered Protein in Canonical Protein Force Fields

Understanding the function of intrinsically disordered proteins is intimately related to our capacity to correctly sample their conformational dynamics. So far, a gap between experimentally and computationally derived ensembles exists, as simulations show overcompacted conformers. Increasing evidence suggests that the solvent plays a crucial role in shaping the ensembles of intrinsically disordered proteins and has led to several attempts to modify water parameters and thereby favor protein-water over protein-protein interactions. This study tackles the problem from a different perspective, which is the use of the Kirkwood-Buff theory of solutions to reproduce the correct conformational ensemble of intrinsically disordered proteins (IDPs). A protein force field recently developed on such a basis was found to be highly effective in reproducing ensembles for a fragment from the FG-rich nucleoporin 153, with dimensions matching experimental values obtained from small-angle X-ray scattering and single molecule FRET experiments. Kirkwood-Buff theory presents a complementary and fundamentally different approach to the recently developed four-site TIP4P-D water model, both of which can rescue the overcollapse observed in IDPs with canonical protein force fields. As such, our study provides a new route for tackling the deficiencies of current protein force fields in describing protein solvation.

Information flow and protein dynamics: the interplay between nuclear magnetic resonance spectroscopy and molecular dynamics simulations

Proteins participate in information pathways in cells, both as links in the chain of signals, and as the ultimate effectors. Upon ligand binding, proteins undergo conformation and motion changes, which can be sensed by the following link in the chain of information. Nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations represent powerful tools for examining the time-dependent function of biological molecules. The recent advances in NMR and the availability of faster computers have opened the door to more detailed analyses of structure, dynamics, and interactions. Here we briefly describe the recent applications that allow NMR spectroscopy and MD simulations to offer unique insight into the basic motions that underlie information transfer within and between cells.

Quantitative biophysical characterization of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are broadly defined as protein regions that do not cooperatively fold into a spatially or temporally stable structure. Recent research strongly supports the hypothesis that a conserved functional role for structural disorder renders IDPs uniquely capable of functioning in biological processes such as cellular signaling and transcription. Recently, the frequency of application of rigorous mechanistic biochemistry and quantitative biophysics to disordered systems has increased dramatically. For example, the launch of the Protein Ensemble Database (pE-DB) demonstrates that the potential now exists to refine models for the native state structure of IDPs using experimental data. However, rigorous assessment of which observables place the strongest and least biased constraints on those ensembles is now needed. Most importantly, the past few years have seen strong growth in the number of biochemical and biophysical studies attempting to connect structural disorder with function. From the perspective of equilibrium thermodynamics, there is a clear need to assess the relative significance of hydrophobic versus electrostatic forces in IDP interactions, if it is possible to generalize at all. Finally, kinetic mechanisms that invoke conformational selection and/or induced fit are often used to characterize coupled IDP folding and binding, although application of these models is typically built upon thermodynamic observations. Recently, the reaction rates and kinetic mechanisms of more intrinsically disordered systems have been tested through rigorous kinetic experiments. Motivated by these exciting advances, here we provide a review and prospectus for the quantitative study of IDP structure, thermodynamics, and kinetics.

The proline-rich region of 18.5 kDa myelin basic protein binds to the SH3-domain of Fyn tyrosine kinase with the aid of an upstream segment to form a dynamic complex in vitro

The intrinsically disordered 18.5 kDa classic isoform of MBP (myelin basic protein) interacts with Fyn kinase during oligodendrocyte development and myelination. It does so primarily via a central proline-rich SH3 (Src homology 3) ligand (T92–R104, murine 18.5 kDa MBP sequence numbering) that is part of a molecular switch due to its high degree of conservation and modification by MAP (mitogen-activated protein) and other kinases, especially at residues T92 and T95. Here, we show using co-transfection experiments of an early developmental oligodendroglial cell line (N19) that an MBP segment upstream of the primary ligand is involved in MBP–Fyn–SH3 association in cellula. Using solution NMR spectroscopy in vitro, we define this segment to comprise MBP residues (T62–L68), and demonstrate further that residues (V83–P93) are the predominant SH3-target, assessed by the degree of chemical shift change upon titration. We show by chemical shift index analysis that there is no formation of local poly-proline type II structure in the proline-rich segment upon binding, and by NOE (nuclear Overhauser effect) and relaxation measurements that MBP remains dynamic even while complexed with Fyn–SH3. The association is a new example first of a non-canonical SH3-domain interaction and second of a fuzzy MBP complex.

MBP interacts with Fyn kinase during oligodendrocyte development and myelination. We show that there is no binding-induced PPII formation in the primary ligand segment, and that a region upstream is required for in vitro interaction.

Disordered structural ensembles of vasopressin and oxytocin and their mutants

Vasopressin and oxytocin are intrinsically disordered cyclic nonapeptides belonging to a family of neurohypophysial hormones. Although unique in their functions, these peptides differ only by two residues and both feature a tocin ring formed by the disulfide bridge between first and sixth cysteine residues. This sequence and structural similarity are experimentally linked to oxytocin agonism at vasopressin receptors and vasopressin antagonism at oxytocin receptors. Yet single- or double-residue mutations in both peptides have been shown to have drastic impacts on their activities at either receptor, and possibly the ability to bind to their neurophysin carrier protein. In this study we perform molecular dynamics simulations of the unbound native and mutant sequences of the oxytocin and vasopressin hormones to characterize their structural ensembles. We classify the subpopulations of these structural ensembles on the basis of the distributions of radius of gyration and secondary structure and hydrogen-bonding features of the canonical tocin ring and disordered tail region. We then relate the structural changes observed in the unbound form of the different hormone sequences to experimental information about peptide receptor binding, and more indirectly, carrier protein binding affinity, receptor activity, and protease degradation. This study supports the hypothesis that the structural characteristics of the unbound form of an IDP can be used to predict structural or functional preferences of its functional bound form.