Re: [AMBER] MCPB.py

From: Pengfei Li <ambermailpengfei.gmail.com>
Date: Tue, 27 Mar 2018 16:13:27 -0500

Hi Abhishek,

I may be because of the atom ID you give. In your NEP_RP_H_fixed_nonprot.pdb file (you did not offer me the NEP_RP_H_fixed.pdb file even it was cited in the NEP_RP.in file), the atom ID of Zinc is 11024 while you have given 11025 in the NEP_RP.in file. This inconsistency will cause an error. Please check about this issue for your modeling.

Kind regards,
Pengfei

> On Mar 26, 2018, at 11:39 AM, Thakur, Abhishek <axt651.miami.edu> wrote:
>
> 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
> <NEP_RP.in><ZN.mol2><MOL.mol2><NEP_RP_H_fixed_nonprot.pdb><HOH.frcmod><HOH.mol2><MOL.frcmod>_______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Mar 27 2018 - 14:30:02 PDT
Custom Search