[gmx-users] Re: question about minimisation
Mark Abraham
Mark.Abraham at anu.edu.au
Thu Jul 26 13:15:52 CEST 2012
On 26/07/2012 9:06 PM, reisingere at rostlab.informatik.tu-muenchen.de wrote:
>>
>> On 7/26/12 6:07 AM, reisingere at rostlab.informatik.tu-muenchen.de wrote:
>>>> On 26/07/2012 6:47 PM, reisingere at rostlab.informatik.tu-muenchen.de
>>>> wrote:
>>>>> Ho,
>>>>> first I minimize my structure. This is the corresponding mdp file:
>>>>>
>>>>> define = -DPOSRES
>>>>> integrator = steep
>>>>> emtol = 10
>>>>> nsteps = 1500
>>>>> nstenergy = 1
>>>>> energygrps = System
>>>>> coulombtype = PME
>>>>> rcoulomb = 0.9
>>>>> rvdw = 0.9
>>>>> rlist = 0.9
>>>>> fourierspacing = 0.12
>>>>> pme_order = 4
>>>>> ewald_rtol = 1e-5
>>>>> pbc = xyz
>>>>>
>>>>>
>>>>> and then I run a md run. This is the corresponding mdp file:
>>>>>
>>>>> define = -DPOSRES
>>>>> integrator = md
>>>>> dt = 0.001
>>>>> nsteps = 5000
>>>>> nstxout = 100
>>>>> nstvout = 0
>>>>> nstfout = 0
>>>>> nstlog = 1000
>>>>> nstxtcout = 500
>>>>> nstenergy = 5
>>>>> energygrps = Protein Non-Protein
>>>>> nstcalcenergy = 5
>>>>> nstlist = 10
>>>>> ns-type = Grid
>>>>> pbc = xyz
>>>>> rlist = 0.9
>>>>> coulombtype = PME
>>>>> rcoulomb = 0.9
>>>>> rvdw = 0.9
>>>>> fourierspacing = 0.12
>>>>> pme_order = 4
>>>>> ewald_rtol = 1e-5
>>>>> gen_vel = yes
>>>>> gen_temp = 200.0
>>>>> gen_seed = 9999
>>>>> constraints = all-bonds
>>>>> tcoupl = V-rescale
>>>>> tc-grps = Protein Non-Protein
>>>>> tau_t = 0.1 0.1
>>>>> ref_t = 298 298
>>>>> pcoupl = no
>>>>>
>>>>>
>>>>>
>>>>> In my topology file I include the restraint files like this:
>>>>>
>>>>> ; Include Position restraint file
>>>>> #ifdef POSRES
>>>>> #include "posre.itp"
>>>>> #endif
>>>>>
>>>>> #ifdef POSRES
>>>>> #include "posre_memb.itp"
>>>>> #endif
>>>> This won't work for multiple [moleculetype] entries. See
>>>> http://www.gromacs.org/Documentation/How-tos/Position_Restraints
>>>>
>>>>>
>>>>> I just recognize that there is a "DPOSRES" in the mdp files and a
>>>>> "POSRES"
>>>>> in my topology file. Is this the problem. Do I have to write it the
>>>>> same
>>>>> way in the several files?
>>>> http://www.gromacs.org/Documentation/Include_File_Mechanism
>>>>
>>>> Mark
>>>
>>> I just see that the [ position_restraints ] part is included under the
>>> [dihedral] and not between the [moleculetype] and the [atom] part. And
>>> according to the site you wrote me this is a problem, right? But this
>>> was
>>> done by gromacs itself. Shell I write it to the [moleculetype] part?
>>>
>> A [position_restraints] directive belongs to the [moleculetype] in which
>> it is
>> declared. The original location of the #include statement produced by
>> Gromacs
>> is correct; it follows sequentially within the protein [moleculetype].
>> Your
>> inclusion of the membrane restraint file within the protein [moleculetype]
>> is,
>> however, incorrect.
>>
>>> But I can not see why it can not work that I have 2 restriction files?
>>> Can
>>> you please explain it to me?
>>>
>> You can have two restraint files for different [moleculetypes] but they
>> must be
>> organized as such.
> okey. So the whole thing with [atoms], [bonds], [dihedral] all contains to
> the entry in [moleculetype] , right?
>
> But directly before I include the membrane restriction file I include the
> membrane definition:
>
> #include "amber03.ff/dum.itp"
> #ifdef POSRES
> #include "posre_memb.itp"
> #endif
>
>
> So I thought that it is directly after the atomtype it belongs to.
Atoms and molecules are different things. So are [atomtype] and
[moleculetype]. Justin, the webpage and I have only been talking about
[moleculetype].
> I
> thought that it is the same in the case with the water, where first the
> water definition is included and after that the restriction of the water.
>
> Or am I wrong?
>
> I include the dummy atoms definition right after the ions. Or is this the
> wrong position?
>
>
> 4935 4952 4950 4951 4
> 4950 4954 4952 4953 4
> 4954 4973 4972 4974 4
>
> ; Include Position restraint file
> #ifdef POSRES
> #include "posre.itp"
> #endif
>
> ; Include water topology
> #include "amber03.ff/tip3p.itp"
>
> #ifdef POSRES_WATER
> ; Position restraint for each water oxygen
> [ position_restraints ]
> ; i funct fcx fcy fcz
> 1 1 1000 1000 1000
> #endif
>
> ; Include topology for ions
> #include "amber03.ff/ions.itp"
>
> #include "amber03.ff/dum.itp"
> #ifdef POSRES
> #include "posre_memb.itp"
> #endif
See new example
http://www.gromacs.org/Documentation/How-tos/Position_Restraints.
Whether posre_memb.itp makes sense depends whether there's an
appropriate [moleculetype] declared above it.
Now might be a good time to have made sure you've read and understood
the examples of 5.7.x of the manual...
Mark
>>> So you mean that it is indeed important if I call it POSRES oder
>>> DPOSRES,
>>> right?
>>>
>> Yes, syntax is crucial. The proper way to make use of #ifdef blocks is
>> described in the manual and in the link Mark posted.
>>
>> -Justin
>>
>> --
>> ========================================
>>
>> Justin A. Lemkul, Ph.D.
>> Research Scientist
>> 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 at gromacs.org
>> http://lists.gromacs.org/mailman/listinfo/gmx-users
>> * Only plain text messages are allowed!
>> * 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-request at gromacs.org.
>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>>
>
More information about the gromacs.org_gmx-users
mailing list