[gmx-users] Restrain specific lipid residues

[Katrina Lexa]

That's a simple enough way to deal with it - thanks!

On Sep 25, 2012, at 2:16 PM, Justin Lemkul wrote:

> 
> 
> On 9/25/12 5:14 PM, Katrina Lexa wrote:
>> Hi everyone,
>> 
>> This is probably a very silly question, but if I want to restrain only certain lipid residues in my bilayer, based on their residue number, is there some other way to do this aside from just having an explicit residue-by-residue list of them & their topologies in my .top? I'm just using Tieleman's lipid parameters, so I cannot define the lipid residues to restrain following the popc.itp file, since it describes just the general topology for any POPC residue & tells me that my index is out of range, since the .top is in atom numbers & I'm trying to think in residue numbers. Does that sort of make sense? Out of my 165 residues, I want to strongly restrain only 33 of them, and then let the rest minimize (I'm trying to add onto an equilibrated box - at this point, it would have just been easier to start over from scratch, but I thought I would ask).
>> 
> 
> Restraints can only be applied per [molecultype], so no, there is no way to have a global lipid topology and individual restraint files for individual molecules.  It is possible to reduce the complexity a bit by defining the 33 restrained lipids as one [moleculetype], though they would have to be consecutive in your coordinate file as well, which may require manual reorganization.
> 
>> Also, thank you to Justin and Mark for your help with  hydrogen building in pdb2gmx - you were right, obviously, I was just not noticing the wrong atom type in my .h2b file (way back Aug 30th).
>> 
> 
> Glad to hear there's at least one problem solved :)
> 
> -Justin
> 
> -- 
> ========================================
> 
> Justin A. Lemkul, Ph.D.
> Research Scientist
> Department of Biochemistry
> Virginia Tech
> Blacksburg, VA
> jalemkul[at]vt.edu | (540) 231-9080
> http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
> 
> ========================================
> -- 
> gmx-users mailing list    gmx-users at gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> * Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request at gromacs.org.
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> 
>