Skip to main content

Amber parameter file generation using prepin or mol2?

I am working on a system that closely resembles this :


The oxygen is protonated so that there are three bonds involving oxygen.

I tried to create the prmtop and inpcrd files for this system using antechamber and xleap.

Using Antechamber, if I try to create to parameter files by using prepin file, I encounter lot of problems.

Initially I used a pdb file of the molecule generated using either vmd or molden. But antechamber finds error in running bondtype.

So, I shifted to gabedit to generate the structure and creating its pdb file. This works fine. Atleast antechamber does not get bondtype error here.

But, then, if you try to create prepin file using command:

$AMBERHOME/exe/antechamber -i F.pdb -fi pdb -o F.prepin -fo prepi -c bcc -s 2
you will face errors related to the charge on the system like:

INFO: Number of electrons is odd: 127
      Please check the total charge (-nc flag) and spin multiplicity (-m flag)


As suggested you need to use -nc and -m flags to take care of charge and multiplicity of the system.

Once you take care of this error, your antechamber runs smoothly and generates corresponding prepin file.

We normally directly jump to parameter check.

If you try to do parameter check using command:

$AMBERHOME/exe/parmchk -i F.prepin -f prepi -o F.frcmod

You will most probably end up with segmentation fault error.

If we open prepin file, we might find some errors corresponding to the extra proton that is added on the oxygen.

35  O1    o     M   33  32  30     1.254   121.207  -171.680 -0.227000
36                  X    1   0   1     2.474       nan -9836.530  0.000000

Atom 36 is the one which corresponds to the extra H on oxygen. The bond type, bond distance, angle, dihedral data reported appears all wrong.

If you try to manually rectify this, you may succeed in generating the frcmod file.
But the catch here is if you visualize the molecule in xleap, you will find that the extra proton is not located on the oxygen as we expect, instead it is located on any neighboring carbon, nitrogen, oxygen(basically sp2 or sp moiety with probability to accept a proton).
After failing in many attempts to locate proton on desired oxygen, I took another approach. I tried to generate prmtop and inpcrd files using mol2 files instead of prepin.

Out of experience, working with prepin files is easier than mol2. But in this case mol2 worked well.

Using command:

antechamber -i F.pdb -fi pdb -o F.mol2 -fo mol2 -c bcc -s 2

along with -nc and -m flags for charge and multiplicity,

the mol2 file generated is more realistic. Parameters for all the atoms including the extra proton are well defined.

Parmchk also works fine and on viewing the molecule in xleap, the proton is located on desired oxygen.










Comments


  1. Really good information to show through this blog. I really appreciate you for all the valuable information that you are providing us through your blog. python Online course

    ReplyDelete

Post a Comment

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-19 Projections -till 22nd July

How would you call the COVID-19 situation in India for the last week. In contrast to what I hear around, I would say, the situation in India is improving. Lets see the reasons why I feel so: - Recovery rate has gone up to 63%. So, there is almost 2/3rd chance of survival for every person getting infected by COVID - The mortality rate has now reduced down to ~2.5%, way down from 3.2% in a 10 days duration.This is happening because of more and more testing. Today's report said, we did record 3.2 lakh testing in a day. This implies that the possibility of death because of COVID-19 is getting low. It may mean,  most of the people around can do fine with a little low panic level. Of course, precautions are necessary, if you don't want to fall in that very low 2.5% mortality set. Coming down to numbers. I am happy to see that India deviated significantly from my mid-July projection of 10.25 lakhs. Good to see the lowered down number. This shows the kind of ...