[gmx-users] Replica Exchange for surface tension

gozde ergin gozdeeergin at gmail.com
Wed Mar 29 17:29:05 CEST 2017


Justin what would you suggest if I simulate the system under constant surface tension since I know the correct surface tension from experiments.
Can I get some information about the force field quality if I compare the conformation of the surface under constant surface tension (N\gammaT) and under NVT ensemble?
 
> On 29 Mar 2017, at 15:21, gozde ergin <gozdeeergin at gmail.com> wrote:
> 
> Exactly SDS is CHARMM36 but cocobetaine not.
> 
>> On 29 Mar 2017, at 15:18, Justin Lemkul <jalemkul at vt.edu <mailto:jalemkul at vt.edu>> wrote:
>> 
>> 
>> 
>> On 3/29/17 9:16 AM, gozde ergin wrote:
>>> Are you sure it is not CHARMM36 because it is written:
>>> 
>>> ;;
>>> ;; Generated by CHARMM-GUI (http://www.charmm-gui.org <http://www.charmm-gui.org/>) v1.7
>>> ;;
>>> ;; psf2itp.py
>>> ;;
>>> ;; Correspondance:
>>> ;; j712l362 at ku.edu <mailto:j712l362 at ku.edu> or wonpil at ku.edu <mailto:wonpil at ku.edu>
>>> ;;
>>> ;; CHARMM36 FF in GROMACS format
>>> ;;
>>> 
>>> on the top of the .itp files.
>>> 
>> 
>> OK, I was thinking too quickly, the workflow goes like this: If a molecule is identified as part of the existing force field, it is used.  If it is not, then CGenFF is called.  So your SDS may be CHARMM36 but your cocobetaine is not. You can check against the CHARMM force field files and also the atom types themselves.  If the second character is "G" then it is CGenFF, otherwise it is stock CHARMM.
>> 
>> Do not rely on header comments; those are often boilerplate.  Check the parameters themselves.
>> 
>> -Justin
>> 
>>> 
>>> 
>>>> On 29 Mar 2017, at 15:06, Justin Lemkul <jalemkul at vt.edu <mailto:jalemkul at vt.edu>> wrote:
>>>> 
>>>> 
>>>> 
>>>> On 3/29/17 4:15 AM, gozde ergin wrote:
>>>>> I used CHARMM GUI Input Generator tool both for SDS and cocobetaine.
>>>> 
>>>> For the sake of clarity, you should understand that you are *not* then using CHARMM36.  CHARMM-GUI calls the CGenFF program to parametrize the molecule, so you are using CGenFF parameters, which are compatible with CHARMM by design, but are not the standard CHARMM36 parameters.  Thus your SDS model will differ from the standard one that we have in the lipids/detergents file.  You should compare it to make sure it is reasonable.  A general force field will sacrifice some accuracy in the name of generality and transferability.
>>>> 
>>>>> So basically the reason of not getting correct surface tension is bad parametrisation of force field.
>>>>> Then I should not spend time on simulation the replica exchange.
>>>> 
>>>> Indeed.  Look at the things I mentioned before.  If the model doesn't work well, then you have to get into the details of why.  Burning through more time hoping it will work out and expensive stuff like REMD is just a waste.
>>>> 
>>>> -Justin
>>>> 
>>>> --
>>>> ==================================================
>>>> 
>>>> Justin A. Lemkul, Ph.D.
>>>> Ruth L. Kirschstein NRSA Postdoctoral Fellow
>>>> 
>>>> Department of Pharmaceutical Sciences
>>>> School of Pharmacy
>>>> Health Sciences Facility II, Room 629
>>>> University of Maryland, Baltimore
>>>> 20 Penn St.
>>>> Baltimore, MD 21201
>>>> 
>>>> jalemkul at outerbanks.umaryland.edu <mailto:jalemkul at outerbanks.umaryland.edu> | (410) 706-7441
>>>> http://mackerell.umaryland.edu/~jalemkul <http://mackerell.umaryland.edu/~jalemkul>
>>>> 
>>>> ==================================================
>>>> --
>>>> Gromacs Users mailing list
>>>> 
>>>> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
>>>> 
>>>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>>>> 
>>>> * For (un)subscribe requests visit
>>>> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.
>>> 
>> 
>> -- 
>> ==================================================
>> 
>> Justin A. Lemkul, Ph.D.
>> Ruth L. Kirschstein NRSA Postdoctoral Fellow
>> 
>> Department of Pharmaceutical Sciences
>> School of Pharmacy
>> Health Sciences Facility II, Room 629
>> University of Maryland, Baltimore
>> 20 Penn St.
>> Baltimore, MD 21201
>> 
>> jalemkul at outerbanks.umaryland.edu <mailto:jalemkul at outerbanks.umaryland.edu> | (410) 706-7441
>> http://mackerell.umaryland.edu/~jalemkul <http://mackerell.umaryland.edu/~jalemkul>
>> 
>> ==================================================
>> -- 
>> Gromacs Users mailing list
>> 
>> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List <http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List> before posting!
>> 
>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists <http://www.gromacs.org/Support/Mailing_Lists>
>> 
>> * For (un)subscribe requests visit
>> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users <https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users> or send a mail to gmx-users-request at gromacs.org <mailto:gmx-users-request at gromacs.org>.
> 



More information about the gromacs.org_gmx-users mailing list