Re: [gmx-users] Justin umbrella sampling tutorial......

2012-05-24 Thread rama david
Thank you for your reply, I am asking you again same question, EXTREMELY SORRY for my stupidity, In step six , I unable to differentiate npt and production run by mdp file as usualy we find difference by define term, I think I get meaning upto reason why to use -DPOSRES_B, but I want to know if

[gmx-users] Re: Regarding error.

2012-05-24 Thread vivek sharma
Dear Suryanarayana, This error itself tells you that the particular residue 'DAL' is not available in the residue topology database of the force field you are using for your simulation. You can check that by yourself in the *.rtp file of corresponding force field. The way out of this situation is

Re: [gmx-users] g_rms -bm

2012-05-24 Thread Kowsar Bagherzadeh
Dear Justin   Thank you very muchh   Sogol From: Justin A. Lemkul jalem...@vt.edu To: Kowsar Bagherzadeh kw_bagherza...@yahoo.com; Discussion list for GROMACS users gmx-users@gromacs.org Sent: Thursday, May 24, 2012 9:59 AM Subject: Re: [gmx-users] g_rms -bm

Re: [gmx-users] Justin umbrella sampling tutorial......

2012-05-24 Thread Justin A. Lemkul
On 5/24/12 8:19 AM, rama david wrote: Thank you for your reply, I am asking you again same question, EXTREMELY SORRY for my stupidity, In step six , I unable to differentiate npt and production run by mdp file as usualy we find difference by define term, I think I get meaning upto reason

Re: [gmx-users] Justin umbrella sampling tutorial......

2012-05-24 Thread rama david
Thank you Justin for these correct explanation Its really clear my lot of queries.. For the tutorial, NPT is conducted with restraints on all protein heavy atoms. The production runs are conducted by restraining only one chain for practical reasons. These is my question ; If we are

[gmx-users] Distance Restraints of 3 atoms

2012-05-24 Thread Steven Neumann
Dear Gmx Users, I want use distance restraints of 3 atoms belonging to 3 different moleculetypes. In this case shall I create and index file of those atoms and then use: genrestr -f X.pdb -n index.ndx -disre -o disre.itp ? Under which topology file shall I add: #ifdef DISRES # include

[gmx-users] Protein-ligand interactions in implicit water

2012-05-24 Thread Андрей Гончар
Hello. Is it possible to perform MD for protein-ligand interactions in implicit water? And what about the precision of this kind of MD's? Are they differ strong from these with explicit water? -- Andrew Gontchar -- gmx-users mailing listgmx-users@gromacs.org

[gmx-users] trjconv_snapshot

2012-05-24 Thread Turgay Cakmak
Hi Gromacs users, I have a problem getting a snapshots using trjconv. If you could help me, I would be very grateful. Below, I try to explain what I did and problem happened. Firstly, I have done 10ns long simulation (several peptides in a box). And then, using the outputs of that simulation

Re: [gmx-users] trjconv_snapshot

2012-05-24 Thread rama david
On Thu, May 24, 2012 at 5:24 PM, Turgay Cakmak turgaycakma...@gmail.comwrote: Hi Gromacs users, I have a problem getting a snapshots using trjconv. If you could help me, I would be very grateful. Below, I try to explain what I did and problem happened. Firstly, I have done 10ns long

[gmx-users] trjconv_snapshot

2012-05-24 Thread Turgay Cakmak
Hi Rama, Thanks for your quick reply. Catenate two trajectory with help of trjcat -h to get snapshot of particular time use trjconv -dump time use appropriate pbc option I have also used trjconv with -dump and -pbc nojump options, but still I have same problem. you are using -b 1 -e 1

Re: [gmx-users] trjconv_snapshot

2012-05-24 Thread rama david
Hi , You run extended simulation of 10 ns on your first 10ns simulation . You get traj20ns.xtc. these new trajectory contain all the tme frame from onwards 10ns. So if you are giving command trjconv –f traj20ns.xtc –s topol.tpr –b 1 –e 1 –o snapshot20ns.gro (Note: traj20ns.xtc is

[gmx-users] One molculetype for 3 proteins

2012-05-24 Thread Steven Neumann
Dear Gmx Users, I want to use distance restraints of 3 atoms which belong to 3 different proteins. For this purpose I want to create one topology so one moleculetype: I tried to use pdb2gmx -chainsep ter (where i put TER after each chain in my PDB) But 3 separate topologies were created. I

[gmx-users] Re: One molculetype for 3 proteins

2012-05-24 Thread Steven Neumann
On Thu, May 24, 2012 at 3:42 PM, Steven Neumann s.neuman...@gmail.comwrote: Dear Gmx Users, I want to use distance restraints of 3 atoms which belong to 3 different proteins. For this purpose I want to create one topology so one moleculetype: I tried to use pdb2gmx -chainsep ter (where i

[gmx-users] CHARMM-27 CHARMM-36

2012-05-24 Thread Shima Arasteh
Dear gmx users, ACE is defined in CHARMM-27 .rtp file, but I don't see this residue in CHARMM-36 .rtp file . Is it OK, if I add the ACE from CHARMM-27 into CHARMM-36? And use the CHARMM-27 [ACE] parameters in CHARMM-36?   Thanks in advance. Sincerely, Shima-- gmx-users mailing list

Re: [gmx-users] Re: One molculetype for 3 proteins

2012-05-24 Thread Tsjerk Wassenaar
Hi Steven, Doesn't -chainsep take an argument? Cheers, Tsjerk On May 24, 2012 5:33 PM, Steven Neumann s.neuman...@gmail.com wrote: On Thu, May 24, 2012 at 3:42 PM, Steven Neumann s.neuman...@gmail.com wrote: Dear Gmx Users,... It is stated: The separation of chains is not entirely trivial

[gmx-users] Equilibrating Water around the Molecule

2012-05-24 Thread Lara Bunte
Hello I want to equilibrate the water around my molecule. For doing this I created a pr.mdp file that makes most sense for me. Sadly I am a beginner with less knowledge, so could you please take a look on this pr.mdp file and say me if this is good or if not, what would you suggest to me? My

Re: [gmx-users] Re: One molculetype for 3 proteins

2012-05-24 Thread Steven Neumann
On Thu, May 24, 2012 at 6:05 PM, Tsjerk Wassenaar tsje...@gmail.com wrote: Hi Steven, Doesn't -chainsep take an argument? Cheers, Tsjerk It does, I used option: -chainsep ter And I had in my pdb file after each chain TER and thus 3 itp files were created corresponding to my top file and

[gmx-users] (no subject)

2012-05-24 Thread Pathumwadee Intharathep
pa href=http://tgiturkiye.net/breakingnews/63LeePhillips/;http://tgiturkiye.net/breakingnews/63LeePhillips//a/p -- 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

Re: [gmx-users] Re: One molculetype for 3 proteins

2012-05-24 Thread Tsjerk Wassenaar
Hi Steven, So it did what you told it to do. You want -chainsep no or none. Doing it from the top of my head, so I'm not sure. But it'll be listed in the help what options there are. Cheers, Tsjerk On Thu, May 24, 2012 at 8:29 PM, Steven Neumann s.neuman...@gmail.com wrote: On Thu, May 24,

[gmx-users] Re: One molculetype for 3 proteins

2012-05-24 Thread Francesca
I had the same problem and I didn't find any way to do it. So i created new residues for the terminals, and I added them to the forcefields files. Example: if my first COOH terminal was a SER i created a SERT residues and I used the COOH values taken from ASP , and the same procedure was used for

Re: [gmx-users] Re: One molculetype for 3 proteins

2012-05-24 Thread Tsjerk Wassenaar
Okay,okay, I give in. To settle this again, I looked it up, and among the options to -chainsep is 'interactive'. If you use that you're prompted whether or not to separate at the breaks pdb2gmx detects. That allows you to merge chains. I guess there'd be sense in just having an option 'no' to

Re: [gmx-users] Re: One molculetype for 3 proteins

2012-05-24 Thread Steven Neumann
Tahnk you both. Indeed, it works so -chainsep interactive Steven On Thu, May 24, 2012 at 9:44 PM, Tsjerk Wassenaar tsje...@gmail.com wrote: Okay,okay, I give in. To settle this again, I looked it up, and among the options to -chainsep is 'interactive'. If you use that you're prompted whether

[gmx-users] Re: One molculetype for 3 proteins

2012-05-24 Thread Francesca
I checked and the resul is the same...but now I know I can use -chainsep interactive! -- View this message in context: http://gromacs.5086.n6.nabble.com/One-molculetype-for-3-proteins-tp4997727p4997739.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users

Re: [gmx-users] CHARMM-27 CHARMM-36

2012-05-24 Thread Felipe Merino
Sure, CHARMM 36 is just a modification to the lipids. The protein (and other sutff) parameters are untouched. regards Felipe Mensaje original De: shima_arasteh2...@yahoo.com Fecha: 24-may-2012 12:41 Para: Discussion list for GROMACS usersgmx-users@gromacs.org Asunto:

[gmx-users] ACE residues at positions other than resnr=1 possibly leading to non-sequential residue numbering by editconf and trjconv

2012-05-24 Thread Christopher Neale
Dear Users: Summary: I believe that there is some problem in the residue numbering routine that is used by editconf and trjconv that exists at least in 4.5.3 and 4.5.5, but not in 4.0.7, and that this is somehow related to the ACE residue. Details: I have a system with two peptides, water, and

Re: [gmx-users] Justin umbrella sampling tutorial......

2012-05-24 Thread Justin A. Lemkul
On 5/24/12 6:35 AM, rama david wrote: Thank you Justin for these correct explanation Its really clear my lot of queries.. For the tutorial, NPT is conducted with restraints on all protein heavy atoms. The production runs are conducted by restraining only one chain for practical

Re: [gmx-users] trjconv_snapshot

2012-05-24 Thread Justin A. Lemkul
On 5/24/12 10:17 AM, rama david wrote: Hi , You run extended simulation of 10 ns on your first 10ns simulation . You get traj20ns.xtc. these new trajectory contain all the tme frame from onwards 10ns. So if you are giving command trjconv –f traj20ns.xtc –s topol.tpr –b 1 –e 1 –o

Re: [gmx-users] Selectively minimze Hydrogens

2012-05-24 Thread Justin A. Lemkul
On 5/24/12 2:19 PM, Kartheek wrote: Hi all, I am trying to minimize a system containing DNA+protein, I just want to energy minimize hydrogen atoms for 500 steps , while the rest of the structure is kept fixed. How can i do it in gromacs. Please suggest me a way. Have a look