[gmx-users] Constructing covalent bond between protein-ligand complex

2012-01-13 Thread Suman Nandy
Dear all,
Guided by gmxusers' list I have used gmx_top_tools.tgz 
(http://www.gromacs.org/Downloads/User_contributions/Other_software) to prepare 
the protein-ligand complex linked by a covalent bond. Unfortunately the tool 
renumtop is detecting a duplicate atom but not updating its number. Only 
gives an error message:
; ERROR: duplicate atom label '4968' for atom #4968 (already used for atom 
#4968)

I used the command ./renumtop topol_Protein_chain_A.itp topol_Chain_B.itp. 
may be I am doing something wrong. Please help.

Suman Nandy


Suman Kumar Nandy
Senior Research Fellow
Department of Biochemistry  Biophysics
University of Kalyani
Kalyani-741235
India.-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] Problem with FE-Cys covalent bond

2012-01-11 Thread Suman Nandy
Dear gromacs users,

I have a problem regarding the simulation of Fe4S4 cluster attached to CYS of a 
protein using opls ff. I editted the aminoacids.rtp file as like this.

[ SF4 ]
... [ atoms ]
FE1 opls_966 2.500 1
FE2 opls_967 2.500 1
FE3 opls_968 2.500 1
FE4 opls_969 2.500 1
S1 opls_970 -2.000 1
S2 opls_971 -2.000 1
S3 opls_972 -2.000 1
S4 opls_973 -2.000 1
[ bonds ]
FE1 S1
FE1 S2
FE1 S3
FE2 S1
FE2 S2
FE2 S4
FE3 S2
FE3 S3
FE3 S4
FE4 S1
FE4 S2
FE4 S3

and also atomtypes.atp as

opls_966 55.84500 ; FE1: HEME
opls_967 55.84500 ; FE2: HEME
opls_968 55.84500 ; FE3: HEME
opls_969 55.84500 ; FE4: HEME
opls_970 32.06000 ; S1: S
opls_971 32.06000 ; S2: S
opls_972 32.06000 ; S3: S
opls_973 32.06000 ; S4: S

I have also added a line to specbond.dat regarding the SG-FE covalent bond. 
But while running the pdb2gmx command, (withouit any error) it just adds a H 
with SG (where the FE to bind)and there is no information on SG-FE bond in 
topology file. Please help.

Suman Kumar Nandy
Senior Research Fellow
Department of Biochemistry  Biophysics
University of Kalyani
Kalyani-741235
India.-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] Problem in Disulfide Bond between different chain!

2012-01-11 Thread Suman Nandy

use -chainsep interactive option in pdb2gmx.

Suman

On Wed, 11 Jan 2012 16:42:12 +0530  wrote
Dear all:
I want to do a simulation using a dimer, but the system is connected by 2
disulfide bonds between A-chain and B-chain. I have produced topology file for 
protein with pdb2gmx using AMBER99SB forcefield. In order to bond disulfide 
bonds I using flag -ss. If I do not delete the TER line between A-chain and 
B-chain in .pdb file, the program pdb2gmx can not recognize the disulfide bonds 
between A-chain and B-chain. But if I delete the TER line, I get the two bonds 
successfully. There would be another problem, the C-terminal of A-chain and the 
N-terminal will be bonded in the .top file. How should I do? 
I am using gromacs 4.5.4!
Sincerely yours,
Xianwei Wang






-- 

gmx-users mailing list  gmx-users@gromacs.org

http://lists.gromacs.org/mailman/listinfo/gmx-users

Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!

Please don't post (un)subscribe requests to the list. Use the 

www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Suman Kumar Nandy
Senior Research Fellow
Department of Biochemistry  Biophysics
University of Kalyani
Kalyani-741235
India.-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists