VMD-L Mailing List

From: Axel Kohlmeyer (akohlmey_at_cmm.chem.upenn.edu)
Date: Wed Jan 10 2007 - 08:39:17 CST

" politr"@huji.ac.il wrote:
> Dear Peter,
> Thank you very much for your fast replay. I want to explain you why I'm
> using only one atom and maybe you can advise me some other way to
> perform the same thing. My target is to create water box with 2M of
> sorbitol (150 mol). I was planning to do it using VMD. First of all I
> create water box and after that I add needed amount of ions (using add
> ions tool from VMD). After that the ions are replaced with atom name and
> residue name of sorbitol and by using psfgen I was planning to create my
> sorbitol molecule. At the moment I did it I can solvate the sorbitol
> molecules and at the end I'm getting water box with needed amount of
> sorbitol molecules. Using this procedure I can't put into water box more
> than one atom for each sorbitol molecule and I don't have a clue how to
> perform it. I will appreciate it very much if you or anyone else can

hi,

an alternate strategy would be to first compute the sorbitol / water
molecule ratio you need, then use the solvate plugin to create a single
solvated sorbitol molecule (you can remove waters manually, if you get
too many), equilibrate this and then just create as many properly translated
copies of that box as you need, merge them into a large file and use a text
editor (or some shell commands) to sort the sorbitols to the front.

i have not done this with the NAMD/VMD toolchain, but with gromacs, and
it worked fine. as a benefit, you get a (too?) evenly distributed mixture
and the equilibration is much faster, as each local sorbitol environment
is pre-equilibrated when you put together the whole system.

cheers,
    axel.

> help me. This procedure worked for glycerol but as I see not for sorbitol.
> Thanks again
> Regina
> Peter Freddolino wrote:
>>
>> ------------------------------------------------------------------------
>>
>> Subject:
>> Re: vmd-l: create new topology file for new molecule
>> From:
>> Peter Freddolino <petefred_at_ks.uiuc.edu>
>> Date:
>> Tue, 09 Jan 2007 11:10:42 -0600
>> To:
>> politr_at_huji.ac.il
>>
>> To:
>> politr_at_huji.ac.il
>>
>>
>> Hi Regina,
>> psfgen does indeed pay attention to internal coordinates, but you need
>> to give it more than one atom as a starting point or else it will just
>> need to guess. Try, for example, giving it C1/C2/C3, or all of the
>> carbons, and you'll see the effects of the IC lines (note that more
>> information is needed in the file you gave -- some bonds and angles
>> would be useful).
>>
>> Peter
>>
>> " politr"@huji.ac.il wrote:
>>
>>> Dear NAMD and VMD users,
>>> I'm trying to build topology file for sorbitol. I have used glycerol
>>> topology file in order to do it. I'm using psfgen with this topology
>>> file in order to create pdb and psf. The problem is that the molecule
>>> created with this topology file is not sorbitol but its isomer. I have
>>> tried to change the Internal coordinate definitions in topology file
>>> but it looks like psfgen ignores this lines at all and creates the
>>> molecule based on its own default parameters. The files i'm using are
>>> attached. I hope somebody can help me with it.
>>> Thank you in advance
>>> Regina
>>> The attached files:
>>> sorbitol_alone.pdb - input file to psfgen in oder to create sorbitol
>>> using only one coordinate (guesscoord can do it)
>>> one_sorbitol.psf/pdb - output files from psfgen
>>> topology file - sorbitol definitions at the end of file
>>> ------------------------------------------------------------------------
>>>
>>> PSF
>>>
>>> 3 !NTITLE
>>> REMARKS original generated structure x-plor psf file
>>> REMARKS topology ../inp/top_all27_prot_lipid_ret_csff_adam_spc.rtf
>>> REMARKS segment SORB { first NONE; last NONE; auto angles dihedrals }
>>>
>>> 26 !NATOM
>>> 1 SORB 1 SORB C1 CTS 0.050000 12.0110 0
>>> 2 SORB 1 SORB H11 HAS 0.100000 1.0080 0
>>> 3 SORB 1 SORB H12 HAS 0.100000 1.0080 0
>>> 4 SORB 1 SORB O1 OHS -0.650000 15.9994 0
>>> 5 SORB 1 SORB H1 HOS 0.400000 1.0080 0
>>> 6 SORB 1 SORB C2 CTS 0.150000 12.0110 0
>>> 7 SORB 1 SORB H21 HAS 0.100000 1.0080 0
>>> 8 SORB 1 SORB O2 OHS -0.650000 15.9994 0
>>> 9 SORB 1 SORB H2 HOS 0.400000 1.0080 0
>>> 10 SORB 1 SORB C3 CTS 0.150000 12.0110 0
>>> 11 SORB 1 SORB H31 HAS 0.100000 1.0080 0
>>> 12 SORB 1 SORB O3 OHS -0.650000 15.9994 0
>>> 13 SORB 1 SORB H3 HOS 0.400000 1.0080 0
>>> 14 SORB 1 SORB C4 CTS 0.150000 12.0110 0
>>> 15 SORB 1 SORB H41 HAS 0.100000 1.0080 0
>>> 16 SORB 1 SORB O4 OHS -0.650000 15.9994 0
>>> 17 SORB 1 SORB H4 HOS 0.400000 1.0080 0
>>> 18 SORB 1 SORB C5 CTS 0.150000 12.0110 0
>>> 19 SORB 1 SORB H51 HAS 0.100000 1.0080 0
>>> 20 SORB 1 SORB O5 OHS -0.650000 15.9994 0
>>> 21 SORB 1 SORB H5 HOS 0.400000 1.0080 0
>>> 22 SORB 1 SORB C6 CTS 0.050000 12.0110 0
>>> 23 SORB 1 SORB H61 HAS 0.100000 1.0080 0
>>> 24 SORB 1 SORB H62 HAS 0.100000 1.0080 0
>>> 25 SORB 1 SORB O6 OHS -0.650000 15.9994 0
>>> 26 SORB 1 SORB H6 HOS 0.400000 1.0080 0
>>>
>>> 25 !NBOND: bonds
>>> 1 2 1 3 1 4 1 6
>>> 4 5 6 7 6 8 6 10
>>> 8 9 10 11 10 12 10 14
>>> 12 13 14 15 14 16 14 18
>>> 16 17 18 19 18 20 18 22
>>> 20 21 22 23 22 24 22 25
>>> 25 26
>>>
>>> 42 !NTHETA: angles
>>> 1 6 10 1 6 8 1 6 7
>>> 1 4 5 2 1 6 2 1 4
>>> 2 1 3 3 1 6 3 1 4
>>> 4 1 6 6 10 14 6 10 12
>>> 6 10 11 6 8 9 7 6 10
>>> 7 6 8 8 6 10 10 14 18
>>> 10 14 16 10 14 15 10 12 13
>>> 11 10 14 11 10 12 12 10 14
>>> 14 18 22 14 18 20 14 18 19
>>> 14 16 17 15 14 18 15 14 16
>>> 16 14 18 18 22 25 18 22 24
>>> 18 22 23 18 20 21 19 18 22
>>> 19 18 20 20 18 22 22 25 26
>>> 23 22 25 23 22 24 24 22 25
>>>
>>> 63 !NPHI: dihedrals
>>> 1 6 8 9 1 6 10 11
>>> 1 6 10 12 1 6 10 14
>>> 2 1 4 5 2 1 6 7
>>> 2 1 6 8 2 1 6 10
>>> 3 1 4 5 3 1 6 7
>>> 3 1 6 8 3 1 6 10
>>> 4 1 6 7 4 1 6 8
>>> 4 1 6 10 5 4 1 6
>>> 6 10 12 13 6 10 14 15
>>> 6 10 14 16 6 10 14 18
>>> 7 6 8 9 7 6 10 11
>>> 7 6 10 12 7 6 10 14
>>> 8 6 10 11 8 6 10 12
>>> 8 6 10 14 9 8 6 10
>>> 10 14 16 17 10 14 18 19
>>> 10 14 18 20 10 14 18 22
>>> 11 10 12 13 11 10 14 15
>>> 11 10 14 16 11 10 14 18
>>> 12 10 14 15 12 10 14 16
>>> 12 10 14 18 13 12 10 14
>>> 14 18 20 21 14 18 22 23
>>> 14 18 22 24 14 18 22 25
>>> 15 14 16 17 15 14 18 19
>>> 15 14 18 20 15 14 18 22
>>> 16 14 18 19 16 14 18 20
>>> 16 14 18 22 17 16 14 18
>>> 18 22 25 26 19 18 20 21
>>> 19 18 22 23 19 18 22 24
>>> 19 18 22 25 20 18 22 23
>>> 20 18 22 24 20 18 22 25
>>> 21 20 18 22 23 22 25 26
>>> 24 22 25 26
>>>
>>> 0 !NIMPHI: impropers
>>>
>>>
>>> 0 !NDON: donors
>>>
>>>
>>> 0 !NACC: acceptors
>>>
>>>
>>> 0 !NNB
>>>
>>> 0 0 0 0 0 0 0 0
>>> 0 0 0 0 0 0 0 0
>>> 0 0 0 0 0 0 0 0
>>> 0 0
>>>
>>> 1 0 !NGRP
>>> 0 0 0
>>>
>>>
>>>
>>
>>
> 

-- 
=======================================================================
Axel Kohlmeyer   akohlmey_at_cmm.chem.upenn.edu   http://www.cmm.upenn.edu
   Center for Molecular Modeling   --   University of Pennsylvania
Department of Chemistry, 231 S.34th Street, Philadelphia, PA 19104-6323
tel: 1-215-898-1582,  fax: 1-215-573-6233,  office-tel: 1-215-898-5425
=======================================================================
If you make something idiot-proof, the universe creates a better idiot.