VMD-L Mailing List
From: Arthur Vale (arthurovale_at_uchicago.edu)
Date: Thu Jul 20 2017 - 15:19:54 CDT
- Next message: Vermaas, Joshua: "Re: Building HMMM using patch files"
- Previous message: Giacomo Fiorin: "Re:"
- Next in thread: Vermaas, Joshua: "Re: Building HMMM using patch files"
- Reply: Vermaas, Joshua: "Re: Building HMMM using patch files"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]
Dear All,
Me and my group are attempting to turn a regular lipid membrane into an HMMM as described in this article: doi: 10.1021/acs.jpcb.6b11378. We have previously built HMMM membranes using CHARMM-GUI but we want to use the enhanced version of it described in the aforementioned paper. We have a pre-equilibrated regular lipid membrane generated through CHARMM-GUI lipid generator, as well as the files provided in the article. From reading the topology files we are under the impression that we can use the patches topology file to both cut the lipid tails and generate the SCSE solvent. Since previously we were successful doing the standard procedure with the autoPSF in VMD, we don't know how to actually proceed to apply these patches in order to achieve this. We tried doing it with the add patches in autoPSF as well as the patch command in psfgen but we were unsuccessful. Could anyone describe the steps to achieve this? I attached the contents of the topology and parameter files provided in the article below for reference.
Thank you in advance,
Arthur Vale
==> patches.top <==
!Patch to terminate C2 tail at carbon 5.
PRES C25T
ATOM C25 CTL3 -0.27
ATOM H5R HAL3 0.09
ATOM H5S HAL3 0.09
ATOM H5T HAL3 0.09
BOND C25 H5T
!Patch to terminate C3 tail at carbon 5.
PRES C35T
ATOM C35 CTL3 -0.27
ATOM H5X HAL3 0.09
ATOM H5Y HAL3 0.09
ATOM H5Z HAL3 0.09
BOND C35 H5Z
==> top_scs.top <==
MASS 1 SCSM 14.027000 SC ! "Simple Carbon Solvent Methane"
MASS 2 SCSE 14.027000 SC ! "Simple Carbon Solvent Ethane"
DEFA FIRS none LAST none
AUTOGENERATE ANGLES DIHEDRALS
RESI SCSM 0.00
ATOM C SCSM 0.00
RESI SCSE !Simple Carbon Solvent ethane
ATOM C1 SCSE 0.00
ATOM C2 SCSE 0.00
BOND C1 C2
END
==> par_scs.prm <==
BONDS
SCSE SCSE 222.500 1.5300 ! CTL2-CTL2 bond parameter.
NONBONDED nbxmod 5 atom cdiel shift vatom vdistance vswitch -
cutnb 14.0 ctofnb 12.0 ctonnb 10.0 eps 1.0 e14fac 1.0 wmin 1.5
!adm jr., 5/08/91, suggested cutoff scheme
!
!V(Lennard-Jones) = Eps,i,j[(Rmin,i,j/ri,j)**12 - 2(Rmin,i,j/ri,j)**6]
!
!epsilon: kcal/mole, Eps,i,j = sqrt(eps,i * eps,j)
!Rmin/2: A, Rmin,i,j = Rmin/2,i + Rmin/2,j
!
!atom ignored epsilon Rmin/2 ignored eps,1-4 Rmin/2,1-4
!
SCSM 0.0 -0.350 2.080
SCSE 0.0 -0.112 2.080
NBFIX
SCSM SCSM -0.625 3.15
SCSE SCSE -0.274 3.15
END
- Next message: Vermaas, Joshua: "Re: Building HMMM using patch files"
- Previous message: Giacomo Fiorin: "Re:"
- Next in thread: Vermaas, Joshua: "Re: Building HMMM using patch files"
- Reply: Vermaas, Joshua: "Re: Building HMMM using patch files"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]