p15PAF is an intrinsically disordered protein with nonrandom structural preferences at sites of interaction with other proteins

Rescaling protein-protein interactions improves Martini 3 for flexible proteins in solution

Multidomain proteins with flexible linkers and disordered regions play important roles in many cellular processes, but characterizing their conformational ensembles is difficult. We have previously shown that the coarse-grained model, Martini 3, produces too compact ensembles in solution, that may in part be remedied by strengthening protein-water interactions. Here, we show that decreasing the strength of protein-protein interactions leads to improved agreement with experimental data on a wide set of systems. We show that the 'symmetry' between rescaling protein-water and protein-protein interactions breaks down when studying interactions with or within membranes; rescaling protein-protein interactions better preserves the binding specificity of proteins with lipid membranes, whereas rescaling protein-water interactions preserves oligomerization of transmembrane helices. We conclude that decreasing the strength of protein-protein interactions improves the accuracy of Martini 3 for IDPs and multidomain proteins, both in solution and in the presence of a lipid membrane.

A structure-based designed small molecule depletes hRpn13Pru and a select group of KEN box proteins

Proteasome subunit hRpn13 is partially proteolyzed in certain cancer cell types to generate hRpn13Pru by degradation of its UCHL5/Uch37-binding DEUBAD domain and retention of an intact proteasome- and ubiquitin-binding Pru domain. By using structure-guided virtual screening, we identify an hRpn13 binder (XL44) and solve its structure ligated to hRpn13 Pru by integrated X-ray crystallography and NMR to reveal its targeting mechanism. Surprisingly, hRpn13Pru is depleted in myeloma cells following treatment with XL44. TMT-MS experiments reveal a select group of off-targets, including PCNA clamp-associated factor PCLAF and ribonucleoside-diphosphate reductase subunit M2 (RRM2), that are similarly depleted by XL44 treatment. XL44 induces hRpn13-dependent apoptosis and also restricts cell viability by a PCLAF-dependent mechanism. A KEN box, but not ubiquitination, is required for XL44-induced depletion of PCLAF. Here, we show that XL44 induces ubiquitin-dependent loss of hRpn13Pru and ubiquitin-independent loss of select KEN box containing proteins.

The pathogen-encoded signalling receptor Tir exploits host-like intrinsic disorder for infection

The translocated intimin receptor (Tir) is an essential type III secretion system (T3SS) effector of attaching and effacing pathogens contributing to the global foodborne disease burden. Tir acts as a cell-surface receptor in host cells, rewiring intracellular processes by targeting multiple host proteins. We investigated the molecular basis for Tir's binding diversity in signalling, finding that Tir is a disordered protein with host-like binding motifs. Unexpectedly, also are several other T3SS effectors. By an integrative approach, we reveal that Tir dimerises via an antiparallel OB-fold within a highly disordered N-terminal cytosolic domain. Also, it has a long disordered C-terminal cytosolic domain partially structured at host-like motifs that bind lipids. Membrane affinity depends on lipid composition and phosphorylation, highlighting a previously unrecognised host interaction impacting Tir-induced actin polymerisation and cell death. Furthermore, multi-site tyrosine phosphorylation enables Tir to engage host SH2 domains in a multivalent fuzzy complex, consistent with Tir's scaffolding role and binding promiscuity. Our findings provide insights into the intracellular Tir domains, highlighting the ability of T3SS effectors to exploit host-like protein disorder as a strategy for host evasion.

Balanced Three-Point Water Model OPC3-B for Intrinsically Disordered and Ordered Proteins

Intrinsically disordered proteins (IDPs) play a critical role in many biological processes. Due to the inherent structural flexibility of IDPs, experimental methods present significant challenges for sampling their conformational information at the atomic level. Therefore, molecular dynamics (MD) simulations have emerged as the primary tools for modeling IDPs whose accuracy depend on force field and water model. To enhance the accuracy of physical modeling of IDPs, several force fields have been developed. However, current water models lack precision and underestimate the interaction between water molecules and proteins. Here, we used Monte-Carlo re-weighting method to re-parameterize a three-point water model based on OPC3 for IDPs (named OPC3-B). We benchmarked the performance of OPC3-B compared with nine different water models for 10 IDPs and three ordered proteins. The results indicate that the performance of OPC3-B is better than other water models for both IDPs and ordered proteins. At the same time, OPC3-B possess the power of transferability with other force field to simulate IDPs. This newly developed water model can be used to insight into the research of sequence-disordered-function paradigm for IDPs.

Improved predictions of phase behaviour of intrinsically disordered proteins by tuning the interaction range

The formation and viscoelastic properties of condensates of intrinsically disordered proteins (IDPs) is dictated by amino acid sequence and solution conditions. Because of the involvement of biomolecular condensates in cell physiology and disease, advancing our understanding of the relationship between protein sequence and phase separation (PS) may have important implications in the formulation of new therapeutic hypotheses. Here, we present CALVADOS 2, a coarse-grained model of IDPs that accurately predicts conformational properties and propensities to undergo PS for diverse sequences and solution conditions. In particular, we systematically study the effect of varying the range of the nonionic interactions and use our findings to improve the temperature scale of the model. We further optimize the residue-specific model parameters against experimental data on the conformational properties of 55 proteins, while also leveraging 70 hydrophobicity scales from the literature to avoid overfitting the training data. Extensive testing shows that the model accurately predicts chain compaction and PS propensity for sequences of diverse length and charge patterning, as well as at different temperatures and salt concentrations.

Improved predictions of phase behaviour of intrinsically disordered proteins by tuning the interaction range

The formation and viscoelastic properties of condensates of intrinsically disordered proteins (IDPs) is dictated by amino acid sequence and solution conditions. Because of the involvement of biomolecular condensates in cell physiology and disease, advancing our understanding of the relationship between protein sequence and phase separation (PS) may have important implications in the formulation of new therapeutic hypotheses. Here, we present CALVADOS 2, a coarse-grained model of IDPs that accurately predicts conformational properties and propensities to undergo PS for diverse sequences and solution conditions. In particular, we systematically study the effect of varying the range of the nonionic interactions and use our findings to improve the temperature scale of the model. We further optimize the residue-specific model parameters against experimental data on the conformational properties of 55 proteins, while also leveraging 70 hydrophobicity scales from the literature to avoid overfitting the training data. Extensive testing shows that the model accurately predicts chain compaction and PS propensity for sequences of diverse length and charge patterning, as well as at different temperatures and salt concentrations.

Performance evaluation of the balanced force field ff03CMAP for intrinsically disordered and ordered proteins

Intrinsically disordered proteins (IDPs) have been found to be closely associated with various human diseases. Because IDPs have no fixed tertiary structure under physiological conditions, current experimental methods, such as X-ray spectroscopy, NMR, and CryoEM, cannot capture all the dynamic conformations. Molecular dynamics simulation is an useful tool that is widely used to study the conformer distributions of IDPs and has become an important complementary tool for experimental methods. However, the accuracy of MD simulations directly depends on utilizing a precise force field. Recently a CMAP optimized force field based on the Amber ff03 force field (termed ff03CMAP herein) was developed for a balanced sampling of IDPs and folded proteins. In order to further evaluate the performance, more types of disordered and ordered proteins were used to test the ability for conformer sampling. The results showed that simulated chemical shifts, J-coupling, and R_g distribution with the ff03CMAP force field were in better agreement with NMR measurements and were more accurate than those with the ff03 force field. The sampling conformations by ff03CMAP were more diverse than those of ff03. At the same time, ff03CMAP could stabilize the conformers of the ordered proteins. These findings indicate that ff03CMAP can be widely used to sample diverse conformers for proteins, including the intrinsically disordered regions.

Improving Martini 3 for Disordered and Multidomain Proteins

Coarse-grained molecular dynamics simulations are a useful tool to determine conformational ensembles of proteins. Here, we show that the coarse-grained force field Martini 3 underestimates the global dimensions of intrinsically disordered proteins (IDPs) and multidomain proteins when compared with small-angle X-ray scattering (SAXS) data and that increasing the strength of protein-water interactions favors more expanded conformations. We find that increasing the strength of interactions between protein and water by ca. 10% results in improved agreement with the SAXS data for IDPs and multidomain proteins. We also show that this correction results in a more accurate description of self-association of IDPs and folded proteins and better agreement with paramagnetic relaxation enhancement data for most IDPs. While simulations with this revised force field still show deviations to experiments for some systems, our results suggest that it is overall a substantial improvement for coarse-grained simulations of soluble proteins.

Accurate model of liquid-liquid phase behavior of intrinsically disordered proteins from optimization of single-chain properties

Many intrinsically disordered proteins (IDPs) may undergo liquid-liquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalization of intracellular biochemical reactions. The phase behavior of IDPs is sequence dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intramolecular and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.

Accurate model of liquid-liquid phase behavior of intrinsically disordered proteins from optimization of single-chain properties

Many intrinsically disordered proteins (IDPs) may undergo liquid-liquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalization of intracellular biochemical reactions. The phase behavior of IDPs is sequence dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intramolecular and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.

Balanced Solvent Model for Intrinsically Disordered and Ordered Proteins

Intrinsically disordered proteins (IDPs) have no fixed three-dimensional (3D) structures under physiological conditions, with the content being about 51% in human proteomics. IDPs are associated with many human diseases, such as cancer, diabetes, and neurodegenerative diseases. Because IDPs do not crystallize and have diverse conformers, traditional experimental methods such as crystallization and NMR can hardly capture their conformation ensemble and just provide average structural characters of IDPs. Therefore, molecular dynamics (MD) simulations become a valuable complement to the experimental data. However, the accuracy of molecular dynamics simulation for IDPs depends on the combination of force fields and solvent models. Recently, we released an environment-specific force field (ESFF1) for IDPs, which can well reproduce the local structural properties (such as J-coupling and secondary chemical shifts). However, there is still a large deviation for the radius of gyration (R_g). Therefore, a solvent model combined with ESFF1 is necessary to capture the local and global characters for IDPs and ordered proteins. Here, we investigated the underestimation or overestimation of the solvent interaction for four solvent models (TIP3P, TIP4P-Ew, TIP4P-D, OPC) under ESFF1 and found the important ε parameter of the solvent model to play a key role in scaling R_g. A near-linear relationship between the simulation R_g and the ε parameter was used to develop the new solvent model, named TIP4P-B. The results indicate that the simulated R_g with TIP4P-B is in better agreement with the experimental observations than the other four solvent models. Simultaneously, TIP4P-B can also maintain the advantages of the ESFF1 force field for the local structural properties. Additionally, TIP4P-B can successfully sample the conformation of ordered proteins. These findings confirm that TIP4P-B is a balanced solvent model and can improve sampling R_g performance for folded proteins and IDPs.

SANS (USH1G) regulates pre-mRNA splicing by mediating the intra-nuclear transfer of tri-snRNP complexes

Splicing is catalyzed by the spliceosome, a compositionally dynamic complex assembled stepwise on pre-mRNA. We reveal links between splicing machinery components and the intrinsically disordered ciliopathy protein SANS. Pathogenic mutations in SANS/USH1G lead to Usher syndrome—the most common cause of deaf-blindness. Previously, SANS was shown to function only in the cytosol and primary cilia. Here, we have uncovered molecular links between SANS and pre-mRNA splicing catalyzed by the spliceosome in the nucleus. We show that SANS is found in Cajal bodies and nuclear speckles, where it interacts with components of spliceosomal sub-complexes such as SF3B1 and the large splicing cofactor SON but also with PRPFs and snRNAs related to the tri-snRNP complex. SANS is required for the transfer of tri-snRNPs between Cajal bodies and nuclear speckles for spliceosome assembly and may also participate in snRNP recycling back to Cajal bodies. SANS depletion alters the kinetics of spliceosome assembly, leading to accumulation of complex A. SANS deficiency and USH1G pathogenic mutations affects splicing of genes related to cell proliferation and human Usher syndrome. Thus, we provide the first evidence that splicing dysregulation may participate in the pathophysiology of Usher syndrome.

Extensive evaluation of environment-specific force field for ordered and disordered proteins

Intrinsically disordered proteins (IDPs) have no fixed tertiary structure under physiological conditions and are associated with many human diseases. Because IDPs have the characteristic of possessing diverse conformations, current experimental methods cannot capture all the conformations of IDPs. However, molecular dynamics simulation can sample these atomistically diverse conformations as a valuable complement to experimental data. To accurately describe the properties of IDPs, the environment-specific precise force field (ESFF1) was successfully released to reproduce the conformer character of ordered and disordered proteins. Here, three typical IDPs and thirteen folded proteins were used to further evaluate the performance of this force field. The results indicate that the NMR observables of ESFF1 better approach experimental data than do those of ff14SB for IDPs. The sampling conformations by ESFF1 are more diverse than those of ff14SB. For folded proteins, these force fields have comparable performances for reproducing conformers. Therefore, ESFF1 can be used to reveal the model of sequence-disorder-function for IDPs.

Comment on the Optimal Parameters to Derive Intrinsically Disordered Protein Conformational Ensembles from Small-Angle X-ray Scattering Data Using the Ensemble Optimization Method

The Ensemble Optimization Method (EOM) is a popular approach to describe small-angle X-ray scattering (SAXS) data from highly disordered proteins. The EOM algorithm selects subensembles of coexisting states from large pools of randomized conformations to fit the SAXS data. Based on the unphysical bimodal radius of gyration (R_g) distribution of conformations resulting from the EOM analysis, a recent article (Fagerberg et al. J. Chem. Theory Comput. 2019, 15 (12), 6968-6983) concluded that this approach inadequately described the SAXS data measured for human Histatin 5 (Hst5), a peptide with antifungal properties. Using extensive experimental and synthetic data, we explored the origin of this observation. We found that the one-bead-per-residue coarse-grained representation with averaged scattering form factors (provided in the EOM as an add-on to represent disordered missing loops or domains) may not be appropriate for EOM analyses of scattering data from short (below 50 residues) proteins/peptides. The method of choice for these proteins is to employ atomistic models (e.g., from molecular dynamics simulations) to sample the protein conformational landscape. As a convenient alternative, we have also improved the coarse-grained approach by introducing amino acid specific form factors in the calculations. We also found that, for small proteins, the search for relatively large subensembles of 20-50 conformers (as implemented in the original EOM version) more adequately describes the conformational space sampled in solution than the procedures optimizing the ensemble size. Our observations have been added as recommendations into the information for EOM users to promote the proper utilization of the program for ensemble-based modeling of SAXS data for all types of disordered systems.

MMM: Integrative ensemble modeling and ensemble analysis

Proteins and their complexes can be heterogeneously disordered. In ensemble modeling of such systems with restraints from several experimental techniques the following problems arise: (a) integration of diverse restraints obtained on different samples under different conditions; (b) estimation of a realistic ensemble width; (c) sufficient sampling of conformational space; (d) representation of the ensemble by an interpretable number of conformers; (e) recognition of weak order with site resolution. Here, I introduce several tools that address these problems, focusing on utilization of distance distribution information for estimating ensemble width. The RigiFlex approach integrates such information with high-resolution structures of ordered domains and small-angle scattering data. The EnsembleFit module provides moderately sized ensembles by fitting conformer populations and discarding conformers with low population. EnsembleFit balances the loss in fit quality upon combining restraint subsets from different techniques. Pair correlation analysis for residues and local compaction analysis help in feature detection. The RigiFlex pipeline is tested on data simulated from the structure 70 kDa protein-RNA complex RsmE/RsmZ. It recovers this structure with ensemble width and difference from ground truth both being on the order of 4.2 Å. EnsembleFit reduces the ensemble of the proliferating-cell-nuclear-antigen-associated factor p15PAF from 4,939 to 75 conformers while maintaining good fit quality of restraints. Local compaction analysis for the PaaA2 antitoxin from E. coli O157 revealed correlations between compactness and enhanced residual dipolar couplings in the original NMR restraint set.

Towards Decoding the Sequence-Based Grammar Governing the Functions of Intrinsically Disordered Protein Regions

A substantial portion of the proteome consists of intrinsically disordered regions (IDRs) that do not fold into well-defined 3D structures yet perform numerous biological functions and are associated with a broad range of diseases. It has been a long-standing enigma how different IDRs successfully execute their specific functions. Further putting a spotlight on IDRs are recent discoveries of functionally relevant biomolecular assemblies, which in some cases form through liquid-liquid phase separation. At the molecular level, the formation of biomolecular assemblies is largely driven by weak, multivalent, but selective IDR-IDR interactions. Emerging experimental and computational studies suggest that the primary amino acid sequences of IDRs encode a variety of their interaction behaviors. In this review, we focus on findings and insights that connect sequence-derived features of IDRs to their conformations, propensities to form biomolecular assemblies, selectivity of interaction partners, functions in the context of physiology and disease, and regulation of function. We also discuss directions of future research to facilitate establishing a comprehensive sequence-function paradigm that will eventually allow prediction of selective interactions and specificity of function mediated by IDRs.

Intrinsically Disordered Bacterial Polar Organizing Protein Z, PopZ, Interacts with Protein Binding Partners Through an N-terminal Molecular Recognition Feature

The polar organizing protein Z (PopZ) is necessary for the formation of three-dimensional microdomains at the cell poles in Caulobacter crescentus, where it functions as a hub protein that recruits multiple regulatory proteins from the cytoplasm. Although a large portion of the protein is predicted to be natively unstructured, in reconstituted systems PopZ can self-assemble into a macromolecular scaffold that directly binds to at least ten different proteins. Here we report the solution NMR structure of PopZ^Δ134-177, a truncated form of PopZ that does not self-assemble but retains the ability to interact with heterologous proteins. We show that the unbound form of PopZ^Δ134-177 is unstructured in solution, with the exception of a small amphipathic α-helix in residues M10-I17, which is included within a highly conserved region near the N-terminal. In applying NMR techniques to map the interactions between PopZ^Δ134-177 and one of its binding partners, RcdA, we find evidence that the α-helix and adjoining amino acids extending to position E23 serve as the core of the binding motif. Consistent with this, a point mutation at position I17 severely compromises binding. Our results show that a partially structured Molecular Recognition Feature (MoRF) within an intrinsically disordered domain of PopZ contributes to the assembly of polar microdomains, revealing a structural basis for complex network assembly in Alphaproteobacteria that is analogous to those formed by intrinsically disordered hub proteins in other kingdoms.

Analysis of Protein Disorder Predictions in the Light of a Protein Structural Alphabet

Intrinsically-disordered protein (IDP) characterization was an amazing change of paradigm in our classical sequence-structure-function theory. Moreover, IDPs are over-represented in major disease pathways and are now often targeted using small molecules for therapeutic purposes. This has had created a complex continuum from order-that encompasses rigid and flexible regions-to disorder regions; the latter being not accessible through classical crystallographic methodologies. In X-ray structures, the notion of order is dictated by access to resolved atom positions, providing rigidity and flexibility information with low and high experimental B-factors, while disorder is associated with the missing (non-resolved) residues. Nonetheless, some rigid regions can be found in disorder regions. Using ensembles of IDPs, their local conformations were analyzed in the light of a structural alphabet. An entropy index derived from this structural alphabet allowed us to propose a continuum of states from rigidity to flexibility and finally disorder. In this study, the analysis was extended to comparing these results to disorder predictions, underlying a limited correlation, and so opening new ideas to characterize and predict disorder.

Human PCNA Structure, Function and Interactions

Proliferating cell nuclear antigen (PCNA) is an essential factor in DNA replication and repair. It forms a homotrimeric ring that embraces the DNA and slides along it, anchoring DNA polymerases and other DNA editing enzymes. It also interacts with regulatory proteins through a sequence motif known as PCNA Interacting Protein box (PIP-box). We here review the latest contributions to knowledge regarding the structure-function relationships in human PCNA, particularly the mechanism of sliding, and of the molecular recognition of canonical and non-canonical PIP motifs. The unique binding mode of the oncogene p15 is described in detail, and the implications of the recently discovered structure of PCNA bound to polymerase δ are discussed. The study of the post-translational modifications of PCNA and its partners may yield therapeutic opportunities in cancer treatment, in addition to illuminating the way PCNA coordinates the dynamic exchange of its many partners in DNA replication and repair.

Multivalent interaction of ESCO2 with the replication machinery is required for sister chromatid cohesion in vertebrates

The tethering together of sister chromatids by the cohesin complex ensures their accurate alignment and segregation during cell division. In vertebrates, sister chromatid cohesion requires the activity of the ESCO2 acetyltransferase, which modifies the Smc3 subunit of cohesin. It was shown recently that ESCO2 promotes cohesion through interaction with the MCM replicative helicase. However, ESCO2 does not significantly colocalize with the MCM complex, suggesting there are additional interactions important for ESCO2 function. Here we show that ESCO2 is recruited to replication factories, sites of DNA replication, through interaction with PCNA. We show that ESCO2 contains multiple PCNA-interaction motifs in its N terminus, each of which is essential to its ability to establish cohesion. We propose that multiple PCNA-interaction motifs embedded in a largely flexible and disordered region of the protein underlie the unique ability of ESCO2 to establish cohesion between sister chromatids precisely as they are born during DNA replication.

The PCNA interaction motifs revisited: thinking outside the PIP-box

Proliferating cell nuclear antigen (PCNA) is a cellular hub in DNA metabolism and a potential drug target. Its binding partners carry a short linear motif (SLiM) known as the PCNA-interacting protein-box (PIP-box), but sequence-divergent motifs have been reported to bind to the same binding pocket. To investigate how PCNA accommodates motif diversity, we assembled a set of 77 experimentally confirmed PCNA-binding proteins and analyzed features underlying their binding affinity. Combining NMR spectroscopy, affinity measurements and computational analyses, we corroborate that most PCNA-binding motifs reside in intrinsically disordered regions, that structure preformation is unrelated to affinity, and that the sequence-patterns that encode binding affinity extend substantially beyond the boundaries of the PIP-box. Our systematic multidisciplinary approach expands current views on PCNA interactions and reveals that the PIP-box affinity can be modulated over four orders of magnitude by positive charges in the flanking regions. Including the flanking regions as part of the motif is expected to have broad implications, particularly for interpretation of disease-causing mutations and drug-design, targeting DNA-replication and -repair.

Double Monoubiquitination Modifies the Molecular Recognition Properties of p15PAF Promoting Binding to the Reader Module of Dnmt1

The proliferating cell nuclear antigen (PCNA)-associated factor p15^PAF is a nuclear protein that acts as a regulator of DNA repair during DNA replication. The p15^PAF gene is overexpressed in several types of human cancer, and its function is regulated by monoubiquitination of two lysines (K15 and K24) at the protein N-terminal region. We have previously shown that p15^PAF is an intrinsically disordered protein which partially folds upon binding to PCNA and independently contacts DNA through its N-terminal tail. Here we present an NMR conformational characterization of p15^PAF monoubiquitinated at both K15 and K24 via a disulfide bridge mimicking the isopeptide bond. We show that doubly monoubiquitinated p15^PAF is monomeric, intrinsically disordered, and binds to PCNA as nonubiquitinated p15^PAF does but interacts with DNA with reduced affinity. Our SAXS-derived conformational ensemble of doubly monoubiquitinated p15^PAF shows that the ubiquitin moieties, separated by eight disordered residues, form transient dimers because of the high local effective ubiquitin concentration. This observation and the sequence similarity with histone H3 N-terminal tail suggest that doubly monoubiquitinated p15^PAF is a binding target of DNA methyl transferase Dnmt1, as confirmed by calorimetry. Therefore, doubly monoubiquitinated p15^PAF directly interacts with PCNA and recruits Dnmt1 for maintenance of DNA methylation during replication.

Dynamics of the intrinsically disordered protein NUPR1 in isolation and in its fuzzy complexes with DNA and prothymosin α

Intrinsically disordered proteins (IDPs) explore diverse conformations in their free states and, a few of them, also in their molecular complexes. This functional plasticity is essential for the function of IDPs, although their dynamics in both free and bound states is poorly understood. NUPR1 is a protumoral multifunctional IDP, activated during the acute phases of pancreatitis. It interacts with DNA and other IDPs, such as prothymosin α (ProTα), with dissociation constants of ~0.5 μM, and a 1:1 stoichiometry. We studied the structure and picosecond-to-nanosecond (ps-ns) dynamics by using both NMR and SAXS in: (i) isolated NUPR1; (ii) the NUPR1/ProTα complex; and (iii) the NUPR1/double stranded (ds) GGGCGCGCCC complex. Our SAXS findings show that NUPR1 remained disordered when bound to either partner, adopting a worm-like conformation; the fuzziness of bound NUPR1 was also pinpointed by NMR. Residues with the largest values of the relaxation rates (R₁, R_1ρ, R₂ and η_xy), in the free and bound species, were mainly clustered around the 30s region of the sequence, which agree with one of the protein hot-spots already identified by site-directed mutagenesis. Not only residues in this region had larger relaxation rates, but they also moved slower than the rest of the molecule, as indicated by the reduced spectral density approach (RSDA). Upon binding, the energy landscape of NUPR1 was not funneled down to a specific, well-folded conformation, but rather its backbone flexibility was kept, with distinct motions occurring at the hot-spot region.

p15PAF binding to PCNA modulates the DNA sliding surface

p15^PAF is an oncogenic intrinsically disordered protein that regulates DNA replication and lesion bypass by interacting with the human sliding clamp PCNA. In the absence of DNA, p15^PAF traverses the PCNA ring via an extended PIP-box that contacts the sliding surface. Here, we probed the atomic-scale structure of p15^PAF–PCNA–DNA ternary complexes. Crystallography and MD simulations show that, when p15^PAF occupies two subunits of the PCNA homotrimer, DNA within the ring channel binds the unoccupied subunit. The structure of PCNA-bound p15^PAF in the absence and presence of DNA is invariant, and solution NMR confirms that DNA does not displace p15^PAF from the ring wall. Thus, p15^PAF reduces the available sliding surfaces of PCNA, and may function as a belt that fastens the DNA to the clamp during synthesis by the replicative polymerase (pol δ). This constraint, however, may need to be released for efficient DNA lesion bypass by the translesion synthesis polymerase (pol η). Accordingly, our biochemical data show that p15^PAF impairs primer synthesis by pol η–PCNA holoenzyme against both damaged and normal DNA templates. In light of our findings, we discuss the possible mechanistic roles of p15^PAF in DNA replication and suppression of DNA lesion bypass.

Effect of Weakly Interacting Cosolutes on Lysozyme Conformations

Exposure of a protein to cosolutes, like denaturants, changes its folding equilibrium. To determine the ensemble of protein conformations at equilibrium, in the presence of weakly interacting cosolutes, we present a two-stage analysis of solution X-ray scattering data. In the first stage, Guinier analysis and Kratky plot revealed information about the compactness and flexibility of the protein. In the second stage, elastic network contact model and coarse-grained normal mode analysis were used to generate an ensemble of conformations. The scattering curves of the conformations were computed and fitted to the measured scattering curves to get insights into the dominating folding states at equilibrium. Urea and guanidine hydrochloride (GuHCl) behaved as preferentially included weakly interacting cosolutes and induced denaturation of hen egg-white lysozyme, which served as our test case. The computed models adequately fit the data and gave ensembles of conformations that were consistent with our measurements. The analysis suggests that in the presence of urea, lysozyme retained its compactness and assumed molten globule characteristics, whereas in the presence of GuHCl lysozyme adopted random coiled conformations. Interestingly, no equilibrium intermediate states were observed in both urea and GuHCl.

Steric interference from intrinsically disordered regions controls dynamin-related protein 1 self-assembly during mitochondrial fission

The self-assembling, mechanoenzymatic dynamin superfamily GTPase, dynamin-related protein 1 (Drp1), catalyzes mitochondrial and peroxisomal fission. Distinct intrinsically disordered regions (IDRs) in Drp1 substitute for the canonical pleckstrin homology (PH) domain and proline-rich domain (PRD) of prototypical dynamin, which cooperatively regulate endocytic vesicle scission. Whether the Drp1 IDRs function analogously to the corresponding dynamin domains however remains unknown. We show that an IDR unique to the Drp1 GTPase (G) domain, the 'extended 80-loop', albeit dissimilar in location, structure, and mechanism, functions akin to the dynamin PRD by enabling stable Drp1 mitochondrial recruitment and by suppressing Drp1 cooperative GTPase activity in the absence of specific partner-protein interactions. Correspondingly, we find that another IDR, the Drp1 variable domain (VD), in conjunction with the conserved stalk L1N loop, functions akin to the dynamin PH domain; first, in an 'auto-inhibitory' capacity that restricts Drp1 activity through a long-range steric inhibition of helical inter-rung G-domain dimerization, and second, as a 'fulcrum' for Drp1 self-assembly in the proper helical register. We show that the Drp1 VD is necessary and sufficient for specific Drp1-phospholipid interactions. We further demonstrate that the membrane-dependent VD conformational rearrangement essential for the alleviation of Drp1 auto-inhibition is contingent upon the basal GTP hydrolysis-dependent generation of Drp1 dimers from oligomers in solution. IDRs thus conformationally couple the enzymatic and membrane activities of Drp1 toward membrane fission.

Developing a molecular dynamics force field for both folded and disordered protein states

Many proteins that perform important biological functions are completely or partially disordered under physiological conditions. Molecular dynamics simulations could be a powerful tool for the structural characterization of such proteins, but it has been unclear whether the physical models (force fields) used in simulations are sufficiently accurate. Here, we systematically compare the accuracy of a number of different force fields in simulations of both ordered and disordered proteins, finding that each force field has strengths and limitations. We then describe a force field that substantially improves on the state-of-the-art accuracy for simulations of disordered proteins without sacrificing accuracy for folded proteins, thus broadening the range of biological systems amenable to molecular dynamics simulations.

Molecular dynamics (MD) simulation is a valuable tool for characterizing the structural dynamics of folded proteins and should be similarly applicable to disordered proteins and proteins with both folded and disordered regions. It has been unclear, however, whether any physical model (force field) used in MD simulations accurately describes both folded and disordered proteins. Here, we select a benchmark set of 21 systems, including folded and disordered proteins, simulate these systems with six state-of-the-art force fields, and compare the results to over 9,000 available experimental data points. We find that none of the tested force fields simultaneously provided accurate descriptions of folded proteins, of the dimensions of disordered proteins, and of the secondary structure propensities of disordered proteins. Guided by simulation results on a subset of our benchmark, however, we modified parameters of one force field, achieving excellent agreement with experiment for disordered proteins, while maintaining state-of-the-art accuracy for folded proteins. The resulting force field, a99SB-disp, should thus greatly expand the range of biological systems amenable to MD simulation. A similar approach could be taken to improve other force fields.

Disentangling polydispersity in the PCNA-p15PAF complex, a disordered, transient and multivalent macromolecular assembly

The intrinsically disordered p15PAF regulates DNA replication and repair when interacting with the Proliferating Cell Nuclear Antigen (PCNA) sliding clamp. As many interactions between disordered proteins and globular partners involved in signaling and regulation, the complex between p15PAF and trimeric PCNA is of low affinity, forming a transient complex that is difficult to characterize at a structural level due to its inherent polydispersity. We have determined the structure, conformational fluctuations, and relative population of the five species that coexist in solution by combining small-angle X-ray scattering (SAXS) with molecular modelling. By using explicit ensemble descriptions for the individual species, built using integrative approaches and molecular dynamics (MD) simulations, we collectively interpreted multiple SAXS profiles as population-weighted thermodynamic mixtures. The analysis demonstrates that the N-terminus of p15PAF penetrates the PCNA ring and emerges on the back face. This observation substantiates the role of p15PAF as a drag regulating PCNA processivity during DNA repair. Our study reveals the power of ensemble-based approaches to decode structural, dynamic, and thermodynamic information from SAXS data. This strategy paves the way for deciphering the structural bases of flexible, transient and multivalent macromolecular assemblies involved in pivotal biological processes.

SAS profile correlations reveal SAS hierarchical nature and information content

In structural biology, Small-Angle Scattering experiments (SAS) are unique, because although they provide low resolution data, they can be performed in closer-to-native conditions than those arising in X-Ray crystallography. A number of questions on SAS, however, remain unsolved, particularly in the light of modelling ensembles of conformers in solution. In this article, we study the ensemble average and covariance of SAS profiles analytically. Using this ensemble covariance, we demonstrate the hierarchical nature of SAS profiles. Furthermore, we show that the information content is not uniform and reaches its maximum in the intermediate q range. The arguments are generalized using microsecond-scale molecular dynamics trajectories of the lysozyme and on an ensemble of the intrinsically disordered protein p15PAF. We show that for highly flexible systems, the SAS profile is a representation of the ensemble of conformers in solution, and not that of one conformer in particular.

The metastasis suppressor KISS1 is an intrinsically disordered protein slightly more extended than a random coil

The metastasis suppressor KISS1 is reported to be involved in the progression of several solid neoplasias, making it a promising molecular target for controlling their metastasis. The KISS1 sequence contains an N-terminal secretion signal and several dibasic sequences that are proposed to be the proteolytic cleavage sites. We present the first structural characterization of KISS1 by circular dichroism, multi-angle light scattering, small angle X-Ray scattering and NMR spectroscopy. An analysis of the KISS1 backbone NMR chemical shifts does not reveal any preferential conformation and deviation from a random coil ensemble. The backbone ¹⁵N transverse relaxation times indicate a mildly reduced mobility for two regions that are rich in bulky residues. The small angle X-ray scattering curve of KISS1 is likewise consistent with a predominantly random coil ensemble, although an ensemble optimization analysis indicates some preference for more extended conformations possibly due to positive charge repulsion between the abundant basic residues. Our results support the hypothesis that KISS1 mostly samples a random coil conformational space, which is consistent with its high susceptibility to proteolysis and the generation of Kisspeptin fragments.

Small-angle scattering studies of intrinsically disordered proteins and their complexes

Intrinsically Disordered Proteins (IDPs) perform a broad range of biological functions. Their relevance has motivated intense research activity seeking to characterize their sequence/structure/function relationships. However, the conformational plasticity of these molecules hampers the application of traditional structural approaches, and new tools and concepts are being developed to address the challenges they pose. Small-Angle Scattering (SAS) is a structural biology technique that probes the size and shape of disordered proteins and their complexes with other biomolecules. The low-resolution nature of SAS can be compensated with specially designed computational tools and its combined interpretation with complementary structural information. In this review, we describe recent advances in the application of SAS to disordered proteins and highly flexible complexes and discuss current challenges.

Quarterly intrinsic disorder digest (January-February-March, 2014)

This is the 5^th issue of the Digested Disorder series that represents a reader's digest of the scientific literature on intrinsically disordered proteins. We continue 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 first quarter of 2014; i.e., during the period of January, February, and March 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.

There is Diversity in Disorder-"In all Chaos there is a Cosmos, in all Disorder a Secret Order"

The protein universe consists of a continuum of structures ranging from full order to complete disorder. As the structured part of the proteome has been intensively studied, stably folded proteins are increasingly well documented and understood. However, proteins that are fully, or in large part, disordered are much less well characterized. Here we collected NMR chemical shifts in a small database for 117 protein sequences that are known to contain disorder. We demonstrate that NMR chemical shift data can be brought to bear as an exquisite judge of protein disorder at the residue level, and help in validation. With the help of secondary chemical shift analysis we demonstrate that the proteins in the database span the full spectrum of disorder, but still, largely segregate into two classes; disordered with small segments of order scattered along the sequence, and structured with small segments of disorder inserted between the different structured regions. A detailed analysis reveals that the distribution of order/disorder along the sequence shows a complex and asymmetric distribution, that is highly protein-dependent. Access to ratified training data further suggests an avenue to improving prediction of disorder from sequence.

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.

Binding cavities and druggability of intrinsically disordered proteins

To assess the potential of intrinsically disordered proteins (IDPs) as drug design targets, we have analyzed the ligand-binding cavities of two datasets of IDPs (containing 37 and 16 entries, respectively) and compared their properties with those of conventional ordered (folded) proteins. IDPs were predicted to possess more binding cavity than ordered proteins at similar length, supporting the proposed advantage of IDPs economizing genome and protein resources. The cavity number has a wide distribution within each conformation ensemble for IDPs. The geometries of the cavities of IDPs differ from the cavities of ordered proteins, for example, the cavities of IDPs have larger surface areas and volumes, and are more likely to be composed of a single segment. The druggability of the cavities was examined, and the average druggable probability is estimated to be 9% for IDPs, which is almost twice that for ordered proteins (5%). Some IDPs with druggable cavities that are associated with diseases are listed. The optimism versus obstacles for drug design for IDPs is also briefly discussed.