Hi, Everything should be updated correctly. I added in the DPPC 128 underneath the protein in [molecules] but I still get the same results. What do you mean by " Further manual modifications are necessary based on how many lipids InflateGRO removes"? After I generated the new position restrain file, I updated the minim.mdp withe the correct info. I then scaled by a factor of 4, updated [molecules], ran energy minimization (mdrun -v -deffnm em), and scaled by a factor 0.95. I have no idea what could possibly going wrong. Maybe something with the manipulation of InflateGRO????
Erica ________________________________ From: Justin A. Lemkul [via GROMACS] [ml-node+s5086n4998242...@n6.nabble.com] Sent: Friday, June 08, 2012 4:41 PM To: Hicks, Erica Subject: RE: Energy Minimization - not getting correct lipid area On 6/8/12 5:32 PM, Erica Hicks wrote: > Hi, > > I went back and found that a possible error could be that the coordinate file > was not updated correctly. After concatenated the Kalp_newbox.gro and > dppc128_whole.gro, it said to update the coordinate file with the correct > number of atoms. This should be the topology.top, right? But when I look > through this, all I find is: > [ molecules ] > ; Compound #mols > Protein 1 > Where should the number of atoms be and is topology.top the correct file? Make sure you added the correct topology information (step 2 in the tutorial, note the green text): http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/membrane_protein/02_topology.html Step 3.2 tells you: "Remove unnecessary lines (the box vectors from the KALP structure, the header information from the DPPC structure) and update the second line of the coordinate file (total number of atoms) accordingly." The tutorial assumes you're familiar with such operations, so in reality what it's implying is that, since you've now concatenated a protein and a 128-lipid membrane, you need to add: DPPC 128 to the [molecules] directive (after the Protein line to keep the correct order). Further manual modifications are necessary based on how many lipids InflateGRO removes. -Justin > > Erica > ________________________________ > From: Justin A. Lemkul [via GROMACS] [[hidden email]<UrlBlockedError.aspx>] > Sent: Friday, June 08, 2012 4:02 PM > To: Hicks, Erica > Subject: Re: Energy Minimization - not getting correct lipid area > > > > On 6/8/12 4:57 PM, Erica Hicks wrote: >> It was after just one iteration. The number of lipids removed after scaling >> the lipids by a factor of 4 was two. I am not sure how to update this as the >> [molecules] directive in the topology.top shows only the number of moles, >> i.e. 1. >> > > The number shown in [molecules] is the number of molecules of a given species. > If you started with 128 lipids, a removal of 4 means you need to list 124 > lipids. > > I don't know how you can get such a small area with only one iteration of > shrinking. Your box should be enormous at that point. > > -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 [hidden email]<UrlBlockedError.aspx> > 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 [hidden email]<UrlBlockedError.aspx>. > Can't post? Read http://www.gromacs.org/Support/Mailing_Lists > > > ________________________________ > If you reply to this email, your message will be added to the discussion > below: > > NAML<<UrlBlockedError.aspx>http://gromacs.5086.n6.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml> > > > -- > View this message in context: > http://gromacs.5086.n6.nabble.com/Energy-Minimization-not-getting-correct-lipid-area-tp4998235p4998241.html > Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- ======================================== 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 [hidden email]<UrlBlockedError.aspx> 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 [hidden email]<UrlBlockedError.aspx>. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists ________________________________ If you reply to this email, your message will be added to the discussion below: http://gromacs.5086.n6.nabble.com/Energy-Minimization-not-getting-correct-lipid-area-tp4998235p4998242.html To unsubscribe from Energy Minimization - not getting correct lipid area, click here<http://gromacs.5086.n6.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=4998235&code=ZXJpY2EuaGlja3NAbXkudW5kLmVkdXw0OTk4MjM1fC0xMTYzMzMyODk1>. NAML<http://gromacs.5086.n6.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml> -- View this message in context: http://gromacs.5086.n6.nabble.com/Energy-Minimization-not-getting-correct-lipid-area-tp4998235p4998243.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing list gmx-users@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-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
