Dear Justin, <br> <br><br>Thanks.<br><br>The POPC bilayer i am using is with berger lipids, corrected for dihedrals so as to be compatible with the OPLS FF for aminoacids.<br><br>While searching for the literature on compatibility of lipid FF and protein FF, I found few references where similar modification was done for DOPC lipid bilayer and were suitable with various FF for proteins and also with CHARMM FF:<br>
<br>1. Membrane protein simulations with a united-atom lipid and all-atom protein model: lipid–protein interactions, side chain transfer free energies and model proteins.J. Phys.: Condens. Matter 18 (2006) S1221–S1234<br>
<br>2. Combination of the CHARMM27 Force Field with United-Atom Lipid Force Fields.J Comput Chem 32: 1400–1410, 2011.<br><br>I don't have the lipid bilayer with their itp files with CHARMM FF parameterization. Please could you inform me where to obtain them, so that i can use the lipid bilayer structure for embedding the protein and use the related CHARMM FF parameterised itp in the topology file in gromacs for MD simulation.<br>
<br>Thanks in advance,<br><br>Pramod<br><br><br><div class="gmail_quote">On Wed, Oct 12, 2011 at 7:51 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<br>
<br>
ram bio wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Justin,<br>
<br><div class="im">
Thanks, and I accept your suggestions;<br>
<br>
If SwissParam was designed to be used with CHARMM, the most intuitive next step is to use CHARMM for the MD, is it not?I understand the point about trying to keep the force fields consistent between docking and MD, but it may not be feasible (i.e., there may not be suitable parameters in OPLS for the bizarre functional groups you're dealing with).<br>
<br>
Yes, I also tried CHARMM FF to generate the topology file of the protein using pdb2gmx (without ligand), and as per the swissparam and gromacs tutorial i could build the protein-ligand-lipid bilayer and minimize it using mdrun and and i am at the NPT equilibration step, everything is ok with this procedure and without errors, but my lipid bilayer is made up of POPC and the POPC itp file has OPLS FF topologies. So, i was wondering whether the POPC itp file i am using for MD simulations can be used with the protein and ligand topology file generated by CHARMM.<br>
<br>
</div></blockquote>
<br>
You shouldn't mix and match force fields. Suitable CHARMM lipid parameters are widely available.<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
and as per the swissparam tutorial the command to generate topology file for protein is:<br>
<br>
pdb2gmx -f protein.pdb -ff charmm27 -water tip3p -ignh -o conf.pdb -nochargegrp<br>
<br>
<br>
<br>
in the gromacs 4.5.4 version the option to select Charmm FF from the pdb2gmx command is available, but i could not understand the usage of -nochargegp flag as per the tutorial, is this flag still valid while generating toplogies.<br>
<br>
</blockquote>
<br></div>
CHARMM does not use charge groups. Therefore, each atom should be its own "group" in the topology. Using -nochargegrp overrides the default behavior of the .rtp files (which has multi-atom charge groups, although I think this was changed somewhere along the way, but I don't remember if it was before or after 4.5.4).<br>
<font color="#888888">
<br>
-Justin</font><div><div></div><div class="h5"><br>
<br>
-- <br>
==============================<u></u>==========<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | <a href="tel:%28540%29%20231-9080" value="+15402319080" target="_blank">(540) 231-9080</a><br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a><br>
<br>
==============================<u></u>==========<br>
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>