Structural Characterization of N-WASP Domain V Using MD Simulations with NMR and SAXS Data

Because of their large conformational heterogeneity, structural characterization of intrinsically disordered proteins (IDPs) is very challenging using classical experimental methods alone. In this study, we use NMR and small-angle x-ray scattering (SAXS) data with multiple molecular dynamics (MD) simulations to describe the conformational ensemble of the fully disordered verprolin homology domain of the neural Aldrich syndrome protein involved in the regulation of actin polymerization. First, we studied several back-calculation software of SAXS scattering intensity and optimized the adjustable parameters to accurately calculate the SAXS intensity from an atomic structure. We also identified the most appropriate force fields for MD simulations of this IDP. Then, we analyzed four conformational ensembles of neural Aldrich syndrome protein verprolin homology domain, two generated with the program flexible-meccano with or without NMR-derived information as input and two others generated by MD simulations with two different force fields. These four conformational ensembles were compared to available NMR and SAXS data for validation. We found that MD simulations with the AMBER-03w force field and the TIP4P/2005s water model are able to correctly describe the conformational ensemble of this 67-residue IDP at both local and global level.     

On the Calculation of SAXS Profiles of Folded and Intrinsically Disordered Proteins from Computer Simulations

Solution techniques such as small-angle X-ray scattering (SAXS) play a central role in structural studies of intrinsically disordered proteins (IDPs); yet, due to low resolution, it is generally necessary to combine SAXS with additional experimental sources of data and to use molecular simulations. Computational methods for the calculation of theoretical SAXS intensity profiles can be separated into two groups, depending on whether the solvent is modeled implicitly as continuous electron density or considered explicitly. The former offers reduced computational cost but requires the definition of a number of free parameters to account for, for example, the excess density of the solvation layer. Overfitting can thus be an issue, particularly when the structural ensemble is unknown. Here, we investigate and show how small variations of the contrast of the hydration shell, δρ, severely affect the outcome, analysis and interpretation of computed SAXS profiles for folded and disordered proteins. For both the folded and disordered proteins studied here, using a default δρ may, in some cases, result in the calculation of non-representative SAXS profiles, leading to an overestimation of their size and a misinterpretation of their structural nature. The solvation layer of the different IDP simulations also impacts their size estimates differently, depending on the protein force field used. The same is not true for the folded protein simulations, suggesting differences in the solvation of the two classes of proteins, and indicating that different force fields optimized for IDPs may cause expansion of the polypeptide chain through different physical mechanisms.     

Effects of Molecular Crowding on the Dynamics of Intrinsically Disordered Proteins

Inside cells, the concentration of macromolecules can reach up to 400 g/L. In such crowded environments, proteins are expected to behave differently than in vitro. It has been shown that the stability and the folding rate of a globular protein can be altered by the excluded volume effect produced by a high density of macromolecules. However, macromolecular crowding effects on intrinsically disordered proteins (IDPs) are less explored. These proteins can be extremely dynamic and potentially sample a wide ensemble of conformations under non-denaturing conditions. The dynamic properties of IDPs are intimately related to the timescale of conformational exchange within the ensemble, which govern target recognition and how these proteins function. In this work, we investigated the macromolecular crowding effects on the dynamics of several IDPs by measuring the NMR spin relaxation parameters of three disordered proteins (ProTα, TC1, and α-synuclein) with different extents of residual structures. To aid the interpretation of experimental results, we also performed an MD simulation of ProTα. Based on the MD analysis, a simple model to correlate the observed changes in relaxation rates to the alteration in protein motions under crowding conditions was proposed. Our results show that 1) IDPs remain at least partially disordered despite the presence of high concentration of other macromolecules, 2) the crowded environment has differential effects on the conformational propensity of distinct regions of an IDP, which may lead to selective stabilization of certain target-binding motifs, and 3) the segmental motions of IDPs on the nanosecond timescale are retained under crowded conditions. These findings strongly suggest that IDPs function as dynamic structural ensembles in cellular environments.     

Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding

The sequence–structure–function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well‐defined three‐dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.     

NMR illuminates intrinsic disorder

Nuclear magnetic resonance (NMR) has long been instrumental in the characterization of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs). This method continues to offer rich insights into the nature of IDPs in solution, especially in combination with other biophysical methods such as small-angle scattering, single-molecule fluorescence, electron paramagnetic resonance (EPR), and mass spectrometry. Substantial advances have been made in recent years in studies of proteins containing both ordered and disordered domains and in the characterization of problematic sequences containing repeated tracts of a single or a few amino acids. These sequences are relevant to disease states such as Alzheimer's, Parkinson's, and Huntington's diseases, where disordered proteins misfold into harmful amyloid. Innovative applications of NMR are providing novel insights into mechanisms of protein aggregation and the complexity of IDP interactions with their targets. As a basis for understanding the solution structural ensembles, dynamic behavior, and functional mechanisms of IDPs and IDRs, NMR continues to prove invaluable.     

Uncoupled binding and folding of immune signaling-related intrinsically disordered proteins

The classical protein structure-function paradigm has been challenged by the emergence of intrinsically disordered proteins (IDPs), the proteins that do not adopt well-defined three-dimensional structures under physiological conditions. This development was accompanied by the introduction of a “coupled binding and folding” paradigm that suggests folding of IDPs upon binding to their partners. However, our recent studies challenge this general view by revealing a novel, previously unrecognized phenomenon – uncoupled binding and folding. This biologically important mechanism is characteristic of members of a new family of IDPs involved in immune signaling and underlies their unusual properties including: (1) specific homodimerization, (2) the lack of folding upon binding to a well-folded protein, another IDP molecule, or to lipid bilayer membranes, and (3) the “scissors-cut paradox”. The third phenomenon occurs in diverse IDP interactions and suggests that properties of IDP fragments are not necessarily additive in the context of the entire protein. The “no disorder-to-order transition” type of binding is distinct from known IDP interactions and is characterized by an unprecedented observation of the lack of chemical shift and peak intensity changes in multidimensional NMR spectra, a fingerprint of proteins, upon complex formation. Here, I focus on those interactions of IDPs with diverse biological partners where the binding phase driven by electrostatic interactions is not be necessarily followed by the hydrophobic folding phase. I also review new multidisciplinary knowledge about immune signaling-related IDPs and show how it expands our understanding of cell function with multiple applications in biology and medicine.     

The Conformational Ensembles of α-Synuclein and Tau: Combining Single-Molecule FRET and Simulations

Intrinsically disordered proteins (IDPs) are increasingly recognized for their important roles in a range of biological contexts, both in normal physiological function and in a variety of devastating human diseases. However, their structural characterization by traditional biophysical methods, for the purposes of understanding their function and dysfunction, has proved challenging. Here, we investigate the model IDPs α-Synuclein (αS) and tau, that are involved in major neurodegenerative conditions including Parkinson’s and Alzheimer’s diseases, using excluded volume Monte Carlo simulations constrained by pairwise distance distributions from single-molecule fluorescence measurements. Using this, to our knowledge, novel approach we find that a relatively small number of intermolecular distance constraints are sufficient to accurately determine the dimensions and polymer conformational statistics of αS and tau in solution. Moreover, this method can detect local changes in αS and tau conformations that correlate with enhanced aggregation. Constrained Monte Carlo simulations produce ensembles that are in excellent agreement both with experimental measurements on αS and tau and with all-atom, explicit solvent molecular dynamics simulations of αS, with much lower configurational sampling requirements and computational expense.Abbreviations used: AAMD, all-atom molecular dynamics; ECMC, experimentally constrained Monte Carlo; ET_eff, energy transfer efficiency; (sm)FRET, (single molecule) Förster resonance energy transfer; IDP, intrinsically disordered protein; LJ, Lennard-Jones; MC, Monte Carlo; MD, molecular dynamics; PRE, paramagnetic relaxation enhancement; SAX(N)S, small-angle x-ray (neutron) scattering; UMC, unconstrained Monte Carlo     

The origin and impact of bound water around intrinsically disordered proteins

Proteins and water couple dynamically over a wide range of time scales. Motivated by their central role in protein function, protein-water dynamics and thermodynamics have been extensively studied for structured proteins, where correspondence to structural features has been made. However, properties controlling intrinsically disordered protein (IDP)-water dynamics are not yet known. We report results of megahertz-to-terahertz dielectric spectroscopy and molecular dynamics simulations of a group of IDPs with varying charge content along with structured proteins of similar size. Hydration water around IDPs is found to exhibit more heterogeneous rotational and translational dynamics compared with water around structured proteins of similar size, yielding on average more restricted dynamics around individual residues of IDPs, charged or neutral, compared with structured proteins. The on-average slower water dynamics is found to arise from excess tightly bound water in the first hydration layer, which is related to greater exposure to charged groups. The more tightly bound water to IDPs correlates with the smaller hydration shell found experimentally, and affects entropy associated with protein-water interactions, the contribution of which we estimate based on the dielectric measurements and simulations. Water-IDP dynamic coupling at terahertz frequencies is characterized by the dielectric measurements and simulations.     

Intrinsically disordered proteins: modes of binding with emphasis on disordered domains

Our notions of protein function have long been determined by the protein structure–function paradigm. However, the idea that protein function is dictated by a prerequisite complementarity of shapes at the binding interface is becoming increasingly challenged. Interactions involving intrinsically disordered proteins (IDPs) have indicated a significant degree of disorder present in the bound state, ranging from static disorder to complete disorder, termed ‘random fuzziness’. This review assesses the anatomy of an IDP and relates how its intrinsic properties permit promiscuity and allow for the various modes of interaction. Furthermore, a mechanistic overview of the types of disordered domains is detailed, while also relating to a recent example and the kinetic and thermodynamic principles governing its formation.     

Alzheimer's-disease-associated conformation of intrinsically disordered tau protein studied by intrinsically disordered protein liquid-phase competitive enzyme-linked immunosorbent assay

Tau protein, the major constituent of paired helical filaments in Alzheimer's disease, belongs to the intrinsically disordered proteins (IDPs). IDPs are an emerging group in the protein kingdom characterized by the absence of a rigid three-dimensional structure. Disordered proteins usually acquire a "functional fold" upon binding to their interaction partner(s). This property of IDPs implies the need for innovative approaches to measure their binding affinity. We have mapped and measured the Alzheimer's-disease-associated epitope on intrinsically disordered tau protein with a novel two-step sandwich competitive enzyme-linked immunosorbent assay (ELISA). This approach allowed us to determine the binding affinity of disordered tau protein in liquid phase without any disturbance to the competitive equilibrium and without any need for covalent or noncovalent modification of tau protein. Furthermore, the global fitting method, used for the reconstruction of tau binding curves, significantly improved the assay readout. The proposed novel competitive ELISA allowed us to determine the changes in the standard Gibbs energy of binding, thus enabling measurement of tau protein conformation in the core of paired helical filaments. IDP competitive ELISA results showed, for the first time, that the tau protein C terminus of the Alzheimer's-disease-derived paired helical filaments core subunit adopts beta-turn type I' fold and is accessible from solution.     

Chain collapse of an amyloidogenic intrinsically disordered protein

Natively unfolded or intrinsically disordered proteins (IDPs) are under intense scrutiny due to their involvement in both normal biological functions and abnormal protein misfolding disorders. Polypeptide chain collapse of amyloidogenic IDPs is believed to play a key role in protein misfolding, oligomerization, and aggregation leading to amyloid fibril formation, which is implicated in a number of human diseases. In this work, we used bovine κ-casein, which serves as an archetypal model protein for amyloidogenic IDPs. Using a variety of biophysical tools involving both prediction and spectroscopic techniques, we first established that monomeric κ-casein adopts a collapsed premolten-globule-like conformational ensemble under physiological conditions. Our time-resolved fluorescence and light-scattering data indicate a change in the mean hydrodynamic radius from ∼4.6 nm to ∼1.9 nm upon chain collapse. We then took the advantage of two cysteines separated by 77 amino-acid residues and covalently labeled them using thiol-reactive pyrene maleimide. This dual-labeled protein demonstrated a strong excimer formation upon renaturation from urea- and acid-denatured states under both equilibrium and kinetic conditions, providing compelling evidence of polypeptide chain collapse under physiological conditions. The implication of the IDP chain collapse in protein aggregation and amyloid formation is also discussed.     

Sequence characteristics responsible for protein-protein interactions in the intrinsically disordered regions of caseins,amelogenins, and small heat-shock proteins

Milk caseins and dental amelogenins are intrinsically disordered proteins (IDPs) that associate with themselves and others. Paradoxically, they are also described as hydrophobic proteins, which is difficult to reconcile with a solvent-exposed conformation. We attempt to resolve this paradox. We show that caseins and amelogenins are not hydrophobic proteins but they are more hydrophobic than most IDPs. Remarkably, uncharged residues from different regions of these mature proteins have a nearly constant average hydropathy but these regions exhibit different charged residue frequencies. A novel sequence analysis method was developed to identify hydrophobic and order-promoting regions that would favor conformational collapse. We found that such regions were uncommon; most hydrophobic and order-promoting residues were adjacent to hydrophilic or disorder-promoting residues. A further reason why caseins and amelogenins do not collapse is their high proportion of disorder-promoting proline residues. We conclude that in these proteins the hydrophobic effect is not large enough to cause conformational collapse but it can contribute, along with polar interactions, to protein-protein interactions. This behaviour is similar to the interaction of the disordered N-terminal region of small heat-shock proteins with either themselves during oligomer formation or other, unfolding, proteins during chaperone action.     

Biophysical Characterisation of Calumenin as a Charged F508del-CFTR Folding Modulator

The cystic fibrosis transmembrane regulator (CFTR) is a cyclic-AMP dependent chloride channel expressed at the apical surface of epithelial cells lining various organs such as the respiratory tract. Defective processing and functioning of this protein caused by mutations in the CFTR gene results in loss of ionic balance, defective mucus clearance, increased proliferation of biofilms and inflammation of human airways observed in cystic fibrosis (CF) patients. The process by which CFTR folds and matures under the influence of various chaperones in the secretory pathway remains incompletely understood. Recently, calumenin, a secretory protein, belonging to the CREC family of low affinity calcium binding proteins has been identified as a putative CFTR chaperone whose biophysical properties and functions remain uncharacterized. We compared hydropathy, instability, charge, unfoldability, disorder and aggregation propensity of calumenin and other CREC family members with CFTR associated chaperones and calcium binding proteins, wild-type and mutant CFTR proteins and intrinsically disordered proteins (IDPs). We observed that calumenin, along with other CREC proteins, was significantly more charged and less folded compared to CFTR associated chaperones. Moreover like IDPs, calumenin and other CREC proteins were found to be less hydrophobic and aggregation prone. Phylogenetic analysis revealed a close link between calumenin and other CREC proteins indicating how evolution might have shaped their similar biophysical properties. Experimentally, calumenin was observed to significantly reduce F508del-CFTR aggregation in a manner similar to AavLEA1, a well-characterized IDP. Fluorescence microscopy based imaging analysis also revealed altered trafficking of calumenin in bronchial cells expressing F508del-CFTR, indicating its direct role in the pathophysiology of CF. In conclusion, calumenin is characterized as a charged protein exhibiting close similarity with IDPs and is hypothesized to regulate F508del-CFTR folding by electrostatic effects. This work provides useful insights for designing optimized synthetic structural correctors of CFTR mutant proteins in the future.     

Sequence Determinants of Compaction in Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs), which lack folded structure and are disordered under nondenaturing conditions, have been shown to perform important functions in a large number of cellular processes. These proteins have interesting structural properties that deviate from the random-coil-like behavior exhibited by chemically denatured proteins. In particular, IDPs are often observed to exhibit significant compaction. In this study, we have analyzed the hydrodynamic radii of a number of IDPs to investigate the sequence determinants of this compaction. Net charge and proline content are observed to be strongly correlated with increased hydrodynamic radii, suggesting that these are the dominant contributors to compaction. Hydrophobicity and secondary structure, on the other hand, appear to have negligible effects on compaction, which implies that the determinants of structure in folded and intrinsically disordered proteins are profoundly different. Finally, we observe that polyhistidine tags seem to increase IDP compaction, which suggests that these tags have significant perturbing effects and thus should be removed before any structural characterizations of IDPs. Using the relationships observed in this analysis, we have developed a sequence-based predictor of hydrodynamic radius for IDPs that shows substantial improvement over a simple model based upon chain length alone.     

Using a FRET Library with Multiple Probe Pairs To Drive Monte Carlo Simulations of α-Synuclein

We describe a strategy for experimentally-constraining computational simulations of intrinsically disordered proteins (IDPs), using α-synuclein, an IDP with a central role in Parkinson’s disease pathology, as an example. Previously, data from single-molecule Förster Resonance Energy Transfer (FRET) experiments have been effectively utilized to generate experimentally constrained computational models of IDPs. However, the fluorophores required for single-molecule FRET experiments are not amenable to the study of short-range (<30 Å) interactions. Using ensemble FRET measurements allows one to acquire data from probes with multiple distance ranges, which can be used to constrain Monte Carlo simulations in PyRosetta. To appropriately employ ensemble FRET data as constraints, we optimized the shape and weight of constraining potentials to afford ensembles of structures that are consistent with experimental data. We also used this approach to examine the structure of α-synuclein in the presence of the compacting osmolyte trimethylamine-N-oxide. Despite significant compaction imparted by 2 M trimethylamine-N-oxide, the underlying ensemble of α-synuclein remains largely disordered and capable of aggregation, also in agreement with experimental data. These proof-of-concept experiments demonstrate that our modeling protocol enables one to efficiently generate experimentally constrained models of IDPs that incorporate atomic-scale detail, allowing one to study an IDP under a variety of conditions.     

Length-dependent compaction of intrinsically disordered proteins

This work investigates the effect of chain length on the degree of compaction of intrinsically disordered proteins (IDPs). The three main IDP types, native coil (NC), pre-molten globule (PMG) and molten globule (MG), are compared by means of a compaction index (CI) normalized for chain length. The results point out a strong variability of compactness as a function of chain length within each group, with larger proteins populating more compact states. While qualitative sequence features are responsible for the main differences among groups, chain length seems to have an unspecific effect modulating the extent of compaction within each group. The results are consistent with a cooperative character of the weak interactions responsible for chain collapse.     

The Disordered Region of the HCV Protein NS5A: Conformational Dynamics,SH3 Binding,and Phosphorylation

Intrinsically disordered proteins (IDPs) perform their physiological role without possessing a well-defined three-dimensional structure. Still, residual structure and conformational dynamics of IDPs are crucial for the mechanisms underlying their functions. For example, regions of transient secondary structure are often involved in molecular recognition, with the structure being stabilized (or not) upon binding. Long-range interactions, on the other hand, determine the hydrodynamic radius of the IDP, and thus the distance over which the protein can catch binding partners via so-called fly-casting mechanisms. The modulation of long-range interactions also presents a convenient way of fine-tuning the protein’s interaction network, by making binding sites more or less accessible. Here we studied, mainly by nuclear magnetic resonance spectroscopy, residual secondary structure and long-range interactions in nonstructural protein 5A (NS5A) from hepatitis C virus (HCV), a typical viral IDP with multiple functions during the viral life cycle. NS5A comprises an N-terminal folded domain, followed by a large (∼250-residue) disordered C-terminal part. Comparing nuclear magnetic resonance spectra of full-length NS5A with those of a protein construct composed of only the C-terminal residues 191–447 (NS5A-D2D3) allowed us to conclude that there is no significant interaction between the globular and disordered parts of NS5A. NS5A-D2D3, despite its overall high flexibility, shows a large extent of local residual (α-helical and β-turn) structure, as well as a network of electrostatic long-range interactions. Furthermore, we could demonstrate that these long-range interactions become modulated upon binding to the host protein Bin1, as well as after NS5A phosphorylation by CK2. As the charged peptide regions involved in these interactions are well conserved among the different HCV genotypes, these transient long-range interactions may be important for some of the functions of NS5A over the course of the HCV life cycle.