[gmx-users] DPPC replication error
Priya Das
priyadas001 at gmail.com
Tue Nov 25 06:31:20 CET 2014
Dear Sir,
I took 128 dppc.pdb, dppc.itp and lipid.itp from Tielemans site (
http://people.ucalgary.ca/~tieleman/download.html). For my work i need to
replicate 128 dppc.pdb four times to get 512 lipid molecules.
The sequence of commands are as follows:
# replicate 4 boxes of 128 dppc with SOL
$ genconf -f dppc128.gro -o dppcnew.gro -box 2 2 1
# generated .tpr file for new dppc with 512 lipid
$ grompp -f minim.mdp -c dppcnew.gro -p topol_dppc.top -o em.tpr
#removed periodicity
$ trjconv -s em.tpr -f dppcnew.gro -o dppcnew_whole.gro -pbc mol -ur compact
# inserted my protein in new lipid bilayer
$editconf -f 1AKL.gro -o 1AKL_new.gro -c -box 12.41840 12.44350 6.59650
#concatenate protein and lipid
cat 1AKL.gro dppcnew_whole.gro > system.gro
#inflate.pl
$ perl inflategro.pl system.gro 4 DPPC 0 system_inflated.gro 5 area.dat
*After running inflate.pl <http://inflate.pl>, the protein is centered
in the empty box ; with only 128 lipid. *
--
*Let us all join hands to save our " Mother Earth"*
------------------------------------------------------------------------
Regards,
Priya Das
Research Scholar
Dept. of Computational Biology and Bioinformatics,
University of Kerala
More information about the gromacs.org_gmx-users
mailing list