Question: Automation for Energy Minimizing protien structure
gravatar for yashwanthakumar
17 months ago by
yashwanthakumar10 wrote:

I have nearly 100 pdb structures to be energy minimized. I need a tool or script to run multiple energy minimization at a stretch. My actual work is to predict the changes in function of a protein with few mutations. So I m in search of automating this process. This is ligand-protein structure.

ADD COMMENTlink modified 17 months ago by fracarb850 • written 17 months ago by yashwanthakumar10
gravatar for fracarb8
17 months ago by
fracarb850 wrote:

Have a look at gromacs ( You can either use bash to wrap the gromacs calls (as I did) and batch process the pdbs or you can have a look at one of the wrappers on github. There is a bit of set up to do, but their documentation is quite good.

You may want to consider running some docking or MD simulations if you want to full study the ligand-protein bind over time.

During my Phd I study the effect of single point mutation of a protein structure, and gromacs was my tool of choice.

ADD COMMENTlink modified 17 months ago • written 17 months ago by fracarb850

I have been using gromacs for energy minimization following some tutorial (Tutorial link: ) But i have to manipulate the topology files manually for that reason only I m searching for some GUI based or any scripts.Please go through the step in the link above and give me a suggestion to automate that tutorial.

ADD REPLYlink written 17 months ago by yashwanthakumar10

It depends a bit where your mutations are. The topology of your 100 proteins are likely to be exactly the same, but for few residues. You can generate the topology for the wild-type and then replace the coordinates of the mutate residues, to obtain all the mutants you want.

The quickest way could be to use the "mutagenesis" function in pymol to mutate your residues programmatically (see A few lines of python should do the trick.

You can then replace the new coordinates using sed

ADD REPLYlink modified 17 months ago • written 17 months ago by fracarb850

Thanks for your reply. i have all the mutated pdb files now. I have to perform energy minimization for all the mutated structures. I cant run gromacs for all mutated structures one by one so I have to find a way for automating that.

ADD REPLYlink written 16 months ago by yashwanthakumar10

bash is your friend here :D I would suggest to have a look at github as there are several gromacs wrappers (e.g. GromacsWrapper).

If you need inspiration to write your own piece of code, you can have a look at the project I used.

Note that it is massively outdated: doitGROMACS_equilib

ADD REPLYlink modified 16 months ago • written 16 months ago by fracarb850

does the code will minimize the structure of ur git hub code

ADD REPLYlink written 4 months ago by sagarhm0320000
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1888 users visited in the last hour