IC Edit not modifying coordinates

Print Thread

#38367 11/09/21 06:36 AM

Joined: Dec 2017

Posts: 3

Anmol

Forum Member

Anmol

Forum Member

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,650 Likes: 26 39 59 56 N, 77 06 07 W rmv Forum Member
rmv Forum Member Joined: Sep 2003 Posts: 8,650 Likes: 26 39 59 56 N, 77 06 07 W	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 OP Forum Member
OP Anmol 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

CHARMM Development Project Forums User Discussion & Questions Setup, I/O, and Basic questions IC Edit not modifying coordinates