A structural model for unfolded proteins from residual dipolar couplings and small-angle x-ray scattering

Structural disorder within Henipavirus nucleoprotein and phosphoprotein: from predictions to experimental assessment

Henipaviruses are newly emerged viruses within the Paramyxoviridae family. Their negative-strand RNA genome is packaged by the nucleoprotein (N) within alpha-helical nucleocapsid that recruits the polymerase complex made of the L protein and the phosphoprotein (P). To date structural data on Henipaviruses are scarce, and their N and P proteins have never been characterized so far. Using both computational and experimental approaches we herein show that Henipaviruses N and P proteins possess large intrinsically disordered regions. By combining several disorder prediction methods, we show that the N-terminal domain of P (PNT) and the C-terminal domain of N (NTAIL) are both mostly disordered, although they contain short order-prone segments. We then report the cloning, the bacterial expression, purification and characterization of Henipavirus PNT and NTAIL domains. By combining gel filtration, dynamic light scattering, circular dichroism and nuclear magnetic resonance, we show that both NTAIL and PNT belong to the premolten globule sub-family within the class of intrinsically disordered proteins. This study is the first reported experimental characterization of Henipavirus P and N proteins. The evidence that their respective N-terminal and C-terminal domains are highly disordered under native conditions is expected to be invaluable for future structural studies by helping to delineate N and P protein domains amenable to crystallization. In addition, following previous hints establishing a relationship between structural disorder and protein interactivity, the present results suggest that Henipavirus PNT and NTAIL domains could be involved in manifold protein-protein interactions.

Influence of side chain conformations on local conformational features of amino acids and implication for force field development

Statistical analysis of coil regions in protein structures has been used to obtain the local backbone phi, psi preferences of amino acids, which agree well with the NMR experiments of unfolded peptides and proteins. We analyzed the conformational features of amino acid residues in a restricted coil library of 4220 high-resolution protein crystal structures. In addition to Gly, Ala, and Pro, the phi, psi distribution (Ramachandran plot) of each amino acid is analyzed with respect to three side chain conformers: g+ (chi(1) approximately -60 degrees), g- (chi(1) approximately 60 degrees), and t (chi(1) approximately 180 degrees). The statistical study indicates that the effect of side chain conformations on phi, psi distributions is even greater than the effect of amino acid types. On the basis of the chi(1), phi, psi conformational preferences, the amino acids in addition to Gly, Pro, and Ala can be divided into five types: (1) ordinary amino acids, (2) Ser, (3) Asp and Asn, (4) Val and Ile, and (5) Thr, each with distinguished chi(1) rotamers. The alpha-helix, beta-sheet, and type-I beta-turn preferences of the different rotamers of various amino acid types can be captured by their intrinsic phi, psi preferences from our coil library. Molecular dynamics simulations of dipeptide Ac-X-NHMe and tetrapeptide Ac-A-X-A-NHMe models give nearly the same side chain rotamer distributions. However, for many amino acids, both OPLS-AA/L and AMBER-FF03 force fields give very different chi(1) rotamer distributions from the coil library. This may partially explain why dipeptide models sometimes cannot reproduce those of protein structures well. The current coil library analysis may be valuable in improving the force field for protein simulations.

The cold denatured state of the C-terminal domain of protein L9 is compact and contains both native and non-native structure

Cold denaturation is a general property of globular proteins, and the process provides insight into the origins of the cooperativity of protein folding and the nature of partially folded states. Unfortunately, studies of protein cold denaturation have been hindered by the fact that the cold denatured state is normally difficult to access experimentally. Special conditions such as addition of high concentrations of denaturant, encapsulation into reverse micelles, the formation of emulsified solutions, high pressure, or extremes of pH have been applied, but these can perturb the unfolded state of proteins. The cold denatured state of the C-terminal domain of the ribosomal protein L9 can be populated under native-like conditions by taking advantage of a destabilizing point mutation which leads to cold denaturation at temperatures above 0 degrees C. This state is in slow exchange with the native state on the NMR time scale. Virtually complete backbone (15)N, (13)C, and (1)H as well as side-chain (13)C(beta) and (1)H(beta) chemical shift assignments were obtained for the cold denatured state at pH 5.7, 12 degrees C. Chemical shift analysis, backbone N-H residual dipolar couplings, amide proton NOEs, and R(2) relaxation rates all indicate that the cold denatured state of CTL9 (the C-terminal domain of the ribosomal protein L9) not only contains significant native-like secondary structure but also non-native structure. The regions corresponding to the two native alpha-helices show a strong tendency to populate helical Phi and Psi angles. The segment which connects alpha-helix 2 and beta-strand 2 (residues 107-124) in the native state exhibits a significant preference to form non-native helical structure in the cold denatured state. The structure observed in the cold denatured state of the I98A mutant is similar to that observed in the pH 3.8 unfolded state of wild type CTL9 at 25 degrees C, suggesting that it is a robust feature of the denatured state ensemble of this protein. The implications for protein folding and for studies of cold denatured states are discussed.

Small-angle X-ray and neutron scattering as a tool for structural systems biology

Small-angle scattering (SAS) of X-rays and neutrons reveals low-resolution structures of biological macromolecules in solution. With the recent experimental and methodological advances, SAS became a unique tool for characterising biological systems. The method covers an extremely broad range of molecule sizes (from a few kDa to hundreds of MDa) and experimental conditions (temperature, pH, salinity, ligand addition, etc.), which is of primary importance for a systemic approach in structural biology. The method provides unique information about the overall structure and conformational changes of native individual proteins, functional complexes, flexible macromolecules and hierarchical systems. New developments in small-angle X-ray and neutron scattering studies of biological macromolecules in solution are briefly reviewed, with a special emphasis on technical and methodological approaches useful for structural systems biology. Possibilities of synergistic use of the method with other techniques are considered.

Structural characterization of proteins and complexes using small-angle X-ray solution scattering

Small-angle scattering of X-rays (SAXS) is an established method for the low-resolution structural characterization of biological macromolecules in solution. The technique provides three-dimensional low-resolution structures, using ab initio and rigid body modeling, and allow one to assess the oligomeric state of proteins and protein complexes. In addition, SAXS is a powerful tool for structure validation and the quantitative analysis of flexible systems, and is highly complementary to the high resolution methods of X-ray crystallography and NMR. At present, SAXS analysis methods have reached an advanced state, allowing for automated and rapid characterization of protein solutions in terms of low-resolution models, quaternary structure and oligomeric composition. In this communication, main approaches to the characterization of proteins and protein complexes using SAXS are reviewed. The tools for the analysis of proteins in solution are presented, and the impact that these tools have made in modern structural biology is discussed.

Side-chain chi(1) conformations in urea-denatured ubiquitin and protein G from (3)J coupling constants and residual dipolar couplings

Current NMR information on side-chain conformations of unfolded protein states is sparse due to the poor dispersion particularly of side-chain proton resonances. We present here optimized schemes for the detection of (3)J(HalphaHbeta), (3)J(NHbeta), and (3)J(C'Hbeta) scalar and (1)D(CbetaHbeta) residual dipolar couplings (RDCs) in unfolded proteins. For urea-denatured ubiquitin and protein G, up to six (3)J-couplings to (1)H(beta) are detected, which define the chi(1) angle at very high precision. Interpretation of the (3)J couplings by a model of mixed staggered chi(1) rotamers yields excellent agreement and also provides stereoassignments for (1)H(beta) methylene protons. For all observed amino acids with the exception of leucine, the chemical shift of (1)H(beta3) protons was found downfield from (1)H(beta2). For most residues, the precision of individual chi(1) rotamer populations is better than 2%. The experimental chi(1) rotamer populations are in the vicinity of averages obtained from coil regions in folded protein structures. However, individual variations from these averages of up to 40% are highly significant and indicate sequence- and residue-specific interactions. Particularly strong deviations from the coil average are found for serine and threonine residues, an effect that may be explained by a weakening of side-chain to backbone hydrogen bonds in the urea-denatured state. The measured (1)D(CbetaHbeta) RDCs correlate well with predicted RDCs that were calculated from a sterically aligned coil model ensemble and the (3)J-derived chi(1) rotamer populations. This agreement supports the coil model as a good first approximation of the unfolded state. Deviations between measured and predicted values at certain sequence locations indicate that the description of the local backbone conformations can be improved by incorporation of the RDC information. The ease of detection of a large number of highly precise side-chain RDCs opens the possibility for a more rigorous characterization of both side-chain and backbone conformations in unfolded proteins.

NMR characterization of hydrophobic collapses in amyloidogenic unfolded states and their implications for amyloid formation

NMR spectroscopy was used to characterize hydrophobic clusters in amyloidogenic unfolded states of a protein and their implications for amyloid formation. Three local hydrophobic clusters were observed in the amyloidogenic state of the phosphatidylinositol 3-kinase (PI3K) SH3 domain. Our NMR studies showed that residues with high average area buried upon folding (AABUF) parameter collapsed to form the clusters. Interestingly, the hydrophobic collapses were not stabilized by long-range tertiary interactions among the clusters that were typically observed in non-amyloidogenic unfolded states of various proteins. The lack of the long-range interactions may be a critical property of the amyloidogenic unfolded state. The SH3 domain was also engineered to disrupt one of the clusters by a single-point mutagenesis (W55G), which allowed us to investigate the effect of the clustering on folding and misfolding. The mutant form of the SH3 domain was not able to fold under folding conditions of the wild type protein (pH 3.6-4.0), supporting the cooperative folding hypothesis. However, aggregation properties of the mutant form were not influenced by the mutation, suggesting the SH3 domain forms amyloid via a non-cooperative process.

Design and structure of an equilibrium protein folding intermediate: a hint into dynamical regions of proteins

Partly unfolded protein conformations close to the native state may play important roles in protein function and in protein misfolding. Structural analyses of such conformations which are essential for their fully physicochemical understanding are complicated by their characteristic low populations at equilibrium. We stabilize here with a single mutation the equilibrium intermediate of apoflavodoxin thermal unfolding and determine its solution structure by NMR. It consists of a large native region identical with that observed in the X-ray structure of the wild-type protein plus an unfolded region. Small-angle X-ray scattering analysis indicates that the calculated ensemble of structures is consistent with the actual degree of expansion of the intermediate. The unfolded region encompasses discontinuous sequence segments that cluster in the 3D structure of the native protein forming the FMN cofactor binding loops and the binding site of a variety of partner proteins. Analysis of the apoflavodoxin inner interfaces reveals that those becoming destabilized in the intermediate are more polar than other inner interfaces of the protein. Natively folded proteins contain hydrophobic cores formed by the packing of hydrophobic surfaces, while natively unfolded proteins are rich in polar residues. The structure of the apoflavodoxin thermal intermediate suggests that the regions of natively folded proteins that are easily responsive to thermal activation may contain cores of intermediate hydrophobicity.

Why Hofmeister effects of many salts favor protein folding but not DNA helix formation

The majority ( approximately 70%) of surface buried in protein folding is hydrocarbon, whereas in DNA helix formation, the majority ( approximately 65%) of surface buried is relatively polar nitrogen and oxygen. Our previous quantification of salt exclusion from hydrocarbon (C) accessible surface area (ASA) and accumulation at amide nitrogen (N) and oxygen (O) ASA leads to a prediction of very different Hofmeister effects on processes that bury mostly polar (N, O) surface compared to the range of effects commonly observed for processes that bury mainly nonpolar (C) surface, e.g., micelle formation and protein folding. Here we quantify the effects of salts on folding of the monomeric DNA binding domain (DBD) of lac repressor (lac DBD) and on formation of an oligomeric DNA duplex. In accord with this prediction, no salt investigated has a stabilizing Hofmeister effect on DNA helix formation. Our ASA-based analyses of model compound data and estimates of the surface area buried in protein folding and DNA helix formation allow us to predict Hofmeister effects on these processes. We observe semiquantitative to quantitative agreement between these predictions and the experimental values, obtained from a novel separation of coulombic and Hofmeister effects. Possible explanations of deviations, including salt-dependent unfolded ensembles and interactions with other types of surface, are discussed.

Ensemble calculations of unstructured proteins constrained by RDC and PRE data: a case study of urea-denatured ubiquitin

The detailed, quantitative characterization of unfolded proteins is a largely unresolved task due to the enormous experimental and theoretical difficulties in describing the highly dimensional space of their conformational ensembles. Recently, residual dipolar coupling (RDC) and paramagnetic relaxation enhancement (PRE) data have provided large numbers of experimental parameters on unfolded states. To obtain a minimal model of the unfolded state according to such data we have developed new modules for the use of steric alignment RDCs and PREs as constraints in ensemble structure calculations by the program XPLOR-NIH. As an example, ensemble calculations were carried out on urea-denatured ubiquitin using a total of 419 previously obtained RDCs and 253 newly determined PREs from eight cysteine mutants coupled to MTSL. The results show that only a small number of about 10 conformers is necessary to fully reproduce the experimental RDCs, PREs and average radius of gyration. C(alpha) contacts determined on a large set (400) of 10-conformer ensembles show significant (10-20%) populations of conformations that are similar to ubiquitin's A-state, i.e. corresponding to an intact native first beta-hairpin and alpha-helix as well as non-native alpha-helical conformations in the C-terminal half. Thus, methanol/acid (A-state) and urea denaturation lead to similar low energy states of the protein ensemble, presumably due to the weakening of the hydrophobic core. Similar contacts are obtained in calculations using solely RDCs or PREs. The sampling statistics of the C(alpha) contacts in the ensembles follow a simple binomial distribution. It follows that the present RDC, PRE, and computational methods allow the statistically significant detection of subconformations in the unfolded ensemble at population levels of a few percent.

Quantitative description of backbone conformational sampling of unfolded proteins at amino acid resolution from NMR residual dipolar couplings

An atomic resolution characterization of the structural properties of unfolded proteins that explicitly invokes the highly dynamic nature of the unfolded state will be extremely important for the development of a quantitative understanding of the thermodynamic basis of protein folding and stability. Here we develop a novel approach using residual dipolar couplings (RDCs) from unfolded proteins to determine conformational behavior on an amino acid specific basis. Conformational sampling is described in terms of ensembles of structures selected from a large pool of conformers. We test this approach, using extensive simulation, to determine how well the fitting of RDCs to reduced conformational ensembles containing few copies of the molecule can correctly reproduce the backbone conformational behavior of the protein. Having established approaches that allow accurate mapping of backbone dihedral angle conformational space from RDCs, we apply these methods to obtain an amino acid specific description of ubiquitin denatured in 8 M urea at pH 2.5. Cross-validation of data not employed in the fit verifies that an ensemble size of 200 structures is appropriate to characterize the highly fluctuating backbone. This approach allows us to identify local conformational sampling properties of urea-unfolded ubiquitin, which shows that the backbone sampling of certain types of charged or polar amino acids, in particular threonine, glutamic acid, and arginine, is affected more strongly by urea binding than amino acids with hydrophobic side chains. In general, the approach presented here establishes robust procedures for the study of all denatured and intrinsically disordered states.

Determination of the free energy landscape of alpha-synuclein using spin label nuclear magnetic resonance measurements

Natively unfolded proteins present a challenge for structure determination because they populate highly heterogeneous ensembles of conformations. A useful source of structural information about these states is provided by paramagnetic relaxation enhancement measurements by nuclear magnetic resonance spectroscopy, from which long-range interatomic distances can be estimated. Here we describe a method for using such distances as restraints in molecular dynamics simulations to obtain a mapping of the free energy landscapes of natively unfolded proteins. We demonstrate the method in the case of alpha-synuclein and validate the results by a comparison with electron transfer measurements. Our findings indicate that our procedure provides an accurate estimate of the relative statistical weights of the different conformations populated by alpha-synuclein in its natively unfolded state.

Structural characterization of alpha-synuclein in an aggregation prone state

The relation of alpha-synuclein (alphaS) aggregation to Parkinson's disease has long been recognized, but the pathogenic species and its molecular properties have yet to be identified. To obtain insight into the properties of alphaS in an aggregation-prone state, we studied the structural properties of alphaS at acidic pH using NMR spectroscopy and computation. NMR demonstrated that alphaS remains natively unfolded at lower pH, but secondary structure propensities were changed in proximity to acidic residues. The ensemble of conformations of alphaS at acidic pH is characterized by a rigidification and compaction of the Asp and Glu-rich C-terminal region, an increased probability for proximity between the NAC-region and the C-terminal region and a lower probability for interactions between the N- and C-terminal regions.

Comparing models of evolution for ordered and disordered proteins

Most models of protein evolution are based upon proteins that form relatively rigid 3D structures. A significant fraction of proteins, the so-called disordered proteins, do not form rigid 3D structures and sample a broad conformational ensemble. Disordered proteins do not typically maintain long-range interactions, so the constraints on their evolution should be different than ordered proteins. To test this hypothesis, we developed and compared models of evolution for disordered and ordered proteins. Substitution matrices were constructed using the sequences of putative homologs for sets of experimentally characterized disordered and ordered proteins. Separate matrices, at three levels of sequence similarity (>85%, 85–60%, and 60–40%), were inferred for each type of protein structure. The substitution matrices for disordered and ordered proteins differed significantly at each level of sequence similarity. The disordered matrices reflected a greater likelihood of evolutionary changes, relative to the ordered matrices, and these changes involved nonconservative substitutions. Glutamic acid and asparagine were interesting exceptions to this result. Important differences between the substitutions that are accepted in disordered proteins relative to ordered proteins were also identified. In general, disordered proteins have fewer evolutionary constraints than ordered proteins. However, some residues like tryptophan and tyrosine are highly conserved in disordered proteins. This is due to their important role in forming protein–protein interfaces. Finally, the amino acid frequencies for disordered proteins, computed during the development of the matrices, were compared with amino acid frequencies for different categories of secondary structure in ordered proteins. The highest correlations were observed between the amino acid frequencies in disordered proteins and the solvent-exposed loops and turns of ordered proteins, supporting an emerging structural model for disordered proteins.

A self-consistent description of the conformational behavior of chemically denatured proteins from NMR and small angle scattering

Characterization of the conformational properties of unfolded proteins is essential for understanding the mechanisms of protein folding and misfolding. This information is also fundamental to determining the relationship between flexibility and function in the highly diverse families of intrinsically disordered proteins. Here we present a self-consistent model of conformational sampling of chemically denatured proteins in agreement with experimental data reporting on long-range distance distributions in unfolded proteins using small-angle x-ray scattering and nuclear magnetic resonance pulse-field gradient-based measurements. We find that standard statistical coil models, selected from folded protein databases with secondary structural elements removed, need to be refined to correct backbone dihedral angle sampling of denatured proteins, although they appear to be appropriate for intrinsically disordered proteins. For denatured proteins, pervasive increases in the sampling of more-extended regions of Ramachandran space {50 degrees Collapse

Key Words Collapse

MESH Headings

• Computer Simulation
• Databases, Protein
• Models, Chemical
• Nuclear Magnetic Resonance, Biomolecular
• Protein Conformation
• Protein Denaturation
• Protein Folding
• Protein Structure, Secondary
• Scattering, Small Angle

Collapse

Grants Collapse

Affiliation(s)

Pau Bernadó Institute for Research in Biomedicine, c/ Baldiri Reixac, Barcelona, Spain
Martin Blackledge Protein Dynamics and Flexibility, Institut de Biologie Structurale, UMR 5075, Commissariat à l'Énergie Atomique-Centre National de la Recherche Scientifique-Université Joseph Fourier, Grenoble, France

Collapse

Improvement and analysis of computational methods for prediction of residual dipolar couplings

We describe a new, computationally efficient method for computing the molecular alignment tensor based on the molecular shape. The increase in speed is achieved by re-expressing the problem as one of numerical integration, rather than a simple uniform sampling (as in the PALES method), and by using a convex hull rather than a detailed representation of the surface of a molecule. This method is applicable to bicelles, PEG/hexanol, and other alignment media that can be modeled by steric restrictions introduced by a planar barrier. This method is used to further explore and compare various representations of protein shape by an equivalent ellipsoid. We also examine the accuracy of the alignment tensor and residual dipolar couplings (RDC) prediction using various ab initio methods. We separately quantify the inaccuracy in RDC prediction caused by the inaccuracy in the orientation and in the magnitude of the alignment tensor, concluding that orientation accuracy is much more important in accurate prediction of RDCs.

Structure and flexibility within proteins as identified through small angle X-ray scattering

Flexibility between domains of proteins is often critical for function. These motions and proteins with large scale flexibility in general are often not readily amenable to conventional structural analysis such as X-ray crystallography, nuclear magnetic resonance spectroscopy (NMR) or electron microscopy. A common evolution of a crystallography project, once a high resolution structure has been determined, is to postulate possible sights of flexibility. Here we describe an analysis tool using relatively inexpensive small angle X-ray scattering (SAXS) measurements to identify flexibility and validate a constructed minimal ensemble of models, which represent highly populated conformations in solution. The resolution of these results is sufficient to address the questions being asked: what kinds of conformations do the domains sample in solution? In our rigid body modeling strategy BILBOMD, molecular dynamics (MD) simulations are used to explore conformational space. A common strategy is to perform the MD simulation on the domains connections at very high temperature, where the additional kinetic energy prevents the molecule from becoming trapped in a local minimum. The MD simulations provide an ensemble of molecular models from which a SAXS curve is calculated and compared to the experimental curve. A genetic algorithm is used to identify the minimal ensemble (minimal ensemble search, MES) required to best fit the experimental data. We demonstrate the use of MES in several model and in four experimental examples.

Paramagnetic relaxation enhancements in unfolded proteins: theory and application to drkN SH3 domain

Site-directed spin labeling in combination with paramagnetic relaxation enhancement (PRE) measurements is one of the most promising techniques for studying unfolded proteins. Since the pioneering work of Gillespie and Shortle (J Mol Biol 1997;268:158), PRE data from unfolded proteins have been interpreted using the theory that was originally developed for rotational spin relaxation. At the same time, it can be readily recognized that the relative motion of the paramagnetic tag attached to the peptide chain and the reporter spin such as (1)H(N) is best described as a translation. With this notion in mind, we developed a number of models for the PRE effect in unfolded proteins: (i) mutual diffusion of the two tethered spheres, (ii) mutual diffusion of the two tethered spheres subject to a harmonic potential, (iii) mutual diffusion of the two tethered spheres subject to a simulated mean-force potential (Smoluchowski equation); (iv) explicit-atom molecular dynamics simulation. The new models were used to predict the dependences of the PRE rates on the (1)H(N) residue number and static magnetic field strength; the results are appreciably different from the Gillespie-Shortle model. At the same time, the Gillespie-Shortle approach is expected to be generally adequate if the goal is to reconstruct the distance distributions between (1)H(N) spins and the paramagnetic center (provided that the characteristic correlation time is known with a reasonable accuracy). The theory has been tested by measuring the PRE rates in three spin-labeled mutants of the drkN SH3 domain in 2M guanidinium chloride. Two modifications introduced into the measurement scheme-using a reference compound to calibrate the signals from the two samples (oxidized and reduced) and using peak volumes instead of intensities to determine the PRE rates-lead to a substantial improvement in the quality of data. The PRE data from the denatured drkN SH3 are mostly consistent with the model of moderately expanded random-coil protein, although part of the data point toward a more compact structure (local hydrophobic cluster). At the same time, the radius of gyration reported by Choy et al. (J Mol Biol 2002;316:101) suggests that the protein is highly expanded. This seemingly contradictory evidence can be reconciled if one assumes that denatured drkN SH3 forms a conformational ensemble that is dominated by extended conformations, yet also contains compact (collapsed) species. Such behavior is apparently more complex than predicted by the model of a random-coil protein in good solvent/poor solvent.

Distribution of conformations sampled by the central amino acid residue in tripeptides inferred from amide I band profiles and NMR scalar coupling constants

The conformational preference of individual amino acid residues in the unfolded state of peptides and proteins is the subject of a continuous debate. Research has mostly been focused on alanine, owing to its abundance in proteins and its relevance for the understanding of helix <----> coil transitions. In the current study, we have analyzed the amide I band profiles of the IR, isotropic and anisotropic Raman, and VCD profiles of trialanine in terms of a conformational model which, for the first time, explicitly considers the entire ensemble of possible conformations rather than representative structures. The distribution function utilized for a satisfactory simulation of the amide I band profiles was found to also reproduce a set of five J coupling constants reported by Graf et al. (Graf, J.; et al. J. Am. Chem. Soc. 2007, 129, 1179). The results of our analysis reveal a PPII fraction of approximately 0.84 for the central alanine residue, which strongly corroborates the notion that alanine has a very high PPII propensity, exceeding the values obtained from restricted coil libraries. We performed a similar analysis for trivaline and found that the dominant fraction of its central residue is a beta-strand. The fraction of the respective distribution is 0.68. The remaining fraction contains contributions from helical and PPII conformations. The results of our analysis enable us to decide on the suitability of force fields used for MD simulations of short alanine-containing peptides. The paper establishes vibrational spectroscopy as a suitable method to explore the energy landscape of amino acid residues.

Quantitative modelfree analysis of urea binding to unfolded ubiquitin using a combination of small angle X-ray and neutron scattering

Characterization of the conformational properties of denatured proteins is essential to our understanding of the molecular basis of protein folding and stability. Here we combine small angle neutron and X-ray scattering to study the interaction of urea with the protein ubiquitin. Comparing coherent intensities scattered at zero angle, and exploiting the scattering densities of H(2)O, D(2)O, ubiquitin, and urea for X-rays and neutrons, we quantitatively determine the number of urea molecules preferentially bound during unfolding of ubiquitin. We find that a pH change from 6.5 to 2.5 triggers recruitment of approximately 20 urea molecules from bulk solution per ubiquitin molecule during the unfolding process.

Probing the urea dependence of residual structure in denatured human alpha-lactalbumin

Backbone (15)N relaxation parameters and (15)N-(1)H(N) residual dipolar couplings (RDCs) have been measured for a variant of human alpha-lactalbumin (alpha-LA) in 4, 6, 8 and 10 M urea. In the alpha-LA variant, the eight cysteine residues in the protein have been replaced by alanines (all-Ala alpha-LA). This protein is a partially folded molten globule at pH 2 and has been shown previously to unfold in a stepwise non-cooperative manner on the addition of urea. (15)N R(2) values in some regions of all-Ala alpha-LA show significant exchange broadening which is reduced as the urea concentration is increased. Experimental RDC data are compared with RDCs predicted from a statistical coil model and with bulkiness, average area buried upon folding and hydrophobicity profiles in order to identify regions of non-random structure. Residues in the regions corresponding to the B, D and C-terminal 3(10) helices in native alpha-LA show R(2) values and RDC data consistent with some non-random structural propensities even at high urea concentrations. Indeed, for residues 101-106 the residual structure persists in 10 M urea and the RDC data suggest that this might include the formation of a turn-like structure. The data presented here allow a detailed characterization of the non-cooperative unfolding of all-Ala alpha-LA at higher concentrations of denaturant and complement previous studies which focused on structural features of the molten globule which is populated at lower concentrations of denaturant.

Prediction of the rotational tumbling time for proteins with disordered segments

For well-structured, rigid proteins, the prediction of rotational tumbling time (tau(c)) using atomic coordinates is reasonably accurate, but is inaccurate for proteins with long unstructured sequences. Under physiological conditions, many proteins contain long disordered segments that play important regulatory roles in fundamental biological events including signal transduction and molecular recognition. Here we describe an ensemble approach to the boundary element method that accurately predicts tau(c) for such proteins by introducing two layers of molecular surfaces whose correlated velocities decay exponentially with distance. Reliable prediction of tau(c) will help to detect intra- and intermolecular interactions and conformational switches between more ordered and less ordered states of the disordered segments. The method has been extensively validated using 12 reference proteins with 14 to 103 disordered residues at the N- and/or C-terminus and has been successfully employed to explain a set of published results on a system that incorporates a conformational switch.

Structural characterization of the natively unfolded N-terminal domain of human c-Src kinase: insights into the role of phosphorylation of the unique domain

The N-terminal regions of the members of Src family of non-receptor protein tyrosine kinases are intrinsically unfolded and contain the maximum sequence divergence among them. In this study, we have addressed the structural characterization by nuclear magnetic resonance of this region of 84 residues that encompasses the SH4 and the unique domains (USrc) of the human c-Src. With this aim, the backbone assignment was performed using (13)C-detected experiments that overcome the spectral resolution problems and the large number of prolines that are typical for intrinsically unfolded proteins. The analysis of the residual dipolar couplings measured for the USrc indicates the presence of a low populated helical structure in the 60-75 region. No long-range contacts between remote fragments of the chain were detected with paramagnetic relaxation enhancement experiments. The structural characterization was extended to two different phosphorylation states of USrc that encompassed three different phosphorylated sites, Ser17, Thr37, and Ser75. The structural and conformational changes upon phosphorylation were monitored through chemical shift perturbations and residual dipolar couplings, indicating that modifications occur at local level and no global rearrangements were apparent. These results suggest a scenario where phosphorylation induces a global electrostatic perturbation that could be involved in the membrane unbinding of c-Src and that could be related with the localization of the enzyme. These observations suggest the unique domain of Src kinases as a source of selectivity and reinforce the relevant role of intrinsically disordered proteins in biological processes.

Structure and disorder in an unfolded state under nondenaturing conditions from ensemble models consistent with a large number of experimental restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include (1)H-(15)N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O(2)-induced (13)C paramagnetic shifts, hydrogen-exchange protection factors, and (15)N R(2) data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu(2+)-Ni(2+) binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative alpha-helical structure and both native-like and nonnative tertiary structure.

ProtSA: a web application for calculating sequence specific protein solvent accessibilities in the unfolded ensemble

Background

The stability of proteins is governed by the heat capacity, enthalpy and entropy changes of folding, which are strongly correlated to the change in solvent accessible surface area experienced by the polypeptide. While the surface exposed in the folded state can be easily determined, accessibilities for the unfolded state at the atomic level cannot be obtained experimentally and are typically estimated using simplistic models of the unfolded ensemble. A web application providing realistic accessibilities of the unfolded ensemble of a given protein at the atomic level will prove useful.

Results

ProtSA, a web application that calculates sequence-specific solvent accessibilities of the unfolded state ensembles of proteins has been developed and made freely available to the scientific community. The input is the amino acid sequence of the protein of interest. ProtSA follows a previously published calculation protocol which uses the Flexible-Meccano algorithm to generate unfolded conformations representative of the unfolded ensemble of the protein, and uses the exact analytical software ALPHASURF to calculate atom solvent accessibilities, which are averaged on the ensemble.

Conclusion

ProtSA is a novel tool for the researcher investigating protein folding energetics. The sequence specific atom accessibilities provided by ProtSA will allow obtaining better estimates of the contribution of the hydrophobic effect to the free energy of folding, will help to refine existing parameterizations of protein folding energetics, and will be useful to understand the influence of point mutations on protein stability.

Modular organization of rabies virus phosphoprotein

A phosphoprotein (P) is found in all viruses of the Mononegavirales order. These proteins form homo-oligomers, fulfil similar roles in the replication cycles of the various viruses, but differ in their length and oligomerization state. Sequence alignments reveal no sequence similarity among proteins from viruses belonging to the same family. Sequence analysis and experimental data show that phosphoproteins from viruses of the Paramyxoviridae contain structured domains alternating with intrinsically disordered regions. Here, we used predictions of disorder of secondary structure, and an analysis of sequence conservation to predict the domain organization of the phosphoprotein from Sendai virus, vesicular stomatitis virus (VSV) and rabies virus (RV P). We devised a new procedure for combining the results from multiple prediction methods and locating the boundaries between disordered regions and structured domains. To validate the proposed modular organization predicted for RV P and to confirm that the putative structured domains correspond to autonomous folding units, we used two-hybrid and biochemical approaches to characterize the properties of several fragments of RV P. We found that both central and C-terminal domains can fold in isolation, that the central domain is the oligomerization domain, and that the C-terminal domain binds to nucleocapsids. Our results suggest a conserved organization of P proteins in the Rhabdoviridae family in concatenated functional domains resembling that of the P proteins in the Paramyxoviridae family.

Biophysical characterization of intrinsically disordered proteins

The challenges associated with the structural characterization of disordered proteins have resulted in the application of a host of biophysical methods to such systems. NMR spectroscopy is perhaps the most readily suited technique for providing high-resolution structural information on disordered protein states in solution. Optical methods, solid state NMR, ESR and X-ray scattering can also provide valuable information regarding the ensemble of conformations sampled by disordered states. Finally, computational studies have begun to assume an increasingly important role in interpreting and extending the impact of experimental data obtained for such systems. This article discusses recent advances in the applications of these methods to intrinsically disordered proteins.

Linking folding and binding

Many cellular proteins are intrinsically disordered and undergo folding, in whole or in part, upon binding to their physiological targets. The past few years have seen an exponential increase in papers describing characterization of intrinsically disordered proteins, both free and bound to targets. Although NMR spectroscopy remains the favored tool, a number of new biophysical techniques are proving exceptionally useful in defining the limits of the conformational ensembles. Advances have been made in prediction of the recognition elements in disordered proteins, in elucidating the kinetics and mechanism of the coupled folding and binding process, and in understanding the role of post-translational modifications in tuning the biological response. Here we review these and other recent advances that are providing new insights into the conformational propensities and interactions of intrinsically disordered proteins and are beginning to reveal general principles underlying their biological functions.

A robust approach for analyzing a heterogeneous structural ensemble

Intrinsically unstructured proteins (IUP) are widespread in eukaryotes and participate in numerous cellular processes, but a structural explanation of the mechanisms they use to recognize and bind their diverse targets has proved elusive. Transcriptional activator domains are one class of IUP that function by recruiting other factors into basal transcription complexes. Transcriptional activator domains are known to use electrostatic interactions for recognition, but it is unclear how this could be accomplished by a structurally heterogeneous ensemble. To investigate this question, we performed principal component analysis on the atomic contact maps of an experimentally restrained ensemble of the human p53 transcriptional activator domain. This analysis revealed that the ensemble is conspicuously nonrandom and permitted a straightforward identification of persistent structural features and their relative probabilities. It was observed that six predominant long-range contacts are combinatorially arranged in 13 clusters of structures. Potential surfaces of the aligned clusters showed that these contacts uniformly organize the negative charges of the highly acidic p53 transcriptional activator domain on one face of the clusters. This observation provides a structural basis for the recruitment of other factors into basal transcription complexes and further supports the hypothesis that the structural ensembles of IUPs are not random and instead have evolved under selection to maintain specific structural features.

Model-independent interpretation of NMR relaxation data for unfolded proteins: the acid-denatured state of ACBP

We have investigated the acid-unfolded state of acyl-coenzyme A binding protein (ACBP) using 15N laboratory frame nuclear magnetic resonance (NMR) relaxation experiments at three magnetic field strengths. The data have been analyzed using standard model-free fitting and models involving distribution of correlation times. In particular, a model-independent method of analysis that does not assume any analytical form for the correlation time distribution is proposed. This method explains correlations between model-free parameters and the analytical distribution parameters found by other authors. The analysis also shows that the relaxation data are consistent with and complementary to information obtained from other parameters, especially secondary chemical shifts and residual dipolar couplings, and strengthens the conclusions of previous observations that three out of the four regions that form helices in the native structure appear to contain residual secondary structure also in the acid-denatured state.

Nucleocapsid structure and function

Measles virus belongs to the Paramyxoviridae family within the Mononegavirales order. Its nonsegmented, single-stranded, negative-sense RNA genome is encapsidated by the nucleoprotein (N) to form a helical nucleocapsid. This ribonucleoproteic complex is the substrate for both transcription and replication. The RNA-dependent RNA polymerase binds to the nucleocapsid template via its co-factor, the phosphoprotein (P). This chapter describes the main structural information available on the nucleoprotein, showing that it consists of a structured core (N(CORE)) and an intrinsically disordered C-terminal domain (N(TAIL)). We propose a model where the dynamic breaking and reforming of the interaction between N(TAIL) and P would allow the polymerase complex (L-P) to cartwheel on the nucleocapsid template. We also propose a model where the flexibility of the disordered N and P domains allows the formation of a tripartite complex (No-P-L) during replication, followed by the delivery of N monomers to the newly synthesized genomic RNA chain. Finally, the functional implications of structural disorder are also discussed in light of the ability of disordered regions to establish interactions with multiple partners, thus leading to multiple biological effects.

NMR: prediction of molecular alignment from structure using the PALES software

Orientational restraints such as residual dipolar couplings promise to overcome many of the problems that traditionally limited liquid-state nuclear magnetic resonance spectroscopy. Recently, we developed methods to predict a molecular alignment tensor and thus residual dipolar couplings for a given molecular structure. This provides many new opportunities for the study of the structure and dynamics of proteins, nucleic acids, oligosaccharides and small molecules. This protocol details the use of the software PALES (Prediction of AlignmEnt from Structure) for prediction of an alignment tensor from a known three-dimensional (3D) coordinate file of a solute. The method is applicable to alignment of molecules in many neutral and charged orienting media and takes into account the molecular shape and 3D charge distribution of the molecule.

Structural biology by NMR: structure, dynamics, and interactions

The function of bio-macromolecules is determined by both their 3D structure and conformational dynamics. These molecules are inherently flexible systems displaying a broad range of dynamics on time-scales from picoseconds to seconds. Nuclear Magnetic Resonance (NMR) spectroscopy has emerged as the method of choice for studying both protein structure and dynamics in solution. Typically, NMR experiments are sensitive both to structural features and to dynamics, and hence the measured data contain information on both. Despite major progress in both experimental approaches and computational methods, obtaining a consistent view of structure and dynamics from experimental NMR data remains a challenge. Molecular dynamics simulations have emerged as an indispensable tool in the analysis of NMR data.

Structural and thermodynamic characterization of T4 lysozyme mutants and the contribution of internal cavities to pressure denaturation

Using small-angle X-ray scattering (SAXS) and tryptophan fluorescence spectroscopy, we have identified multiple compact denatured states of a series of T4 lysozyme mutants that are stabilized by high pressures. Recent studies imply that the mechanism of pressure denaturation is the penetration of water into the protein rather than the transfer of hydrophobic residues into water. To investigate water penetration and the volume change associated with pressure denaturation, we studied the solution behavior of four T4 lysozyme mutants having different cavity volumes at low and neutral pH up to a pressure of 400 MPa (0.1 MPa = 0.9869 atm). At low pH, L99A T4 lysozyme expanded from a compact folded state to a partially unfolded state with a corresponding change in radius of gyration from 17 to 32 A. The volume change upon denaturation correlated well with the total cavity volume, indicating that all of the molecule's major cavities are hydrated with pressure. As a direct comparison to high-pressure crystal structures of L99A T4 lysozyme solved at neutral pH [Collins, M. D., Hummer, G., Quillin, M. L., Matthews, B. W., and Gruner, S. M. (2005) Proc. Natl. Acad. Sci. U.S.A. 102, 16668-16671], pressure denaturation of L99A and the structurally similar L99G/E108V mutant was studied at neutral pH. The pressure-denatured state at neutral pH is even more compact than at low pH, and the small volume changes associated with denaturation suggest that the preferential filling of large cavities is responsible for the compactness of the pressure-denatured state. These results confirm that pressure denaturation is characteristically distinct from thermal or chemical denaturation.

Domain conformation of tau protein studied by solution small-angle X-ray scattering

Tau is one of the two main proteins involved in the pathology of Alzheimer's disease via formation of beta-sheet rich intracellular aggregates named paired helical filaments (PHFs). Given that tau is a natively unfolded protein with no folded core (even upon binding to physiological partners such as microtubules), its structural analysis by high-resolution techniques has been difficult. In this study, employing solution small-angle X-ray scattering from the full length isoforms and from a variety of deletion and point mutants the conformation of tau in solution is structurally characterized. A recently developed ensemble optimization method was employed to generate pools of random models and to select ensembles of coexisting conformations, which fitted simultaneously the scattering data from the full length protein and deletion mutants. The analysis of the structural properties of these selected ensembles allowed us to extract information about residual structure in different domains of the native protein. The short deletion mutants containing the repeat domain (considered the core constituent of the PHFs) are significantly more extended than random coils, suggesting an extended conformation of the repeat domain. The longer tau constructs are comparable in size with the random coils, pointing to long-range contacts between the N- and C-termini compensating for the extension of the repeat domain. Moreover, most of the aggregation-promoting mutants did not show major differences in structure from their wild-type counterparts, indicating that their increased pathological effect is triggered only after an aggregation core has been formed.

The effect of a DeltaK280 mutation on the unfolded state of a microtubule-binding repeat in Tau

Tau is a natively unfolded protein that forms intracellular aggregates in the brains of patients with Alzheimer's disease. To decipher the mechanism underlying the formation of tau aggregates, we developed a novel approach for constructing models of natively unfolded proteins. The method, energy-minima mapping and weighting (EMW), samples local energy minima of subsequences within a natively unfolded protein and then constructs ensembles from these energetically favorable conformations that are consistent with a given set of experimental data. A unique feature of the method is that it does not strive to generate a single ensemble that represents the unfolded state. Instead we construct a number of candidate ensembles, each of which agrees with a given set of experimental constraints, and focus our analysis on local structural features that are present in all of the independently generated ensembles. Using EMW we generated ensembles that are consistent with chemical shift measurements obtained on tau constructs. Thirty models were constructed for the second microtubule binding repeat (MTBR2) in wild-type (WT) tau and a DeltaK280 mutant, which is found in some forms of frontotemporal dementia. By focusing on structural features that are preserved across all ensembles, we find that the aggregation-initiating sequence, PHF6*, prefers an extended conformation in both the WT and DeltaK280 sequences. In addition, we find that residue K280 can adopt a loop/turn conformation in WT MTBR2 and that deletion of this residue, which can adopt nonextended states, leads to an increase in locally extended conformations near the C-terminus of PHF6*. As an increased preference for extended states near the C-terminus of PHF6* may facilitate the propagation of beta-structure downstream from PHF6*, these results explain how a deletion at position 280 can promote the formation of tau aggregates.

Fluorescence characterization of denatured proteins

Characterization of unfolded states, while critical to a complete understanding of protein folding, is inherently difficult due to structural heterogeneity and dynamic interchange between states. The growing body of work focusing on single molecule fluorescence techniques for the study of protein folding, also highlights their potential for studies of unfolded proteins. These methods can obtain conformational information about individual subpopulations of molecules in an ensemble, and measure dynamics without the need for synchronization. The studies highlighted here demonstrate the promise of these techniques for obtaining novel information about unfolded states in vitro and in more physiologically relevant milieu.

Intramolecular cohesion of coils mediated by phenylalanine--glycine motifs in the natively unfolded domain of a nucleoporin

The nuclear pore complex (NPC) provides the sole aqueous conduit for macromolecular exchange between the nucleus and the cytoplasm of cells. Its diffusion conduit contains a size-selective gate formed by a family of NPC proteins that feature large, natively unfolded domains with phenylalanine–glycine repeats (FG domains). These domains of nucleoporins play key roles in establishing the NPC permeability barrier, but little is known about their dynamic structure. Here we used molecular modeling and biophysical techniques to characterize the dynamic ensemble of structures of a representative FG domain from the yeast nucleoporin Nup116. The results showed that its FG motifs function as intramolecular cohesion elements that impart order to the FG domain and compact its ensemble of structures into native premolten globular configurations. At the NPC, the FG motifs of nucleoporins may exert this cohesive effect intermolecularly as well as intramolecularly to form a malleable yet cohesive quaternary structure composed of highly flexible polypeptide chains. Dynamic shifts in the equilibrium or competition between intra- and intermolecular FG motif interactions could facilitate the rapid and reversible structural transitions at the NPC conduit needed to accommodate passing karyopherin–cargo complexes of various shapes and sizes while simultaneously maintaining a size-selective gate against protein diffusion.

The nuclear pore complex is a molecular filter that gates macromolecular exchange between the cytoplasm and the nucleoplasm of cells. It contains a size-selective diffusion barrier at its center composed of proteins named FG nucleoporins. These nucleoporins feature large, structurally disordered domains that are highly decorated with phenylalanine–glycine (FG) sequence motifs. The dynamic structure of these disordered FG domains excludes them from classical structural biology analyses such as X-ray crystallography; thus, new approaches are needed to characterize their shape. Here computational and biophysical approaches were used to elucidate the ensemble of structures adopted by the FG domain of a nucleoporin. The analyses showed that the FG motifs function as intramolecular cohesion elements that compact the shape of the FG domain, forcing it to adopt loosely knit globular configurations that are constantly reconfiguring. Within the nuclear pore complex, dozens of these nucleoporin FG domains may stack as loosely knit globules forming a porous sieve that gates molecular diffusion by size exclusion.

A structure refinement protocol combining NMR residual dipolar couplings and small angle scattering restraints

We present the implementation of a target function based on Small Angle Scattering data (Gabel et al. Eur Biophys J 35(4):313-327, 2006) into the Crystallography and NMR Systems (CNS) and demonstrate its utility in NMR structure calculations by simultaneous application of small angle scattering (SAS) and residual dipolar coupling (RDC) restraints. The efficiency and stability of the approach are demonstrated by reconstructing the structure of a two domain region of the 31 kDa nuclear export factor TAP (TIP-associated protein). Starting with the high resolution X-ray structures of the two individual TAP domains, the translational and orientational domain arrangement is refined simultaneously. We tested the stability of the protocol against variations of the SAS target parameters and the number of RDCs and their uncertainties. The activation of SAS restraints results in an improved translational clustering of the domain positions and lifts part of the fourfold degeneracy of their orientations (associated with a single alignment tensor). The resulting ensemble of structures reflects the conformational space that is consistent with the experimental SAS and RDC data. The SAS target function is computationally very efficient. SAS restraints can be activated at different levels of precision and only a limited SAS angular range is required. When combined with additional data from chemical shift perturbation, paramagnetic relaxation enhancement or mutational analysis the SAS refinement is an efficient approach for defining the topology of multi-domain and/or multimeric biomolecular complexes in solution based on available high resolution structures (NMR or X-ray) of the individual domains.

Calculation of residual dipolar couplings from disordered state ensembles using local alignment

Residual dipolar couplings (RDCs) have been observed in disordered states of several proteins. While their nonuniform values were initially surprising, it has been shown that reasonable approximation of experimental RDCs can be obtained using simple statistical coil models and assuming global alignment of each structure, provided that many thousands of conformers are averaged. Here we show that, by using short local alignment tensors, we can achieve good agreement between experimental and simulated RDCs with far fewer structures than required when using global alignment. This makes the possibility of using RDCs as direct restraints in structural calculations of disordered proteins much more feasible. In addition, it provides insight into the nature of RDCs in disordered states, suggesting that they are primarily reporting on local structure.

Kinetic barriers and the role of topology in protein and RNA folding

This review compares the folding behavior of proteins and RNAs. Topics covered include the role of topology in the determination of folding rates, major folding events including collapse, properties of denatured states, pathway heterogeneity, and the influence of the mode of initiation on the folding pathway.