[gmx-users] empty velocity *.xvg file from g_traj
Good morning! Firstly, I'd like to thank congratulate everyone involved for ver 4.0 Now for business: I'm getting empty *.xvg files from g_traj... this doesn't always happen, and I can't find what I'm doing differently... I'm running version 3.3.3 on Ubuntu: grompp_d -f eq -c b4eq -n -t mdrun_d -c b4md -v grompp_d -f md -c b4md -n -t mdrun_d -c md -v g_traj_d -n -ox x_ndxgrp1 -ov v_ndxgrp1 gives me two *.xvg files: the x_ndxgrp1 file is fine, but v_ndxgrp1 has an xmgrace formatting header, and then no data afterward. This formatting information agrees with the number of atoms in the group ndxgrp1. what am I doing wrong? -- Inon Sharony ינון שרוני +972(3)6407634 atto.TAU.ac.IL/~inonshar Please consider your environmental responsibility before printing this e-mail. ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] density graph ( density vs time)
huan wrote: Dear all gmx-users and developers. I wish to get a density graph with density versus time instead of nm, so how can i obtained this? Please ask a more specific question. We can't tell whether your problem is knowing what tool to use, how to use it, how to label axes of plots or something else. Mark ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] density graph ( density vs time)
huan wrote: Dear all gmx-users and developers. I wish to get a density graph with density versus time instead of nm, so how can i obtained this? You could extract the box vectors with g_energy and compute the density from the volume. Best, Jochen Thanks. ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php . -- Dr. Jochen Hub Max Planck Institute for Biophysical Chemistry Computational biomolecular dynamics group Am Fassberg 11 D-37077 Goettingen, Germany Email: jhub[at]gwdg.de Tel.: +49 (0)551 201-2312 ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
[gmx-users] dielectric coffiecient
Deat justin thanks for your help. I delete cl- from (peptide,tfe ,water) in my system and I see about his charge not get message. when I add tfe in peptide (pdb2gmx -f name.pdb -ter -ignh )the lysin amino acid is protonated (2 lysin).another charge aminoacid is D,E.do you think it related to adding tfe to system.because my peptide charge is +1.please advise me. and another ,do I need to define e (dielectric E) for water and tfe in my mdp files ?or has the program default? best -- sh-karbalaee ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
[gmx-users] dihedrals again
Hi all, I did some work on thise dihedrals for polystyrene, and I simplified my version now, using also the inspiration from PHE residue. My problem now is that for the dihedral around the CH2-CH bond in the inner monomer (PS thus) I define in the rtp -CH CH2 CH +CH2gd_17 but in the top fileinstead of 2 9 10 17 gd_17 I get this: 2 9 10 11 (without any parameter). This make me think that pdb2gmx does not count the dihedral as I defined it (why?) but puts a default dihedral there, without any parameter. I am right? What would be a solution for that? My dihedrals are now defined as follows (it helps if you really draw the molecule): the entries in the rtp file for first (PS1), inner (PS) and last (PSN) monomer: [PS1] [dihedrals] ;in the chain - 3-fold (n=3) CH3 CH +CH2+CH gd_17 ;around CH-CB - 2-fold (n=2) CH3 CH CB CG5 gd_1 [PS] [dihedrals] ;in the chain - 3-fold CH2 CH +CH2+CH gd_17 -CH CH2 CH +CH2gd_17 ;around CH-CB - 2-fold CH2 CH CB CG5 gd_1 [PSN] [dihedrals] ;in the chain - 3-fold (n=3) -CH CH2 CH CH3 gd_17 ;around CH-CB - 2-fold (n=2) CH2 CH CB CG5 gd_1 The results for trimer in the top file. [ dihedrals ] ; aiajakal functc0c1 c2c3c4c5 1 2 3 8 1gd_1 1 2 910 1gd_17 2 91011 1 9101116 1gd_1 9101718 1gd_17 10171825 1gd_17 17181924 1gd_1 Best regards, Andrea ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
[gmx-users] Re: dihedrals again
Remark: if I change that dihedral with -CHCH2 CH CBgd_17, that it gives again 2 9 10 11 , but this time with gd_17, as it should do. But why it does not work with the other definition? Does it have to do with the fact that in the dihedral I have atoms belonging to 3 residues? And are the 2 definitions (with +CH2 and CB at the end) equivalent, taking in acount that we define also an improper of type gi_2 around the CH atom? Cheers, Andrea ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] dielectric coffiecient
shahrbanoo karbalaee wrote: Deat justin thanks for your help. I delete cl- from (peptide,tfe ,water) in my system and I see about his charge not get message. when I add tfe in peptide (pdb2gmx -f name.pdb -ter -ignh )the lysin amino acid is protonated (2 lysin).another charge aminoacid is D,E.do you think it related to adding tfe to system.because my peptide charge is +1.please advise me. and another ,do I need to define e (dielectric E) for water and tfe in my mdp files ?or has the program default? If your peptide has a net charge of +1, then you need a counterion to neutralize the charge. But I'm confused. If you have 2 Lys and an Asp and Glu, then the net charge should be 0, unless you've done something strange with your termini. There is no need to set a dielectric, at least not with PME for electrostatics. If you're using reaction field, then you would specify a dielectric beyond a certain cut-off. -Justin best -- Justin A. Lemkul Graduate Research Assistant Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] Adding Na ions to DPPG
On Mon, Oct 27, 2008 at 10:21 AM, Justin A. Lemkul [EMAIL PROTECTED] wrote: Peyman Yamin wrote: On Monday 27 October 2008 15:07, Justin A. Lemkul wrote: Peyman Yamin wrote: On Monday 27 October 2008 14:03, Justin A. Lemkul wrote: Actually I'm doing a funny mistake telling I should replace an O with an Na! (probably because of high adrenaline level in my blood for the deadline is near!) I should just put one Na beside the O I have! I simply have a PO4 where one O has no H and thus is negative. an Na+ should accompany this O. Well that makes more sense. I'm assuming you have a DPPG bilayer with solvent? Or are you simulating a single DPPG in vacuo? If it's a bilayer in solution, just use genion; you can even specify exactly how many Na+ you want using -np. If the topology organization is confusing you, refer to Chapter 5 of the manual. If you make a .top for DPPG, you can specify: [ molecules ] DPPGx NA+ 1 (or whatever) after including the appropriate force field call and ions.itp; the parameters for all atom types are taken from these files. Well this looks great! I simulate one single DPPG in solvent. and you know, this Na+ is not in solvent but a part of DPPG. So you mean I just make a DPPG with an O which has no H so is minus and then use genion to add Na? will it place the Na aroung O hopefully? I think if I put no H arounf O the prodrg will force me to have it! so I should delete it and replace it with Na? or I should never use the word replace at all??? Well, the Na+ may be associated with DPPG electrostatically, but they are still distinct chemical entities. You can have a DPPG that has a net -1 charge; certainly all chemical species are not net neutral! Indeed, genion will place Na+ in whatever solvent you choose, but may not necessarily be bound to DPPG. Some simulation may show the association, however. Yes, keep the DPPG. The charge distribution around anionic lipid has been studied. http://scitation.aip.org/getabs/servlet/GetabsServlet?prog=normalid=PRLTAO0001010303810301idtype=cvipsgifs=yes And no, don't replace any part of your DPPG. Because then you don't really have DPPG any more do you? -Justin Thanks again ;) Peyman -Justin Peyman Yamin wrote: Dear GMX users, I want to add sodium ions (Na+) to a phospholipid (DPPG) I have the structure with the Na replaced with O. I get a .top with this from prodrg and try to add the Na manually at the right place by removing the corresponding O. I use GROMOS 96.1 forcefield. GMX version 3.3.1.I have a problem: Are you trying to replace the O in DPPG with Na+? Or are you simply trying to add ions to the surrounding solvent (water)? If you just need ions in the solvent, use genion; you don't have to do it manually. * In .top file, as far as I can understand, I should just replace OA with NA+ and set the mass and charge. DPPG has a centre: [ O-P(O)=O ] - The Na+ should be around this with +1 charge. should I put charge +1 ? Again, I'm not clear on what you're trying to accomplish. If you replace one of these oxygen atoms, you will likely not have an integer charge within the DPPG molecule itself. * which parameters should I place in the top as the Na+ is not bonding. I need LJ and Coulomb params. I have C6 and C12 them from ffG43a1nb.itp. (NA+ OA) but I'm not sure where should I put it in the top file. What should I insert manually in top file at all? Simply #include ffG43a1.itp at the top of your .top file; it will include nonbonded and bonded parameters for all the atoms in the system. -Justin Thanks in advance for the time, I've been away for some time from GMX and might sound Cheers, Peyman -- Justin A. Lemkul Graduate Research Assistant Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php -- Best wishes, Myunggi Yi == KLB 419 Institute of Molecular Biophysics Florida State University Tallahassee, FL 32306 Office: (850) 645-1334 http://www.scs.fsu.edu/~myunggi [EMAIL PROTECTED] ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL
Re: [gmx-users] density graph ( density vs time)
huan wrote: Dear all, Previously i obtained a density (y-axis) versus box length (x axis) using g_density.. So now i wonder is it there is another way that we can use other way to obtain a density (y axis) versus time length (x-axis). Use g_energy. -Justin Thanks.. --- On Tue, 10/28/08, Mark Abraham [EMAIL PROTECTED] wrote: From: Mark Abraham [EMAIL PROTECTED] Subject: Re: [gmx-users] density graph ( density vs time) To: Discussion list for GROMACS users gmx-users@gromacs.org Date: Tuesday, October 28, 2008, 4:46 PM huan wrote: Dear all gmx-users and developers. I wish to get a density graph with density versus time instead of nm, so how can i obtained this? Please ask a more specific question. We can't tell whether your problem is knowing what tool to use, how to use it, how to label axes of plots or something else. Mark ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php -- Justin A. Lemkul Graduate Research Assistant Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] density graph ( density vs time)
I used g_energy before but i just obtained density (y-axis) versus box length (x axis). now i wish to obtain density vs time thanks --- On Tue, 10/28/08, Justin A. Lemkul [EMAIL PROTECTED] wrote: From: Justin A. Lemkul [EMAIL PROTECTED] Subject: Re: [gmx-users] density graph ( density vs time) To: [EMAIL PROTECTED], Discussion list for GROMACS users gmx-users@gromacs.org Date: Tuesday, October 28, 2008, 11:31 PM huan wrote: Dear all, Previously i obtained a density (y-axis) versus box length (x axis) using g_density.. So now i wonder is it there is another way that we can use other way to obtain a density (y axis) versus time length (x-axis). Use g_energy. -Justin Thanks.. --- On Tue, 10/28/08, Mark Abraham [EMAIL PROTECTED] wrote: From: Mark Abraham [EMAIL PROTECTED] Subject: Re: [gmx-users] density graph ( density vs time) To: Discussion list for GROMACS users gmx-users@gromacs.org Date: Tuesday, October 28, 2008, 4:46 PM huan wrote: Dear all gmx-users and developers. I wish to get a density graph with density versus time instead of nm, so how can i obtained this? Please ask a more specific question. We can't tell whether your problem is knowing what tool to use, how to use it, how to label axes of plots or something else. Mark ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php -- Justin A. Lemkul Graduate Research Assistant Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] density graph ( density vs time)
huan wrote: I used g_energy before but i just obtained density (y-axis) versus box length (x axis). huan wrote: Previously i obtained a density (y-axis) versus box length (x axis) using g_density.. Sounds like a contradition to me... :-) thanks --- On Tue, 10/28/08, Justin A. Lemkul [EMAIL PROTECTED] wrote: From: Justin A. Lemkul [EMAIL PROTECTED] Subject: Re: [gmx-users] density graph ( density vs time) To: [EMAIL PROTECTED], Discussion list for GROMACS users gmx-users@gromacs.org Date: Tuesday, October 28, 2008, 11:31 PM huan wrote: Dear all, Previously i obtained a density (y-axis) versus box length (x axis) using g_density.. So now i wonder is it there is another way that we can use other way to obtain a density (y axis) versus time length (x-axis). Use g_energy. -Justin Thanks.. --- On Tue, 10/28/08, Mark Abraham [EMAIL PROTECTED] wrote: From: Mark Abraham [EMAIL PROTECTED] Subject: Re: [gmx-users] density graph ( density vs time) To: Discussion list for GROMACS users gmx-users@gromacs.org Date: Tuesday, October 28, 2008, 4:46 PM huan wrote: Dear all gmx-users and developers. I wish to get a density graph with density versus time instead of nm, so how can i obtained this? Please ask a more specific question. We can't tell whether your problem is knowing what tool to use, how to use it, how to label axes of plots or something else. Mark ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php -- Justin A. Lemkul Graduate Research Assistant Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php . -- Dr. Jochen Hub Max Planck Institute for Biophysical Chemistry Computational biomolecular dynamics group Am Fassberg 11 D-37077 Goettingen, Germany Email: jhub[at]gwdg.de Tel.: +49 (0)551 201-2312 ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] density graph ( density vs time)
Jochen Hub wrote: huan wrote: I used g_energy before but i just obtained density (y-axis) versus box length (x axis). huan wrote: Previously i obtained a density (y-axis) versus box length (x axis) using g_density.. Sounds like a contradition to me... :-) Indeed; the g_energy output gives the x-axis as time... -Justin thanks --- On Tue, 10/28/08, Justin A. Lemkul [EMAIL PROTECTED] wrote: From: Justin A. Lemkul [EMAIL PROTECTED] Subject: Re: [gmx-users] density graph ( density vs time) To: [EMAIL PROTECTED], Discussion list for GROMACS users gmx-users@gromacs.org Date: Tuesday, October 28, 2008, 11:31 PM huan wrote: Dear all, Previously i obtained a density (y-axis) versus box length (x axis) using g_density.. So now i wonder is it there is another way that we can use other way to obtain a density (y axis) versus time length (x-axis). Use g_energy. -Justin Thanks.. --- On Tue, 10/28/08, Mark Abraham [EMAIL PROTECTED] wrote: From: Mark Abraham [EMAIL PROTECTED] Subject: Re: [gmx-users] density graph ( density vs time) To: Discussion list for GROMACS users gmx-users@gromacs.org Date: Tuesday, October 28, 2008, 4:46 PM huan wrote: Dear all gmx-users and developers. I wish to get a density graph with density versus time instead of nm, so how can i obtained this? Please ask a more specific question. We can't tell whether your problem is knowing what tool to use, how to use it, how to label axes of plots or something else. Mark ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php -- Justin A. Lemkul Graduate Research Assistant Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php . -- Justin A. Lemkul Graduate Research Assistant Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
[gmx-users] table use
Hi all users, When I input the mdrun -table table_nonbond.xvg it always shows that : Fatal error: Library file tablep.xvg not found in current dir nor in default directories. (You can set the directories to search with the GMXLIB path variable) Has anyone encountered the same problem? Can you share me your experience about how to solve it ? Thank you very much. Yang ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
[gmx-users] Is there any realation between structure and RMSD?
Hi all I have simulated two systems(mutated and unmutated)of protein for 7ns and I plotted RMSD. here my doubt is that when I see these structures(superimpose)in VMD mutated final simulated structure with respective initial one drastic changes(means helix-coil and Beta-alpha transitions) but there is no change in unmutated structure with corresponding initial one. But when I plot RMSD, mutated structure showing low fluctuations(0.19nm)while unmutated(RMSD 0.3nm)high fluctuations. I have mentined RMSD command like this g_rms -f 7ns.xtc -s min.tpr -o 7ns_rmsd -pbc Selected c-alpha for both least square fit and RMSD calculation What I mean to say that the structure which shows changes(coil-helix and beta-alpha)it is supposed to has RMSD high viceversa but I am getting in opposite. Why its happening? Any suggestions would be appreciated Thanks in advance. ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] cannot generate .top file
Hi. Thanks for your enlightenment. Now, I'm able to run pdb2gmx with the output files produced. I just added 'sudo' before ./pdb2gmx and it allows me to write in the directory. But unfortunately, this step seem does not work for ./gmxtest.pl. I'm having the same output even after i had added 'sudo' in front of it. Joyce On Fri, Oct 24, 2008 at 8:52 PM, Mark Abraham [EMAIL PROTECTED]wrote: Kwee Hong wrote: Hi. I seem having some problem in running ./pdb2gmx. Here is the output file. Can you help me in giving me some idea to deal with this? bcsb09s-imac52:bin bcsb09$ pdb2gmx -f 1aml.pdb -p 1aml.top -o 1aml.gro snip --- Program pdb2gmx, VERSION 3.3.3 Source code file: futil.c, line: 313 File input/output error: 1aml.top --- This means it can't write 1aml.top, which likely means you don't have write permission to the directory in which pdb2gmx resides. This is a Good Thing. You should do daily work in a suitable subdirectory of (say) your home directory, and adjust your environment variables to give you ready access to GROMACS stuff by using source /path/to/GMXRC using the path in which you installed the GROMACS binaries. I guess this is your attempted working directory above. I seem to be having problem in running the test tun as well. Below is its output: bcsb09s-imac52:bin bcsb09$ ./gmxtest.pl Usage: ./gmxtest.pl [ -np N ] [-verbose ] [ -double ] [ simple | complex | kernel | pdb2gmx | all ] or: ./gmxtest.pl clean | refclean | dist Again this won't work in a directory in which you don't have write permissions. The above message tells you that an empty command line isn't acceptable to gmxtest.pl. *Technically* it ought to be acceptable per the usage line, since square brackets denote optional things and |-symbols denote a series of mutually-exclusive possibilities. I think the brackets should be missing around simple | complex | kernel | pdb2gmx | all. Mark ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] cannot generate .top file
Kwee Hong wrote: Hi. Thanks for your enlightenment. Now, I'm able to run pdb2gmx with the output files produced. I just added 'sudo' before ./pdb2gmx and it allows me to write in the directory. But unfortunately, this step seem does not work for ./gmxtest.pl. I'm having the same output even after i had added 'sudo' in front of it. If you do this kind of thing, your machine is a ticking time bomb. I gave you good advice last time. I'm not going to give you more :-) Mark ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
[gmx-users] distance restraint
Dear gmx's user, May I asked you some technicl Questions? Do you know how to make distance restraint in gromacs? what I did I have created file name disre.itp as shown below [ distance_restraints ] ; ai aj type index type low up1 up2 fac 489 1133 1 1 1 0.25 0.35 0.40 1.0 489 1164 1 2 1 0.25 0.35 0.40 1.0 490 1127 1 3 1 0.25 0.35 0.40 1.0 490 1167 1 4 1 0.25 0.35 0.40 1.0 in md.mdp file: title = M2_amand in POPC cpp = /lib/cpp constraints = all-bonds define = -DDISRES integrator = md dt = 0.002 ; ps ! nsteps = 50 ; total 1000 ps = 1 ns ! nstcomm = 1 nstxout = 250 nstvout = 1000 nstfout = 0 nstlog = 100 nstenergy = 100 nstlist = 10 ns_type = grid rlist = 1.2 coulombtype = PME rcoulomb = 1.2 rvdw = 1.2 pbc = xyz ; Berendsen temperature coupling is on in two groups Tcoupl = berendsen tc-grps = Protein POPC SOL Cl tau_t = 0.1 0.1 0.1 0.1 ref_t = 310 310 310 310 ; Energy monitoring energygrps = Protein POPC SOL Cl ; Isotropic pressure coupling is now on Pcoupl = berendsen Pcoupltype = isotropic tau_p = 1.0 compressibility = 4.5e-5 ref_p = 1.0 ; Generate velocites is off at 300 K. gen_vel = no gen_temp = 310.0 disres = simple in topology.top file I added... ; Include distance restraint file #ifdef DISRES #include disre.itp #endif My questions are Q1:What is column index mean in disre.itp file? and how is it benefit? Q2:Where can I find in the output to make sure that my distances has already restrainted? Thanks so much pathum ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] cannot generate .top file
Oo.. I'm so sorry for I'm not aware of this. But I'm really new to linux environment and don't really have a good mastering skill on the commandline that it works with. As you said that I'm not given the permission to write in that directory, sudo was the only thing that I can think of out from my knowledge. I don't really know how to do my work at the subdirectory of my hoe directory. And does the meaning of adjust your environment variables to give you ready access to GROMACS stuff by using source /path/to/GMXRC using the path in which you installed the GROMACS binaries requires me to make some changes in the source code or just by typing source /path/to/GMXRC at the terminal would have open the access for me in the directory? I'm sorry for my ignorance... I think I would better get some books on linux environment and read on it first. Joyce On Wed, Oct 29, 2008 at 12:11 PM, Mark Abraham [EMAIL PROTECTED]wrote: Kwee Hong wrote: Hi. Thanks for your enlightenment. Now, I'm able to run pdb2gmx with the output files produced. I just added 'sudo' before ./pdb2gmx and it allows me to write in the directory. But unfortunately, this step seem does not work for ./gmxtest.pl. I'm having the same output even after i had added 'sudo' in front of it. If you do this kind of thing, your machine is a ticking time bomb. I gave you good advice last time. I'm not going to give you more :-) Mark ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED] Can't post? Read http://www.gromacs.org/mailing_lists/users.php