Previous Thread
Next Thread
Print Thread
Page 2 of 2 1 2
Re: NMA parameters
alex #7999 08/29/05 02:53 PM
Joined: Oct 2004
Posts: 71
F
Forum Member
OP Offline
Forum Member
F
Joined: Oct 2004
Posts: 71
Dear Alex,

Thanks for your reply.
Our group is now working on QM calculations for NMA. I will report something on this issue later.

Re: NMA parameters
alex #8000 09/01/05 03:47 PM
Joined: Oct 2004
Posts: 71
F
Forum Member
OP Offline
Forum Member
F
Joined: Oct 2004
Posts: 71
Dear Alex,

I have a question to modify the dihedral force parameters. As you know, I want to change the dihedral parameters for the methyl group in NMA (or any peptide), and I changed the parameter file from

HA CT3 NH1 H 0.00 3 0.0

to

HA CT3 NH1 H 0.10 3 0.0

and I calculated the dihedral energy as a function of the methyl group rotation angle (while other degrees of freedom fixed). I got the three fold dihedral interaction as expected. The barrier height is 0.6 kcal/mol because there are three hydrogens.

Then I changed the delta=0.0 to, say, 180.0, that is, from

HA CT3 NH1 H 0.10 3 0.0

to

HA CT3 NH1 H 0.10 3 180.0

and I expected 60 degree shift for the potential curve, but the result is about 12.5 degree shift. How can I understand this?

I appreciate your comments.

Re: NMA parameters
fujisaki #8001 09/03/05 06:28 PM
Joined: Oct 2003
Posts: 782
A
Forum Member
Offline
Forum Member
A
Joined: Oct 2003
Posts: 782

1) be sure that there is not an additional dihedral parameter with a different multiplicity acting on the dihedral you are studying.

2) the combination of electrostatic, LJ and dihedral terms are contributing to the energy surface, that the obtained value of the minimum is due to the combination of these terms. do an analysis of the individual energy terms as a function of the dihedral angle to understand in better detail what is going on. the substitution options can write out the individual terms as a function of the dihedral.

alex


School of Pharmacy
University of Maryland
20 Penn Street
Baltimore, MD, 21201
Re: NMA parameters
alex #8002 09/04/05 12:57 AM
Joined: Oct 2004
Posts: 71
F
Forum Member
OP Offline
Forum Member
F
Joined: Oct 2004
Posts: 71
Dear Alex,

I checked what you said, but I still can't solve this problem.
Could you look at my input code, and give me some advice?

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

* NMA in vacuum
*

BOMBLEV -1

! Topology and parameters

open unit 1 card read name top_all22_model.inp
read RTF card unit 1
close unit 1

open unit 1 card read name par_all22_prot.inp
read PARA card unit 1
close unit 1

!
! read in sequence
read sequ card
* sequence for NMA
*
1
NMA

generate nma setup

! build in missing coordinates using values in
! the parameter set
ic para
ic seed 1 CR 1 N 1 C
ic build

coor init sele type H* end
hbuild

print coor

mini conj nstep 1000 Inbfrq 10 CUTNB 13.0 CTONNB 8.0 CTOFNB 12.0 -
group vdw fswitch vswitch

coor axis sele type CR end sele type N end

set j 1
label loop

coor rotate axis phi @j sele type HR1 .or. type HR2 .or. type HR3 end
energy

!open write formatted unit 10 name scratch/NMA-rot@j.pdb
!write coor pdb unit 10
!* Rotated structure
!*
!close unit 10

coor rotate axis phi -@j sele type HR1 .or. type HR2 .or. type HR3 end

incr j by 1
if j lt 181 goto loop

stop

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

I just changed the parameter files as follows:

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

!HA CT3 NH1 C 0.0000 3 0.00 ! ALLOW PEP
HA CT3 NH1 C 0.1000 3 120.00 ! ALLOW PEP
! LK for autogenerate dihe, sp2-methyl, no dihedral potential

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

I expected that the dihedral interaction energy displayed in "ENER INTERN>" as a function of the dihedral angle should shift by 60 degrees, but it results in 16.5 degrees shift.

I appreciate your reply.

Re: NMA parameters
fujisaki #8003 09/04/05 07:02 PM
Joined: Oct 2003
Posts: 782
A
Forum Member
Offline
Forum Member
A
Joined: Oct 2003
Posts: 782

below is an example for calculating the surface; note the use of "cons dihe" to sample the dihedral surface. keep the phase of the dihedral at 0 or 180.

alex

set n 999.
set 3 @n ! cutim
decr n by 1
set 4 @n ! cutnb
decr n by 3
set 5 @n ! ctonnb
incr n by 2
set 6 @n ! ctofnb
set 7 switch
set 8 atom
set 9 vatom

update inbfrq 20 ihbfrq 0 imgfrq 20 -
@7 @8 vswitch @9 cutim @3 cutnb @4 ctonnb @5 ctofnb @6

mini abnr nstep 100 nprint 10
mini nraph nstep 50 nprint 5 tolgrd 0.0

! set minimimum values to offset surface to 0.0
set a ?ener
set b ?elec
set c ?vdw
set d ?dihe
set e ?angl
set f ?bond
set g ?urey

open unit 99 write form name butane_cccc_@id2.surf

set dihe 180.

label loop

cons dihe but 1 c1 but 1 c2 but 1 c3 but 1 c4 forc 10000.0 min @dihe

mini abnr nstep 100 nprint 10 inbfrq 100
mini nraph nstep 20 nprint 10 tolgrd 0.0
cons cldh

energy
! offset by global minimum energy
set 5 ?ener
decr 5 by @a
set 6 ?elec
decr 6 by @b
set 7 ?vdw
decr 7 by @c
set 8 ?dihe
decr 8 by @d
set 9 ?angl
decr 9 by @e
set 0 ?bond
decr 0 by @f
set 4 ?urey
decr 4 by @g

write title unit 99
* @dihe @5 @6 @7 @8 @9 @0 @4
*

decr dihe by 5.

if dihe ge 0.0 goto loop


School of Pharmacy
University of Maryland
20 Penn Street
Baltimore, MD, 21201
Re: NMA parameters
alex #8004 09/05/05 12:05 AM
Joined: Oct 2004
Posts: 71
F
Forum Member
OP Offline
Forum Member
F
Joined: Oct 2004
Posts: 71
Dear Alex,

Thanks for your reply and information.

First off I have to say that I misunderstood in the previous post. I looked at the minimized structure with the following parameter file change

!HA CT3 NH1 C 0.0000 3 0.00 ! ALLOW PEP
HA CT3 NH1 C 0.1000 3 120.00 ! ALLOW PEP
! LK for autogenerate dihe, sp2-methyl, no dihedral potential

and the methyl group in NMA certainly shifts by 120/3=40 degrees. But the problem is if I shift further the methyl group rotation angle, this is not the stable structure with respect to the dihedral interaction. I think this is strange. This does not occur if I choose 0 degree or 180 degrees, but it does for other dihedral shift angles. How can I understand this?

Re: NMA parameters
fujisaki #8005 09/05/05 12:29 PM
Joined: Oct 2003
Posts: 782
A
Forum Member
Offline
Forum Member
A
Joined: Oct 2003
Posts: 782

look at the individual energy components and you will figure it out. note that you can use the INTER command with a selection of 4 atoms to output the dihedral energy for those 4 atoms.

alex


School of Pharmacy
University of Maryland
20 Penn Street
Baltimore, MD, 21201
Re: NMA parameters
alex #8006 09/05/05 06:39 PM
Joined: Oct 2004
Posts: 71
F
Forum Member
OP Offline
Forum Member
F
Joined: Oct 2004
Posts: 71
Dear Alex,

THanks for your reply.

I used

INTER sele type N .or. type CR end

instead of just "energy" command in the code of post #7789, and I looked at INTE INTERN lines, but the result turns out to be the same as the ENER INTERN lines with the energy command. So this minimized structure is not stable when I rotate the methyl group. I would like to know what happens when the dihedral parameter in the CHARMM parameter file changes from 0 or 180. It seems that only 0 and 180 degrees are allowed in CHARMM, and if we use others, some error might occur ...

I appreciate your comment.

Re: NMA parameters
fujisaki #8007 09/06/05 12:12 PM
Joined: Oct 2003
Posts: 782
A
Forum Member
Offline
Forum Member
A
Joined: Oct 2003
Posts: 782

at this stage I don't think I can be of additional help. you have to go carefully through the documentation and really think about what you are trying to do. our contribution to the forum is to point you in the right direction to solve your problem; we can't solve all the details of the problem.

alex


School of Pharmacy
University of Maryland
20 Penn Street
Baltimore, MD, 21201
Re: NMA parameters
alex #8008 09/08/05 05:39 PM
Joined: Oct 2004
Posts: 71
F
Forum Member
OP Offline
Forum Member
F
Joined: Oct 2004
Posts: 71
Hi Alex,

Thanks for your reply.
I think I have solved my problem! Thanks for your guidance!

Page 2 of 2 1 2

Moderated by  alex, lennart, rmv 

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

PHP: 5.6.33-0+deb8u1 Page Time: 0.016s Queries: 35 (0.008s) Memory: 0.9895 MB (Peak: 1.1232 MB) Data Comp: Off Server Time: 2020-10-25 02:10:31 UTC
Valid HTML 5 and Valid CSS