Software
How to compare two pairwise alignments using Modeller-10.1?
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=“.
Files required
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:
/home/user/Downloads/modeller-10.1/examples/commands/compare_alignments.py
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.
Comparing alignments
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
or
‘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.
References
- 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.
Autodock suite is used for docking small molecules [1]. Recently, Autodock-GPU [2] is developed to accelerate the docking process. Its installation is described in this article. In this tutorial, we will install Autodock 4.2.6 on Ubuntu.
Discovery Studio (DS) Visualizer (from BIOVIA) is a visualization tool for viewing, sharing, and analyzing proteins [1]. Here are some uses and applications of DS Visualizer.
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…)
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 [1].
AlphaFold is a popular artificial intelligence based protein prediction tool [1]. 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.,[2]. (more…)
Electron density maps are available for most of the protein structures in PDB. Therefore, in this article, we are using PDB to generate electron density maps in Pymol.
PyVOL [1] is an excellent plugin of Pymol [2] for pocket visualization of proteins. In this article, we will install the PyVOL plugin in Pymol on Ubuntu. (more…)
While trying to install Matlab [1], it generally gives an error stating “matlab: command not found.”. Here, we provide a solution to this error.
Kpax is a bioinformatics program to search and align protein structures [1]. 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…)
In the latest version of GROMACS (2022) [1], 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…)
SMINA is a fork of AutoDock Vina supporting better scoring function and high-performance energy minimization [1]. In this article, we are going to install SMINA on Ubuntu (Linux). (more…)
ClusCo is a bioinformatics software to perform clustering and comparison of protein models [1]. In this article, we are going to install ClusCo on Ubuntu (Linux). (more…)
Running a plugin on an old version of Pymol [1] can give you multiple errors that are not easy to troubleshoot. For example, LigAlign plugin [2] runs on an old version of Pymol. Previously, we explained how to install LigAlign on Pymol. In this article, we will run the LigAlign command on Pymol. (more…)
Few errors appear when we try to run the LigAlign plugin [1] in Pymol [2]. For example, if you try to run the ligand_alignment plugin, it will give you multiple errors including “Unable to initialize LigAlign v1.00“, or “can’t run LigAlign v1.00” or “incorrect Python syntax” or “Plugin has been installed but initialization failed“. In this article, we explain the reason for this issue and how you can rectify these errors. (more…)
Sometimes we need to access old versions of Pymol [1]. Running some plugins on Pymol is difficult due to Python incompatibility. New Pymol versions require Python 3.x whereas older versions run on Python 2.x. Therefore, we need to maintain multiple versions of Pymol on a single system. In this article, we will install Pymol 1.7.x along with the latest version (Pymol 2.5.2) on Ubuntu. Later, we will create shortcuts for them.
Previously, we provided a tutorial for Pymol installation on Ubuntu. In this article, we are going to install Pymol on Mac OS. (more…)
In this article, we are going to install VMD [1] on Ubuntu. (more…)
We have provided a few articles on GROMACS installation on Ubuntu. In this article, we are going to install GROMACS [1] on Mac OS. (more…)
PyMol is one of the most widely used bioinformatics software. Generally, it is used as a molecular viewer to visualize macromolecules and small molecules. In this article, we are discussing several uses and applications of PyMol. (more…)
You must be logged in to post a comment Login