Open Zkkkkkui opened 1 year ago
Yes. The input feature generation is similar to AF2. A sample shell script is here for individual sequences.
What about this .py file: run_af2c_fea.py, which is also said to be used to get features?
Yes. The input feature generation is similar to AF2. A sample shell script is here for individual sequences.
What about this .py file: run_af2c_fea.py, which is also said to be used to get features?
The shell script calls this python script you referred to do the job.
Yes. The input feature generation is similar to AF2. A sample shell script is here for individual sequences.
What about this .py file: run_af2c_fea.py, which is also said to be used to get features?
The shell script calls this python script you referred to do the job.
Could you please give an example of the script on colab just to generate the features of a protein sequence from uniprot?
CalledProcessError Traceback (most recent call last)
The Colab notebook we provided only takes features.pkl
files of individual monomers. If you use other AlphaFold notebook to generate input features, make sure that you use the monomer, not multimer, pipeline to generate the input. And tar these pickle files into one single tar ball.
In this example, you have a heterodimer HgcAB composed of two monomers, HgcA and HgcB. Organizes the feature input as the following:
hgc.tar
├── hgc
│ ├── HgcA
│ │ └── features.pkl.gz
│ └── HgcB
│ └── features.pkl.gz
Then tar this folder into a single tarball and upload it to our notebook. Note that our code can take gzipped pickle files directly. It is up to you whether or not to gzip the pickle files before you make the tarball.
After you upload the attached tarball, you may run a test to predict a heterodimer using the target syntax: HgcA/HgcB 433 HgcAB
Thank you for the instruction! I have successfully got the features from other AF notebook(it used 0 sequence template and I am not sure whether it would matter compared to your examples) and did prediction on some protein complexes. However it seems to have a very high false negative rate(the proteins were supposed to be interacting but the output was not), is there any way to improve that?
Thank you for the instruction! I have successfully got the features from other AF notebook(it used 0 sequence template and I am not sure whether it would matter compared to your examples) and did prediction on some protein complexes. However it seems to have a very high false negative rate(the proteins were supposed to be interacting but the output was not), is there any way to improve that?
Many things to try, such as:
@FreshAirTonight In AlphaFold, MSAs are built using jackhmmer and HHblits. To avoid the extensive data downloads and CPU processing, precomputed MSAs and the feature.pkl file can be used instead.
Due to my limited resources, I'm focusing on the GPU-based second step and considering tools like ColabFold to create the feature.pkl files.
What files and steps are needed to create the feature.pkl file? Is there a tool available on Google Colab or Kaggle for this?
@fereidoon27 An example of feature generation script run_fea_gen.sh
can be found under the example
folder. If you have limited resources, consider using the uniprot
mode, under which MSA construction uses only the UniProt library (creating dummy files for other sequence library to get around file checking). With this option, you can generate features for hundreds of protein sequences of moderate lengths with a decent workstation (e.g., with 4TB nvme, 24 cores).
You may use precomputed MSAs as well. Under the intended output folder of a protein, create a subfolder named msas
and place the MSA files under that folder. The only MSAs required in the uniprot
mode is uniprot_hits.sto
. Add pdb_hits.sto
for templates if you need.
Yes. The input feature generation is similar to AF2. A sample shell script is here for individual sequences.