Re: [AMBER] MCPB.py

From: Thakur, Abhishek <axt651.miami.edu>
Date: Mon, 26 Mar 2018 16:39:25 +0000

Hi Pengfei,


Thank you for your response. I have checked all those points in the system but still, it is giving the same error.

Will you mind looking into my system once?


I am sorry for disturbing you, but I have tried everything and still, it seems to give me a problem.

I will really appreciate your help.



Thanking you,

-Abhishek

________________________________
From: Pengfei Li <ambermailpengfei.gmail.com>
Sent: Monday, March 26, 2018 5:24:03 AM
To: AMBER Mailing List
Subject: Re: [AMBER] MCPB.py

Hi Abhishek,

It is an atom name error. Please check the following section in the tutorial: https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fambermd.org%2Ftutorials%2Fadvanced%2Ftutorial20%2Fmcpbpy.htm&data=02%7C01%7Caxt651%40miami.edu%7C362d4f5a07c84e4f91e708d593361541%7C2a144b72f23942d48c0e6f0f17c48e33%7C0%7C0%7C636576782858525722&sdata=1K3%2FGd03Dxwv%2FRQ2SEScrsDS2ebpJE5544G6G%2BwCols%3D&reserved=0 <https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fambermd.org%2Ftutorials%2Fadvanced%2Ftutorial20%2Fmcpbpy.htm&data=02%7C01%7Caxt651%40miami.edu%7C362d4f5a07c84e4f91e708d593361541%7C2a144b72f23942d48c0e6f0f17c48e33%7C0%7C0%7C636576782858525722&sdata=1K3%2FGd03Dxwv%2FRQ2SEScrsDS2ebpJE5544G6G%2BwCols%3D&reserved=0>.

“MCPB.py now supports more than 80 metal ions (see Figure 3 in the reference paper <https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fpubs.acs.org%2Fdoi%2Fabs%2F10.1021%2Facs.jcim.5b00674&data=02%7C01%7Caxt651%40miami.edu%7C362d4f5a07c84e4f91e708d593361541%7C2a144b72f23942d48c0e6f0f17c48e33%7C0%7C0%7C636576782858525722&sdata=EDCCarKbkMKR7BibB4R8ksoY%2FHyGkgpriP1tJbvcDPk%3D&reserved=0>). For all of these metal ions and halide ions (if in the ion format but not neutral format), capitalized element name are suggested to use for both the residue name and atom name. Please make sure the residue name and atom name of the metal ion (not only for zinc used here) are all captilized (e.g. herein all equals "ZN", not "Zn", "zN", or "zn"), only in this way MCPB.py will recognize it as a metal ion. Meanwhile, please make sure each metal ion or halide ion (if it is in its ion format but not neutral format) is treated seperately as an indepedent residue in the PDB file. If you have metal site which has metal ion embeded in the ligand residue, for example, the HEME group in PDB file, please seperate the metal ion into an independent residue (also, with an unique residue number and atom number)in the PDB file. Meanwhile, in order to make MCPB.py recognize the atoms well, if you have atoms have atom names capitalized with numbers in your original PDB file (applies for protein or ligand), please change these atom names to atom names capitalized with their element symbols (for example, change "2HA1" to "HA12", change "1HAA" to "HAA1") in the PDB file before performing following steps.”

Hope it helps,
Pengfei

> On Mar 22, 2018, at 2:12 PM, Thakur, Abhishek <axt651.miami.edu> wrote:
>
> Hi, I am having a small problem in generating gaussian input file using MCPB.py.
>
> Kindly help me in solving this problem.
>
> The error that I am getting is.
>
>
> ***Creating the small model...
> It contains the residue 530-HID as sidechain coordinated.
> It contains the residue 534-HID as sidechain coordinated.
> It contains the residue 593-GLU as sidechain coordinated.
> It contains the residue 696-ZN as normal.
> It contains the residue 697-MOL as normal.
> It contains the residue 698-HOH as normal.
> Totally there are 83 atoms in the small model.
> Traceback (most recent call last):
> File "/nethome/axt651/apps/amber16_tools/bin/MCPB.py", line 644, in <module>
> addred, lgchg, lgspin)
> File "/nethome/axt651/apps/amber16_tools/lib/python2.7/site-packages/pymsmt/mcpb/gene_model_files.py", line 1895, in gene_model_files
> smspin, addred, outf, sqmopt)
> File "/nethome/axt651/apps/amber16_tools/lib/python2.7/site-packages/pymsmt/mcpb/gene_model_files.py", line 1380, in build_small_model
> AtNum = Atnum[gatm.element]
> KeyError: '1'
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&data=02%7C01%7Caxt651%40miami.edu%7C362d4f5a07c84e4f91e708d593361541%7C2a144b72f23942d48c0e6f0f17c48e33%7C0%7C0%7C636576782858525722&sdata=eEbhtS24%2BGpmYRCfR0%2Br8EoAccjcupsoRoyK47Erb9A%3D&reserved=0

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&data=02%7C01%7Caxt651%40miami.edu%7C362d4f5a07c84e4f91e708d593361541%7C2a144b72f23942d48c0e6f0f17c48e33%7C0%7C0%7C636576782858525722&sdata=eEbhtS24%2BGpmYRCfR0%2Br8EoAccjcupsoRoyK47Erb9A%3D&reserved=0


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber

Received on Mon Mar 26 2018 - 10:00:06 PDT
Custom Search