- Posts: 8
Reverse transformation- pdb2gmx error
- c.poojari
- Topic Author
- Offline
- Fresh Boarder
I was successfully able to coarse grain, self assembly of helix structure within the bilayer.
Now i want to convert this CG structure to fine grained structure. I went through the tutorial for reverse transformation and tried the first step:
pdb2gmx -f *.gro -missing
When prompted, i choose GROMOS96 53a6 force field with MAPPING option.
Now i get a error message as such:
Fatal error:
Atom BB in residue MET 1 not found in rtp entry with 9 atoms while sorting atoms
Looking into rtp file, they dont have CG atoms in it.
Please can I know how do i fix this error..
Kind Regards,
Chetan.
Please Log in or Create an account to join the conversation.
- andrzej
- Visitor
did you use coarse grained file as pdb2gmx input ? It should be the atomistic one.
Andrzej
Please Log in or Create an account to join the conversation.
- c.poojari
- Topic Author
- Offline
- Fresh Boarder
- Posts: 8
Yes i had used coarse grained str., after using atomistic it worked fine.
I have another question regarding reconstruction simulation.
Before releasing restraints, my mdp parameters for reconstruction
simulation are;
Restrained SA:
cap_force = 15000
cap_a = 100
fc_restr = 12000
r_CGW = 0.21
fc_restrW = 400
rel_steps = 0
rel_water = no
For Release of restraints,
cap_force = 15000
cap_a = 100
fc_restr = 0
r_CGW = 0.21
fc_restrW = 0
rel_steps = 5000
rel_water = yes
i will be setting rel_steps to 5000,
Should i set rel_water = yes, and set fc_restr = 0, fc_restrW = 0 ??
Kind Regards,
Chetan.
Please Log in or Create an account to join the conversation.
- xavier
- Offline
- Admin
- Posts: 416
c.poojari wrote: Dear Andrzej,
Yes i had used coarse grained str., after using atomistic it worked fine.
I have another question regarding reconstruction simulation.
Before releasing restraints, my mdp parameters for reconstruction
simulation are;
Restrained SA:
cap_force = 15000
cap_a = 100
fc_restr = 12000
r_CGW = 0.21
fc_restrW = 400
rel_steps = 0
rel_water = no
For Release of restraints,
cap_force = 15000
cap_a = 100
fc_restr = 0
r_CGW = 0.21
fc_restrW = 0
rel_steps = 5000
rel_water = yes
i will be setting rel_steps to 5000,
Should i set rel_water = yes, and set fc_restr = 0, fc_restrW = 0 ??
Kind Regards,
Chetan.
Sorry for the delay in answering but we had a problem with the server ...
did you manage to fix the problem you had?
XAvier.
Please Log in or Create an account to join the conversation.
- c.poojari
- Topic Author
- Offline
- Fresh Boarder
- Posts: 8
I am still stuck with it.
Regards,
Chetan
Please Log in or Create an account to join the conversation.
- xavier
- Offline
- Admin
- Posts: 416
XAvier.
Please Log in or Create an account to join the conversation.
- c.poojari
- Topic Author
- Offline
- Fresh Boarder
- Posts: 8
I am running protein-membrane reverse transformation simulation.
I ran restrained simulated annealing with below parameters, run completed successfully. I ran this at 1 fs time step for 40 ps.
cap_force = 15000
cap_a = 100
fc_restr = 12000
r_CGW = 0.21
fc_restrW = 400
rel_steps = 0
rel_water = no
Now for release of restraints, i have used the following parameters. Even this run completed without any errors. I ran this at 1 fs time step for 40 ps.
cap_force = 15000
cap_a = 100
fc_restr = 12000
r_CGW = 0.21
fc_restrW = 400
rel_steps = 5000
rel_water = no
In total i ran reverse tranformation for 80 ps. Please can i know if my settings for release of restraints is alright or should i be setting rel_water= yes.
Kind Regards,
Chetan.
Please Log in or Create an account to join the conversation.