Site-specific 13C isotopically edited infrared (IR) spectroscopy combined with an Ising-like (Wako–Saitô–Muñoz–Eaton) statistical-mechanical model showed in 2018 that two point mutations of the de novo designed helix-turn-helix motif αtα — a stabilizing E15L near the turn and a destabilizing L32V in the hydrophobic core — change the stability of all probed local segments in the same direction, leaving the relative order of local thermal unfolding essentially unchanged. This supported the conclusion that the native structure, not the sequence details, dominates the folding mechanism of this elementary motif, and that only a full re-parameterization of the coarse-grained model (contact energy ε\varepsilon and per-residue entropy loss Δs\Delta s) can fit the mutant data. The open problem is to establish the robustness and the physical-chemical origin of this 'global rather than local' mutation effect: whether it survives a more detailed description of the mutant structures and interactions, whether it generalizes to other positions and more disruptive substitutions, whether the folding pathway itself remains unchanged, and whether the mutation-induced changes of the coarse-grained parameters can be given a physical-chemical basis.

Background

How a linear amino-acid sequence folds into its native three-dimensional structure, and how the folding mechanism is distributed among local segments of that structure, remains a central question in biophysical chemistry. Elementary structural motifs — the smallest protein fragments that fold autonomously — are the natural test systems, and the helix-turn-helix (HTH) motif, two α-helices packed against each other and connected by a turn, is one of the simplest in which secondary structure (helices) and tertiary structure (interhelical contacts) can be studied together. The peptide αtα is a de novo designed HTH motif whose NMR solution structure is known (Protein Data Bank entry 1ABZ); its two isolated helix fragments are unfolded, so its folding is driven by interhelical tertiary contacts, making it a minimal model for cooperativity between secondary and tertiary structure.

The experimental methodology that motivates this problem is site-specific isotopically edited IR spectroscopy: selectively substituting 13C at chosen amide carbonyls shifts their amide I' (C=O stretch) absorption to a distinct frequency band, so the folding of several individual residues or segments can be monitored separately in a single thermal unfolding experiment, alongside circular dichroism (CD), which reports average α-helical content. Applied to αtα at multiple labeled positions, this revealed a heterogeneous, non-two-state unfolding pattern: the central parts of both helices (the hydrophobic core) are most stable, while the turn and the helix termini are least stable.

These multi-probe data were analyzed with an Ising-like statistical-mechanical model of folding — in this context equivalent to the Wako–Saitô–Muñoz–Eaton (WSME) model, in which each residue is a binary variable (native-like or disordered), only native (structure-derived) contacts contribute favorable energy, and the partition function can be evaluated exactly or under approximations. The variant used here is the double sequence approximation with loops (DSA/L), which allows up to two independently folded segments that can interact across a disordered loop. Remarkably, a parametrization with a single contact energy ε\varepsilon and a single conformational-entropy cost Δs\Delta s per residue, with the native contact map weighted by its fractional occurrence in the NMR ensemble, quantitatively reproduced both the CD data and the 13C IR data of multiple isotopologues — not only for αtα but also for a second, naturally occurring HTH motif (the P22 repressor subdomain) — without any sequence-specific residue–residue interactions (Miyazawa–Jernigan statistical potentials) or explicit hydrogen-bond terms. This suggested that native topology, not sequence detail, determines the local stability pattern. Independently, all-atom folding simulations of wild-type αtα with residue-specific force fields (RSFF1/RSFF2) reproduced the NMR structure, site-specific helix stabilities, and a 'zipping-out' folding mechanism initiated at the inner helix segments — an all-atom confirmation available for the wild type only.

The 2018 study then asked how this picture responds to perturbing the sequence itself. After screening seven substitutions, two were studied in site-specific detail: E15L (Glu→Leu near the turn, net stabilizing) and L32V (Leu→Val in the hydrophobic core, net destabilizing), each monitored at four 13C-labeled positions. Fitting the Ising-like model showed that adjusting only the mutated residue's parameters could not reproduce the data — all interaction parameters had to be re-optimized — and that both mutants retained the wild-type ranking of local stabilities, with all probed segments shifting up (E15L) or down (L32V) together. The folding free-energy surfaces, and hence the inferred folding pathway, were essentially unchanged. The authors concluded that the mutation effect is global rather than local, evidencing strong interconnectivity of the stabilizing interaction network even in this small, apparently non-cooperative motif, with the native structure dominating the folding mechanism.

Problem Statement

For the de novo designed helix-turn-helix motif αtα, determine whether the experimentally observed effect of a point mutation on site-specific stability and folding is genuinely global (position-independent) rather than local, and explain this behavior physically. The 2018 study of the E15L and L32V mutants concluded that mutations act globally and that the native structure dominates the folding mechanism; that conclusion rests on a coarse-grained Ising-like (DSA/L) model with a single contact energy ε\varepsilon and entropy cost Δs\Delta s per residue, on the unverified assumption that the mutant folded structures and contact maps are identical to wild type, on only two mutations, and on equilibrium data alone. The problem is to test and ground this conclusion; the following interconnected aspects all belong to the same question: (1) whether a more detailed description of the same mutants reproduces the global effect or reveals local perturbations that the single-parameter coarse-grained model cannot capture; (2) what the actual structures of the αtα mutants are, and how any side-chain rearrangements or contact-map changes propagate into the model predictions; (3) whether the global effect generalizes to other positions in αtα and to more disruptive substitutions; (4) whether the folding pathway itself is unchanged by the mutations; and (5) whether the mutation-induced changes of the coarse-grained parameters ε\varepsilon and Δs\Delta s have a physical-chemical basis in the substitution, so that the parametrization becomes transferable and predictive rather than fitted anew for each mutant. A resolution is a defensible answer to the question posed above — whether the effect of a point mutation is genuinely global or local, and the physical reason for it — supported by evidence beyond the original coarse-grained (DSA/L) fit

References

  1. Jason K. Lai, Ginka S. Kubelka, Jan Kubelka, 'Effect of Mutations on the Global and Site-Specific Stability and Folding of an Elementary Protein Structural Motif', J. Phys. Chem. B 122 (2018) 11083–11094, DOI: 10.1021/acs.jpcb.8b05280, https://doi.org/10.1021/acs.jpcb.8b05280
  2. Jason K. Lai, Ginka S. Kubelka, Jan Kubelka, 'Sequence, structure, and cooperativity in folding of elementary protein structural motifs', Proc. Natl. Acad. Sci. USA 112 (2015) 9890–9895, DOI: 10.1073/pnas.1506309112, https://doi.org/10.1073/pnas.1506309112
  3. Ginka S. Kubelka, Jan Kubelka, 'Site-specific thermodynamic stability and unfolding of a de novo designed protein structural motif mapped by 13C isotopically edited IR spectroscopy', J. Am. Chem. Soc. 136 (2014) 6037–6048, DOI: 10.1021/ja500918k, https://doi.org/10.1021/ja500918k
  4. Youcef Fezoui, Peter J. Connolly, John J. Osterhout, 'Solution structure of αtα, a helical hairpin peptide of de novo design', Protein Sci. 6 (1997) 1869–1877, DOI: 10.1002/pro.5560060907, https://doi.org/10.1002/pro.5560060907
  5. Juan Zeng, Fan Jiang, Yun-Dong Wu, 'Folding Simulations of an α-Helical Hairpin Motif αtα with Residue-Specific Force Fields', J. Phys. Chem. B 120 (2016) 33–41, DOI: 10.1021/acs.jpcb.5b09027, https://doi.org/10.1021/acs.jpcb.5b09027
  6. Koji Ooka, Runjing Liu, Munehito Arai, 'The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics', Molecules 27 (2022) 4460, DOI: 10.3390/molecules27144460, https://doi.org/10.3390/molecules27144460
  7. Koji Ooka, Munehito Arai, 'Accurate prediction of protein folding mechanisms by simple structure-based statistical mechanical models', Nat. Commun. 14 (2023) 6338, DOI: 10.1038/s41467-023-41664-1, https://doi.org/10.1038/s41467-023-41664-1
  8. Soundhararajan Gopi, Devanshu Devanshu, Nandakumar Rajasekaran, Sathvik Anantakrishnan, Athi N. Naganathan, 'pPerturb: A Server for Predicting Long-Distance Energetic Couplings and Mutation-Induced Stability Changes in Proteins via Perturbations', ACS Omega 5 (2020) 1142–1146, DOI: 10.1021/acsomega.9b03371, https://doi.org/10.1021/acsomega.9b03371
  9. Ginka S. Kubelka, Jan Kubelka, 'Multi-Probe Equilibrium Analysis of Gradual (Un)Folding Processes', Methods Mol. Biol. 2376 (2022) 161–171, DOI: 10.1007/978-1-0716-1716-8_9, https://doi.org/10.1007/978-1-0716-1716-8_9