francesca vitalini wrote:
Thanks Tsjerk, Florian,
However I need the .gro file for the reverse transformation and I
cannot simply obtain it through pdb2gmx as it won't recognize the
atoms


Use editconf to convert between .pdb and .gro (as well as other formats). pdb2gmx is for producing topologies, not interconverting coordinate files.

-Justin

Fatal error:
Atom BB in residue ALA 1 was not found in rtp entry ALA with 8 atoms
while sorting atoms.
.

Any help with that?
Thanks a lot



2012/2/20 Dommert Florian <domm...@icp.uni-stuttgart.de>:
On Mon, 2012-02-20 at 16:05 +0100, francesca vitalini wrote:
Just  one last question Tsjerk, I was trying to load the cg.gro file
obtained with martinize in vmd and the program says it is unable to
load the molecule. the result of the martinize script looks like:

MODEL        1
TITLE     Protein in water
CRYST1   29.094   29.094   29.094  90.00  90.00  90.00 P 1           1
ATOM      1   BB ALA     1      15.063  17.501  13.232  1.00  0.00
ATOM      2   BB GLU     2      14.568  13.732  14.199  1.00  0.00
ATOM      3  SC1 GLU     2      17.373  14.208  15.223  1.00  0.00
TER
ENDMDL

This looks like a PDB file. If you load cg.gro into VMD it will
recognize it as GRO file, but you provide the wrong format. Rename
cg.gro to cg.pdb and it should be loaded successfully into VMD.

/Flo

Does it sound right to you?
Thanks



2012/2/20 francesca vitalini <francesca.vitalin...@gmail.com>:
Thank you so much.. Now it is working. I was trying to do it without
specifying the secondary structure as I'm having problems with the
do_dssp command. It is not working basically and I have to use pdb2gmx
to convert into pdb then go to the dssp webpage and create the dssp
file from there through the pdb. Do you know a fastest way to make it
work? Thanks

However, my Martini account never worked, don't know why, so I cannot
post it directly.

2012/2/20 Tsjerk Wassenaar <tsje...@gmail.com>:
Hey Francesca,

Now there's a small bug in the program. Sorry about it. We'll put the
fixed version on in a bit. The problem arises because you don't
specify the secondary structure and pymol is not yet available for
doing so. The workaround is to explicitly set the secondary structure
to loop, by adding "-ss LLL" on the command line. Thanks for pointing
that out.

By the way, you might want to post these issues on the Martini forum:
http://md.chem.rug.nl/cgmartini/index.php/user-platform/forum
It's not exactly Gromacs...

Cheers,

Tsjerk

On Mon, Feb 20, 2012 at 2:34 PM, francesca vitalini
<francesca.vitalin...@gmail.com> wrote:
Done while waiting for your e-mail. the error message now says

INFO       Chain termini will be charged
INFO       Residues at chain brakes will not be charged
INFO       Local elastic bonds will be used for extended regions.
INFO       Position restraints will be generated.
WARNING    Position restraints are only enabled if -DPOSRES is set in
the MDP file
INFO       Read input structure from file.
INFO       Input structure is a PDB file.
INFO       Found 1 chains:
INFO          1:     (Protein), 25 atoms in 3 residues.
INFO       Total size of the system: 3 residues.
Traceback (most recent call last):
 File "./martinize-1.0.py", line 2415, in <module>
   elif options["-pymol"]:
KeyError: '-pymol'

and the pdb looks like

TITLE     Protein in water
REMARK    THIS IS A SIMULATION BOX
CRYST1   29.094   29.094   29.094  90.00  90.00  90.00 P 1           1
MODEL        1
ATOM      1  N   ALA     1      15.650  18.630  12.470  1.00  0.00
ATOM      2  H1  ALA     1      16.142  19.465  12.716  1.00  0.00
ATOM      3  H2  ALA     1      16.062  18.229  11.652  1.00  0.00
ATOM      4  H3  ALA     1      14.692  18.849  12.287  1.00  0.00
ATOM      5  CA  ALA     1      15.730  17.670  13.580  1.00  0.00
ATOM      6  CB  ALA     1      15.240  18.220  14.910  1.00  0.00
ATOM      7  C   ALA     1      14.970  16.400  13.180  1.00  0.00
ATOM      8  O   ALA     1      13.880  16.420  12.610  1.00  0.00
ATOM      9  N   GLU     2      15.500  15.310  13.740  1.00  0.00
ATOM     10  H   GLU     2      16.295  15.467  14.326  1.00  0.00
ATOM     11  CA  GLU     2      15.070  13.910  13.610  1.00  0.00
ATOM     12  CB  GLU     2      16.230  12.950  13.880  1.00  0.00
ATOM     13  CG  GLU     2      16.930  13.120  15.230  1.00  0.00
ATOM     14  CD  GLU     2      17.660  14.430  15.540  1.00  0.00
ATOM     15  OE1 GLU     2      18.650  14.820  14.880  1.00  0.00
ATOM     16  OE2 GLU     2      17.070  15.190  16.330  1.00  0.00
ATOM     17  C   GLU     2      13.910  13.480  14.500  1.00  0.00
ATOM     18  O   GLU     2      13.760  12.300  14.810  1.00  0.00
ATOM     19  N   GLU     3      13.000  14.410  14.790  1.00  0.00
ATOM     20  H   GLU     3      13.192  15.329  14.446  1.00  0.00
ATOM     21  CA  GLU     3      11.750  14.260  15.550  1.00  0.00
ATOM     22  CB  GLU     3      11.920  13.690  16.960  1.00  0.00
ATOM     23  C   GLU     3      11.190  15.670  15.730  1.00  0.00
ATOM     24  O1  GLU     3      11.907  16.693  15.192  1.00  0.00
ATOM     25  O2  GLU     3      10.020  15.787  16.413  1.00  0.00
TER


can you please explain which the problem is to me?

thanks



2012/2/20 Tsjerk Wassenaar <tsje...@gmail.com>:
Hi Francesca,

Given the error, it seems there's still a mismatch between the number
of atoms and the number indicated. Try converting your structure to
PDB and use that for coarsegraining.

Cheers,

Tsjerk

On Mon, Feb 20, 2012 at 2:08 PM, francesca vitalini
<francesca.vitalin...@gmail.com> wrote:
I've changed that but it is still complaining...

INFO       Chain termini will be charged
INFO       Residues at chain brakes will not be charged
INFO       Local elastic bonds will be used for extended regions.
INFO       Position restraints will be generated.
WARNING    Position restraints are only enabled if -DPOSRES is set in
the MDP file
INFO       Read input structure from file.
INFO       Input structure is a GRO file. Chains will be labeled consecutively.
Traceback (most recent call last):
 File "./martinize-1.0.py", line 2306, in <module>
   for title,atoms,box in frameIterator(inStream):
 File "./martinize-1.0.py", line 1183, in groFrameIterator
   atoms  = [groAtom(streamIterator.next())  for i in range(natoms)]
 File "./martinize-1.0.py", line 1170, in groAtom
   return (S(a[10:15]), S(a[5:10]),   I(a[:5]), " ",
10*F(a[20:28]),10*F(a[28:36]),10*F(a[36:44]))
ValueError: invalid literal for int() with base 10: '2.'



2012/2/20 Tsjerk Wassenaar <tsje...@gmail.com>:
Hi Francesca,

The problem is that the second line of your gro file indicates there
are 2410 atoms in the file, while there are only 25. Did you manually
remove water? In that case you have to update the number of atoms in
the second line. The error message should be more explanatory though.

Cheers,

Tsjerk

On Mon, Feb 20, 2012 at 1:57 PM, francesca vitalini
<francesca.vitalin...@gmail.com> wrote:
I was using the one available from the tutorial but now I have
downloaded the new one and it gives me still an error message like

INFO       Chain termini will be charged
INFO       Residues at chain brakes will not be charged
INFO       Local elastic bonds will be used for extended regions.
INFO       Position restraints will be generated.
WARNING    Position restraints are only enabled if -DPOSRES is set in
the MDP file
INFO       Read input structure from file.
INFO       Input structure is a GRO file. Chains will be labeled consecutively.
Traceback (most recent call last):
 File "./martinize-1.0.py", line 2306, in <module>
   for title,atoms,box in frameIterator(inStream):
 File "./martinize-1.0.py", line 1183, in groFrameIterator
   atoms  = [groAtom(streamIterator.next())  for i in range(natoms)]
 File "./martinize-1.0.py", line 1170, in groAtom
   return (S(a[10:15]), S(a[5:10]),   I(a[:5]), " ",
10*F(a[20:28]),10*F(a[28:36]),10*F(a[36:44]))
ValueError: invalid literal for int() with base 10: '2.'

I have used the command line ./martinize-1.0.py -f npt_strip_wat.gro
-o frag_npt_cg_pdb.top -x frag_npt_cg_pdb.gro -v  -p backbone

were my input file is very simply

Protein in water
 2410
   1ALA      N    1   1.497   1.862   1.290  0.1206 -0.3069 -0.2102
   1ALA     H1    2   1.409   1.866   1.338  0.8590 -0.1748  1.1978
   1ALA     H2    3   1.547   1.948   1.296  2.0372 -1.4201  0.8359
   1ALA     H3    4   1.483   1.837   1.194  1.8924 -0.0993 -0.5621
   1ALA     CA    5   1.585   1.753   1.335 -0.1231 -0.2165  0.4955
   1ALA     CB    6   1.625   1.781   1.480  0.3554 -0.5405  0.4275
   1ALA      C    7   1.517   1.617   1.323  0.1234 -0.2807 -0.1917
   1ALA      O    8   1.407   1.624   1.268 -0.0033  0.4302  0.1378
   2GLU      N    9   1.566   1.515   1.393 -0.2655 -0.2542  0.1191
   2GLU      H   10   1.655   1.536   1.434  0.1561  0.1079 -0.9482
   2GLU     CA   11   1.501   1.384   1.402  0.3096 -0.5394  0.1823
   2GLU     CB   12   1.601   1.272   1.427  0.7527 -0.2033 -0.0775
   2GLU     CG   13   1.694   1.280   1.548 -0.0279  0.1789  0.5017
   2GLU     CD   14   1.779   1.405   1.571  0.4664 -0.1532  0.4921
   2GLU    OE1   15   1.882   1.412   1.500  0.2904 -0.0189  0.2469
   2GLU    OE2   16   1.740   1.489   1.654  0.2746  0.1001  0.1509
   2GLU      C   17   1.378   1.355   1.489  0.5085 -0.1385  0.5964
   2GLU      O   18   1.352   1.237   1.510  0.3837 -0.1792  0.2200
   3GLU      N   19   1.307   1.458   1.535 -0.1887 -0.4392  0.2009
   3GLU      H   20   1.315   1.547   1.489 -1.1961  0.1298  1.0835
   3GLU     CA   21   1.180   1.436   1.604 -0.3118 -0.2722  0.0307
   3GLU     CB   22   1.192   1.426   1.757 -0.3397  0.1329  0.0602
   3GLU      C   23   1.066   1.533   1.572 -0.3117 -0.4554 -0.5274
   3GLU     O1   24   1.072   1.606   1.471  0.0488  0.2651  0.0047
   3GLU     O2   25   0.965   1.542   1.646 -0.1745 -0.2606 -0.3621
  2.91477   2.91477   2.91477

Can you help me on that?
Thanks

2012/2/20 Tsjerk Wassenaar <tsje...@gmail.com>:
Hi Francesca,

Is this the latest version
(http://md.chem.rug.nl/cgmartini/index.php/downloads/tools/204-martinize)?
If it is, please send me the input file and I'll fix the bug. Note
that the previous version that was available online was one used in a
workshop, while the script was still in beta.

Cheers,

Tsjerk


On Mon, Feb 20, 2012 at 10:43 AM, francesca vitalini
<francesca.vitalin...@gmail.com> wrote:
Hi all,
I'm trying to coarsegrain my structure using the script martinize.py
and using my gro file as inmput and the dssp file with the second
structure downloaded from http://swift.cmbi.ru.nl/gv/dssp using the
pdb structure as input, I get the following error message that I
really don't understand.

INFO       Chain termini will be charged
INFO       Residues at chain brakes will not be charged
INFO       Local elastic bonds will be used for extended regions.
INFO       Position restraints will be generated.
WARNING    Position restraints are only enabled if -DPOSRES is set in
the MDP file
INFO       Read input structure from file.
INFO       Input structure is a GRO file. Chains will be labeled consecutively.
Traceback (most recent call last):
 File "./martinize.py", line 2037, in <module>
   for title,atoms,box in frameIterator(inStream):
 File "./martinize.py", line 1303, in groFrameIterator
   atoms  = [groAtom(streamIterator.next())  for i in range(natoms)]
 File "./martinize.py", line 1290, in groAtom
   return (S(a[10:15]), S(a[5:10]),   I(a[:5]), " ",
10*F(a[20:28]),10*F(a[28:36]),10*F(a[36:44]))
ValueError: invalid literal for int() with base 10: '2.'

 Can anyone help me on that?
Thanks



--
Francesca
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


--
Tsjerk A. Wassenaar, Ph.D.

post-doctoral researcher
Molecular Dynamics Group
* Groningen Institute for Biomolecular Research and Biotechnology
* Zernike Institute for Advanced Materials
University of Groningen
The Netherlands
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


--
Francesca Vitalini

PhD student at Computational Molecular Biology Group,
Department of Mathematics and Informatics, FU-Berlin
Arnimallee 6 14195 Berlin

vital...@zedat.fu-berlin.de
francesca.vital...@fu-berlin.de

+49 3083875776
+49 3083875412
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


--
Tsjerk A. Wassenaar, Ph.D.

post-doctoral researcher
Molecular Dynamics Group
* Groningen Institute for Biomolecular Research and Biotechnology
* Zernike Institute for Advanced Materials
University of Groningen
The Netherlands
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


--
Francesca Vitalini

PhD student at Computational Molecular Biology Group,
Department of Mathematics and Informatics, FU-Berlin
Arnimallee 6 14195 Berlin

vital...@zedat.fu-berlin.de
francesca.vital...@fu-berlin.de

+49 3083875776
+49 3083875412
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


--
Tsjerk A. Wassenaar, Ph.D.

post-doctoral researcher
Molecular Dynamics Group
* Groningen Institute for Biomolecular Research and Biotechnology
* Zernike Institute for Advanced Materials
University of Groningen
The Netherlands
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


--
Francesca Vitalini

PhD student at Computational Molecular Biology Group,
Department of Mathematics and Informatics, FU-Berlin
Arnimallee 6 14195 Berlin

vital...@zedat.fu-berlin.de
francesca.vital...@fu-berlin.de

+49 3083875776
+49 3083875412
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


--
Tsjerk A. Wassenaar, Ph.D.

post-doctoral researcher
Molecular Dynamics Group
* Groningen Institute for Biomolecular Research and Biotechnology
* Zernike Institute for Advanced Materials
University of Groningen
The Netherlands
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


--
Francesca Vitalini

PhD student at Computational Molecular Biology Group,
Department of Mathematics and Informatics, FU-Berlin
Arnimallee 6 14195 Berlin

vital...@zedat.fu-berlin.de
francesca.vital...@fu-berlin.de

+49 3083875776
+49 3083875412


--
Francesca Vitalini

PhD student at Computational Molecular Biology Group,
Department of Mathematics and Informatics, FU-Berlin
Arnimallee 6 14195 Berlin

vital...@zedat.fu-berlin.de
francesca.vital...@fu-berlin.de

+49 3083875776
+49 3083875412
--
Florian Dommert
Dipl. - Phys.

Institute for Computational Physics
University Stuttgart

Pfaffenwaldring 27
70569 Stuttgart

EMail: domm...@icp.uni-stuttgart.de
Homepage: http://www.icp.uni-stuttgart.de/~icp/Florian_Dommert

Tel.: +49 - (0)711 - 68563613
Fax.: +49 - (0)711 - 68563658

--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists




--
========================================

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
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 list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to