[gmx-users] Re:Re:Why REMD simulation becomes so slow when the number of replicas becomes large?

2011-02-08 Thread Qiong Zhang
Hi Mark,   Your analyses are quite reasonable. The low-temperature replicas are indeed doing much more work than the high-temperature replicas. As you said, the lowest temperature replica in the 24-replica should take an amount of time comparable to that of the lowest in the 42-replica. So for

[gmx-users] Increasing machine precision

2011-02-08 Thread abdullah ahmed
Hello, Is there a way to continue minimizing after reaching machine precision? Emtol and the number of iterations are sufficient to continue. I am assuming that reaching machine precision means that the gradient of change from one iteration to another has become so small that further

Re: [gmx-users] Increasing machine precision

2011-02-08 Thread Erik Marklund
Hello, Is there a way to continue minimizing after reaching machine precision? Emtol and the number of iterations are sufficient to continue. I am assuming that reaching machine precision means that the gradient of change from one iteration to another has become so small that further

[gmx-users] Re: No residue type for 'ARG' as a starting terminus

2011-02-08 Thread bipin singh
Hi all, I am getting following during while running pdb2gmx for a RNA moleculei am using amber99sb force field parameters The details of the error is as: Fatal error: In the chosen force field there is no residue type for 'ARG' as a starting terminus For more information and tips for

[gmx-users] OPLS and RB parameters in GROMACS

2011-02-08 Thread sulatha M. S
Hi all, Hi I've converted the OPLS-AA torsional potential parameters for the alkane C-C-C-C (CT-CT-CT-CT in the gromacs parameter file notation), C-C-C-H (CT-CT-CT-HC), and H-C-C-H (HC-CT-CT-HC) torsions from the OPLS format given in Jorgensen et al, JACS 118, 11225 (1996) to the

[gmx-users] Re:gmx-users Digest, Vol 82, Issue 56

2011-02-08 Thread gromacs
Hi,guys I obtained surface tension of water and salt solution. But the value RMSD was always very big although the value of 65 is reasonable. Is it all right about the RMSD?? Energy Average Err.Est. RMSD Tot-Drift

[gmx-users] sufrace tension RMSD too big

2011-02-08 Thread gromacs
Hi,guys I obtained surface tension of water and salt solution. But the value RMSD was always very big although the value of 65 is reasonable. Is it all right about the RMSD?? Energy Average Err.Est. RMSD Tot-Drift

Re: [gmx-users] OPLS and RB parameters in GROMACS

2011-02-08 Thread Andrew Paluch
There is no error. The alkane dihedral parameters were updated in 1999, and differ from those originally published in 1996. Andrew On Tue, Feb 8, 2011 at 5:20 AM, sulatha M. S mssula...@gmail.com wrote: Hi all, Hi I've converted the OPLS-AA torsional potential parameters for the alkane

[gmx-users] Reduced Units

2011-02-08 Thread Tomy van Batis
Dear all I would like to do coarse-grained simulations using reduced LJ units. As I can see in the manual, this is possible, but I don't understand how... I use the Gromos96 53a6 force field. In the *ffG53a6nb.itp *the* σ* and *ε* are in *(nm*) and *(KJ/mol*) Also in the *ffG53a6.atp* the

Re: [gmx-users] Re: free energy calculation , grompp crash

2011-02-08 Thread Da-Wei Li
hello Now it is almost clear what happened. When couple-intramol is no (default), all pairwise vdm and charge interaction becomes bonded interaction. All intra-molecular non-bonded interactions for moleculetype couple-moltype are replaced by exclusions and explicit pair interactions. In this

[gmx-users] Nice Present

2011-02-08 Thread boaz kan-tor
This was very useful for me http://www.sunsetstripbillboards.com/info.html Wanna share with you. -- Boaz Kan-Tor I am a humanist, which means, in part, that I have tried to behave decently without expectations of rewards or punishments after I am dead. -- gmx-users mailing list

[gmx-users] ask for help on MARTINI SIMULATION OF POPC lipids with polarisable water

2011-02-08 Thread wezhao
Dear All, I have performed a simulation of POPC with polarisable water model in presence of 0.2 mol CaCl2 based on MARTINI CG model. I used shift for the electrostatic interactions and the job was done on a 8-core node with domain decomposition 2 2 2. The lipid bilayer consists of 512 lipid and

[gmx-users] Help with B2AR within the POPC membrane

2011-02-08 Thread Aldo Segura
Dear gmxusers, I need to perform molecular dynamics simulation of a B2AR within the POPC membrane. I have downloaded the 128b.pdb, popc.itp and lipid.itp files from Prof.Tieleman's group. My protein of interest is 343 residues. Also, I aligned the protein and membrane. I followed the Justin

Re: [gmx-users] Help with B2AR within the POPC membrane

2011-02-08 Thread Justin A. Lemkul
Aldo Segura wrote: Dear gmxusers, I need to perform molecular dynamics simulation of a B2AR within the POPC membrane. I have downloaded the 128b.pdb, popc.itp and lipid.itp files from Prof.Tieleman's group. My protein of interest is 343 residues. Also, I aligned the protein and membrane. I

Re: [gmx-users] Help with B2AR within the POPC membrane

2011-02-08 Thread Aldo Segura
Thanks for your answer. You're right in the procedure for the packing of the protein and lipids. However, after several iterations (~30) the lipids are packaged to form the bilayer and the protein is outside of it. I can send you a couple of pictures for a better explanation of my problem.

Re: [gmx-users] Help with B2AR within the POPC membrane

2011-02-08 Thread Justin A. Lemkul
Aldo Segura wrote: Thanks for your answer. You're right in the procedure for the packing of the protein and lipids. However, after several iterations (~30) the lipids are packaged to form the bilayer and the protein is outside of it. I can send you a couple of pictures for a better explanation

Re: [gmx-users] Help with B2AR within the POPC membrane

2011-02-08 Thread Aldo Segura
Dear Justin, You're right, I corrected the box vectors , and it works! Thanks, -- 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

[gmx-users] Re: g_rmsf

2011-02-08 Thread bharat gupta
Hi, I want to calculate the RMSF of residues and not of protein ... how can this be done with g_rmsf.. Also I want to see the rmsf of certain residues .. for which I created the .ndx file containint those residues only .. and after using g_rmsf with index file gives the RMSF for whole protein

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread Justin A. Lemkul
bharat gupta wrote: Hi, I want to calculate the RMSF of residues and not of protein ... how can this be done with g_rmsf.. Also I want to see the rmsf of certain residues .. for which I created the .ndx file containint those residues only .. and after using g_rmsf with index file gives

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread bharat gupta
Actually after loop incorporation I want to check which region of the protein shows much deviation , which I think can be done by plotting rmsf values from both proteins.. but the problem here is that one structure which contains loops has more no. of atoms as compared to other str. without loop

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread Justin A. Lemkul
bharat gupta wrote: Actually after loop incorporation I want to check which region of the protein shows much deviation , which I think can be done by plotting rmsf values from both proteins.. but the problem here is that one structure which contains loops has more no. of atoms as compared to

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread bipin singh
try -res option On Wed, Feb 9, 2011 at 08:51, Justin A. Lemkul jalem...@vt.edu wrote: bharat gupta wrote: Actually after loop incorporation I want to check which region of the protein shows much deviation , which I think can be done by plotting rmsf values from both proteins.. but the

Re: [gmx-users] solvation_box_preparation

2011-02-08 Thread shahid nayeem
Hi Justin Thanks a lot. I tried doing energy minimization and th lowering emtotal to 200 and the system converged to Steepest Descents converged to Fmax 200 in 1411 steps Potential Energy = -3.9063050e+05 Maximum force = 1.3442458e+02 on atom 927 Norm of force = 1.4758101e+01 For

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread bharat gupta
I used the -res option ... and I got the rmsf in terms of residues but still the problem is that the two structures contain different amount of residues due to loop replacement in one structure.. In that case how shall proceed to check the effect of loop insertion on the overall topology of the

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread Tsjerk Wassenaar
Hi Bharat, You can do it with post-processing the data you obtain from g_rmfs, if it's okay that the fit uses all residues in either case. Otherwise, you can make an two index files, including only the residues that are common to both. Hope it helps, Tsjerk On Wed, Feb 9, 2011 at 5:58 AM,

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread bipin singh
You calculate rmsf of both proteins separately and then plot them together and look the region of your interest.. On Wed, Feb 9, 2011 at 10:28, bharat gupta bharat.85.m...@gmail.com wrote: I used the -res option ... and I got the rmsf in terms of residues but still the problem is that the

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread bharat gupta
Thanks for the advice and while creating the index file for first 100 common residues it found that both structures shows different no.of atoms. .. how is that possible ?? On Tue, Feb 8, 2011 at 9:00 PM, Tsjerk Wassenaar tsje...@gmail.com wrote: Hi Bharat, You can do it with post-processing

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread Tsjerk Wassenaar
We can't know that. We don't have your files (and we don't want them). Check what you did, and what index groups you have. Write out the structure for the index group and have a look. We're not a substitute for your brain here... :p Cheers, Tsjerk On Wed, Feb 9, 2011 at 6:20 AM, bharat gupta

Re: [gmx-users] Re: No residue type for 'ARG' as a starting terminus

2011-02-08 Thread Mark Abraham
On 8/02/2011 9:27 PM, bipin singh wrote: Hi all, I am getting following during while running pdb2gmx for a RNA moleculei am using amber99sb force field parameters The details of the error is as: Fatal error: In the chosen force field there is no residue type for 'ARG' as a starting

Re: [gmx-users] OPLS and RB parameters in GROMACS

2011-02-08 Thread sulatha M. S
I looked at the the paper published in 1999, ( Jorgensen et al JACS, 121, 20, 4831, 1999), the aliphatic torsional parameters are the same as those from 1996. here are the values, V1 V2V3 (kcal/mol) CT-CT-CT-CT 1.740

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread bharat gupta
Tsjerk, Sorry for asking that .. actually I made a silly mistake while selecting residues.. After plotting the graphs for common regions I found that first 100 amino acids shows a lot of fluctuations (compared to the one without any loop insertion) ... Does the insertion caused a great change in

Re: [gmx-users] Re: No residue type for 'ARG' as a starting terminus

2011-02-08 Thread bipin singh
Sir, Actually ARG is present as a ligand bound to RNA molecule On Wed, Feb 9, 2011 at 11:16, Mark Abraham mark.abra...@anu.edu.au wrote: On 8/02/2011 9:27 PM, bipin singh wrote: Hi all, I am getting following during while running pdb2gmx for a RNA moleculei am using amber99sb force

Re: [gmx-users] OPLS and RB parameters in GROMACS

2011-02-08 Thread Mark Abraham
On 9/02/2011 4:48 PM, sulatha M. S wrote: I looked at the the paper published in 1999, ( Jorgensen et al JACS, 121, 20, 4831, 1999), the aliphatic torsional parameters are the same as those from 1996. here are the values, V1 V2

Re: [gmx-users] Re: g_rmsf

2011-02-08 Thread Mark Abraham
On 9/02/2011 4:52 PM, bharat gupta wrote: Tsjerk, Sorry for asking that .. actually I made a silly mistake while selecting residues.. After plotting the graphs for common regions I found that first 100 amino acids shows a lot of fluctuations (compared to the one without any loop insertion)

Re: [gmx-users] Re: No residue type for 'ARG' as a starting terminus

2011-02-08 Thread bipin singh
Thanks for your help. On Wed, Feb 9, 2011 at 11:44, Mark Abraham mark.abra...@anu.edu.au wrote: On 9/02/2011 4:56 PM, bipin singh wrote: Sir, Actually ARG is present as a ligand bound to RNA molecule Then you've got work to do. pdb2gmx copes well with linear polymers of

Re: [gmx-users] OPLS and RB parameters in GROMACS

2011-02-08 Thread sulatha M. S
Sorry for not giving the manual version earlier. P. 63 of manual 4.0 is what I am referring to. The equations relate to converting the OPLS torsional parameters to RB parameters in GROMACS with the OPLS ff. The equations are C0 = V0+V2+0.5(V1+V3) C1= 0.5(3 * V3-V1) C2= -V2 + 4 * V4 C3= -2 * V3

[gmx-users] add force field

2011-02-08 Thread lina zhao
Hi, Anyone has successfully tried to add some new force field in the default force field (in 4.5) so pdb2gmx can show up the (newly-added) force field which will also be recognized by gromacs, Thanks, lina -- gmx-users mailing listgmx-users@gromacs.org