Search

Intrinsically disordered proteins (IDPs) have fluctuating heterogeneous conformations, which makes structural characterization challenging. Although challenging, characterizing the conformational ensembles of IDPs is of great interest, since their conformational ensembles are the link between their sequences and functions. An accurate description of IDP conformational ensembles depends crucially on the amount and quality of the experimental data, how it is integrated, and if it supports a consistent structural picture. We used integrative modeling and validation to apply conformational restraints and assess agreement with the most common structural techniques for IDPs: Nuclear Magnetic Resonance (NMR) spectroscopy, Small-angle X-ray Scattering (SAXS), and single-molecule Förster Resonance Energy Transfer (smFRET). Agreement with such a diverse set of experimental data suggests that details of the generated ensembles can now be examined with a high degree of confidence. Using the disordered N-terminal region of the Sic1 protein as a test case, we examined relationships between average global polymeric descriptions and higher-moments of their distributions. To resolve apparent discrepancies between smFRET and SAXS inferences, we integrated SAXS data with non-smFRET (NMR) data and reserved the smFRET data as an independent validation. Consistency with smFRET, which was not guaranteed a priori, indicates that, globally, the perturbative effects of NMR or smFRET labels on the Sic1 ensemble are minimal. Analysis of the ensembles revealed distinguishing features of Sic1, such as overall compactness and large end-to-end distance fluctuations, which are consistent with biophysical models of Sic1’s ultrasensitive binding to its partner Cdc4. Our results underscore the importance of integrative modeling and validation in generating and drawing conclusions from IDP conformational ensembles.

Many cell signaling events are coordinated by intrinsically disordered protein regions (IDRs) that undergo multisite Serine/Threonine phosphorylation. The conformational properties of these IDRs prior to and following multi-site phosphorylation are directly relevant to understanding their functions. Here, we present results from biophysical studies and molecular simulations that quantify the conformational properties of an 81-residue IDR from the S. cerevisiae transcription factor Ash1. We show that the unphosphorylated Ash1 IDR adopts coil-like conformations that are expanded and well-solvated. This result contradicts inferences regarding global compaction that are derived from heuristics based on amino acid compositions for IDRs with low proline contents. Upon phosphorylation at ten distinct sites, the global conformational properties of pAsh1 are indistinguishable from those of unphosphorylated Ash1. This insensitivity derives from compensatory changes to the pattern of local and long-range intra-chain contacts. We show that the conformational properties of Ash1 and pAsh1 can be explained in terms of the linear sequence patterning of proline and charged residues vis-à-vis all other residues. The sequence features of the Ash1 IDR are shared by many other IDRs that undergo multisite phosphorylation. Accordingly, we propose that our findings might be generalizable to other IDRs involved in cell signaling.

In the photosynthetic reaction center from Rhodobacter sphaeroides, the primary (Q(A)) and secondary (Q(B)) electron acceptors are both ubiquinone-10, but with very different properties and functions. To investigate the protein environment that imparts these functional differences, we have applied X-band HYSCORE, a 2D pulsed EPR technique, to characterize the exchangeable protons around the semiquinone (SQ) in the Q(A) and Q(B) sites, using samples of (15)N-labeled reaction centers, with the native high spin Fe(2+) exchanged for diamagnetic Zn(2+), prepared in (1)H(2)O and (2)H(2)O solvent. The powder HYSCORE method is first validated against the orientation-selected Q-band ENDOR study of the Q(A) SQ by Flores et al. (Biophys. J.2007, 92, 671-682), with good agreement for two exchangeable protons with anisotropic hyperfine tensor components, T, both in the range 4.6-5.4 MHz. HYSCORE was then applied to the Q(B) SQ where we found proton lines corresponding to T ≈ 5.2, 3.7 MHz and T ≈ 1.9 MHz. Density functional-based quantum mechanics/molecular mechanics (QM/MM) calculations, employing a model of the Q(B) site, were used to assign the observed couplings to specific hydrogen bonding interactions with the Q(B) SQ. These calculations allow us to assign the T = 5.2 MHz proton to the His-L190 N(δ)H···O(4) (carbonyl) hydrogen bonding interaction. The T = 3.7 MHz spectral feature most likely results from hydrogen bonding interactions of O1 (carbonyl) with both Gly-L225 peptide NH and Ser-L223 hydroxyl OH, which possess calculated couplings very close to this value. The smaller 1.9 MHz coupling is assigned to a weakly bound peptide NH proton of Ile-L224. The calculations performed with this structural model of the Q(B) site show less asymmetric distribution of unpaired spin density over the SQ than seen for the Q(A) site, consistent with available experimental data for (13)C and (17)O carbonyl hyperfine couplings. The implications of these interactions for Q(B) function and comparisons with the Q(A) site are discussed.