.. _`getInitialModel`: Get the initial model with AlphaFold2 ===================================== In this section of the AlphaFold tutorial we cover the first steps in the small workflow shown before (:numref:`model_building_fig16_workflow`) addressed to get the AlphaFold2 prediction structure of the human TACAN protein (isoform CRA_a). .. figure:: Images/Fig17_workflow_1.svg :alt: Part of the general Scipion workflow performed to predict the AlphaFold2 structure of the human TACAN complex (isoform CRA_a). :name: model_building_fig17_workflow :align: center :width: 100.0% Part of the general Scipion workflow performed to predict the AlphaFold2 structure of the human TACAN complex (isoform CRA_a). Predicting the three-dimensional (3D) structure of a protein from its sequence alone remains an unsolved problem. However, some programs such as AlphaFold2 are able to predict the 3D atomic coordinates of a large collection of folded protein structures with remarcable accuracy :cite:p:`jumper2021`, :cite:p:`AlphaFold-Multimer2021`. Alphafold2 execution can take up to hours for a single protein, it needs more than 2 TB of disk space and the use of a GPU card is recommendable to speed up the execution of processes. The *Scipion* protocol **chimerax-alphafold prediction** has implemented the possibility of running Alphafold2 in your own computer selecting the option *Local AlphaFold* (:numref:`model_building_fig1_alphafold_prediction` (4), see :ref:`Appendix ` for details). .. figure:: Images/Fig1_protocol_prediction.svg :alt: Protocol **chimerax-alphafold prediction** detailing the 4 different options to available in *Scipion* to retrieve AlphaFold predicted models. :name: model_building_fig1_alphafold_prediction :align: center :width: 100.0% Protocol **chimerax-alphafold prediction** detailing the 4 different options to available in *Scipion* to retrieve AlphaFold predicted models. There are, nevertheless, different options that enable researchers retrieve AlphaFold2 predicted models if they lack of these computational resources and, with this aim, the *Scipion* protocol **chimerax-alphafold prediction** has also implemented three options to retrieve AlphaFold2 predicted models. The first two options take advantage of the AlphaFold protein structure `database `_ and give users the possibility of retrieving AlphaFold predicted models of the protein in which users are interested (option *EBI Database (identical sequence)*, :numref:`model_building_fig1_alphafold_prediction` (1)) or of a homologous protein, if that protein is not included in the database yet (option *EBI Database (homologous sequence)* :numref:`model_building_fig1_alphafold_prediction` (2)). Look at :ref:`Appendix ` for :ref:`use cases `. Take into account, however, that although the database is constantly growing with new additional predictions, the protein relevant for the user might not have close homologous. Executing Alphafold2 is thus required. The third option implemented in the *Scipion* protocol **chimerax-alphafold prediction** enables the user run the program taking advantage of two independent solutions based on Google Colaboratory (:numref:`model_building_fig1_alphafold_prediction` (3), option *Google Colab*) involving *Phenix* and *Chimera21* Colab Notebooks, respectively. In this tutorial we will use one of the solutions involving the option *Google Colab* to create an initial atomic model of the protein *TACAN* (isoform CRA_a). As mentioned before, *TACAN* isoform induced by tumor necrosis factor alpha probably folds as a homodimer of the protein with Uniprot ID `A0A024R4K9 `_. Let us start by importing the sequence *A0A024R4K9* (:numref:`model_building_fig2_alphafold_prediction`) using the protocol **pwem-import-sequence**. Complete the protocol form with your prefered name of the sequence (1), the type of sequence (2, aminoacids), the database source (3, UniProt), and the UniProt ID (4). Finally, execute the protocol (5). Look at protocol :ref:`appendix ` for details. .. figure:: Images/Fig2_protocol_prediction.svg :alt: Completing the *Scipion* protocol **pwem-import-sequence** to import the sequence of TACAN protein, isoform CRA_a, from UniProt database (workflow step 1; :numref:`model_building_fig17_workflow`). :name: model_building_fig2_alphafold_prediction :align: center :width: 80.0% Completing the *Scipion* protocol **pwem-import-sequence** to import the sequence of TACAN protein, isoform CRA_a, from UniProt database (workflow step 1; :numref:`model_building_fig17_workflow`). Protein sequences are the only input required to predict the atomic structure of a complex running AlphaFold2. Templates from PDB or a particular user's template can also be included as inputs depending on the Colab Notebook. In this tutorial we are going to execute two different examples using *Phenix* and *Chimera21* Colab Notebooks, respectively, in the *Scipion* protocol **chimerax-alphafold prediction**. AlphaFold initial model of monomer ---------------------------------- In the first one we will execute AlphaFold2 based on *Phenix* Colab Notebook and a single copy of the sequence as input. Fill in the protocol form as indicated in :numref:`model_building_fig3_alphafold_prediction`. Select *Google Colab* option (1), *Phenix* Colab Notebook (2), include the sequence of TACAN protein (isoform CRA_a) (3) and execute the protocol (4). For details go to the protocol :ref:`appendix `. .. figure:: Images/Fig3_protocol_prediction.svg :alt: **A**: Completing the *Scipion* protocol **chimerax-alphafold prediction** to predict the structure of TACAN protein (isoform CRA_a) with AlphaFold2 using *Phenix* Colab Notebook (workflow step 2; :numref:`model_building_fig17_workflow`). **B**: Warning message about the requirement of RAM. **C**: Warning message about the authorship of the Notebook. :name: model_building_fig3_alphafold_prediction :align: center :width: 100.0% **A**: Completing the *Scipion* protocol **chimerax-alphafold prediction** to predict the structure of TACAN protein (isoform CRA_a) with AlphaFold2 using *Phenix* Colab Notebook (workflow step 2; :numref:`model_building_fig17_workflow`). **B**: Warning message about the requirement of RAM. **C**: Warning message about the authorship of the Notebook. After executing the protocol a couple of warning messages usually appear (:numref:`model_building_fig3_alphafold_prediction` (B and C)). Accept them with *OK* and *Run anyway*, respectively, to run the AlphaFold2 prediction. The browser of *Phenix* Colab Notebook them will start to run, as you can observe in each of the independent running cells (:numref:`model_building_fig4_alphafold_prediction` (3, 4)), although control running and stopping of all of them is possible openning the *Runtime* menu (:numref:`model_building_fig4_alphafold_prediction` (1)). Take into account that you need a Google account (2) and the sequence of the protein should be written in the form (4). .. figure:: Images/Fig4_protocol_prediction.svg :alt: Browser of *Phenix* Colab Notebook completed with the TACAN protein (isoform CRA_a) sequence and running. :name: model_building_fig4_alphafold_prediction :align: center :width: 100.0% Browser of *Phenix* Colab Notebook completed with the TACAN protein (isoform CRA_a) sequence and running. When the execution finishes, since the default option *Show results in ChimeraX* has been set (:numref:`model_building_fig3_alphafold_prediction` (A)), the prediction structure of TACAN protein (isoform CRA_a) will be opened in the *ChimeraX* GUI. The user can close this *ChimeraX* session if (s)he doesn't save any modification in the structure. After clicking ``Analyze Results`` in the *Scipion* framework, the output result of AlphaFold2 execution can be visualized anyway (:numref:`model_building_fig5_alphafold_prediction`). .. figure:: Images/Fig5_protocol_prediction.svg :alt: Prediction of TACAN protein (isoform CRA_a) structure obtained using the *Phenix* Colab Notebook. **A**: Protein structure in *ChimeraX* GUI. **B**: PAE plot. :name: model_building_fig5_alphafold_prediction :align: center :width: 100.0% Prediction of TACAN protein (isoform CRA_a) structure obtained using the *Phenix* Colab Notebook. **A**: Protein structure in *ChimeraX* GUI. **B**: PAE plot. The unique prediction structure of TACAN protein, isoform CRA_a, (:numref:`model_building_fig5_alphafold_prediction` (A, model #2)) is colored according to the *Model Confidence* score. The reoriented model observed in Fig. :numref:`model_building_fig6_alphafold_prediction` (A) helps to observe in detail colors of the two helices diametrically opposed (1 and 2) and connected by linker 3. As you can check in *"How confident should I be in a prediction?"* in AlphaFold Protein Structure Database `FAQ `_, only the residues of alpha helices (1 and 2) display acceptable confidence pLDDT values (color blue and close to blue), although the confidence of helix 1 residues seems to be better than the observed for helix 2. The residues included in ribbons (yellow and orange) are disordered and can't be interpreted in general. .. figure:: Images/Fig6_protocol_prediction.svg :alt: Prediction of the isoform CRA_a of TACAN protein structure obtained using the *Phenix* Colab Notebook. **A**: Protein structure in *ChimeraX* GUI reoriented compared to :numref:`model_building_fig5_alphafold_prediction`. **B**: PAE plot split in homogeneous squares. :name: model_building_fig6_alphafold_prediction :align: center :width: 100.0% Prediction of the isoform CRA_a of TACAN protein structure obtained using the *Phenix* Colab Notebook. **A**: Protein structure in *ChimeraX* GUI reoriented compared to :numref:`model_building_fig5_alphafold_prediction`. **B**: PAE plot split in homogeneous squares. Regarding the relative position of domains inferred from PAE plot (Predicted Aligned Error plot reporting AlphaFold’s expected position error at residue x, when the predicted and true structures are aligned on residue y), several conclusions can be extracted. Look at *"How should I interpret the relative positions of domains?"* in AlphaFold Protein Structure Database `FAQ `_ for a detailed explanation. The big dark green square of the upper left corner (:numref:`model_building_fig5_alphafold_prediction` (B)) corresponds to the N-terminal part of the protein where two sub-domains are included (:numref:`model_building_fig6_alphafold_prediction` (B, 1 and 2)). These domains are the two alpha helices 1 and 2 of :numref:`model_building_fig6_alphafold_prediction` (A). The relative orientation of these two sub-domains is assessed by the squares with orange frame (:numref:`model_building_fig6_alphafold_prediction` (B, 3)). Although the error values of these squares are worse than the error values of sub-domains 1 and 2, they are still much better than the background, indicating that the relative position of the sub-domains might be close to the real one. The orientation of the yellow linker represented in :numref:`model_building_fig6_alphafold_prediction` (A, 3) could be reliable. A similar prediction could have been obtained using the *Chimera21* Colab Notebook. We have completed the protocol with this aim, as shown in :numref:`model_building_fig7_alphafold_prediction` (A). Again, the option *Google Colab* has been selected (1), as well as the *Chimera21* Colab Notebook (2). To avoid using more computational resources than are available from Google Colab, we maintain the default option to *Skip Amber relaxation* (3). The last step of molecular dynamics performed with the best prediction structure is thus skipped, which also save processing time. Take into account that surpassing the available resources might drive to the prediction failing. And finally, the sequence of the isoform CRA_a of TACAN protein has been included (4). .. figure:: Images/Fig7_protocol_prediction.svg :alt: **A**: Completing the *Scipion* protocol **chimerax-alphafold prediction** to predict the structure of TACAN protein (isoform CRA_a) with AlphaFold2 using *Chimera21* Colab Notebook (workflow step 3; :numref:`model_building_fig17_workflow`). **B**: Browser of *Chimera21* Colab Notebook completed with the TACAN protein sequence and running. :name: model_building_fig7_alphafold_prediction :align: center :width: 100.0% **A**: Completing the *Scipion* protocol **chimerax-alphafold prediction** to predict the structure of TACAN protein (isoform CRA_a) with AlphaFold2 using *Chimera21* Colab Notebook (workflow step 3; :numref:`model_building_fig17_workflow`). **B**: Browser of *Chimera21* Colab Notebook completed with the TACAN protein sequence and running. Compared with the browser of *Phenix* Colab Notebook, the browser of *Chimera21* Colab Notebook contains only one running cell (:numref:`model_building_fig7_alphafold_prediction` (B)). Check that the protein sequence has been included before starting the running. In this particular case, upstream the sequence appears the command *"dont_minimize"* referring the option *Skip Amber relaxation*. When the run time finishes, the **5** resulting unrelaxed models of the isoform CRA_a of TACAN will appear in the *ChimeraX* GUI (:numref:`model_building_fig8_alphafold_prediction` (A)). ``Remark`` that *Phenix* Colab Notebook generates a unique structure prediction. Observing one by one the predicted structures two conclusions can be inferred: a) Residues from the horizontal helix show on average better confidence score values than residues from vertical helices. b) Comparing the five predictions, the agreement in the prediction of the horizontal helix seems to be higher that in the prediction of the vertical helix. c) There are accuracy differences among predictions and one displays the highest confidence score values. This model (number #3 in *Models* panel) is the best model. Rotate, show and hide models to confirm the result. .. figure:: Images/Fig8_protocol_prediction.svg :alt: Prediction of the isoform CRA_a of TACAN protein structure obtained using the *Chimera21* Colab Notebook. **A**: 5 structure predictions visualized in *ChimeraX GUI*. **B**: Best structure prediction of TACAN protein displayed in the *ChimeraX viewer*. **C**: Plot *Number of Aligned Sequences with no Gap for each Residue Position*. **D**: PAE plot. :name: model_building_fig8_alphafold_prediction :align: center :width: 100.0% Prediction of the isoform CRA_a of TACAN protein structure obtained using the *Chimera21* Colab Notebook. **A**: 5 structure predictions visualized in *ChimeraX GUI*. **B**: Best structure prediction of TACAN protein displayed in the *ChimeraX viewer*. **C**: Plot *Number of Aligned Sequences with no Gap for each Residue Position*. **D**: PAE plot. The best model is individually displayed in the *ChimeraX* GUI viewer opened adter clicking ``Analyze Results`` (:numref:`model_building_fig8_alphafold_prediction` (B)). A couple of additional plots are opened simultaneously. The first one is the plot *Number of Aligned Sequences with no Gap for each Residue Position* that provides the map of the multiple alignment coverage, that is, given a multiple alignment containing N proteins how many of these N proteins contain each one of the problem protein aminoacids (:numref:`model_building_fig8_alphafold_prediction` (C)). This plot shows on the left the residues with high coverage of aligned sequences. This part of the sequence coincides with the N-terminal end of the protein for which AlphaFold2 achieved a quite fair prediction. The C-terminal part of the isoform CRA_a of TACAN doesn't show a minimum coverage of aligned sequences, thus avoiding the structure prediction of that sequence. PAE plot is the second plot opened to visualize the results (:numref:`model_building_fig8_alphafold_prediction` (D)). Quite similar to the same plot obtained from *Phenix* Colab Notebook (:numref:`model_building_fig6_alphafold_prediction` (B)), the PAE plot obtained with *Chimera21* Colab Notebook shows a slightly better prediction of the inter-domain linker. If you try to superpose the models generated both with *Phenix* and *Chimera21* Colab Notebooks with the *Scipion* protocol :ref:`Superpose PDBs ` you will see that there is not an obvious superposition (workflow step 4; :numref:`model_building_fig17_workflow`). This indicates that both prediction models are quite different in the spatial orientation of helices, which might be expected due to the low accuracy of residues that are not included in the two helices. AlphaFold initial model of multimer ----------------------------------- One interesting advantage of *Chimera21* Colab Notebook compared with *Phenix* Colab Notebook is that the first one enables structure prediction of multimers. As it was mentioned :ref:`before `, the isoform CRA_a of TACAN protein oligomerizes probably as a homodimer and we can predict directly the structure of this protein complex using AlphaFold2. Let us execute again AlphaFold with the protocol **chimerax-alphafold prediction** completing it with two identical sequences of TACAN that should be part of the homodimer. ``WARNING:`` Take into account that you have to import twice the same sequence in the *Scipion* framework (workflow step 5; :numref:`model_building_fig17_workflow`). Observe that this is the only difference between the prediction of monomer (:numref:`model_building_fig7_alphafold_prediction` (4)) and dimer (:numref:`model_building_fig9_alphafold_prediction` (1)). .. figure:: Images/Fig9_protocol_prediction.svg :alt: Completing the *Scipion* protocol **chimerax-alphafold prediction** to predict the structure of the isoform CRA_a of TACAN homodimer with AlphaFold2 using *Chimera21* Colab Notebook (workflow step 6; :numref:`model_building_fig17_workflow`). :name: model_building_fig9_alphafold_prediction :align: center :width: 100.0% Completing the *Scipion* protocol **chimerax-alphafold prediction** to predict the structure of the isoform CRA_a of TACAN homodimer with AlphaFold2 using *Chimera21* Colab Notebook (workflow step 6; :numref:`model_building_fig17_workflow`). When the run time finishes, the **5** resulting unrelaxed models of TACAN homodimer (isoform CRA_a) will appear in the *ChimeraX* GUI. This window can be closed if no modifications are saved. Then, clicking on ``Analyze Results`` three windows will pop-up again allowing the visualization of the best homodimer structure prediction (:numref:`model_building_fig10_alphafold_prediction`). .. figure:: Images/Fig10_protocol_prediction.svg :alt: Prediction of the isoform CRA_a of TACAN homodimer structure obtained using the *Chimera21* Colab Notebook. **A**: Best structure prediction of TACAN protein displayed in the *ChimeraX viewer*. **B**: PAE plot. **C**: Plot *Number of Aligned Sequences with no Gap for each Residue Position*. :name: model_building_fig10_alphafold_prediction :align: center :width: 100.0% Prediction of the isoform CRA_a of TACAN homodimer structure obtained using the *Chimera21* Colab Notebook. **A**: Best structure prediction of TACAN protein displayed in the *ChimeraX viewer*. **B**: PAE plot. **C**: Plot *Number of Aligned Sequences with no Gap for each Residue Position*. Concerning the confidence score values of each residue (:numref:`model_building_fig10_alphafold_prediction` (A)), we can see that most of the residues included in the N-terminal part (helix-linker-helix) of the each protein (1 and 2) seem to show better confidence scores in the homodimer than in any of monomer predictions. The PAE plot (:numref:`model_building_fig10_alphafold_prediction` (B)) shows a small error for this part of the protein in each monomer (boxes 1 and 2), as well as the relative position between the N-terminal parts of each monomer (boxes 3). Of course, the disordered region continues disordered in the dimer, mainly because no additional sequences have been found to align in the C-terminal end of the TACAN protein, as can be observed in :numref:`model_building_fig10_alphafold_prediction` (C). ----------- NOTE: Google Colab notebooks have an idle timeout of 90 minutes and absolute timeout of 12 hours. This means, if user does not interact with his Google Colab notebook for more than 90 minutes, its instance is automatically terminated. So, is there any way to programmatically prevent Google Colab from disconnecting on a timeout? In the past we have been able to implement the solution provided by seyon-seyon `here `_.