Re: [AMBER] sum of charges from parm file

From: Pengfei Li <ambermailpengfei.gmail.com>
Date: Thu, 9 Apr 2020 12:12:11 -0400

Hi Lenka,

I guess Clorice implied that 0.15 M NaCl solution would not retain this concentration after you running the NPT simulations. This is because when adding the water box tleap usually generates a system with density smaller than 1.0 g/ml. So if you want to have 0.15 M concentration of NaCl in your production MD simulations (which is usually done after NPT equilibration that corrects the system density), you need to recalculate how many Na+ and Cl- ions you needed based on the volume of the system after the NPT equilibration, and then start your whole simulations from the tleap modeling all over again, in this way you would likely to have the ion concentration as you want in your production MD simulations.

Hope it helps,
Pengfei

> On Mar 31, 2020, at 6:47 PM, Lenka <Vellryba.seznam.cz> wrote:
>
>
> Dear Clorice,
>
>
>
>
>
> Of course! Duh! Not my best day:)
>
>
>
>
>
> Thank you for answering my question, yes that all makes sense now and I dont
> have to dive into Ewald calculations.
>
>
>
>
>
> L.
>
>
>
>
>
>
>
> ---------- Původní e-mail ----------
> Od: Clorice Reinhardt <clorice.reinhardt.yale.edu>
> Komu: AMBER Mailing List <amber.ambermd.org>
> Datum: 31. 3. 2020 23:42:06
> Předmět: Re: [AMBER] sum of charges from parm file
> "Hi Lenka,
>
> There are ions in the selected file yes, but because its NaCl, there are
> equal amounts Na+ as Cl-, so it will be neutral. I think the answer to your
> first question will resolve the second, but i am unsure whether you are
> also referring to building the protein in tleap or something with the
> actual MD.
>
> Hopefully that helped.
> Clorice
>
> On Tue, Mar 31, 2020 at 6:28 PM Lenka <Vellryba.seznam.cz> wrote:
>
>>
>> Hi,
>>
>>
>>
>>
>>
>> I am to simulate protein in 0.15M NaCl. It is a membrane protein so I
>> build
>> the membrane in Charmm and generate files for Amber.
>>
>>
>>
>>
>>
>> Despite the fact I select 0.15M NaCl, I get this message in the min_mdout
>> file:
>>
>>
>>
>>
>>
>> Sum of charges from parm topology file = 0.00000000
>> Forcing neutrality...
>>
>>
>>
>>
>>
>> Shouldn't there be a charge present in the parm file if the ions add to
> 150
>> mM NaCl? There are ions present in the generated pdb.
>>
>>
>>
>>
>>
>> My second question is, if I select the charge will the Amber PME
>> calculation
>> neutralize the charge anyway? I was under the impression that this happens
>
>> so how can you run a simulation at 0.15M?
>>
>>
>>
>>
>>
>> Thank you for the explanation
>>
>>
>>
>>
>>
>> L.
>> _______________________________________________
>> AMBER mailing list
>> AMBER.ambermd.org
>> http://lists.ambermd.org/mailman/listinfo/amber
>>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
> "
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Apr 09 2020 - 09:30:01 PDT
Custom Search