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

Scrutinising the Conformational Ensemble of the Intrinsically Mixed-Folded Protein Galectin-3

Galectin-3 is a protein involved in many intra- and extra-cellular processes. It has been identified as a diagnostic or prognostic biomarker for certain types of heart disease, kidney disease and cancer. Galectin-3 comprises a carbohydrate recognition domain (CRD) and an N-terminal domain (NTD), which is unstructured and contains eight collagen-like Pro-Gly-rich tandem repeats. While the structure of the CRD has been solved using protein crystallography, current knowledge about conformations of full-length galectin-3 is limited. To fill in this knowledge gap, we performed molecular dynamics (MD) simulations of full-length galectin-3. We systematically re-scaled the solute-solvent interactions in the Martini 3 force field to obtain the best possible agreement between available data from SAXS experiments and the ensemble of conformations generated in the MD simulations. The simulation conformations were found to be very diverse, as reflected, e.g., by (i) large fluctuations in the radius of gyration, ranging from about 2 to 5 nm, and (ii) multiple transient contacts made by amino acid residues in the NTD. Consistent with evidence from NMR experiments, contacts between the CRD and NTD were observed to not involve the carbohydrate-binding site on the CRD surface. Contacts within the NTD were found to be made most frequently by aromatic residues. Formation of fuzzy complexes with unspecific stoichiometry was observed to be mediated mostly by the NTD. Taken together, we offer a detailed picture of the conformational ensemble of full-length galectin-3, which will be important for explaining the biological functions of this protein at the molecular level.

The structural properties of full-length annexin A11

Annexin A11 (ANXA11) is a calcium-dependent phospholipid-binding protein belonging to the annexin protein family and implicated in the neurodegenerative amyotrophic lateral sclerosis. Structurally, ANXA11 contains a conserved calcium-binding C-terminal domain common to all annexins and a putative intrinsically unfolded N-terminus specific for ANXA11. Little is known about the structure and functions of this region of the protein. By analogy with annexin A1, it was suggested that residues 38 to 59 within the ANXA11 N-terminus could form a helical region that would be involved in interactions. Interestingly, this region contains residues that, when mutated, may lead to clinical manifestations. In the present study, we have studied the structural features of the full-length protein with special attention to the N-terminal region using a combination of biophysical techniques which include nuclear magnetic resonance and small angle X-ray scattering. We show that the N-terminus is intrinsically disordered and that the overall features of the protein are not markedly affected by the presence of calcium. We also analyzed the 38-59 helix hypothesis using synthetic peptides spanning both the wild-type sequence and clinically relevant mutations. We show that the peptides have a remarkable character typical of a native helix and that mutations do not alter the behaviour suggesting that they are required for interactions rather than being structurally important. Our work paves the way to a more thorough understanding of the ANXA11 functions.

Phosphorylation motif dictates GPCR C-terminal domain conformation and arrestin interaction

Arrestin-dependent G protein-coupled receptor (GPCR) signaling pathway is regulated by the phosphorylation state of GPCR's C-terminal domain, but the molecular bases of arrestin:receptor interaction are to be further illuminated. Here we investigated the impact of phosphorylation on the conformational features of the C-terminal region from three rhodopsin-like GPCRs, the vasopressin V2 receptor (V2R), the growth hormone secretagogue or ghrelin receptor type 1a (GHSR), and the β2-adernergic receptor (β2AR). Using phosphomimetic variants, we identified pre-formed secondary structure elements, or short linear motifs (SLiMs), that undergo specific conformational transitions upon phosphorylation. Of importance, such conformational transitions appear to favor arrestin-2 binding. Hence, our results suggest a model in which the phosphorylation-dependent structuration of the GPCR C-terminal regions would modulate arrestin binding and therefore signaling outcomes in arrestin-dependent pathways.

The Action of Chemical Denaturants: From Globular to Intrinsically Disordered Proteins

Proteins perform their many functions by adopting either a minimal number of strictly similar conformations, the native state, or a vast ensemble of highly flexible conformations. In both cases, their structural features are highly influenced by the chemical environment. Even though a plethora of experimental studies have demonstrated the impact of chemical denaturants on protein structure, the molecular mechanism underlying their action is still debated. In the present review, after a brief recapitulation of the main experimental data on protein denaturants, we survey both classical and more recent interpretations of the molecular basis of their action. In particular, we highlight the differences and similarities of the impact that denaturants have on different structural classes of proteins, i.e., globular, intrinsically disordered (IDP), and amyloid-like assemblies. Particular attention has been given to the IDPs, as recent studies are unraveling their fundamental importance in many physiological processes. The role that computation techniques are expected to play in the near future is illustrated.

Structural studies of protein-nucleic acid complexes: A brief overview of the selected techniques

Protein-nucleic acid complexes are involved in all vital processes, including replication, transcription, translation, regulation of gene expression and cell metabolism. Knowledge of the biological functions and molecular mechanisms beyond the activity of the macromolecular complexes can be determined from their tertiary structures. Undoubtably, performing structural studies of protein-nucleic acid complexes is challenging, mainly because these types of complexes are often unstable. In addition, their individual components may display extremely different surface charges, causing the complexes to precipitate at higher concentrations used in many structural studies. Due to the variety of protein-nucleic acid complexes and their different biophysical properties, no simple and universal guideline exists that helps scientists chose a method to successfully determine the structure of a specific protein-nucleic acid complex. In this review, we provide a summary of the following experimental methods, which can be applied to study the structures of protein-nucleic acid complexes: X-ray and neutron crystallography, nuclear magnetic resonance (NMR) spectroscopy, cryogenic electron microscopy (cryo-EM), atomic force microscopy (AFM), small angle scattering (SAS) methods, circular dichroism (CD) and infrared (IR) spectroscopy. Each method is discussed regarding its historical context, advancements over the past decades and recent years, and weaknesses and strengths. When a single method does not provide satisfactory data on the selected protein-nucleic acid complex, a combination of several methods should be considered as a hybrid approach; thus, specific structural problems can be solved when studying protein-nucleic acid complexes.

Interaction of the tau fibrils with the neuronal membrane

Tau proteins are recently gaining a lot of interest due to their active role in causing a range of tauopathies. Molecular mechanisms underlying the tau interaction with the neuronal membrane are hitherto unknown and difficult to characterize using experimental methods. Using the cryo-EM structure of the tau-fibrils we have used atomistic molecular dynamics simulation to model the tau fibril and neuronal membrane interaction using explicit solvation. The dynamics and structural characteristics of the tau fibril with the neuronal membrane are compared to the tau fibril in the aqueous phase to corroborate the effect of the neuronal membrane in the tau structure. Tau fibrils have been modelled using CHARMM-36m force field and the six component neuronal membrane composition is taken from the earlier simulation results. The timescale conceivable in our molecular dynamics simulations is of the order of microseconds which captures the onset of the interaction of the tau fibrils with the neuronal membrane. This interaction is found to impact the tau pathogenesis that finally causes neuronal toxicity. Our study initiates the understanding of tau conformational ensemble in the presence of neuronal membrane and sheds the light on the significant tau-membrane interactions.

Modeling of flexible membrane-bound biomolecular complexes for solution small-angle scattering

Recent advances in protein expression protocols, sample handling, and experimental set up of small-angle scattering experiments have allowed users of the technique to structurally investigate biomolecules of growing complexity and structural disorder. Notable examples include intrinsically disordered proteins, multi-domain proteins and membrane proteins in suitable carrier systems. Here, we outline a modeling scheme for calculating the scattering profiles from such complex samples. This kind of modeling is necessary for structural information to be refined from the corresponding data. The scheme bases itself on a hybrid of classical form factor based modeling and the well-known spherical harmonics-based formulation of small-angle scattering amplitudes. Our framework can account for flexible domains alongside other structurally elaborate components of the molecular system in question. We demonstrate the utility of this modeling scheme through a recent example of a structural model of the growth hormone receptor membrane protein in a phospholipid bilayer nanodisc which is refined against experimental SAXS data. Additionally we investigate how the scattering profiles from the complex would appear under different scattering contrasts. For each contrast situation we discuss what structural information is contained and the related consequences for modeling of the data.

Conformational Search for the Building Block of Proteins Based on the Gradient Gravitational Search Algorithm (ConfGGS) Using Force Fields: CHARMM, AMBER, and OPLS-AA

Proteins are linear polymers built from a repertoire of 20 different amino acids, which are considered building blocks of proteins. The diversity and versatility of these 20 building blocks with regard to their conformations are key to adopting three-dimensional structures that facilitate proteins to undergo important mechanistic biological processes in living systems. The present investigation reports a conformational search of 20 different amino acids, building blocks of proteins, using three different force fields, CHARMM, AMBER, and OPLS-AA, implemented in the gradient gravitational search algorithm. The search technique (ConfGGS) includes the contribution from both bonded and nonbonded terms using Cartesian coordinates. The efficiency of such conformational searches has also been compared with other optimization algorithms: DE/Best, DE/Rand, and PSO algorithms with respect to computational time and accuracy based on the minimum number of iteration steps and computed lowest mean absolute error (MAE) and mean standard deviation (MSD) values for dihedral angles of respective near-optimal structures. Moreover, the ConfGGS technique has also been extended to an ordered protein fragment (PQITL) extracted from HIV-1 protease (PDB ID: 1YTH), an intrinsically disordered protein fragment, i.e., an amyloid-forming segment (AVVTGVTAV), from the NAC domain of Parkinson's disease protein α-synuclein, residues 69-77 (PDB ID: 4RIK), the experimental NMR atomic-resolution structure of α-synuclein fibrils (PDB ID: 2N0A), and a disulfide bond-containing protein fragment sequence (PCYGWPVCY), residues 59-67 (PDB ID: 6Y4F) toward structure prediction as a close homologue compared with experimental accuracy, using the CHARMM force field. The MolProbity validation results for the protein fragment (PQITL) obtained by ConfGGS/CHARMM are in better agreement with the native protein fragment structure of HIV-1 protease (PDB ID: 1YTH). Furthermore, the computed results have also been compared with the coordinates obtained from the AlphaFold network.

How phosphorylation impacts intrinsically disordered proteins and their function

Phosphorylation is the most common post-translational modification (PTM) in eukaryotes, occurring particularly frequently in intrinsically disordered proteins (IDPs). These proteins are highly flexible and dynamic by nature. Thus, it is intriguing that the addition of a single phosphoryl group to a disordered chain can impact its function so dramatically. Furthermore, as many IDPs carry multiple phosphorylation sites, the number of possible states increases, enabling larger complexities and novel mechanisms. Although a chemically simple and well-understood process, the impact of phosphorylation on the conformational ensemble and molecular function of IDPs, not to mention biological output, is highly complex and diverse. Since the discovery of the first phosphorylation site in proteins 75 years ago, we have come to a much better understanding of how this PTM works, but with the diversity of IDPs and their capacity for carrying multiple phosphoryl groups, the complexity grows. In this Essay, we highlight some of the basic effects of IDP phosphorylation, allowing it to serve as starting point when embarking on studies into this topic. We further describe how recent complex cases of multisite phosphorylation of IDPs have been instrumental in widening our view on the effect of protein phosphorylation. Finally, we put forward perspectives on the phosphorylation of IDPs, both in relation to disease and in context of other PTMs; areas where deep insight remains to be uncovered.

PYK2 senses calcium through a disordered dimerization and calmodulin-binding element

Multidomain kinases use many ways to integrate and process diverse stimuli. Here, we investigated the mechanism by which the protein tyrosine kinase 2-beta (PYK2) functions as a sensor and effector of cellular calcium influx. We show that the linker between the PYK2 kinase and FAT domains (KFL) encompasses an unusual calmodulin (CaM) binding element. PYK2 KFL is disordered and engages CaM through an ensemble of transient binding events. Calcium increases the association by promoting structural changes in CaM that expose auxiliary interaction opportunities. KFL also forms fuzzy dimers, and dimerization is enhanced by CaM binding. As a monomer, however, KFL associates with the PYK2 FERM-kinase fragment. Thus, we identify a mechanism whereby calcium influx can promote PYK2 self-association, and hence kinase-activating trans-autophosphorylation. Collectively, our findings describe a flexible protein module that expands the paradigms for CaM binding and self-association, and their use for controlling kinase activity.

Protein tyrosine kinase 2-beta is shown to function as a sensor and effector of cellular calcium influx through self-association.

SAXS-guided unbiased coarse-grained Monte Carlo simulation for identification of self-assembly nanostructures and dimensions

Recent studies have shown that solvated amphiphiles can form nanostructured self-assemblies called dynamic binary complexes (DBCs) in the presence of ions. Since the nanostructures of DBCs are directly related to their viscoelastic properties, it is important to understand how the nanostructures change under different solution conditions. However, it is challenging to obtain a three-dimensional molecular description of these nanostructures by utilizing conventional experimental characterization techniques or thermodynamic models. To this end, we combined the structural data from small angle X-ray scattering (SAXS) experiments and thermodynamic knowledge from coarse-grained Monte Carlo (CGMC) simulations to identify the detailed three-dimensional nanostructure of DBCs. Specifically, unbiased CGMC simulations are performed with SAXS-guided initial conditions, which aids us to sample accurate nanostructures in a computationally efficient fashion. As a result, an elliptical bilayer nanostructure is obtained as the most probable nanostructure of DBCs whose dimensions are validated by scanning electron microscope (SEM) images. Then, utilizing the obtained molecular model of DBCs, we could also explain the pH tunability of the system. Overall, our results from SAXS-guided unbiased CGMC simulations highlight that using potential energy combined with SAXS data, we can distinguish otherwise degenerate nanostructures resulting from the inherent ambiguity of SAXS patterns.

Conformational ensemble of the full-length SARS-CoV-2 nucleocapsid (N) protein based on molecular simulations and SAXS data

The nucleocapsid protein of the SARS-CoV-2 virus comprises two RNA-binding domains and three regions that are intrinsically disordered. While the structures of the RNA-binding domains have been solved using protein crystallography and NMR, current knowledge of the conformations of the full-length nucleocapsid protein is rather limited. To fill in this knowledge gap, we combined coarse-grained molecular simulations with data from small-angle X-ray scattering (SAXS) experiments using the ensemble refinement of SAXS (EROS) method. Our results show that the dimer of the full-length nucleocapsid protein exhibits large conformational fluctuations with its radius of gyration ranging from about 4 to 8 nm. The RNA-binding domains do not make direct contacts. The disordered region that links these two domains comprises a hydrophobic α-helix which makes frequent and nonspecific contacts with the RNA-binding domains. Each of the intrinsically disordered regions adopts conformations that are locally compact, yet on average, much more extended than Gaussian chains of equivalent lengths. We offer a detailed picture of the conformational ensemble of the nucleocapsid protein dimer under near-physiological conditions, which will be important for understanding the nucleocapsid assembly process.

Prion-like strain effects in tauopathies

Tau is a microtubule-associated protein that plays crucial roles in physiology and pathophysiology. In the realm of dementia, tau protein misfolding is associated with a wide spectrum of clinicopathologically diverse neurodegenerative diseases, collectively known as tauopathies. As proposed by the tau strain hypothesis, the intrinsic heterogeneity of tauopathies may be explained by the existence of structurally distinct tau conformers, “strains”. Tau strains can differ in their associated clinical features, neuropathological profiles, and biochemical signatures. Although prior research into infectious prion proteins offers valuable lessons for studying how a protein-only pathogen can encompass strain diversity, the underlying mechanism by which tau subtypes are generated remains poorly understood. Here we summarize recent advances in understanding different tau conformers through in vivo and in vitro experimental paradigms, and the implications of heterogeneity of pathological tau species for drug development.

Structural Insights into the Intrinsically Disordered GPCR C-Terminal Region, Major Actor in Arrestin-GPCR Interaction

Arrestin-dependent pathways are a central component of G protein-coupled receptor (GPCRs) signaling. However, the molecular processes regulating arrestin binding are to be further illuminated, in particular with regard to the structural impact of GPCR C-terminal disordered regions. Here, we used an integrated biophysical strategy to describe the basal conformations of the C-terminal domains of three class A GPCRs, the vasopressin V2 receptor (V2R), the growth hormone secretagogue or ghrelin receptor type 1a (GHSR) and the β2-adernergic receptor (β2AR). By doing so, we revealed the presence of transient secondary structures in these regions that are potentially involved in the interaction with arrestin. These secondary structure elements differ from those described in the literature in interaction with arrestin. This suggests a mechanism where the secondary structure conformational preferences in the C-terminal regions of GPCRs could be a central feature for optimizing arrestins recognition.

Disordered regions tune order in chromatin organization and function

Intrinsically disordered proteins or hybrid proteins with ordered domains and disordered regions (both collectively designated as IDP(R)s) defy the well-established structure-function paradigm due to their ability to perform multiple biological functions even in the absence of a well-defined 3D structure. IDP(R)s have a unique ability to exist as a functional heterogeneous ensemble, where they adopt multiple thermodynamically stable conformations with low energy barriers between states. The resultant structural plasticity or conformational adaptability provides them with a high functional diversity and ease of regulation. Hence, IDP(R)s are highly efficient biological machinery to mediate intricate cellular functions such as signaling, gene expression, and assembly of complex structures. One such structure is the nucleoprotein complex known as Chromatin. Interestingly, the proteins involved in shaping up the structure and function of chromatin are abundant in disordered regions, which serve more than just as mere flexible linkers. The disordered regions are involved in crucial processes such as gene expression regulation, chromatin architecture maintenance, and liquid-liquid phase separation initiation. This review is an attempt to explore the advantages and the functional and regulatory roles of intrinsic disorder in several Chromatin Associated Proteins from a mechanistic standpoint.

Recent advances in small-angle scattering and its expanding impact in structural biology

Applications of small-angle scattering (SAS) in structural biology have benefited from continuing developments in instrumentation, tools for data analysis, modeling capabilities, standards for data and model presentation, and data archiving. The interplay of these capabilities has enabled SAS to contribute to advances in structural biology as the field pushes the boundaries in studies of biomolecular complexes and assemblies as large as whole cells, membrane proteins in lipid environments, and dynamic systems on time scales ranging from femtoseconds to hours. This review covers some of the important advances in biomolecular SAS capabilities for structural biology focused on over the last 5 years and presents highlights of recent applications that demonstrate how the technique is exploring new territories.

Biophysical applications in structural and molecular biology

The main objective of structural biology is to model proteins and other biological macromolecules and link the structural information to function and dynamics. The biological functions of protein molecules and nucleic acids are inherently dependent on their conformational dynamics. Imaging of individual molecules and their dynamic characteristics is an ample source of knowledge that brings new insights about mechanisms of action. The atomic-resolution structural information on most of the biomolecules has been solved by biophysical techniques; either by X-ray diffraction in single crystals or by nuclear magnetic resonance (NMR) spectroscopy in solution. Cryo-electron microscopy (cryo-EM) is emerging as a new tool for analysis of a larger macromolecule that couldn't be solved by X-ray crystallography or NMR. Now a day's low-resolution Cryo-EM is used in combination with either X-ray crystallography or NMR. The present review intends to provide updated information on applications like X-ray crystallography, cryo-EM and NMR which can be used independently and/or together in solving structures of biological macromolecules for our full comprehension of their biological mechanisms.

Identification and characterization of metamorphic proteins: Current and future perspectives

Proteins that can reversibly alternate between distinctly different folds under native conditions are described as being metamorphic. The "metamorphome" is the collection of all metamorphic proteins in the proteome, but it remains unknown the extent to which the proteome is populated by this class of proteins. We propose that uncovering the metamorphome will require a synergy of computational screening of protein sequences to identify potential metamorphic behavior and validation through experimental techniques. This perspective discusses computational and experimental approaches that are currently used to predict and characterize metamorphic proteins as well as the need for developing improved methodologies. Since metamorphic proteins act as molecular switches, understanding their properties and behavior could lead to novel applications of these proteins as sensors in biological or environmental contexts.

Conformational landscape of multidomain SMAD proteins

SMAD transcription factors, the main effectors of the TGFβ (transforming growth factor β) network, have a mixed architecture of globular domains and flexible linkers. Such a complicated architecture precluded the description of their full-length (FL) structure for many years. In this study, we unravel the structures of SMAD4 and SMAD2 proteins through an integrative approach combining Small-angle X-ray scattering, Nuclear Magnetic Resonance spectroscopy, X-ray, and computational modeling. We show that both proteins populate ensembles of conformations, with the globular domains tethered by disordered and flexible linkers, which defines a new dimension of regulation. The flexibility of the linkers facilitates DNA and protein binding and modulates the protein structure. Yet, SMAD4FL is monomeric, whereas SMAD2FL is in different monomer-dimer-trimer states, driven by interactions of the MH2 domains. Dimers are present regardless of the SMAD2FL activation state and concentration. Finally, we propose that SMAD2FL dimers are key building blocks for the quaternary structures of SMAD complexes.

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.

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.

Full structural ensembles of intrinsically disordered proteins from unbiased molecular dynamics simulations

Molecular dynamics (MD) simulation is widely used to complement ensemble-averaged experiments of intrinsically disordered proteins (IDPs). However, MD often suffers from limitations of inaccuracy. Here, we show that enhancing the sampling using Hamiltonian replica-exchange MD (HREMD) led to unbiased and accurate ensembles, reproducing small-angle scattering and NMR chemical shift experiments, for three IDPs of varying sequence properties using two recently optimized force fields, indicating the general applicability of HREMD for IDPs. We further demonstrate that, unlike HREMD, standard MD can reproduce experimental NMR chemical shifts, but not small-angle scattering data, suggesting chemical shifts are insufficient for testing the validity of IDP ensembles. Surprisingly, we reveal that despite differences in their sequence, the inter-chain statistics of all three IDPs are similar for short contour lengths (< 10 residues). The results suggest that the major hurdle of generating an accurate unbiased ensemble for IDPs has now been largely overcome.

PED in 2021: a major update of the protein ensemble database for intrinsically disordered proteins

The Protein Ensemble Database (PED) (https://proteinensemble.org), which holds structural ensembles of intrinsically disordered proteins (IDPs), has been significantly updated and upgraded since its last release in 2016. The new version, PED 4.0, has been completely redesigned and reimplemented with cutting-edge technology and now holds about six times more data (162 versus 24 entries and 242 versus 60 structural ensembles) and a broader representation of state of the art ensemble generation methods than the previous version. The database has a completely renewed graphical interface with an interactive feature viewer for region-based annotations, and provides a series of descriptors of the qualitative and quantitative properties of the ensembles. High quality of the data is guaranteed by a new submission process, which combines both automatic and manual evaluation steps. A team of biocurators integrate structured metadata describing the ensemble generation methodology, experimental constraints and conditions. A new search engine allows the user to build advanced queries and search all entry fields including cross-references to IDP-related resources such as DisProt, MobiDB, BMRB and SASBDB. We expect that the renewed PED will be useful for researchers interested in the atomic-level understanding of IDP function, and promote the rational, structure-based design of IDP-targeting drugs.

Molecular simulations of IDPs: From ensemble generation to IDP interactions leading to disorder-to-order transitions

Intrinsically disordered proteins (IDPs) lack a well-defined three-dimensional structure but do exhibit some dynamical and structural ordering. The structural plasticity of IDPs indicates that entropy-driven motions are crucial for their function. Many IDPs undergo function-related disorder-to-order transitions upon by their interaction with specific binding partners. Approaches that are based on both experimental and theoretical tools enable the biophysical characterization of IDPs. Molecular simulations provide insights into IDP structural ensembles and disorder-to-order transition mechanisms. However, such studies depend strongly on the chosen force field parameters and simulation techniques. In this chapter, we provide an overview of IDP characteristics, review all-atom force fields recently developed for IDPs, and present molecular dynamics-based simulation methods that allow IDP ensemble generation as well as the characterization of disorder-to-order transitions. In particular, we introduce metadynamics, replica exchange molecular dynamics simulations, and also kinetic models resulting from Markov State modeling, and provide various examples for the successful application of these simulation methods to IDPs.

Solution NMR insights into dynamic supramolecular assemblies of disordered amyloidogenic proteins

The extraordinary flexibility and structural heterogeneity of intrinsically disordered proteins (IDP) make them functionally versatile molecules. We have now begun to better understand their fundamental role in biology, however many aspects of their behaviour remain difficult to grasp experimentally. This is especially true for the intermolecular interactions which lead to the formation of transient or highly dynamic supramolecular self-assemblies, such as oligomers, aggregation intermediates and biomolecular condensates. Both the emerging functions and pathogenicity of these structures have stimulated great efforts to develop methodologies capable of providing useful insights. Significant progress in solution NMR spectroscopy has made this technique one of the most powerful to describe structural and dynamic features of IDPs within such assemblies at atomic resolution. Here, we review the most recent works that have illuminated key aspects of IDP assemblies and contributed significant advancements towards our understanding of the complex conformational landscape of prototypical disease-associated proteins. We also include a primer on some of the fundamental and innovative NMR methods being used in the discussed studies.

Dynamics of the intrinsically disordered inhibitor IF7 of glutamine synthetase in isolation and in complex with its partner

Glutamine synthetase (GS) catalyzes the ATP-dependent formation of glutamine from glutamate and ammonia. The activity of Synechocystis sp. PCC 6803 GS is regulated, among other mechanisms, by protein-protein interactions with a 65-residue-long, intrinsically disordered protein (IDP), named IF7. IDPs explore diverse conformations in their free states and, in some cases, in their molecular complexes. We used both nuclear magnetic resonance (NMR) at 11.7 T and small angle X-ray scattering (SAXS) to study the size and the dynamics in the picoseconds-to-nanosecond (ps-ns) timescale of: (i) isolated IF7; and (ii) the IF7/GS complex. Our SAXS findings, together with MD results, show: (i) some of the possible IF7 structures in solution; and, (ii) that the presence of IF7 affected the structure of GS in solution. The joint use of SAXS and NMR shows that movements of each amino acid of IF7 were uncorrelated with those of its neighbors. Residues of IF7 with the largest values of the relaxation rates (R₁, R₂ and η_xy), in the free and bound species, were mainly clustered around: (i) the C terminus of the protein; and (ii) Ala30. These residues, together with Arg8 (which is a hot-spot residue in the interaction with GS), had a restricted mobility in the presence of GS. The C-terminal region, which appeared more compact in our MD simulations of isolated IF7, seemed to be involved in non-native contacts with GS that help in the binding between the two macromolecules.

New technologies to analyse protein function: an intrinsic disorder perspective

Functions of intrinsically disordered proteins do not require structure. Such structure-independent functionality has melted away the classic rigid "lock and key" representation of structure-function relationships in proteins, opening a new page in protein science, where molten keys operate on melted locks and where conformational flexibility and intrinsic disorder, structural plasticity and extreme malleability, multifunctionality and binding promiscuity represent a new-fangled reality. Analysis and understanding of this new reality require novel tools, and some of the techniques elaborated for the examination of intrinsically disordered protein functions are outlined in this review.

Structural Modeling Using Solution Small-Angle X-ray Scattering (SAXS)

Small-angle X-ray scattering (SAXS) offers a way to examine the overall shape and oligomerization state of biological macromolecules under quasi native conditions in solution. In the past decades, SAXS has become a standard tool for structure biologists due to the availability of high brilliance X-ray sources and the development of data analysis/interpretation methods. Sample handling robots and software pipelines have significantly reduced the time necessary to conduct SAXS experiments. Presently, most synchrotrons feature beamlines dedicated to biological SAXS, and the SAXS-derived models are deposited into dedicated and accessible databases. The size of macromolecules that may be analyzed ranges from small peptides or snippets of nucleic acids to gigadalton large complexes or even entire viruses. Compared to other structural methods, sample preparation is straightforward, and the risk of inducing preparation artefacts is minimal. Very importantly, SAXS is a method of choice to study flexible systems like unfolded or disordered proteins, providing the structural ensembles compatible with the data. Although it may be utilized stand-alone, SAXS profits a lot from available experimental or predicted high-resolution data and information from complementary biophysical methods. Here, we show the basic principles of SAXS and review latest developments in the fields of hybrid modeling and flexible systems.

Structural Analyses of Intrinsically Disordered Proteins by Small-Angle X-Ray Scattering

Small-angle X-ray scattering (SAXS) is a low-resolution method for the structural characterization of biological macromolecules in solution. Information about the overall structural features provided by SAXS is highly complementary to X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy, and cryo-electron microscopy, which are high-resolution methods. SAXS not only provides the shape, oligomeric state, and quaternary structure of folded proteins and protein complexes but also allows for quantitative analysis of flexible biomolecules. In this chapter, the most relevant SAXS procedures for structural characterization of flexible macromolecules, including intrinsically disordered proteins (IDPs), are presented. The sample requirements for SAXS experiments on protein solutions and the sequence of steps in data collection and processing are described. The use of the advanced data analysis tools to quantitatively characterize flexible proteins is presented in detail. Typical experimental issues and potential problems encountered during SAXS data measurements and analyses are discussed.

Structural Characterization of Protein-Protein Interactions with pyDockSAXS

Structural characterization of protein-protein interactions can provide essential details to understand biological functions at the molecular level and to facilitate their manipulation for biotechnological and biomedical purposes. Unfortunately, the 3D structure is available for only a small fraction of all possible protein-protein interactions, due to the technical limitations of high-resolution structural determination methods. In this context, low-resolution structural techniques, such as small-angle X-ray scattering (SAXS), can be combined with computational docking to provide structural models of protein-protein interactions at large scale. In this chapter, we describe the pyDockSAXS web server ( https://life.bsc.es/pid/pydocksaxs ), which uses pyDock docking and scoring to provide structural models that optimally satisfy the input SAXS data. This server, which is freely available to the scientific community, provides an automatic pipeline to model the structure of a protein-protein complex from SAXS data.

Extensive tests and evaluation of the CHARMM36IDPSFF force field for intrinsically disordered proteins and folded proteins

Intrinsically disordered proteins (IDPs) have received increasing attention in recent studies due to their structural heterogeneity and critical biological functions. To fully understand the structural properties and determine accurate ensembles of IDPs, molecular dynamics (MD) simulation was widely used to sample diverse conformations and reveal the structural dynamics. However, the classical state-of-the-art force fields perform well for folded proteins while being unsatisfactory for the simulations of disordered proteins reported in many previous studies. Thus, improved force fields were developed to precisely describe both folded proteins and disordered proteins. Preliminary tests show that our newly developed CHARMM36IDPSFF (C36IDPSFF) force field can well reproduce the experimental observables of several disordered proteins, but more tests on different types of proteins are needed to further evaluate the performance of C36IDPSFF. Here, we extensively simulate short peptides, disordered proteins, and fast-folding proteins as well as folded proteins, and compare the simulated results with the experimental observables. The simulation results show that C36IDPSFF could substantially reproduce the experimental observables for most of the tested proteins but some limitations are also found in the radius of gyration of large disordered proteins and the stability of fast-folding proteins. This force field will facilitate large scale studies of protein structural dynamics and functions using MD simulations.

Gas-phase conformations of intrinsically disordered proteins and their complexes with ligands: Kinetically trapped states during transfer from solution to the gas phase

Flexible structures of intrinsically disordered proteins (IDPs) are crucial for versatile functions in living organisms, which involve interaction with diverse partners. Electrospray ionization ion mobility mass spectrometry (ESI-IM-MS) has been widely applied for structural characterization of apo-state and ligand-associated IDPs via two-dimensional separation in the gas phase. Gas-phase IDP structures have been regarded as kinetically trapped states originated from conformational features in solution. However, an implication of the states remains elusive in the structural characterization of IDPs, because it is unclear what structural property of IDPs is preserved. Recent studies have indicated that the conformational features of IDPs in solution are not fully reproduced in the gas phase. Nevertheless, the molecular interactions captured in the gas phase amplify the structural differences between IDP conformers. Therefore, an IDP conformational change that is not observed in solution is observable in the gas-phase structures obtained by ESI-IM-MS. Herein, we have presented up-to-date researches on the key implications of kinetically trapped states in the gas phase with a brief summary of the structural dynamics of IDPs in ESI-IM-MS.

Conformational selection of the intrinsically disordered plant stress protein COR15A in response to solution osmolarity - an X-ray and light scattering study

The plant stress protein COR15A stabilizes chloroplast membranes during freezing. COR15A is an intrinsically disordered protein (IDP) in aqueous solution, but acquires an α-helical structure during dehydration or the increase of solution osmolarity. We have used small- and wide-angle X-ray scattering (SAXS/WAXS) combined with static and dynamic light scattering (SLS/DLS) to investigate the structural and hydrodynamic properties of COR15A in response to increasing solution osmolarity. Coarse-grained ensemble modelling allowed a structure-based interpretation of the SAXS data. Our results demonstrate that COR15A behaves as a biomacromolecule with polymer-like properties which strongly depend on solution osmolarity. Biomacromolecular self-assembly occurring at high solvent osmolarity is initiated by the occurrence of two specific structural subpopulations of the COR15A monomer. The osmolarity dependent structural selection mechanism is an elegant way for conformational regulation and assembly of COR15A. It highlights the importance of the polymer-like properties of IDPs for their associated biological function.

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.

Contrasting Assemblies of Oppositely Charged Proteins

Oppositely charged proteins can form soluble assemblies that under specific physical chemical conditions lead to liquid-liquid phase separation, also called heteroprotein coacervation. Increasing evidence suggests that surface charge anisotropy plays a key role in heteroprotein complexation, and coacervation. Here, we investigated complexation of an acidic protein, β-lactoglobulin (BLG), with two basic proteins, rapeseed napin (NAP) and lysozyme (LYS), of similar net charge and size but differing in surface charge distribution. Using turbidity measurements and isothermal titration calorimetry, we confirmed that LYS binds BLG as expected from previous studies. This interaction leads to two types of phase separation phenomena, depending on pH: liquid-solid phase separation in the case of strong electrostatic attraction and liquid-liquid phase separation for weaker attraction. More interestingly, we showed using dynamic light scattering that NAP interacts with BLG, resulting in formation of assemblies in the nanometer size range. The formation of assemblies was also evident when modeling the interactions using Brownian dynamics for both BLG + NAP and BLG + LYS. Similarly, to DLS, BLG and NAP formed smaller assemblies than BLG with LYS. The molecular details rather than the net charge of BLG and NAP may therefore play a role in their assembly. Furthermore, simulated BLG + NAP assemblies were larger than those experimentally detected by DLS. We discuss the discrepancy between experiments and simulations in relation to the limitations of modelling precisely the molecular characteristics of proteins.

An arsenal of methods for the experimental characterization of intrinsically disordered proteins - How to choose and combine them?

In this review, we detail the most common experimental approaches to assess and characterize protein intrinsic structural disorder, with the notable exception of NMR and EPR spectroscopy, two ideally suited approaches that will be described in depth in two other reviews within this special issue. We discuss the advantages, the limitations, as well as the caveats of the various methods. We also describe less common and more demanding approaches that enable achieving further insights into the conformational properties of IDPs. Finally, we present recent developments that have enabled assessment of structural disorder in living cells, and discuss the currently available methods to model IDPs as conformational ensembles.

Phosphatidylinositol 4-kinase IIIβ (PI4KB) forms highly flexible heterocomplexes that include ACBD3, 14-3-3, and Rab11 proteins

Phosphatidylinositol 4-kinase IIIβ (PI4KB) is a key enzyme of the Golgi system because it produces its lipid hallmark - the phosphatidylinositol 4-phosphate (PI4P). It is recruited to Golgi by the Golgi resident ACBD3 protein, regulated by 14-3-3 proteins and it also serves as an adaptor because it recruits the small GTPase Rab11. Here, we analyzed the protein complexes formed by PI4KB in vitro using small angle x-ray scattering (SAXS) and we discovered that these protein complexes are highly flexible. The 14-3-3:PI4KB:Rab11 protein complex has 2:1:1 stoichiometry and its different conformations are rather compact, however, the ACBD3:PI4KB protein complex has both, very compact and very extended conformations. Furthermore, in vitro reconstitution revealed that the membrane is necessary for the formation of ACBD3:PI4KB:Rab11 protein complex at physiological (nanomolar) concentrations.

SAXS-guided Enhanced Unbiased Sampling for Structure Determination of Proteins and Complexes

Molecular simulations can be utilized to predict protein structure ensembles and dynamics, though sufficient sampling of molecular ensembles and identification of key biologically relevant conformations remains challenging. Low-resolution experimental techniques provide valuable structural information on biomolecule at near-native conditions, which are often combined with molecular simulations to determine and refine protein structural ensembles. In this study, we demonstrate how small angle x-ray scattering (SAXS) information can be incorporated in Markov state model-based adaptive sampling strategy to enhance time efficiency of unbiased MD simulations and identify functionally relevant conformations of proteins and complexes. Our results show that using SAXS data combined with additional information, such as thermodynamics and distance restraints, we are able to distinguish otherwise degenerate structures due to the inherent ambiguity of SAXS pattern. We further demonstrate that adaptive sampling guided by SAXS and hybrid information can significantly reduce the computation time required to discover target structures. Overall, our findings demonstrate the potential of this hybrid approach in predicting near-native structures of proteins and complexes. Other low-resolution experimental information can be incorporated in a similar manner to collectively enhance unbiased sampling and improve the accuracy of structure prediction from simulation.

Characterization of intrinsically disordered proteins and their dynamic complexes: From in vitro to cell-like environments

Over the last two decades, it has become increasingly clear that a large fraction of the human proteome is intrinsically disordered or contains disordered segments of significant length. These intrinsically disordered proteins (IDPs) play important regulatory roles throughout biology, underlining the importance of understanding their conformational behavior and interaction mechanisms at the molecular level. Here we review recent progress in the NMR characterization of the structure and dynamics of IDPs in various functional states and environments. We describe the complementarity of different NMR parameters for quantifying the conformational propensities of IDPs in their isolated and phosphorylated states, and we discuss the challenges associated with obtaining structural models of dynamic protein-protein complexes involving IDPs. In addition, we review recent progress in understanding the conformational behavior of IDPs in cell-like environments such as in the presence of crowding agents, in membrane-less organelles and in the complex environment of the human cell.

Intrinsic disorder in the regulatory N-terminal domain of diacylglycerol acyltransferase 1 from Brassica napus

Proteins with multifunctional regulatory domains often demonstrate structural plasticity or protein disorder, allowing the binding of multiple regulatory factors and post-translational modifications. While the importance of protein disorder is clear, it also poses a challenge for in vitro characterization. Here, we report protein intrinsic disorder in a plant molecular system, which despite its prevalence is less studied. We present a detailed biophysical characterization of the entire cytoplasmic N-terminal domain of Brassica napus diacylglycerol acyltransferase, (DGAT1), which includes an inhibitory module and allosteric binding sites. Our results demonstrate that the monomeric N-terminal domain can be stabilized for biophysical characterization and is largely intrinsically disordered in solution. This domain interacts with allosteric modulators of DGAT1, CoA and oleoyl-CoA, at micromolar concentrations. While solution scattering studies indicate conformational heterogeneity in the N-terminal domain of DGAT1, there is a small gain of secondary structure induced by ligand binding.

Role of Computational Methods in Going beyond X-ray Crystallography to Explore Protein Structure and Dynamics

Protein structural biology came a long way since the determination of the first three-dimensional structure of myoglobin about six decades ago. Across this period, X-ray crystallography was the most important experimental method for gaining atomic-resolution insight into protein structures. However, as the role of dynamics gained importance in the function of proteins, the limitations of X-ray crystallography in not being able to capture dynamics came to the forefront. Computational methods proved to be immensely successful in understanding protein dynamics in solution, and they continue to improve in terms of both the scale and the types of systems that can be studied. In this review, we briefly discuss the limitations of X-ray crystallography in studying protein dynamics, and then provide an overview of different computational methods that are instrumental in understanding the dynamics of proteins and biomacromolecular complexes.

Elucidating binding mechanisms and dynamics of intrinsically disordered protein complexes using NMR spectroscopy

Advances in characterizing complexes of intrinsically disordered proteins (IDPs) have led to the discovery of a remarkably diverse interaction landscape that includes folding-upon-binding, highly dynamic complexes, multivalent interactions as well as regulatory switches controlled by post-translational modifications. Nuclear magnetic resonance (NMR) spectroscopy has in recent years made significant contributions to this field by describing the binding mechanisms and mapping conformational dynamics on multiple time scales. Importantly, this progress has been associated with specific methodological developments in NMR, for example in exchange techniques, allowing challenging biological systems to be studied at atomic resolution. In general, the level of dynamics observed in IDP complexes does not correlate with binding affinities, demonstrating the intricate relationship between conformational dynamics and IDP regulatory function.