Structural and Functional Characterization of the ABA-Water Deficit Stress Domain from Wheat and Barley: An Intrinsically Disordered Domain behind the Versatile Functions of the Plant Abscissic Acid, Stress and Ripening Protein Family

The ASR protein family has been discovered thirty years ago in many plant species and is involved in the tolerance of various abiotic stresses such as dehydration, salinity and heat. Despite its importance, nothing is known about the conserved ABA-Water Deficit Stress Domain (ABA-WDS) of the ASR gene family. In this study, we characterized two ABA-WDS domains, isolated from durum wheat (TtABA-WDS) and barley (HvABA-WDS). Bioinformatics analysis shows that they are both consistently predicted to be intrinsically disordered. Hydrodynamic and circular dichroism analysis indicate that both domains are largely disordered but belong to different structural classes, with HvABA-WDS and TtABA-WDS adopting a PreMolten Globule-like (PMG-like) and a Random Coil-like (RC-like) conformation, respectively. In the presence of the secondary structure stabilizer trifluoroethanol (TFE) or of increasing glycerol concentrations, which mimics dehydration, the two domains acquire an α-helical structure. Interestingly, both domains are able to prevent heat- and dehydration-induced inactivation of the enzyme lactate dehydrogenase (LDH). Furthermore, heterologous expression of TtABA-WDS and HvABA-WDS in the yeast Saccharomyces cerevisiae improves its tolerance to salt, heat and cold stresses. Taken together our results converge to show that the ABA-WDS domain is an intrinsically disordered functional domain whose conformational plasticity could be instrumental to support the versatile functions attributed to the ASR family, including its role in abiotic stress tolerance. Finally, and after validation in the plant system, this domain could be used to improve crop tolerance to abiotic stresses.


Introduction
Crop production must increase to meet population feed global demand that is estimated to increase continuously. It is well known that abiotic stresses such as extremes in temperature, drought, salinity, heavy metal concentrations, and radiation represent the most limiting factors for agricultural productivity worldwide [1]. To cope with such challenges, plants developed some adaptation strategies to escape the changing environmental conditions. Abiotic stress responses in plants are complex due to the interrelationship of mechanisms and variability of genes involved. Gene discovery, with the genomic resources developed recently for wheat, has enabled rapid assessment of gene function and directly accelerates efficient wheat breeding. previously reported [9]) argues for a functional significance. Noteworthy, in some cases, as for instance in the ASR protein from rice (Os ASR1), the domain represents more than 65% of the full-length protein [27].
Domains represent one of the most useful levels to understand protein function, and domain family-based analysis has a profound impact on the functional characterization and annotation of recently released genome data. To the best of our knowledge, the conserved ABA-WDS domain of this protein family has never been experimentally investigated. Although the ASR gene family has been extensively explored in many plant species, the precise physiological function, mechanism of action and classification of ASR proteins remain largely unknown. In addition, experimental and functional data on the ABA-WDS domain in isolation are lacking, while they are of paramount importance to illuminate the function of the ASR family in plant stress tolerance breeding. In order to fill in this gap in knowledge, and with the aim of ascertaining whether the ABA-WDS domain is functional on its own, we herein report the identification, structural and functional characterization of two ABA-WDS domains, TtABA-WDS and HvABA-WDS, isolated respectively, from salt-stressed durum wheat and barley leaves. The rationale for focusing on durum wheat and barley resides in pursuing and extending our previous studies that reported the characterization of two ASR proteins from these two cereals [17] and that showed the role of the ASR family in salt and drought tolerance in durum wheat [9].
Our results reveal that these domains are consistently predicted to be intrinsically disordered and further confirmed their disordered nature experimentally. In addition, we investigated the conformational plasticity of these two domains under external environmental factors, such as temperature and glycerol-induced water deficit. Finally, we show that the two domains are able to prevent the inactivation of the enzyme lactate dehydrogenase (LDH) under heat and desiccation stress and we provide evidence that their heterologous expression in yeast improves its tolerance to salt and heat-cold stress.

Cloning of TtABA-WDS and HvABA-WDS cDNA and Disorder Prediction of the Encoded Protein Domains
Using primers designed from a previously reported multiple sequence alignment of plant full-length ASR proteins and conceived to amplify ABA-WDS domains [9], two cDNAs were amplified, one from salt-stressed durum wheat (Tt) and one from barley (Hv) leaves. The amplified cDNAs were subsequently cloned and sequenced. Analysis of the two deduced amino acid sequences against InterPro and Conserved Domains Database (CDD) confirmed that the two sequences (hereafter referred to as TtABA-WDS and HvABA-WDS) correspond to ABA-WDS domains.
TtABA-WDS and HvABA-WDS cDNA encode polypeptides of 66 amino acids in length with a predicted molecular mass of~7.2 kDa. The TtABA-WDS and HvABA-WDS sequences are very close to each other (94% of identity) and only differ in four amino acids (10H-> R; 17A-> V; 45V-> L; 63Q-> T, Figure 1A). The two domains' acidic isoelectric points are 6.12 and 6.24, respectively (Table 1). Their predicted grand average of hydropathy index (GRAVY) values (as obtained from http://www.gravy-calculator.de/ (accessed on 25 January 2021)) are negative, underscoring their hydrophilic nature (Table 1). : residue number; MM: molecular mass in Daltons; f+: fraction of positively charged residues; f−: fraction of negatively charged residues; FCR: fraction of charged residues (f+ + f−); NCPR: net charge per residue, value of the difference between the fraction of positively charged and negatively positively residues, (|f+ − f−|); κ: linear distribution of opposite charges; GRAVY: grand average of hydropathy index (sum of hydropathy values of all amino acids divided by the protein length); PDR: phase diagram region. Basic residues are in blue, acidic residues in red, prolines in pink, histidines in green and aromatic residues in orange. Amino acid substitutions are shown in bold on a yellow background. The alignment was generated with Clustal Omega and then edited manually. (B) Deviation in amino acid composition from the Swiss-PROT database of TtABA-WDS and HvABA-WDS. The relative enrichment in disorder promoting (pink and green bars) and depletion in order-promoting (grey bars) residues is shown. Residues have been ordered on the x-axis according to the TOP-IDP flexibility index as described in [28].
The sequence composition of TtABA-WDS and HvABA-WDS domains were compared to that of proteins within the SWISS-PROT database ( Figure 1B). Both domains have a biased sequence composition, being enriched in the most disorder-promoting amino acids (A, G, R, D, H, Q, K, S, E and P) and depleted in order-promoting residues (W, F, Y, I, M, L, V, C and T) [28]. Moreover, Ala (A), Glu (E) and His (H) are the most abundant amino acids in both domains. In line with these observations, TtABA-WDS and HvABA-WDS domains are predicted to be intrinsically disordered by the mean hydrophobicity/mean net charge ratio [29], as judged from their location in the left-hand side of the RH-plot (Supplementary Figure S1A,B).
The group of Pappu showed that sequence polarity is a determinant of IDP compaction, with polar IDPs having been found to favor collapsed ensembles in water de- Basic residues are in blue, acidic residues in red, prolines in pink, histidines in green and aromatic residues in orange. Amino acid substitutions are shown in bold on a yellow background. The alignment was generated with Clustal Omega and then edited manually. (B) Deviation in amino acid composition from the Swiss-PROT database of TtABA-WDS and HvABA-WDS. The relative enrichment in disorder promoting (pink and green bars) and depletion in order-promoting (grey bars) residues is shown. Residues have been ordered on the x-axis according to the TOP-IDP flexibility index as described in [28]. The sequence composition of TtABA-WDS and HvABA-WDS domains were compared to that of proteins within the SWISS-PROT database ( Figure 1B). Both domains have a biased sequence composition, being enriched in the most disorder-promoting amino acids (A, G, R, D, H, Q, K, S, E and P) and depleted in order-promoting residues (W, F, Y, I, M, L, V, C and T) [28]. Moreover, Ala (A), Glu (E) and His (H) are the most abundant amino acids in both domains. In line with these observations, TtABA-WDS and HvABA-WDS domains are predicted to be intrinsically disordered by the mean hydrophobicity/mean net charge ratio [29], as judged from their location in the left-hand side of the RH-plot (Supplementary Figure S1A,B).
The group of Pappu showed that sequence polarity is a determinant of IDP compaction, with polar IDPs having been found to favor collapsed ensembles in water despite the absence of hydrophobic groups [30]. To assess sequence polarity, the net charge per residue (NCPR, defined as |f + − f − |) [30,31], the total fraction of charged residues (FCR, defined as f + + f − ), and linear distribution of opposite charges (κ value) [32] were calculated (Table 1). For both domains, the NCPR is negative, and the FCR values of TtABA-WDS and HvABA-WDS domains are close (0.333 and 0.348, respectively). According to the predictive diagram of states developed by Pappu and colleagues [33], both domains fall in phase diagram region (PDR) 2, with both being however very close to the boundary separating PDR2 and PDR3 (Figure 2A,C). PDR 2 corresponds to proteins that adopt conformations likely representing a continuum of possibilities between PDR 1, which accommodates weak polyampholytes or polyelectrolytes adopting a globular conformation, and PDR 3, which embeds strong polyampholytes adopting non-globular, swollen coil-like conformations. Both TtABA-WDS and HvABA-WDS domains have very low κ values indicating that opposite charges are well mixed and suggesting that these two domains adopt preferentially extended and swollen, coil-like conformation. spite the absence of hydrophobic groups [30]. To assess sequence polarity, the net charge per residue (NCPR, defined as |f+ − f−|) [30,31], the total fraction of charged residues (FCR, defined as f+ + f−), and linear distribution of opposite charges (κ value) [32] were calculated (Table 1). For both domains, the NCPR is negative, and the FCR values of TtABA-WDS and HvABA-WDS domains are close (0.333 and 0.348, respectively). According to the predictive diagram of states developed by Pappu and colleagues [33], both domains fall in phase diagram region (PDR) 2, with both being however very close to the boundary separating PDR2 and PDR3 (Figure 2A,C). PDR 2 corresponds to proteins that adopt conformations likely representing a continuum of possibilities between PDR 1, which accommodates weak polyampholytes or polyelectrolytes adopting a globular conformation, and PDR 3, which embeds strong polyampholytes adopting non-globular, swollen coil-like conformations. Both TtABA-WDS and HvABA-WDS domains have very low κ values indicating that opposite charges are well mixed and suggesting that these two domains adopt preferentially extended and swollen, coil-like conformation.  [34] ouputs of TtABA-WDS (B) and HvABA-WDS (D) featuring the amino acid sequence above the HCA plot. The regions predicted as disordered by IUPRED [35] are represented by bi-directional red arrows.
The amino acid sequences of HvABA-WDS and TtABA-WDS were analyzed also using the MeDor server for predicting protein disorder [34]. Hydrophobic Cluster Analysis (HCA) [36] (see Supplementary text) of the two domains shows that they are depleted in hydrophobic clusters ( Figure 2B,D), a property indicative of protein disorder. The HCA plots of the two domains are very similar except for the size of the first hy-  [34] ouputs of TtABA-WDS (B) and HvABA-WDS (D) featuring the amino acid sequence above the HCA plot. The regions predicted as disordered by IUPRED [35] are represented by bi-directional red arrows.
The amino acid sequences of HvABA-WDS and TtABA-WDS were analyzed also using the MeDor server for predicting protein disorder [34]. Hydrophobic Cluster Analysis (HCA) [36] (see Supplementary text) of the two domains shows that they are depleted in hydrophobic clusters ( Figure 2B,D), a property indicative of protein disorder. The HCA plots of the two domains are very similar except for the size of the first hydrophobic cluster that encompasses two hydrophobic residues (LeuMet) in the case of TtABA-WDS, and three hydrophobic residues (LeuMetVal) in the case of HvABA-WDS ( Figure 2B,D). This difference arises from the 17A-> V substitution. Additionally, the number of disordered residues predicted by IUPRED [35] is larger in TtABA-WDS than in HvABA-WDS ( Figure 2B,D).
Analysis of the amino acid sequence of the two domains with the Phyre2 protein fold recognition server [37] returned a very low-confidence (~16%) model encompassing residues A 30 -V 52 , in which residues H 40 -A 51 adopt an α-helical conformation (data not shown). Failure of the server to predict a model over most of the sequence, together with the very low confidence of the predicted α-helix, is consistent with the expectedly predominantly disordered nature of the domains.

Expression and Purification of HvABA-WDS and TtABA-WDS Domains
TtABA-WDS and HvABA-WDS cDNAs were cloned into the pGEX4-T1 expression vector to generate N-terminally glutathione S transferase (GST) tagged proteins. The two recombinant plasmids were then used to transform the E. coli strain BL21, where they were expressed as a major protein product (data not shown).
They were purified from the soluble fraction by GST affinity chromatography, followed by thrombin cleavage to remove the GST tag and size-exclusion chromatography (SEC). The final purified products ( Figure 3A) exhibit an abnormally slow migration in SDS-PAGE: the apparent molecular mass (MM) of the recombinant domains, is estimated to be about 12 kDa, i.e., 1.3 times larger than that predicted by the amino acid sequence (~9.2 kDa, including the vector-encoded residues). Note that for both proteins, MALDI-TOF mass spectrometry (MS) analysis revealed the presence of a unique species with a molecular mass very close to the value calculated from the amino acid sequence ( Figure 3B,C). In addition, MS analysis of peptides resulting from GluC digestion confirmed the identity of the purified proteins (Supplementary Figure S2). From MS analysis the purity of the proteins was estimated to be >95%.
IDPs often display an aberrant electrophoretic migration that is mainly due to their typical compositional bias, e.g., high net charge and low hydrophobicity that causes them to bind less SDS than globular proteins [38,39]. The anomalous electrophoretic behavior of TtABA-WDS and HvABA-WDS constitutes a hint of their likely disordered nature.

Hydrodynamiques Properties of HvABA-WDS and TtABA-WDS Domains from Size Exclusion Chromatography (SEC)
The Stokes radii (R S of TtABA-WDS and HvABA-WDS domains were inferred using SEC. The elution profile of proteins in fact reflects their hydrodynamic properties. SEC was carried out using PBS pH 7.3 as elution buffer. TtABA-WDS and HvABA-WDS domains were eluted from the gel filtration column as sharp peaks (data not shown). From their elution volume (75.5 mL for HvABA-WDS and 72.1 mL for TtABA-WDS), their R S were estimated to be 20.3 Å for HvABA-WDS and 21.5 Å for TtABA-WDS (Table 2).  IDPs often display an aberrant electrophoretic migration that is mainly due to their typical compositional bias, e.g., high net charge and low hydrophobicity that causes them to bind less SDS than globular proteins [38,39]. The anomalous electrophoretic behavior of TtABA-WDS and HvABA-WDS constitutes a hint of their likely disordered nature.

Hydrodynamiques Properties of HvABA-WDS and TtABA-WDS Domains from Size Exclusion Chromatography (SEC)
The Stokes radii (RS of TtABA-WDS and HvABA-WDS domains were inferred using SEC. The elution profile of proteins in fact reflects their hydrodynamic properties. SEC was carried out using PBS pH 7.3 as elution buffer. TtABA-WDS and HvABA-WDS domains were eluted from the gel filtration column as sharp peaks (data not shown). From their elution volume (75.5 mL for HvABA-WDS and 72.1 mL for TtABA-WDS), their RS were estimated to be 20.3 Å for HvABA-WDS and 21.5 Å for TtABA-WDS (Table 2).

Differential Scanning Fluorimetry of HvABA-WDS and TtABA-WDS Domains
TtABA-WDS and HvABA-WDS were further analyzed by differential scanning fluorimetry (DSF), a technique used to determine the conformational stability of proteins that relies on the use of a fluorescent dye whose fluorescence is enhanced when protein hydrophobic cavities become accessible as a result of temperature-induced protein unfolding. IDPs, which do not possess hydrophobic cavities being devoid of a stable 3D structure, are characterized by a fluorescence profile that is rather flat and temperature-independent [17]. The experimentally observed profile obtained for HvABA-WDS is flat, consistent with the lack of a stable 3D structure ( Figure 4). Albeit the profile of TtABA-WDS is not flat and rather features a progressive decrease in the fluorescence with increasing temperature, it still lacks the typical transition peak observed for a folded protein undergoing a cooperative transition ( Figure 4). As such, the observed profile is indicative of the lack of 3D structure like in the case of HvABA-WDS. Therefore these results strengthen the conclusion that in silico analyses and with the above-described experimental lines of evidence.

Circular Dichroism (CD) Studies of HvABA-WDS and TtABA-WDS
To assess the secondary structure content of TtABA-WDS and HvABA-WDS, circular dichroism (CD) measurements were performed in the far-ultraviolet (UV) region. Under native conditions and neutral pH, both spectra display a large negative peak centered at 200 nm, low amplitude in the 220-230 nm region, and low ellipticity at 190 nm ( Figure 5A). Both domains present spectra typical of disordered proteins lacking any stable organized secondary structure. Spectral deconvolution revealed a high content (about 50%) of unordered structure in both domains ( Figure 5A, inset).
Previous studies unveiled that far-UV CD spectroscopy enables discriminating between IDPs adopting a random coil-like (RC) state and a premolten globule-like (PMG) state, based on the ratio of the ellipticity values at 200 and 222 nm [40]. PMGs possess some transiently populated secondary structure and are thus more structured than RCs. According to their ellipticity values at 200 and 222 nm, TtABA-WDS and HvABA-WDS domains can be classified in the RC-like and PMG-like region, respectively ( Figure 5B). Thus, HvABA-WDS is more compact than TtABA-WDS, in agreement with SEC data that highlighted a larger compaction index (CI) for HvABA-WDS compared to TtABA-WDS (Table 2).
We next assessed whether the two domains undergo temperature-induced folding. To this end, we monitored the ellipticity at 230 nm as a function of temperature ( Figure  5C). While TtABA-WDS undergoes temperature-induced folding, the profile of HvABA-WDS is roughly linear and devoid of a clear temperature-dependence ( Figure  5C). To rule out any possible temperature-induced aggregation of the proteins, we monitored the HT voltage associated with the transmission of light through the protein in the CD spectropolarimeter during the heating of the proteins. The HT voltage is in fact extremely sensitive to changes in the level of scattering of transmitted light. Aggregation of protein at high temperature will lead to an increase in the HT voltage. We observed no significant changes in the HT voltage associated with heating indicating a lack of aggre-

Circular Dichroism (CD) Studies of HvABA-WDS and TtABA-WDS
To assess the secondary structure content of TtABA-WDS and HvABA-WDS, circular dichroism (CD) measurements were performed in the far-ultraviolet (UV) region. Under native conditions and neutral pH, both spectra display a large negative peak centered at 200 nm, low amplitude in the 220-230 nm region, and low ellipticity at 190 nm ( Figure 5A). Both domains present spectra typical of disordered proteins lacking any stable organized secondary structure. Spectral deconvolution revealed a high content (about 50%) of unordered structure in both domains ( Figure 5A, inset).
Previous studies unveiled that far-UV CD spectroscopy enables discriminating between IDPs adopting a random coil-like (RC) state and a premolten globule-like (PMG) state, based on the ratio of the ellipticity values at 200 and 222 nm [40]. PMGs possess some transiently populated secondary structure and are thus more structured than RCs. According to their ellipticity values at 200 and 222 nm, TtABA-WDS and HvABA-WDS domains can be classified in the RC-like and PMG-like region, respectively ( Figure 5B). Thus, HvABA-WDS is more compact than TtABA-WDS, in agreement with SEC data that highlighted a larger compaction index (CI) for HvABA-WDS compared to TtABA-WDS (Table 2).
We next assessed whether the two domains undergo temperature-induced folding. To this end, we monitored the ellipticity at 230 nm as a function of temperature ( Figure 5C). While TtABA-WDS undergoes temperature-induced folding, the profile of HvABA-WDS is roughly linear and devoid of a clear temperature-dependence ( Figure 5C). To rule out any possible temperature-induced aggregation of the proteins, we monitored the HT voltage associated with the transmission of light through the protein in the CD spectropolarimeter during the heating of the proteins. The HT voltage is in fact extremely sensitive to changes in the level of scattering of transmitted light. Aggregation of protein at high temperature will lead to an increase in the HT voltage. We observed no significant changes in the HT voltage associated with heating indicating a lack of aggregation during heating (data not shown). In addition, we also monitored the absorbance at 340 nm (indicative of the presence of aggregates [41]) before and after heating the samples at 80 • C and observed no differences. We also monitored the ellipticity at 230 nm upon cooling the sample from 80 • C to 20 • C. The obtained profiles are almost indistinguishable from the heating profiles ( Figure 5C) thereby attesting to the full reversibility of the process and hence definitively excluding possible irreversible aggregation-related processes.
To provide more evidence on the conformational state of TtABA-WDS and HvABA-WDS domains, we analyzed the environment of the aromatic amino acid side chains, which provides information about the tertiary structure, by measuring the CD spectra in the near-UV region (250-350 nm) at 20 • C [42]. Although the near-UV CD spectra of TtABA-WDS and HvABA-WDS are not superimposable, they are both characterized by a weak intensity and by the presence of a broad negative band ( Figure 5D), once again consistent with the lack of a stable tertiary structure and of a hydrophobic core containing oriented aromatic residues as typically found in globular proteins [43]. The differences observed between the two near-UV CD spectra likely reflect differences between the two domains that are also pinpointed by far-UV CD studies. gation during heating (data not shown). In addition, we also monitored the absorbance at 340 nm (indicative of the presence of aggregates [41]) before and after heating the samples at 80 °C and observed no differences. We also monitored the ellipticity at 230 nm upon cooling the sample from 80 °C to 20 °C. The obtained profiles are almost indistinguishable from the heating profiles ( Figure 5C) thereby attesting to the full reversibility of the process and hence definitively excluding possible irreversible aggregation-related processes.
To provide more evidence on the conformational state of TtABA-WDS and HvABA-WDS domains, we analyzed the environment of the aromatic amino acid side chains, which provides information about the tertiary structure, by measuring the CD spectra in the near-UV region (250-350 nm) at 20 °C [42]. Although the near-UV CD spectra of TtABA-WDS and HvABA-WDS are not superimposable, they are both characterized by a weak intensity and by the presence of a broad negative band ( Figure 5D), once again consistent with the lack of a stable tertiary structure and of a hydrophobic core containing oriented aromatic residues as typically found in globular proteins [43]. The differences observed between the two near-UV CD spectra likely reflect differences between the two domains that are also pinpointed by far-UV CD studies.

TtABA-WDS and HvABA-WDS Domains Have an Intrinsic Propensity to Fold into an α-Helical Conformation
Some of the most striking functional advantages of IDPs come from their ability to contextually change their conformation due to binding or environmental changes [44]. In order to assess the folding potential of TtABA-WDS and HvABA-WDS domains, their far-UV CD spectra were recorded in the presence of TFE. This organic solvent is a secondary structure stabilizer that is used to mimic the hydrophobic environment experienced by proteins during interactions [45,46] It is thus widely used as a probe of hidden structural propensities to unveil protein regions having a propensity to undergo induced folding [39,47]. For both domains, the addition of TFE induces a gain of α-helicity indicated by the characteristic maximum at 190 nm and double minima at 208 and 222 nm ( Figure 6A,B). The estimated α-helical content gradually increases upon increasing the TFE concentration from 0 to 50% ( Figure 6A,B and insets). For both proteins, the most important transitions take place at 30% of TFE. We noted also that for the two domains, the spectra display an isodichroic point at 202 nm indicative of a two-state transition ( Figure 6A,B). We therefore plotted the percentage of α-helix as a function of the TFE concentration and fitted the data to a sigmoidal curve, corresponding to a two-state transition ( Figure 6C). The fitting yielded a midpoint of transition at a TFE concentration of~27% for both domains. From the fitting procedure, the m value could also be obtained. The latter is a measure of the osmolyte efficacy in folding or unfolding a protein. It also reflects the cooperativity of folding/unfolding in the presence of the osmolyte, and arises from the free energy contributions of protein groups that either become exposed upon unfolding or buried upon folding the protein.
The m value is typically positive for protecting osmolytes, such as TFE, that drive the equilibrium towards the folded state. In agreement, m values of 0.049 kcal mol −1 M −1 and of 0.042 kcal mol −1 M −1 were obtained for TtABA-WDS and HvABA-WDS, respectively. The two domains therefore feature midpoints and m values very close to each other, indicating that they follow a very similar folding pathway.

TtABA-WDS and HvABA-WDS Domains Gain α-Helicity under Conditions Mimicking Dehydration
Glycerol is known as a simplistic model for mimicking cellular dehydration that induces changes in the water surface tension around proteins. Here, we sought at assessing whether the two cereal domains could gain structure under increasing glycerol concentrations. The CD spectra of TtABA-WDS and HvABA-WDS domains show that the two domains fold in the presence of glycerol (Figure 7). Increasing the glycerol concentration triggers an increase in the α-helical content of both domains. Upon increasing the glycerol concentration from 0 to 50% the α-helical content augments from about 8% to 21% for HvABA-WDS and from about 6% to 20% for TtABA-WDS (Figure 7). The CD spectra of TtABA-WDS display an isodichroic point at around 208 nm. In the case of HvABA-WDS, the spectra, with the notable exception of the spectrum at 50% glycerol, exhibit an isodichroic point at 204 nm. The peculiar behavior of HvABA-WDS observed at 50% glycerol might reflect the fact that additional, different structural transitions take place under these conditions suggesting that the glycerol-induced structural transition experienced by this domain is more complex than a two-state transition.  (6)). MRE: molar residue ellipticity (Θ) in deg cm 2 dmol −1 . (C) Transition diagrams of HvABA-WDS and TtABA-WDS. The α-helical content of the two domains, as derived from the CD spectra, was plotted as a function of the TFE concentration and data were fitted to sigmoidal curve, corresponding to a two-state transition. The inset shows the midpoint of transition and m values, as derived from the fitting, along with the quality of the fitting.

TtABA-WDS and HvABA-WDS Domains Gain α-Helicity under Conditions Mimicking Dehydration
Glycerol is known as a simplistic model for mimicking cellular dehydration that induces changes in the water surface tension around proteins. Here, we sought at assessing whether the two cereal domains could gain structure under increasing glycerol concentrations. The CD spectra of TtABA-WDS and HvABA-WDS domains show that the two domains fold in the presence of glycerol (Figure 7). Increasing the glycerol concentration triggers an increase in the α-helical content of both domains. Upon increasing the glycerol concentration from 0 to 50% the α-helical content augments from about 8% to 21% for HvABA-WDS and from about 6% to 20% for TtABA-WDS (Figure 7). The CD spectra of TtABA-WDS display an isodichroic point at around 208 nm. In the case of HvABA-WDS, the spectra, with the notable exception of the spectrum at 50% glycerol, exhibit an isodichroic point at 204 nm. The peculiar behavior of HvABA-WDS observed at 50% glycerol might reflect the fact that additional, different structural transitions take The gain of structure observed in the presence of glycerol could arise from dehydration or from crowding effects (i.e., increased excluded volume effects). To discriminate between these two hypotheses, we also recorded the far-UV CD spectra of both domains in the presence of 30% of sucrose. The addition of sucrose does not induce a pronounced folding in the two domains and triggers only a slight increase in the α-helical content of TtABA-WDS (Supplementary Figure S3). Upon addition of sucrose, the overall content in disorder decreases in TtABA-WDS domain, as observed from the decrease in the amplitude of the negative peak at 200 nm. However, the structural transition is not comparable to that observed with glycerol, supporting that these two domains probably fold because hydrogen bonds with water molecules are replaced by intramolecular ones and not because of excluded volume effects. place under these conditions suggesting that the glycerol-induced structural transition experienced by this domain is more complex than a two-state transition. The gain of structure observed in the presence of glycerol could arise from dehydration or from crowding effects (i.e., increased excluded volume effects). To discriminate between these two hypotheses, we also recorded the far-UV CD spectra of both domains in the presence of 30% of sucrose. The addition of sucrose does not induce a pronounced folding in the two domains and triggers only a slight increase in the α-helical content of TtABA-WDS (Supplementary Figure S3). Upon addition of sucrose, the overall content in disorder decreases in TtABA-WDS domain, as observed from the decrease in the amplitude of the negative peak at 200 nm. However, the structural transition is not comparable to that observed with glycerol, supporting that these two domains probably fold because hydrogen bonds with water molecules are replaced by intramolecular ones and not because of excluded volume effects.

TtABA-WDS and HvABA-WDS Confer Salt, Heat and Cold Stress Tolerance to Yeast Cells
The ASR gene family is absent in yeast cells [13]. In our previous work, we showed that overexpression of TtASR1 protein (cloned from Tunisian durum wheat landraces) increases Saccharomyces cerevisiae tolerance to salt and osmotic stress [9]. In a context

TtABA-WDS and HvABA-WDS Confer Salt, Heat and Cold Stress Tolerance to Yeast Cells
The ASR gene family is absent in yeast cells [13]. In our previous work, we showed that overexpression of TtASR1 protein (cloned from Tunisian durum wheat landraces) increases Saccharomyces cerevisiae tolerance to salt and osmotic stress [9]. In a context where no functional characterization has yet been carried out for the ABA-WDS domain, the TtABA-WDS and HvABA-WDS were expressed in the yeast Saccharomyces cerevisiae to assess their possible protective effects. The cDNAs encoding the two domains were sub-cloned into the pYES2 expression vector downstream of the inducible GAL1 promoter. Yeast cells were transformed with pYES2 empty vector or with the two recombinant vectors bearing the insert encoding either TtABA-WDS or HvABA-WDS. The recombinant yeast cells were subjected to salt (2M NaCl), cold (4 • C) and heat (50 • C) stress.
After transformation, growth assays of yeast cells on solid and liquid media in the presence of galactose and salt (2M NaCl) were conducted. Under salt stress the recombinant yeast strains expressing HvABA-WDS and TtABA-WDS domains grow better than the control strain bearing an empty vector, with HvABA-WDS conferring a higher resistance to salt stress compared to TtABA-WDS in both liquid and solid media ( Figure 8A,B).
After transformation, growth assays of yeast cells on solid and liquid media in the presence of galactose and salt (2M NaCl) were conducted. Under salt stress the recombinant yeast strains expressing HvABA-WDS and TtABA-WDS domains grow better than the control strain bearing an empty vector, with HvABA-WDS conferring a higher resistance to salt stress compared to TtABA-WDS in both liquid and solid media ( Figure  8A,B). In addition, expression of the two domains in yeast increases cell viability under extreme temperatures ( Figure 8C). The overexpression of the HvABA-WDS and TtABA-WDS domains leads to a 2.5-fold and 2.9-fold increase, respectively, in the number of CFUs under heat stress compared to the control strain bearing the empty vector. Likewise, the two domains exert a protective effect against cold stress. The overexpression of TtABA-WDS and HvABA-WDS confers a 1.3-fold and 2.4-fold increase, respectively, in cell viability compared to the control strain ( Figure 8C).
In conclusion, yeast cells transformed with plasmids driving the expression of TtABA-WDS and HvABA-WDS are less sensitive to abiotic stress (salt, cold and heat) In addition, expression of the two domains in yeast increases cell viability under extreme temperatures ( Figure 8C). The overexpression of the HvABA-WDS and TtABA-WDS domains leads to a 2.5-fold and 2.9-fold increase, respectively, in the number of CFUs under heat stress compared to the control strain bearing the empty vector. Likewise, the two domains exert a protective effect against cold stress. The overexpression of TtABA-WDS and HvABA-WDS confers a 1.3-fold and 2.4-fold increase, respectively, in cell viability compared to the control strain ( Figure 8C).
In conclusion, yeast cells transformed with plasmids driving the expression of TtABA-WDS and HvABA-WDS are less sensitive to abiotic stress (salt, cold and heat) compared to cells transformed with an empty vector, indicating that the ABA-WDS domains are expressed and functional in yeast cells and confer them abiotic stress tolerance.

HvABA-WDS and TtABA-WDS Domains Stabilize the Enzyme LDH against Thermal and Dessication Stress
To ascertain whether the HvABA-WDS and TtABA-WDS domains have a chaperone activity in vitro, we assessed their ability to protect LDH activity under heat and dehydration stress, as already described for ASR proteins from tomato [48], banana [13], lily [49] and durum wheat [9]. The protective effects exerted by the two domains were compared to those observed either in the presence of BSA, as an example of a non-specific protectant, or in the absence of any protein. In the presence of BSA, LDH gradually loses its activity, with its residual activity being about 11% after 60 min at 43 • C ( Figure 9A). By contrast, heat inactivation of LDH is significantly reduced in the presence of HvABA-WDS and TtABA-WDS with 33% and 52% of its initial activity being retained, respectively, after 60 min at 43 • C ( Figure 9A). After dehydration and rehydration, the LDH activity was reduced to 28% of its initial value. TtABA-WDS and HvABA-WDS domains were found to confer a protective effect against dessication-induced enzyme inactivation, with the LDH residual activity reaching values of~70% and 60%, respectively ( Figure 9B). Results indicate that TtABA-WDS provides a higher degree of protection against thermal and dessication-induced denaturation compared to HvABA-WDS.
[49] and durum wheat [9]. The protective effects exerted by the two domains were com-pared to those observed either in the presence of BSA, as an example of a non-specific protectant, or in the absence of any protein. In the presence of BSA, LDH gradually loses its activity, with its residual activity being about 11% after 60 min at 43 °C ( Figure 9A). By contrast, heat inactivation of LDH is significantly reduced in the presence of HvABA-WDS and TtABA-WDS with 33% and 52% of its initial activity being retained, respectively, after 60 min at 43 °C ( Figure 9A). After dehydration and rehydration, the LDH activity was reduced to 28% of its initial value. TtABA-WDS and HvABA-WDS domains were found to confer a protective effect against dessication-induced enzyme inactivation, with the LDH residual activity reaching values of ~70% and 60%, respectively ( Figure 9B). Results indicate that TtABA-WDS provides a higher degree of protection against thermal and dessication-induced denaturation compared to HvABA-WDS.

Discussion
Experimental information on the conserved ABA-WDS domain of the ABA-WDS family (Pfam 02496) in isolation is lacking. To fill this gap in knowledge, we herein isolated from stressed leaves the cDNAs encoding two ABA-WDS domains, one from durum wheat and one from barley. Using various bioinformatics approaches, we showed that these domains are predicted to be intrinsically disordered, and further confirmed their disordered nature experimentally. Analysis of their deduced amino acid sequence revealed a biased composition, i.e., an enrichment in disorder-promoting amino acids (A, H, and E) as reported for most ASR [17] and LEA proteins [21]. Many Pfam domains contain sequences with 100% predicted intrinsic disorder [50]. Noteworthy, the Pfam database recently added "disordered" as a specific entry type, underscoring the im- with respect to the LDH activity recorded before stress treatment. The values are the mean ± SD from three independent samples. The asterisks denote statistically significant differences (** p < 0.01, *** p < 0.001) as evaluated with Student's T test.

Discussion
Experimental information on the conserved ABA-WDS domain of the ABA-WDS family (Pfam 02496) in isolation is lacking. To fill this gap in knowledge, we herein isolated from stressed leaves the cDNAs encoding two ABA-WDS domains, one from durum wheat and one from barley. Using various bioinformatics approaches, we showed that these domains are predicted to be intrinsically disordered, and further confirmed their disordered nature experimentally. Analysis of their deduced amino acid sequence revealed a biased composition, i.e., an enrichment in disorder-promoting amino acids (A, H, and E) as reported for most ASR [17] and LEA proteins [21]. Many Pfam domains contain sequences with 100% predicted intrinsic disorder [50]. Noteworthy, the Pfam database recently added "disordered" as a specific entry type, underscoring the importance of disorder in the definition of protein domains [51]. Illustrative examples of proteins with IDDs playing key regulatory roles in controlling the cell cycle are p27kip1, p21Waf1/Cip1/Sdi, p57kip2 [52,53] and p53 [54]. Many of the signaling pathways and regulatory systems in eukaryotic cells are controlled by proteins bearing IDDs that mediate specific protein-protein and protein-nucleic acid interactions. As an example, TFs contain a high fraction of functionally essential intrinsically disordered regions (IDRs) [55]. In this context, direct and indirect evidence that ASR proteins can function as TFs [6,8] does not come as a surprise.
Although the two ABA-WDS domains are predicted to be disordered, they slightly differ in their predicted disorder content, with HvABA-WDS being predicted to be slightly more ordered than its cognate TtABA-WDS domain. This subtle predicted difference was experimentally confirmed by spectroscopic (i.e., CD) and hydrodynamic (i.e., SEC) studies.
The ellipticity values at 200 and 222 nm, as observed in the far-UV CD spectra, enabled classifying TtABA-WDS and HvABA-WDS within the RC-like and PMG-like subfamilies, respectively. Slight differences between the two domains are also observed in their near-UV CD spectra. In line with spectroscopic data, the CI values, calculated from SEC, spotlight conformational differences between the two domains and showed that HvABA-WDS is more compact than TtABA-WDS. This slight difference in disorder content could result from subtle sequence differences (see Figure 1A). Although the two domains exhibit a high sequence identity (~94%), the presence of a valine in HvABA-WDS that replaces Ala17 in TtABA-WDS entails the formation of a bigger hydrophobic cluster in HvABA-WDS (Figure 2A) that may explain the comparatively higher order-propensity of the latter. It should be pointed out that conformational differences have the potential to also impact function, as well exemplified for instance by ubiquitin E3 ligases whose ubiquitination activity was shown to be reduced by amino acid substitutions inducing a more extended shape [56].
The differences in compactness and in the extent of secondary structure between the two domains mirror differences in their folding abilities. Indeed, while TtABA-WDS undergoes temperature-induced folding, as already observed for other IDPs [57,58], HvABA-WDS does not. The lack of effect of temperature on the secondary structure content of HvABA-WDS might arise from the fact that this domain is already more structured compared to TtABA-WDS. Interestingly, this differential behavior of the two domains vis-à-vis temperature mirrors to some extent the behavior we previously observed for the ASR proteins from durum wheat (TtASR1) and from barley (HvASR1), where the latter was found to undergo a temperature-induced folding less pronounced than that of the cognate ASR protein from durum wheat in spite of the high sequence similarity between the two proteins [17].
Many studies showed that various mechanisms trigger a gain of the structure by IDPs/IDDs, including disulfide bond formation, ion coordination, macromolecular partner binding and many others [59]. Taking into account the role of ABA-WDS domains in plant abiotic stress tolerance, it is conceivable that environmental changes may induce their folding. In line with this, in a previous work, we showed that TtASR1 and HvASR1 gain structure upon the addition of TFE, glycerol and zinc ions [17]. In the presence of TFE, a cosolvent commonly used in folding studies, both HvABA-WDS and TtABA-WDS domains gain α-helical structure. The two domains also undergo a disorder to α-helix transition upon glycerol-induced water deficit stress. Notably, while the CD spectra obtained in presence of TFE show an isodichroic point reflecting a two-state transition ( Figure 6), a more complex scenario is observed in the presence of glycerol. In this case indeed, while the CD spectra of TtABA-WDS display an isodichroic point, in the case of HvABA-WDS only the CD spectra obtained in the 10-40% glycerol range do exhibit an isodichroic point ( Figure 7) suggesting that distinct structural transitions could take place at 50% and that the glycerol-induced structural transition experienced by this domain is more complex than a simple two-state transition.
The subtle differences in terms of disorder content and folding propensity between the two domains parallel those previously observed with the ASR proteins from durum wheat (TtASR1) and barley (HvASR1) [17]. These ASR proteins, which encompass an ABA-WDS domain distinct, though similar to the domain herein reported (see Supplementary Figure  S4), were indeed shown to feature slight conformational differences, with HvASR1 being more ordered than TtASR1 [17].
The fact that even as few amino acid differences as those observed between the two ABA-WDS domains engender measurable differences in their conformational properties and folding abilities, illustrates how the amino acid sequence finely tune protein behavior and properties. In addition, it should be emphasized that the conformational differences between the two domains could be hardly anticipated in light of the conservative nature of the amino acid substitutions that they bear, which further underscores the value of experimental approaches.
The abundance of histidine residues in the HvABA-WDS (12.1%) and TtABA-WDS (13.6%) domains provides a hint of their possible ability to bind bivalent cations such as zinc. We previously showed that the HvASR1 and TtASR1 proteins fold upon binding to zinc and mapped the region that undergoes Zn-induced structural transition [17]. This region contains a motif (PEHAHKHK) that was shown to be responsible for binding to zinc in tomato ASR1 [60]. This motif is also conserved in the ABA-WDS domains herein described (see Figure 1A and Supplementary Figure S4), thus further supporting their expected ability to bind zinc. The metal binding abilities of these domains and their impact on their conformational properties are the focus of ongoing studies that, ultimately, are expected to enable establishing correlations with their functions and possible transactivation capabilities.
Although the precise physiological function of the ASR gene family remains elusive, ASR proteins have been shown to have chaperone-like activity for direct plant protection under stress [3,7,9]. We recently reported that the ASR protein durum wheat (TtASR1) can function as a chaperone-like protein and increases Sacchaomyces cerevisiae tolerance under salt and osmotic stress [9]. In light of the fact that no functional characterization has yet been carried out for the ABA-WDS domain in isolation, we resorted to investigating the functional impact of TtABA-WDS and HvABA-WDS domains both in vitro and in vivo. S. cerevisiae is a widely used model for the functional characterization of stress-associated proteins. Here, we show that the heterologous expression of TtABA-WDS and HvABA-WDS domains increase S. cerevisiae viability under cold, heat and salt stress. The yeast growth assays carried out under stress conditions show that HvABA-WDS has a more pronounced protective effect compared to TtABA-WDS. Although these functional differences may be tied to the subtle differences observed in the structural and folding properties of these two domains, the possibility that the lower efficiency of TtABA-WDS in conferring tolerance to abiotic stress to yeast might arise from a lower expression level compared to HvABA-WDS cannot be ruled out. Definitive conclusions await future studies that are currently in progress in our lab.
LDH activity assays confirmed and extended the results obtained in vivo and showed that the TtABA-WDS and HvABA domains are able to protect the sensitive LDH enzyme from heat and dehydration in vitro. While HvABA-WDS was found to be more effective in protecting yeast against abiotic stress than TtABA-WDS, the opposite scenario is observed in vitro. The discrepancy between the in vitro and the in vivo results may arise from the expectedly higher complexity of chaperone activity in vivo, where the function of chaperones can potentially rely on various co-chaperones (e.g., proteins, membrane, DNA . . . ). Moreover, ASR proteins were shown to function as transcriptional factors in addition to their chaperone activity. The in vivo abiotic stress tolerance phenotype is the result of all these combined activities.
The finding that the ABA-WDS domains have a chaperone activity is in line with the well-established relationship between intrinsic disorder and chaperone function [61]. Many chaperones are either fully disordered or contain IDRs that are involved in the regulation of the chaperone or in the interaction with the substrate itself [62]. Interactions mediated by IDRs of chaperones allow them to assist the folding of a much broader range of substrates. In other words, the lack of structure of chaperones is instrumental for their chaperone activity. In support of this, the chaperone activity of tomato ASR1 is not enhanced by the gain of structure induced by Zn 2+ [19]. The conformational plasticity of the ABA-WDS domain could explain the multiple physiological functions attributed to the proteins of the ASR family. Additional experiments will be needed to unravel the relationships between the conformation adopted by the ABA-WDS domains under various stress conditions and their diverse biological functions.
Although awaiting validation in planta, the in vivo assays in yeast and the in vitro LDH protective assay provide evidence for the implication of the ABA-WDS domain in the known function attributed to the ASR family and advocate for the potential use of the ABA-WDS domain in biotechnological applications.
Finally, the present study contributes to shedding light on the hitherto ambiguous classification of ASR proteins. Their biased amino acid composition constitutes one of the pillars onto which relies on the previously proposed classification of the ASR protein family within the LEA superfamily [23,24]. In light of the lack of sequence similarity with any of the recognized LEA group proteins, Hunault and Jaspard proposed however the exclusion of ASR proteins from the LEA family [26]. While waiting for more extensive studies, the results herein presented would advocate in favor of a reconsideration of the ABA-WDS family as a member of the LEA superfamily. In fact, the structural and functional characterization of the wheat and barley domains unveiled that beyond a biased composition, they share most of the known features of LEA proteins, i.e., (i) they are intrinsically disordered, (ii) they undergo α-helical folding under desiccation, and (iii) they are endowed with a heat, desiccation and cold protective activity and their overexpression in yeast confer abiotic stress tolerance.

Plant Materials and Stress Treatments
The local Tunisian genotype of tetraploid wheat T. turgidum L. subsp. durum (2n = 4x = 28) cv. Mahmoudi and Hordeum vulgare L. (Rihane), were used for ABA_WDS cDNA cloning. All seeds were initially surface sterilized by a 0.5% NaOCl wash for 15 min, rinsed three times with sterile water, and germinated on wet Whatman paper filter placed in Petri dishes after 2 days in the dark. Ten-day-old seedlings grown were subjected to stress. For salinity seedlings were incubated in 200 mM NaCl. Stressed seedlings were sampled at 24 h of treatment. Leaves were harvested and frozen in liquid nitrogen for RNA isolation.

RNA Extraction and Amplification of TtABA-WDS and HvABA-WDS cDNA Domains from Durum Wheat and Barely
RNA was isolated from approximately 200 mg of durum wheat and barley leaves according to the Trizol method (Invitrogen) and following the manufacturer's instructions. RNA was treated with RNase-free DNase to remove any contaminating DNA. Reverse transcription reactions were performed for 1 h at 37 • C using MML-reverse transcriptase (Invitrogen) and oligo-dT. First-strand complementary DNA (cDNA) was used as a template for PCR amplifications using a couple of specific primers (PD1-PD2) designed from the alignment of a set of full ASR protein sequences belonging to monocot plant species [9]. Control amplifications in the absence of the reverse transcriptase were also performed to rule out any amplification caused by the presence of contaminating DNA. The TtABA-WDS and HvABA-WDS cDNA were cloned into the EcoRI site of the E. coli expression vector pGEX-4T-1, resulting in a fusion with glutathione S-transferase (GST). The two domains were expressed and purified as previously described [17]. Removal of the GST tag by thrombin results in a recombinant product bearing a vector-encoded N-terminal GSPEF and a C-terminal EFPGRLERPHRD amino acid extension. Disorder content analysis of the resulting recombinant products (Supplementary Figure S1C,D) indicates that the presence of the vector-encoded residues preserves (actually slightly increases) the predicted disordered nature of the domains, in line with the expectedly disordered nature of vector-encoded linkers.
Protein concentrations were calculated using the theoretical absorption coefficients at 280 nm as obtained using the program ProtParam at the EXPASY server (https://web. expasy.org/protparam/ (accessed on 25 January 2021)).

Intact Protein Mass Analysis
Protein masses were determined on purified solution samples. An amount of 1 µL of protein at~20 µM was mixed with 1 µL of α-Cyano-4-hydroxycinnamic acid matrix (Bruker Daltonics, Wissembourg, France) solution in 0.3% TFA/CH3CN (50:50 v/v). One µL of the mix was spotted on the target and analysed by MALDI-TOF on an Ultraflex III spectrometer (Bruker Daltonics, Wissembourg, France) controlled by the Flexcontrol 3.0 package (Build 51) and operated in the linear mode, using a maximum accelerating potential of 25 kV and a 5000-15,000 m/z range (LP_Protmix_Method). The laser frequency was fixed to 100 Hz and~1000 shots per sample were cumulated. Four external standards (Protein Calibration Standard I, Bruker Daltonics, Wissembourg, France) were used to calibrate each spectrum to a mass accuracy within 100 ppm. Peak picking was performed using the FlexAnalysis 3.0 software with an adapted analysis method. Parameters used were: centroid peak detection algorithm, S/N threshold fixed to 5 and a quality factor threshold of 30.
For peptide ionization in the EASY-Spray source, the spray voltage was set at 1.9 kV and the capillary temperature at 275 • C. The spectrometer was operated in the datadependent mode to switch consistently between MS and MS/MS. MS spectra were acquired in the 300-1700 m/z range at full width at half maximum (FWHM) resolution of 30,000 measured at 400 m/z. For internal mass calibration, the 445.12 m/z ion was used as lock mass. The 10 most abundant precursor ions were selected and collision-induced dissociation fragmentation was performed in the ion trap to insure both maximum sensitivity and a maximum amount of MS/MS data. The signal threshold for an MS/MS event was set to 500 counts. Charge state screening was enabled to exclude precursors with 0 and 1 charge states. Dynamic exclusion was enabled with a repeat count of 1, exclusion list size of 500 MS runs and exclusion duration of 30 s. The acquired raw data were processed using Proteome Discoverer (version 1.4.1.14, Thermo Fisher Scientific). Spectra were searched using SEQUEST (Thermo Fisher Scientific) against a homemade database comprising 20,150 human and 4306 E. coli sequences and the TtABA-WDS and the HvABA-WDS sequences. To detect contaminants, an additional search against the Swiss-Prot database (https://www.uniprot.org/uniprot (accessed on 25 January 2021); version 2014.02; 20,284 entries) supplemented with a set of 245 frequently observed contaminants (e.g., human keratin) was performed using MASCOT (http://www.matrixscience.com (accessed on 25 January 2021)). Search parameters were: (i) GluC; five miscleavage allowed; (ii) mass tolerance of 6 ppm for monoisotopic precursor ions and 0.8 ppm for fragment ions from MS/MS; (iii) Cys carbamidomethylation (+57.02146 Da) as a fixed modification; Met oxidation (+15.99491 Da) and N-terminal acetylation (+42.0106 Da) as variable modifications; (iv) minimum peptide length of four residues. Only high-score peptides were selected. Proteins were identified with a false discovery rate (FDR) of 1%.

Size Exclusion Chromatography and Calculation of Hydrodynamic Radii
The hydrodynamic radii (Stokes radii, R S ) of the TtABA-WDS and HvABA-WDS domains were estimated by analytical SEC as previously described [17]. Typically 2 mg mL-1 of purified protein was injected. The SEC buffer was PBS pH 7.3.
The Stokes radii of domains eluted from the SEC column were deduced from a calibration curve obtained using globular proteins of known molecular mass (MM, in Daltons) and whose RS (in Å) was calculated according to [63]: The R S (in Å) of a natively folded protein (R s N F ) and of a fully unfolded state in urea (R S U) with a molecular mass (MM) (in Daltons) were calculated according to [64]: The R S of an IDP with N residues was also calculated according to [65] using the simple power-law model: where R 0 = 2.49 and ν = 0.509. The compaction index (CI) is expressed as according to [66]: This parameter, which allows comparison between proteins of different lengths, increases with increasing compaction and, in principle, varies between 0 and 1.

Differential Scanning Fluorimetry (DSF)
The DSF analysis of TtABA-WDS and HvABA-WDS domains was carried out in the presence of a fluorescent dye as already described [17]. This experiment was conducted using a PCR instrument (Biorad) and 96-well plates containing 25 µL of mixture per well. Each well contained 21.5 µL of TtABA-WDS domain or HvABA-WDS domain at 1 mg mL −1 in solution and 3.5 µL of SYPRO Orange solution (from a 5000× stock solution). Fluorescent signals were acquired with excitation and emission wavelengths at 485 nm and 625 nm, respectively. Temperature scans were performed from 20 • C to 90 • C.

Circular Dichroism (CD) Measurements
The CD spectra of TtABA-WDS and HvABA-WDS were recorded using a Jasco 810 dichrograph, flushed with N2 and equipped with a Peltier thermoregulation system [17]. One-mm or 1-cm thick quartz cuvettes were used for far-and near-UV CD measurements, respectively. Proteins concentrations were 0.12 mg mL −1 and 1.2 mg mL −1 for far-and near-UV CD studies, respectively. Far-UV CD spectra were measured between 190 and 260 nm, while near-UV CD spectra were recorded between 250 and 350 nm. Unless differently specified, CD spectra were recorded in 10 mM sodium phosphate pH 7 at 20 • C. The scanning speed was 20 nm/min, with data pitch of 0.2 nm. Each spectrum is the average of three acquisitions. The spectrum of buffer was subtracted from the protein spectrum. Spectra were smoothed using the "means-movement" smoothing procedure implemented in the Spectra Manager package.
Mean molar ellipticity values per residue (MRE) were calculated as (Θ) = 3300 M∆A/(lcn), where l is the pathlength in cm, n is the number of residues, M is the molecular mass in Daltons and c is the concentration of the protein in mg mL −1 . Numbers of amino acid residues are 83 for both domains. The molecular mass is 9200 Da and 9234 Da for TtABA-WDS and HvABA-WDS, respectively.
Temperature-induced protein folding events were assessed by measuring the ellipticity at a fixed wavelength of 230 nm in the temperature range of 20-80 • C, with data pitch 2 • C and a temperature slope of 2 • C/min and protein concentrations of 0.12 mg mL −1 . The reversibility of the process was assessed by monitoring the ellipticity during cooling from 80 • C to 20 • C.
The DICHROWEB website (http://dichroweb.cryst.bbk.ac.uk/html/home.shtml (accessed on 25 January 2021)), which was supported by grants to the BBSRC Centre for Protein and Membrane Structure and Dynamics (CPMSD) [67], was used to analyze the experimental data in the 190-260 nm range. The content in the various types of secondary structure was estimated using the CONTIN deconvolution method with the reference protein set 7 [68].
For the estimation of the percentage of residues adopting an α-helical conformation in the presence of various additives (TFE, glycerol, and sucrose) we analyzed the value of molar ellipticity at 220 nm observed under each condition and divided it by the value expected for a protein in which all residues adopt an α-helical conformation (i.e., 100% α-helix). The latter was calculated according to the following empirical relationship [69]: where n is the number of residues.

Salt Tolerance and Growth Assays of Yeast Cells
The effect of salt stress on the growth of yeast strain W303-1B (MATα; ade2-1; ura3-1; his3-11,15; leu2-3112; trp1-1; can1-100) of S. cerevisiae harboring recombinants plasmids (TtABA-WDS or HvABA-WDS) or empty vector (pYES2) was studied. The pYES2 expression vector (Invitrogen) is a 2 µ based multicopy plasmid and contains the URA3 gene and the Gal1 promoter for selection and expression in yeast. Stationary yeast precultures grown on glucose-containing minimal medium (MM) were diluted to 5 × 106 cells mL −1 and used to inoculate (1:100) yeast cultures on MM Gal in the absence or presence of NaCl (2M) under continuous agitation (200 r.p.m.) at 30 • C. The OD at 600 nm was measured every 2 h. For kinetic studies on solid medium, serial dilutions of different recombinant yeast clones were prepared (10 −1 -10 −3 ) from the same precultures, and then 7 µL from each dilution were spotted onto solid MMGal containing NaCl (2M) and incubated at 30 • C for three days. Cell viability was estimated by counting the number of colony-forming units after incubation of the plate overnight at 30 • C. The viability ratio under salt stress was calculated according to the following formula: (mean of colony number on a stressed plate/mean of colony number on a control plate) × 100%.

LDH Protective Assay
The activity of lactate dehydrogenase (LDH, EC1.1.1.27, rabbit muscle lactate dehydrogenase) was measured as previously described [9]. One µL of a solution of LDH at 20 µg/mL (Sigma, Tokyo, Japan) was added to a 20 µL-solution containing 20 µg/mL of either bovine serum albumin (BSA), or TtABA-WDS or HvABA-WDS in 10 mM sodium phosphate, pH 7.4. Samples were then submitted to heat stress treatment. To test the effect of desiccation on LDH activity with or without TtABA-WDS or HvABA-WDS, samples were vacuum-dried in a Speed Vac (Concentrator plus / Vacufuge ® plus of Eppendorf AG, Hamburg, Germany) to a final volume of 6 µL and rehydrated by the addition of 14 µL of water.
To determine the LDH activity, 20 µL of the LDH mixture was added to 980 µL of assay buffer (10 mM sodium phosphate, pH 7.4, 2 mM NADH, and 10 mM pyruvic acid). NADH oxidation was monitored by measuring the absorbance at 340 nm over 3 min, during which the reaction rate was linear. The rate of absorbance was then used to calculate the activity ∆A340/min × 8095 = U/L (Biomaghreb kit, Biomaghreb, Tunis, Tunisia). All samples were assayed in triplicate.

Bioinformatics Analyses
Identification of the conserved domain of TtABA-WDS and HvABA-WDS was carried out by query against the Conserved Domain Database (CDD) (http://www.ncbi.nlm. nih.gov/Structure/cdd/wrpsb.cgi (accessed on 25 January 2021)). The ExPASy server (http://web.expasy.org/protparam/ (accessed on 25 January 2021)) was used to derive the physicochemical properties of the domains. Sequence alignment of the TtABA-WDS and HvABA-WDS domains was performed using Clustal Omega (https://www.ebi.ac. uk/Tools/msa/clustalo/ (accessed on 25 January 2021)). The RH-plot was generated by the PONDR server (http://www.pondr.com (accessed on 25 January 2021)). The phase diagram plots of TtABA-WDS and HvABA-WDS were generated using the CIDER server (http://pappulab.wustl.edu/CIDER/analysis/ (accessed on 25 January 2021)). Deviations in amino acid composition of TtABA-WDS and HvABA-WDS were computed as previously described [17] using the average amino acid frequencies of the SWISS-PROT database (as obtained from https://web.expasy.org/protscale/pscale/A.A.Swiss-Prot.html (accessed on 25 January 2021)) as the reference value. The average amino acid frequencies of the SWISS-PROT database roughly correspond to the mean composition of proteins in nature. If the average composition of an amino acid X in SWISS-PROT proteins is CSPX, and CPX is the composition of X within a protein P, deviation from the composition of X of SWISS-PROT proteins was defined for P as (CPX-CSPX)/CSPX. The HCA plots and the IUPRED disorder prediction were obtained using the MeDor server [34]. Structure prediction was carried out using the Phyre2 protein fold recognition server (http://www.sbg.bio.ic.ac.uk/ phyre2/html/page.cgi?id=index (accessed on 25 January 2021)) which uses advanced remote homology detection methods to build 3D models [37].

Conclusions
In this work, two ABA-WDS domains of the ASR family were isolated from durum wheat and barley. These two domains were shown to be intrinsically disordered and to be able to undergo induced folding in the presence of TFE and glycerol, which mimic interaction and dehydration respectively. Strikingly, both domains confer stress resistance to yeast cells, as well as heat and dessication resistance to LDH in vitro. To the best of our knowledge, this is the first time a molecular function as a chaperone can be experimentally attributed to an ABA-WDS domain. Our results reveal that the ABA-WDS domain is a bona fide domain, capable of exerting its function in isolation, and whose conformational plasticity could underlie the functional versatility of the ASR protein family. The present results provide a foundation for further evolutionary and functional characterization of the ASR gene family, and contribute to enlarge our knowledge of this protein family. Beyond a fundamental interest, the present study also holds promise in the field of plant improvement. After an in planta validation step (which is currently in progress in our lab), the sequence encoding the ABA-WDS could be transferred to various crops to improve their tolerance to abiotic stress.
Lastly, and from a more general perspective, the data herein presented provide evidence that the ABA-WDS domain could be used to preserve enzyme activities under heat and cold stress thus paving the way towards possible biotechnological applications in the field of enzymology.