Modeller-10.1 [1, 2] is the latest version of the software. We can perform multiple tasks using this software including homology modeling, comparison of alignments and structures, and so on. In this article, we are going to compare two pairwise alignments using Modeller-10.1.
To run modeller, you need to enter a license key. If you have entered a license key during its installation, then you can proceed further otherwise you will have to define one as shown below. You can easily obtain a license key from here. It is free for academic users.
- Go to home/user/Downloads/modeller-10.1/modlib/modeller
- Open file ‘config.py’ in a text editor.
- Write the license key in the second line next to “license=“.
For comparing alignments, you need two input PDB files and a Python script, namely, ‘compare_alignments.py’. You can find this script inside your modeller directory. For example, if you have installed it in the Downloads directory, then the path to the script would be:
Remember, here user refers to the username in your system.
Preparing input files
Download PDB files that you want to compare and save them inside /home/user/Downloads/modeller-10.1/examples/atom_files directory. You can also remove extra chains from the file.
Preparing Python script
Open the Python script ‘compare_alignments.py‘ and write the names of the input PDB files that you want to compare as shown below.
# Generate and save sequence-sequence alignment: aln = Alignment(env) for code in ('1fas', '2ctx'): mdl = Model(env, file=code) aln.append_model(mdl=mdl, align_codes=code, atom_files=code) aln.align(gap_penalties_1d=(-600, -400)) aln.write(file='toxin-seq.ali')
Change the names of the files in line 3 shown in the above code snippet, i.e., “for code in (‘1fas’, ‘2ctx’):”.
Here, we are using the default files: 1fas & 2ctx.
After changing the names, save this file.
Use modpy.sh file to run modeller instead of mod10.1. mod10.1 will keep giving you warnings and errors such as
‘import site’ failed
‘Could not find platform independent libraries <prefix>’.
To run modeller using modpy.sh, you need Python 2.3 and 3.9 installed.
Now run the script using the following command.
$ /home/user/Downloads/modeller-10.1/bin/modpy.sh python3 /home/user/Downloads/modeller-10.1/examples/commands/compare_alignments.py
If you are not using Python3 then just write python. It will take a few seconds to finish the job and the output will look like this.
- Webb, B., & Sali, A. (2014). Comparative protein structure modeling using MODELLER. Current protocols in bioinformatics, 47(1), 5-6.
- Martí-Renom, M. A., Stuart, A. C., Fiser, A., Sánchez, R., Melo, F., & Šali, A. (2000). Comparative protein structure modeling of genes and genomes. Annual review of biophysics and biomolecular structure, 29(1), 291-325.
Video Tutorial: Calculating binding pocket volume using PyVol plugin.
How to generate topology from SMILES for MD Simulation?
[Tutorial] Installing jdock on Ubuntu (Linux).
jdock is an extended version of idock . It has the same features as the idock along with some bug fixes. However, the binary name and the GitHub repository names are changed. We are installing jdock on Ubuntu (Linux).
How to upgrade cmake on Ubuntu (Linux)?
In bioinformatics, cmake is used to install multiple software including GROMACS, jdock, and so on. Here is a short tutorial on how to upgrade cmake on Ubuntu and get rid of the previous version. (more…)
How to install GMXPBSA on Ubuntu (Linux)?
[Tutorial] Installing Pyrx on Windows.
How to solve ‘Could NOT find CUDA: Found unsuitable version “10.1”‘ error during GROMACS installation?
Compiling GROMACS  with GPU can be trivial. Previously, we have provided a few articles on the same. In this article, we will solve an error frequently occurring during GROMACS installation.
Installing Autodock4 on MacOS.
How to install Autodock4 on Ubuntu?
DS Visualizer: Uses & Applications
Protein structure & folding information exploited from remote homologs.
Remote homologs are similar protein structures that share similar functions, but there is no easily detectable sequence similarity in them. A new study has revealed that the protein folding information can be exploited from remote homologous structures. A new tool is developed to recognize such proteins and predict their structure and folding pathway. (more…)
Pathonoia- A new tool to detect pathogens in RNA-seq data.
Detecting viruses and bacteria in RNA-seq data with less false positive rate is a difficult task. A new tool is introduced to detect pathogens in RNA-seq data with high precision and recall known as Pathonoia .
AlphaFill- New algorithm to fill ligands in AlphaFold models.
AlphaFold is a popular artificial intelligence based protein prediction tool . Though it predicts good protein structures, it lacks the capability to predict the small molecules present in the structure such as ligands. For this purpose, AlphaFill is introduced by Hekkelman et al.,. (more…)
How to calculate binding pocket volume using PyVol plugin in PyMol?
How to generate electron density map using Pymol?
Installing PyVOL plugin in Pymol on Ubuntu (Linux).
How to execute matlab from terminal in Ubuntu (Linux)?
How to install Kpax on Ubuntu (Linux)?
Kpax is a bioinformatics program to search and align protein structures . It is currently available for Linux platforms only. In this article, we are going to install the latest version of Kpax (5.1.3) on Ubuntu (Linux). (more…)
How to run do_dssp command (mkdssp) in Gromacs 2022?
In the latest version of GROMACS (2022) , there are some issues regarding the gmx do_dssp command. Apparently, this command either does not run displaying a fatal error, or if it runs then it does not read any frame from MD simulation files. In this article, we are going to run the same command for GROMACS 2022. (more…)