CoChat
Sign Up
Cited by 573Protein Structure and DynamicsEnzyme Structure and FunctionReceptor Mechanisms and Signaling

The DOI matches its Crossref record; 90 of the 90 references it cites were found in a registry or an index; none of the 49 sources its claims cite is checked for support yet.

Predicting multiple conformations via sequence clustering and AlphaFold2

Open access
Wayment-Steele, Hannah K., Ojoawo, Adedolapo, Otten, Renee, Apitz, Julia M., Pitsawong, Warintra, Hömberger, Marc, Ovchinnikov, Sergey, Colwell, Lucy, Kern, Dorothee·2024·Nature
DOI 10.1038/s41586-023-06832-9PMC PMC10808063
Open this paper in CoChat — free account

Abstract

AlphaFold2 (ref. 1 ) has revolutionized structural biology by accurately predicting single structures of proteins. However, a protein’s biological function often depends on multiple conformational substates 2 , and disease-causing point mutations often cause population changes within these substates 3,4 . We demonstrate that clustering a multiple-sequence alignment by sequence similarity enables AlphaFold2 to sample alternative states of known metamorphic proteins with high confidence. Using this method, named AF-Cluster, we investigated the evolutionary distribution of predicted structures for the metamorphic protein KaiB 5 and found that predictions of both conformations were distributed in clusters across the KaiB family. We used nuclear magnetic resonance spectroscopy to confirm an AF-Cluster prediction: a cyanobacteria KaiB variant is stabilized in the opposite state compared with the more widely studied variant. To test AF-Cluster’s sensitivity to point mutations, we designed and experimentally verified a set of three mutations predicted to flip KaiB from Rhodobacter sphaeroides from the ground to the fold-switched state. Finally, screening for alternative states in protein families without known fold switching identified a putative alternative state for the oxidoreductase Mpt53 in Mycobacterium tuberculosis . Further development of such bioinformatic methods in tandem with experiments will probably have a considerable impact on predicting protein energy landscapes, essential for illuminating biological function.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

AF-Cluster clusters a protein multiple-sequence alignment by sequence similarity and feeds each cluster to AlphaFold2, enabling sampling of alternative conformational states for metamorphic proteins. The authors experimentally validate predicted KaiB states and designed point mutations, and identify a putative alternative state for Mpt53.

Key findings

  • AF-Cluster predicted both ground and fold-switched states for KaiB, RfaH and MAD2; for RfaH, top-five-model plDDT was 84.2 for the active state and 73.9 for the autoinhibited state, versus 68.6 for the autoinhibited state from the complete MSA.
  • NMR confirmed the AF-Cluster prediction for KaiBTV-4: CS-Rosetta models were within 1.8 ± 0.3 Å r.m.s.d. to the AF-Cluster model, and the NMR structure matched with 1.89 ± 0.13 Å r.m.s.d. and average pairwise r.m.s.d. of 0.7 Å over backbone atoms.
  • Three mutations (I68R, V83D, N84R) switched KaiBRS from the ground to the fold-switched state; NMR showed the mutant occupied the FS state at 86% versus 11% in wild type, and the ground state at 14% versus 89% in wild type.
  • Screening 628 protein families identified a putative alternative state for M. tuberculosis Mpt53; in the GA/GB system, AF-Cluster correctly predicted the most stable fold for 10 out of 12 point-mutation sets, whereas default AF2 predicted 8 out of 12.

Methods

  • MSAs were generated with ColabFold/MMseqs2, sequences with >25% gaps removed, then clustered by edit distance using DBSCAN; each cluster was used as input to AlphaFold2, with MSA Transformer used to probe evolutionary couplings.
  • Experimental validation used NMR spectroscopy, CS-Rosetta, NOESY-based structure calculation, SEC–MALS, and a curated KaiB phylogenetic tree of 487 variants.
  • The fold-switch screen used 628 proteins of 48–150 residues with >1,000 MSA sequences from a database of MSAs associated with crystal structures; Mpt53 homologues were analyzed with DALI and phylogenetic analysis.

Limitations

  • The paper states that pLDDT cannot be used as a measure of free energy; the thermodynamically disfavoured KaiB FS state still had higher pLDDT than the ground state, and the number of models per state reflects construct abundance rather than Boltzmann weight.
  • The paper states that AF-Cluster could not predict the oligomeric states of selecase, lymphotactin or CLIC1, and that conformational substates absent from the evolutionary signal may require alternative methods.
  • The paper notes that novel computationally identified folds may be difficult to validate experimentally if the triggering stimulus is unknown; the Mpt53 alternative state is a prediction without experimental validation in this work.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1038/s41586-023-06832-9

Cites

90 of 90 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Fig. 1 AF2 predictions from MSA clusters for the fold-switching protein KaiB return both known structures. a,b, Crystal structures of KaiB from T. elongatus (KaiBTE) in the ground state (PDB: 2QKE) (a) and the FS state (PDB: 5JYT) (b). c, The default ColabFold prediction of KaiBTE returns the FS state. Using only the closest 50 sequences by sequence distance returned from the MSA returns the ground state. For a–c, the first 50 residues that are identical in both states are coloured grey and the fold-switching elements are coloured the same in both states. d, Overview of the AF-Cluster method. Left, MSA is clustered by sequence similarity. Sequence space is depicted using a t-distributed stochastic neighbour embedding (t-SNE)56 of the one-hot sequence encoding. Right, clusters are used as an input to AF2, resulting in a distribution of predicted structures, coloured by plDDT. e, The top five models for the ground and FS state, ranked by plDDT. f, The r.m.s.d. of AF2 structure predictions for all clusters relative to the ground and FS state. The highest-confidence regions of the AF-Cluster distribution for KaiBTE are within 3 Å r.m.s.d. of crystal structures of both the ground and FS state. By contrast, sampling the MSA uniformly returns only the FS state with high confidence.
Fig. 1 AF2 predictions from MSA clusters for the fold-switching protein KaiB return both known structures. a,b, Crystal structures of KaiB from T. elongatus (KaiBTE) in the ground state (PDB: 2QKE) (a) and the FS state (PDB: 5JYT) (b). c, The default ColabFold prediction of KaiBTE returns the FS state. Using only the closest 50 sequences by sequence distance returned from the MSA returns the ground state. For a–c, the first 50 residues that are identical in both states are coloured grey and the fold-switching elements are coloured the same in both states. d, Overview of the AF-Cluster method. Left, MSA is clustered by sequence similarity. Sequence space is depicted using a t-distributed stochastic neighbour embedding (t-SNE)56 of the one-hot sequence encoding. Right, clusters are used as an input to AF2, resulting in a distribution of predicted structures, coloured by plDDT. e, The top five models for the ground and FS state, ranked by plDDT. f, The r.m.s.d. of AF2 structure predictions for all clusters relative to the ground and FS state. The highest-confidence regions of the AF-Cluster distribution for KaiBTE are within 3 Å r.m.s.d. of crystal structures of both the ground and FS state. By contrast, sampling the MSA uniformly returns only the FS state with high confidence.CC BY · PubMed Central
Fig. 2 The KaiB family contains pockets of sequences predicted to be stabilized for both states. a, AF2 predictions for each variant in a phylogenetic tree using the ten closest sequences as the input MSA. Left, each node is coloured by predicted state (blue, ground state; red, FS state). Right, the same tree, coloured by plDDT. b, Three known fold-switching KaiB variants from R. sphaeroides32 (i), T. elongatus5 (ii) and S. elongatus31 (iii) are predicted in the ground state, and a variant from L. pneumonia37 (iv), crystallized in the FS state, is predicted in the FS state with a high plDDT. c, A KaiB copy present in T. elongatus vestitus, KaiBTV-4, is predicted to favour the FS state. d–f, Experimental testing of KaiBTV-4. d, The secondary structure propensity determined by NMR backbone chemical shifts, calculated using TALOS-N57 for KaiBTV-4, fully agrees with the FS state predicted by AF-Cluster. Unassigned amino acid residues are indicated by stars. e, Structure models calculated using CS-Rosetta39, shown in grey, have 1.8 ± 0.3 Å r.m.s.d. to the AF-Cluster model (magenta). f, NMR structural models calculated from 3D 1H-15N- and 3D 1H-13C-edited NOESY spectra have an average pairwise r.m.s.d. of 0.7 Å, and 1.89 ± 0.13 Å r.m.s.d. to the AF-Cluster model. r.m.s.d. values in e and f were calculated over backbone atoms in secondary structure regions.
Fig. 2 The KaiB family contains pockets of sequences predicted to be stabilized for both states. a, AF2 predictions for each variant in a phylogenetic tree using the ten closest sequences as the input MSA. Left, each node is coloured by predicted state (blue, ground state; red, FS state). Right, the same tree, coloured by plDDT. b, Three known fold-switching KaiB variants from R. sphaeroides32 (i), T. elongatus5 (ii) and S. elongatus31 (iii) are predicted in the ground state, and a variant from L. pneumonia37 (iv), crystallized in the FS state, is predicted in the FS state with a high plDDT. c, A KaiB copy present in T. elongatus vestitus, KaiBTV-4, is predicted to favour the FS state. d–f, Experimental testing of KaiBTV-4. d, The secondary structure propensity determined by NMR backbone chemical shifts, calculated using TALOS-N57 for KaiBTV-4, fully agrees with the FS state predicted by AF-Cluster. Unassigned amino acid residues are indicated by stars. e, Structure models calculated using CS-Rosetta39, shown in grey, have 1.8 ± 0.3 Å r.m.s.d. to the AF-Cluster model (magenta). f, NMR structural models calculated from 3D 1H-15N- and 3D 1H-13C-edited NOESY spectra have an average pairwise r.m.s.d. of 0.7 Å, and 1.89 ± 0.13 Å r.m.s.d. to the AF-Cluster model. r.m.s.d. values in e and f were calculated over backbone atoms in secondary structure regions.CC BY · PubMed Central
Fig. 3 A designed minimal set of mutations switches the predominant fold of KaiBRS from the ground state to the FS state. a, Sequence features enriched in clusters that predict the FS and ground state. b, Three mutations are sufficient to switch the structure prediction for KaiBRS in AF2 from the ground state to the FS state. Top, AF-Cluster models for KaiBRS and KaiBRS-3m, coloured by plDDT. Bottom, three mutation sites are highlighted. c, Overlaid 1H-15N HSQC spectra of KaiBRS (blue) and KaiBRS-3m (red). d, Examples of residues from well-resolved regions in the 1H-15N HSQC assigned in both states are shown for WT KaiBRS and KaiBRS-3m to illustrate the flip in populations through the three mutations. e, Chemical-shift-based secondary structure calculated using TALOS-N57 analysis of the ground and FS states of KaiBRS and the major state of KaiBRS-3m. Unassigned amino acid residues are indicated by stars. The green box indicates the fold-switching region. f, Average of the NMR peak intensity ratio of ground versus FS state for select residues that could be assigned in both states for both variants in well-resolved regions. The error bars represent the s.e.m. n = 5 residues.
Fig. 3 A designed minimal set of mutations switches the predominant fold of KaiBRS from the ground state to the FS state. a, Sequence features enriched in clusters that predict the FS and ground state. b, Three mutations are sufficient to switch the structure prediction for KaiBRS in AF2 from the ground state to the FS state. Top, AF-Cluster models for KaiBRS and KaiBRS-3m, coloured by plDDT. Bottom, three mutation sites are highlighted. c, Overlaid 1H-15N HSQC spectra of KaiBRS (blue) and KaiBRS-3m (red). d, Examples of residues from well-resolved regions in the 1H-15N HSQC assigned in both states are shown for WT KaiBRS and KaiBRS-3m to illustrate the flip in populations through the three mutations. e, Chemical-shift-based secondary structure calculated using TALOS-N57 analysis of the ground and FS states of KaiBRS and the major state of KaiBRS-3m. Unassigned amino acid residues are indicated by stars. The green box indicates the fold-switching region. f, Average of the NMR peak intensity ratio of ground versus FS state for select residues that could be assigned in both states for both variants in well-resolved regions. The error bars represent the s.e.m. n = 5 residues.CC BY · PubMed Central
Fig. 4 AF-Cluster predicts fold switching for the proteins RfaH and MAD2. a, Fold switching in the RfaH transcription factor in E. coli. In RfaH’s autoinhibited state, the CTD (red) forms an α-helix bundle (PDB: 5OND)58. In the active state, the CTD unbinds and forms a β-sheet that is homologous to the transcription factor NusG (CTD PDB: 2LCL)22. b, AF-Cluster returns structure models that include both the autoinhibited and the active state with high confidence. Note that the CTD orientation is not defined due to the flexible linker between the two domains. c, The closed state (PDB: 1S2H)27 and the open state (PDB: 1DUJ)59 of the MAD2 spindle checkpoint in humans with the fold-switching portions coloured. d, Both MAD2 states are predicted by AF-Cluster with high confidence.
Fig. 4 AF-Cluster predicts fold switching for the proteins RfaH and MAD2. a, Fold switching in the RfaH transcription factor in E. coli. In RfaH’s autoinhibited state, the CTD (red) forms an α-helix bundle (PDB: 5OND)58. In the active state, the CTD unbinds and forms a β-sheet that is homologous to the transcription factor NusG (CTD PDB: 2LCL)22. b, AF-Cluster returns structure models that include both the autoinhibited and the active state with high confidence. Note that the CTD orientation is not defined due to the flexible linker between the two domains. c, The closed state (PDB: 1S2H)27 and the open state (PDB: 1DUJ)59 of the MAD2 spindle checkpoint in humans with the fold-switching portions coloured. d, Both MAD2 states are predicted by AF-Cluster with high confidence.CC BY · PubMed Central
Fig. 5 Screening for fold switching in many protein families predicts a putative alternative fold for the M. tuberculosis secreted protein Mpt53. a, Overview of the strategy for detecting novel predicted alternative folds. Screening of 628 families with more than 1,000 sequences in their MSA and residue length 48–150 from ref. 30. After clustering, we ran AF2 predictions using ten randomly selected clusters from each. b, Candidates for further sampling were selected by looking for outlier predictions with a high r.m.s.d. to the reference structure and high plDDT. c, Sampled models for candidate Mpt53, visualized using PCA of the closest heavy-atom contacts. Two states with a higher plDDT than the background were observed. d, The top five models by plDDT for the known state (top) and the putative alternative state (bottom), coloured by plDDT per residue. e, The crystal structure of the reduced state of M. tuberculosis Mpt53 (PDB: 1LU4), which corresponds to state 1 in the sampled landscape (top). In the putative alternative state 2, strand β1 replaces β5 in the five-strand β-sheet. Helix α4 shifts to the other side of the β-sheet and helix α5 is displaced.
Fig. 5 Screening for fold switching in many protein families predicts a putative alternative fold for the M. tuberculosis secreted protein Mpt53. a, Overview of the strategy for detecting novel predicted alternative folds. Screening of 628 families with more than 1,000 sequences in their MSA and residue length 48–150 from ref. 30. After clustering, we ran AF2 predictions using ten randomly selected clusters from each. b, Candidates for further sampling were selected by looking for outlier predictions with a high r.m.s.d. to the reference structure and high plDDT. c, Sampled models for candidate Mpt53, visualized using PCA of the closest heavy-atom contacts. Two states with a higher plDDT than the background were observed. d, The top five models by plDDT for the known state (top) and the putative alternative state (bottom), coloured by plDDT per residue. e, The crystal structure of the reduced state of M. tuberculosis Mpt53 (PDB: 1LU4), which corresponds to state 1 in the sampled landscape (top). In the putative alternative state 2, strand β1 replaces β5 in the five-strand β-sheet. Helix α4 shifts to the other side of the β-sheet and helix α5 is displaced.CC BY · PubMed Central
Extended Data Fig. 1 Investigating two highly-similar sets of sequences in the KaiBTE MSA. a) Sequences 1-50 predict the ground state in all 5 AF2 models, whereas sequences 50–100 predict the FS state in 4 of 5 models. Sequences are ranked by sequence similarity from the ColabFold MSA generation routine. b) MSA Transformer predicted contacts for both sets of sequences. c) Taking the difference of both contact maps highlights that sequences 50–100 contain features for the FS state corresponding to beta-strands (boxed in orange, magenta) and the helix-helix interaction (boxed in red). Right: Structure model (PDB: 5JYT) for the FS state of KaiBTE, features coloured analogously.
Extended Data Fig. 1 Investigating two highly-similar sets of sequences in the KaiBTE MSA. a) Sequences 1-50 predict the ground state in all 5 AF2 models, whereas sequences 50–100 predict the FS state in 4 of 5 models. Sequences are ranked by sequence similarity from the ColabFold MSA generation routine. b) MSA Transformer predicted contacts for both sets of sequences. c) Taking the difference of both contact maps highlights that sequences 50–100 contain features for the FS state corresponding to beta-strands (boxed in orange, magenta) and the helix-helix interaction (boxed in red). Right: Structure model (PDB: 5JYT) for the FS state of KaiBTE, features coloured analogously.CC BY · PubMed Central
Extended Data Fig. 2 Empirically maximizing information content of clustering using DBSCAN35. a) Varying the parameter epsilon, which controls the maximum allowable distance for points to be in a cluster, results in a peak in the number of clusters DBSCAN identifies for a set of sequences. b) For epsilon <epsmax, fewer sequences are clustered, i.e. more are identified as outliers by the DBSCAN algorithm. For epsilon > epsmax, more sequences are clustered but fewer clusters are returned as more clusters are joined. c) Example clusterings of KaiB sequences at different epsilon values (compare to Fig. 1d). d) Corresponding KaiB landscape of predictions for these epsilon values. e) The plDDT values of models within 3 Å RMSD of the ground and FS state from the clustered sampling method are statistically significantly higher than the rest of the models. Box plots depict median and 25/75% interquartile range, whiskers = 1.5 *interquartile range. P-values for sample comparisons with p < 0.05 indicated, calculated via a two-sided test for the null hypothesis that 2 independent samples have identical mean values. n = 500 models for the two Uniform sampling methods, n = 230 for AF-Cluster sampling.
Extended Data Fig. 2 Empirically maximizing information content of clustering using DBSCAN35. a) Varying the parameter epsilon, which controls the maximum allowable distance for points to be in a cluster, results in a peak in the number of clusters DBSCAN identifies for a set of sequences. b) For epsilon <epsmax, fewer sequences are clustered, i.e. more are identified as outliers by the DBSCAN algorithm. For epsilon > epsmax, more sequences are clustered but fewer clusters are returned as more clusters are joined. c) Example clusterings of KaiB sequences at different epsilon values (compare to Fig. 1d). d) Corresponding KaiB landscape of predictions for these epsilon values. e) The plDDT values of models within 3 Å RMSD of the ground and FS state from the clustered sampling method are statistically significantly higher than the rest of the models. Box plots depict median and 25/75% interquartile range, whiskers = 1.5 *interquartile range. P-values for sample comparisons with p < 0.05 indicated, calculated via a two-sided test for the null hypothesis that 2 independent samples have identical mean values. n = 500 models for the two Uniform sampling methods, n = 230 for AF-Cluster sampling.CC BY · PubMed Central
Extended Data Fig. 3 AF-Cluster sampling detects KaiBTE ground state in evolutionary couplings. a) Contacts under 5 Å that correspond uniquely to the FS state (left) or ground state (right). Boxed features correspond to features unique to both states. b) AUC scores to both states of contact maps predicted by MSA transformer, a method trained by unsupervised learning. Randomly-subsampled MSAs have higher score to the FS state, and AF-Cluster contacts have higher score to ground state. c) Contact maps of sampled MSAs with the highest AUC to the ground state from (i) AF-Cluster and (ii) random sampling. The best-scoring random sample does not include the beta-strand unique to the ground state (boxed in blue in i). (iii) Contacts calculated from the whole MSA show features corresponding the FS state: beta-strands (orange, magenta) and the helix-helix interaction (red) boxed in (A). d) Contact map scores for both states correlate to the AF2 prediction RMSD for each state (Ground state: Spearman R = −0.32, p = 2e-09, FS state: Spearman R = −0.34, p = 4e-10 via a two-sided statistical test. No adjustment for multiple comparisons was made). Error bands for the linear trendline are 95% confidence intervals obtained via bootstrapping.
Extended Data Fig. 3 AF-Cluster sampling detects KaiBTE ground state in evolutionary couplings. a) Contacts under 5 Å that correspond uniquely to the FS state (left) or ground state (right). Boxed features correspond to features unique to both states. b) AUC scores to both states of contact maps predicted by MSA transformer, a method trained by unsupervised learning. Randomly-subsampled MSAs have higher score to the FS state, and AF-Cluster contacts have higher score to ground state. c) Contact maps of sampled MSAs with the highest AUC to the ground state from (i) AF-Cluster and (ii) random sampling. The best-scoring random sample does not include the beta-strand unique to the ground state (boxed in blue in i). (iii) Contacts calculated from the whole MSA show features corresponding the FS state: beta-strands (orange, magenta) and the helix-helix interaction (red) boxed in (A). d) Contact map scores for both states correlate to the AF2 prediction RMSD for each state (Ground state: Spearman R = −0.32, p = 2e-09, FS state: Spearman R = −0.34, p = 4e-10 via a two-sided statistical test. No adjustment for multiple comparisons was made). Error bands for the linear trendline are 95% confidence intervals obtained via bootstrapping.CC BY · PubMed Central
Extended Data Fig. 4 Supplemental experimental data for KaiBTV-4 and KaiBRS-3m. a) 1H 15N HSQC spectra of KaiBTV-4 indicates one major folded state. Assignments are shown. b) Strip plot extracted from a 150 ms mixing time 15N-edited NOESY-HSQC spectrum of KaiBTV-4 illustrating the inter-strand NOEs between residues V3-V8; T35-D40; T58-Y63; R68-Y71, used in confirming KaiBTV-4 is in the fold-switched state. c) Summary of NOEs between the parallel β-sheets V3-V8 and T35-D40, and the antiparallel β-sheets T58-V62 and R68-V71. Confirmed NOEs are depicted by dashed lines. NOEs not depicted could not be confirmed unambiguously. SEC-MALS analysis of (d) KaiBTV-4 and (e) KaiBRS-3m at NMR concentration of 500 μM indicate both are monomeric. The profiles on the left show the full SEC-MALS run with the light scattering (LS) profile in blue, normalized UV profile in red, and refractive index (RI) profile in green. On the right is the region of the peak of interest showing the light scattering profile (blue) plotted against elution time, and the protein molar masses are indicated in red. The molar masses of KaiBTV−4 and KaiBRS-3m have been determined from light scattering and refractometry data to be (9.5 +/− 3.0) kDa and (9.4 +/− 1.7) kDa, respectively.
Extended Data Fig. 4 Supplemental experimental data for KaiBTV-4 and KaiBRS-3m. a) 1H 15N HSQC spectra of KaiBTV-4 indicates one major folded state. Assignments are shown. b) Strip plot extracted from a 150 ms mixing time 15N-edited NOESY-HSQC spectrum of KaiBTV-4 illustrating the inter-strand NOEs between residues V3-V8; T35-D40; T58-Y63; R68-Y71, used in confirming KaiBTV-4 is in the fold-switched state. c) Summary of NOEs between the parallel β-sheets V3-V8 and T35-D40, and the antiparallel β-sheets T58-V62 and R68-V71. Confirmed NOEs are depicted by dashed lines. NOEs not depicted could not be confirmed unambiguously. SEC-MALS analysis of (d) KaiBTV-4 and (e) KaiBRS-3m at NMR concentration of 500 μM indicate both are monomeric. The profiles on the left show the full SEC-MALS run with the light scattering (LS) profile in blue, normalized UV profile in red, and refractive index (RI) profile in green. On the right is the region of the peak of interest showing the light scattering profile (blue) plotted against elution time, and the protein molar masses are indicated in red. The molar masses of KaiBTV−4 and KaiBRS-3m have been determined from light scattering and refractometry data to be (9.5 +/− 3.0) kDa and (9.4 +/− 1.7) kDa, respectively.CC BY · PubMed Central
Extended Data Fig. 5 Three mutations are sufficient to switch KaiBRS AF2 prediction to high-confidence FS state prediction. a) plDDT from AF2 (no MSA, 12 recycles, model 1) for all combinations of 8 possible point mutations most enriched from FS state analysis (cf. Fig. 3b). Quadruple-mutants and greater are not labelled by residue mutation, as we searched for the minimal set of mutations to flip the conformational equilibrium. b) Structure models of single mutant V83D, double mutant V83D-I68R, and triple mutant V83D-I68R-N84A demonstrating that V83D switches the C-terminal strand to a helix, and I68R switches the C-terminal helices to a strand. N84A increases the plDDT of the prediction of the FS state. Top row: structures coloured as in Fig. 1a. Bottom row: structures coloured by plDDT.
Extended Data Fig. 5 Three mutations are sufficient to switch KaiBRS AF2 prediction to high-confidence FS state prediction. a) plDDT from AF2 (no MSA, 12 recycles, model 1) for all combinations of 8 possible point mutations most enriched from FS state analysis (cf. Fig. 3b). Quadruple-mutants and greater are not labelled by residue mutation, as we searched for the minimal set of mutations to flip the conformational equilibrium. b) Structure models of single mutant V83D, double mutant V83D-I68R, and triple mutant V83D-I68R-N84A demonstrating that V83D switches the C-terminal strand to a helix, and I68R switches the C-terminal helices to a strand. N84A increases the plDDT of the prediction of the FS state. Top row: structures coloured as in Fig. 1a. Bottom row: structures coloured by plDDT.CC BY · PubMed Central
Extended Data Fig. 6 Results corresponding to testing AF-Cluster for other proteins. a) Predicting the structure of RfaH in AF2 with the complete MSA from ColabFold34 returns the autoinhibited state with a mean plDDT of 68.6 (note low confidence in the first alpha-helix of the CTD.) b) B-factors of PDB 5OND58, indicating that the last helical turn of the second to last helix has high B-factors (arrow). C) AF-Cluster only predicts the monomeric state for proteins that switch between monomeric and oligomeric states.
Extended Data Fig. 6 Results corresponding to testing AF-Cluster for other proteins. a) Predicting the structure of RfaH in AF2 with the complete MSA from ColabFold34 returns the autoinhibited state with a mean plDDT of 68.6 (note low confidence in the first alpha-helix of the CTD.) b) B-factors of PDB 5OND58, indicating that the last helical turn of the second to last helix has high B-factors (arrow). C) AF-Cluster only predicts the monomeric state for proteins that switch between monomeric and oligomeric states.CC BY · PubMed Central
Extended Data Fig. 7 Investigating the source of the AF-Cluster prediction for an alternate state of Mpt53. a) plDDT vs. RMSD for AF-Cluster sampling on oxidoreductase Mpt53. Each prediction coloured by MSA size. b) plDDT values for state 1, corresponding to the known thioredoxin-like state, and an alternate unknown state are significantly higher than background. Box plots depict median and 25/75% interquartile range, whiskers = 1.5 *interquartile range. P-values for sample comparisons with p < 0.05 indicated, calculated via a two-sided test for the null hypothesis that 2 independent samples have identical mean values. n = 1642 models total. c) The conserved CxxC active site is very similar between its conformation in the crystal structure and models for the putative alternate state. d) Workflow for using DALI91 to screen for structure homologues to both Mpt53’s original state and predicted alternate state to search for any similar structures in the PDB that might have been in AF2’s training set. e) Plotting RMSD normalized by alignment length to both structures reveals some structures with lower weighted RMSD to the alternate state than to the original state. f) 7 of 9 DALI hits with lower alternate state RMSD contained an alpha-helix positioned in similar same way as in the Mpt53 alternate state (coloured in green). One structure (3EMX) also contained an N-terminus beta-strand positioned similarly to the alternate state.
Extended Data Fig. 7 Investigating the source of the AF-Cluster prediction for an alternate state of Mpt53. a) plDDT vs. RMSD for AF-Cluster sampling on oxidoreductase Mpt53. Each prediction coloured by MSA size. b) plDDT values for state 1, corresponding to the known thioredoxin-like state, and an alternate unknown state are significantly higher than background. Box plots depict median and 25/75% interquartile range, whiskers = 1.5 *interquartile range. P-values for sample comparisons with p < 0.05 indicated, calculated via a two-sided test for the null hypothesis that 2 independent samples have identical mean values. n = 1642 models total. c) The conserved CxxC active site is very similar between its conformation in the crystal structure and models for the putative alternate state. d) Workflow for using DALI91 to screen for structure homologues to both Mpt53’s original state and predicted alternate state to search for any similar structures in the PDB that might have been in AF2’s training set. e) Plotting RMSD normalized by alignment length to both structures reveals some structures with lower weighted RMSD to the alternate state than to the original state. f) 7 of 9 DALI hits with lower alternate state RMSD contained an alpha-helix positioned in similar same way as in the Mpt53 alternate state (coloured in green). One structure (3EMX) also contained an N-terminus beta-strand positioned similarly to the alternate state.CC BY · PubMed Central
Extended Data Fig. 8 An analogous fold-switch state is predicted for some Mpt53 structure homologues. 6 of the 10 screened homologues from DALI91 with the lowest RMSD to the original state predicted an alternate state similar to that of Mpt53. a) Conformational landscapes, visualized by RMSD to two states of Mpt53, and showing the corresponding known structures (above) and predicted alternate structures (below), coloured analogously to Mpt53 (cf. Fig. 5e). b) Alternate structures in (a), coloured by plDDT. c) Conformational landscapes of 4 structure homologues with no evidence for predicted alternate state.
Extended Data Fig. 8 An analogous fold-switch state is predicted for some Mpt53 structure homologues. 6 of the 10 screened homologues from DALI91 with the lowest RMSD to the original state predicted an alternate state similar to that of Mpt53. a) Conformational landscapes, visualized by RMSD to two states of Mpt53, and showing the corresponding known structures (above) and predicted alternate structures (below), coloured analogously to Mpt53 (cf. Fig. 5e). b) Alternate structures in (a), coloured by plDDT. c) Conformational landscapes of 4 structure homologues with no evidence for predicted alternate state.CC BY · PubMed Central
Extended Data Fig. 9 Phylogenetic tree of closest structure matches for Mpt53 states. Homologues for Mpt53 original state and alternate state are dispersed across a calculated phylogenetic tree of the structure hits for both identified via DALI91 (cf. Extended Data Fig. 7).
Extended Data Fig. 9 Phylogenetic tree of closest structure matches for Mpt53 states. Homologues for Mpt53 original state and alternate state are dispersed across a calculated phylogenetic tree of the structure hits for both identified via DALI91 (cf. Extended Data Fig. 7).CC BY · PubMed Central
Extended Data Fig. 10 MSA clusters enable correct predictions for engineered fold-switching point mutations in the protein GA/GB system. A) Sequences of the 12 sets of GA/GB point mutations tested, from refs. 49–51 Point mutations different from neighbouring sequences in the series are coloured in orange. Right: Representative NMR structures of the GA and GB fold. B) Left: Visualization of sequence identity and coverage of the MSA returned by ColabFold for GA98. Right: Visualization of MSA clusters with more than 10 sequences from the AF-Cluster clustering routine. C) We compared 3 types of MSAs for each point mutation: i) the full MSA returned by ColabFold, ii) MSA clusters returned by AF-Cluster, and iii) MSAs of the wild-type GA and GB proteins in ref. 50 Predictions for which the highest plDDT is incorrect are marked with an X. AF-Cluster has a higher success rate and returns predictions with higher plDDT.
Extended Data Fig. 10 MSA clusters enable correct predictions for engineered fold-switching point mutations in the protein GA/GB system. A) Sequences of the 12 sets of GA/GB point mutations tested, from refs. 49–51 Point mutations different from neighbouring sequences in the series are coloured in orange. Right: Representative NMR structures of the GA and GB fold. B) Left: Visualization of sequence identity and coverage of the MSA returned by ColabFold for GA98. Right: Visualization of MSA clusters with more than 10 sequences from the AF-Cluster clustering routine. C) We compared 3 types of MSAs for each point mutation: i) the full MSA returned by ColabFold, ii) MSA clusters returned by AF-Cluster, and iii) MSAs of the wild-type GA and GB proteins in ref. 50 Predictions for which the highest plDDT is incorrect are marked with an X. AF-Cluster has a higher success rate and returns predictions with higher plDDT.CC BY · PubMed Central

Verification

Claims the paper makes

60 Own result · 32 Not checked yet · 10 Own interpretation · 3 Own method · 3 No basis given

Verify all 49 references — free account
  • It may state two separate claims as one. 33
  • Probably not a substantial claim. 21
  • The quote only partly states this claim. 1
50 of 50

  • The paper relies on this without support.

    Proteins must evolve in the context of the multiple conformational states that they adopt.

    “However, proteins must evolve in the context of the multiple conformational states that they adopt.”

  • The paper relies on this without support.

    There is a need for methods that deconvolve signal from multiple states if they are not already both present at the level of the entire MSA.

    “However, there is a need for methods that deconvolve signal from multiple states if they are not already both present at the level of the entire MSA.”
2 of 2

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

2 of 2

Citations

Cites

Show all 90 references

Research works printed but not linked: 3

  1. Ester, M., Kriegel, H.-P., Sander, J. & Xu, X. A density-based algorithm for discovering clusters in large spatial databases with noise. In Proc. 2nd International Conference on Knowledge Discovery and Data Mining (KDD-96) 226–231 (ACM, 1996).
  2. DeLano, W. L. Pymol: an open-source molecular graphics tool. CCP4 Newsl. Protein Crystallogr. 40, 82–92 (2002).
  3. Pedregosa, F. et al. Scikit-learn: machine learning in Python. J. Mach. Learn. Res. 12, 2825–2830 (2011).

Cited by

Sorted by citations

See all 573

Related

See all

Metadata

Authors
Wayment-Steele, H. K., Ojoawo, A., Otten, R., Apitz, J. M., Pitsawong, W., Hömberger, M., Ovchinnikov, S., Colwell, L., & Kern, D.
Affiliations
  • Howard Hughes Medical Institute (US)
  • Harvard University (US)
  • Google (United States) (US)
Corresponding author
Kern, Dorothee
ORCID
Venue
Nature
Publisher
Springer Science and Business Media LLC
Published
2023-11-13
Volume
625
Issue
7996
Pages
832-839
ISSN
  • 0028-0836 (ISSN-L)
  • 1476-4687
DOI
10.1038/s41586-023-06832-9
PMC
PMC10808063
OpenAlex
W4388627063
Semantic Scholar
6d464c44cad5eb182ae580adfff8624562e457b8
Cited by
545
References
90
Citations by year
2026: 175 · 2025: 202 · 2024: 160 · 2023: 8
FWCI
48.74
Type
journal article
Language
en
Open access
Licence
Locations
  • Nature (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • RePEc: Research Papers in Economics (submittedVersion)
  • Apollo (University of Cambridge) (submittedVersion, cc-by)
Funders
National Institute of General Medical Sciences (P41GM111135, T32GM135126, R24GM141526)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Protein Structure and Dynamics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Structural Biology › Enzyme Structure and Function
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Receptor Mechanisms and Signaling
Keywords
protein conformational switching, sequence clustering, AlphaFold2, KaiB, fold switching, point mutations, protein energy landscape, oxidative enzymes, protein family evolution
MeSH terms
Machine Learning, Antigens, Bacterial, Bacterial Proteins, Mutation, Protein Conformation, Proteins, Rhodobacter sphaeroides, Cluster Analysis, Sequence Alignment, Protein Folding

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-07, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-07 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review