AW: AW: [gmx-users] question about the minimization

Rausch, Felix frausch at ipb-halle.de
Thu Oct 18 16:07:23 CEST 2012


If you got the same problems with the protein only (and without any restraints), I wonder if you have any non-standard residues (i.e. ligands or protein modifications) present in your system?
Because a normal protein should work just fine with the parameters you gave us.

-----Ursprüngliche Nachricht-----
Von: gmx-users-bounces at gromacs.org [mailto:gmx-users-bounces at gromacs.org] Im Auftrag von reisingere at rostlab.informatik.tu-muenchen.de
Gesendet: Donnerstag, 18. Oktober 2012 15:28
An: Discussion list for GROMACS users
Betreff: Re: AW: [gmx-users] question about the minimization

Hi Justin,
I am sorry, I switched the radius and the numbers:
radius=0
charge=1
mass=1

But I already tried it without the membrane and there I had the same problem. And also when I did not restrained anything. Not even the membrane I got the high force and the LINCS warnings.

Of what physical problems do you think?

I already looked at the protein with Pymol and the residue with the high force does not in interact with any other.

>
>
> On 10/18/12 8:51 AM, reisingere at rostlab.informatik.tu-muenchen.de wrote:
>> Hi Justin,
>>
>> Well I already tried the protein without the membrane and I got the 
>> same result. And I also make sure that I had no error messages during 
>> the preparation runs.
>>
>
> If the protein crashes in the absence of anything else, then likely 
> you have problems with its starting geometry.  Applying restraints 
> during minimization is probably causing the clashes to be 
> insufficiently resolved.
>
>> I am using the amber03 force field.
>> My dummy atoms have a radius of 1.00 and a charge of 0.00. And yes 
>> they are non bonded atoms.
>>
>
> Radius of 1.0 nm?  If the crash occurs without these atoms, then it's 
> likely irrelevant, but with particles 10 times the size of normal 
> atoms, you may be asking for trouble.
>
> -Justin
>
>
>>> On 10/18/12 7:05 AM, reisingere at rostlab.informatik.tu-muenchen.de
>>> wrote:
>>>> Hi Justin,
>>>>
>>>> I already read the page you send me.
>>>> And I am sure about the first points. But I am not completely sure 
>>>> about
>>>
>>> Sure about the first five points?  The five suggestions on the page 
>>> I linked are all diagnostic steps.  What did they turn up?
>>>
>>>> the mdp file for the nvt run.
>>>>
>>>> It looks like this:
>>>>
>>>> define               = -DPOSRES
>>>> integrator           = md
>>>> dt                   = 0.002
>>>> nsteps               = 10000
>>>> nstxout              = 0
>>>> nstvout              = 0
>>>> nstfout              = 0
>>>> nstlog               = 1000
>>>> nstxtcout            = 0
>>>> 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
>>>>
>>>
>>> What force field are you using?  What about my questions related to 
>>> your so-called dummy atoms?  What are they and what are their 
>>> parameters?
>>> Your
>>> minimization seemed fine, so there's clearly something wrong with 
>>> the physical model that leads to a crash during MD.  At present, 
>>> there's not enough information to say exactly what.
>>>
>>> -Justin
>>>
>>>>
>>>> Thank you for your help,
>>>>    Eva
>>>>
>>>>
>>>>>
>>>>>
>>>>> On 10/18/12 6:15 AM, reisingere at rostlab.informatik.tu-muenchen.de
>>>>> wrote:
>>>>>> Hi Erik,
>>>>>>
>>>>>> my .mdp file looks like this:
>>>>>>
>>>>>> define                  = -DPOSRES
>>>>>> integrator              = steep
>>>>>> emtol                   = 10
>>>>>> nsteps                  = 30000
>>>>>> 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 yes, I use constraints. I have a membrane represented by 
>>>>>> dummy atoms around my protein. They are fixed and I also fixed 
>>>>>> the backbone. Or is that a problem.
>>>>>>
>>>>>
>>>>> The EM seems to have worked fine, so we need the .mdp file for 
>>>>> NVT, not EM.
>>>>> Your membrane representation could be an issue, but I'm not clear 
>>>>> on the terminology here.  A "dummy atom" (or virtual site) is a 
>>>>> particle with nonbonded parameters but no mass.  Is that what 
>>>>> you're dealing with?  If this membrane environment is fixed and 
>>>>> the backbone is too, you may get clashes in the system, but it 
>>>>> doesn't sound like there's a whole lot that's moving in the first 
>>>>> place.
>>>>>     General troubleshooting advice (that I seem to post almost
>>>>> daily):
>>>>>
>>>>> http://www.gromacs.org/Documentation/Terminology/Blowing_Up#Diagno
>>>>> sing_an_Unstable_System
>>>>>
>>>>> Your simulation is failing fairly early, so finding the source of 
>>>>> the crash should not be terribly difficult.
>>>>>
>>>>> -Justin
>>>>>
>>>>>
>>>>>>> It could be due to your mdp-settings. What's your time step? Do 
>>>>>>> you use constraints and virtual interaction sites?
>>>>>>>
>>>>>>> Erik
>>>>>>>
>>>>>>> 18 okt 2012 kl. 12.04 skrev
>>>>>>> reisingere at rostlab.informatik.tu-muenchen.de:
>>>>>>>
>>>>>>>> Hi Felix,
>>>>>>>> thank you for your answer.
>>>>>>>> So you think that a force of 3.2341940e+02 is okey? To work on?
>>>>>>>>
>>>>>>>> I also tried to work on with a NVT run but then I got many many 
>>>>>>>> LINCS
>>>>>>>> warnings:
>>>>>>>>
>>>>>>>> Step 1114, time 2.228 (ps)  LINCS WARNING relative constraint 
>>>>>>>> deviation after LINCS:
>>>>>>>> rms 0.000169, max 0.001796 (between atoms 3480 and 3482) bonds 
>>>>>>>> that rotated more than 30 degrees:
>>>>>>>> atom 1 atom 2  angle  previous, current, constraint length
>>>>>>>>      3480   3482   30.3    0.1090   0.1092      0.1090
>>>>>>>>
>>>>>>>> Step 1121, time 2.242 (ps)  LINCS WARNING relative constraint 
>>>>>>>> deviation after LINCS:
>>>>>>>> rms 0.000108, max 0.000561 (between atoms 2000 and 2001) bonds 
>>>>>>>> that rotated more than 30 degrees:
>>>>>>>> atom 1 atom 2  angle  previous, current, constraint length
>>>>>>>>      2002   2005   30.4    0.1090   0.1090      0.1090
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> And in the end there is the error:
>>>>>>>>
>>>>>>>> Fatal error:
>>>>>>>> 11 particles communicated to PME node 2 are more than 2/3 times 
>>>>>>>> the cut-off out of the domain decomposit ion cell of their 
>>>>>>>> charge group in dimension x.
>>>>>>>>
>>>>>>>>
>>>>>>>> I looked it up on the GROMACS site and it was written there 
>>>>>>>> that such a error is caused by a not well equilibrated system.
>>>>>>>>
>>>>>>>> So I thought that the high force might be the reason.
>>>>>>>>
>>>>>>>> What do you think?
>>>>>>>>
>>>>>>>> Thanks again!!
>>>>>>>>
>>>>>>>>
>>>>>>>>> Hi Eva,
>>>>>>>>>
>>>>>>>>> The change in energy in every step became smaller than the 
>>>>>>>>> machine precision can represent (as is stated in the output). 
>>>>>>>>> In other
>>>>>>>>> words:
>>>>>>>>> your
>>>>>>>>> system is minimized to a very high degree.
>>>>>>>>> A value of 10 for the Fmax is very low for a system like 
>>>>>>>>> yours, even with a pure protein in water system you probably 
>>>>>>>>> wouldn't hit this value.
>>>>>>>>> Furthermore, 30.000 steps for a minimization seems pretty high 
>>>>>>>>> to me.
>>>>>>>>>
>>>>>>>>> Cheers,
>>>>>>>>> Felix
>>>>>>>>>
>>>>>>>>> -----Ursprüngliche Nachricht-----
>>>>>>>>> Von: gmx-users-bounces at gromacs.org 
>>>>>>>>> [mailto:gmx-users-bounces at gromacs.org]
>>>>>>>>> Im Auftrag von reisingere at rostlab.informatik.tu-muenchen.de
>>>>>>>>> Gesendet: Donnerstag, 18. Oktober 2012 10:56
>>>>>>>>> An: gmx-users at gromacs.org
>>>>>>>>> Betreff: [gmx-users] question about the minimization
>>>>>>>>>
>>>>>>>>> Hi everybody,
>>>>>>>>>
>>>>>>>>> I have a question about the minimization of my protein. For 
>>>>>>>>> this minimization I fixed the backbone and the membrane atoms 
>>>>>>>>> which around my protein.
>>>>>>>>>
>>>>>>>>> I minimized it 30000 steps long but after step 20385 it stops 
>>>>>>>>> with the
>>>>>>>>> output:
>>>>>>>>>
>>>>>>>>> Steepest Descents converged to machine precision in 20385 
>>>>>>>>> steps, but did not reach the requested Fmax < 10.
>>>>>>>>> Potential Energy  = -8.8805600e+05
>>>>>>>>> Maximum force     =  3.2341940e+02 on atom 3050
>>>>>>>>> Norm of force     =  2.0930219e+00
>>>>>>>>>
>>>>>>>>> I do not see why it can not reach a lower force.
>>>>>>>>>
>>>>>>>>> I already looked at the structure with Pymol and there I saw 
>>>>>>>>> that the atom
>>>>>>>>> 3050 is in a loop and do not interact with any other residue.
>>>>>>>>>
>>>>>>>>> My .mdp file looks like this:
>>>>>>>>>
>>>>>>>>> define                  = -DPOSRES
>>>>>>>>> integrator              = steep
>>>>>>>>> emtol                   = 10
>>>>>>>>> nsteps                  = 30000
>>>>>>>>> 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
>>>>>>>>>
>>>>>>>>> Do you see a reason why it can not reach a force lower than 10?
>>>>>>>>>
>>>>>>>>> Thank you,
>>>>>>>>> Eva
>>>>>>>>>
>>>>>>>>> --
>>>>>>>>> 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/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
>>>>>>>>> --
>>>>>>>>> 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/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
>>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> --
>>>>>>>> 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/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
>>>>>>>
>>>>>>> -----------------------------------------------
>>>>>>> Erik Marklund, PhD
>>>>>>> Dept. of Cell and Molecular Biology, Uppsala University.
>>>>>>> Husargatan 3, Box 596,    75124 Uppsala, Sweden
>>>>>>> phone:    +46 18 471 6688        fax: +46 18 511 755
>>>>>>> erikm at xray.bmc.uu.se
>>>>>>> http://www2.icm.uu.se/molbio/elflab/index.html
>>>>>>>
>>>>>>> --
>>>>>>> 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/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
>>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>> --
>>>>> ========================================
>>>>>
>>>>> 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
>>>>> * 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
>>>>>
>>>>
>>>>
>>>
>>> --
>>> ========================================
>>>
>>> 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
>>> * 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
>>>
>>
>>
>
> --
> ========================================
>
> 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
> * 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
>


-- 
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/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