Predicting the sequence-dependent backbone dynamics of intrinsically disordered proteins

eLife assessment

In this useful study, a solid machine learning approach based on a broad set of systems to predict the R2 relaxation rates of residues in intrinsically disordered proteins (IDPs) is described. The ability to predict the patterns of R2 will be helpful to guide experimental studies of IDPs. A potential weakness is that the predicted R2 values may include both fast and slow motions, thus the predictions provide only limited new physical insights into the nature of the underlying protein dynamics, such as the most relevant timescale.

https://doi.org/10.7554/eLife.88958.3.sa0

Significance of the findings:

Useful: Findings that have focused importance and scope

Landmark
Fundamental
Important
Valuable
Useful

Strength of evidence:

Solid: Methods, data and analyses broadly support the claims with only minor weaknesses

Exceptional
Compelling
Convincing
Solid
Incomplete
Inadequate

During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments

Abstract
Introduction
Results
Discussion
Methods
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

How the sequences of intrinsically disordered proteins (IDPs) code for functions is still an enigma. Dynamics, in particular residue-specific dynamics, holds crucial clues. Enormous efforts have been spent to characterize residue-specific dynamics of IDPs, mainly through NMR spin relaxation experiments. Here, we present a sequence-based method, SeqDYN, for predicting residue-specific backbone dynamics of IDPs. SeqDYN employs a mathematical model with 21 parameters: one is a correlation length and 20 are the contributions of the amino acids to slow dynamics. Training on a set of 45 IDPs reveals aromatic, Arg, and long-branched aliphatic amino acids as the most active in slow dynamics whereas Gly and short polar amino acids as the least active. SeqDYN predictions not only provide an accurate and insightful characterization of sequence-dependent IDP dynamics but may also serve as indicators in a host of biophysical processes, including the propensities of IDP sequences to undergo phase separation.

Introduction

Intrinsically disordered proteins (IDPs) or regions (IDRs) do not have the luxury of a three-dimensional structure to help decipher the relationship between sequence and function. Instead, dynamics has emerged as a crucial link between sequence and function for IDPs (Dey et al., 2022). Nuclear magnetic resonance (NMR) spin relaxation is a uniquely powerful technique for characterizing IDP dynamics, capable of yielding residue-specific information (Camacho-Zarco et al., 2022). Backbone ¹⁵N relaxation experiments typically yield three parameters per residue: transverse relaxation rate ( $R_{2}$ ), longitudinal relaxation rate ( $R_{1}$ ), and steady-state heteronuclear Overhauser enhancement (NOE). While all three parameters depend on ps-ns dynamics, $R_{2}$ is the one most affected by slower dynamics (10 s of ns to 1 μs). An increase in either the timescale or the amplitude of slower dynamics results in higher $R_{2}$ . For IDPs, $R_{2}$ is also the parameter that exhibits the strongest dependence on sequence (Dey et al., 2022; Camacho-Zarco et al., 2022).

$R_{2}$ was noted early on as an important indicator of residual structure in the unfolded state of the structured protein lysozyme (Klein-Seetharaman et al., 2002). This property has since been measured for many IDPs to provide insight into various biophysical processes. Just as the residual structure in the unfolded state biases the folding pathway of lysozyme (Klein-Seetharaman et al., 2002), a nascent α-helix in the free state of Sendai virus nucleoprotein C-terminal domain (Sev-NT), as indicated by highly elevated $R_{2}$ (Abyzov et al., 2016), biases the coupled binding and folding pathway in the presence of its target phosphoprotein (Schneider et al., 2015). Local secondary structure preformation also facilitates the binding of yes-associated protein (YAP) with its target transcription factor (Feichtinger et al., 2022). Likewise a correlation has been found between $R_{2}$ in the free state and the membrane binding propensity of synaptobrevin-2: residues with elevated $R_{2}$ have increased propensity for membrane binding (Lakomek et al., 2019). $R_{2}$ in the free state has also been used to uncover factors that promote liquid-liquid phase separation of IDPs. For example, a nascent α-helix (shown by elevated $R_{2}$ ) is important for the phase separation of the TDP-43 low-complexity domain, as both the deletion of the helical region and a helix-breaking mutation (Ala to Pro) abrogates phase separation (Conicella et al., 2016). Similarly, nascent α-helices in the free state of cytosolic abundant heat-soluble 8 (CAHS-8), upon raising concentration and lowering temperature stabilize to form the core of fibrous gels (Malki et al., 2022). For the hnRNPA1 low-complexity domain (A1-LCD), aromatic residues giving rise to local peaks in $R_{2}$ also mediate phase separation (Martin et al., 2020).

Both NMR relaxation data and molecular dynamics (MD) simulations have revealed determinants of $R_{2}$ for IDPs. It has been noted that the flexible Gly tends to lower $R_{2}$ , whereas secondary structure and contact formation tend to raise $R_{2}$ (Cook et al., 2019). This conclusion agrees well with recent MD simulations (Dey et al., 2022; Hicks et al., 2020; Yu and Brüschweiler, 2022; Smrt et al., 2023). These MD studies, using IDP-specific force fields, are able to predict $R_{2}$ in quantitative agreement with NMR measurements, without ad hoc reweighting as done in earlier studies. According to MD, most contact clusters are formed by local sequences, within blocks of up to a dozen or so residues (Dey et al., 2022; Hicks et al., 2020; Smrt et al., 2023). Tertiary contacts can also form but are relatively rare; as such their accurate capture requires extremely extensive sampling and still poses a challenge for MD simulations. Contrary to Gly, aromatic residues have been noted as mediators of contact clusters (Klein-Seetharaman et al., 2002; Martin et al., 2020).

Schwalbe et al., 1997 introduced a mathematical model to describe the $R_{2}$ profile along the sequence for lysozyme in the unfolded state. The $R_{2}$ value of a given residue was expressed as the sum of contributions from this residue and its neighbors. This model yields a mostly flat profile across the sequence, except for a falloff at the termini, resulting in an overall bell shape. Klein-Seetharaman et al., 2002 then fit peaks above this flat profile as a sum of Gaussians. Cho et al., 2007 proposed bulkiness as a qualitative indicator of backbone dynamics. Recently Sekiyama et al., 2022 calculated $R_{2}$ as the geometric mean of ‘indices of local dynamics’; the latter were parameterized by fitting to the measured $R_{2}$ for a single IDP. All these models merely describe the $R_{2}$ profile of a given IDP, and none of them is predictive.

Here, we present a method, SeqDYN, for predicting $R_{2}$ of IDPs. Using a mathematical model introduced by Li et al., 2020 to predict propensities for binding nanoparticles and also adapted for predicting propensities for binding membranes (Qin et al., 2022), we express the $R_{2}$ value of a residue as the product of contributing factors from all residues. The contributing factor attenuates as the neighboring residue becomes more distant from the central residue. The model, after training on a set of 45 IDPs, has prediction accuracy that is competitive against that of the recent MD simulations using IDP-specific force fields (Dey et al., 2022; Hicks et al., 2020; Yu and Brüschweiler, 2022; Smrt et al., 2023). For lysozyme and other structured proteins, the SeqDYN prediction agrees remarkably well with $R_{2}$ measured in their unfolded state.

Results

The data set of IDPs with R₂ rates

We collected R₂ data for a total of 54 nonhomologous IDPs or IDRs (Table 1; Figure 1). According to indicators from NMR properties, including low or negative NOEs, narrow dispersion in backbone amide proton chemical shifts, and small secondary chemical shifts (SCSs), most of the proteins are disordered with at most transient α-helices. A few are partially folded, including Sev-NT with a well-populated (~80%) long helix (residues 478–491; Jensen et al., 2008), CREB-binding protein fourth intrinsically disordered linker (CBP-ID4) with >50% propensities for two long helices (residues 2–25 and 101–128; Piai et al., 2016), HOX transcription factor DFD (HOX-DFD) with a well-folded domain comprising three helices (Maiti et al., 2019), and Hahellin (apo form) as a molten globule (Patel et al., 2014). In Figure 2, we display representative conformations of five IDPs, ranging from fully disordered MAPK kinase 4 (MKK4; Delaforge et al., 2018) and α-synuclein (Sung and Eliezer, 2007) to Measles virus phosphoprotein N-terminal domain (Mev-P_NTD; Milles et al., 2018) with transient short helices to Sev-NT and CBP-ID4 with stable long helices. The sequences of all the IDPs are listed in Appendix 1.

Figure 1

Download asset Open asset

Clock-like tree plot showing lack of homology among the 45 IDPs.

The level of homology between two sequences is measured by the distance from their convergence point to the center of the clock. The highest level of apparent identity is between A1-LCD and TDP-43, at 25%, but these two proteins differ in both secondary structure formation and $R_{2}$ characteristics. There is, however, a 20-residue overlap between the N-terminus of MBP-xα2 and the C-terminus of rmBG21.

Figure 2

Download asset Open asset

Representative conformations of five IDPs.

(**A–E**) MKK4, α-synuclein, Mev-P_NTD, Sev-NT, and CBP-ID4. Conformations were initially generated using TraDES (http://trades.blueprint.org; Feldman and Hogue, 2002), selected to have radius of gyration close to predicted by a scaling function $R_{g} = 2.54 N^{0.522}$ (Å) (Bernadó and Blackledge, 2009). Conformations for residues predicted as helical by PsiPred plus filtering were replaced by an ideal helix. Finally residues are colored according to a scheme ranging from green for low predicted $R_{2}$ to red for high predicted $R_{2}$ .

Table 1

Experimental conditions, mean and standard deviation of measured $R_{2}$ , and SeqDYN prediction RMSE.

Protein name	# of res	Temp(K)	B₀ (MHz)	${\bar{R}}_{2}$ (s^–1)	$σ_{R_{2}}$ (s^–1)	RMSE(s^–1)	PMID; ref
Training set (45 IDPs) ^*
A1-LCD	131	298	800	2.68	0.46	0.60	32029630; Martin et al., 2020
Aβ40	40	278	600	3.40	0.92	0.38	31181936; Rezaei-Ghaleh et al., 2019
Ash1	83	278	800	9.80	1.40	1.41	27807972; Martin et al., 2016
Beclin1	165	288	800	5.37	1.03	1.14	27288992; Yao et al., 2016
CAPRIN1	103	303	600	5.34	0.88	0.72	31898464; Wong et al., 2020
CBP-ID4	207	283	700	5.45	2.55	2.01;1.90^†	29790640; Murrali et al., 2018
GbnD4-DHD	91	280	700	6.81	1.55	1.28	29309054; Jenner et al., 2018
ERD14	185	288	600	3.96	0.87	0.54	21336827; Szalainé Ágoston et al., 2011
ExsE	88	298	600	3.18	0.88	0.76	22138394; Zheng et al., 2012
FCP1	85	298	500	2.94	0.54	0.43	26286791; Lawrence and Showalter, 2012
FUS	163	298	850	3.48	0.51	0.54	26455390; Burke et al., 2015
GAb1	82	298	500	3.99	0.88	0.89	34929201; Gruber et al., 2022
hACTR	69	304	600	3.26	0.47	0.49	18177052; Ebert et al., 2008
Hahellin	92	298	800	9.94	2.69	2.85	24671380; Patel et al., 2014
hCSD1	141	298	500	3.56	0.93	0.99	18537264; Kiss et al., 2008
HOX-DFD	90	298	600	6.98	3.15	1.99	30802457; Maiti et al., 2019
hZIP4-ICL2	100	283	800	9.54	2.37	1.58	30793391; Bafaro et al., 2019
Jaburetox	94	298	800	6.01	2.30	2.27	25605001; Lopes et al., 2015
KRS-NT	72	303	600	3.26	0.93	0.83	24983501; Cho et al., 2014
MBP-xα2	70	295	600	3.83	0.60	0.54	25343306; De Avila et al., 2014
MKK4	86	278	850	4.49	1.42	0.63	29276882; Delaforge et al., 2018
N-Cby	63	298		4.19	1.20	1.25	21182262; Mokhtarzada et al., 2011
Niv-P_NTD	406	288	700	5.41	1.82	1.66	33177626; Schiavina et al., 2020
NS5A-D2D3	268	278	800	8.62	3.85	2.14	26445449; Sólyom et al., 2015
NUPR1	93	298	600	2.98	0.82	0.76	31325636; Neira et al., 2019
OPN	220	310	800	2.59	0.82	0.54	31794728; Mateos et al., 2020
p53TAD	73	298	850	2.72	0.66	0.33	30240067; Xie et al., 2018
PDEγ	87	298		3.96	1.05	0.71	18230733; Song et al., 2008
PKIα	75	300	900	3.41	0.87	0.52	32338601; Olivieri et al., 2020
Mev-P_NTD	304	298	950	2.92	0.59	0.48	30140745; Milles et al., 2018
ProTα	113	283	800	3.40	0.56	0.43	29466338; Borgia et al., 2018
Pup	64	298	850	2.66	0.51	0.43	30240067; Xie et al., 2018
rmBG21	199	300	600	4.06	0.90	0.63	17676872; Ahmed et al., 2007
RPB1	201	277	850	6.48	1.74	1.33	28945358; Janke et al., 2018
securin	202	283	500	5.49	1.13	1.08	19053469; Csizmok et al., 2008
Sev-NT	124	298	600	3.20	1.42	0.76;0.38^†	27112095; Abyzov et al., 2016
Sic1	92	278	500	3.34	0.59	0.48	20399186; Mittag et al., 2010
SKIPN	71	298		5.64	1.05	1.46	20007319; Wang et al., 2010
SLBP-NT	113	298	600	3.96	1.40	1.61	15260482; Thapar et al., 2004
α-synuclein	140	298	600	2.96	0.53	0.44	30184304; Rezaei-Ghaleh et al., 2018
SOCS5-JIR	70	303	800	4.32	2.36	1.91	26173083; Chandrashekaran et al., 2015
tau K18	129	283	700	4.12	0.95	0.83	23740819; Barré and Eliezer, 2013
TC1	106	298	600	4.65	1.61	1.24	23189168; Cino et al., 2012
TDP-43	151	283	500	4.07	1.51	0.96	27545621; Conicella et al., 2016
γ-tubulin-CT	39	288	500	2.23	0.35	0.27	29127738; Harris et al., 2018
Test set (9 IDPs)
AMOTL1	207	283	800	8.45	2.55	2.04	35481651; Vogel et al., 2022
CAHS-8	233	303	850	4.43	3.25	2.36;1.92^†	34750927; Malki et al., 2022
ChiZ	64	298	800	4.33	0.89	0.74	32585849; Hicks et al., 2020
α-endosulfine	121	298	800	3.21	0.81	0.48	34346186; Thapa et al., 2022
FtsQ	99	305	800	6.44	3.78	2.32;1.71^†	36959324; Smrt et al., 2023
Pdx1	83	298	500	2.98	0.70	0.76	30525611; Cook et al., 2019
synaptobrevin-2	96	278	600	5.54	1.80	0.72	30975750; Lakomek et al., 2019
TIA-1	91	310	800	4.01	0.89	0.55	36112647; Sekiyama et al., 2022
YAP	122	298	800	3.19	1.44	1.23	35378854; Feichtinger et al., 2022

*

For training set, RMSE is calculated for prediction based on leave-one-out training (using 44 IDPs).
†

First number is for SeqDYN prediction; second number is after applying a helix boost.

We used 45 of the 54 IDPs to train and validate SeqDYN and reserved the remaining 9 for testing. The sequence lengths of the training set range from 39 to 406 residues, with an average of 125.3 residues. Altogether $R_{2}$ data are available for 3966 residues. A large majority (35 out of 45) of the 45 IDPs have mean $R_{2}$ values ( ${\bar{R}}_{2}$ , calculated among all the residues in a protein) between 2.5 and 5.5 s^–1 (Table 1 and Figure 3A). This ${\bar{R}}_{2}$ range is much lower than that of structured proteins with similar sequence lengths. The low ${\bar{R}}_{2}$ values and lack of dependence on sequence length (Figure 3—figure supplement 1A) suggest that $R_{2}$ of the IDPs is mostly dictated by local sequence instead of tertiary interaction.

Figure 3 with 1 supplement see all

Download asset Open asset

Properties of the 45 IDPs in the training set.

(A) Histograms of means and standard deviations, calculated for individual proteins. Curves are drawn to guide the eye. Inset: correlation between ${\bar{R}}_{2}$ and $σ_{R_{2}}$ . (B) Experimental mean scaled $R_{2}$ ( $m s R_{2}$ ) and SeqDYN $q$ parameters, for the 20 types of amino acids. Note that Pro residues have low $m s R_{2}$ for the lack of backbone amide proton. Amino acids are in descending order of $q$ .

Figure 3—source data 1 Source data for Figure 3.: https://cdn.elifesciences.org/articles/88958/elife-88958-fig3-data1-v1.xlsx
Download elife-88958-fig3-data1-v1.xlsx

The most often used temperature for acquiring the $R_{2}$ data was 298 K, but low temperatures (277–280 K) were used in a few cases (Table 1 and Figure 3—figure supplement 1B). Of the seven IDPs with ${\bar{R}}_{2}$ > 6.4 s^–1, four can be attributed to low temperatures (Sólyom et al., 2015; Martin et al., 2016; Janke et al., 2018; Jenner et al., 2018), one is due to a relatively low temperature (283 K) as well as the presence of glycerol (20% v/v; Bafaro et al., 2019), and two can be explained by tertiary structure formation [a folded domain (Maiti et al., 2019) or molten globule (Patel et al., 2014)]. A simple reason for higher $R_{2}$ values at lower temperatures is the higher water viscosity, resulting in a slowdown in molecular tumbling; a similar effect is achieved by adding glycerol. In some cases, $R_{2}$ was measured at both low and room temperatures (Abyzov et al., 2016; Martin et al., 2020). To a good approximation, the effect of lowering temperature is a uniform scaling of $R_{2}$ across the IDP sequence. For Sev-NT, downscaling of the $R_{2}$ values at 278 K by a factor of 2.0 brings them into close agreement with those at 298 K (Figure 3—figure supplement 1C), with a root-mean-square-deviation (RMSD) of 0.5 s^–1 among all the residues. Likewise, for A1-LCD, downscaling by a factor of 2.4 brings the $R_{2}$ values at 288 K into good match with those at 298 K (Figure 3—figure supplement 1D), with an RMSD of 0.4 s^–1. Because SeqDYN is concerned with the sequence dependence of $R_{2}$ , a uniform scaling has no effect on model parameter or prediction; therefore mixing the data from different temperatures is justified. The same can be said about the different magnetic fields in acquiring the $R_{2}$ data (Table 1 and Figure 3—figure supplement 1E). Increasing the magnetic field raises $R_{2}$ values, and the effect is also approximated well by a uniform scaling (Abyzov et al., 2016; Conicella et al., 2016; Janke et al., 2018).

One measure on the level of sequence dependence of $R_{2}$ is the standard deviation, $σ_{R_{2}}$ , calculated among the residues of an IDP. Among the training set, the $R_{2}$ values of 30 IDPs have moderate sequence variations, with $σ_{R_{2}}$ ranging from 0.5 to 1.5 s^–1 (Table 1); the histogram of $σ_{R_{2}}$ calculated for the entire training set peaks around 0.75 s^–1 (Figure 3A). There is a moderate correlation between $σ_{R_{2}}$ and ${\bar{R}}_{2}$ (Figure 3A, inset), reflecting in part the fact that $σ_{R_{2}}$ can be raised simply by a uniform upscaling, for example as a result of lowering temperature. Still, only two of the five IDPs with high ${\bar{R}}_{2}$ attributable to lower temperature or presence of glycerol are among the seven IDPs with high sequence variations ( $σ_{R_{2}}$ >2 s^–1). Therefore, the sequence variation of $R_{2}$ as captured by $σ_{R_{2}}$ manifests mostly the intrinsic effect of the IDP sequence, not the influence of external factors such as temperature or magnetic field strength. The mean $σ_{R_{2}}$ value among the training set is 1.24 s^–1.

One way to eliminate the influence of external factors is to scale the $R_{2}$ values of each IDP by its ${\bar{R}}_{2}$ ; we refer to the results as scaled $R_{2}$ , or $s R_{2}$ . We then pooled the $s R_{2}$ values for all residues in the training set, and separated them according to amino-acid types. The amino acid type-specific mean $s R_{2}$ values, or $m s R_{2}$ , are displayed in Figure 3B. The seven amino acids with the highest $m s R_{2}$ in descending order are Trp, Arg, Tyr, Phe, Ile, His, and Leu. The presence of all the four aromatic amino acids in this “high-end” group immediately suggests π-π stacking as important for raising $m s R_{2}$ ; the presence of Arg further implicates cation-π interactions. In the other extreme, the seven amino acids with the lowest $m s R_{2}$ in ascending order are Gly, Cys, Val, Asp, Ser, Thr, and Asn. Gly is well-known as a flexible residue; it is also interesting that all the four amino acids with short polar sidechains are found in this “low-end” group. Pro has an excessively low $m s R_{2}$ [with data from only two IDPs (Murrali et al., 2018; Wong et al., 2020)], but that is due to the absence of an amide proton.

The SeqDYN model and parameters

The null model is to assume a uniform $R_{2}$ for all the residues in an IDP. The root-mean-square-error (RMSE) of the null model is equal to the standard deviation, $σ_{R_{2}}$ , of the measured $R_{2}$ values. The mean RMSE, $\bar{R M S E}$ , of the null model, equal to 1.24 s^–1 for the training set, serves as the upper bound for evaluating the errors of $R_{2}$ predictors. The next improvement is a one-residue predictor, where first each residue (with index n) assumes its amino acid-specific mean $s R_{2}$ ( $m s R_{2}$ ) and then a uniform scaling factor $Υ$ is applied:

R_{2} (n) = Υ ∙ {m s R}_{2} (n)

This one-residue model does only minutely better than the null model, with a $\bar{R M S E}$ of 1.22 s^–1.

In SeqDYN, we account for the influence of neighboring residues. Specifically, each residue i contributes a factor $f (i; n)$ to the $R_{2}$ value of residue n. Therefore,

R_{2} (n) = Υ \prod_{i = 1}^{N} f (i; n)

where N is the total number of residues in the IDP. The contributing factor depends on the sequence distance $s = | i - n |$ and the amino-acid type of residue $i$ :

f (i; n) = 1 + \frac{q (i) - 1}{1 + b s^{2}}

There are 21 global parameters. The first 20 are the $q$ values, one for each of the 20 types of amino acids; the last parameter is $b$ , appearing in the Lorentzian form of the sequence-distance dependence. We define the correlation length, $L_{c o r r}$ , as the sequence distance at which the contributing factor is midway between the values at $s = 0$ and $\infty$ . It is easy to verify that $L_{c o r r} = b^{- 1 / 2}$ . Note that the single-residue model can be seen as a special case of SeqDYN, with $L_{c o r r}$ set to 0 and $q$ set to $m s R_{2}$ .

The functional forms of Equation 2a and Equation 2b were adapted from Li et al., 2020; we also used them for predicting residue-specific membrane association propensities of IDPs (Qin et al., 2022). In these previous applications, a linear term was also present in the denominator of Equation 2b. In our initial training of SeqDYN, the coefficient of the linear term always converged to near zero. We thus eliminated the linear term. In addition to the Lorentzian form, we also tested a Gaussian form for the sequence-distance dependence and found somewhat worse performance. The more gradual attenuation of the Lorentzian form with increasing sequence distance evidently provides an overall better model for the $R_{2}$ data in the entire training set. Others (Cho et al., 2007; Sekiyama et al., 2022; Delaforge et al., 2018) have modeled $R_{2}$ as the average of some parameters over a window; a window has an extremely abrupt sequence-distance dependence (1 for $s < L_{c o r r}$ and 0 for $s > L_{c o r r}$ ).

We parametrized the SeqDYN model represented by Equation 2a and Equation 2b on the training set of 45 IDPs. In addition to the 21 global parameters noted above, there are also 45 local parameters, namely one uniform scaling factor ( $Υ$ ) per IDP. The parameter values were selected to minimize the sum of the mean-square-errors for the IDPs in the training set, calculated on $R_{2}$ data for a total of 3924 residues. We excluded the 42 Pro residues in the training set because, as already noted, their $R_{2}$ values are lower for chemical reasons. We will present validation and test results below, but first let us look at the parameter values.

The $q$ values are displayed in Figure 3B alongside $m s R_{2}$ . In descending order, the seven amino acids with the highest $q$ values are Trp, Ile, Tyr, Arg, His, Phe, and Leu. These are exactly the same amino acids in the high-end group for $m s R_{2}$ , though their order there is somewhat different. In ascending order, the seven amino acids (excluding Pro) with the lowest $q$ values are Gly, Asn, Ser, Asp, Val, Thr, and Cys. The composition of the low-end group is also identical to that for $m s R_{2}$ . The $q$ values thus also suggest that π-π and cation-π interactions in local sequences may raise $R_{2}$ , whereas Gly and short-polar residues may lower $R_{2}$ .

Given the common amino acids at both the high and low ends for $m s R_{2}$ and $q$ , it is not surprising that these two properties exhibit a strong correlation, with a coefficient of determination (R²; excluding Pro) at 0.92 (Figure 4A). Also, because the high-end group contains the largest amino acids (e.g. Trp and Tyr) whereas the low-end group contains the smallest amino acids (e.g. Gly and Ser), we anticipated some correlation of $m s R_{2}$ and $q$ with amino-acid size. We measure the latter property by the molecular mass (m). As shown in Figure 4B, both $m s R_{2}$ and $q$ indeed show a medium correlation with m, with R²=0.67 (excluding Pro) and 0.61, respectively. A bulkiness parameter was proposed as an indicator of sequence-dependent backbone dynamics of IDPs (Cho et al., 2007; Delaforge et al., 2018). Bulkiness was defined as the sidechain volume-to-length ratio, and identified amino acids with aromatic or branched aliphatic sidechains as bulky (Zimmerman et al., 1968). We found only modest correlations between either $m s R_{2}$ or $q$ and bulkiness, with R² just below 0.4 (Figure 4C).

Figure 4 with 1 supplement see all

Download asset Open asset

SeqDYN model parameters.

(A) Correlation between $m s R_{2}$ and $q$ . The values are also displayed as bars in Figure 3B. (B) Correlation of $m s R_{2}$ and $q$ with amino-acid molecular mass. (C) Correlation of $m s R_{2}$ and $q$ with bulkiness. (D) The optimal correlation length and deterioration of SeqDYN prediction as the correlation length is moved away from the optimal value.

Figure 4—source data 1 Source data for Figure 4.: https://cdn.elifesciences.org/articles/88958/elife-88958-fig4-data1-v1.xlsx
Download elife-88958-fig4-data1-v1.xlsx

The optimized value of $b$ is 3.164 × 10^–2, corresponding to an $L_{c o r r}$ of 5.6 residues. The resulting optimized $\bar{R M S E}$ is 0.95 s^–1, a clear improvement over the value 1.24 s^–1 of the null model. To check the sensitivity of prediction accuracy to $b$ , we set $b$ to values corresponding to $L_{c o r r}$ = 0, 1, 2,…, and retrained SeqDYN for $b$ fixed at each value (Figure 4D). Note that the null-model $\bar{R M S E}$ , 1.24 s^–1, sets an upper bound. This upper bound is slowly reached when $L_{c o r r}$ is increased from the optimal value. In the opposite direction, when $L_{c o r r}$ is decreased from the optimal value, $\bar{R M S E}$ rises quickly, reaching 1.22 s^–1 at $L_{c o r r}$ = 0. The latter $\bar{R M S E}$ is the same as that of the single-residue model. Lastly we note that there is a strong correlation between the uniform scaling factors and ${\bar{R}}_{2}$ values among the 45 IDPs (R²=0.77), as to be expected. For 39 of the 45 IDPs, $Υ$ values fall in the range of 0.8–2.0 s^–1.

As presented next, we evaluate the performance of SeqDYN by leave-one-out cross validation, where each IDP in turn was left out of the training set and the model was trained on the remaining 44 IDPs to predict $R_{2}$ for the IDP that was left out. The parameters from the leave-one-out (also known as jackknife) training sessions allow us to assess the potential bias of the training set. For this purpose, we compare the values of the 21 global parameters, either from the full training set or from taking the averages of the jackknife training sessions. For each of the $q$ parameters, the values from these two methods differ only in the fourth digit; for example for Leu, they are both 1.1447 from full training and from jackknife training. The values for $b$ are 3.164×10^–2 from full training as stated above and 3.163×10^–2 from jackknife training. The close agreement in parameter values between full training and jackknife training suggests no significant bias in the training set.

Another question of interest is whether the difference between the $q$ parameters of two amino acids is statistically significant. To answer this question, we carried out fivefold cross-validation training, resulting in five independent estimates for each parameter. For example, the mean ±standard deviation of the $q$ parameter is 1.1405 ± 0.0066 for Leu and 1.2174 ± 0.0211 for Ile. A t-test shows that their difference is extremely statistically significant (P<0.0001). In contrast, the difference between Leu and Phe ( $q$ =1.1552 ± 0.0304) is not significant. t-test results for other pairs of amino acids are found in Figure 4—figure supplement 1.

Validation of SeqDYN predictions

We now present leave-one-out cross-validation results. We denote the RMSE of the $R_{2}$ prediction for the left-out IDP as RMSE(–1). As expected, RMSE(–1) is higher than the RMSE obtained with the IDP kept in the training set, but the increases are generally slight. Specifically, all but eight of the IDPs have increases <0.1 s^–1; the largest increase is 0.35 s^–1, for CBP-ID4. The mean RMSE(–1), or $\bar{R M S E} (- 1)$ , for the 45 IDPs is increased by 0.05 s^–1 over $\bar{R M S E}$ , to 1.00 s^–1. The latter value is still a distinct improvement over the mean RMSE 1.24 s^–1 of the null model. The histogram of RMSE(–1) for the 45 IDPs is shown in Figure 5A. It peaks at 0.5 s^–1, which is a substantial downshift from the corresponding peak at 0.75 s^–1 for $σ_{R_{2}}$ (Figure 3A). Thirty-four of the 45 IDPs have RMSE(–1) values lower than the corresponding $σ_{R_{2}}$ .

Figure 5 with 1 supplement see all

Download asset Open asset

Quality of SeqDYN predictions.

(A) Histogram of RMSE(–1). Letters indicate RMSE(–1) values of the IDPs to be presented in panels (**B–F**). (**B–F**) Measured (bars) and predicted (curves) $R_{2}$ profiles for MKK4, α-synuclein, Mev-P_NTD, Sev-NT, and CBP-ID4. In (E) and (F), green curves are SeqDYN predictions and red curves are obtained after a helix boost.

Figure 5—source data 1 Source data for Figure 5.: https://cdn.elifesciences.org/articles/88958/elife-88958-fig5-data1-v1.xlsx
Download elife-88958-fig5-data1-v1.xlsx

To further illustrate the performance of SeqDYN, we present the comparison of predicted and measured $R_{2}$ values for five IDPs: MKK4, α-synuclein, Mev-P_NTD, Sev-NT, and CBP-ID4 (Figure 5B–F). A simple common feature is the falloff of $R_{2}$ at the N- and C-termini, resulting from missing upstream or downstream residues that otherwise would be coupled to the terminal residues, as first recognized by Schwalbe et al., 1997. Representative conformations of the five IDPs are displayed in Figure 2, with residues colored according to the predicted $R_{2}$ values. For four of these IDPs, the RMSE(–1) values range from 0.44 to 0.76 s^–1 and are scattered around the peak of the histogram, while the RMSE(–1) for the fifth IDP, namely CBP-ID4, the RMSE(–1) value is 2.01 s^–1 and falls on the tail of the histogram (Figure 5A). Figure 5B displays the measured and predicted $R_{2}$ for MKK4. SeqDYN correctly predicts higher $R_{2}$ values in the second half of the sequence than in the first half. It even correctly predicts the peak around residue Arg75. The sequence in this region is H₇₂IERLRTH₇₉; six of these eight residues belong to the high-end group. In contrast, the lowest $R_{2}$ values occur in the sequence S₇GGGGSGGGSGSG₁₉, comprising entirely of two amino acids in the low-end group.

$R_{2}$ values for α-synuclein are shown in Figure 5C. Here, SeqDYN correctly predicts higher $R_{2}$ near the C-terminus and a dip around Gly68. However, it misses the $R_{2}$ peaks around Tyr39 and Asp121. MD simulations Dey et al., 2022 have found that these $R_{2}$ peaks can be explained by a combination of secondary structure formation (β-sheet around Tyr39 and polyproline II helix around Asp121) and local (between Tyr39 and Ser42) or long-range (between Asp121 and Lys96) interactions. SeqDYN cannot account for long-range interactions (e.g. between β-strands and between Asp121 and Lys96). Figure 5D shows that SeqDYN gives excellent $R_{2}$ predictions for Mev-P_NTD. It correctly predicts the high peaks around Arg17, Glu31, Leu193, and lower peaks around Arg235 and Trp285, but does underpredict the narrow peak around Tyr113.

The overall $R_{2}$ profile of Sev-NT is predicted well by SeqDYN, but the peak in the long helical region (residues 478–491) is severely underestimated (green curve in Figure 5E). A similar situation occurs for CBP-ID4, where the peak in the second long helical region (around Glu113) is underpredicted (green curve in Figure 5F). While the measured $R_{2}$ exhibits a higher peak in the second helical region than in the first helical region (around Arg16), the opposite is predicted by SeqDYN. When the $R_{2}$ data were included in the training set (i.e., full training), the second peak is higher than the first one, but that is not a real prediction because the $R_{2}$ data themselves were used for training the model. It merely means that the SeqDYN functions can be parameterized to produce any prescribed $R_{2}$ profile along the sequence. Indeed, when the $R_{2}$ data of CBP-ID4 alone were used to parameterize SeqDYN, the measured $R_{2}$ profile is closely reproduced (Figure 5—figure supplement 1). The reversal in $R_{2}$ peak heights between the two helical regions is the reason for the aforementioned unusual increase in RMSE when CBP-ID4 was left out of the training set.

R2 boost in long helical regions

It is apparent that SeqDYN underestimates the $R_{2}$ of stable long helices. Transient short helices does not seem to be a problem, since these are present, for example in Mev-P_NTD, where transient helix formation in the first 37 residues and between residues 189–198 (Milles et al., 2018) coincides with $R_{2}$ peaks that are correctly predicted by SeqDYN. SeqDYN can treat coupling between residues within the correlation length of 5.6 residues, but a much longer helix would tumble more slowly than implied by an $L_{c o r r}$ of 5.6, and thus it makes sense that SeqDYN would underestimate $R_{2}$ in that case.

Our solution then is to apply a boost factor to the long helical region. To do so, we have to know whether an IDP does form long helices and if so what the constituent residues are. Secondary structure predictors tend to overpredict α-helices and β-strands for IDPs, as they are trained on structured proteins. One way to counter that tendency is to make the criteria for α-helices and β-strands stricter. We found that, by filtering PsiPred (http://bioinf.cs.ucl.ac.uk/psipred; McGuffin et al., 2000) helix propensity scores ( $p_{H l x}$ ) with a very high cutoff of 0.99, the surviving helix predictions usually correspond well with residues identified by NMR as having high helix propensities. For example, for Mev-P_NTD, PsiPred plus filtering predicts residues 14–17, 28–33, and 191–193 as helical; all of them are in regions that form transient helices according to chemical shifts (Milles et al., 2018). Likewise long helices are also correctly predicted for Sev-NT (residues 477–489) and CBP-ID4 (residues 6–17 and 105–116; Jensen et al., 2008; Piai et al., 2016).

We apply a boost factor, $B_{H l x}$ , to helices with a threshold length of 12:

B_{H l x} = 1 + α p_{H l x} Θ (p_{H l x} \geq 0.99; L_{H l x} \geq 12)

The $Θ$ function is 1 if the helix propensity score is above the filtering cutoff and the helix length ( $L_{H l x}$ ) is above the threshold, and 0 otherwise With a boost amplitude $α$ at 0.5, the boosted SeqDYN prediction for Sev-NT reaches excellent agreement with the measured $R_{2}$ (Figure 5E, red curve). The RMSE(–1) is reduced from 0.76 s^–1 to 0.38 s^–1 upon boosting. Applying the same helix boost to CBP-ID4 also results in a modest reduction in RMSE(–1), from 2.01 to 1.90 s^–1 (Figure 5F, red curve). The only other IDP for which PsiPred plus filtering predicts a long helix is the N-terminal region of lysyl-tRNA synthetase (KRS-NT). The authors who studied this protein did not report on secondary structure (Cho et al., 2014), but feeding their reported chemical shifts to the TALOS +server (https://spin.niddk.nih.gov/bax/nmrserver/talos/; Shen et al., 2009) found only short stretches of residues that fall into the helical region of the Ramachandran map. The SeqDYN prediction for KRS-NT is already good [RMSE(–1)=0.83 s^–1]; applying a helix boost would deteriorate the RMSE(−1) to 1.16 s^–1.

Further test on a set of nine IDPs

We have reserved nine IDPs for testing SeqDYN (parameterized on the training set of 45 IDPs). The level of disorder in these test proteins also spans the full range, from absence of secondary structures [ChiZ N-terminal region (Hicks et al., 2020), Pdx1 C-terminal region (Cook et al., 2019), and TIA-1 prion-like domain (Sekiyama et al., 2022)] to presence of transient short helices [synaptobrevin-2 (Lakomek et al., 2019), α-endosulfine (Thapa et al., 2022), YAP (Feichtinger et al., 2022), angiomotin-like 1 (AMOTL1) (Vogel et al., 2022)] to formation of stable long helices [FtsQ Smrt et al., 2023 and CAHS-8 Malki et al., 2022]. For eight of the nine test IDPs, the RMSEs of SeqDYN predictions are lower than the experimental $σ_{R_{2}}$ values, by an average of 0.66 s^–1. For the ninth IDP (Pdx1), the SeqDYN RMSE is slightly higher, by 0.06 s^–1, than the experimental $σ_{R_{2}}$ . Together, the nine test IDPs have a mean RMSE of 1.13 s^–1, close to the $\bar{R M S E} (- 1)$ of 1.00 s^–1 for the training set in the leave-one-out cross-validation.

The comparison of predicted and measured $R_{2}$ profiles along the sequence is presented in Figure 6A–I. For ChiZ, SeqDYN correctly predicts the major peak around Arg25 and the minor peak around Arg46 (Figure 6A). The $R_{2}$ profile of Pdx1 is largely featureless, except for a dip around Gly216, which is correctly predicted by SeqDYN (Figure 6B). Correct prediction is also obtained for the higher $R_{2}$ in the first half of TIA-1 prion-like domain than in the second half (Figure 6C). SeqDYN gives an excellent prediction for synaptobrevin-2, including a linear increase up to Arg56 and the major peak around Trp89 (Figure 6D).

Figure 6

Download asset Open asset

Measured (bars) and predicted (curves) $R_{2}$ profiles for ChiZ N-terminal region, TIA1 prion-like domain, Pdx1 C-terminal region, synaptobrevin-2, α-endosulfine, YAP, AMOTL1, FtsQ, and CAHS-8.

In (C), $R_{2}$ does not fall off at the N-terminus because the sequence is preceded by an expression tag MGSSHHHHHHHHHHHHS. In (H) and (I), green curves are SeqDYN predictions and red curves are obtained after a helix boost.

Figure 6—source data 1 Source data for Figure 6.: https://cdn.elifesciences.org/articles/88958/elife-88958-fig6-data1-v1.xlsx
Download elife-88958-fig6-data1-v1.xlsx

The prediction is also very good for α-endosulfine, including elevated $R_{2}$ around Glu34, which coincides with the presence of a transient helix, and depressed $R_{2}$ in the last 40 residues (Figure 6E). The only miss is an underprediction for the peak around Lys74. SeqDYN also predicts well the overall shape of the $R_{2}$ profile for YAP, including peaks around Asn70, Leu91, Arg124, and Arg161, but severely underestimates the peak height around Asn70 (Figure 6F). NOE signals indicate contacts between Met86, Leu91, Fhe95, and Fhe96 (Feichtinger et al., 2022); evidently this type of local contacts is captured well by SeqDYN. The $R_{2}$ elevation around Asn70 is mostly due to helix formation: residues 61–74 have helix propensities up to 40% (Feichtinger et al., 2022). PsiPred predicts helix for residues 62–73, but only residues 65–68 survive the filtering that we impose, resulting in a helix that is too short to apply a helix boost. The prediction for AMOTL1 is mostly satisfactory, including peaks around Phe200 and Arg264 and a significant dip around Gly292 (Figure 6G). However, whereas the two peaks have approximately equal heights in the measured $R_{2}$ profile, the predicted peak height around Phe200 is too low. SCSs indicate helix propensity around both $R_{2}$ peaks (Vogel et al., 2022). PsiPred also predicts helix in both regions, but only five and two residues, respectively, survive after filtering, and are too short for applying a helix boost.

For FtsQ, SeqDYN correctly predicts elevated $R_{2}$ for the long helix [residues 46–74 Smrt et al., 2023] but underestimates the magnitude (RMSE = 2.32 s^–1; green curve in Figure 6H). PsiPred plus filtering predicts a long helix formed by residues 47–73. Applying the helix boost substantially improves the agreement with the measured $R_{2}$ , with RMSE reducing to 1.71 s^–1 (red curve in Figure 6H). SeqDYN also gives a qualitatively correct $R_{2}$ profile for CAHS-8, with higher $R_{2}$ for the middle section (residues 95–190; RMSE = 2.36 s^–1; green curve in Figure 6I). However, it misses the extra elevation in $R_{2}$ for the first half of the middle section (residues 95–145). According to SCS, the first and second halves have helix propensities of 60% and 30%, respectively (Malki et al., 2022). PsiPred plus filtering predicts helices for residues 96–121, 124–141, 169–173, and 179–189. Only the first two helices, both in the first half of the middle section, are considered long according to our threshold. Once again, applying the helix boost leads to marked improvement in the predicted in $R_{2}$ , with RMSE reducing to 1.92 s^–1 (red curve in Figure 6I).

Inputting the sequences of structured proteins predicts $R_{2}$ in the unfolded state

SeqDYN is trained on IDPs, what if we feed it with the sequence of a structured protein? The prediction using the sequence of hen egg white lysozyme, a well-studied single-domain protein, is displayed in Figure 7A. It shows remarkable agreement with the $R_{2}$ profile measured by Klein-Klein-Seetharaman et al., 2002 in the unfolded state (denatured by 8 M urea at pH 2 and reduced to break disulfide bridges), including a major peak around Trp62, a second peak around Trp111, and a third peak around Trp123. Klein-Seetharaman et al. mutated Trp62 to Gly and the major peak all but disappeared. This result is also precisely predicted by SeqDYN with the mutant sequence (Figure 7B).

Figure 7

Download asset Open asset

$R_{2}$ profiles predicted (curves) by SeqDYN show close agreement with those measured (bars) on structured proteins in the unfolded state.

(A) Wild-type lysozyme (8 M urea; pH 2; cysteine-methylated). (B) Lysozyme with Trp62 to Gly mutation (pH 2). Methylated cysteines were treated as Ala in the SeqDYN predictions. (C) Apomyoglogin (8 M urea; pH 2.3). (D) Ubiquitin (8 M urea; pH 2).

Figure 7—source data 1 Source data for Figure 7.: https://cdn.elifesciences.org/articles/88958/elife-88958-fig7-data1-v1.xlsx
Download elife-88958-fig7-data1-v1.xlsx

SeqDYN also predicts well the $R_{2}$ profiles of other proteins in the unfolded state. For unfolded apomyoglobin (8 M urea; pH 2.3), Schwarzinger et al., 2002 claimed that depressed $R_{2}$ corresponded to stretches of small amino acids (Gly and Ala), whereas elevated corresponded to local hydrophobic interactions. SeqDYN reproduces all the observed peaks and valleys in the $R_{2}$ profile (Figure 7C). The deepest valley indeed occurs over a Gly/Ala-rich stretch, G₁₂₅ADAQGA₁₃₁, but the highest peak occurs over a stretch, I₁₀₂KYLEFI₁₀₈, that contains both hydrophobic and charged residues, all of which are on the high end of the $q$ parameters (Figure 3B). The $R_{2}$ profile of unfolded ubiquitin (8 M urea; pH 2) is relatively flat, which Wirmer et al., 2006 attributed to lack of residual secondary structure, based on the assumption that β-sheets (major elements of folded ubiquitin) are less resistant to denaturation than α-helices. SeqDYN predicts a relatively flat $R_{2}$ profile (Figure 7D), but the reason is that the ubiquitin sequence lacks a contiguous stretch of high- $q$ amino acids.

Discussion

We have developed a powerful method, SeqDYN, that predicts the backbone amide transverse relaxation rates ( $R_{2}$ ) of IDPs. The method is based on IDP sequences, is extremely fast, and available as a web server at https://zhougroup-uic.github.io/SeqDYNidp/ (Qin and Zhou, 2024). The excellent performance supports the notion that the ns-dynamics reported by $R_{2}$ is coded by the local sequence, comprising up to 6 residues on either side of a given residue. The amino-acid types that contribute the most to coupling within a local sequence are aromatic (Trp, Tyr, Phe, and His), Arg, and long branched aliphatic (Ile and Leu), suggesting the importance of π-π, cation-π, and hydrophobic interactions in raising $R_{2}$ . These interactions are interrupted by Gly and amino acids with short polar sidechains (Ser, Thr, Asn, and Asp), leading to reduced $R_{2}$ . Transient short helices produce moderate elevation in $R_{2}$ , whereas stable long helices result in a big boost in $R_{2}$ . Tertiary contacts can also raise $R_{2}$ , but appear to be infrequent in most IDPs (Dey et al., 2022).

It is also possible that $R_{2}$ reported by backbone amide ¹⁵N relaxation (as is the case for most of the IDPs studied here) may not be particularly sensitive to exchange effects, which likely involve tertiary contact formation. For the D2 domain of p27^Kip1, the exchange contributions measured using ¹⁵N relaxation were small (<2.5 s^–1) but were as large as 25 s^–1 when measured by high-power ¹H relaxation dispersion (Ban et al., 2017). This experiment measures the effective transverse relation rate, $R_{2, e f f}$ , over a range of effective radiofrequency $ω_{e f f}$ . The exchange contribution is maximal for the value $R_{2, e f f}^{l o w ω_{e f f}}$ at low $ω_{e f f}$ but is largely quenched for the value $R_{2,0}^{a p p}$ in the high- $ω_{e f f}$ limit. The SeqDYN prediction for this IDP matches much better with $R_{2,0}^{a p p}$ than with $R_{2, e f f}^{l o w ω_{e f f}}$ (Figure 8). It is not clear whether this IDP is unique in forming persistent tertiary contacts that give rise to substantial exchange contributions or the ¹H relaxation dispersion experiment is unique in reporting the exchange contributions. At the minimum, SeqDYN yields the exchange-free portion of the transverse relaxation rate, enabling easy identification of residues that potentially participate in tertiary contacts. For the D2 domain of p27^Kip1, SeqDYN correctly predicts the $R_{2,0}^{a p p}$ local maxima at W76 and Y88. It is these same two residues that show substantial exchange contributions and putatively participate in tertiary contact (Ban et al., 2017). Therefore local contacts may seed tertiary contacts. If $R_{2, e f f}$ data with substantial exchange contributions become available for more IDPs, SeqDYN may be retrained to make predictions for IDPs forming persistent tertiary contacts.

Figure 8

Download asset Open asset

Comparison between SeqDYN prediction (curves) and effective transverse relaxation rate (bars) from ¹H dispersion relaxation experiment.

(A) $R_{2, e f f}$ in the high- $ω_{e f f}$ limit. (B) $R_{2, e f f}$ at low $ω_{e f f}$ .

Figure 8—source data 1 Source data for Figure 8.: https://cdn.elifesciences.org/articles/88958/elife-88958-fig8-data1-v1.xlsx
Download elife-88958-fig8-data1-v1.xlsx

The $q$ parameters, while introduced here to characterize the propensities of amino acids to participate in local interactions, appear to correlate with the tendencies of amino acids to drive liquid-liquid phase separation. Consistent with the rank order of $q$ , Trp, Tyr, and Arg have been reported to be strong drivers of phase separation, Lys is a moderate driver, whereas Gly and Ser suppress phase separation (Martin et al., 2020; Wong et al., 2020; Wang et al., 2018). Recent measurements of the threshold concentration produced the following order for the propensity of phase separation by eight nonpolar amino acids in homotetrapeptides of the form XXssXX (ss: backbone disulfide bond): Trp > Phe > Leu>Met > Ile>Val > Ala>Pro (Zhang et al., 2024). This order is the same as that of the $q$ parameters, except that the $q$ values of Ile and Val are in the second and last places, respectively. Threshold concentrations of IDPs are now predicted reasonably well by coarse-grained simulations where each amino acid is modeled by a single bead with a Lennard-Jones diameter $d_{0}$ and a stickiness parameter $λ$ (Tesei and Lindorff-Larsen, 2022). Our $q$ parameter shows a good correlation (R²=0.59) with the compound parameter $d_{0}^{3} λ$ (Figure 9). Therefore, the $q$ parameter may serve as a predictor for the tendency of an amino acid to drive phase separation. In essence, the same ability of an amino acid, for example Trp, to form interactions with neighboring residues of an IDP in the free state also applies when it comes to interactions with residues on neighboring chains in a dense phase.

Figure 9

Download asset Open asset

Correlation between the stickiness parameters (λ) and the NMR relaxation parameters (q).

The regression line is shown as dashes.

Figure 9—source data 1 Source data for Figure 9.: https://cdn.elifesciences.org/articles/88958/elife-88958-fig9-data1-v1.xlsx
Download elife-88958-fig9-data1-v1.xlsx

Our method incorporates ideas from a number of previous efforts at describing $R_{2}$ . The first serious effort was by Schwalbe et al., 1997, who accounted for contributions from neighboring residues as additive terms, instead of multiplicative factors as in SeqDYN. Cho et al., 2007 and Delaforge et al., 2018 used the running average of the bulkiness parameter over a window of five to nine residues as a qualitative indicator of $R_{2}$ . Here again the calculation was based on an additive model. Sekiyama et al., 2022 employed a multiplicative model, with $R_{2}$ calculated as a geometric mean of ‘indices of local dynamics’ over a five-residue window. These indices, akin to our $q$ parameters, were trained on a single IDR (TIA-1 prion-like domain) and used to reproduce the measured $R_{2}$ for the same IDR. As we have illustrated on CBP-ID4 (Figure 5—figure supplement 1), training on a single protein merely biases the parameters to that model and has little value in predicting $R_{2}$ for other proteins. In comparison, SeqDYN is trained on 45 IDPs and its predictions are robust and achieve quantitative agreement with measured $R_{2}$ .

Ten of the IDPs tested here have been studied recently by MD simulations using IDP-specific force fields (Dey et al., 2022; Hicks et al., 2020; Yu and Brüschweiler, 2022; Smrt et al., 2023). In Table 2, we compare the RMSEs of SeqDYN predictions with those for $R_{2}$ calculations from MD simulations. For five of these IDPs: A1-LCD, Aβ40, α-synuclein, tau K18, and FtsQ, RMSEs of SeqDYN and MD are remarkably similar. Four of these IDPs lack significant population of α-helices or β-sheets, but FtsQ forms a stable long helix. For one other IDP, namely HOX-DFD, MD, by explicitly modeling its folded domain, does a much better job in predicting $R_{2}$ than SeqDYN (RMSEs of 1.40 s^–1 vs 1.99 s^–1). However, for the four remaining IDPs: p53TAD, Pup, Sev-NT, and ChiZ, SeqDYN significantly outperforms MD, with RMSEs averaging only 0.47 s^–1, compared to the MD counterpart of 1.14 s^–1. Overall, SeqDYN is very competitive against MD in predicting $R_{2}$ , but without the significant computational cost. While MD simulations can reveal details of local interactions, as noted for α-synuclein, and capture tertiary interactions if they occur, they still suffer from perennial problems of force-field imperfection and inadequate sampling. SeqDYN provides an accurate description of IDP dynamics at a ‘mean-field’ level, but could miss idiosyncratic behaviors of specific local sequences.

Table 2

RMSEs (s^–1) of $R_{2}$ predictions by SeqDYN and MD for 10 IDPs.

IDP name	SeqDYN	MD
A1-LCD	0.60*	0.59 ^§^{, ¶}
Aβ40	0.38*	0.38 ^§
HOX-DFD	1.99*	1.40 ^§
α-synuclein	0.44*	0.50 ^§
p53TAD	0.33*	1.04 ^**
Pup	0.43*	1.00**
Sev-NT	0.38*^,†	1.10 ^§^,††
tau K18	0.83*	0.80 ^§
ChiZ	0.74 ^‡	1.40 ^{‡ ‡}
FtsQ	1.71 ^§^,†	1.70 ^{§ §}

*

Based on leave-one-out training (using 44 IDPs).
†

Helix boost applied.
‡

Based on training by the full training set (45 IDPs).
§

From Dey et al., 2022.
¶

RMSE is scaled down by a factor of 2.39, to correct for the effect of temperature (MD at 288 K; see Figure 3—figure supplement 1C).
**

From Yu and Brüschweiler, 2022.
††

RMSE is scaled down by a factor of 2.99, to correct for the effects of temperature and magnetic field (MD at 274 K and 850 MHz; see Figure 3—figure supplement 1B).
‡ ‡

Originally calculated in Hicks et al., 2020 with correction in Hicks et al., 2021.
§ §

From Smrt et al., 2023.

Deep-learning models have become very powerful, but they usually have millions of parameters and require millions of protein sequences for training (Rives et al., 2021). In contrast, SeqDYN employs a mathematical model with dozens of parameters and requires only dozens of proteins for training. Reduced models (by collapsing amino acids into a small number of distinct types) have even been trained on <10 IDPs to predict propensities for binding nanoparticles (Li et al., 2020) or membranes (Qin et al., 2022). The mathematical model-based approach may be useful in other applications where data, similar to $R_{2}$ , are limited, including predictions of IDP secondary chemical shifts or residues that bind drug molecules (Robustelli et al., 2022) or protein targets, or even in protein design, for example for recognizing an antigenic site or a specific DNA site.

Methods

Collection of IDPs with measured $R_{2}$

Starting from six nonhomologous IDPs in our previous MD study (Dey et al., 2022), we obtained $R_{2}$ data for eight IDPs from the Bimolecular Magnetic Resonance Data Bank (BMRB; https://bmrb.io); data for two other IDPs were from our collaborators (Hicks et al., 2020; Smrt et al., 2023). Most of the 54 IDPs studied here were from searching the literature. Disorder was judged by dispersion in backbone amide proton chemical shifts, NOE, and SCS. $R_{2}$ data that were not available from the authors or BMRB were obtained by digitizing $R_{2}$ plots presented in figures of published papers, using WebPlotDigitizer (https://automeris.io/WebPlotDigitizer; Rohatgi, 2022) and further inspected visually.

Homology of IDPs was checked by sequence alignment using Clustal W (http://www.clustal.org/clustal2; Larkin et al., 2007), and presented as a clock-like tree using the ‘ape’ package (http://ape-package.ird.fr; Paradis et al., 2004). IDPs that had discernible homology with the selected training set were removed. Removed IDPs included HOX-SCR and β-synuclein from our previous MD study (Dey et al., 2022), due to homology with HOX-DFD and α-synuclein, respectively.

Coding for SeqDYN

The training of SeqDYN was coded in python, similar to our previous work for predicting residue-specific membrane association propensities (ReSMAP; https://zhougroup-uic.github.io/ReSMAPidp/; Qin et al., 2022). The cost function was the sum of mean-squared-errors for the IDPs in the training set. We used the least_squares function in scipy.optimize, with Trust Region Reflective as the minimization algorithm and all parameters restricted to the positive range. For the web server (https://zhougroup-uic.github.io/SeqDYNidp/; Qin and Zhou, 2024), we rewrote the prediction code javascript.

Appendix 1

Sequences of 54 IDPs

Terminal tags and other insertions are underlined.

Training set (45 IDPs)

>A1-LCD (Uniprot ID P04256, residues 186-320; deletion of 258-263)
GSMASASSSQRGRSGSGNFGGGRGGGFGGNDNFGRGGNFSGRGGFGGSRGGGGYGGSGDGYNGFGNDGSNFGGGGNYNNQSSNFGPMKGGNFGGRSSGPYGGGGQYFAKPRNQGGYGGSSSSSSYGSGRRF
>Aβ40 (Uniprot ID Q28053, residues 7-46)
DAEFRHDSGYEVHHQKLVFFAEDVGSNKGAIIGLMVGGV
>Ash1 (Uniprot ID P34233, residues 420-500)
GASASSSPSPSTPTKSGKMRSRSSSPVRPKAYTPSPRSPNYHRFALDSPPQSPRRSSNSSITKKGSRRSSGSSPTRHTTRVCV
>Beclin1 (Uniprot ID Q14457, residues 1-150)
MGSSHHHHHHSQDPMEGSKTSNNSTMQVSFVSQRSSQPLKLDTSFKILDRVTIQELTAPLLTTAQAKPGETQEEETNSGEEPFIETPRQDGVSRRFIPPARMMSTESANSFTLIGEASDGGTMENLSRRLKVTGDLFDIMSGQTDVDHPLSEESTDTLLDQLDTY
>CAPRIN1 (Uniprot ID Q14444, residues 607-709)
SRGVSRGGSRGARGLMNGYRGPANGFRGGYDGYRPSFSNTPNSGYTQSQFSAPRDYSGYQRDGYQQNFKRGSGQSGPRGAPRGRGGPPRPNRGMPQMNTQQVN
>CBP-ID4 (Uniprot ID Q92793, residues 1852-2057)
MQQQIQHRLQQAQLMRRRMATMNTRNVPQQSLPSPTSAPPGTPTQQPSTPQTPQPPAQPQPSPVSMSPAGFPSVARTQPPTTVSTGKPTSQVPAPPPPAQPPPAAVEAARQIEREAQQQQHLYRVNINNSMPPGRTGMGTPGSQMAPVSLNVPRPNQVSGPVMPSMPPGQWQQAPLPQQQPMPGLPRPVISMQAQAAVAGPRMPSVQ
>GbnD4-DHD (Uniprot ID A0A808VWJ6, residues 412-482)
MKHHHHHHHHGGLVPRGSHGSDEGVPDALRADTVPRAGPVRYARRRYWIGEARSDALAPAAPLEREPLPAEAMGAYFAIRRTDADDTVAAH
>ERD14 (Uniprot ID P42763, residues 1-185)
MAEEIKNVPEQEVPKVATEESSAEVTDRGLFDFLGKKKDETKPEETPIASEFEQKVHISEPEPEVKHESLLEKLHRSDSSSSSSSEEEGSDGEKRKKKKEKKKPTTEVEVKEEEKKGFMEKLKEKLPGHKKPEDGSAVAAAPVVVPPPVEEAHPVEKKGILEKIKEKLPGYHPKTTVEEEKKDKE
>ExsE (Uniprot ID Q9I322, residues 1-81)
MKIESIPPVQPSQDAGAEAVGHFEGRSVTRAAVRGDDRSSVAGLARWLARNVAGDPRSEQALQRLADGDGTPLEARTVRRREFLEGSS
>FCP1 (Uniprot ID Q9Y5B0, residues 879-961)
PGPEEQEEEPQPRKPGTRRERTLGAPASSERSAAGGRGPRGHKRKLNEEDAASESSRESSNEDEGSSSEADEMAKALEAELNDLM
>FUS (Uniprot ID P35637, residues 1-163)
MASNDYTQQATQSYGAYPTQPGQGYSQQSSQPYGQQSYSGYSQSTDTSGYGQSSYSSYGQSQNTGYGTQSTPQGYGSTGGYGSSQSSQSSYGQQSSYPGYGQQPAPSSTSGSYGSSSQSSSYGQPQSGSYSQQPSYGGQQQSYGQQQSYNPPQGYGQQNQYNS
>GAb1 (Uniprot ID B7Z3B9, residues 510-591)
SSPMIKPKGDKQVEYLDLDLDSGKSTPPRKQKSSGSGSSVADERVDYVVVDQQKTLALKSTREAWTDGRQSTESETPAKSVK
>hACTR (Uniprot ID Q9Y6Q9, residues 1023-1091)
GTQNRPLLRNSLDDLVGPPSNLEGQSDERALLDQLHTLLSNTDATGLEEIDRALGIPELVNQGQALEPK>Hahellin (Uniprot ID Q2SHN6, residues 162-252)
MGEKTVKLYEDTHFKGYSVELPVGDYNLSSLISRGALNDDLSSARVPSGLRLEVFQHNNFKGVRDFYTSDAAELSRDNDASSVRVSKMETTN
>hCSD1 (Uniprot ID P20810, residues 137-277)
AVPVESKPDKPSGKSGMDAALDDLIDTLGGPEETEEENTTYTGPEVSDPMSSTYIEELGKREVTIPPKYRELLAKKEGITGPPADSSKPIGPDDAIDALSSDFTCGSPTAAGKKTEKEESTEVLKAQSAGTVRSAAPPQEK
>HOX-DFD (Uniprot ID P07548, residues 337-426)
TDGERIIYPWMKKIHVAGVANGSYQPGMEPKRQRTAYTRHQILELEKEFHYNRYLTRRRRIEIAHTLVLSERQIKIWFQNRRMKWKKDNK
>hZIP4-ICL2 (Uniprot ID Q6P5W5, residues 424-498)
GDRGPEFELGTLPRDPEDLEDGPCGHSSHSHGGHSHGVSLQLAPSELRQPKPPHEGSRADLVAEESPELLNPEPRRLSPELRLLPYGHGLSAWSHPQFEK
>Jaburetox (Uniprot ID I1K3K3, residues 230-320)
MGPVNEANCKAAMEIVCRREFGHKEEEDASEGVTTGDPDCPFTKAIPREEYANKYGPTIGDKIRLGDTDLIAEIEKDFALYGDESVFGGGKVIH
>KRS-NT (Uniprot ID Q15046, residues 1-72)
MAAVQAAEVKVDGSEPKLSKNELKRRLKAEKKVAEKEAKQKELSEKQLSQATAAATNHTTDNGVGPEEESVD
>MBP-xα2 (Uniprot ID P04370, residues 172-237)
SIGRFFSGDRGAPKRGSGKDSHTRTTHYGSLPQKSQHGRTQDENPVVHFFKNIVTPRTPPPSQGKGRGLS
>MKK4 (Uniprot ID P45985, residues 1-86)
MAAPSPSGGGGSGGGSGSGTPGPVGSPAPGHPAVSSMQGKRKALKLNFANPPFKSTARFTLNPNPTGVQNPHIERLRTHSIESSGK
>N-Cby (Uniprot ID B0QY54, residues 1-63)
MPFFGNTFSPKKTPPRKSASLSNLHSLDRSTREVELGLEYGSPTMNLAGQSLKFENGQWIAET
>Niv-PNTD (Uniprot ID P0C1C7, residues 1-406)
MDKLELVNDGLNIIDFIQKNQKEIQKTYGRSSIQQPSIKDQTKAWEDFLQCTSGESEQVEGGMSKDDGDVERRNLEDLSSTSPTDGTIGKRVSNTRDWAEGSDDIQLDPVVTDVVYHDHGGECTGYGFTSSPERGWSDYTSGANNGNVCLVSDAKMLSYAPEIAVSKEDRETDLVHLENKLSTTGLNPTAVPFTLRNLSDPAKDSPVIAEHYYGLGVKEQNVGPQTSRNVNLDSIKLYTSDDEEADQLEFEDEFAGSSSEVIVGISPEDEEPSSVGGKPNESIGRTIEGQSIRDNLQAKDNKSTDVPGAGPKDSAVKEEPPQKRLPMLAEEFECSGSEDPIIRELLKENSLINCQQGKDAQPPYHWSIERSISPDKTEIVNGAVQTADRQRPGTPMPKSRGIPIKK
>NS5A-D2D3 (Uniprot ID O92972, residues 2163-2419)
GHMASGSLRGGEPEPDVTVLTSMLTDPSHITAETAKRRLARGSPPSLASSSASQLSAPSLKATCTTHHDSPDADLIEANLLWRQEMGGNITRVESENKVVILDSFEPLHADGDEREISVAAEILRKSRKFPSALPIWARPDYNPPLLESWKDPDYVPPVVHGCPLPPTKAPPIPPPRRKRTVVLTESNVSSALAELATKTFGSSGSSAVDSGTATALPDQASDDGDKGSDVESYSSMPPLEGEPGDPDLSDGSWSTVSEEASEDVVCC
>NUPR1 (Uniprot ID O60356, residues 2-82)
MRGSHHHHHHGSATFPPATSAPQQPPGPEDEDSSLDESDLYSLAHSYLGGGGRKGRTKREAAANTNRPSPGGHERKLVTKLQNSERKKRGARR
OPN (Uniprot ID F1NSM8, residues 46-264)
MHQDHVDSQSQEHLQQTQNDLASLQQTHYSSEENADVPEQPDFPDVPSKSQETVDDDDDDDNDSNDTDESDEVFTDFPTEAPVAPFNRGDNAGRGDSVAYGFRAKAHVVKASKIRKAARKLIEDDATTEDGDSQPAGLWWPKESREQNSRELPQHQSVENDSRPKFDSREVDGGDSKASAGVDSRESQGSVPAVDASNQTLESAEDAEDRHSIENNEVTR
>p53TAD (Uniprot ID P04637, residues 1-71)
MEEPQSDPSVEPPLSQETFSDLWKLLPENNVLSPLPSQAMDDLMLSPDDIEQWFTEDPGPDEAPRMPEAAPRV
>PDEγ (Uniprot ID P61248, residues 1-87)
MNLEPPKAEIRSATRVMGGPVTPRKGPPKFKQRQTRQFKSKPPKKGVQGFGDDIPGMEGLGTDITVIAPWEAFNHLELHELAQYGII
>PKIα (Uniprot ID P61925, residues 2-76)
TDVETTYADFIASGRTGRRNAIHDILVSSASGNSNELALKLAGLDINKTEGEEDAQRSSTEQSGEAQGEAAKSES
>Mev-PNTD (Uniprot ID P03422, residues 1-304)
MAEEQARHVKNGLECIRALKAEPIGSLAIEEAMAAWSEISDNPGQERATCREEKAGSSGLSKPCLSAIGSTEGGAPRIRGQGPGESDDDAETLGIPPRNLQASSTGLQCYYVYDHSGEAVKGIQDADSIMVQSGLDGDSTLSGGDNESENSDVDIGEPDTEGYAITDRGSAPISMGFRASDVETAEGGEIHELLRLQSRGNNFPKLGKTLNVPPPPDPGRASTSGTPIKKGTERRLASFGTEIASLLTGGATQCARKSPSEPSGPGAPAGNVPECVSNAALIQEWTPESGTTISPRSQNNEEGG
>ProTα (Uniprot ID P06454, residues 1-111)
GPMSDAAVDTSSEITTKDLKEKKEVVEEAENGRDAPANGNAENEENGEQEADNEVDEEEEEGGEEEEEEEEGDGEEEDGDEDEEAESATGKRAAEDDEDDDVDTKKQKTDEDD
>Pup (Uniprot ID P9WHN4, residues 1-64)
MAQEQTKRGGGGGDDDDIAGSTAAGQERREKLTEETDDLLDEIDDVLEENAEDFVRAYVQKGGQ
>rmBG21 (Uniprot ID P04370, residues 2-190)
GNHSGKRELSAEKASKDGEIHRGEAGKKRSVGKLSQTASEDSDVFGEADAIQNNGTSAEDTAVTDSKHTADPKNNWQGAHPADPGNRPHLIRLFSRDAPGREDNTFKDRPSESDELQTIQEDPTAASGGLDVMASQKRPSQRSKYLATASTMDHARHGFLPRHRDTGILDSIGRFFSGDRGAPKRGSGKVSLEHHHHHH
>RPB1 (Uniprot ID P24928, residues 1773-1970)
GHMSPNYTPTSPNYSPTSPSYSPTSPSYSPTSPSYSPSSPRYTPQSPTYTPSSPSYSPSSPSYSPASPKYTPTSPSYSPSSPEYTPTSPKYSPTSPKYSPTSPKYSPTSPTYSPTTPKYSPTSPTYSPTSPVYTPTSPKYSPTSPTYSPTSPKYSPTSPTYSPTSPKGSTYSPTSPGYSPTSPTYSLTSPAISPDDSDEEN
>securin (Uniprot ID O95997, residues 1-202)
MATLIYVDKENGEPGTRVVAKDGLKLGSGPSIKALDGRSQVSTPRFGKTFDAPPALPKATRKALGTVNRATEKSVKTKGPLKQKQPSFSAKKMTEKTVKAKSSVPASDDAYPEIEKFFPFNPLDFESFDLPEEHQIAHLPLSGVPLMILDEERELEKLFQLGPPSPVKMPSPPWESNLLQSPSSILSTLDVELPPVCCDIDI
>Sev-NT (Uniprot ID Q07097, residues 401-524)
LSGGDGAYHEPTGGGAIEVALDNADIDLETEAHADQDARGWGGESGERWARQVSGGHFVTLHGAERLEEETNDEDVSDIERRIAMRLAERRQEDSATHGDEGRNNGVDHDEDDDAAAVAGIGGI
>Sic1 (Uniprot ID P38634, residues 1-90)
GSMTPSTPPRSRGTRYLAQPSGNTSSSALMQGQKTPQKPSQNLVPVTPSTTKSFKNAPLLAPPNSNMGMTSPFNGLTSPQRSPFPKSSVKRT
>SKIPN (Uniprot ID G3V5R3, residues 59-129)
GDGGAFPEIHVAQYPLDMGRKKKMSNALAIQVDSEGKIKYDAIARQGQSKDKVIYSKYTDLVPKEVMNADD
>SLBP-NT (Uniprot ID Q9VAN6, residues 17-108)
MGSSHHHHHHSSGLVPRGSHMGSGSLNSSASSISIDVKPTMQSWAQEVRAEFGHSDEASSSLNSSAASCGSLAKKETADGNLESKDGEGREMAFEFLDGVNEVKFERLVKEEK
>α-synuclein (Uniprot ID P37840, residues 1-140)
MDVFMKGLSKAKEGVVAAAEKTKQGVAEAAGKTKEGVLYVGSKTKEGVVHGVATVAEKTKEQVTNVGGAVVTGVTAVAQKTVEGAGSIAAATGFVKKDQLGKNEEGAPQEGILEDMPVDPDNEAYEMPSEEGYQDYEPEA
>SOCS5-JIR (Uniprot ID A0A5E4BAI0, residues 12-81)
RSLRQRLQDTVGLCFPMRTYSKQSKPLFSNKRKIHLSELMLEKCPFPAGSDLAQKWHLIKQHTAPVSPHS
>tau K18 (Uniprot ID Q9MYX8, residues 186-314)
QTAPVPMPDLKNVKSKIGSTENLKHQPGGGKVQIINKKLDLSNVQSKCGSKDNIKHVPGGGSVQIVYKPVDLSKVTSKAGSLGNIHHKPGGGQVEVKSEKLDFKDRVQSKIGSLDNITHVPGGGNKKIE
>TC1 (Uniprot ID Q9NR00, residues 1-106)
MKAKRSHQAIIMSTSLRVSPSIHGYHFDTASRKKAVGNIFENTDQESLERLFRNSGDKKAEERAKIIFAIDQDVEEKTRALMALKKRTKDKLFQFLKLRKYSIKVH
>TDP-43 (Uniprot ID Q13148, residues 267-414)
GHMNRQLERSGRFGGNPGGFGNQGGFGNSRGGGAGLGNNQGSNMGGGMNFGAFSINPAMMAAAQAALQSSWGMMGMLASQQNQSGPSGNNQNQGNMQREPNQAFGSGNNSYSGSNSGAAIGWGSASNAGSGSGFNGGFGSSMDSKSSGWGM
>γ-tubulin-CT (Uniprot ID P53378, residues 439-473)
LLRGAAEQDSYLDDVLVDDENMVGELEEDLDADGDHKLV

Test set (9 IDPs)

>AMOTL1 (Uniprot ID Q8IY63, residues 178-384)
STQPQQNNEELPTYEEAKAQSQFFRGQQQQQQQQGAVGHGYYMAGGTSQKSRTEGRPTVNRANSGQAHKDEALKELKQGHVRSLSERIMQLSLERNGAKQHLPGSGNGKGFKVGGGPSPAQPAGKVLDPRGPPPEYPFKTKQMMSPVSKTQEHGLFYGDQHPGMLHEMVKPYPAPQPVRTDVAVLRYQPPPEYGVTSRPCQLPFPST
>CAHS-8 (Uniprot ID P0CU50, residues 1-227)
MSGRNVESHMERNEKVVVNNSGHADVKKQQQQVEHTEFTHTEVKAPLIHPAPPIISTGAAGLAEEIVGQGFTASAARISGGTAEVHLQPSAAMTEEARRDQERYRQEQESIAKQQEREMEKKTEAYRKTAEAEAEKIRKELEKQHARDVEFRKDLIESTIDRQKREVDLEAKMAKRELDREGQLAKEALERSRLATNVEVNFDSAAGHTVSGGTTVSTSDKMEIKRNENLYFQ
>ChiZ (Uniprot ID I6YA32, residues 1-64)
MTPVRPPHTPDPLNLRGPLDGPRWRRAEPAQSRRPGRSRPGGAPLRYHRTGVGMSRTGHGSRPV
>-endosulfine (Uniprot ID O43768, residues 1-121)
MSQKQEEENPAEETGEEKQDTQEKEGILPERAEEAKLKAKYPSLGQKPGGSDFLMKRLQKGQKYFDSGDYNMAKAKMKNKQLPSAGPDKNLVTGDHIPTPQDLPQRKSSLVTSKLAGGQVE
>FtsQ (Uniprot ID Q8IY63, residues 1-99)
MTEHNEDPQIERVADDAADEEAVTEPLATESKDEPAEHPEFEGPRRRARRERAERRAAQARATAIEQARRAAKRRARGQIVSEQNPAKPAARGVVRGLK
>Pdx1 (Uniprot ID P52945, residues 204-283)
GPGEEDKKRGGGTAVGGGGVAEPEQDCAVTSGEELLALPPPPPPGGAVPPAAPVAAREGRLPPGLSASPQPSSVAPRRPQEPR
>synaptobrevin-2 (Uniprot ID P63027, residues 1-96)
MSATAATAPPAAPAGEGGPPAPPPNLTSNRRLQQTQAQVDEVVDIMRVNVDKVLERDQKLSELDDRADALQAGASQFETSAAKLKRKYWWKNLKMM
>TIA-1 (Uniprot ID P31483, residues 320-386)
MGSSHHHHHHHHHHHHSENLYFQGGQYVPNGWQVPAYGVYGQPWSQQGFNQTQSSAPWMGPNYSVPPPQGQNGSMLPSQPAGYRVAGYETQ
>YAP (Uniprot ID P46937, residues 50-171)
AGHQIVHVRGDSETDLEALFNAVMNPKTANVPQTVPMRLRKLPDSFFKPPEPKSHSRQASTDAGTAGALTPQHVRAHSSPASLQLGAVSPGTLTPTGVVSGPAATPTAQHLRQSSFEIPDDV

Data availability

All data generated or analyzed during this study are included in the manuscript and supplementary files; source data have been provided for Figures 3—9, Figure 3—figure supplement 1, Figure 4—figure supplement 1, and Figure 5—figure supplement 1.

References

1. Abyzov A
2. Salvi N
3. Schneider R
4. Maurin D
5. Ruigrok RWH
6. Jensen MR
7. Blackledge M
(2016) Identification of dynamic modes in an intrinsically disordered protein using temperature-dependent NMR relaxation
Journal of the American Chemical Society 138:6240–6251.

https://doi.org/10.1021/jacs.6b02424
- PubMed
- Google Scholar
(2007) The BG21 isoform of Golli myelin basic protein is intrinsically disordered with a highly flexible amino-terminal domain
Biochemistry 46:9700–9712.

https://doi.org/10.1021/bi700632x
- PubMed
- Google Scholar
(2019) Concomitant disorder and high-affinity zinc binding in the human zinc- and iron-regulated transport protein 4 intracellular loop
Protein Science 28:868–880.

https://doi.org/10.1002/pro.3591
- PubMed
- Google Scholar
1. Ban D
2. Iconaru LI
3. Ramanathan A
4. Zuo J
5. Kriwacki RW
(2017) A small molecule causes a population shift in the conformational landscape of an intrinsically disordered protein
Journal of the American Chemical Society 139:13692–13700.

https://doi.org/10.1021/jacs.7b01380
- PubMed
- Google Scholar
1. Barré P
2. Eliezer D
(2013) Structural transitions in tau k18 on micelle binding suggest a hierarchy in the efficacy of individual microtubule-binding repeats in filament nucleation
Protein Science 22:1037–1048.

https://doi.org/10.1002/pro.2290
- PubMed
- Google Scholar
1. Bernadó P
2. Blackledge M
(2009) A self-consistent description of the conformational behavior of chemically denatured proteins from NMR and small angle scattering
Biophysical Journal 97:2839–2845.

https://doi.org/10.1016/j.bpj.2009.08.044
- PubMed
- Google Scholar
1. Borgia A
2. Borgia MB
3. Bugge K
4. Kissling VM
5. Heidarsson PO
6. Fernandes CB
7. Sottini A
8. Soranno A
9. Buholzer KJ
10. Nettels D
11. Kragelund BB
12. Best RB
13. Schuler B
(2018) Extreme disorder in an ultrahigh-affinity protein complex
Nature 555:61–66.

https://doi.org/10.1038/nature25762
- PubMed
- Google Scholar
1. Burke KA
2. Janke AM
3. Rhine CL
4. Fawzi NL
(2015) Residue-by-residue view of in vitro FUS granules that bind the C-terminal domain of RNA polymerase II
Molecular Cell 60:231–241.

https://doi.org/10.1016/j.molcel.2015.09.006
- PubMed
- Google Scholar
1. Camacho-Zarco AR
2. Schnapka V
3. Guseva S
4. Abyzov A
5. Adamski W
6. Milles S
7. Jensen MR
8. Zidek L
9. Salvi N
10. Blackledge M
(2022) NMR provides unique insight into the functional dynamics and interactions of intrinsically disordered proteins
Chemical Reviews 122:9331–9356.

https://doi.org/10.1021/acs.chemrev.1c01023
- PubMed
- Google Scholar
(2015) Structure and functional characterization of the conserved jak interaction region in the intrinsically disordered N-terminus of SOCS5
Biochemistry 54:4672–4682.

https://doi.org/10.1021/acs.biochem.5b00619
- PubMed
- Google Scholar
(2007) Amino acid bulkiness defines the local conformations and dynamics of natively unfolded alpha-synuclein and tau
Journal of the American Chemical Society 129:3032–3033.

https://doi.org/10.1021/ja067482k
- PubMed
- Google Scholar
1. Cho HY
2. Ul Mushtaq A
3. Lee JY
4. Kim DG
5. Seok MS
6. Jang M
7. Han B-W
8. Kim S
9. Jeon YH
(2014) Characterization of the interaction between lysyl-tRNA synthetase and laminin receptor by NMR
FEBS Letters 588:2851–2858.

https://doi.org/10.1016/j.febslet.2014.06.048
- PubMed
- Google Scholar
(2012) Effects of molecular crowding on the dynamics of intrinsically disordered proteins
PLOS ONE 7:e49876.

https://doi.org/10.1371/journal.pone.0049876
- PubMed
- Google Scholar
(2016) ALS mutations disrupt phase separation mediated by α-Helical Structure in the TDP-43 Low-Complexity C-terminal domain
Structure 24:1537–1549.

https://doi.org/10.1016/j.str.2016.07.007
- PubMed
- Google Scholar
(2019) Solution ensemble of the C-terminal domain from the transcription factor Pdx1 resembles an excluded volume polymer
The Journal of Physical Chemistry. B 123:106–116.

https://doi.org/10.1021/acs.jpcb.8b10051
- PubMed
- Google Scholar
1. Csizmok V
2. Felli IC
3. Tompa P
4. Banci L
5. Bertini I
(2008) Structural and dynamic characterization of intrinsically disordered human securin by NMR spectroscopy
Journal of the American Chemical Society 130:16873–16879.

https://doi.org/10.1021/ja805510b
- PubMed
- Google Scholar
1. De Avila M
2. Vassall KA
3. Smith GST
4. Bamm VV
5. Harauz G
(2014) The proline-rich region of 18.5 kDa myelin basic protein binds to the SH3-domain of Fyn tyrosine kinase with the aid of an upstream segment to form a dynamic complex in vitro
Bioscience Reports 34:e00157.

https://doi.org/10.1042/BSR20140149
- PubMed
- Google Scholar
1. Delaforge E
2. Kragelj J
3. Tengo L
4. Palencia A
5. Milles S
6. Bouvignies G
7. Salvi N
8. Blackledge M
9. Jensen MR
(2018) Deciphering the dynamic interaction profile of an intrinsically disordered protein by nmr exchange spectroscopy
Journal of the American Chemical Society 140:1148–1158.

https://doi.org/10.1021/jacs.7b12407
- PubMed
- Google Scholar
1. Dey S
2. MacAinsh M
3. Zhou HX
(2022) Sequence-dependent backbone dynamics of intrinsically disordered proteins
Journal of Chemical Theory and Computation 18:6310–6323.

https://doi.org/10.1021/acs.jctc.2c00328
- PubMed
- Google Scholar
1. Ebert MO
2. Bae SH
3. Dyson HJ
4. Wright PE
(2008) NMR relaxation study of the complex formed between CBP and the activation domain of the nuclear hormone receptor coactivator ACTR
Biochemistry 47:1299–1308.

https://doi.org/10.1021/bi701767j
- PubMed
- Google Scholar
1. Feichtinger M
2. Beier A
3. Migotti M
4. Schmid M
5. Bokhovchuk F
6. Chène P
7. Konrat R
(2022) Long-range structural preformation in yes-associated protein precedes encounter complex formation with TEAD
iScience 25:104099.

https://doi.org/10.1016/j.isci.2022.104099
- Google Scholar
1. Feldman HJ
2. Hogue CWV
(2002)
Probabilistic sampling of protein conformations: new hope for brute force?

Proteins 46:8–23.
- PubMed
- Google Scholar
1. Gruber T
2. Lewitzky M
3. Machner L
4. Weininger U
5. Feller SM
6. Balbach J
(2022) Macromolecular crowding induces a binding competent transient structure in intrinsically disordered Gab1
Journal of Molecular Biology 434:167407.

https://doi.org/10.1016/j.jmb.2021.167407
- PubMed
- Google Scholar
(2018) Concerted millisecond timescale dynamics in the intrinsically disordered carboxyl terminus of γ-tubulin induced by mutation of a conserved tyrosine residue
Protein Science 27:531–545.

https://doi.org/10.1002/pro.3345
- PubMed
- Google Scholar
1. Hicks A
2. Escobar CA
3. Cross TA
4. Zhou HX
(2020) Sequence-dependent correlated segments in the intrinsically disordered region of ChiZ
Biomolecules 10:1–23.

https://doi.org/10.3390/biom10060946
- PubMed
- Google Scholar
(2021) Removing thermostat distortions of protein dynamics in constant-temperature molecular dynamics simulations
Journal of Chemical Theory and Computation 17:5920–5932.

https://doi.org/10.1021/acs.jctc.1c00448
- PubMed
- Google Scholar
1. Janke AM
2. Seo DH
3. Rahmanian V
4. Conicella AE
5. Mathews KL
6. Burke KA
7. Mittal J
8. Fawzi NL
(2018) Lysines in the RNA polymerase II C-terminal domain contribute to TAF15 fibril recruitment
Biochemistry 57:2549–2563.

https://doi.org/10.1021/acs.biochem.7b00310
- PubMed
- Google Scholar
1. Jenner M
2. Kosol S
3. Griffiths D
4. Prasongpholchai P
5. Manzi L
6. Barrow AS
7. Moses JE
8. Oldham NJ
9. Lewandowski JR
10. Challis GL
(2018) Mechanism of intersubunit ketosynthase-dehydratase interaction in polyketide synthases
Nature Chemical Biology 14:270–275.

https://doi.org/10.1038/nchembio.2549
- PubMed
- Google Scholar
(2008) Quantitative conformational analysis of partially folded proteins from residual dipolar couplings: application to the molecular recognition element of Sendai virus nucleoprotein
Journal of the American Chemical Society 130:8055–8061.

https://doi.org/10.1021/ja801332d
- PubMed
- Google Scholar
1. Kiss R
2. Kovács D
3. Tompa P
4. Perczel A
(2008) Local structural preferences of calpastatin, the intrinsically unstructured protein inhibitor of calpain
Biochemistry 47:6936–6945.

https://doi.org/10.1021/bi800201a
- PubMed
- Google Scholar
1. Klein-Seetharaman J
2. Oikawa M
3. Grimshaw SB
4. Wirmer J
5. Duchardt E
6. Ueda T
7. Imoto T
8. Smith LJ
9. Dobson CM
10. Schwalbe H
(2002) Long-range interactions within a nonnative protein
Science 295:1719–1722.

https://doi.org/10.1126/science.1067680
- PubMed
- Google Scholar
(2019) Structural dynamics and transient lipid binding of synaptobrevin-2 tune SNARE assembly and membrane fusion
PNAS 116:8699–8708.

https://doi.org/10.1073/pnas.1813194116
- PubMed
- Google Scholar
1. Larkin MA
2. Blackshields G
3. Brown NP
4. Chenna R
5. McGettigan PA
6. McWilliam H
7. Valentin F
8. Wallace IM
9. Wilm A
10. Lopez R
11. Thompson JD
12. Gibson TJ
13. Higgins DG
(2007) Clustal W and clustal X version 2.0
Bioinformatics 23:2947–2948.

https://doi.org/10.1093/bioinformatics/btm404
- PubMed
- Google Scholar
1. Lawrence CW
2. Showalter SA
(2012) Carbon-detected (15)N NMR spin relaxation of an intrinsically disordered protein: FCP1 dynamics unbound and in complex with RAP74
The Journal of Physical Chemistry Letters 3:1409–1413.

https://doi.org/10.1021/jz300432e
- PubMed
- Google Scholar
(2020) Quantitative cooperative binding model for intrinsically disordered proteins interacting with nanomaterials
Journal of the American Chemical Society 142:10730–10738.

https://doi.org/10.1021/jacs.0c01885
- PubMed
- Google Scholar
1. Lopes FC
2. Dobrovolska O
3. Real-Guerra R
4. Broll V
5. Zambelli B
6. Musiani F
7. Uversky VN
8. Carlini CR
9. Ciurli S
(2015) Pliable natural biocide: Jaburetox is an intrinsically disordered insecticidal and fungicidal polypeptide derived from jack bean urease
The FEBS Journal 282:1043–1064.

https://doi.org/10.1111/febs.13201
- PubMed
- Google Scholar
1. Maiti S
2. Acharya B
3. Boorla VS
4. Manna B
5. Ghosh A
6. De S
(2019) Dynamic studies on intrinsically disordered regions of two paralogous transcription factors reveal rigid segments with important biological functions
Journal of Molecular Biology 431:1353–1369.

https://doi.org/10.1016/j.jmb.2019.02.021
- PubMed
- Google Scholar
(2022) Intrinsically disordered tardigrade proteins self-assemble into fibrous gels in response to environmental stress
Angewandte Chemie 61:e202109961.

https://doi.org/10.1002/anie.202109961
- PubMed
- Google Scholar
1. Martin EW
2. Holehouse AS
3. Grace CR
4. Hughes A
5. Pappu RV
6. Mittag T
(2016) Sequence determinants of the conformational properties of an intrinsically disordered protein prior to and upon multisite phosphorylation
Journal of the American Chemical Society 138:15323–15335.

https://doi.org/10.1021/jacs.6b10272
- PubMed
- Google Scholar
1. Martin EW
2. Holehouse AS
3. Peran I
4. Farag M
5. Incicco JJ
6. Bremer A
7. Grace CR
8. Soranno A
9. Pappu RV
10. Mittag T
(2020) Valence and patterning of aromatic residues determine the phase behavior of prion-like domains
Science 367:694–699.

https://doi.org/10.1126/science.aaw8653
- PubMed
- Google Scholar
(2020) The ambivalent role of proline residues in an intrinsically disordered protein: from disorder promoters to compaction facilitators
Journal of Molecular Biology 432:3093–3111.

https://doi.org/10.1016/j.jmb.2019.11.015
- PubMed
- Google Scholar
(2000) The PSIPRED protein structure prediction server
Bioinformatics 16:404–405.

https://doi.org/10.1093/bioinformatics/16.4.404
- PubMed
- Google Scholar
1. Milles S
2. Jensen MR
3. Lazert C
4. Guseva S
5. Ivashchenko S
6. Communie G
7. Maurin D
8. Gerlier D
9. Ruigrok RWH
10. Blackledge M
(2018) An ultraweak interaction in the intrinsically disordered replication machinery is essential for measles virus function
Science Advances 4:eaat7778.

https://doi.org/10.1126/sciadv.aat7778
- PubMed
- Google Scholar
1. Mittag T
2. Marsh J
3. Grishaev A
4. Orlicky S
5. Lin H
6. Sicheri F
7. Tyers M
8. Forman-Kay JD
(2010) Structure/function implications in a dynamic complex of the intrinsically disordered Sic1 with the Cdc4 subunit of an SCF ubiquitin ligase
Structure 18:494–506.

https://doi.org/10.1016/j.str.2010.01.020
- PubMed
- Google Scholar
(2011) Structural characterization of partially disordered human Chibby: insights into its function in the Wnt-signaling pathway
Biochemistry 50:715–726.

https://doi.org/10.1021/bi101236z
- PubMed
- Google Scholar
1. Murrali MG
2. Piai A
3. Bermel W
4. Felli IC
5. Pierattelli R
(2018) Proline fingerprint in intrinsically disordered proteins
Chembiochem 19:1625–1629.

https://doi.org/10.1002/cbic.201800172
- PubMed
- Google Scholar
(2019) Dynamics of the intrinsically disordered protein NUPR1 in isolation and in its fuzzy complexes with DNA and prothymosin α
Biochimica et Biophysica Acta - Proteins and Proteomics 1867:140252.

https://doi.org/10.1016/j.bbapap.2019.07.005
- Google Scholar
1. Olivieri C
2. Wang Y
3. Li GC
4. V S M
5. Kim J
6. Stultz BR
7. Neibergall M
8. Porcelli F
9. Muretta JM
10. Thomas DD
11. Gao J
12. Blumenthal DK
13. Taylor SS
14. Veglia G
(2020) Multi-state recognition pathway of the intrinsically disordered protein kinase inhibitor by protein kinase A
eLife 9:e55607.

https://doi.org/10.7554/eLife.55607
- PubMed
- Google Scholar
(2004) APE: analyses of phylogenetics and evolution in R language
Bioinformatics 20:289–290.

https://doi.org/10.1093/bioinformatics/btg412
- PubMed
- Google Scholar
(2014) Conformational propensities and dynamics of a βγ-crystallin, an intrinsically disordered protein
Physical Chemistry Chemical Physics 16:12703.

https://doi.org/10.1039/c3cp53558d
- Google Scholar
1. Piai A
2. Calçada EO
3. Tarenzi T
4. del Grande A
5. Varadi M
6. Tompa P
7. Felli IC
8. Pierattelli R
(2016) Just a flexible linker? the structural and dynamic properties of CBP-ID4 Revealed by NMR Spectroscopy
Biophysical Journal 110:372–381.

https://doi.org/10.1016/j.bpj.2015.11.3516
- Google Scholar
1. Qin S
2. Hicks A
3. Dey S
4. Prasad R
5. Zhou HX
(2022) ReSMAP: web server for predicting residue-specific membrane-association propensities of intrinsically disordered proteins
Membranes 12:773.

https://doi.org/10.3390/membranes12080773
- PubMed
- Google Scholar
Website
1. Qin S
2. Zhou HX
(2024) SeqDYN: Predictor for the Sequence-Dependent Backbone Dynamics of Intrinsically Disordered Proteins
Accessed October 23, 2024.

https://zhougroup-uic.github.io/SeqDYNidp/
(2018) Local and global dynamics in intrinsically disordered synuclein
Angewandte Chemie 57:15262–15266.

https://doi.org/10.1002/anie.201808172
- PubMed
- Google Scholar
(2019) reorientational dynamics of amyloid-β from NMR spin relaxation and molecular simulation
The Journal of Physical Chemistry Letters 10:3369–3375.

https://doi.org/10.1021/acs.jpclett.9b01050
- PubMed
- Google Scholar
1. Rives A
2. Meier J
3. Sercu T
4. Goyal S
5. Lin Z
6. Liu J
7. Guo D
8. Ott M
9. Zitnick CL
10. Ma J
11. Fergus R
(2021) Biological structure and function emerge from scaling unsupervised learning to 250 million protein sequences
PNAS 118:e2016239118.

https://doi.org/10.1073/pnas.2016239118
- PubMed
- Google Scholar
(2022) Molecular basis of small-molecule binding to α-Synuclein
Journal of the American Chemical Society 144:2501–2510.

https://doi.org/10.1021/jacs.1c07591
- PubMed
- Google Scholar
Software
1. Rohatgi A
(2022) Webplotdigitizer, version 4.6
WebPlotDigitizer.

https://apps.automeris.io/wpd4/
1. Schiavina M
2. Salladini E
3. Murrali MG
4. Tria G
5. Felli IC
6. Pierattelli R
7. Longhi S
(2020) Ensemble description of the intrinsically disordered N-terminal domain of the Nipah virus P/V protein from combined NMR and SAXS
Scientific Reports 10:19574.

https://doi.org/10.1038/s41598-020-76522-3
- PubMed
- Google Scholar
1. Schneider R
2. Maurin D
3. Communie G
4. Kragelj J
5. Hansen DF
6. Ruigrok RWH
7. Jensen MR
8. Blackledge M
(2015) Visualizing the molecular recognition trajectory of an intrinsically disordered protein using multinuclear relaxation dispersion NMR
Journal of the American Chemical Society 137:1220–1229.

https://doi.org/10.1021/ja511066q
- PubMed
- Google Scholar
1. Schwalbe H
2. Fiebig KM
3. Buck M
4. Jones JA
5. Grimshaw SB
6. Spencer A
7. Glaser SJ
8. Smith LJ
9. Dobson CM
(1997) Structural and dynamical properties of a denatured protein. Heteronuclear 3D NMR experiments and theoretical simulations of lysozyme in 8 M urea
Biochemistry 36:8977–8991.

https://doi.org/10.1021/bi970049q
- PubMed
- Google Scholar
(2002) Molecular hinges in protein folding: the urea-denatured state of apomyoglobin
Biochemistry 41:12681–12686.

https://doi.org/10.1021/bi020381o
- PubMed
- Google Scholar
1. Sekiyama N
2. Takaba K
3. Maki-Yonekura S
4. Akagi K-I
5. Ohtani Y
6. Imamura K
7. Terakawa T
8. Yamashita K
9. Inaoka D
10. Yonekura K
11. Kodama TS
12. Tochio H
(2022) ALS mutations in the TIA-1 prion-like domain trigger highly condensed pathogenic structures
PNAS 119:e2122523119.

https://doi.org/10.1073/pnas.2122523119
- PubMed
- Google Scholar
1. Shen Y
2. Delaglio F
3. Cornilescu G
4. Bax A
(2009) TALOS+: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts
Journal of Biomolecular NMR 44:213–223.

https://doi.org/10.1007/s10858-009-9333-z
- PubMed
- Google Scholar
1. Smrt ST
2. Escobar CA
3. Dey S
4. Cross TA
5. Zhou HX
(2023) An Arg/Ala-rich helix in the N-terminal region of M. tuberculosis FtsQ is a potential membrane anchor of the Z-ring
Communications Biology 6:311.

https://doi.org/10.1038/s42003-023-04686-5
- PubMed
- Google Scholar
1. Sólyom Z
2. Ma P
3. Schwarten M
4. Bosco M
5. Polidori A
6. Durand G
7. Willbold D
8. Brutscher B
(2015) The disordered region of the HCV Protein NS5A: conformational dynamics, SH3 binding, and phosphorylation
Biophysical Journal 109:1483–1496.

https://doi.org/10.1016/j.bpj.2015.06.040
- PubMed
- Google Scholar
1. Song J
2. Guo L-W
3. Muradov H
4. Artemyev NO
5. Ruoho AE
6. Markley JL
(2008) Intrinsically disordered gamma-subunit of cGMP phosphodiesterase encodes functionally relevant transient secondary and tertiary structure
PNAS 105:1505–1510.

https://doi.org/10.1073/pnas.0709558105
- PubMed
- Google Scholar
1. Sung YH
2. Eliezer D
(2007) Residual structure, backbone dynamics, and interactions within the synuclein family
Journal of Molecular Biology 372:689–707.

https://doi.org/10.1016/j.jmb.2007.07.008
- PubMed
- Google Scholar
(2011) Full backbone assignment and dynamics of the intrinsically disordered dehydrin ERD14
Biomolecular NMR Assignments 5:189–193.

https://doi.org/10.1007/s12104-011-9297-2
- Google Scholar
1. Tesei G
2. Lindorff-Larsen K
(2022) Improved predictions of phase behaviour of intrinsically disordered proteins by tuning the interaction range
Open Research Europe 2:94.

https://doi.org/10.12688/openreseurope.14967.2
- PubMed
- Google Scholar
1. Thapa C
2. Roivas P
3. Haataja T
4. Permi P
5. Pentikäinen U
(2022) Interaction mechanism of endogenous PP2A inhibitor protein ENSA with PP2A
The FEBS Journal 289:519–534.

https://doi.org/10.1111/febs.16150
- PubMed
- Google Scholar
(2004) The N-terminal domain of the Drosophila histone mRNA binding protein, SLBP, Is Intrinsically disordered with nascent helical structure ^,
Biochemistry 43:9390–9400.

https://doi.org/10.1021/bi036314r
- Google Scholar
(2022) Multivalent Angiomotin-like 1 and Yes-associated protein form a dynamic complex
Protein Science 31:e4295.

https://doi.org/10.1002/pro.4295
- PubMed
- Google Scholar
1. Wang X
2. Zhang S
3. Zhang J
4. Huang X
5. Xu C
6. Wang W
7. Liu Z
8. Wu J
9. Shi Y
(2010) A large intrinsically disordered region in SKIP and its disorder-order transition induced by PPIL1 binding revealed by NMR
The Journal of Biological Chemistry 285:4951–4963.

https://doi.org/10.1074/jbc.M109.087528
- PubMed
- Google Scholar
1. Wang J
2. Choi J-M
3. Holehouse AS
4. Lee HO
5. Zhang X
6. Jahnel M
7. Maharana S
8. Lemaitre R
9. Pozniakovsky A
10. Drechsel D
11. Poser I
12. Pappu RV
13. Alberti S
14. Hyman AA
(2018) A molecular grammar governing the driving forces for phase separation of prion-like RNA binding proteins
Cell 174:688–699.

https://doi.org/10.1016/j.cell.2018.06.006
- PubMed
- Google Scholar
(2006) Motional properties of unfolded ubiquitin: a model for a random coil protein
Journal of Biomolecular NMR 35:175–186.

https://doi.org/10.1007/s10858-006-9026-9
- PubMed
- Google Scholar
1. Wong LE
2. Kim TH
3. Muhandiram DR
4. Forman-Kay JD
5. Kay LE
(2020) NMR experiments for studies of dilute and condensed protein phases: application to the phase-separating protein CAPRIN1
Journal of the American Chemical Society 142:2471–2489.

https://doi.org/10.1021/jacs.9b12208
- PubMed
- Google Scholar
1. Xie M
2. Li D
3. Yuan J
4. Hansen AL
5. Brüschweiler R
(2018) Quantitative binding behavior of intrinsically disordered proteins to nanoparticle surfaces at individual residue level
Chemistry – A European Journal 24:16997–17001.

https://doi.org/10.1002/chem.201804556
- Google Scholar
(2016) Characterisation of the conformational preference and dynamics of the intrinsically disordered N-terminal region of Beclin 1 by NMR spectroscopy
Biochimica et Biophysica Acta 1864:1128–1137.

https://doi.org/10.1016/j.bbapap.2016.06.005
- PubMed
- Google Scholar
1. Yu L
2. Brüschweiler R
(2022) Quantitative prediction of ensemble dynamics, shapes and contact propensities of intrinsically disordered proteins
PLOS Computational Biology 18:e1010036.

https://doi.org/10.1371/journal.pcbi.1010036
- Google Scholar
1. Zhang Y
2. Prasad R
3. Su S
4. Lee D
5. Zhou HX
(2024) Amino acid-dependent phase equilibrium and material properties of tetrapeptide condensates
Cell Reports Physical Science 5:102218.

https://doi.org/10.1016/j.xcrp.2024.102218
- Google Scholar
1. Zheng Z
2. Ma D
3. Yahr TL
4. Chen L
(2012) The transiently ordered regions in intrinsically disordered ExsE are correlated with structural elements involved in chaperone binding
Biochemical and Biophysical Research Communications 417:129–134.

https://doi.org/10.1016/j.bbrc.2011.11.070
- PubMed
- Google Scholar
(1968) The characterization of amino acid sequences in proteins by statistical methods
Journal of Theoretical Biology 21:170–201.

https://doi.org/10.1016/0022-5193(68)90069-6
- PubMed
- Google Scholar

Article and author information

Author details

Sanbo Qin

Department of Chemistry, University of Illinois Chicago, Chicago, United States

Contribution
Resources, Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Methodology

Competing interests
No competing interests declared
Huan-Xiang Zhou
1. Department of Chemistry, University of Illinois Chicago, Chicago, United States
2. Department of Physics, University of Illinois Chicago, Chicago, United States
Contribution
Conceptualization, Resources, Data curation, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing

For correspondence
hzhou43@uic.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-9020-0302

Funding

National Institutes of Health (GM118091)

Huan-Xiang Zhou

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work was supported by Grant GM118091 from the National Institutes of Health.

Version history

Preprint posted: February 3, 2023
Sent for peer review: May 9, 2023
Reviewed Preprint version 1: July 18, 2023
Reviewed Preprint version 2: September 17, 2024
Version of Record published: October 30, 2024

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.88958. This DOI represents all versions, and will always resolve to the latest one.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

757

views
49

downloads
2

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Sanbo Qin
Huan-Xiang Zhou

(2024)

Predicting the sequence-dependent backbone dynamics of intrinsically disordered proteins

eLife 12:RP88958.

https://doi.org/10.7554/eLife.88958.3

Share this article

Cite this article

Clock-like tree plot showing lack of homology among the 45 IDPs.

Representative conformations of five IDPs.

Experimental conditions, mean and standard deviation of measured R2, and SeqDYN prediction RMSE.

Properties of the 45 IDPs in the training set.

Figure 3—source data 1

SeqDYN model parameters.

Figure 4—source data 1

Quality of SeqDYN predictions.

Figure 5—source data 1

Measured (bars) and predicted (curves) R2 profiles for ChiZ N-terminal region, TIA1 prion-like domain, Pdx1 C-terminal region, synaptobrevin-2, α-endosulfine, YAP, AMOTL1, FtsQ, and CAHS-8.

Figure 6—source data 1

R2 profiles predicted (curves) by SeqDYN show close agreement with those measured (bars) on structured proteins in the unfolded state.

Figure 7—source data 1

Comparison between SeqDYN prediction (curves) and effective transverse relaxation rate (bars) from 1H dispersion relaxation experiment.

Figure 8—source data 1

Correlation between the stickiness parameters (λ) and the NMR relaxation parameters (q).

Figure 9—source data 1

RMSEs (s–1) of R2 predictions by SeqDYN and MD for 10 IDPs.

Author details

Sanbo Qin

Contribution

Competing interests

Huan-Xiang Zhou

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading

Experimental conditions, mean and standard deviation of measured $R_{2}$ , and SeqDYN prediction RMSE.

Measured (bars) and predicted (curves) $R_{2}$ profiles for ChiZ N-terminal region, TIA1 prion-like domain, Pdx1 C-terminal region, synaptobrevin-2, α-endosulfine, YAP, AMOTL1, FtsQ, and CAHS-8.

$R_{2}$ profiles predicted (curves) by SeqDYN show close agreement with those measured (bars) on structured proteins in the unfolded state.

Comparison between SeqDYN prediction (curves) and effective transverse relaxation rate (bars) from ¹H dispersion relaxation experiment.

RMSEs (s^–1) of $R_{2}$ predictions by SeqDYN and MD for 10 IDPs.