Skip to main content

Using Packmol

Packmol is a software that allows you to generate a single combined PDB file for all the molecules you are considering.
It is a free software and can be downloaded from this link:

http://www.ime.unicamp.br/~martinez/packmol/

The installation guide is also provided on the page and its really easy to follow.http://www.ime.unicamp.br/~martinez/packmol/

You can also find some examples but I tried on my own.
Anyway what you need is an input file which looks like:

------------------
tolerance 2.0
output trial.pdb
add_amber_ter
filetype pdb
structure 2a.pdb
  number 1
  inside cube 0. 0. 0. 10.
end structure
structure 2d.pdb
  number 1
  inside cube 0. 0. 0. 10.
end structure
---------------------------------

IT includes some information like:
how close the molecules considered can be.(tolerance)

name of output pdb file
"add_amber_ter" adds term "TER" after every molecule thus defining the start and end of every molecule.
After mentioning the filetype it has information of the input pdb files.
You can add as many different molecules here.
If you want repeating unit of same molecule then just mention information of one molecule and change the number setting from 1 to n, n is the number of units of that molecule.

I am showing here and example with pdb files of 2 molecules to be merged into one single PDB file.Name of the input file is trial.inp and output will be trial.pdb

you can execute packmol using command:

./packmol < trial.inp

You might require to launch packmol from its home directory packmol.
Sometimes I use command like:

../../../packmol/./packmol < trial.inp

The output is trial.pdb file that contains information of both the molecules.
Here is the output:

-----------------
HEADER
TITLE     Built with Packmol
REMARK   Packmol generated pdb file
REMARK   Home-Page: http://www.ime.unicamp.br/~martinez/packmol
REMARK
ATOM      1  C   NIT A   1       1.042   7.542   8.196  1.00  0.00
ATOM      2  N   NIT A   1       2.161   7.815   8.002  1.00  0.00
ATOM      3  C1  NIT A   1       3.530   8.137   7.761  1.00  0.00
ATOM      4  H   NIT A   1       4.155   7.317   8.121  1.00  0.00
ATOM      5  H1  NIT A   1       3.785   9.052   8.302  1.00  0.00
ATOM      6  H2  NIT A   1       3.683   8.275   6.691  1.00  0.00
TER
ATOM      7  OH  HDX B   1       0.527   9.826   9.789  1.00  0.00
ATOM      8 1HH  HDX B   1       0.526   8.877   9.941  1.00  0.00
ATOM      9 2HH  HDX B   1       0.536   9.998   8.845  1.00  0.00
TER
END

--------------------

You will get such a pdb file if you have generated your input pdb files by the method mentioned in previous posts. Check the unit names "NIT" and "HDX".
If your final pdb file appears like this, then it will compatible to the lib file you generated earlier. Otherwise you will landup with errors while try to load it on xleap later on.

Comments

Popular posts from this blog

Material Studio - selective molecule movement

So, working with Material Studio is awesome. You get to handle bigger, complicated and even unimaginable system. But I am facing problems in very minute functioning. I recently faced two problems, but managed to resolve on my own. Somehow, tutorials doesn't help in resolving them and you end up doing some exercise on your own to get the solution. Among these, the problem of bringing two molecular structures in one window and then adjusting one molecule by moving it around, while other structure remains untouched, as it is. Actually, this should not be a problem to do. There is a simple way of doing this. Select the molecule to be moved and then there is a symbol (button) to select, which allows you to selectively move that molecule only. This is the same way it can be done in Material studio too.  But among all the dozens of symbols, the new person handling MS, kind of gets lost. So, what he does is check the tutorials, for any help. But the tutorials don't talk about suc...

How would COVID bahave in August - few projections

Entering August with COVID by side   Entering in the month of August, COVID is still here by our side. Situation seems grim, worse than earlier, as number of active cases as well as the count of deceased patients are increasing day by day. Yesterday we crossed the mark of 50,000 cases in a day. Closely competing with Brazil with US counts also in range now. On the other side, few other numbers like reduced mortality rate to 2.18%, recovery rate of 64.6% and slightly increased doubling rate of 21 days can reduce the anxiety. My projections of 16.75 lakhs for India might be crossed by coupe of thousands of more cases after today’s count is added. This number would reach to a huge count of 28.5 lakhs by August 15. Maharashtra just doesn’t like competition. Just like what Delhi tried few days back, Andhra Pradesh also managed to cross Maharashtra’s per day count for2 consecutive days. Seems like not a good idea. Maharashtra immediately went to the next zone. The state which w...

COVID projections till July end

This last week seems terrible. The kind of spurt we saw in last many days was unexpected. Yesterday's count was 45000+. I did not expect. Are coming days going to report bigger numbers than this? We need to understand the cycle it follows. Mumbai is a good example. Since beginning Mumbai was heading the race. While other cities were reporting the count in 10's or 100's, Mumbai was reporting in thousands. back then, as much as 60-70% of Mahrashtra's total count was coming from Mumbai. Now after a period of almost 3 months, Mumbai is lagging behind. Pune, Kalyan, Thane have overtaken Mumbai. The famous hotspots of Dharavi, Dadar, Sion, Worli are longer uttered in the news. Of course, the counts are still there, but they do not attract any more attention. This is a cycle which will be followed in every new red zone area. The point is how big spike we allow the zone to get and how quickly we nab it. Delhi at one point started competing with Mumbai/Maharashtra. But, it h...