[gmx-users] Re: protein in SDS/water
Mark Abraham
Mark.Abraham at anu.edu.au
Wed Jun 17 15:18:40 CEST 2009
Anna Marabotti wrote:
> Dear Mark,
> thank you very much for suggestions. I'm pasting here an "extract" of the .top file I used during the grompp
> process (the complete one is too big to be sent to the mailing list, as I told yesterday - the dots indicate
> that I delete the information here, but they are present in the original file). As you can see, I added
> manually the #include for SDS.itp before [ molecules ] and "SDS" under [ moleculetype ] and [ molecules ]
> section. Maybe I made a trivial error, but I can't find it. All lines under [ atoms], [bonds], [pairs],
> [angles] and [dihedral] sections refer only to the protein.
> Moreover, it is quite strange to me that when I'm doing genbox, the top file is not updated automatically. I
> copied the SDS.itp file from PRODRG in the directory ../gromacs/share/top containing all others .itp files,
> but it seems that the program does not find information for SDS to include automatically in the .top file
>
> Many thanks for help
> Kind regards
> Anna
>
> Here the my_prot.top file (extract):
>
> ;
> ; File 'my_prot.top' was generated
> ; By user: anna (62867)
> ; On host: bioserv1.isa.cnr.it
> ; At date: Tue Jun 16 16:15:08 2009
> ;
> ; This is your topology file
> ; "The World is a Friendly Place" (Magnapop)
> ;
> ; Include forcefield parameters
> #include "ffG43a1.itp"
>
> [ moleculetype ]
> ; Name nrexcl
> Protein 3
> SDS 3
This is wrong. [moleculetype] declares the existence of a new type of
molecule, whose name follows. A subsequent [moleculetype] ends the
previous molecule type and starts a new one. Thus only one type of
molecule can be described in a [moleculetype] section. (My previous post
erroneously referred to [molecule] sections, sorry!) The entries in
[molecules] at the end refer to the names given in [moleculetype] sections.
Mark
> [ atoms ]
> ; nr type resnr residue atom cgnr charge mass typeB chargeB massB
> 1 NL 1 SER N 1 0.129 14.0067 ; qtot 0.129
> 2 H 1 SER H1 1 0.248 1.008 ; qtot 0.377
> 3 H 1 SER H2 1 0.248 1.008 ; qtot 0.625
> .........................................................................................
> .........................................................................................
> 2280 C 229 GLU C 980 0.27 12.011 ; qtot -6.73
> 2281 OM 229 GLU O1 980 -0.635 15.9994 ; qtot -7.365
> 2282 OM 229 GLU O2 980 -0.635 15.9994 ; qtot -8
>
> [ bonds ]
> ; ai aj funct c0 c1 c2 c3
> 1 2 2 gb_2
> 1 3 2 gb_2
> 1 4 2 gb_2
> 1 5 2 gb_20
> .......................................................
> .......................................................
> 2275 2276 2 gb_26
> 2276 2277 2 gb_26
> 2277 2278 2 gb_5
> 2277 2279 2 gb_5
> 2280 2281 2 gb_5
> 2280 2282 2 gb_5
>
> [ pairs ]
> ; ai aj funct c0 c1 c2 c3
> 1 7 1
> 1 10 1
> 1 11 1
> ...................................
> ...................................
> 2273 2275 1
> 2273 2280 1
> 2274 2277 1
> 2275 2278 1
> 2275 2279 1
> 2275 2281 1
> 2275 2282 1
> 2276 2280 1
>
> [ angles ]
> ; ai aj ak funct c0 c1 c2 c3
> 2 1 3 2 ga_9
> 2 1 4 2 ga_9
> 2 1 5 2 ga_10
> 3 1 4 2 ga_9
> 3 1 5 2 ga_10
> .....................................
> .....................................
> 2275 2276 2277 2 ga_14
> 2276 2277 2278 2 ga_21
> 2276 2277 2279 2 ga_21
> 2278 2277 2279 2 ga_37
> 2274 2280 2281 2 ga_21
> 2274 2280 2282 2 ga_21
> 2281 2280 2282 2 ga_37
>
> [ dihedrals ]
> ; ai aj ak al funct c0 c1 c2 c3 c4
> c5
> 2 1 5 9 1 gd_14
> 1 5 6 7 1 gd_17
> 1 5 9 11 1 gd_20
> 5 6 7 8 1 gd_12
> 5 9 11 13 1 gd_4
> ..............................................
> ..............................................
> 2270 2272 2274 2280 1 gd_19
> 2272 2274 2275 2276 1 gd_17
> 2272 2274 2280 2282 1 gd_20
> 2274 2275 2276 2277 1 gd_17
> 2275 2276 2277 2279 1 gd_20
>
> [ dihedrals ]
> ; ai aj ak al funct c0 c1 c2 c3
> 5 1 9 6 2 gi_2
> 9 5 11 10 2 gi_1
> 11 9 13 12 2 gi_1
> 14 13 16 15 2 gi_1
> 16 14 18 17 2 gi_1
> 18 16 20 19 2 gi_2
> 20 18 22 21 2 gi_1
> ..........................................
> ..........................................
> 2260 2256 2258 2261 2 gi_1
> 2262 2249 2264 2263 2 gi_1
> 2264 2262 2266 2265 2 gi_1
> 2266 2264 2270 2267 2 gi_2
> 2270 2266 2272 2271 2 gi_1
> 2272 2270 2274 2273 2 gi_1
> 2274 2272 2280 2275 2 gi_2
> 2276 2279 2278 2277 2 gi_1
> 2280 2274 2282 2281 2 gi_1
>
> ; Include Position restraint file
> #ifdef POSRES
> #include "posre.itp"
> #endif
>
> ; Include water topology
> #include "spc.itp"
>
> #ifdef POSRES_WATER
> ; Position restraint for each water oxygen
> [ position_restraints ]
> ; i funct fcx fcy fcz
> 1 1 1000 1000 1000
> #endif
>
> ; Include generic topology for ions
> #include "ions.itp"
>
> ; Include topology for SDS
> #include "SDS.itp"
>
> [ system ]
> ; Name
> Protein in water with SDS
>
> [ molecules ]
> ; Compound #mols
> Protein 1
> SDS 2000
> SOL 36173
>
>> Dear gmx-users,
>> I would like to setup a simulation of a protein in a mixture SDS/water using GROMACS 3.3.1.
>> I proceeded in this way:
>>
>> - after pdb2gmx, I created with editconf a cubic box of 10 nm per side centered on my protein:
>> editconf -f my_prot.gro -o my_prot_box.gro -bt cubic -box 10 -c
>>
>> - then I used genbox to fill the box with water and to add the desired quantity of SDS (to simulate 3.5 M
> SDS
>> I added approx. 2000 molecules of SDS to the box). I obtained the coordinates of SDS from a pdb file and
> then
>> used PRODRG server to create the topology and the coordinate in .gro format. The command was:
>> genbox -cp my_prot_box.gro -cs spc216.gro -ci SDS.gro -nmol 2000 -o my_prot_boxsolv.gro -p my_prot.top
>>
>> The output of this command is a .gro file in which the protein is at the centre of the cubic box, filled
> with
>> approx. 36000 water molecules and 2000 molecules of SDS. I had a look at it with VMD, and all seems to be
> OK.
>> The only "strange" thing is that the SDS molecules are not considered as "solvent" molecules.
>
> Shrug... that depends on the definition of solvent.
>
>> The problem comes when I'm starting to minimize this box. When I use grompp:
>> grompp -f em.mdp -c my_prot_boxsolv.gro -o my_prot_mini.tpr -p my_prot.top
>>
>> the program returns the error:
>> Fatal error:
>> number of coordinates in coordinate file (my_prot_boxsolv.gro, 144801) does not match topology (my_prot.top,
>> 0)
>
> OK, so probably something is wrong with your [molecules] section.
>
>> I tried to manually edit the .top file by adding the number of SDS molecules under the section [ molecules
> ],
>> I also added manually "SDS.itp" to the topology file, but all was useless. I don't really know especially
> why
>> GROMACS claims that the number of coordinates in the topology is 0, since at least the protein is present.
>
> Not if [molecules] is mangled. The order of the directives is quite
> important. The #include for SDS.itp must come before [molecules] and
> after the final subsection of any other [molecule] section. See parts of
> chapter 5 of the manual.
>
> Mark
>
>> What's wrong? I saw many procedures to create mixed solvents for a simulation, but I don't know what is the
>> better one. However, genbox -nmol -ci seems perfect for me and apparently it functions in creating the
> system.
>> So why the topology file is not updated by genbox? Or in case, could you suggest me a better procedure to
> add
>> 3.5 M SDS to my protein in water?
>>
>> Thank you and best regards
>> Anna
>>
>> I tried to copy the .top file at the end of the message, but it is too big, so please if you need to see it
>> let me know how to send to you.
>>
>> ______________________________________________
>> Anna Marabotti, Ph.D.
>> Laboratory of Bioinformatics and Computational Biology
>> Institute of Food Science, CNR
>> Via Roma 64
>> 83100 Avellino (Italy)
>> Phone: +39 0825 299651
>> Fax: +39 0825 781585
>> Skype: annam1972
>> E-mail: amarabotti at isa.cnr.it
>> Web page: http://bioinformatica.isa.cnr.it/anna/anna.htm
>> ____________________________________________________
>> "If you think you are too small to make a difference, try sleeping with a mosquito"
>>
>
> _______________________________________________
> gmx-users mailing list gmx-users at gromacs.org
> http://lists.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 gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>
More information about the gromacs.org_gmx-users
mailing list