Previous Thread
Next Thread
Print Thread
Page 1 of 2 1 2
Joined: Jan 2013
Posts: 79
T
Tarsis Offline OP
Forum Member
OP Offline
Forum Member
T
Joined: Jan 2013
Posts: 79
Hi all,

After setting up the patch for adding the sulfamate to the N- glucosamine, I got the folloing error:

CHARMM> patch nhshb S 2 setup
ATOM S AGLCNA 2 S ADDED.
ATOM S AGLCNA 2 O2 ADDED.
ATOM S AGLCNA 2 O7 ADDED.
ATOM S AGLCNA 2 O8 ADDED.

Message from MAPIC: Atom numbers are changed.
*** LEVEL -3 WARNING *** BOMLEV IS -3
BOMLEV HAS BEEN SATISFIED.supposedly TERMINATING.

I cannot explain which atom numbers changed and why. The patch supposedly adds SO3 to the AGLCNA already deacetylated by another patch.

Thank in advance

Tarsis

Attached Images
patch.tar.gz (56.45 KB, 772 downloads)
nhshb.out.txt (7.67 KB, 801 downloads)
NHSSHB.txt (1.47 KB, 843 downloads)
Joined: Dec 2005
Posts: 1,535
Forum Member
Offline
Forum Member
Joined: Dec 2005
Posts: 1,535
WRNLEV 1 BOMLEV -3 = . Sorry, you're on your own. Please review the READ THIS FIRST topic.

Edit: in case you wondered why we're so fanatic about our users not touching WRNLEV or BOMLEV, your fatal error is not caused by the change in atom numbers. CHARMM would likely tell you the real cause of the problem if you wouldn't suppress it by changing WRNLEV .

Last edited by Kenno; 08/21/13 05:34 PM.
Joined: Jan 2013
Posts: 79
T
Tarsis Offline OP
Forum Member
OP Offline
Forum Member
T
Joined: Jan 2013
Posts: 79
I'm sorry,

I ran it before with BOMLEV 0 but wanted to test whether it would work with -3.

That's how the script ends with BOMLEV 0

CHARMM> patch amin S 2 setup
ATOM S 2 AGLCNA C AND ALL REFERENCES TO IT DELETED.
ATOM S 2 AGLCNA O AND ALL REFERENCES TO IT DELETED.
ATOM S 2 AGLCNA CT AND ALL REFERENCES TO IT DELETED.
ATOM S 2 AGLCNA HT1 AND ALL REFERENCES TO IT DELETED.
ATOM S 2 AGLCNA HT2 AND ALL REFERENCES TO IT DELETED.
ATOM S 2 AGLCNA HT3 AND ALL REFERENCES TO IT DELETED.
ATOM S AGLCNA 2 S ADDED.
ATOM S AGLCNA 2 O2 ADDED.
ATOM S AGLCNA 2 O7 ADDED.
ATOM S AGLCNA 2 O8 ADDED.

Message from MAPIC: Atom numbers are changed.
DELTIC: 6 bonds deleted
DELTIC: 11 angles deleted
DELTIC: 13 dihedrals deleted
DELTIC: 2 improper dihedrals deleted
*** LEVEL -3 WARNING *** BOMLEV IS 0
BOMLEV HAS BEEN SATISFIED. TERMINATING.

ABNORMAL TERMINATION
MOST SEVERE WARNING WAS AT LEVEL -3

Attached Images
create_disa_nhshb.inp.txt (594 Bytes, 769 downloads)
nhshb.out.txt (6.82 KB, 790 downloads)
Joined: Dec 2005
Posts: 1,535
Forum Member
Offline
Forum Member
Joined: Dec 2005
Posts: 1,535
As mentioned in my last post, you need to remove the WRNLEV command as well. There's still information missing from your output.

Joined: Jan 2013
Posts: 79
T
Tarsis Offline OP
Forum Member
OP Offline
Forum Member
T
Joined: Jan 2013
Posts: 79
Yep, as you said it before, haste makes waste frown

I removed WRNLEV and got what the error was. I will try to fix it by myself and then post the results/questions when they appear.

Thanks again for the patience Kenno.

Joined: Jan 2013
Posts: 79
T
Tarsis Offline OP
Forum Member
OP Offline
Forum Member
T
Joined: Jan 2013
Posts: 79
I think I got it shocked (except by some nasty clashes between the sulfur S1 and O7 and HN nitrogen).

Does anyone know how to fix that?

Thanks all again,

Tarsis

Attached Images
patch.tar.gz (94.67 KB, 760 downloads)
disa.pdb.txt (3.74 KB, 836 downloads)
nhshb.out.txt (8.83 KB, 801 downloads)
Joined: Dec 2005
Posts: 1,535
Forum Member
Offline
Forum Member
Joined: Dec 2005
Posts: 1,535
Now you're talking! The clash can be relieved by editing the values (especially the dihedral) in the IC C1 C2 N HN entry. I speculate it arose after editing IC C1 C2 N S1 without doing the above. To make the hydrogen "follow" automatically when editing the dihedral associated with the heavy atoms, it is recommended to use an improper IC entry, like this:
Code:
IC  C2   S1   *N   HN         0.0000    0.00  120.00    0.00   0.0000
Where the 120 can be changed to -120 depending on the pyramid flip, and the zeroes can be substituted manually with the correct coordinates, or automatically with IC PARAM.

On an unrelated note:
  • You can suppress the charge warning at the beginning of your output by changing the total charge in your topology file like this: PRES AMIN -1.00 .
  • It is good practice to leave the main force field files untouched, and make a private self-contained stream file for any additions.

See the CGenFF files and tutorials for examples illustrating all of the above.

Joined: Jan 2013
Posts: 79
T
Tarsis Offline OP
Forum Member
OP Offline
Forum Member
T
Joined: Jan 2013
Posts: 79
Hi Kenno,

I managed to fix the clash by editing the IC entry as you said, tks by the way for all your help. But now two questions arose.

1) There seems to be a problem with the coordinates for the BGLC atoms (this error does not reflect in anything, just annoys me).

CHARMM> read coor univ unit 1 offset 0
SPATIAL COORDINATES BEING READ FROM UNIT 1
TITLE> NONE
TITLE> *
** WARNING ** NO COORDINATES FOR 1 1 BGLCA C1
** WARNING ** NO COORDINATES FOR 2 1 BGLCA H1
** WARNING ** NO COORDINATES FOR 3 1 BGLCA O1
** WARNING ** NO COORDINATES FOR 4 1 BGLCA HO1
** WARNING ** NO COORDINATES FOR 5 1 BGLCA C5

2) Where can I find impropers that may match those entries?

: No improper parameters for 2 ( SC OC2DP NG2S3 OC2DP )
: No improper parameters for 3 ( NG2S3 SC CC3161 HCP1 )
: No improper parameters for 4 ( SC OC2DP NG2S3 OC2DP )
: A TOTAL OF 3 MISSING PARAMETERS


I will, when done with this parametrization, generate the stream file. Thanks for reminding me that! smile

All the best

Joined: Sep 2003
Posts: 8,602
Likes: 21
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,602
Likes: 21
The "NO COORDINATES" error occurs because you may not using READ COOR UNIV correctly, and because the PDB format does not support the 5- and 6-char residue names used by the all36_carb force field.


Rick Venable
computational chemist

Joined: Dec 2005
Posts: 1,535
Forum Member
Offline
Forum Member
Joined: Dec 2005
Posts: 1,535
Concerning the impropers, you should never define them on nonplanar/tetrahedral centers (and not even on all planar centers; see this FAQ). So just remove the 3 superfluous IMPR lines from your topology file and you're done (I'm assuming the sulfamate nitrogen is not planar).

Joined: Jan 2013
Posts: 79
T
Tarsis Offline OP
Forum Member
OP Offline
Forum Member
T
Joined: Jan 2013
Posts: 79
Got it, thanks Rick.

Joined: Jan 2013
Posts: 79
T
Tarsis Offline OP
Forum Member
OP Offline
Forum Member
T
Joined: Jan 2013
Posts: 79
Yes, I removed those improper lines but after running a simple minimization and that group clash again. Isn't that cause by the missing impropers?

Attached Images
disa-initialmini.crd.txt (6.58 KB, 763 downloads)
disa-initialmini.pdb.txt (3.74 KB, 764 downloads)
disa-initialmini.psf.txt (15.43 KB, 2677 downloads)
in_glc.out.txt (314.42 KB, 778 downloads)
Initial-minimization_glc.inp.txt (817 Bytes, 770 downloads)
Joined: Dec 2005
Posts: 1,535
Forum Member
Offline
Forum Member
Joined: Dec 2005
Posts: 1,535
Nope, it's a (frequently made) mistake in the script that generated the psf file. You need to issue the command "AUTOGENereate ANGLes DIHEdrals" after you're finished patching.

Joined: Jan 2013
Posts: 79
T
Tarsis Offline OP
Forum Member
OP Offline
Forum Member
T
Joined: Jan 2013
Posts: 79
Thanks for the answer and sorry for being such a novice, but now when trying to minimize I get this error:

CHARMM> energy

NONBOND OPTION FLAGS:
ELEC VDW ATOMs CDIElec FSWItch VATOm VSWItch
BYGRoup NOEXtnd NOEWald
CUTNB = 16.000 CTEXNB =999.000 CTONNB = 10.000 CTOFNB = 12.000
CGONNB = 0.000 CGOFNB = 10.000
WMIN = 1.500 WRNMXD = 0.500 E14FAC = 1.000 EPS = 1.000
NBXMOD = 5
There are 946 atom pairs and 264 atom exclusions.
There are 0 group pairs and 23 group exclusions.
Generating nonbond list with Exclusion mode = 5
== PRIMARY == SPACE FOR 1082 ATOM PAIRS AND 0 GROUP PAIRS

General atom nonbond list generation found:
946 ATOM PAIRS WERE FOUND FOR ATOM LIST
0 GROUP PAIRS REQUIRED ATOM SEARCHES

: No torsion parameters for 72 ( CC3162 CC3161 NG2S3 SC )
: No torsion parameters for 101 ( CC3161 NG2S3 SC OC2DP )
: No torsion parameters for 102 ( CC3161 NG2S3 SC OC2DP )
: No torsion parameters for 103 ( CC3161 NG2S3 SC OC2DP )
: No torsion parameters for 108 ( HCA1 CC3161 NG2S3 SC )
: No torsion parameters for 116 ( HCP1 NG2S3 SC OC2DP )
: No torsion parameters for 117 ( HCP1 NG2S3 SC OC2DP )
: No torsion parameters for 118 ( HCP1 NG2S3 SC OC2DP )
: No torsion parameters for 119 ( SC NG2S3 CC3161 CC3161 )
: A TOTAL OF 9 MISSING PARAMETERS

***** LEVEL -1 WARNING FROM *****
***** CODES> MISSING PARAMETERS
******************************************
BOMLEV ( 0) IS REACHED - TERMINATING. WRNLEV IS 5

May I add those torsions?


Thanks again

Attached Images
in_glc.out.txt (8.72 KB, 770 downloads)
create_disa_nhshb.inp.txt (627 Bytes, 745 downloads)
create.out.txt (16.13 KB, 900 downloads)
Initial-minimization_glc.inp.txt (817 Bytes, 768 downloads)
disa.tar.gz (58.95 KB, 719 downloads)
Joined: Dec 2005
Posts: 1,535
Forum Member
Offline
Forum Member
Joined: Dec 2005
Posts: 1,535
I already thought this went a little bit too easy. Yes, all these torsions are required. And validation/optimization would be highly recommended. Especially the C-C-N-S are quite critical.

Page 1 of 2 1 2

Moderated by  alex, lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.5
(Release build 20201027)
Responsive Width:

PHP: 7.3.31-1~deb10u1 Page Time: 0.013s Queries: 49 (0.008s) Memory: 0.8205 MB (Peak: 0.9402 MB) Data Comp: Off Server Time: 2022-08-16 13:53:10 UTC
Valid HTML 5 and Valid CSS