User:Eric Martz/AlphaFold3 case studies
From Proteopedia
| Line 2: | Line 2: | ||
This page is under construction. This notice will be removed when it is completed. [[User:Eric Martz|Eric Martz]] 02:03, 13 November 2024 (UTC) | This page is under construction. This notice will be removed when it is completed. [[User:Eric Martz|Eric Martz]] 02:03, 13 November 2024 (UTC) | ||
</td></tr></table> | </td></tr></table> | ||
| - | The following case studies of predictions by the [https://alphafoldserver.com AlphaFold3 Server] were done in November 2024. At that time, the [https://alphafold.ebi.ac.uk/ AlphaFold Database] had been generated with AlphaFold2. Superpositions and RMSD values were obtained from [https://fatcat.godziklab.org/ FATCAT]. Rocking animations were generated by [http://firstglance.jmol.org FirstGlance in Jmol], and morphs were [http://firstglance.jmol.org/videocapture.htm captured] from FATCAT. [[pLDDT]] averages for specified sequence ranges | + | The following case studies of predictions by the [https://alphafoldserver.com AlphaFold3 Server] were done in November 2024. At that time, the [https://alphafold.ebi.ac.uk/ AlphaFold Database] had been generated with AlphaFold2. Superpositions and RMSD values were obtained from [https://fatcat.godziklab.org/ FATCAT]. Rocking animations were generated by [http://firstglance.jmol.org FirstGlance in Jmol], and morphs were [http://firstglance.jmol.org/videocapture.htm captured] from FATCAT. [[pLDDT]] averages for specified sequence ranges [[FirstGlance/How to get average pLDDT from AlphaFold models|were obtained using FirstGlance in Jmol]]. |
==AlphaFold3 Example 8AW3== | ==AlphaFold3 Example 8AW3== | ||
Revision as of 21:08, 16 November 2024
|
This page is under construction. This notice will be removed when it is completed. Eric Martz 02:03, 13 November 2024 (UTC) |
The following case studies of predictions by the AlphaFold3 Server were done in November 2024. At that time, the AlphaFold Database had been generated with AlphaFold2. Superpositions and RMSD values were obtained from FATCAT. Rocking animations were generated by FirstGlance in Jmol, and morphs were captured from FATCAT. pLDDT averages for specified sequence ranges were obtained using FirstGlance in Jmol.
Contents |
AlphaFold3 Example 8AW3
Summary: This is one of three examples featured in the AlphaFold3 Server. AlphaFold3 predicted the fold of the mid-341 residues of tRNA A34 deaminase (UniProt Q381Q7) accurately, when compared to the empirical model in 8aw3. It made no confident predictions for what appear to be 90 residues in two disordered loops, and a 25-residue disordered C terminus. In contrast, the prediction from the AlphaFold2 Database was bent relative to the empirical model, but otherwise largely accurate.
UniProt Q381Q7 (tRNA A34 deaminase) is the longest chain in the example 8aw3 provided by the AlphaFold3 Server. In 8aw3, its chain name is "3". 8aw3 is a cryo-EM structure with resolution 3.6 Å. 8aw3 indicates that the experimental material included full-length chain 3 of 369 residues, but 28 end-residues lack coordinates, so the coordinates run from UniProt 4-344, length 341. Within this range, there are two loops lacking coordinates of lengths 49 and 13, leaving 279 amino acids with coordinates. Neither missing loop is predicted to be disordered at UniProt, although the estimated disorder propensity is high for both.
Q381Q7: AlphaFold3 vs Empirical Structure
The overall 4-chain prediction (for the complex in 8aw3) has pTM of 0.75, at the low end of the range 0.7-0.9 deemed "confident" by the server. Both of the missing loops have "no-confidence" (red) coordinates. The longer one is modeled as a large circle protruding from the surface of the protein, typical for loops with intrinsic disorder. The average pLDDT for chain 3 is 68 (weak confidence), but with the large disordered loop removed, this increases to 74[1].
The AlphaFold3 prediction superposes (FATCAT rigid) with all 279 alpha carbons of chain 3 of 8aw3 with RMSD 2.2 Å. The morph shows the largest discrepancy in the red (no confidence) loop at center bottom, which has sequence SNSGCRKSNR (238-247 in the numbering of both UniProt and 8aw3). These residues have coordinates in 8aw3, but 3 residues, CRK, have incomplete sidechains.
|
|
|
|
AlphaFold3 prediction for Q38107. |
FATCAT morph between AlphaFold3 prediction and chain 3 of 8aw3. |
Q381Q7: AlphaFold3 vs AlphaFold2 Database
|
The long disordered loop was in different positions in the AlphaFold3 vs. AlphaFold2 Database models. When it was included, FATCAT superposed it between models, causing the RMSD to be 4.8-4.9 Å (rigid or flexible superposition). Deleting that loop from the models resulted in a much better flexible superposition RMSD value of 2.0 Å for 312 alpha carbons (98% of the 317 present). This close superposition required flexibility. The rigid FATCAT superposition still had RMSD 4.5 Å. The morph shows bending and twisting between the two models. Since the AlphaFold2 Database prediction is bent/twisted relative to the AlphaFold3 prediction, and since the latter superposes well with the empirical model in a rigid superposition, the AlphaFold2 Database prediction should not superpose as well with the empirical model in a rigid superposition. Indeed, the rigid superposition has RMSD 4.1 Å for 270 (of 279) alpha carbons (morph not shown). (Flexible superposition achieved RMSD 2.4 Å.) |
|
|
FATCAT morph between AlphaFold3 prediction and AlphaFold2 Database prediction for Q381Q7. |
8JRP Chain A: Partially Untemplated
|
AlphaFold3 Server FAQ says that when empirical templates exist in the wwPDB, they will be used unconditionally, adding that “… the server searches PDB for template structures with a cutoff date of 30th September 2021 …”. According to this, the N-terminal half of chain A of 8jrp[2] (120-519) should be untemplated, since a sequence search at RCSB finds no matching entries prior to 2023. Being untemplated should make structure prediction more challenging. Of course, a structure need not be sequence-identical to serve as a template. A search for similar chain structures finds no good matches (top hits "not significantly similar" according to FATCAT), consistent with this N-terminal half being untemplated. The C-terminal half (520-869) is templated by 2 entries, e.g. 1c4z[3], a 1999 X-ray 2.6 Å resolution structure. |
|
|
Chain A of 8jrp: putatively untemplated half, templated half. |
Missing/disordered
N-terminal residues 1-119, missing in 8jrp, are predicted to be partially disordered. Loop 170-230, missing from 8jrp, is predicted to be disordered at RCSB. Sequence range 388-424 is predicted to be disordered, but has coordinates in 8jrp.
|
|
RCSB sequence graphic (simplified) for chain A of 8jrp. Gray bars are missing residues. *Brown bars are residues missing atoms in their sidechains. |
AlphaFold 3 prediction
Chain A in 8jrp is human ubiquitin-protein ligase E3A. The full-length sequence, Q05086, has 875 amino acids. The 3.6 Å resolution cryo-EM structure indicates that the full-length protein was imaged, but the N-terminal end of the model lacks coordinates for 119 residues. The ends of the model are sequence numbers 120 and 869 (length 750), but one loop of length 61 (170-230) is missing coordinates, leaving 689 amino acids with coordinates. Conveniently, the sequence numbers in UniProt and 8jrp are the same.
|
|
|
|
|
Chain A of 8jrp. Missing ends and loops shown as "empty baskets" by FirstGlance in Jmol.
|
AlphaFold3 prediction for Q05086[4]. pLDDT > 70 is confident.
|
Morph of FATCAT rigid superposition between AlphaFold3 prediction and chain A of 8jrp, RMSD 2.1 Å. |
AlphaFold3 was unable to predict with confidence the loop 170-230 (length 61) missing in 8jrp, and predicted to be disordered. The average pLDDT for the predicted loop was 40.2 (no confidence). Similarly, the N-terminal residues 1-119 missing from 8jrp were predicted with average pLDDT of 55.5 (no confidence). The 27-residue loop 397-423, predicted to be disordered but with coordinates in 8jrp, appears red in the predicted model with average pLDDT 47.1 (no confidence). In contrast, the average pLDDT for the residues having coordinates in 8jrp (120-169, 231-869) is 86.1 (high confidence).
The AlphaFold3 prediction superposed 688/689 alpha carbons onto chain A of 8jrp with RMSD 2.1 Å. Thus, the prediction is very close to the empirical model.
See Also
- AlphaFold3 Server alphafoldserver.com
- How to predict structures with AlphaFold
- AlphaFold/Index lists pages in Proteopedia that relate to AlphaFold.
Notes
- ↑ Average pLDDT values were obtained with FirstGlance in Jmol, which reports the minimum, average, and maximum in the upper right as "Reliability".
- ↑ 8jrp was chosen despite its poor resolution of 3.6 Å because the better-resolution 8jrn, 2.6 Å, is not available in PDB format. A rigid superposition by FATCAT matched all 689 residues with coordinates, giving RMSD 0.7 Å, a near perfect match.
- ↑ A FATCAT rigid superposition of 8jrp chain A with 1cz4 chain C shows a twist between the two domains. A FATCAT flexible superposition with one twist gives RMSD 1.3 Å for all 350 alpha carbons. Thus, 1c4z would be an excellent template for predicting the C-terminal half of 8jrp chain A.
- ↑ Predictions were done for the full length Q05086 sequence (shown above), and with 1-119 removed. Rigid FATCAT superposition between the two predictions (only the 689 alpha carbons present in 8jrp) gave RMSD 1.4 Å. Thus, inclusion of 1-119 had little if any effect on the prediction of the remainder of the structure.
