Global Implications of Local Unfolding Phenomena, Probed by Cysteine Reactivity in Human Frataxin

Local events that affect specific regions of proteins are of utmost relevance for stability and function. The aim of this study is to quantitatively assess the importance of locally-focused dynamics by means of a simple chemical modification procedure. Taking human Frataxin as a working model, we investigated local fluctuations of the C-terminal region (the last 16 residues of the protein) by means of three L → C replacement mutants: L98C, L200C and L203C. The conformation and thermodynamic stability of each variant was assessed. All the variants exhibited native features and high stabilities: 9.1 (wild type), 8.1 (L198C), 7.0 (L200C) and 10.0 kcal mol−1 (L203C). In addition, kinetic rates of Cys chemical modification by DTNB and DTDPy were measured, conformational dynamics data were extracted and free energy for the local unfolding of the C-terminal region was estimated. The analysis of these results indicates that the conformation of the C-terminal region fluctuates with partial independence from global unfolding events. Additionally, numerical fittings of the kinetic model of the process suggest that the local transition occurs in the seconds to minutes timescale. In fact, standard free energy differences for local unfolding were found to be significantly lower than those of the global unfolding reaction, showing that chemical modification results may not be explained in terms of the global unfolding reaction alone. These results provide unequivocal experimental evidence of local phenomena with global effects and contribute to understanding how global and local stability are linked to protein dynamics.

Even though proteins tend to fold into compact conformations of roughly equal stability that minimize their energy and maximize the system's entropy, they are not rigid, but flexible objects with dynamics that allow them to interact with the surroundings and accomplish their biological function. Several weak intramolecular interactions that stabilize the native conformation exhibit half-lives in the order of the picoseconds to the nanoseconds, implying that those bonds are constantly being formed and broken, giving place to the faster events of protein dynamics 1 . Other common events take place in longer time scales and might have local or global effects over the molecule. Those events, which are a consequence of conformational flexibility, taken together give place to the so-called protein breathing.
Such fluctuations are characterized by the temporal time-scale in which they occur (kinetical component) and by the magnitude and directionality of the change taking place (structural component). It is worth noting that dynamics are closely dependent on the system's conditions (temperature, pH, ionic strength, buffer components, intermolecular protein−protein interactions, etc.), and that modifying those conditions will result in the alteration of relative populations of states and their interconversion kinetics. The perturbation of global stability, as a consequence of the alteration of the local stability of a region, might result in a change in a molecule's internal motions. The prediction of the amplitude or the time-scale of the movements that characterize this change is not obvious if the analysis is performed over the 3D structure of the protein alone. To meaningfully address this issue, molecule is unfolded; in the latter, only a region loses its structure-the CTR in our naive model-while the rest of the molecule retains its native fold. Thiol modification with different reactants enables the study of the mobility (conformational exchange) of diverse positions and the determination of the contribution to stability of the interaction of the CTR with the rest of the molecule. In this way, by applying different stability-modulating conditions, we gained a detailed understanding of the open-close mechanism, strictly understood as the process by means of which each position fluctuates between an exposed (reactive) and a hidden (non-reactive) conformation.
We found that, despite the lack of periodic structure, the modification rate significantly differs for Cys residues located in different positions along the CTR and the reactivity is significantly lower than the expected for free thiols. Besides, we described the effect over local dynamics of stabilizing or denaturing agents. Based on our results, we propose a reaction mechanism consistent with the existence of at least two native substates.

Materials and Methods
Please note that while as much methods description as is necessary to follow the logic of the text is detailed below, the bulk of technical description of materials and methods is included in Supplementary materials. Cysteine modification kinetics. In order to assess the dynamics of specific positions of FXN, we determined proteins' reactivity towards 5,5′-dithiobis-2-nitrobenzoic acid (DTNB)-also known as Ellman's reagentand 4,4′-dithiodipyridine (DTDPy). β-mercaptoethanol (BME) was used to quantify DTNB and DTDPy stock solutions. In order to find the intrinsic modification rate coefficient (k mod ), we decided to use short peptides corresponding to the stretch comprised between T196 and S206, in which each Leu residue, one at a time, was replaced by a Cys residue. Peptides pL198C: TKCDLSSLAYS, pL200C: TKLDCSSLAYS and pL203C: TKLDLSSCAYS were synthesized (GenScript, Piscataway, NJ, USA) and their labeling rates were measured at varying concentrations of DTNB and DTDPy. Proteins or peptides were used at a 30 µM concentration, in a final volume of 500 µl. Solutions were prepared in 20 mM Tris·HCl, 100 mM NaCl, pH 7.0, with varying concentrations of DTNB or DTDPy (0-3.0 mM). Once the probe was added, the solution was gently homogenized and transferred to a cuvette to follow the absorbance change. Time courses were carried out in a Jasco V550 UV/VIS spectrophotometer with a thermostated cell holder connected to a circulating water bath set at 25 °C. For DTNB modification, absorbance was measured at 425 nm, with band-pass of 4 nm, using a 1.0 cm path length cell. For DTDPy modification, absorbance was measured at 324 nm, with band-pass of 4 nm, using a 0.3 cm path length cell. For each kinetic trace, between 1000 and 3000 data points were obtained. When reactions occured in the timescale of seconds, an RX2000 stopped flow instrument (Applied Photophysics Ltd, UK) was used.
Analysis of local dynamics. In order to describe the local unfolding process, we considered a model that assumes that for a given position, the open (reactive) and closed (non-reactive) conformations are at equilibrium 14,15 . In such a condition, the concentration of species depends on an open microscopic rate coefficient (k open ) and a close microscopic rate coefficient (k close ). In turn, the open conformation may react with an appropriate probe (P), with an intrinsic modification rate coefficient (k mod ) that describes the reaction rate of the thiolate with the reactant. The value of k mod is the product of the probe's concentration and the bimolecular rate constant for the reaction of an unprotected free thiol (k mod = k [P]), such as small Cys-containing peptides, BME, cysteine or glutathione. The model can be summarized in the following way: To rule out in which regime the reaction occurs in a certain experimental condition, it suffices to find the dependence of the value of k label with the concentration of the reactant. If conditions determine an EX1 regime, the value of k label will be constant given that it only depends on the value of k open , which is independent of the reactant concentration. In the case that the regime is of the EX2 type, the magnitude of k label will depend on the value of k mod , which varies linearly with the reactant concentration. Given that in EX2 k label depends on the values of  K C O and k mod -only dependent on the chemical nature of the reactant and the thiolate-it is possible to determine the local unfolding free energy just by finding the value of the apparent rate coefficient for the labeling reaction. Under steady state conditions k label may be written in terms of the microscopic rate constants 17 : On the other hand, a reaction's equilibrium free energy can be expressed in terms of the constant that defines that equilibrium: The combination of Eqs 5 and 6 gives a relation between the experimentally-measured (apparent) rate coefficient (k label ), the intrinsic microscopic modification coefficient (k mod ) and Δ

Results
Design of CTR mutants. The alteration of the C-terminal stretch of different FXN homologues has been shown to affect the protein's global stability 4 . Besides, the effect of point mutations in the CTR and of the complete deletion of this region over the stability, dynamics and folding have been previously studied 5,6,18 . Here we focus on residues L198, L200 and L203, which form a cluster of hydrophobic interactions among themselves and with residues from the core of the protein (Fig. 1A). Variant FXN90-210, nine residues shorter in its N-terminal region than the mature form FXN81-210 (currently considered the functional form in vivo), was the wild-type control sample for conformation and stability, and the base variant over which all mutants were constructed. Note that residues 81-89 have been shown to be disordered 19 . We prepared point-mutants FXN L198C, FXN L200C and FXN L203C, which locate a thiol in different positions of the CTR. Those constructions allowed us to investigate the molecular motions by assessing thiol-exchange reactivity. A fourth mutant, FXN V134C, with its Cys completely protected in the folded ensemble served as a hidden-Cys control.
The four variants are soluble and monomeric, the same as wild-type FXN. Far and near-UV circular dichroism (CD) spectra showed that both the secondary structure content and the tertiary packing are native-like (Fig. S2A,B). Besides, variants' resistance to proteolysis probed their rigidity (Fig. S3). Moreover, wild-type and Cys point mutants are able to activate desulfurase NFS1 alike, indicating that this function remains intact upon mutation (Fig. S4A). On the other hand, iron binding capability was slightly altered, especially in the case of L200C and in a minor degree in the case of L198C (Fig. S4B).
Solvent accessible surface area and pK a predictions performed over computational models based on the X-ray structure of FXN indicated that the mutated residue in each mutant may expose less than 10% of its surface ( Fig. S1 and Table S1) and suggested that cysteines would be protonated in a broad range of pH values, particularly in our working conditions. This means that cysteines would mainly behave as apolar residues, capable of establishing most of the hydrophobic interactions that a leucine would establish in the context of the wild-type protein ( Fig. S1 and Table S1). However, given the dual polar-apolar nature of cysteine residue, the existence of polar interactions cannot be ruled-out.

Thermodynamic stability of FXN variants.
A severe thermodynamic destabilization would preclude us from applying the proposed experimental strategy. For instance, if mutations give rise to very altered equilibrium unfolding constants, a partially folded or unfolded state might be so populated that the labeling reaction could significantly proceed from these states rather that from a substate of N. To investigate whether Cys point mutations affect the native stability of FXN90-210, chemical and temperature equilibrium unfolding experiments were performed (Fig. 2). Although folding kinetic experiments showed that FXN mainly folds via an intermediate state 18 , the shape and overlapping of the unfolding profiles followed by CD and tryptophan fluorescence revealed that all mutants studied in this paper exhibit unfolding curves compatible with a two-state model (  N U). This fact was previously described for equilibrium unfolding of several full-length FXN variants 5,6,20 , and it is consequence of the very low population of the I state in equilibrium. Besides, the unfolding reaction is reversible, as judged by signal recovery after temperature-induced unfolding and refolding experiments 18 .
A global fitting of a two-state model, allowing for a unique value for the  m N U parameter (the dependence of unfolding free energy with denaturant concentration), was performed to experimental data of the three Cys mutants and wild-type FXN ( Fig. 2 and Table 1). Our results revealed that Leu → Cys mutations disturb the stability, though not in a great extent. This result is compatible with the fact that the side-chain of cysteine can establish interactions through its −SH group 21,22 . Both C m and T m decrease in FXN L198C and FXN L200C. On the contrary, a significant increase of C m and T m is observed for FXN L203C. Mutations L198C and L200C destabilize FXN by ~1 and 2 kcal mol −1 , respectively, while mutation L203C stabilizes FXN by ~1 kcal mol −1 .
The thermodynamic stabilization observed for the L203C mutant was unexpected. A possible explanation deduced by visually inspecting the crystallographic structure of FXN is that the introduced cysteine might establish polar interactions with the side chains of S105 and H183, and with atoms from the backbone of L200 (Fig. S5). Besides, the sulfhydryl group may also be involved in stabilizing electrostatic interactions such as hydrogen bonds, possibly leading to the alteration of the pK a of that group 21 , in agreement with molecular dynamic simulation results (data not shown).
The dependence of stability on pH was also analyzed. No significant effect of the pH in the range 6.0-8.0 was detected, suggesting that the population of molecules with their cysteine's thiol in the thiolate form is negligible, even at pH 8.0 ( Fig. S6 and Table S2). This result is consistent with in silico pK a predictions (Table S1), and allows us infer that in our regular working conditions the cysteine of the three variants is protonated and would consequently behave as a hydrophobic residue. This supports the idea that cysteine residues have the capability of satisfying to a good extent the original interactions established by leucine in wild-type FXN.
Cysteine residue modification. In order to estimate the dynamics of the C-terminal region, we used two reactive probes with different size and charge: 5,5′-dithiobis-2-nitrobenzoic acid (DTNB) and 4,4′-dithiodipyridine (DTDPy). Each probe is composed by two identical aromatic rings, linked by a disulfide (Fig. 3). DTNB consists of two molecules of 2-nitro-5-thiobenzoic acid, linked through position 5. DTDPy consists of two pyridines linked through position 4. A free sulfhydryl group, in its thiolate form can react with a molecule of DTNB or DTDPy to produce a mixed disulfide (labeled protein) and a molecule of 5-thio-2-nitrobenzoic acid ( = TNB ) or 4-thiopyridine ( − TPy ), respectively: DTDPy R S S TPy TPy (9) Note that in the case of FXN, the formation of dimers of the form FXN-S-S-FXN seems rather improbable, given we have never observed dimers despite not having treated the protein with any reducing agent. Besides, the effective concentration of free thiols that should be consider in order to establish the magnitude of probe's excess is that of the protein in which the Cys residue is in the open (solvent accessible) conformation, which is small in all our variants, as discussed later. The reaction progress is readily followed given that the hemimolecules of both probes possess different absorption spectra than the original molecule (Fig. 3). Besides, at the specific wavelengths that we used to follow the reactions (425 nm for DTNB and 324 nm for DTDPy), there is no significant protein absorption.  Table 1. Thermodynamic parameters obtained from equilibrium unfolding experiments. * A two-state model was simultaneously fitted to the data obtained in urea-induced unfolding experiments followed by CD and Trp fluorescence for all variants ( Fig. 2A,B). The value of the  m N U parameter was assumed the same for all variants, and found to be 1.8 ± 0.1 kcal mol −1 M −1 , considerably larger than the value inferred by considering the protein's length (1.5 . † A two-state model was simultaneously fitted to the data obtained in temperature-induced unfolding experiments followed by CD at 220 nm (Fig. 2C). The value of the Δ  C P N U parameter-the difference in the heat capacity between the native and unfolded states-was assumed the same for all variants and found to be 1. Kinetics of cysteine modification. Reactivity toward DTNB and DTDPy was analyzed. Cys mutants show significant differences in the observed rate coefficients values (k label ), being FXN L200C the most reactive and FXN L203C the one with the lowest modification rate (Fig. 4). The three mutants are much less reactive than their respective peptides (pL198C, pL200C or pL203C) or BME, indicating that there exists some degree of protection of the thiol toward chemical modification, irrespectively of the position of the Cys residue. Furthermore, given that the modification reaction effectively occurs, there must exist some dynamic process that allows the sulfhydryl group to be exposed only during a fraction of the time. In addition, given the molecular size of DTNB and DTDPy (van der Waals surface areas 434 and 277 Å 2 , respectively), a significant extent of Cys exposition would be necessary for the probe to react. Interestingly, the three peptides are significantly more reactive than BME and significant differences in reactivity are observed among them. Given that the peptides are only 11-residue long, this result suggests that the delicate modulation of the thiol reactivity depends on the very local molecular environment of the Cys residue. In particular, Lys2 gives the peptide a net positive charge in the N-terminal region; considering the proximity between the Lys and the Cys residue, it may be expected that reactivity of pL198C > pL200C > pL203C 22 , coincidently with our experimental observations. Further experiments will be carried out in order to establish the differences in pKa and/or nucleophilicity of the thiol for each peptide. It is worthy to mention that the use of small molecules like BME, glutathione or even cysteine should be used with caution for this kind of studies.
Values corresponding to the observed labeling coefficients, k label , were obtained for each Cys mutant at varying probe concentration. Although knowing k mod and k label at an identical reagent concentration would suffice to calculate the free energy associated to the CTR's open-close process (Equation 7), k label dependence on the reagent's concentration is necessary to establish the regime under which the reaction occurs. Moreover, given this is a bimolecular reaction, the rate of change of k label with the change in the concentration of the probe provides its concentration-independent value, which is a more accurate value as it is the result of the regression to multiple experimental data. The same criteria apply for the determination of the intrinsic modification coefficient, k mod for BME and the three model peptides pL198C, pL200C and pL203C (Fig. 4C,F). Obtained concentration-independent values of k label and k mod were used to calculate the value of ΔG° for the open-close reaction (Fig. 5 and Table 2). Additionally, the comparison between ΔG° values for the open-close process obtained using the k mod values for BME or peptide models is shown Fig. S7 (Supplementary Material). ΔG° values for the open-close reaction indicate that local conformational changes are occurring, given that the difference in free energy associated to the global unfolding process (  N U) is considerably larger (Fig. 5). In the same vein, FXN L134C shows a lower reactivity than any of the CTR mutants, despite being even less stable than FXN L200C. In fact, the difference in ΔG° for the open-close reaction is 5.3 ± 0.2 kcal mol −1 in this case (calculated by considering the k mod corresponding to the BME thiol), a value that is compatible with the complete unfolding of the protein.  (B and E). The reaction intrinsic constant, k mod , was obtained from a similar fitting to the Cys-peptides' curves (C and F). We obtained the concentration-independent second-order kinetic rate constants from the linear regression of pseudo-first order constants found each probe concentration (B, C, E and F).

Discussion
The complete understanding of protein function cannot be obtained from the mere analyses of the static three-dimensional structure. Instead, it requires a deep knowledge of the native state structural dynamics, given that function tightly depends on the molecule's internal mobility. Although the analyses based on X-ray diffraction or NMR solution structures may provide good information-particularly about spatial localization of reactive residues, interaction regions y ligand-binding sites-the study of the temporal courses, frequencies and amplitudes of intramolecular movements and fluctuations is essential in order to understand biological processes and minutely characterize function-related native substates. Experimental access to the temporal dimension and high-resolution study of substates at the secondary structure and amino acid residue level is a challenging task. One of the most popular tools used to investigate internal mobility and activated states is NMR relaxation dispersion. Besides, hydrogen-deuterium exchange allows the study of slow processes (occurring in the milliseconds and seconds time scales) and the determination of equilibrium constants of local events, but relies on auxiliary techniques such as NMR or mass spectroscopy for the identification of the residues involved in exchange processes 1 .
The fact that certain regions may experience local unfolding events that do not imply the global unfolding of the protein, suggests the existence of several minimal-energy substates: the native state ensemble. That kind of events have been described in several proteins that depend on local unfolding for a ligand to enter a cavity, such as retinol-binding protein 23 , the interaction with the active site, like in transglutaminase 24 , enzymatic catalysis itself, such as in peroxirredoxines 25 , and even as a determinant of the degradation rate 26 or of the capability to aggregate as amyloid fibrils 27 . It was also suggested that defencins-a group of short, cationic, cysteine-rich, amphiphilic peptides, capable of neutralizing several different exotoxins-may exert their action mechanism by inducing local unfolding of certain regions of several microbial and viral proteins 28 .
In this work, we studied the reactivity of sulfhydryl groups introduced in positions 198, 200 and 203 of FXN to explore the effect of the alteration of the interaction network established by the CTR of the molecule. By means of the analysis of the differential reactivity we characterized the dynamics of each position. Besides, we could calculate the change in free energy associated with the closed  open reaction in order to determine the energetics of the process involving those fluctuations.
A cysteine's reactivity depends on its thiol group nucleophilicity, its accessibility to the solvent (exposition rate) and its ionization state, which is modulated by the chemical environment (e.g. ionic interactions with neighbor residues), being the thiolate state (R − S − ) the one that may react. Naturally, a hidden, completely buried cysteine will have no access to the probe dampening the reactivity of the closed conformation (Fig. 1B). The apolar environment provided by the hydrophobic core favors the less reactive thiol state (R − SH). Instead, when the sulfhydryl group is exposed to the solvent (open conformation) its deprotonation will be a more probable event.
In wild-type FXN, L198, L200 and L203 contribute to the stability of the interaction between the core and the CTR 3,4 . The substitution of leucine for cysteine produces the shortening of the side-chain by two methyl groups and replaces the Cγ for an S atom. Such a modification necessarily implies the loss of interactions, what may destabilize the molecule 6 . However, given the strong apolar character of a protonated thiol 29 , it is possible for the original interactions established by leucine to partially prevail. In fact, pK a values are predicted higher than expected for a free solvated thiol (Table S1), and in agreement with this prediction, the conformational stability of none the Cys mutants changed significantly in response to a change in pH (Fig. S6).
The intrinsic microscopic modification coefficient (k mod ) for each mutant was evaluated using the synthetic peptides pL198C, pL200C and pL203C, which mimic the chemical environment of the corresponding Cys residue when the C-terminal region is unfolded (open state). Our results showed that BME-a completely-exposed thiol of low molecular weight-reacts significantly slower than the peptides, suggesting that reactivity is modulated by the peptidic environment under our experimental conditions. For model peptides, significant shifts of thiol pKa values may be explained by the presence of positively charged residues near the Cys side chain 22 . On the other hand, we cannot rule out a potential overestimation of the reactivity of these peptides as the result of the proximity of the N-terminus amine group. In this regard, the N-acetylation of the peptides might contribute to solve this issue.  Table 2. Local (  C O) and global (  N U) unfolding equilibrium constants. † The closed (C) to open (O) ratio and the native (N) to unfolded (U) ratio are the inverse value of the corresponding equilibrium constant. * Equilibrium constants were obtained by fitting a two-state model to the urea-induced unfolding experiments followed by CD ( Fig. 2A and Table 1).
SCIENTIFIC REPORTS | (2019) 9:1731 | https://doi.org/10.1038/s41598-019-39429-2 DTNB exhibited significantly lower modification rates than DTDPy. Consequently, Δ  G C O o values obtained for DTDPy were lower than for DTNB. However, the modification rates of the three Cys mutants followed the same pattern for both reagents. DTNB is charged and much voluminous a molecule than DTDPy, which consists in just two pyridine rings linked by a disulfide bond (Fig. 3). Therefore, DTDPy may have an easier access to partially exposed thiols, being able to detect fluctuations that imply a smaller structural change. In the same fashion, Stratton, et al. 13 showed that the kinetic of thiol-disulfide exchange involving Cys residues located at specific positions in myoglobin was highly dependent on the size of the thiosulfonate-probe reagent.
Cys exposition rate is only explainable by local fluctuations. The accessibility of probes to the thiolate might well be explained by global unfolding events. In the case of a protein that unfolds via a two-state mechanism, those events occur according to the equilibrium constant of the  N U reaction,  K N U . Such is the case of variant FXN V134C, which has a Cys residue located deep in the core of the protein, only accessible when the molecule unfolds. In native conditions, the fraction of molecules in the U state is several orders of magnitude lower than in the N state. Given that its reactivity relies on global unfolding, and is thus related to the global conformational stability, FXN V134C reacted extremely slowly.
The case of CTR Cys-mutants is contrasting. As derived from the calculated equilibrium constants (Table 2), it is not possible to explain cysteine modification just by global unfolding events, given that  K N U is two, three and three orders of magnitude lower than the equilibrium constant of the  C O reaction,  K C O , for the L198C, L200C and L203C, respectively.
For a ballpark analysis, let us consider that the dynamics of the CTR in wild-type FXN is somewhere in between the behavior of mutants FXN L198C and FXN L203C. The difference in free energy between the open and closed forms dictates that every 10 4 molecules, less than 2 exist in the open state (when the mechanism is studied by means of DTNB reactivity). This means that just a fraction of native molecules would exhibit the open conformation, while a much higher number of molecules would have the CTR packed against the rest of the protein ( Table 2).
A comprehensive model for cysteine exposition. In order to challenge the robustness of our simplified analytical approach, and to account for contributions of global unfolding processes to the apparent modification rate, we proposed a more complex model which includes the intermediate state (Fig. 6), only detectable by means of rapid-mixing kinetic studies 18 . This consideration is particularly relevant, given that the CTR is not structured in that state 18 . Therefore, a cysteine located in that region would be just as reactive in the intermediate (I) state as it would be in the unfolded (U) state. Consequently, the observed modification rate may well be explained by the  Table 3. The equilibrium constant was systematically reproduced in several fittings with varying starting values of the parameters. The values obtained for  K C O coincided with those calculated using our simplified procedure, reinforcing its validity.
The numerical fitting of the model allows us to analyze the contribution of the rate coefficients to the  K C O . It is interesting that the most notorious variation is observed for the values of k close . Not only are k close values more distinct between different mutants, but also, they are between one and four orders of magnitude higher than k open . This suggests that the local fluctuation and the correlated global stability are heavily dependent on the local refolding process, adding a complementary explanation to the conformational lock hypothesis 6 .
Simulations were useful to determine if it is correct to assume a monoexponential behavior (governed by a unique observed rate coefficient) of the labeling process (Fig. S8). We have described that if experimental conditions are those that guarantee the reaction to occur under an EX2 regime, it is possible to simplify equations enough so that the global reaction might be consider a first-order one. Given that a monoexponential function of time can be correctly fitted to the result of the simulations of the labeling reaction (still obeying the model shown in Fig. 6), it is evident that such an equation accurately explains the data (Fig. S8B).
It is noteworthy that our analysis shed light on the timescale in which the conformational changes of the studied process occurs. The inferred values for the microscopic coefficients k open indicate that these conformational changes are seemingly infrequent. Assuming a good correlation between global free energy of unfolding and k close , one may infer the corresponding k close value for wild-type FXN to be ∼47 s −1 . Cys-modification analysis allowed us to detect motions that were not previously detected by NMR. The T 1 /T 2 ratio values suggest no enhanced relaxation in the CTR under native conditions and CPMG relaxation dispersion experiments show that this stretch of wild-type FXN does not significantly contribute to the protein chemical/conformational exchange 5,30 . Although such observations are apparently in conflict with Cys-modification analysis and with results that suggest that Y205 (located in the CTR) is the first proteolytic site for wild-type FXN after a long incubation with chymotrypsin 6 , they are not unexpected given that those NMR experiments explores conformational dynamics in the pico to milliseconds time scales. Remarkably, the fluctuations involving the CTR in our working model are consistent with processes occurring in the order of seconds to minutes (Table 3).
One may consider that the differences in reactivity found for our mutants might be the result of very large perturbations of the hydrophobic interactions occurring between the CTR and the core of FXN; consequently, the behavior observed for this set of mutants might not be illustrative of the motions and local unfolding events that exist in wild-type FXN. However, these mutants are quite stable from the thermodynamic viewpoint, revealing that the interactions that CTR establishes are substantially preserved. This fact, added to the finding of stabilizing (L203C) and destabilizing mutations (L198C and L200C), and to the existence of a trend between global and local stabilities, might allow us to estimate local unfolding and CTR dynamics of wild-type FXN, therefore gaining significant information about its native motions.
Finally, our results indicate that there is a relationship between local stability of the CTR and global stability of FXN and this is consistent with the strong loss of stability and increase of mobility of the CTR seen in the pathogenic variant FXN L198R 5 . It has been shown that a number of pathogenic mutations that produce Friedreich's ataxia involves protein variants that are highly destabilized showing significantly lower melting temperatures and lower free energy of unfolding by comparison to the wild-type FXN, among them: G130V (ΔT m = 23 °C and ΔΔG°N U ∼2.9 kcal mol −1 ) 31 , G137V (ΔT m ~20 °C) 32 and L198R (ΔT m = 15.3 °C and ΔΔG°N U ∼4.2 kcal mol −1 ) 5 . Local stabilization of the CTR in these cases may produce an increase in global stability of the mutant proteins. This means that the CTR may be a promising target for drug design and therapeutic strategies.
Concluding remarks. In this paper, we addressed the study of local fluctuations, derived from native state dynamics, by means of the analysis of the reactivity of cysteine in a set of rationally-designed point mutants.  Table 3. Best-fitting values of the parameters of the model. The model presented in Fig. 6 was fitted simultaneously to DTNB modification time courses of each Cys-mutant and their corresponding Cys-peptide using COPASI 34 . Initial values were those corresponding to the equilibrium distribution of species in the absence of probe. Parameters are presented as value ± 1 S.E. Fixed parameter values were the following 18  We proved the existence of local conformational changes, which imply the exposition of the side-chains of residues that are relevant for global stability and are hidden from the solvent in the most populated native substate. Furthermore, we could get some insight into the energetics of that local process and quantified the local stability of FXN's CTR for each of the mutants. Our results show that chemical modification can be efficiently used to assess the existence of native substates that characterize the dynamics and affect stability of protein molecules.