# Identification and Structural Characterization of an Intermediate in the Folding of the Measles Virus X Domain

Daniela Bonetti, Carlo Camilloni, Lorenzo Visconti, Sonia Longhi, Maurizio Brunori, Michele Vendruscolo, Stefano Gianni

<figure><img src="image_1.png" /></figure>

## To cite this version:

Daniela Bonetti, Carlo Camilloni, Lorenzo Visconti, Sonia Longhi, Maurizio Brunori, et al.. Identification and Structural Characterization of an Intermediate in the Folding of the Measles Virus X Domain. Journal of Biological Chemistry, 2016, 291 (20), pp.10886+. 10.1074/jbc.M116.721126. hal-01439066

HAL Id: hal-01439066

https://hal.science/hal-01439066v1

Submitted on 11 Jul 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# Identification and Structural Characterization of an Intermediate in the Folding of the Measles Virus X Domain*§

Received for publication, February 10, 2016, and in revised form, March 17, 2016 Published, JBC Papers in Press, March 21, 2016, DOI 10.1074/jbc.M116.721126

Daniela Bonetti‡, Carlo Camilloni§¶, Lorenzo Visconti‡, Sonia Longhi||**, Maurizio Brunori‡, Michele Vendruscolo§, and Stefano Gianni‡§<fnref n="1" />

From the ‡Istituto Pasteur-Fondazione Cenci Bolognetti and Istituto di Biologia e Patologia Molecolari del CNR, Dipartimento di Scienze Biochimiche "A. Rossi Fanelli," Sapienza Università di Roma, 00185 Rome, Italy, ¶Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom, ‡Department of Chemistry and Institute for Advanced Study, Technische Universität München, Lichtenbergstrasse 4, D-85747 Garching, Germany, ||Aix-Marseille Université, Architecture et Fonction des Macromolécules Biologiques (AFMB), CNRS, UMR 7257, 13288 Marseille, France, and **CNRS, AFMB UMR 7257, 13288 Marseille, France

Although most proteins fold by populating intermediates, the transient nature of such states makes it difficult to characterize their structures. In this work we identified and characterized the structure of an intermediate of the X domain of phosphoprotein (P) of measles virus. We obtained this result by a combination of equilibrium and kinetic measurements and NMR chemical shifts used as structural restraints in replica-averaged metadynamics simulations. The structure of the intermediate was then validated by rationally designing four mutational variants predicted to affect the stability of this state. These results provide a detailed view of an intermediate state and illustrate the opportunities offered by a synergistic use of experimental and computational methods to describe non-native states at atomic resolution.

A major goal in protein folding studies lies in identifying intermediates and in characterizing their structure (1–8). This task is often challenging because, as marginally stable, folding intermediates are elusive and can be detected only transiently during the folding reaction. The combination of experimental and theoretical approaches represents a powerful strategy to address this problem, providing the possibility to unveil the molecular structure of folding intermediates and to define their role in the reaction mechanism.

The viral polymerase complex of paramyxoviruses is composed by the protein L and by the phosphoprotein P (9–11). In measles virus (a member of the Morbillivirus genus), the C-terminal region of phosphoprotein P is a globular domain of 49 amino acids, called X domain (XD),<fnref n="2" /> composed of three α-helices organized as an anti-parallel bundle (12). Previous investigations on measles virus XD have suggested this domain to be structurally heterogeneous, populating at least two alternative conformations under native conditions (13). This feature is consistent with the findings of Kingston *et al.* (14), who suggested that the native state

of XD from mumps virus (a *Rubulavirus* member) represents an example of folding induced by crystal packing effects. Notably, the lack of a unique stable three-dimensional structure is not a feature unique to mumps virus XD, being also shared by the corresponding domains from other rubulaviruses that were found to populate in solution a structural continuum ranging from stable α-helical bundles to largely disordered (15). Because of such heterogeneity, XD represents an ideal candidate for folding studies, as it is likely to populate partially structured states at equilibrium, displaying an overall stability comparable with that of the native conformation.

In this work we characterized the folding pathway of XD, and present the high resolution structure of an intermediate while characterizing its mechanistic role in time-resolved kinetic experiments. The occurrence of a folding intermediate is supported by the observation of bi-phasic (un)folding transitions in discharge-induced temperature-jump experiments as well as by all-atom replica-averaged metadynamics simulations with NMR chemical shift restraints. In an effort to further address the structure of the intermediate state, we designed rationally and characterized experimentally mutational variants that selectively tune the stability of this state with respect to that of the native state.

## Experimental Procedures

Experiments were performed on a fluorescent pseudo-wild type XD variant Y480W, which was previously produced and characterized. All proteins were expressed and purified as described (16). All reagents were of analytical grade.

*Circular Dichroism Equilibrium Experiments*—Circular dichroism equilibrium denaturation experiments were carried out on a JASCO circular dichroism (CD) spectropolarimeter (Jasco, Inc., Easton, MD). CD spectra were recorded between 200 and 250 nm using a quartz cuvette with a light path of 1 cm, at different urea concentrations. For each urea concentration three spectra were averaged. Protein concentration was typically 10 µM, and the buffer used was 50 mM sodium phosphate, 300 mM NaCl at pH 7.2. The experiments were performed at 25 °C.

*Fluorescence Equilibrium Experiments*—Fluorescence equilibrium denaturation experiments were carried out using a Fluoromax single photon counting spectrofluorometer (Jobin-Yvon, Edison, NJ). Emission spectra were recorded between 300 and 400 nm using a quartz cuvette with a light path of 1 cm, at different urea concentrations. The excitation wavelength was

<footnote>* This work was supported by the Italian Ministero dell'Istruzione dell'Università e della Ricerca (Progetto di Interesse "Invecchiamento"; to S. G.) and by Sapienza Università di Roma (C26A155S48; to S. G.). The authors declare that they have no conflicts of interest with the contents of this article.</footnote>

<footnote>§ This article contains supplemental Figs. S1 and S2.</footnote>

<footnote n="1">¹ To whom correspondence should be addressed. E-mail: stefano.gianni@uniroma1.it.</footnote>

<footnote n="2">² The abbreviations used are: XD, X domain; D, denatured state; I, intermediate state; N, native state.</footnote>
An Intermediate in the Folding of the Measles Virus X Domain

280 nm. Protein concentration was typically 5 µM, and the buffer used was 50 mM sodium phosphate, 300 mM NaCl at pH 7.2. The experiments were performed at 25 °C.

*Temperature-jump Fluorescence Spectroscopy*—Kinetic (un)folding experiments were performed using a Hi-Tech PTJ-64 capacitor-discharge T-jump apparatus (Hi-Tech, Salisbury, UK). Temperature was rapidly changed with a jump-size of 9 °C, from 16 °C to 25 °C. The fluorescence change of N-acetyltryptophanamide (NATA) was used in control measurements. Degassed samples were slowly pumped through the 0.5 × 2-mm quartz flow cell before data acquisition. Usually 10 individual traces were averaged. The excitation wavelength was 296 nm, and the fluorescence emission was collected using a 320-nm cut-off glass filter. Protein concentration was typically 100 µM, and the buffer used was 50 mM sodium phosphate, 300 mM NaCl at pH 7.2.

*Replica-averaged Metadynamics Simulations*—Replica-averaged metadynamics simulations (17, 18) were performed using GROMACS compiled with PLUMED and ALMOST (19). The system was simulated using the Amber03W force field (20) in explicit TIP4P05 water (21). Van der Waals and short-range electrostatic interactions were cut off at 0.9 nm, whereas long range electrostatic interaction was treated with the Particle Mesh Ewald method and a mesh size of 0.12 nm (22). The isothermal-isobaric ensemble was enforced using the Bussi thermostat (23) at 300 K and the Parrinello-Rahman barostat (24) at 1 bar. The starting conformation for XD was taken from the 2K9D NMR solution structure (25). The structure was solvated with 5845 water molecules and 3 chloride ions.

Replica-averaged metadynamic simulations were performed using chemical shifts as replica-averaged restraints and bias-exchange metadynamics (18). Four replicas of the system were simulated in parallel with a restraint applied on the square difference between the CamShift (26) back-calculated NMR chemical shifts for the Hα and HN hydrogens and the experimental data using a force constant of 6 and 3 kJ/(mol ppm²), respectively. Each of the four replicas is bias along one of the following four collective variables: α-helical content, radius of gyration, root mean square deviation from the crystal structure calculated over the Cα carbons, and number of contacts among the heavy atoms of the hydrophobic residues. Gaussian deposition was performed with an initial rate of 0.125 kJ/mol/ps, a bias factor of 10, and with σ values set to 0.2, 0.005, 0.1, and 0.25, respectively. Each replica was been evolved for 180 ns, with exchange trials every 50 ps.

Results and Discussion

*Detection of an Intermediate in the Folding Pathway of XD*—Urea-induced denaturation of XD was monitored both by far-UV CD and by Trp fluorescence emission (Fig. 1A). In each data set the apparent equilibrium transition was consistent with a two-state sigmoidal transition (27); however, a comparison between these different techniques highlighted an additional complexity. In fact, the denaturation midpoint shows a small difference, 5.4 ± 0.1 M from CD experiments and 5.0 ± 0.1 M from fluorescence, with a similar *m* value of ~0.7 ± 0.1 kcal mol⁻¹ M⁻¹. These differences, which are just larger than the sensitivity of the experiments, would suggest the presence

<figure><img src="image_4.png" /><figcaption>FIGURE 1. Equilibrium and time resolved unfolding of XD, the X domain of the phosphoprotein P, as follows at pH 7.2, 50 mM sodium phosphate, 300 mM NaCl, and 25 °C. Panel A, urea-induced equilibrium denaturation of XD, monitored by fluorescence (close circles) and CD (open circles). Both fluorescence and CD, recorded at multiple wavelengths, were fitted to a two-state (continuous lines). Panel B, temperature jump induced unfolding transition of XD, induced by a rapid jump of 9 °C to a final temperature of 25 °C.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>FIGURE 2. Chevron plot of XD. The urea dependence of the fast λ₁ (closed circles) and slow λ₂ (open circles) folding phases of XD is reported. The lines are the best global fit to an on-pathway (black) and off-pathway (gray) scenario. It must be noted that the off-pathway model would imply the microscopic rate constant k₁D (for the transition from the intermediate to the denatured state) to be necessarily faster than the observed rate constant λ₂, accounting for a maximal stability of I versus D of 0.4 ± 0.3 kcal mol⁻¹.</figcaption></figure>

of an equilibrium folding intermediate (28) and demand additional investigation.

The (un)folding time course of XD in the absence of denaturant, measured by T-jump, is reported in Fig. 1B. The time-resolved fluorescence conforms to a double exponential decay, the two kinetic phases displaying, in the absence of denaturant, similar relative amplitude but opposite sign. This finding confirms unequivocally that the folding of XD implies the presence of an intermediate, which is remarkably populated at equilibrium also in the absence of denaturant. Importantly, because the observed transition displays two distinct kinetic phases, it follows that the intermediate state displays a fluorescence emission, which is different from that of the denatured and native states.

The urea dependence of the two kinetic rate constants λ₁ and λ₂ is reported in Fig. 2 (chevron plot). Given the relatively low
An Intermediate in the Folding of the Measles Virus X Domain

<table><caption>TABLE 1<br>Kinetic folding parameters of XD</caption><thead><tr><th></th><th>k<sub>DI</sub><sup>a</sup><br>s<sup>-1</sup></th><th>k<sub>ID</sub><sup>a</sup><br>s<sup>-1</sup></th><th>k<sub>IN</sub><sup>a</sup><br>s<sup>-1</sup></th><th>k<sub>IN</sub><sup>a</sup><br>s<sup>-1</sup></th><th>k<sub>DN</sub><sup>a</sup><br>s<sup>-1</sup></th><th>k<sub>DN</sub><sup>a</sup><br>s<sup>-1</sup></th><th>ΔG<sub>DI</sub><br>kcal mol<sup>-1</sup></th><th>ΔG<sub>DN</sub><br>kcal mol<sup>-1</sup></th></tr></thead><tbody><tr><td>On pathway</td><td>5700 ± 500 (0.2 ± 0.05)</td><td>2.9 ± 0.7 (0.8 ± 0.05)</td><td>480 ± 70 (0.04 ± 0.05)</td><td>94 ± 10 (0.2 ± 0.03)</td><td>4.2 ± 0.9 (1.0 ± 0.07)</td><td>0.9 ± 0.1 (0.24 ± 0.07)</td><td></td><td></td></tr><tr><td>Off pathway</td><td>1400 ± 300 (0.3 ± 0.05)</td><td>683 ± 60 (0.16 ± 0.05)</td><td></td><td></td><td>0.4 ± 0.1 (0.46 ± 0.07)</td><td>2.3 ± 0.5 (0.30 ± 0.15)</td><td></td><td></td></tr></tbody><tfoot><tr><td colspan="9"><sup>a</sup> The microscopic rate constants were calculated by globally fitting the two observed relaxation rate constants to the analytical solution of a three-state mechanism involving an off-pathway intermediate. The m value associated to each rate constant, expressed in kcal mol<sup>-1</sup>, is reported in parentheses.</td></tr></tfoot></table>

<figure><img src="image_6.png" /><figcaption>FIGURE 3. Free energy surface (in kJ/mol) of the replica-averaged meta-dynamic XD ensemble as a function of α-helical content and number of hydrophobic contacts. The projection highlights a 70% populated ground state and a 30% populated intermediate state. The two states are represented as an ensemble of structures (top) as well as an average structure (bottom).</figcaption></figure>

stability of the intermediate, we could not observe the fast phase $\lambda_1$ at urea concentrations higher than 3 M. A folding intermediate may represent either an on-pathway species en route to the native state or a misfolded off-pathway kinetic trap competing with productive folding, with the two different scenarios displaying characteristic signatures in the corresponding chevron plots (2, 3, 7, 29–34). In fact, although the off-pathway kinetic model implies the microscopic rate constant $k_{\text{ID}}$ (i.e. the rate constant for the unfolding of the intermediate to the denatured state) as necessarily faster than the observed rate constant $\lambda_2$, the on-pathway model does not imply this constraint. To quantitatively describe the folding pathway of XD, we globally fitted the chevron plots for the two observed relaxation rate constants. As shown in Fig. 2, both the on- and off-pathway models fit fairly well with the experimental data and, therefore, remain indistinguishable at this stage. The main difference between the two models lies in the lower stability of the intermediate state predicted by the off-pathway model, corresponding to a maximal difference in stability between I and D being 0.4 ± 0.1 kcal
An Intermediate in the Folding of the Measles Virus X Domain

<figure><img src="image_7.png" /><figcaption>FIGURE 4. **Urea-induced denaturation of the XD variants.** Data were recorded in 50 mm sodium phosphate at pH 7.2, 300 mm NaCl, and 25 °C, monitored by fluorescence (*closed circles*) and CD (*open circles*). All the lines are the best fit to a two-state model.</figcaption></figure>

mol⁻¹, which in turn suggests a negligible population of the intermediate state in the absence of denaturant (≪1% of the molecules). The calculated rate constants from the global analysis arising from the two curve-fitting are listed in Table 1 together with the predicted stability of the intermediate and native states for the two different scenarios. Analysis of the calculated *m* values, reported in Table 1, suggests that the intermediate state is relatively compact, the β~T~ values being 0.8 ± 0.1 and 0.6 ± 0.1, as calculated from the on-pathway and the off-pathway scenarios, respectively.

*Structure of the Folding Intermediate—*To provide a structural characterization of the folding intermediate, we performed all-atom replica-averaged metadynamics simulations with the previously published NMR chemical shift of XD as structural restraints (17, 18, 35, 36). This approach is based on the observation that, given that the relaxation kinetics suggests that the intermediate is significantly populated even in the absence of denaturant, it may be possible to detect the structure of the intermediate in the structural ensemble of XD, obtained under native conditions.

Simulations were performed at 300 K using chemical shifts as replica-averaged structural restraints over 4 replicas and bias-exchange metadynamics to enhance the sampling, as previously described (17, 18, 35, 36). Metadynamics was applied on four collective variables, namely the α-helical content, the radius of gyration (*R*~g~), the root mean square deviation from the crystal structure calculated over the Cα carbons, and the number of

contacts among the heavy atoms of the hydrophobic residues. The conformational search was limited up to 6 Å from the crystal structure to avoid a complete unfolding of the protein. After 140 ns of simulation per replica the free energies as a function of the collective variables were converged, and the simulations were extended up to 180 ns per replicas.

The effect of the bias on the simulations was removed using a weighted-histogram analysis technique (37), and the resulting ensemble was clustered in the four-dimensional space defined by the selected collective variables. We thus determined the presence of two major states populated at ~70 and 30%, respectively (Fig. 3). The most relevant population of the ensemble is structurally homogeneous and superposes well onto the native state of XD, as determined by x-ray crystallography (PDB code 1OKS) (12). On the other hand, the alternative conformational ensemble of XD, although characterized by a native-like content of secondary structure, displays a different structural organization of the three-helix bundle, with a much stronger heterogeneity. This less populated state was taken to reflect the structure of the intermediate state.

To validate the structural ensembles described above, the chemical shifts for both the H~N~ and Hα atoms were calculated and compared with those experimentally determined by NMR and those calculated from the crystal structure (supplemental Fig. S1). It is evident that the values from the ensemble obtained by metadynamics resemble very closely those obtained by NMR and are in significantly better agreement than those calculated
An Intermediate in the Folding of the Measles Virus X Domain

<figure><img src="image_8.png" /><figcaption>FIGURE 5. Temperature jump induced unfolding transition of the XD variants measured in 50 mM sodium phosphate at pH 7.2, 300 mM NaCl, and at a final temperature of 25 °C. Unfolding was induced by a rapid increase of 9 °C. Although the R465A and the D486N mutants display a double exponential time course, in the case of I504A and I505A, data were consistent with a single exponential time course.</figcaption></figure>

from the crystal structure. This finding suggests that some effects of the dynamics of the ensemble are lost in the crystallization, as previously suggested for the XD of mumps virus (14).

Validation of the Three-dimensional Structure of the Folding Intermediate by Tuning Its Stability—To further support the structural model, we adopted the strategy to rationally design site-directed mutants aimed at tuning its thermodynamic stability. A comparison between the intermediate and the native conformations highlights a reorganization of the α-helix bundle while maintaining a similar α-helical content. In particular, we observed that positions Arg-465 and Asp-486, located at the interface between helix1 and helix2, are involved in a salt bridge in the native state, whereas they are not in contact in the intermediate state; furthermore, positions I504 and I505, located at the C-terminal end of the domain, constitute a structural link between the three α-helices. Therefore, we produced the D486N, R465A, I504A, and I505A variants. Fig. 3 highlights the different environment of these four residues in the two states.

The urea-induced equilibrium (un)folding transitions followed by CD and fluorescence for the four variants are reported in Fig. 4. In all cases each individual transition conforms to an apparent two-state process similar to that reported above for the wild type protein. However, in the case of the R465A and D486N variants there was a large difference between the observed midpoints of urea denaturation as observed by fluorescence and CD and the values of [urea]₁/₂, 6.8 ± 0.3 M from CD and 4.8 ± 0.1 M from fluorescence for R465A, and 6.3 ± 0.2 M

from CD and 4.6 ± 0.1 M from fluorescence in the case of D486N. By contrast, in the case of the I504A and I505A variants, the overall transitions and urea midpoints followed by CD and fluorescence are superimposable within experimental error. These findings indicate that, although the intermediate state is stabilized in the D486N and R465A variants, the I504A and I505A variants conform to a two-state transition. Because Asp-486 and Arg-465 are engaged in a salt bridge in the native but not in the intermediate state, it is likely that the observed stabilization of the intermediate in the D486N and R465A variants arises from a relative destabilization of the native state relative to I.

In light of these experimental observations, it is of interest to compare the emission spectra of the different variants to those of wild type XD (supplemental Fig. S2). In fact, both the I504A and I505A mutants, which display a two-state folding transition, display a different λmax compared with the other variants. This finding supports the view that in the case of wild type XD, R465A, and D486N, the intermediate state is partially populated, whereas in the case of the I504A and I505A there is a shift toward the native state only.

Additional evidence for the stabilization of the intermediate in the D486N and R465A variants and its destabilization in the case of the I504A and I505A mutants was then provided by assessing the (un)folding kinetics of these variants. As shown by time-resolved temperature jump transition of the four different mutants (Fig. 5), in the case of D486N and R465A, the relative
An Intermediate in the Folding of the Measles Virus X Domain

<figure><img src="image_9.png" /></figure>

FIGURE 6. **Chevron plot of XD variants.** Urea dependence of the fast $\lambda_1$ (closed circles) and slow $\lambda_2$ (open circles) unfolding phases of XD. The lines are the best global fit to an on-pathway scenario. For simplicity, data were fitted by constraining the *m* values to those observed in wild type XD. As described under "Results and Discussion," in the case of the I504A and I505A variants, only the slow phase $\lambda_2$ could be observed. Accordingly, in these cases, the chevron was fitted to a two-state model.

amplitude of the fast *versus* the slow unfolding phase increases with respect to that of the native state. On the other hand, both I504A and I505A (un)fold via an apparent single-exponential time course, consistent with a two-state transition. The chevron plots for these four variants are shown in Fig. 6.

*Conclusions*—We have identified and characterized an intermediate state of the X domain of phosphoprotein P, which is populated for ~30% of the time. To obtain this result we exploited the recent development of NMR-restrained all-atom replica-averaged metadynamics simulations, which allow the structures of weakly populated states under-native conditions to be determined (18, 38).

Through a comparison with the thermodynamic parameters obtained from chevron analysis, we found that the intermediate state is an on-pathway species en route to the native state rather than an off-pathway kinetic trap. Although in the case of the former model, the intermediate is estimated to be populated in the absence of denaturant by ~30% of the molecules, the latter model would predict the intermediate to be only weakly populated under native conditions (<<1%) due to its intrinsic constraint(s) (i.e. $k_{1D} > \lambda_2$). The calculated structure of the intermediate is fairly compact and displays a native-like secondary structure content. The α-helix bundle, however, is characterized by a different overall topological organization. This finding is fully consistent with equilibrium denaturation experiments, which indicate that the intermediate displays a native-like

$\alpha$-helical content, as probed by CD, while characterized by a different fluorescence emission. Moreover, mutations designed with the specific aim of perturbing the interfaces between the $\alpha$-helices modulated the stability of the intermediate as predicted such that folding can be tuned toward a simple two-state reaction. Future work based on site-directed mutagenesis will elucidate the mechanistic details of the three-state folding transition of XD and its possible biological significance.

*Author Contributions*—D. B. and L. V. conducted the experimental work. C. C. determined the structure of the intermediate. D. B, C. C., L. V., S. L., M. B., M. V., and S. G, analyzed and discussed the data. S. G. conceived the idea of the project and wrote the main draft of the manuscript, which was later edited by all the authors.

References

1. Capaldi, A. P., Kleanthous, C., and Radford, S. E. (2002) Im7 folding mechanism: misfolding on a path to the native state. *Nat. Struct. Biol.* **9**, 209–216

2. Capaldi, A. P., Shastry, M. C., Kleanthous, C., Roder, H., and Radford, S. E. (2001) Ultrarapid mixing experiments reveal that Im7 folds via an on-pathway intermediate. *Nat. Struct. Biol.* **8**, 68–72

3. Gianni, S., Ivarsson, Y., De Simone, A., Travaglini-Allocatelli, C., Brunori, M., and Vendruscolo, M. (2010) Structural characterization of a misfolded intermediate populated during the folding process of a PDZ domain. *Nat. Struct. Mol. Biol.* **17**, 1431–1437

4. Matouschek, A., Kellis, J. T. Jr., Serrano, L., Bycroft, M., Fersht, A. R. (1990)
An Intermediate in the Folding of the Measles Virus X Domain

Transient folding intermediates characterized by protein engineering. *Nature* **346**, 440–445

5. Parker, M. J., Spencer, J., and Clarke, A. R. (1995) An integrated kinetic analysis of intermediates and transition states in protein folding reactions. *J. Mol. Biol.* **253**, 771–786

6. Religa, T. L., Markson, J. S., Mayor, U., Freund, S. M., and Fersht, A. R. (2005) Solution structure of a protein denatured state and folding intermediate. *Nature* **437**, 1053–1056

7. Wildegger, G., and Kiefhaber, T. (1997) Three-state model for lysozyme folding: triangular folding mechanism with an energetically trapped intermediate. *J. Mol. Biol.* **270**, 294–304

8. Gruebele, M. (2002) An intermediate seeks instant gratification. *Nat. Struct. Biol.* **9**, 154–155

9. Bourhis, J. M., Canard, B., and Longhi, S. (2006) Structural disorder within the replicative complex of measles virus: functional implications. *Virology* **344**, 94–110

10. Bourhis, J. M., Johansson, K., Receveur-Bréchot, V., Oldfield, C. J., Dunker, K. A., Canard, B., and Longhi, S. (2004) The C-terminal domain of measles virus nucleoprotein belongs to the class of intrinsically disordered proteins that fold upon binding to their physiological partner. *Virus Res.* **99**, 157–167

11. Sedlmeier, R., and Neubert, W. (1998) The replicative complex of paramyxoviruses: structure and function. *Adv. Virus. Res.* **50**, 101–139

12. Johansson, K., Bourhis, J.-M., Campanacci, V., Cambillau, C., Canard, B., and Longhi, S. (2003) Crystal structure of the measles virus phosphoprotein domain responsible for the induced folding of the C-terminal domain of the nucleoprotein. *J. Biol. Chem.* **278**, 44567–44573

13. D'Urzo, A., Konijnenberg, A., Rossetti, G., Habchi, J., Li, J., Carloni, P., Sobott, F., Longhi, S., and Grandori, R. (2015) Molecular basis for structural heterogeneity of an intrinsically disordered protein bound to a partner by combined ESI-IM-MS and modeling. *J. Am. Soc. Mass. Spectrom.* **26**, 472–481

14. Kingston, R. L., Gay, L. S., Baase, W. S., and Matthews, B. W. (2008) Structure of the nucleocapsid-binding domain from the mumps virus polymerase; an example of protein folding induced by crystallization. *J. Mol. Biol.* **379**, 719–731

15. Yegambaram, K., Bulloch, E. M., and Kingston, R. L. (2013) Protein domain definition should allow for conditional disorder. *Protein Sci.* **22**, 1502–1518

16. Dosnon, M., Bonetti, D., Morrone, A., Erales, J., di Silvio, E., Longhi, S., and Gianni, S. (2015) Demonstration of a folding after binding mechanism in the recognition between the measles virus NTAIL and X domains. *ACS Chem. Biol.* **10**, 795–802

17. Camilloni, C., Cavalli, A., and Vendruscolo, M. (2013) Replica-averaged metadynamics. *J. Chem. Theory Comput.* **9**, 5610–5617

18. Camilloni, C., and Vendruscolo, M. (2014) Statistical mechanics of the denatured state of a protein using replica-averaged metadynamics. *J. Am. Chem. Soc.* **136**, 8982–8991

19. Fu, B., Sahakyan, A. B., Camilloni, C., Tartaglia, G. G., Paci, E., Caflisch, A., Vendruscolo, M., and Cavalli, A. (2014) ALMOST: an all atom molecular simulation toolkit for protein structure determination. *J. Comput. Chem.* **35**, 1101–1105

20. Best, R. B., and Mittal, J. (2010) Protein simulations with an optimized water model: cooperative helix formation and temperature-induced unfolded state collapse. *J. Phys. Chem. B* **114**, 14916–14923

21. Abascal, J. L., and Vega, C. (2005) A general purpose model for the condensed phases of water: TIP4P/2005. *J. Chem. Phys.* **123**, 234505

22. Essmann, U., Perera, L., Berkowitz, M. L., Darden, T., Lee, H., and PederSEN, L. G. (1995) A smooth particle mesh Ewald method. *J. Chem. Phys.* **103**, 8577–8593

23. Bussi, G., Donadio, D., and Parrinello, M. (2007) Canonical sampling through velocity rescaling. *J. Chem. Phys.* **126**, 014101–014105

24. Parrinello, M., and Rahman, A. (1981) Polymorphic transitions in single crystals: a new molecular dynamics method. *J. App. Phys.* **52**, 7182–7190

25. Bernard, C., Gely, S., Bourhis, J. M., Morelli, X., Longhi, S., and Darbon, H. (2009) Interaction between the C-terminal domains of N and P proteins of measles virus investigated by NMR. *FEBS Lett.* **583**, 1084–1089

26. Kohlhoff, K. J., Robustelli, P., Cavalli, A., Salvatella, X., and Vendruscolo, M. (2009) Fast and accurate predictions of protein NMR chemical shifts from interatomic distances. *J. Am. Chem. Soc.* **131**, 13894–13895

27. Pace, C. N., and Shaw, K. L. (2000) Linear extrapolation method of analyzing solvent denaturation curves. *Proteins* **4**, 1–7

28. Jackson, S. E., Fersht, A. R. (1991) Folding of chymotrypsin inhibitor 2. 1: evidence for a two-state transition. *Biochemistry* **30**, 10428–10435

29. Ivarsson, Y., Travaglini-Allocatelli, C., Jemth, P., Malatesta, F., Brunori, M., and Gianni, S. (2007) An on-pathway intermediate in the folding of a PDZ domain. *J. Biol. Chem.* **282**, 8568–8572

30. Jemth, P., Gianni, S., Day, R., Li, B., Johnson, C. M., Daggett, V., Fersht, A. R. (2004) Demonstration of a low-energy on-pathway intermediate in a fast-folding protein by kinetics, protein engineering, and simulation. *Proc. Natl. Acad. Sci. U.S.A.* **101**, 6450–6455

31. Khorasanizadeh, S., Peters, I. D., and Roder, H. (1996) Evidence for a three-state model of protein folding from kinetic analysis of ubiquitin variants with altered core residues. *Nat. Struct. Biol.* **3**, 193–205

32. Travaglini-Allocatelli, C., Gianni, S., Dubey, V. K., Borgia, A., Di Matteo, A., Bonivento, D., Cutruzzolà, F., Bren, K. L., and Brunori, M. (2005) An obligatory intermediate in the folding pathway of cytochrome c552 from *Hydrogenobacter thermophilus*. *J. Biol. Chem.* **280**, 25729–25734

33. Travaglini-Allocatelli, C., Gianni, S., Morea, V., Tramontano, A., Soulimane, T., and Brunori, M. (2003) Exploring the cytochrome c folding mechanism: cytochrome c552 from thermus thermophilus folds through an on-pathway intermediate. *J. Biol. Chem.* **278**, 41136–41140

34. Gianni, S., Ivarsson, Y., Jemth, P., Brunori, M., and Travaglini-Allocatelli, C. (2007) Identification and characterization of protein folding intermediates. *Biophys. Chem.* **128**, 105–113

35. Camilloni, C., Robustelli, P., De Simone, A., Cavalli, A., and Vendruscolo, M. (2012) Characterisation of the conformational equilibrium between the two major substates of RNase A using NMR chemical shifts. *J. Am. Chem. Soc.* **134**, 3968–3971

36. Cavalli, A., Camilloni, C., and Vendruscolo, M. (2013) Molecular dynamics simulations with replica-averaged structural restraints generate structural ensembles according to the maximum entropy principle. *J. Chem. Phys.* **138**, 094112

37. Kumar, S., Rosenberg, J. M., Bouzida, D., Swendsen, R. H., and Kollman, P. A. (1992) The weighted histogram analysis method for free-energy calculations on biomolecules. I. The method. *J. Comput. Chem.* **13**, 1011–1021

38. Lindorff-Larsen, K., Best, R. B., Depristo, M. A., Dobson, C. M., and Vendruscolo, M. (2005) Simultaneous determination of protein structure and dynamics. *Nature* **433**, 128–132