[gmx-users] methanol as a solvent

Nuno R. L. Ferreira nunolf at ci.uc.pt
Wed Sep 29 10:31:35 CEST 2004


Hi Alok

IIRC, now you just have to edit your topology file by hand.
Right under the #include describing the peptide, put the .itp of methanol.

Like for e.g.:

; Include forcefield parameters
#include "ffG43a1.itp"
#include "meoh_wfg.itp"  ; check the .itp name

and at the tail of your .top file, add one line with specifying the number
of methanol molecules added to the system.

Methanol 458 ; check the name of [ moleculetype ] under .itp
And you're done.

Best regards.
Nuno


----- Original Message ----- 
From: "Alok" <alokjain at iitk.ac.in>
To: <gmx-users at gromacs.org>
Sent: Wednesday, September 29, 2004 9:15 AM
Subject: [gmx-users] methanol as a solvent


Hello,

         i was trying to solvate  speptide(given in tutotorial of gromacs)
by methanol in place of SPC .
         i did first two steps same as given in tutorial for speptide
(pdb2gmx,editconf).
     For genbox i only change the solvent input file. i specified
meoh_WFG.gro in place of spc216.gro (which is uploaded by "Christoph
Freudenberger" in gromacs topogies. thanks for it )

 genbox -cp out -cs meoh_WFG.gro -p speptide -o b4em

my problem is that solvent(methanol) is not  added  in speptide.top file.
can someone point out my mistake.
Thanks in advance.

i am attaching the output of genbox :-

Opening library file /usr/local/share/gromacs/top/aminoacids.dat
Opening library file /usr/local/share/gromacs/top/atommass.dat
Opening library file /usr/local/share/gromacs/top/vdwradii.dat
Opening library file /usr/local/share/gromacs/top/dgsolv.dat
#Entries in atommass.dat: 82 vdwradii.dat: 26 dgsolv.dat: 7
Reading solute configuration
Go Rough, Oppose Many Angry Chinese Serial killers
Containing 191 atoms in 19 residues
Initialising van der waals distances...
Reading solvent configuration
"MeOH solvent box, 400 molecules, CF 13/03/03"
solvent configuration contains 1200 atoms in 400 residues

Initialising van der waals distances...
Will generate new solvent configuration of 1x2x1 boxes
Generating configuration
Sorting configuration
Found 1 molecule type:
   MeOH (   3 atoms):   800 residues
Calculating Overlap...
box_margin = 0.315
Removed 609 atoms that were outside the box
Table routines are used for coulomb: FALSE
Table routines are used for vdw:     FALSE
Cut-off's:   NS: 0.48   Coulomb: 0.48   LJ: 0.48
Generated table with 0 data points for COUL.
Tabscale = 500 points/nm
Generated table with 0 data points for LJ6.
Tabscale = 500 points/nm
Generated table with 0 data points for LJ12.
Tabscale = 500 points/nm
Going to determine what solvent types we have.
There are 0 molecules, 2591 charge groups and 2591 atoms
There are 0 optimized solvent molecules on node 0
There are 0 optimized water molecules on node 0
Grid: 12 x 19 x 10 cells
Succesfully made neighbourlist
nri = 4450, nrj = 50057
Checking Protein-Solvent overlap: tested 1864 pairs, removed 162 atoms.
Checking Solvent-Solvent overlap: tested 11554 pairs, removed 255 atoms.
Added 458 molecules
Generated solvent containing 1374 atoms in 458 residues
Writing generated configuration to b4em.gro
Go Rough, Oppose Many Angry Chinese Serial killers

Output configuration contains 1565 atoms in 477 residues
Volume                 :     34.9038 (nm^3)
Density                :     731.971 (g/l)
Number of SOL molecules:      0



                                           Alok jain
                                    Indian Institute of Tecnology
                                        Kanpur, India




_______________________________________________
gmx-users mailing list
gmx-users at gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-request at gromacs.org.





More information about the gromacs.org_gmx-users mailing list