Re: [gmx-users] Salt Bridge

2014-06-12 Thread bharat gupta
Hi,

I want to know the default cutoff distance g_saltbr tool uses in
identifying the salt-bridges..


Regards
--
Bharat
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] Hbonds per residue

2014-06-12 Thread Natalia Alveal F.
Dear gromacs users,
I have used g_hbond with -hbn option to generate a .ndx file. There is a
method to plot the number of hbonds per residue in gromacs 4.5.5???
This issue is very confuse.
Thanks for any suggestion!!!

-- 
Natalia Alveal Fuentealba
Ingeniera en Bioinformática
Laboratory of Membranes Proteins Structure and Signalling
Departamento de Fisiología
Pontificia Universidad Católica de Chile
Alameda 340, Santiago, Chile
Fono: 56-2-3542877
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Interaction energy between helices

2014-06-12 Thread Justin Lemkul



On 6/11/14, 5:36 PM, Natalia Alveal F. wrote:

Dear Gromacs users,
I want to calculate the energy of interaction between two helices of a
protein, but g_energy dont have the option for -ndx file to select these
specific residues of a helix. How I can do this?



That's what energygrps in the .mdp file are for.  Re-calculate the energies with 
mdrun -rerun using the new .tpr file.


-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] how many groups shall I use to control the temperature?

2014-06-12 Thread 陈功
Dear justin,
 I have doubts that whether I should use two groups to control the 
temperature or just one since there are just 604 hundred Protein beads and 
6+ solvent beads. It is said that two groups are more accurate but when 
when i split them into two groups Protein nonProtein, I found that 
fluctulation of temperature in group Protein are obvious, Im afraid that 
would have an effect on accuracy on my simulation results. So how many groups 
shall I use to control the temperature?
Thanks in advance
regards,
chen
 



-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Why does the list keep unsubscribing me?

2014-06-12 Thread X Rules
This is the 4th time I am thrown out of the list by some random phishing email 
and I am done with this. If I am thrown out this time (I know I will be and 
no-one cares to fix any issues), I am not coming back. 
I just came back the last time to thank you folks for all your comments and 
helps regarding gromacs.have fun !!

 From: xru...@live.com
 To: gmx-us...@gromacs.org
 Date: Mon, 2 Jun 2014 04:16:24 +
 Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
 
 Hi Administrators,   I have been unsubscribed now for the third time (twice 
 in the same month). This is getting increasingly annoying now. Can there be a 
 fix to mailman to look for genuine email, or to send a confirmation link to 
 unsubscribe (the same as for subscription) ?
 Thanks, Abhi 
 
 
  From: xru...@live.com
  To: gmx-us...@gromacs.org
  Date: Thu, 15 May 2014 03:51:37 +
  Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
  
  Thanks for the explanation Rossen.I hope mailman gets some more artificial 
  intelligence in coming future :).cheers,Abhi
  
   Date: Wed, 14 May 2014 14:33:00 +0200
   From: ros...@kth.se
   To: gmx-us...@gromacs.org
   Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
   
   Hi,
   
   Someone else mentioned the same recently.
   
   Automatic removal can happen if someone, e.g. a spammer, sends a mail to 
   gmx-users-unsubscr...@gromacs.org while faking your email address in the 
   From: field.
   
   Mailman is now maintained centrally by the university, so I can't dig in 
   the logs for more information.
   
   Rossen
   
   On 14/05/14 10:29, X Rules wrote:
This is the second time in 2 months, that I am unsubscribed from the 
list and i have to subscribe again. I just get an email that you have 
been unsubscribed from the list?
How can I get unsubscribed without my intervention is beyond 
understanding?
   

   
   -- 
   Gromacs Users mailing list
   
   * Please search the archive at 
   http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before 
   posting!
   
   * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
   
   * For (un)subscribe requests visit
   https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or 
   send a mail to gmx-users-requ...@gromacs.org.

  -- 
  Gromacs Users mailing list
  
  * Please search the archive at 
  http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
  
  * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
  
  * For (un)subscribe requests visit
  https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send 
  a mail to gmx-users-requ...@gromacs.org.
 
 -- 
 Gromacs Users mailing list
 
 * Please search the archive at 
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
 
 * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 
 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
 mail to gmx-users-requ...@gromacs.org.
  
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Why does the list keep unsubscribing me?

2014-06-12 Thread Mark Abraham
Hi,

One simple solution is to subscribe with a different email address.

Mark


On Thu, Jun 12, 2014 at 9:54 AM, X Rules xru...@live.com wrote:

 This is the 4th time I am thrown out of the list by some random phishing
 email and I am done with this. If I am thrown out this time (I know I will
 be and no-one cares to fix any issues), I am not coming back.
 I just came back the last time to thank you folks for all your comments
 and helps regarding gromacs.have fun !!

  From: xru...@live.com
  To: gmx-us...@gromacs.org
  Date: Mon, 2 Jun 2014 04:16:24 +
  Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
 
  Hi Administrators,   I have been unsubscribed now for the third time
 (twice in the same month). This is getting increasingly annoying now. Can
 there be a fix to mailman to look for genuine email, or to send a
 confirmation link to unsubscribe (the same as for subscription) ?
  Thanks, Abhi
 
 
   From: xru...@live.com
   To: gmx-us...@gromacs.org
   Date: Thu, 15 May 2014 03:51:37 +
   Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
  
   Thanks for the explanation Rossen.I hope mailman gets some more
 artificial intelligence in coming future :).cheers,Abhi
  
Date: Wed, 14 May 2014 14:33:00 +0200
From: ros...@kth.se
To: gmx-us...@gromacs.org
Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
   
Hi,
   
Someone else mentioned the same recently.
   
Automatic removal can happen if someone, e.g. a spammer, sends a
 mail to
gmx-users-unsubscr...@gromacs.org while faking your email address
 in the
From: field.
   
Mailman is now maintained centrally by the university, so I can't
 dig in
the logs for more information.
   
Rossen
   
On 14/05/14 10:29, X Rules wrote:
 This is the second time in 2 months, that I am unsubscribed from
 the list and i have to subscribe again. I just get an email that you have
 been unsubscribed from the list?
 How can I get unsubscribed without my intervention is beyond
 understanding?


   
--
Gromacs Users mailing list
   
* Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
 posting!
   
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
   
* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users
 or send a mail to gmx-users-requ...@gromacs.org.
  
   --
   Gromacs Users mailing list
  
   * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
 posting!
  
   * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
  
   * For (un)subscribe requests visit
   https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
 send a mail to gmx-users-requ...@gromacs.org.
 
  --
  Gromacs Users mailing list
 
  * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
 posting!
 
  * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 
  * For (un)subscribe requests visit
  https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
 send a mail to gmx-users-requ...@gromacs.org.

 --
 Gromacs Users mailing list

 * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
 posting!

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

 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
 send a mail to gmx-users-requ...@gromacs.org.

-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Hbonds per residue

2014-06-12 Thread Erik Marklund
Hi Natalia,

It's a feature that has been requested a number of times, but has not been 
implemented yet unfortunately.

Erik

On 12 Jun 2014, at 00:50, Natalia Alveal F. nalv...@bio.puc.cl wrote:

 Dear gromacs users,
 I have used g_hbond with -hbn option to generate a .ndx file. There is a
 method to plot the number of hbonds per residue in gromacs 4.5.5???
 This issue is very confuse.
 Thanks for any suggestion!!!
 
 -- 
 Natalia Alveal Fuentealba
 Ingeniera en Bioinformática
 Laboratory of Membranes Proteins Structure and Signalling
 Departamento de Fisiología
 Pontificia Universidad Católica de Chile
 Alameda 340, Santiago, Chile
 Fono: 56-2-3542877
 -- 
 Gromacs Users mailing list
 
 * Please search the archive at 
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
 
 * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 
 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
 mail to gmx-users-requ...@gromacs.org.

-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] Energy minimisation step in Lipid protein simulation

2014-06-12 Thread Balasubramanian Suriyanarayanan
Dear friends, In protein lipid simulation , when I do energy minimisation
ster I get an error as  Ek! No confout.gro at all!

   Died at inflategro.pl line 104.

What is confout.gro. If am right it is the output file that we get from
the minimisation step.

I used the command  perl inflategro.pl confout.gro 0.95 POPC 0
system_shrink1.gro 5 area_shrink1.dat.


Please help me.

regards,
Suriyanarayanan.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] energy minimisation step in lipid protein simulation

2014-06-12 Thread Balasubramanian Suriyanarayanan
Dear friends, In protein lipid simulation , when I do energy minimisation
ster I get an error as  Ek! No confout.gro at all!

   Died at inflategro.pl line 104.

What is confout.gro. If am right it is the output file that we get from
the minimisation step.

I used the command  perl inflategro.pl confout.gro 0.95 POPC 0
system_shrink1.gro 5 area_shrink1.dat.


Please help me.

regards,
Suriyanarayanan.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] energy minimisation step in protein lipid simulation

2014-06-12 Thread Balasubramanian Suriyanarayanan
Dear friends, In protein lipid simulation , when I do energy minimisation
ster I get an error as  Ek! No confout.gro at all!

   Died at inflategro.pl line 104.

What is confout.gro. If am right it is the output file that we get from
the minimisation step.

I used the command  perl inflategro.pl confout.gro 0.95 POPC 0
system_shrink1.gro 5 area_shrink1.dat.


Please help me.

regards,
Suriyanarayanan.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] pH simulation

2014-06-12 Thread Lovika Moudgil
Hi everyone...

I want to do pH simulation in gromacs . I didn't found any specific link
that I can follow for calculation regarding pH . Can you please guide me
..Any link that can help me...


Thanks in advance

Regards
Lovika
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] pH simulation

2014-06-12 Thread rajat desikan
Hi,
http://www.gromacs.org/Documentation/How-tos/Constant_pH_Simulation

Gromacs users list:
http://comments.gmane.org/gmane.science.biology.gromacs.user/40950
https://www.mail-archive.com/gmx-users@gromacs.org/msg05430.html



On Thu, Jun 12, 2014 at 2:44 PM, Lovika Moudgil lovikamoud...@gmail.com
wrote:

 Hi everyone...

 I want to do pH simulation in gromacs . I didn't found any specific link
 that I can follow for calculation regarding pH . Can you please guide me
 ..Any link that can help me...


 Thanks in advance

 Regards
 Lovika
 --
 Gromacs Users mailing list

 * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
 posting!

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

 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
 send a mail to gmx-users-requ...@gromacs.org.




-- 
Rajat Desikan (Ph.D Scholar)
Prof. K. Ganapathy Ayappa's Lab (no 13),
Dept. of Chemical Engineering,
Indian Institute of Science, Bangalore
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] energy minimisation in lipid protein simulation

2014-06-12 Thread Balasubramanian Suriyanarayanan
Dear friends, In protein lipid simulation , when I do energy minimisation
ster I get an error as  Ek! No confout.gro at all!

   Died at inflategro.pl line 104.

What is confout.gro. If am right it is the output file that we get from
the minimisation step.

I used the command  perl inflategro.pl confout.gro 0.95 POPC 0
system_shrink1.gro 5 area_shrink1.dat.


Please help me.

regards,
Suriyanarayanan.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] (no subject)

2014-06-12 Thread Balasubramanian Suriyanarayanan
Dear friends, In protein lipid simulation , when I do energy minimisation
ster I get an error as  Ek! No confout.gro at all!

   Died at inflategro.pl line 104.

What is confout.gro. If am right it is the output file that we get from
the minimisation step.

I used the command  perl inflategro.pl confout.gro 0.95 POPC 0
system_shrink1.gro 5 area_shrink1.dat.


Please help me.

regards,
Suriyanarayanan.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] (no subject)

2014-06-12 Thread Justin Lemkul



On 6/12/14, 5:31 AM, Balasubramanian Suriyanarayanan wrote:

Why is my questions not appearing in the forum.

Do I not follow the rules.'
Please clarify.



Your messages are being posted just fine.  Please stop spamming the list.

-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Energy minimisation step in Lipid protein simulation

2014-06-12 Thread Justin Lemkul



On 6/12/14, 4:54 AM, Balasubramanian Suriyanarayanan wrote:

Dear friends, In protein lipid simulation , when I do energy minimisation
ster I get an error as  Ek! No confout.gro at all!

Died at inflategro.pl line 104.

What is confout.gro. If am right it is the output file that we get from
the minimisation step.

I used the command  perl inflategro.pl confout.gro 0.95 POPC 0
system_shrink1.gro 5 area_shrink1.dat.



You're specifying confout.gro as the source of coordinates.  Apparently that 
file doesn't exist in the working directory, so InflateGRO is giving you an error.


-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] how many groups shall I use to control the temperature?

2014-06-12 Thread Justin Lemkul



On 6/12/14, 2:37 AM, 陈功 wrote:

Dear justin,
  I have doubts that whether I should use two groups to control the temperature or just one 
since there are just 604 hundred Protein beads and 6+ solvent beads. It is said that two groups 
are more accurate but when when i split them into two groups Protein nonProtein, I 
found that fluctulation of temperature in group Protein are obvious, Im afraid that 
would have an effect on accuracy on my simulation results. So how many groups shall I use to 
control the temperature?


What protocol is used in published studies similar to what you're doing?  Small 
groups will have larger fluctuations in temperature, that's a given.  I don't do 
CG simulations, so I don't know how such groups are commonly treated.


-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] -inf potential with TPIC

2014-06-12 Thread Justin Lemkul



On 6/11/14, 8:55 PM, Rafael I. Silverman y de la Vega wrote:

I mistyped amber, when I should have typed charmm, amber is totally
uninvolved in  my calculations. I used the script from
http://www.gromacs.org/Downloads/User_contributions/Other_software titled
charmm2gromacs. The output looked fine after running it.

And, I doubt its the cofactor,  just ran the insertion with double
precision mdrun, it gave reasonable looking results. There seems to be some
sort of issue with unreasonable interaction energies not being discarded
properly.



It's entirely possible that there's something wrong with the tpic code; you can 
certainly file a bug report on Redmine with input files to reproduce it.


-Justin



On Tue, Jun 10, 2014 at 3:41 PM, Justin Lemkul jalem...@vt.edu wrote:




On 6/10/14, 6:14 PM, Rafael I. Silverman y de la Vega wrote:


I meant that I only changed to cofactor to gromacs chamm27 format, not
amber format. I used charmm27 for the apoprotein as well as the cofactor.



What script did you use?  I don't understand why Amber99SB is involved
here, so I don't know exactly what you were doing to create the topology.
  It may not end up being useful information, but it's important to
establish the validity of the underlying physical model.


-Justin



On Tue, Jun 10, 2014 at 3:09 PM, Justin Lemkul jalem...@vt.edu wrote:




On 6/10/14, 6:04 PM, Rafael I. Silverman y de la Vega wrote:

  you think it was the cofactor model?



  If you're using a CHARMM-parametrized cofactor with an Amber99SB

protein,
I'd say the simulations are wholly unreliable.  It's certainly the #1
likely cause for nonsensical results.  If you want to get to the root of
the issue, you need to start with a solid foundation.

Unless that's not what your last message meant?  When you say only
changed the cofactor, that implied to me that you didn't consistently
convert everything in the system.


-Justin


  On Tue, Jun 10, 2014 at 2:23 PM, Justin Lemkul jalem...@vt.edu wrote:





On 6/10/14, 5:01 PM, Rafael I. Silverman y de la Vega wrote:

   wait, I only changed  the cofactor to gromacs charmm27.ff format




   Pretty amazing that the simulation even ran with that mash-up of


parameters :)

I'd say that's the most likely cause of your problems; nonsensical
physical model leads to nonsensical output.

-Justin



   On Tue, Jun 10, 2014 at 1:58 PM, Rafael I. Silverman y de la Vega 


rsilv...@ucsc.edu wrote:

I suppose something might be wrong with the simulations. I am using

  gromacs 4.6.5, I am simulating a protein with a cofactor, someone

else
parametrized it, I changed from charmm to gromacs amber99sb.ff format
with
a script on the gromacs website.  I did a typical em/nvt/npt
equilibration
followed by a 5 ns production run. I picked 250 frames out of the
trajectory, and altered them to the TPIC input format with a script I
wrote, to insert the particle half way between two atoms that I
chose.
I
am
inserting a water molecule. When I run the actual insertion some of
the
frames have a reasonable looking mu, and some have -inf mu. I think
it
may
have to do with cavity size, some frames have too small of a cavity,
so
the
interaction energy of the inserted water diverges...


On Mon, Jun 9, 2014 at 6:40 PM, Justin Lemkul jalem...@vt.edu
wrote:



  On 6/9/14, 6:41 PM, Rafael I. Silverman y de la Vega wrote:


Hi all, I am having a problem with -inf mu coming out of TPIC

  calculations.

Anyone have some suggestions on how to avoid this? Tabulated
potentials?
Or
can I modify the source code to reject the absurd values that cause
this? I
read a few old discussions on this, I didnt see much on how to
handle
this.


I'd be much more concerned that something is going wrong in the

  simulations.  I have no experience with TPIC, but if you post more

detailed
information about what you're doing and which Gromacs version you're
using,
you're more likely to get some suggestions.

-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at http://www.gromacs.org/
Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users
or
send a mail to gmx-users-requ...@gromacs.org.





   --


==


Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, 

Re: [gmx-users] (no subject)

2014-06-12 Thread Vedat Durmaz

i've recieved your email concernig energy minimisation step in protein
lipid simulation already 3 times ..


maybe you've deactivated the delivery of mails? if so, do the following:

go to https://mailman-1.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users

and log in where it says To unsubscribe from gromacs.org_gmx-users, get
a password reminder, or change your subscription options enter your
subscription email address:

there you can change your personal settings.

vedat



Am 12.06.2014 11:31, schrieb Balasubramanian Suriyanarayanan:
 Why is my questions not appearing in the forum.

 Do I not follow the rules.'
 Please clarify.


 regards

 Suriyanarayanan

-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] (no subject)

2014-06-12 Thread Vedat Durmaz
then, i don't know what to do. may be your email provider sorts them out
or something.

anyway, did you get justin's answer he sent recently? see below ..

vedat



On 6/12/14, 4:54 AM, Balasubramanian Suriyanarayanan wrote:
 Dear friends, In protein lipid simulation , when I do energy minimisation
 ster I get an error as  Ek! No confout.gro at all!

 Died at inflategro.pl line 104.

 What is confout.gro. If am right it is the output file that we get from
 the minimisation step.

 I used the command  perl inflategro.pl confout.gro 0.95 POPC 0
 system_shrink1.gro 5 area_shrink1.dat.


You're specifying confout.gro as the source of coordinates.  Apparently
that file doesn't exist in the working directory, so InflateGRO is
giving you an error.

-Justin



Am 12.06.2014 11:46, schrieb Balasubramanian Suriyanarayanan:
 but I get other mails.

 even the last one which i have sent with out title has appeared.

 regards
 Suriyanarayanan


 On Thu, Jun 12, 2014 at 3:14 PM, Vedat Durmaz dur...@zib.de
 mailto:dur...@zib.de wrote:


 i've recieved your email concernig energy minimisation step in
 protein
 lipid simulation already 3 times ..


 maybe you've deactivated the delivery of mails? if so, do the
 following:

 go to
 https://mailman-1.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users

 and log in where it says To unsubscribe from
 gromacs.org_gmx-users, get
 a password reminder, or change your subscription options enter your
 subscription email address:

 there you can change your personal settings.

 vedat



 Am 12.06.2014 11:31, schrieb Balasubramanian Suriyanarayanan:
  Why is my questions not appearing in the forum.
 
  Do I not follow the rules.'
  Please clarify.
 
 
  regards
 
  Suriyanarayanan



-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] (no subject)

2014-06-12 Thread Vedat Durmaz

i've recieved your email concernig energy minimisation step in protein
lipid simulation already 3 times ..


maybe you've deactivated the delivery of mails? if so, do the following:

go to https://mailman-1.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users

and log in where it says To unsubscribe from gromacs.org_gmx-users, get
a password reminder, or change your subscription options enter your
subscription email address:

there you can change your personal settings.

vedat



Am 12.06.2014 11:31, schrieb Balasubramanian Suriyanarayanan:
 Why is my questions not appearing in the forum.

 Do I not follow the rules.'
 Please clarify.


 regards

 Suriyanarayanan

-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Tc grps, one group or two?

2014-06-12 Thread Mark Abraham
On Thu, Jun 12, 2014 at 4:11 AM, 陈功 gchen...@gmail.com wrote:

 Dear gromacs users,
  I have doubts that whether I should use two groups to control the
 temperature or just one since there are just 604 hundred Protein beads and
 6+ solvent beads. It is said that two groups are more accurate but when
 when i split them into two groups Protein nonProtein, I found that
 fluctulation of temperature in group Protein are obvious, Im afraid that
 would have an effect on accuracy on my simulation results.
   Can anyone help me ?  thanks in advance


This is a complex issue. Even with a proper thermostat (i.e not Berendsen),
making a sound observation of temperature is tricky. See
http://www.deshawresearch.com/publications/Equipartition%20and%20the%20Calculation%20of%20Temperature%20in%20Biomolecular%20Simulations.pdf
Measuring and managing temperature in the presence of other possible
problems (invalid or non-transferable force-field parameters,
non-conservative integration schemes) is even harder. Whether any of this
is relevant to your intended observables is usually unknown, because
typically no previous work was based on simulations sufficiently converged
to be definitive.

With orders of magnitude more solvent than solute, I would generally
suggest a single temperature-coupling group, but there may be issues with
how your coarse-grained force field was parameterized of which I am unaware.

Mark

regards,
 chen




 --
 Gromacs Users mailing list

 * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
 posting!

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

 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
 send a mail to gmx-users-requ...@gromacs.org.

-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Energy minimisation step in Lipid protein simulation

2014-06-12 Thread Balasubramanian Suriyanarayanan
But is not confout.gro the output of the compression step.

the command is
perl inflategro.pl system.gro 4 POPC 0 system_inflated.gro 5 area.dat

or else what is the output of above command.

regards

suriyanarayanan





On Thu, Jun 12, 2014 at 3:13 PM, Justin Lemkul jalem...@vt.edu wrote:



 On 6/12/14, 4:54 AM, Balasubramanian Suriyanarayanan wrote:

 Dear friends, In protein lipid simulation , when I do energy minimisation
 ster I get an error as  Ek! No confout.gro at all!

 Died at inflategro.pl line 104.

 What is confout.gro. If am right it is the output file that we get from
 the minimisation step.

 I used the command  perl inflategro.pl confout.gro 0.95 POPC 0
 system_shrink1.gro 5 area_shrink1.dat.


 You're specifying confout.gro as the source of coordinates.  Apparently
 that file doesn't exist in the working directory, so InflateGRO is giving
 you an error.

 -Justin

 --
 ==

 Justin A. Lemkul, Ph.D.
 Ruth L. Kirschstein NRSA Postdoctoral Fellow

 Department of Pharmaceutical Sciences
 School of Pharmacy
 Health Sciences Facility II, Room 601
 University of Maryland, Baltimore
 20 Penn St.
 Baltimore, MD 21201

 jalem...@outerbanks.umaryland.edu | (410) 706-7441
 http://mackerell.umaryland.edu/~jalemkul

 ==
 --
 Gromacs Users mailing list

 * Please search the archive at http://www.gromacs.org/
 Support/Mailing_Lists/GMX-Users_List before posting!

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

 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
 send a mail to gmx-users-requ...@gromacs.org.

-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Why does the list keep unsubscribing me?

2014-06-12 Thread X Rules
Thanks for the solution Mark, but from what I feel the right way to do is to 
fix mailman (add another option of clicking a link sent to your email to 
confirm unsubscribe, the same way subscription works). I might not be the only 
one, or may be I am the only one complaining :). 

Lets see how long this lasts,Abhi

 Date: Thu, 12 Jun 2014 10:22:48 +0200
 From: mark.j.abra...@gmail.com
 To: gmx-us...@gromacs.org
 Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
 
 Hi,
 
 One simple solution is to subscribe with a different email address.
 
 Mark
 
 
 On Thu, Jun 12, 2014 at 9:54 AM, X Rules xru...@live.com wrote:
 
  This is the 4th time I am thrown out of the list by some random phishing
  email and I am done with this. If I am thrown out this time (I know I will
  be and no-one cares to fix any issues), I am not coming back.
  I just came back the last time to thank you folks for all your comments
  and helps regarding gromacs.have fun !!
 
   From: xru...@live.com
   To: gmx-us...@gromacs.org
   Date: Mon, 2 Jun 2014 04:16:24 +
   Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
  
   Hi Administrators,   I have been unsubscribed now for the third time
  (twice in the same month). This is getting increasingly annoying now. Can
  there be a fix to mailman to look for genuine email, or to send a
  confirmation link to unsubscribe (the same as for subscription) ?
   Thanks, Abhi
  
  
From: xru...@live.com
To: gmx-us...@gromacs.org
Date: Thu, 15 May 2014 03:51:37 +
Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
   
Thanks for the explanation Rossen.I hope mailman gets some more
  artificial intelligence in coming future :).cheers,Abhi
   
 Date: Wed, 14 May 2014 14:33:00 +0200
 From: ros...@kth.se
 To: gmx-us...@gromacs.org
 Subject: Re: [gmx-users] Why does the list keep unsubscribing me?

 Hi,

 Someone else mentioned the same recently.

 Automatic removal can happen if someone, e.g. a spammer, sends a
  mail to
 gmx-users-unsubscr...@gromacs.org while faking your email address
  in the
 From: field.

 Mailman is now maintained centrally by the university, so I can't
  dig in
 the logs for more information.

 Rossen

 On 14/05/14 10:29, X Rules wrote:
  This is the second time in 2 months, that I am unsubscribed from
  the list and i have to subscribe again. I just get an email that you have
  been unsubscribed from the list?
  How can I get unsubscribed without my intervention is beyond
  understanding?
 
 

 --
 Gromacs Users mailing list

 * Please search the archive at
  http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
  posting!

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

 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users
  or send a mail to gmx-users-requ...@gromacs.org.
   
--
Gromacs Users mailing list
   
* Please search the archive at
  http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
  posting!
   
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
   
* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
  send a mail to gmx-users-requ...@gromacs.org.
  
   --
   Gromacs Users mailing list
  
   * Please search the archive at
  http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
  posting!
  
   * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
  
   * For (un)subscribe requests visit
   https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
  send a mail to gmx-users-requ...@gromacs.org.
 
  --
  Gromacs Users mailing list
 
  * Please search the archive at
  http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
  posting!
 
  * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 
  * For (un)subscribe requests visit
  https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
  send a mail to gmx-users-requ...@gromacs.org.
 
 -- 
 Gromacs Users mailing list
 
 * Please search the archive at 
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
 
 * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 
 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
 mail to gmx-users-requ...@gromacs.org.
  
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Energy minimisation step in Lipid protein simulation

2014-06-12 Thread Justin Lemkul



On 6/12/14, 7:14 AM, Balasubramanian Suriyanarayanan wrote:

But is not confout.gro the output of the compression step.

the command is
perl inflategro.pl system.gro 4 POPC 0 system_inflated.gro 5 area.dat

or else what is the output of above command.



Read the InflateGRO documentation (top of the script).

-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Why does the list keep unsubscribing me?

2014-06-12 Thread Mark Abraham
Hi,

I've no idea if mailman even has that option, or the KTH IT admins would be
willing to set it up, but that's not worth exploring until automatic
unsubscription is known to be a problem for more than one email address.

Mark


On Thu, Jun 12, 2014 at 1:34 PM, X Rules xru...@live.com wrote:

 Thanks for the solution Mark, but from what I feel the right way to do is
 to fix mailman (add another option of clicking a link sent to your email to
 confirm unsubscribe, the same way subscription works). I might not be the
 only one, or may be I am the only one complaining :).

 Lets see how long this lasts,Abhi

  Date: Thu, 12 Jun 2014 10:22:48 +0200
  From: mark.j.abra...@gmail.com
  To: gmx-us...@gromacs.org
  Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
 
  Hi,
 
  One simple solution is to subscribe with a different email address.
 
  Mark
 
 
  On Thu, Jun 12, 2014 at 9:54 AM, X Rules xru...@live.com wrote:
 
   This is the 4th time I am thrown out of the list by some random
 phishing
   email and I am done with this. If I am thrown out this time (I know I
 will
   be and no-one cares to fix any issues), I am not coming back.
   I just came back the last time to thank you folks for all your comments
   and helps regarding gromacs.have fun !!
  
From: xru...@live.com
To: gmx-us...@gromacs.org
Date: Mon, 2 Jun 2014 04:16:24 +
Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
   
Hi Administrators,   I have been unsubscribed now for the third time
   (twice in the same month). This is getting increasingly annoying now.
 Can
   there be a fix to mailman to look for genuine email, or to send a
   confirmation link to unsubscribe (the same as for subscription) ?
Thanks, Abhi
   
   
 From: xru...@live.com
 To: gmx-us...@gromacs.org
 Date: Thu, 15 May 2014 03:51:37 +
 Subject: Re: [gmx-users] Why does the list keep unsubscribing me?

 Thanks for the explanation Rossen.I hope mailman gets some more
   artificial intelligence in coming future :).cheers,Abhi

  Date: Wed, 14 May 2014 14:33:00 +0200
  From: ros...@kth.se
  To: gmx-us...@gromacs.org
  Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
 
  Hi,
 
  Someone else mentioned the same recently.
 
  Automatic removal can happen if someone, e.g. a spammer, sends a
   mail to
  gmx-users-unsubscr...@gromacs.org while faking your email
 address
   in the
  From: field.
 
  Mailman is now maintained centrally by the university, so I can't
   dig in
  the logs for more information.
 
  Rossen
 
  On 14/05/14 10:29, X Rules wrote:
   This is the second time in 2 months, that I am unsubscribed
 from
   the list and i have to subscribe again. I just get an email that you
 have
   been unsubscribed from the list?
   How can I get unsubscribed without my intervention is beyond
   understanding?
  
  
 
  --
  Gromacs Users mailing list
 
  * Please search the archive at
   http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
   posting!
 
  * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 
  * For (un)subscribe requests visit
 
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users
   or send a mail to gmx-users-requ...@gromacs.org.

 --
 Gromacs Users mailing list

 * Please search the archive at
   http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
   posting!

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

 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users
 or
   send a mail to gmx-users-requ...@gromacs.org.
   
--
Gromacs Users mailing list
   
* Please search the archive at
   http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
   posting!
   
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
   
* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users
 or
   send a mail to gmx-users-requ...@gromacs.org.
  
   --
   Gromacs Users mailing list
  
   * Please search the archive at
   http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
   posting!
  
   * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
  
   * For (un)subscribe requests visit
   https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
   send a mail to gmx-users-requ...@gromacs.org.
  
  --
  Gromacs Users mailing list
 
  * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
 posting!
 
  * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 
  * For (un)subscribe requests visit
  https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
 send a mail to gmx-users-requ...@gromacs.org.

 

Re: [gmx-users] how many groups shall I use to control the temperature?

2014-06-12 Thread 陈功
Thanks a lot, Justin. I do CG simulation using Martini force field. So maybe I 
can find some useful infomation on the forum of martini.
Thanks :)




 
From: Justin Lemkul
Date: 2014-06-12 17:44
To: gmx-users
Subject: Re: [gmx-users] how many groups shall I use to control the temperature?
 
 
On 6/12/14, 2:37 AM, 陈功 wrote:
 Dear justin,
   I have doubts that whether I should use two groups to control the 
 temperature or just one since there are just 604 hundred Protein beads and 
 6+ solvent beads. It is said that two groups are more accurate but when 
 when i split them into two groups Protein nonProtein, I found that 
 fluctulation of temperature in group Protein are obvious, Im afraid that 
 would have an effect on accuracy on my simulation results. So how many groups 
 shall I use to control the temperature?
 
What protocol is used in published studies similar to what you're doing?  Small 
groups will have larger fluctuations in temperature, that's a given.  I don't 
do 
CG simulations, so I don't know how such groups are commonly treated.
 
-Justin
 
-- 
==
 
Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow
 
Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201
 
jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul
 
==
-- 
Gromacs Users mailing list
 
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
 
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 
* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] Acpyte with opls-aa

2014-06-12 Thread Saman Shahriyari
Dear users



i wanted to know if Acpyte generated topologies for a ligand are compatible 
with OPLS-AA.ff? i am searching for a reliable topology builder for OPLS-aa 
like ATP for gromos53a6.ff so that i can avoid manual approached like using 
gaussian .
 regards,
Saman---BeginMessage---
Dear users
i wanted to know if Acpyte generated topologies for a ligand are compatible 
with OPLS-AA.ff? i am searching for a reliable topology builder for OPLS-aa 
like ATP for gromos53a6.ff so that i can avoid manual approached like using 
gaussian .
 regards,
Saman---End Message---
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] g_clustsize

2014-06-12 Thread David van der Spoel

On 2014-06-12 15:33, Raisa Kociurzynski wrote:

Dear Prof. Van der Spoel,

I want to calculate the clustersize of GM1 clusters in a DPPC membrane. I used 
g_clustersize but when I use the option -mol the index file is ignored and the 
whole system is used for the calculation. Is there any way to specify only GM1 
molecules for the cluster calculation? I also tried to modify my xtc and tpr 
files, so that my system consists only of GM1 molecules, but then I only get 
one big cluster with g_clustsize. I could not find any solution for this 
problem in the gromacs forum.

Best regards,

Raisa Kociurzynski
If you have a membrane there is only one cluster unless the membrane is 
fragmented into micelles.



--
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.
sp...@xray.bmc.uu.sehttp://folding.bmc.uu.se

--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] -inf potential with TPIC

2014-06-12 Thread João M . Damas
Hello Rafael,

High energies (energy differences) are not discarded. When summing them
through the exponential of their negative, they contribute very little to
the final sum (lower energies dominate the sum). Hence, they are
discarded.

The -inf value on the mu standard error output is normal because, unlike
the -tpi output, there is no ternary operator checking if the sum of
exponentials is equal to zero (when it is, the log of zero gives you the
inf value).

The change to the double precision not giving you the -inf also makes
sense, because the sum of exponentials, which on single-precision was zero,
can now be non-zero due to the higher number of decimal places on
double-precision. I am guessing with double-precision you are getting very
high positive numbers (energies) for the mu, around 1000 maybe?

My interpretation is that it is hard to find a favorable place to insert
the water molecule in that cavity, at least with a few frames. I am
guessing that with more and more frames, that may change, which is natural.

I hope this has helped.

Cheers,
João



On Thu, Jun 12, 2014 at 10:45 AM, Justin Lemkul jalem...@vt.edu wrote:



 On 6/11/14, 8:55 PM, Rafael I. Silverman y de la Vega wrote:

 I mistyped amber, when I should have typed charmm, amber is totally
 uninvolved in  my calculations. I used the script from
 http://www.gromacs.org/Downloads/User_contributions/Other_software titled
 charmm2gromacs. The output looked fine after running it.

 And, I doubt its the cofactor,  just ran the insertion with double
 precision mdrun, it gave reasonable looking results. There seems to be
 some
 sort of issue with unreasonable interaction energies not being discarded
 properly.


 It's entirely possible that there's something wrong with the tpic code;
 you can certainly file a bug report on Redmine with input files to
 reproduce it.


 -Justin


 On Tue, Jun 10, 2014 at 3:41 PM, Justin Lemkul jalem...@vt.edu wrote:



 On 6/10/14, 6:14 PM, Rafael I. Silverman y de la Vega wrote:

  I meant that I only changed to cofactor to gromacs chamm27 format, not
 amber format. I used charmm27 for the apoprotein as well as the
 cofactor.


  What script did you use?  I don't understand why Amber99SB is involved
 here, so I don't know exactly what you were doing to create the topology.
   It may not end up being useful information, but it's important to
 establish the validity of the underlying physical model.


 -Justin


  On Tue, Jun 10, 2014 at 3:09 PM, Justin Lemkul jalem...@vt.edu wrote:



 On 6/10/14, 6:04 PM, Rafael I. Silverman y de la Vega wrote:

   you think it was the cofactor model?



   If you're using a CHARMM-parametrized cofactor with an Amber99SB

 protein,
 I'd say the simulations are wholly unreliable.  It's certainly the #1
 likely cause for nonsensical results.  If you want to get to the root
 of
 the issue, you need to start with a solid foundation.

 Unless that's not what your last message meant?  When you say only
 changed the cofactor, that implied to me that you didn't consistently
 convert everything in the system.


 -Justin


   On Tue, Jun 10, 2014 at 2:23 PM, Justin Lemkul jalem...@vt.edu
 wrote:




  On 6/10/14, 5:01 PM, Rafael I. Silverman y de la Vega wrote:

wait, I only changed  the cofactor to gromacs charmm27.ff format



Pretty amazing that the simulation even ran with that mash-up of

  parameters :)

 I'd say that's the most likely cause of your problems; nonsensical
 physical model leads to nonsensical output.

 -Justin



On Tue, Jun 10, 2014 at 1:58 PM, Rafael I. Silverman y de la Vega
 

  rsilv...@ucsc.edu wrote:

 I suppose something might be wrong with the simulations. I am
 using

   gromacs 4.6.5, I am simulating a protein with a cofactor, someone

 else
 parametrized it, I changed from charmm to gromacs amber99sb.ff
 format
 with
 a script on the gromacs website.  I did a typical em/nvt/npt
 equilibration
 followed by a 5 ns production run. I picked 250 frames out of the
 trajectory, and altered them to the TPIC input format with a
 script I
 wrote, to insert the particle half way between two atoms that I
 chose.
 I
 am
 inserting a water molecule. When I run the actual insertion some of
 the
 frames have a reasonable looking mu, and some have -inf mu. I think
 it
 may
 have to do with cavity size, some frames have too small of a
 cavity,
 so
 the
 interaction energy of the inserted water diverges...


 On Mon, Jun 9, 2014 at 6:40 PM, Justin Lemkul jalem...@vt.edu
 wrote:



   On 6/9/14, 6:41 PM, Rafael I. Silverman y de la Vega wrote:


 Hi all, I am having a problem with -inf mu coming out of TPIC

   calculations.

 Anyone have some suggestions on how to avoid this? Tabulated
 potentials?
 Or
 can I modify the source code to reject the absurd values that
 cause
 this? I
 read a few old discussions on this, I didnt see much on how to
 handle
 this.


 I'd be much more concerned that something is going wrong in
 the

   simulations.  I 

[gmx-users] question about topolbuild

2014-06-12 Thread Hui Wen Ng


 Dear gmxusers,
 
 I have a question about using topolbuild to generate
 topology for my ligand (intend to use OPLS AA forcefield).
 However, I am facing a problem in using the command, which
 according to the README text, requires the basic input as
 below
 
  topolbuild -n modelname -dir path_to_interp. -ff
 forcefield
 
 I have a ligand.mol2 and I typed in the command as below
 
 topolbuild1_3/src/topolbuild -n ligand -ff oplsaa  -dir
 topolbuild1_3/dat/gromacs/
 
 I got this error: Segmentation fault
 
 Could anyone helpme with this please? Thanks a lot.
 
 Best,
 HW
 
 
 
 
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] how many groups shall I use to control the temperature?

2014-06-12 Thread Dr. Vitaly Chaban
There is no rule regarding how many thermostats per system one needs.

Often people choose separate thermostats for parts of system, which
exhibit different mobility. I do not see any connection to functional
forms of the interactions and, therefore, to force fields.


Dr. Vitaly V. Chaban


On Thu, Jun 12, 2014 at 2:33 PM, 陈功 gchen...@gmail.com wrote:
 Thanks a lot, Justin. I do CG simulation using Martini force field. So maybe 
 I can find some useful infomation on the forum of martini.
 Thanks :)





 From: Justin Lemkul
 Date: 2014-06-12 17:44
 To: gmx-users
 Subject: Re: [gmx-users] how many groups shall I use to control the 
 temperature?


 On 6/12/14, 2:37 AM, 陈功 wrote:
 Dear justin,
   I have doubts that whether I should use two groups to control the 
 temperature or just one since there are just 604 hundred Protein beads and 
 6+ solvent beads. It is said that two groups are more accurate but when 
 when i split them into two groups Protein nonProtein, I found that 
 fluctulation of temperature in group Protein are obvious, Im afraid that 
 would have an effect on accuracy on my simulation results. So how many 
 groups shall I use to control the temperature?

 What protocol is used in published studies similar to what you're doing?  
 Small
 groups will have larger fluctuations in temperature, that's a given.  I don't 
 do
 CG simulations, so I don't know how such groups are commonly treated.

 -Justin

 --
 ==

 Justin A. Lemkul, Ph.D.
 Ruth L. Kirschstein NRSA Postdoctoral Fellow

 Department of Pharmaceutical Sciences
 School of Pharmacy
 Health Sciences Facility II, Room 601
 University of Maryland, Baltimore
 20 Penn St.
 Baltimore, MD 21201

 jalem...@outerbanks.umaryland.edu | (410) 706-7441
 http://mackerell.umaryland.edu/~jalemkul

 ==
 --
 Gromacs Users mailing list

 * Please search the archive at 
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
 mail to gmx-users-requ...@gromacs.org.
 --
 Gromacs Users mailing list

 * Please search the archive at 
 http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

 * For (un)subscribe requests visit
 https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
 mail to gmx-users-requ...@gromacs.org.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Interaction energy between helices

2014-06-12 Thread Justin Lemkul



On 6/12/14, 2:02 PM, Natalia Alveal F. wrote:

Thanks Justin, but for how long do this new simulation?


It's not a new simulation.  You're re-processing your existing trajectory. 
Usually the energy evaluations only take a few minutes.  The practical impact of 
energygrps during a run is usually fairly insignificant, but that's also a 
function of how many groups there are and how often you're writing to disk.


-Justin


Natalia


2014-06-11 18:14 GMT-04:00 Justin Lemkul jalem...@vt.edu:




On 6/11/14, 5:36 PM, Natalia Alveal F. wrote:


Dear Gromacs users,
I want to calculate the energy of interaction between two helices of a
protein, but g_energy dont have the option for -ndx file to select these
specific residues of a helix. How I can do this?



That's what energygrps in the .mdp file are for.  Re-calculate the
energies with mdrun -rerun using the new .tpr file.

-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at http://www.gromacs.org/
Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
send a mail to gmx-users-requ...@gromacs.org.







--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] Adding residue at two end of DNA

2014-06-12 Thread Mehdi Bagherpour
Dear Gromacs User,

Do you know how can I add some residue to ends of DNA?

I am looking forward hearing of you,=
Mahdi Bagherpoor
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] hardware problem of GPU?

2014-06-12 Thread Albert

thanks a lot.

Finally both GPU are working.

Albert

On 06/11/2014 01:21 PM, Mark Abraham wrote:

This is not a hardware issue. See
http://www.gromacs.org/Documentation/Acceleration_and_parallelization#Heterogenous_parallelization.3a_using_GPUs

Mark


--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Using specbond.dat in pdb2gmx

2014-06-12 Thread Matthew Stancea
On 6/11/14, 2:33 PM, Matthew Stancea wrote:
 ? Hello,

 I have been having a bit of issues generating an accurate gromacs topology
 file (topol.top) utilizing a pdb of a cyclic peptide in pdb2gmx. I have
 been able to generate topologies that are almost identical to the original
 pdb using the option -ignh at the end of my command, but doing so deletes
 the bond between the first nitrogen and the last carbon (which should be
 connected for this peptide to be cyclic) and adds two additional hydrogens
 and a positive charge to the first nitrogen and an additional oxygen and a
 negative charge to the last carbon, rendering this peptide as non-cyclic.
 After searching around for quite a while, I found out that many others on
 this mailing list were having the same issues as myself, and some replies to
 their messages including the usage of a file known as specbond.dat which
 may be helpful in retaining that bond.


The -ignh flag is not relevant to those observations.  

In my working directory, I have a file named specbond.dat and it contains the 
following information:

3
CYS N   1   VAL C   1   0.25CYS2VAL2
CYS SG  1   CYS SG  1   0.25CYS2CYS2
CYM SG  1   CYM SG  1   0.25CYS2CYS2


When I input the command pdb2gmx_mpi -ff amber99sb -water tip3p -f 1NB1.pdb -o 
conf.gro -p topol.top -i kalata.itp while specbond.dat is in the working 
directory (Here is a link for the structure of peptide 1NB1 for reference if 
needed: http://www.ncbi.nlm.nih.gov/protein/1NB1_A ), I get the following 
message: 

---
Program pdb2gmx_mpi, VERSION 4.6.2
Source code file: /home/msaum/apps/gromacs-4.6.2/src/kernel/pdb2gmx.c, line: 727

Fatal error:
Atom HB3 in residue CYS 1 was not found in rtp entry NCYS with 13 atoms
while sorting atoms.

For a hydrogen, this can be a different protonation state, or it
might have had a different number in the PDB file and was rebuilt
(it might for instance have been H3, and we only expected H1  H2).
Note that hydrogens might have been added to the entry for the N-terminus.
Remove this hydrogen or choose a different protonation state to solve it.
Option -ignh will ignore all hydrogens in the input.
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors
---


(HB3 would be the hydrogen of a cysteine residue that has not formed a 
disulfide bond with another amino acid residue, however my pdb does not contain 
an HB3).

I then tried adding -ignh as per the suggestion in the twelfth line in that 
message. The following is a part of the message I received:

Checking for duplicate atoms
Generating any missing hydrogen atoms and/or adding termini.
Now there are 29 residues with 385 atoms
Making bonds...
Number of bonds was 390, now 389
Generating angles, dihedrals and pairs...
Before cleaning: 1010 pairs
Before cleaning: 1030 dihedrals
Keeping all generated dihedrals
Making cmap torsions...There are 1030 dihedrals,   73 impropers,  699 angles
  1007 pairs,  389 bonds and 0 virtual sites
Total mass 2916.376 a.m.u.
Total charge -0.000 e
Writing topology

Writing coordinate file...
- PLEASE NOTE 
You have successfully generated a topology from: 1NB1.pdb.
The Amber99sb force field and the tip3p water model are used.
- ETON ESAELP 


On the fifth line of that excerpt, I read Number of bonds was 390, now 389 
and questioned what that could mean. Because of that, I wanted to verify that 
the topology that was successfully generated contains all the necessary 
bonds, so I looked through the topol.top file and saw 3 total hydrogens on the 
terminal nitrogen on the beginning cysteine and two total oxygens on the 
terminal carbon on the ending valine. Also, I did not see a bond between that 
nitrogen and that carbon (atoms 1 and 383); these data lead me to believe that 
this peptide is no longer cyclic.

Construction of termini and assignment of ionization 
state is done with -ter.  For your case, you likely
need to be using -ter and selecting None.


When I input the command pdb2gmx_mpi -ff amber99sb -water tip3p -f 1NB1.pdb -o 
conf.gro -p topol.top -i kalata.itp -ter, I received the same message as 
before without any kind of interactivity for selecting options. Because of 
this, I tried the command pdb2gmx_mpi -ff amber99sb -water tip3p -f 1NB1.pdb 
-o conf.gro -p topol.top -i kalata.itp -inter -ter, which asked about the 
protonation of residue 24 (arginine), and when I answer none, I get this 
message:

Processing chain 1 'A' (376 atoms, 29 residues)
Which ARGININE type do you want for residue 24
0. Not protonated (charge 0) (-)
1. Protonated (charge +1) (ARG)

Type a number:none

---
Program 

Re: [gmx-users] Using specbond.dat in pdb2gmx

2014-06-12 Thread Justin Lemkul



On 6/12/14, 4:53 PM, Matthew Stancea wrote:

On 6/11/14, 2:33 PM, Matthew Stancea wrote:

? Hello,

I have been having a bit of issues generating an accurate gromacs topology
file (topol.top) utilizing a pdb of a cyclic peptide in pdb2gmx. I have
been able to generate topologies that are almost identical to the original
pdb using the option -ignh at the end of my command, but doing so deletes
the bond between the first nitrogen and the last carbon (which should be
connected for this peptide to be cyclic) and adds two additional hydrogens
and a positive charge to the first nitrogen and an additional oxygen and a
negative charge to the last carbon, rendering this peptide as non-cyclic.
After searching around for quite a while, I found out that many others on
this mailing list were having the same issues as myself, and some replies to
their messages including the usage of a file known as specbond.dat which
may be helpful in retaining that bond.




The -ignh flag is not relevant to those observations.


In my working directory, I have a file named specbond.dat and it contains the 
following information:

3
CYS N   1   VAL C   1   0.25CYS2VAL2
CYS SG  1   CYS SG  1   0.25CYS2CYS2
CYM SG  1   CYM SG  1   0.25CYS2CYS2


When I input the command pdb2gmx_mpi -ff amber99sb -water tip3p -f 1NB1.pdb -o 
conf.gro -p topol.top -i kalata.itp while specbond.dat is in the working directory 
(Here is a link for the structure of peptide 1NB1 for reference if needed: 
http://www.ncbi.nlm.nih.gov/protein/1NB1_A ), I get the following message:

---
Program pdb2gmx_mpi, VERSION 4.6.2
Source code file: /home/msaum/apps/gromacs-4.6.2/src/kernel/pdb2gmx.c, line: 727

Fatal error:
Atom HB3 in residue CYS 1 was not found in rtp entry NCYS with 13 atoms
while sorting atoms.

For a hydrogen, this can be a different protonation state, or it
might have had a different number in the PDB file and was rebuilt
(it might for instance have been H3, and we only expected H1  H2).
Note that hydrogens might have been added to the entry for the N-terminus.
Remove this hydrogen or choose a different protonation state to solve it.
Option -ignh will ignore all hydrogens in the input.
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors
---


(HB3 would be the hydrogen of a cysteine residue that has not formed a 
disulfide bond with another amino acid residue, however my pdb does not contain 
an HB3).

I then tried adding -ignh as per the suggestion in the twelfth line in that 
message. The following is a part of the message I received:

Checking for duplicate atoms
Generating any missing hydrogen atoms and/or adding termini.
Now there are 29 residues with 385 atoms
Making bonds...
Number of bonds was 390, now 389
Generating angles, dihedrals and pairs...
Before cleaning: 1010 pairs
Before cleaning: 1030 dihedrals
Keeping all generated dihedrals
Making cmap torsions...There are 1030 dihedrals,   73 impropers,  699 angles
   1007 pairs,  389 bonds and 0 virtual sites
Total mass 2916.376 a.m.u.
Total charge -0.000 e
Writing topology

Writing coordinate file...
 - PLEASE NOTE 
You have successfully generated a topology from: 1NB1.pdb.
The Amber99sb force field and the tip3p water model are used.
 - ETON ESAELP 


On the fifth line of that excerpt, I read Number of bonds was 390, now 389 and 
questioned what that could mean. Because of that, I wanted to verify that the topology that was 
successfully generated contains all the necessary bonds, so I looked through the 
topol.top file and saw 3 total hydrogens on the terminal nitrogen on the beginning cysteine and two 
total oxygens on the terminal carbon on the ending valine. Also, I did not see a bond between that 
nitrogen and that carbon (atoms 1 and 383); these data lead me to believe that this peptide is no 
longer cyclic.



This is related to removing duplicate bonds after merging .rtp and .tdb entries. 
 It's normal output; there are almost always duplicates.



Construction of termini and assignment of ionization
state is done with -ter.  For your case, you likely
need to be using -ter and selecting None.



When I input the command pdb2gmx_mpi -ff amber99sb -water tip3p -f 1NB1.pdb -o conf.gro -p topol.top -i 
kalata.itp -ter, I received the same message as before without any kind of interactivity for selecting 
options. Because of this, I tried the command pdb2gmx_mpi -ff amber99sb -water tip3p -f 1NB1.pdb -o 
conf.gro -p topol.top -i kalata.itp -inter -ter, which asked about the protonation of residue 24 
(arginine), and when I answer none, I get this message:

Processing chain 1 'A' (376 atoms, 29 residues)
Which ARGININE type 

[gmx-users] about thermal intergation calculations

2014-06-12 Thread mircial
Hi All:

I am trying to do free energy calculations by thermal integration method and I 
have read the tutorials on the gromacs web site, they are really nice tutorials 
and help me a lot. But I still have two questions here:

1, After we define the necessary paramteres in the mdp file (such as the lambda 
values), when we run MD simulation using these mdp files, we get a file 
(default name dhdl.xvg) of the derivatives  of Hamilton relative to the lambda 
values (the dH/dLAMBDA) as a function of simulation time, and the free energy 
difference is the integration of this value, right?
My question is, how does this value calculated? Since in the simulations, we 
modified some interactions by using a specific lambda value, I can understand 
how to calculate the Hamilton or energy of the system under this specific 
lambda value, but how does the derivatives of Hamilton calculated?

2, General we use different lambda to modify the vdw or electrostatic 
interactions between ligand and it's environment, i understand how this work in 
gromacs. However, I am planning to restrain distances between specific atoms 
and use different lambda to change the equilibrium distances between these 
atoms. Is there any code implemented in gromacs to do this? Or I need to run 
simulations at each lambda value and calculated the potential energy as well as 
the derivatives and integrations by myself?

Thanks very much for your help!

RXG
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Why does the list keep unsubscribing me?

2014-06-12 Thread X Rules
Hi Mark,
 Thanks for your email. Firstly we don't even know if the exact problem is. It 
is just an assumption that this is happening.  I understand that maybe I am the 
only one (or maybe others never got reported). Its just my attitude to try to 
remove bugs when I find and not wait for enough incidents to happen, but that 
might just be difference of opinions.
 PS. A quick look at Mailman http://www.list.org/mailman-member/node41.html; 
tells that confirmation options do exist in Mailman for unsubscribe and it has 
to be switched on (in fact it seems confirmation is by default, but GROMACS 
mailman at KTH has been configured otherwise). But again, I understand that for 
just a single case, it might to too much work to do. 
Have a nice day,cheers,Abhi



 Date: Thu, 12 Jun 2014 13:57:12 +0200
 From: mark.j.abra...@gmail.com
 To: gmx-us...@gromacs.org
 Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
 
 Hi,
 
 I've no idea if mailman even has that option, or the KTH IT admins would be
 willing to set it up, but that's not worth exploring until automatic
 unsubscription is known to be a problem for more than one email address.
 
 Mark
 
 
 On Thu, Jun 12, 2014 at 1:34 PM, X Rules xru...@live.com wrote:
 
  Thanks for the solution Mark, but from what I feel the right way to do is
  to fix mailman (add another option of clicking a link sent to your email to
  confirm unsubscribe, the same way subscription works). I might not be the
  only one, or may be I am the only one complaining :).
 
  Lets see how long this lasts,Abhi
 
   Date: Thu, 12 Jun 2014 10:22:48 +0200
   From: mark.j.abra...@gmail.com
   To: gmx-us...@gromacs.org
   Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
  
   Hi,
  
   One simple solution is to subscribe with a different email address.
  
   Mark
  
  
   On Thu, Jun 12, 2014 at 9:54 AM, X Rules xru...@live.com wrote:
  
This is the 4th time I am thrown out of the list by some random
  phishing
email and I am done with this. If I am thrown out this time (I know I
  will
be and no-one cares to fix any issues), I am not coming back.
I just came back the last time to thank you folks for all your comments
and helps regarding gromacs.have fun !!
   
 From: xru...@live.com
 To: gmx-us...@gromacs.org
 Date: Mon, 2 Jun 2014 04:16:24 +
 Subject: Re: [gmx-users] Why does the list keep unsubscribing me?

 Hi Administrators,   I have been unsubscribed now for the third time
(twice in the same month). This is getting increasingly annoying now.
  Can
there be a fix to mailman to look for genuine email, or to send a
confirmation link to unsubscribe (the same as for subscription) ?
 Thanks, Abhi


  From: xru...@live.com
  To: gmx-us...@gromacs.org
  Date: Thu, 15 May 2014 03:51:37 +
  Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
 
  Thanks for the explanation Rossen.I hope mailman gets some more
artificial intelligence in coming future :).cheers,Abhi
 
   Date: Wed, 14 May 2014 14:33:00 +0200
   From: ros...@kth.se
   To: gmx-us...@gromacs.org
   Subject: Re: [gmx-users] Why does the list keep unsubscribing me?
  
   Hi,
  
   Someone else mentioned the same recently.
  
   Automatic removal can happen if someone, e.g. a spammer, sends a
mail to
   gmx-users-unsubscr...@gromacs.org while faking your email
  address
in the
   From: field.
  
   Mailman is now maintained centrally by the university, so I can't
dig in
   the logs for more information.
  
   Rossen
  
   On 14/05/14 10:29, X Rules wrote:
This is the second time in 2 months, that I am unsubscribed
  from
the list and i have to subscribe again. I just get an email that you
  have
been unsubscribed from the list?
How can I get unsubscribed without my intervention is beyond
understanding?
   
   
  
   --
   Gromacs Users mailing list
  
   * Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
posting!
  
   * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
  
   * For (un)subscribe requests visit
  
  https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users
or send a mail to gmx-users-requ...@gromacs.org.
 
  --
  Gromacs Users mailing list
 
  * Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
posting!
 
  * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 
  * For (un)subscribe requests visit
  https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users
  or
send a mail to gmx-users-requ...@gromacs.org.

 --
 Gromacs Users mailing list

 * Please search the archive at