Connect with us

Bioinformatics Programming

How to calculate center of mass of a protein structure using Python script?

Published

on

How to calculate center of mass of a protein structure using Python script?

Here is a Python script that helps you calculate the center of mass of a protein using the Pymol [1].

Make sure that you have installed Pymol on your system otherwise, you will get an error stating “No Pymol module found”. If Pymol is not in your path, then run this script from the Pymol directory itself.


import __main__
__main__.pymol_argv = [ 'pymol', '-qc']            #quiet and no GUI
import sys, time, os
import pymol

pymol.finish_launching()

##### Reading user input #####
spath = os.path.abspath(sys.argv[1])
filename = spath.split('/')[-1].split('.')[0]

##### Loading structures #####
pymol.cmd.load(spath, filename)
pymol.cmd.disable("all")
pymol.cmd.enable(filename)
pymol.cmd.centerofmass(filename)

##### Getting output in a file #####
output = pymol.cmd.centerofmass(filename)
sys.stdout = open('output.txt','w')
print ((filename), (output))

##### Quiting #####
pymol.cmd.quit()

If you don’t want to write into a file, then comment out the “Getting output in a file” section. It will display the output on the terminal.

Now, save this file. We are saving it as com.py. Run this file as shown below:

$ python com.py <pdb_file>

OR

$ python3 com.py <pdb_file>

If your input PDB file is in different directory, then provide full path to it.


References

  1. The PyMOL Molecular Graphics System, Version 1.2r3pre, Schrödinger, LLC.

Dr. Muniba is a Bioinformatician based in New Delhi, India. She has completed her PhD in Bioinformatics from South China University of Technology, Guangzhou, China. She has cutting edge knowledge of bioinformatics tools, algorithms, and drug designing. When she is not reading she is found enjoying with the family. Know more about Muniba

Bioinformatics Programming

vs_analysis_compound.py: Python script to search for binding affinities based on compound names.

Published

on

vs_analysis_compound.py: Python script to search for binding affinities based on compound names.

Previously, we have provided the vs_analysis.py script to analyze virtual screening (VS) results obtained from Autodock Vina. In this article, we have provided another script to search for binding affinity associated with a compound. (more…)

Continue Reading

Bioinformatics Programming

How to download files from an FTP server using Python?

Published

on

How to download files from an FTP server using Python?

In this article, we provide a simple Python script to download files from an FTP server using Python. (more…)

Continue Reading

Bioinformatics Programming

How to convert the PDB file to PSF format?

Published

on

How to convert the PDB file to PSF format?

VMD allows converting PDB to PSF format but sometimes it gives multiple errors. Therefore, in this article, we are going to convert PDB into PSF format using a different method. (more…)

Continue Reading

LATEST ISSUE

ADVERT