Hydrodynamic radii of native and denatured proteins measured by pulse field gradient NMR techniques

Pulse field gradient NMR methods have been used to determine the effective hydrodynamic radii of a range of native and nonnative protein conformations. From these experimental data, empirical relationships between the measured hydrodynamic radius (R(h)) and the number of residues in the polypeptide chain (N) have been established; for native folded proteins R(h) = 4.75N (0.29)A and for highly denatured states R(h) = 2.21N (0.57)A. Predictions from these equations agree well with experimental data from dynamic light scattering and small-angle X-ray or neutron scattering studies reported in the literature for proteins ranging in size from 58 to 760 amino acid residues. The predicted values of the hydrodynamic radii provide a framework that can be used to analyze the conformational properties of a range of nonnative states of proteins. Several examples are given here to illustrate this approach including data for partially structured molten globule states and for proteins that are unfolded but biologically active under physiological conditions. These reveal evidence for significant coupling between local and global features of the conformational ensembles adopted in such states. In particular, the effective dimensions of the polypeptide chain are found to depend significantly on the level of persistence of regions of secondary structure or features such as hydrophobic clusters within a conformational ensemble.     

Dihedral angle preferences of amino acid residues forming various non-local interactions in proteins

In theory, a polypeptide chain can adopt a vast number of conformations, each corresponding to a set of backbone rotation angles. Many of these conformations are excluded due to steric overlaps. Ramachandran and coworkers were the first to look into this problem by plotting backbone dihedral angles in a two-dimensional plot. The conformational space in the Ramachandran map is further refined by considering the energetic contributions of various non-bonded interactions. Alternatively, the conformation adopted by a polypeptide chain may also be examined by investigating interactions between the residues. Since the Ramachandran map essentially focuses on local interactions (residues closer in sequence), out of interest, we have analyzed the dihedral angle preferences of residues that make non-local interactions (residues far away in sequence and closer in space) in the folded structures of proteins. The non-local interactions have been grouped into different types such as hydrogen bond, van der Waals interactions between hydrophobic groups, ion pairs (salt bridges), and ππ-stacking interactions. The results show the propensity of amino acid residues in proteins forming local and non-local interactions. Our results point to the vital role of different types of non-local interactions and their effect on dihedral angles in forming secondary and tertiary structural elements to adopt their native fold.     

Intrinsically Disordered Protein Exhibits Both Compaction and Expansion under Macromolecular Crowding

Conformational malleability allows intrinsically disordered proteins (IDPs) to respond agilely to their environments, such as nonspecifically interacting with in vivo bystander macromolecules (or crowders). Previous studies have emphasized conformational compaction of IDPs due to steric repulsion by macromolecular crowders, but effects of soft attraction are largely unexplored. Here we studied the conformational ensembles of the IDP FlgM in both polymer and protein crowders by small-angle neutron scattering. As crowder concentrations increased, the mean radius of gyration of FlgM first decreased but then exhibited an uptick. Ensemble optimization modeling indicated that FlgM conformations under protein crowding segregated into two distinct populations, one compacted and one extended. Coarse-grained simulations showed that compacted conformers fit into an interstitial void and occasionally bind to a surrounding crowder, whereas extended conformers snake through interstitial crevices and bind multiple crowders simultaneously. Crowder-induced conformational segregation may facilitate various cellular functions of IDPs.     

Kinetics of Contact Formation and End-to-End Distance Distributions of Swollen Disordered Peptides

Unstructured polypeptide chains are subject to various degrees of swelling or compaction depending on the combination of solvent condition and amino acid sequence. Highly denatured proteins generally behave like random-coils with excluded volume repulsion, whereas in aqueous buffer more compact conformations have been observed for the low-populated unfolded state of globular proteins as well as for naturally disordered sequences. To quantitatively account for the different mechanisms inducing the swelling of polypeptides, we have examined three 14-residues peptides in aqueous buffer and in denaturant solutions, including the well characterized AGQ repeat as a reference and two variants, in which we have successively introduced charged side chains and removed the glycines. Quenching of the triplet state of tryptophan by close contact with cysteine has been used in conjunction with Förster resonance energy transfer to study the equilibrium and kinetic properties of the peptide chains. The experiments enable accessing end-to-end root mean-square distance, probability of end-to-end contact formation and intrachain diffusion coefficient. The data can be coherently interpreted on the basis of a simple chain model with backbone angles obtained from a library of coil segments of proteins and hard sphere repulsion at each C_α position. In buffered water, we find that introducing charges in a glycine-rich sequence induces a mild chain swelling and a significant speed-up of the intrachain dynamics, whereas the removal of the glycines results in almost a two-fold increase of the chain volume and a drastic slowing down. In denaturants we observe a pronounced swelling of all the chains, with significant differences between the effect of urea and guanidinium chloride.     

Specificity of the initial collapse in the folding of the cold shock protein

The two-state folding reaction of the cold shock protein from Bacillus caldolyticus (Bc-Csp) is preceded by a rapid chain collapse. A fast shortening of intra-protein distances was revealed by F?rster resonance energy transfer (FRET) measurements with protein variants that carried individual pairs of donor and acceptor chromophores at various positions along the polypeptide chain. Here we investigated the specificity of this rapid compaction. Energy transfer experiments that probed the stretching of strand beta2 and the close approach between the strands beta1 and beta2 revealed that the beta1-beta2 hairpin is barely formed in the collapsed form, although it is native-like in the folding transition state of Bc-Csp. The time course of the collapse could not be resolved by pressure or temperature jump experiments, indicating that the collapsed and extended forms are not separated by an energy barrier. The co-solute (NH4)2SO4 stabilizes both native Bc-Csp and the collapsed form, which suggests that the large hydrated SO4(2-) ions are excluded from the surface of the collapsed form in a similar fashion as they are excluded from folded Bc-Csp. Ethylene glycol increases the stability of proteins because it is excluded preferentially from the backbone, which is accessible in the unfolded state. The collapsed form of Bc-Csp resembles the unfolded form in its interaction with ethylene glycol, suggesting that in the collapsed form the backbone is still accessible to water and small molecules. Our results thus rule out that the collapsed form is a folding intermediate with native-like chain topology. It is better described as a mixture of compact conformations that belong to the unfolded state ensemble. However, some of its structural elements are reminiscent of the native protein.     

Non-Sequence-Specific Interactions Can Account for the Compaction of Proteins Unfolded under “Native” Conditions

Proteins unfolded by high concentrations of chemical denaturants adopt expanded, largely structure-free ensembles of conformations that are well approximated as random coils. In contrast, globular proteins unfolded under less denaturing conditions (via mutations, or transiently unfolded after a rapid jump to native conditions) and molten globules (arising due to mutations or cosolvents) are often compact. Here we explore the origins of this compaction using a truncated equilibrium-unfolded variant of the 57-residue FynSH3 domain. As monitored by far-UV circular dichroism, NMR spectroscopy, and hydrogen-exchange kinetics, CΔ4 (a 4-residue carboxy-terminal deletion variant of FynSH3) appears to be largely unfolded even in the absence of denaturant. Nevertheless, CΔ4 is quite compact under these conditions, with a hydrodynamic radius only slightly larger than that of the native protein. In order to understand the origins of this molten-globule-like compaction, we have characterized a random sequence polypeptide of identical amino acid composition to CΔ4. Notably, we find that the hydrodynamic radius of this random sequence polypeptide also approaches that of the native protein. Thus, while native-like interactions may contribute to the formation of compact “unfolded” states, it appears that non-sequence-specific monomer-monomer interactions can also account for the dramatic compaction observed for molten globules and the “physiological” unfolded state.     

Reconciling observations of sequence-specific conformational propensities with the generic polymeric behavior of denatured proteins

Radii of gyration of denatured proteins vary with chain length and are insensitive to details of amino acid sequence. Observations of sequence independence in polymeric properties conflict with results from spectroscopic experiments, which suggest the presence of sequence-specific residual structure in denatured states. Can we reconcile the two apparently conflicting sets of observations? To answer this question, we need knowledge of the ensemble of conformations accessible to proteins in good solvents. The excluded-volume limit provides an ideal mimic of polymers in good solvents. Therefore, we attempt to solve the "reconciliation problem" by simulating conformational ensembles accessible to peptides and proteins in the excluded-volume limit. Analysis of these ensembles for a variety of polypeptide sequences leads to results that are consistent with experimental observations of sequence-specific conformational preferences in short peptides and the scaling behavior of polymeric quantities for denatured proteins. Reconciliation in the excluded-volume limit comes about due to a tug of war between two factors, namely, minimization of steric overlap and the competing effects of conformational entropy. Minimization of steric overlap promotes chain stretching and leads to experimentally observed sequence-dependent preferences for locally extended segments such as polyproline II helices, beta-strands, and very short stretches of alpha-helix. Conformational entropy opposes chain stretching, and the calculated persistence length for sequence-dependent conformational preferences is less than five amino acids. This estimate does not vary with amino acid sequence. The short persistence lengths lead directly to experimental observations of generic sequence-independent behavior of radii of gyration for denatured proteins.     

Native and nonnative conformational preferences in the urea-unfolded state of barstar

The refolding of barstar from its urea-unfolded state has been studied extensively using various spectroscopic probes and real-time NMR, which provide global and residue-specific information, respectively, about the folding process. Here, a preliminary structural characterization by NMR of barstar in 8 M urea has been carried out at pH 6.5 and 25 degrees C. Complete backbone resonance assignments of the urea-unfolded protein were obtained using the recently developed three-dimensional NMR techniques of HNN and HN(C)N. The conformational propensities of the polypeptide backbone in the presence of 8 M urea have been estimated by examining deviations of secondary chemical shifts from random coil values. For some residues that belong to helices in native barstar, 13C(alpha) and 13CO secondary shifts show positive deviations in the urea-unfolded state, indicating that these residues have propensities toward helical conformations. These residues are, however, juxtaposed by residues that display negative deviations indicative of propensities toward extended conformations. Thus, segments that are helical in native barstar are unlikely to preferentially populate the helical conformation in the unfolded state. Similarly, residues belonging to beta-strands 1 and 2 of native barstar do not appear to show any conformational preferences in the unfolded state. On the other hand, residues belonging to the beta-strand 3 segment show weak nonnative helical conformational preferences in the unfolded state, indicating that this segment may possess a weak preference for populating a helical conformation in the unfolded state.     

Observation of residual dipolar couplings in short peptides

Residual dipolar couplings provide information on the orientation of individual bond vectors with respect to a unique set of molecular axes. We report that short peptides from 2 to 15 amino acids in length of arbitrary sequence exhibit a modest range of residual dipolar couplings when aligned in either strained polyacrylamide gels or alkyl-PEG bicelles. The absence of significant line broadening in gels suggests peptides align predominantly through steric interactions with the polyacrylamide matrix. However, broadening of NMR lines for a subset of residues aligned in bicelles indicates some peptides bind weakly to these lipid disks, yet a weak negative correlation between the couplings measured in gels and bicelles is consistent with steric hindrance playing a role in both media. The observation of dipolar couplings for peptides of length 10-15 suggests the statistical segment lengths of polypeptide chains must often be >10-15 residues, with data from denatured proteins indicating even larger values. Presumably, local side-chain backbone interactions severely restrict chain flexibility, with the cumulative effect of many such restrictions giving rise to biases in chain direction that may persist for the entire length of a protein chain. Comparison of experimental dipolar couplings for peptides with couplings calculated for ensembles of conformations generated by molecular dynamics should permit evaluation of the accuracy of molecular mechanics potentials in reproducing sequence-specific preferences for phi and psi angles.     

Distance-Based Metrics for Comparing Conformational Ensembles of Intrinsically Disordered Proteins

Intrinsically disordered proteins are proteins whose native functional states represent ensembles of highly diverse conformations. Such ensembles are a challenge for quantitative structure comparisons because their conformational diversity precludes optimal superimposition of the atomic coordinates necessary for deriving common similarity measures such as the root mean-square deviation of these coordinates. Here, we introduce superimposition-free metrics that are based on computing matrices of the Cα-Cα distance distributions within ensembles and comparing these matrices between ensembles. Differences between two matrices yield information on the similarity between specific regions of the polypeptide, whereas the global structural similarity is captured by the root mean-square difference between the medians of the Cα-Cα distance distributions of two ensembles. Together, our metrics enable rigorous investigations of structure-function relationships in conformational ensembles of intrinsically disordered proteins derived using experimental restraints or by molecular simulations and for proteins containing both structured and disordered regions.     

Small-angle X-ray scattering of reduced ribonuclease A: effects of solution conditions and comparisons with a computational model of unfolded proteins

The disulfide-reduced form of bovine ribonuclease A, with the Cys thiols irreversibly blocked, was characterized by small-angle x-ray scattering. To help resolve the conflicting results and interpretations from previous studies of this model unfolded protein, we measured scattering profiles using a range of solution conditions and compared them with the profiles predicted by a computational model for a random-coil polypeptide. Analysis of the simulated and experimental profiles reveals that scattering intensities at intermediate angles, corresponding to interatomic distances in the range of 5-20 Å, are particularly sensitive to changes in solvation and can be used to assess the internal scaling behavior of the polypeptide chain, expressed as a mass fractal dimension, D_m. This region of the scattering curve is also much less sensitive to experimental artifacts than is the very small angle regime (the Guinier region) that has been more typically used to characterize unfolded proteins. The experimental small-angle x-ray scattering profiles closely matched those predicted by the computational model assuming relatively small solvation energies. The scaling behavior of the polypeptide approaches that of a well-solvated polymer under conditions where it has a large net charge and at high urea concentrations. At lower urea concentrations and neutral pH, the behavior of the chain approaches that expected for θ-conditions, where the effects of slightly unfavorable interactions with solvent balance those of excluded volume, leading to scaling behavior comparable to that of an idealized random walk chain. Though detectable, the shift toward more compact conformations at lower urea concentrations does not correspond to a transition to a globule state and is associated with little or no reduction in conformational entropy. This type of collapse, therefore, is unlikely to greatly reduce the conformational search for the native state.     

Tensile Mechanics of Alanine-Based Helical Polypeptide: Force Spectroscopy versus Computer Simulations

In nature, an α-helix is commonly used to build thermodynamically stable and mechanically rigid protein conformations. In view of growing interest in the mechanical rigidity of proteins, we measured the tensile profile of an alanine-based α-helical polypeptide on an atomic-force microscope to investigate the basic mechanics of helix extension with minimal interference from side-chain interactions. The peptide was extended to its maximum contour length with much less force than in reported cases of poly-L-Glu or poly-L-Lys, indicating that chain stiffness strongly depended on the physicochemical properties of side chains, such as their bulkiness. The low tensile-force extension originated presumably in locally unfolded parts because of spontaneous structural fluctuations. In 50% trifluoroethanol, the well-known helix-promoting agent, the rigidity of the sample polypeptide was markedly increased. Computer simulations of the peptide-stretching process showed that a majority of constituent residues underwent a transition from an α-helical to an extended conformation by overcoming an energy barrier around ψ ∼0° on the Ramachandran plot. The observed lability of an isolated helix signified the biological importance of the lateral bundling of helices to maintain a rigid protein structure.     

Dynamics of unfolded polypeptide chains in crowded environment studied by fluorescence correlation spectroscopy

Proteins have evolved to fold and function within a cellular environment that is characterized by high macromolecular content. The earliest step of protein folding represents intrachain contact formation of amino acid residues within an unfolded polypeptide chain. It has been proposed that macromolecular crowding can have significant effects on rates and equilibria of biomolecular processes. However, the kinetic consequences on intrachain diffusion of polypeptides have not been tested experimentally, yet. Here, we demonstrate that selective fluorescence quenching of the oxazine fluorophore MR121 by the amino acid tryptophan (Trp) in combination with fast fluorescence correlation spectroscopy (FCS) can be used to monitor end-to-end contact formation rates of unfolded polypeptide chains. MR121 and Trp were incorporated at the terminal ends of polypeptides consisting of repetitive units of glycine (G) and serine (S) residues. End-to-end contact formation and dissociation result in "off" and "on" switching of MR121 fluorescence and underlying kinetics can be revealed in FCS experiments with nanosecond time resolution. We revisit previous experimental studies concerning the dependence of end-to-end contact formation rates on polypeptide chain length, showing that kinetics can be described by Gaussian chain theory. We further investigate effects of solvent viscosity and temperature on contact formation rates demonstrating that intrachain diffusion represents a purely diffusive, entropy-controlled process. Finally, we study the influence of macromolecular crowding on polypeptide chain dynamics. The data presented demonstrate that intrachain diffusion is fast in spite of hindered diffusion caused by repulsive interactions with macromolecules. Findings can be explained by effects of excluded volume reducing chain entropy and therefore accelerating the loop search process. Our results suggest that within a cellular environment the early formation of structural elements in unfolded proteins can still proceed quite efficiently in spite of hindered diffusion caused by high macromolecular content.     

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.     

Dynamics in the unfolded state of beta2-microglobulin studied by NMR

Many proteins form amyloid-like fibrils in vitro under conditions that favour the population of partially folded conformations or denatured state ensembles. Characterising the structural and dynamic properties of these states is crucial towards understanding the mechanisms of self-assembly in amyloidosis. The aggregation of beta2-microglobulin (beta2m) into amyloid fibrils in vivo occurs in the condition known as dialysis-related amyloidosis (DRA) and the protein has been shown to form amyloid-like fibrils under acidic conditions in vitro. We have used a number of 1H-15N nuclear magnetic resonance (NMR) experiments in conjunction with site-directed mutagenesis to study the acid-unfolded state of beta2m. 15N NMR transverse relaxation experiments reveal that the acid-denatured ensemble, although predominantly unfolded at the N and C termini, contains substantial non-native structure in the central region of the polypeptide chain, stabilised by long-range interactions between aromatic residues and by the single disulphide bond. Relaxation dispersion studies indicate that the acid-unfolded ensemble involves two or more distinct species in conformational equilibrium on the micro- to millisecond time-scale. One of these species appears to be hydrophobically collapsed, as mutations in an aromatic-rich region of the protein, including residues that are solvent-exposed in the native protein, disrupt this structure and cause a consequent decrease in the population of this conformer. Thus, acid-unfolded beta2m consists of a heterogeneous ensemble of rapidly fluctuating species, some of which contain stable, non-native hydrophobic clusters. Given that amyloid assembly of beta2m proceeds with lag kinetics under the conditions of this study, a rarely populated species such as a conformer with non-native aromatic clustering could be key to the initiation of amyloidosis.     

Characterization of two distinct beta2-microglobulin unfolding intermediates that may lead to amyloid fibrils of different morphology

The self-assembly of beta(2)-microglobulin into fibrils leads to dialysis-related amyloidosis. pH-mediated partial unfolding is required for the formation of the amyloidogenic intermediate that then self-assembles into amyloid fibrils. Two partially folded intermediates of beta(2)-microglobulin have been identified experimentally and linked to the formation of fibrils of distinct morphology, yet it remains difficult to characterize these partially unfolded states at high resolution using experimental approaches. Consequently, we have performed molecular dynamics simulations at neutral and low pH to determine the structures of these partially unfolded amyloidogenic intermediates. In the low-pH simulations, we observed the formation of alpha-sheet structure, which was first proposed by Pauling and Corey. Multiple simulations were performed, and two distinct intermediate state ensembles were identified that may account for the different fibril morphologies. The predominant early unfolding intermediate was nativelike in structure, in agreement with previous NMR studies. The late unfolding intermediate was significantly disordered, but it maintained an extended elongated structure, with hydrophobic clusters and residual alpha-extended chain strands in specific regions of the sequence that map to amyloidogenic peptides. We propose that the formation of alpha-sheet facilitates self-assembly into partially unfolded prefibrillar amyloidogenic intermediates.     

Kinetics of internal-loop formation in polypeptide chains: a simulation study

The speed of simple diffusional motions, such as the formation of loops in the polypeptide chain, places one physical limit on the speed of protein folding. Many experimental studies have explored the kinetics of formation of end-to-end loops in polypeptide chains; however, protein folding more often requires the formation of contacts between interior points on the chain. One expects that, for loops of fixed contour length, interior loops will form more slowly than end-to-end loops, owing to the additional excluded volume associated with the "tails". We estimate the magnitude of this effect by generating ensembles of randomly coiled, freely jointed chains, and then using the theory of Szabo, Schulten, and Schulten to calculate the corresponding contact formation rates for these ensembles. Adding just a few residues, to convert an end-to-end loop to an internal loop, sharply decreases the contact rate. Surprisingly, the relative change in rate increases for a longer loop; sufficiently long tails, however, actually reverse the effect and accelerate loop formation slightly. Our results show that excluded volume effects in real, full-length polypeptides may cause the rates of loop formation during folding to depart significantly from the values derived from recent loop-formation experiments on short peptides.