Previous Thread
Next Thread
Print Thread
IC Edit not modifying coordinates
#38367 11/09/21 06:36 AM
Joined: Dec 2017
Posts: 3
A
Anmol Offline OP
Forum Member
OP Offline
Forum Member
A
Joined: Dec 2017
Posts: 3
I am trying to change a particular dihedral for a known set of coordinates by editing the IC. The IC changes but coordinate remains the same.
I could not figure out the mistake in following code.

Code
stream toppar.str
stream @str
set residue @resi

read sequence @resi 1
generate @resi first none last none setup warn

read coor card name @resicrd

ic fill ! preserve

set value -180

label repeat
ic edit
  dihedral 1 N1 1 C3 1 C4 1 C5 @value
end
ic print
coor init sele (type N1 .or. type C3 or. type C4 .or. type C5) end
ic build

open unit 30 write form name @resi_@value.pdb
write coor pdb unit 30
close unit 30
incr value by 15
if @value eq 180 then
   stop
else
   goto repeat
endif


Re: IC Edit not modifying coordinates
Anmol #38368 11/09/21 03:21 PM
Joined: Sep 2003
Posts: 8,582
Likes: 11
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,582
Likes: 11
The script is not sufficient for diagnosis; there are probably some major hints in the output log file and .pdb files.

For changing a dihedral via IC EDIT, COOR INIT, and IC BUILD, one way is to apply COOR INIT to all atoms on only one side of the dihedral, e.g. N1, C3, and all of atoms connected to those atoms. For a small molecule, you could probably use COOR INIT on all of the atoms. Also, there could be issues if the 4 atoms listed in IC EDIT are not already part of an existing IC entry from the topology file.


Rick Venable
computational chemist

Re: IC Edit not modifying coordinates
rmv #38369 11/10/21 02:26 AM
Joined: Dec 2017
Posts: 3
A
Anmol Offline OP
Forum Member
OP Offline
Forum Member
A
Joined: Dec 2017
Posts: 3
Thanks for your reply. The same script works now. As you pointed out, the dihedral mentioned had wrong atom order, thus was not part of existing IC and was not being used to generate the coordinates.


Moderated by  lennart, rmv 

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

PHP: 7.3.31-1~deb10u1 Page Time: 0.009s Queries: 20 (0.006s) Memory: 0.7230 MB (Peak: 0.7633 MB) Data Comp: Off Server Time: 2021-12-07 23:44:59 UTC
Valid HTML 5 and Valid CSS