Improved structural characterizations of the drkN SH3 domain unfolded state suggest a compact ensemble with native-like and non-native structure

Due to their dynamic ensemble nature and a deficiency of experimental restraints, disordered states of proteins are difficult to characterize structurally. Here, we have expanded upon our previous work on the unfolded state of the Drosophila drk N-terminal (drkN) SH3 domain with our program ENSEMBLE, which assigns population weights to pregenerated conformers in order to calculate ensembles of structures whose properties are collectively consistent with experimental measurements. The experimental restraint set has been enlarged with newly measured paramagnetic relaxation enhancements from Cu(2+) bound to an amino terminal Cu(2+)-Ni(2+) binding (ATCUN) motif as well as nuclear Overhauser effect (NOE) and hydrogen exchange data from recent studies. In addition, two new pseudo-energy minimization algorithms have been implemented that have dramatically improved the speed of ENSEMBLE population weight assignment. Finally, we have greatly improved our conformational sampling by utilizing a variety of techniques to generate both random structures and structures that are biased to contain elements of native-like or non-native structure. Although it is not possible to uniquely define a representative structural ensemble, we have been able to assess various properties of the drkN SH3 domain unfolded state by performing ENSEMBLE minimizations of different conformer pools. Specifically, we have found that the experimental restraint set enforces a compact structural distribution that is not consistent with an overall native-like topology but shows preference for local non-native structure in the regions corresponding to the diverging turn and the beta5 strand of the folded state and for local native-like structure in the region corresponding to the beta6 and beta7 strands. We suggest that this approach could be generally useful for the structural characterization of disordered states.     

Ensemble modeling of protein disordered states: Experimental restraint contributions and validation

Disordered states of proteins include the biologically functional intrinsically disordered proteins and the unfolded states of normally folded proteins. In recent years, ensemble‐modeling strategies using various experimental measurements as restraints have emerged as powerful means for structurally characterizing disordered states. However, these methods are still in their infancy compared with the structural determination of folded proteins. Here, we have addressed several issues important to ensemble modeling using our ENSEMBLE methodology. First, we assessed how calculating ensembles containing different numbers of conformers affects their structural properties. We find that larger ensembles have very similar properties to smaller ensembles fit to the same experimental restraints, thus allowing a considerable speed improvement in our calculations. In addition, we analyzed the contributions of different experimental restraints to the structural properties of calculated ensembles, enabling us to make recommendations about the experimental measurements that should be made for optimal ensemble modeling. The effects of different restraints, most significantly from chemical shifts, paramagnetic relaxation enhancements and small‐angle X‐ray scattering, but also from other data, underscore the importance of utilizing multiple sources of experimental data. Finally, we validate our ENSEMBLE methodology using both cross‐validation and synthetic experimental restraints calculated from simulated ensembles. Our results suggest that secondary structure and molecular size distribution can generally be modeled very accurately, whereas the accuracy of calculated tertiary structure is dependent on the number of distance restraints used. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Structural Interpretation of Paramagnetic Relaxation Enhancement-Derived Distances for Disordered Protein States

Paramagnetic relaxation enhancement (PRE) is a powerful technique for studying transient tertiary organizations of unfolded and partially folded proteins. The heterogeneous and dynamic nature of disordered protein states, together with the r⁻⁶ dependence of PRE, presents significant challenges for reliable structural interpretation of PRE-derived distances. Without additional knowledge of accessible conformational substates, ensemble-simulation-based protocols have been used to calculate structure ensembles that appear to be consistent with the PRE distance restraints imposed on the ensemble level with the proper r⁻⁶ weighting. However, rigorous assessment of the reliability of such protocols has been difficult without intimate knowledge of the true nature of disordered protein states. Here we utilize sets of theoretical PRE distances derived from simulated structure ensembles that represent the folded, partially folded and unfolded states of a small protein to investigate the efficacy of ensemble-simulation-based structural interpretation of PRE distances. The results confirm a critical limitation that, due to r⁻⁶ weighting, only one or a few members need to satisfy the distance restraints and the rest of the ensemble are essentially unrestrained. Consequently, calculated structure ensembles will appear artificially heterogeneous no matter whether the PRE distances are derived from the folded, partially unfolded or unfolded state. Furthermore, the nature of the heterogeneous ensembles is largely determined by the protein model employed in structure calculation and reflects little on the true nature of the underlying disordered state. These findings suggest that PRE measurements on disordered protein states alone generally do not contain enough information for a reliable structural interpretation and that the latter will require additional knowledge of accessible conformational substates. Interestingly, when a very large number of PRE measurements is available, faithful structural interpretation might be possible with intermediate ensemble sizes under ideal conditions.     

Backbone 1H and 15N resonance assignments of the N-terminal SH3 domain of drk in folded and unfolded states using enhanced-sensitivity pulsed field gradient NMR techniques

Summary The backbone ¹H and ¹⁵N resonances of the N-terminal SH3 domain of the Drosophila signaling adapter protein, drk, have been assigned. This domain is in slow exchange on the NMR timescale between folded and predominantly unfolded states. Data were collected on both states simultaneously, on samples of the SH3 in near physiological buffer exhibiting an approximately 1:1 ratio of the two states. NMR methods which exploit the chemical shift dispersion of the ¹⁵N resonances of unfolded states and pulsed field gradient water suppression approaches for avoiding saturation and dephasing of amide protons which rapidly exchange with solvent were utilized for the assignment.Abbreviations 2D, 3D two-, three-dimensional - drkN SH3 N-terminal SH3 domain of Drosophila drk - HSQC heteronuclear single-quantum spectroscopy - NOE nuclear Overhauser enhancement - SH3 Src homology domain 3 - TOCSY total correlation spectroscopy     

Exploring the Minimally Frustrated Energy Landscape of Unfolded ACBP

The unfolded state of globular proteins is not well described by a simple statistical coil due to residual structural features, such as secondary structure or transiently formed long-range contacts. The principle of minimal frustration predicts that the unfolded ensemble is biased toward productive regions in the conformational space determined by the native structure. Transient long-range contacts, both native-like and non-native-like, have previously been shown to be present in the unfolded state of the four-helix-bundle protein acyl co-enzyme binding protein (ACBP) as seen from both perturbations in nuclear magnetic resonance (NMR) chemical shifts and structural ensembles generated from NMR paramagnetic relaxation data. To study the nature of the contacts in detail, we used paramagnetic NMR relaxation enhancements, in combination with single-point mutations, to obtain distance constraints for the acid-unfolded ensemble of ACBP. We show that, even in the acid-unfolded state, long-range contacts are specific in nature and single-point mutations affect the free-energy landscape of the unfolded protein. Using this approach, we were able to map out concerted, interconnected, and productive long-range contacts. The correlation between the native-state stability and compactness of the denatured state provides further evidence for native-like contact formation in the denatured state. Overall, these results imply that, even in the earliest stages of folding, ACBP dynamics are governed by native-like contacts on a minimally frustrated energy landscape.     

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.     

Atomic-level characterization of disordered protein ensembles

The roles of unfolded states of proteins in normal folding and in diseases involving aggregation, as well as the prevalence and regulatory functions of intrinsically disordered proteins, have become increasingly recognized. The structural representation of these disordered states as ensembles of interconverting conformers can therefore provide critical insights. Experimental methods can be used to probe ensemble-averaged structural properties of disordered states and computational approaches generate representative ensembles of conformers using experimental restraints. In particular, NMR and small-angle X-ray scattering provide quantitative data that can readily be incorporated into calculations. These techniques have gleaned structural information about denatured, unfolded and intrinsically disordered proteins. The use of experimental data in different computational approaches, including ensemble molecular dynamics simulations and algorithms that assign populations to pregenerated conformers, has highlighted the presence of both local and long-range structure, and the occurrence of native-like and non-native interactions in unfolded and denatured states. Analysis of the resulting ensembles has suggested important implications of this fluctuating structure for folding, aggregation and binding.     

Specific 15N, NH correlations for residues in15 N, 13C and fractionally deuterated proteins that immediately follow methyl-containing amino acids

A triple-resonance pulse scheme is described which records¹⁵N, NH correlations of residues that immediately follow amethyl-containing amino acid. The experiment makes use of a¹⁵N, ¹³C and fractionally deuterated proteinsample and selects for CH₂D methyl types. The experiment isthus useful in the early stages of the sequential assignment process as wellas for the confirmation of backbone ¹⁵N, NH chemical shiftassignments at later stages of data analysis. A simple modification of thesequence also allows the measurement of methyl side-chain dynamics. This isparticularly useful for studying side-chain dynamic properties in partiallyunfolded and unfolded proteins where the resolution of aliphatic carbon andproton chemical shifts is limited compared to that of amide nitrogens.     

α-Lactalbumin, Engineered to be Nonnative and Inactive, Kills Tumor Cells when in Complex with Oleic Acid: A New Biological Function Resulting from Partial Unfolding

HAMLET (human alpha-lactalbumin made lethal to tumor cells) is a tumoricidal complex consisting of partially unfolded protein and fatty acid and was first identified in casein fractions of human breast milk. The complex can be produced from its pure components through a modified chromatographic procedure where preapplied oleic acid binds with partially unfolded α-lactalbumin on the stationary phase in situ. Because native α-lactalbumin itself cannot trigger cell death, HAMLET's remarkable tumor-selective cytotoxicity has been strongly correlated with the conformational change of the protein upon forming the complex, but whether a recovery to the native state subsequently occurs upon entering the tumor cell is yet unclear. To this end, we utilize a recombinant variant of human α-lactalbumin in which all eight cysteine residues are substituted for alanines (rHLA^all-Ala), rendering the protein nonnative and biologically inactive under all conditions. The HAMLET analogue formed from the complex of rHLA^all-Ala and oleic acid (rHLA^all-Ala-OA) exhibited equivalent strong tumoricidal activity against lymphoma and carcinoma cell lines and was shown to accumulate within the nuclei of tumor cells, thus reproducing the cellular trafficking pattern of HAMLET. In contrast, the fatty acid-free rHLA^all-Ala protein associated with the tumor cell surface but was not internalized and lacked any cytotoxic activity. Structurally, whereas HAMLET exhibited some residual native character in terms of NMR chemical shift dispersion, rHLA^all-Ala-OA showed significant differences to HAMLET and, in fact, was found to be devoid of any tertiary packing. The results identify α-lactalbumin as a protein with strikingly different functions in the native and partially unfolded states. We posit that partial unfolding offers another significant route of functional diversification for proteins within the cell.     

Insight into polyproline II helical bundle stability in an antifreeze protein denatured state

The use of polyproline II (PPII) helices in protein design is currently hindered by limitations in our understanding of their conformational stability and folding. Recent studies of the snow flea antifreeze protein (sfAFP), a useful model system composed of six PPII helices, suggested that a low denatured state entropy contributes to folding thermodynamics. Here, circular dichroism spectroscopy revealed minor populations of PPII like conformers at low temperature. To get atomic level information on the conformational ensemble and entropy of the reduced, denatured state of sfAFP, we have analyzed its chemical shifts and {¹H}-¹⁵N relaxation parameters by NMR spectroscopy at four experimental conditions. No significant populations of stable secondary structure were detected. The stiffening of certain N-terminal residues at neutral versus acidic pH and shifted pK_a values leads us to suggest that favorable charge-charge interactions could bias the conformational ensemble to favor the formation the C1-C28 disulfide bond during nascent folding, although no evidence for preferred contacts between these positions was detected by paramagnetic relaxation enhancement under denaturing conditions. Despite a high content of flexible glycine residues, the mobility of the sfAFP denatured ensemble is similar for denatured α/β proteins both on fast ps/ns as well as slower μs/ms timescales. These results are in line with a conformational entropy in the denatured ensemble resembling that of typical proteins and suggest that new structures based on PPII helical bundles should be amenable to protein design.     

Charge sequence coding in statistical modeling of unfolded proteins

Unfolded proteins recently attracted attention due to accumulation of experimental evidences for their significant role in different life processes. Modeling of electrostatic interactions (EI) in unfolded state of proteins is becoming increasingly important as well. In this paper, we stress on the importance of how the sequence of charged residues of a given protein is incorporated into the models for calculation of EI in the unfolded state. On the basis of the distributions of distances between titratable sites of charged residues calculated for polypeptide chains of various compositions, it was found that the distance distribution for a pair of residues, located close to each other along the sequence of a protein, depends on what residues constitute the pair in question. It was concluded that the consideration of these residue-specific distributions is essential for a statistical model to be accurate from the physical point of view. It was suggested that use of distance intervals in the spherical model of unfolded proteins accounts better for the charge sequence than the set of single distance values. This was illustrated by comparison of the pK values of the titratable groups of the unfolded N-terminal SH3 domain of the Drosophila protein drk to the available experimental data.     

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 ¹³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.     

Characterization of conformational and dynamic properties of natively unfolded human and mouse alpha-synuclein ensembles by NMR: implication for aggregation

Conversion of human α-synuclein (aS) from the free soluble state to the insoluble fibrillar state has been implicated in the etiology of Parkinson's disease. Human aS is highly homologous in amino acid sequence to mouse aS, which contains seven substitutions including the A53T that has been linked to familial Parkinson's disease, and including five substitutions in the C-terminal region. It has been shown that the rate of fibrillation is highly dependent on the exact sequence of the protein, and mouse aS is reported to aggregate more rapidly than human aS in vitro. Nuclear magnetic resonance experiments of mouse and human aS at supercooled temperatures (263 K) are used to understand the effect of sequence on conformational fluctuations in the disordered ensembles and to relate these to differences in propensities to aggregate. We show that both aS are natively unfolded at low temperature with different propensities to secondary structure, backbone dynamics and long-range contacts across the protein. Mouse aS exhibits a higher propensity to helical conformation around the C-terminal substitutions as well as the loss of transient long-range contacts from the C- to the N-terminal end and hydrophobic central regions of the protein relative to human aS. Lack of back-folding from the C-terminal end of mouse aS exposes the N-terminal region, which is shown, by ¹⁵N relaxation experiments, to be very restricted in mobility relative to human aS. We propose that the restricted mobility in the N-terminal region may arise from transient interchain interactions, suggesting that the N-terminal KTK(E/Q)GV repeats may serve as initiation sites for aggregation in mouse aS. These transient interchain interactions coupled with a non-Aβ amyloid component (NAC) region that is both more exposed and has a higher propensity to β structure may accelerate the rate of fibril formation of aS.     

Single-Molecule Studies of the Im7 Folding Landscape

Under appropriate conditions, the four-helical Im7 (immunity protein 7) folds from an ensemble of unfolded conformers to a highly compact native state via an on-pathway intermediate. Here, we investigate the unfolded, intermediate, and native states populated during folding using diffusion single-pair fluorescence resonance energy transfer by measuring the efficiency of energy transfer (or proximity or P ratio) between pairs of fluorophores introduced into the side chains of cysteine residues placed in the center of helices 1 and 4, 1 and 3, or 2 and 4. We show that while the native states of each variant give rise to a single narrow distribution with high P values, the distributions of the intermediates trapped at equilibrium (denoted I^eqm) are fitted by two Gaussian distributions. Modulation of the folding conditions from those that stabilize the intermediate to those that destabilize the intermediate enabled the distribution of lower P value to be assigned to the population of the unfolded ensemble in equilibrium with the intermediate state. The reduced stability of the I^eqm variants allowed analysis of the effect of denaturant concentration on the compaction and breadth of the unfolded state ensemble to be quantified from 0 to 6 M urea. Significant compaction is observed as the concentration of urea is decreased in both the presence and absence of sodium sulfate, as previously reported for a variety of proteins. In the presence of Na₂SO₄ in 0 M urea, the P value of the unfolded state ensemble approaches that of the native state. Concurrent with compaction, the ensemble displays increased peak width of P values, possibly reflecting a reduction in the rate of conformational exchange among iso-energetic unfolded, but compact conformations. The results provide new insights into the initial stages of folding of Im7 and suggest that the unfolded state is highly conformationally constrained at the outset of folding.     

Ensemble models of proteins and protein domains based on distance distribution restraints

Conformational ensembles of intrinsically disordered peptide chains are not fully determined by experimental observations. Uncertainty due to lack of experimental restraints and due to intrinsic disorder can be distinguished if distance distributions restraints are available. Such restraints can be obtained from pulsed dipolar electron paramagnetic resonance (EPR) spectroscopy applied to pairs of spin labels. Here, we introduce a Monte Carlo approach for generating conformational ensembles that are consistent with a set of distance distribution restraints, backbone dihedral angle statistics in known protein structures, and optionally, secondary structure propensities or membrane immersion depths. The approach is tested with simulated restraints for a terminal and an internal loop and for a protein with 69 residues by using sets of sparse restraints for underlying well‐defined conformations and for published ensembles of a premolten globule‐like and a coil‐like intrinsically disordered protein. Proteins 2016; 84:544–560. © 2016 Wiley Periodicals, Inc.     

Determination of the Individual Roles of the Linker Residues in the Interdomain Motions of Calmodulin Using NMR Chemical Shifts

Many protein molecules are formed by two or more domains whose structures and dynamics are closely related to their biological functions. It is thus important to develop methods to determine the structural properties of these multidomain proteins. Here, we characterize the interdomain motions in the calcium-bound state of calmodulin (Ca^2 +-CaM) using NMR chemical shifts as replica-averaged structural restraints in molecular dynamics simulations. We find that the conformational fluctuations of the interdomain linker, which are largely responsible for the overall interdomain motions of CaM, can be well described by exploiting the information provided by chemical shifts. We thus identify 10 residues in the interdomain linker region that change their conformations upon substrate binding. Five of these residues (Met76, Lys77, Thr79, Asp80 and Ser81) are highly flexible and cover the range of conformations observed in the substrate-bound state, while the remaining five (Arg74, Lys75, Asp78, Glu82 and Glu83) are much more rigid and do not populate conformations typical of the substrate-bound form. The ensemble of conformations representing the Ca^2 +-CaM state obtained in this study is in good agreement with residual dipolar coupling, paramagnetic resonance enhancement, small-angle X-ray scattering and fluorescence resonance energy transfer measurements, which were not used as restraints in the calculations. These results provide initial evidence that chemical shifts can be used to characterize the conformational fluctuations of multidomain proteins.     

Pressure-dependent <Superscript>13</Superscript>C chemical shifts in proteins: origins and applications

Pressure-dependent ¹³C chemical shifts have been measured for aliphatic carbons in barnase and Protein G. Up to 200 MPa (2 kbar), most shift changes are linear, demonstrating pressure-independent compressibilities. CH₃, CH₂ and CH carbon shifts change on average by +0.23, −0.09 and −0.18 ppm, respectively, due to a combination of bond shortening and changes in bond angles, the latter matching one explanation for the γ-gauche effect. In addition, there is a residue-specific component, arising from both local compression and conformational change. To assess the relative magnitudes of these effects, residue-specific shift changes for protein G were converted into structural restraints and used to calculate the change in structure with pressure, using a genetic algorithm to convert shift changes into dihedral angle restraints. The results demonstrate that residual ¹³Cα shifts are dominated by dihedral angle changes and can be used to calculate structural change, whereas ¹³Cβ shifts retain significant dependence on local compression, making them less useful as structural restraints.