Alphafold
Alphafold
AlphaFold is an artificial intelligence (AI) program developed by Alphabets's/Google's DeepMind which performs predictions of protein structure.
How to use
use run_alphafold.sh script located at /home/alphafold_folder/alphafold_multimer_non_docker
script reference:
Usage: run_alphafold.sh <OPTIONS> Required Parameters: -d <data_dir> Path to directory with supporting data: AlphaFold parameters and genetic and template databases. Set to the target of download_all_databases.sh. -o <output_dir> Path to a directory that will store the results. -f <fasta_path> Path to a FASTA file containing a single sequence. -t <max_template_date> Maximum template release date to consider (ISO-8601 format: YYYY-MM-DD). Important if folding historical test sets. Optional Parameters: -n <openmm_threads> OpenMM threads (default: all available cores) -b <benchmark> Run multiple JAX model evaluations to obtain a timing that excludes the compilation time, which should be more indicative of the time required for inferencing many proteins (default: false) -g <use_gpu> Enable NVIDIA runtime to run with GPUs (default: true) -a <gpu_devices> Comma separated list of devices to pass to 'CUDA_VISIBLE_DEVICES' (default: 0) -m <model_preset> Choose preset model configuration - the monomer model (monomer), the monomer model with extra ensembling (monomer_casp14), monomer model with pTM head (monomer_ptm), or multimer model (multimer) (default: monomer) -p <db_preset> Choose preset MSA database configuration - smaller genetic database config (reduced_dbs) or full genetic database config (full_dbs) (default: full_dbs) -u <use_precomputed_msas> Whether to read MSAs that have been written to disk. WARNING: This will not check if the sequence, database or configuration have changed. (default: false) -r <remove_msas_after_use> Whether, after structure prediction(s), to delete MSAs that have been written to disk to significantly free up storage space. (default: false) -i <is_prokaryote> Optional for multimer system, not used by the single chain system. This should contain a boolean specifying true where the target complex is from a prokaryote, and false where it is not, or where the origin is unknown. These values determine the pairing method for the MSA (default: false)
Databases
We downloaded the databases to /home/alphafold_folder/alphafold_data on compute-0-300 you may use it, or copy it to your own storage and point to it with -d flag of the run script. also, you may download the databases to your own storage via the script download_all_data.sh located at /home/alphafold_folder/alphafold_multimer_non_docker/scripts/
Sample qsub script
Note:
- Bulleted list item
create folder for output in your home dir mkdir ~/alphafold_output then run the script
- you may download dummy_test folder from this github as well for the output
https://github.com/kalininalab/alphafold_non_docker
- /home/alphafold_folder/alphafold_multimer_non_docker/example/query.fasta = this is sample data, please point to the data you need to query.
#!/bin/bash #PBS -l select=1:ncpus=4:ngpus=1 #PBS -q gpu # load miniconda module load miniconda/miniconda3-4.7.12-environmentally # activate relevant venv conda activate /powerapps/share/centos7/miniconda/miniconda3-4.7.12-environmentally/envs/alphafold_non_docker # run alphafold bash /home/alphafold_folder/alphafold_multimer_non_docker/run_alphafold.sh -d /home/alphafold_folder/alphafold_data -o ~/alphafold_output -f /home/alphafold_folder/alphafold_multimer_non_docker/example/query.fasta -t 2020-05-14