[gmx-users] Regarding indexing the atom labelled same in peptide
Dilip H N
cy16f01.dilip at nitk.edu.in
Tue Aug 22 18:49:49 CEST 2017
Hello,
I have a Diglycine peptide in .gro file as:-
1GLY N 1 1.312 0.940 1.100 0.5173 0.1852 0.2665
1GLY H1 2 1.375 1.015 1.064 0.3308 -0.2009 -0.8897
1GLY H2 3 1.376 0.860 1.115 2.0625 0.7372 -2.9245
1GLY H3 4 1.242 0.907 1.031 0.2385 0.6317 0.3326
1GLY CA 5 1.240 0.963 1.228 0.3801 -0.0854 0.2381
1GLY HA1 6 1.155 0.898 1.233 1.2980 -1.3807 -0.5280
1GLY HA2 7 1.321 0.963 1.300 -1.1296 -1.6137 1.9682
1GLY C 8 1.185 1.101 1.232 0.0318 -0.2122 -0.1516
1GLY O 9 1.171 1.174 1.133 0.1199 0.2166 0.1513
2GLY N 10 1.148 1.148 1.352 -1.0777 -0.3721 -0.4264
2GLY HN 11 1.142 1.091 1.434 -0.1054 0.0667 -0.0362
2GLY CA 12 1.092 1.277 1.377 0.5476 0.2162 0.2355
2GLY HA1 13 1.170 1.351 1.365 0.3836 0.3441 -0.0427
2GLY HA2 14 1.012 1.294 1.306 -0.5941 1.1037 1.7136
2GLY C 15 1.031 1.298 1.514 0.2354 -0.4348 0.1971
2GLY OT1 16 0.923 1.238 1.541 0.2601 -0.2386 0.7420
2GLY OT2 17 1.091 1.369 1.598 0.1878 -0.0650 -0.0813
So now i want to index of the Nitrogen (N) of 1GLY only, but not N of
2GLY...
I tried the the following command as:-
gmx make_ndx -f abc.gro -o def.ndx
1 & a N
But it is giving 2 Protein_&_N : 2 atoms
which is taking the N of 2GLY too into consideration....which i don't
required, I need the N of ...1GLY only....
So what is the correct command for this..??
Thank you...
With Best Regards,
DILIP.H.N
Ph.D Student
<https://mailtrack.io/> Sent with Mailtrack
<https://mailtrack.io/install?source=signature&lang=en&referral=cy16f01.dilip@nitk.edu.in&idSignature=22>
More information about the gromacs.org_gmx-users
mailing list