Well yes I guess I am getting the same thing what I expected but now I am little confused here and I don't know how to explain you this but I will try. <div>Initially I aligned my protein using TMDET results and centered that at ORIGIN and then I centered the lipid at origin using geom_center script provided in VMD and when I visualized that, the protein was properly inserted in the lipid the way MARTINI folks have shown in their website.</div>
<div><br></div><div>After that I followed your tutorial and made the system.gro file by catting lipid and protein and its visualization gives me the same thing as I was getting earlier. Now after that when I run the inflategro script for the system.gro file and scale the lipid by the factor of 4, the protein goes apart from lipid (I dnt know why) and after our recent conversation I did the same thing but now the protein is shifted little upwards but still it is at the center of lipid. Now I don't know whether this factor will matter a lot or not? I hope you got my question.</div>
<div><br></div><div>In summary, earlier the upper portion of protein was inserted fully same as shown in the martini website but after doing the inflategro without running EM it shifted little upwards but still at the center so I dnt know whether this will matter a lot or not?<br>
<div><br></div><div><br><br><div class="gmail_quote">On Thu, Oct 29, 2009 at 5:26 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;">
<div class="im"><br>
<br>
sunny mishra wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Alright. I have done this....I have made another lipid_newbox.gro file and wrote same box vectors as of lipid.gro file<br>
<br>
For LIPID :<br>
<br>
editconf -f lipid.gro -o lipid_newbox.gro -c -box 13.29820 13.29820 6.59160<br>
<br>
For Protein :<br>
<br>
editconf -f protein.gro -o protein_newbox.gro -c -box 13.29820 13.29820 6.59160<br>
<br>
Then cat the system and made system.gro file and did inflategro of the system.gro file. Do you think this is the correct way and now when I visualize this in VMD the protein is shifted little upwards but it is at the center of the lipid.<br>
<br>
</blockquote>
<br></div>
Well, did you get what you expected when you ran InflateGRO? The protein and lipid should have the same center; that's what the editconf commands did, as I've told you. If you need the protein to be placed differently, you must specify a different center (providing the coordinates to editconf -center).<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Sunny<div class="im"><br>
<br>
On Thu, Oct 29, 2009 at 4:59 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
<br>
<br>
<br>
sunny mishra wrote:<br>
<br>
No. But in order to do that with lipid.gro what co-ordinates<br>
should I have to take?<br>
<br>
<br>
I think you need to differentiate between "coordinates" and "box<br>
vectors." You should center the molecule in the box (defined by the<br>
vectors in the last line of the .gro file). I assumed you were<br>
specifying these box vectors with -box.<br>
<br>
So, for both the protein and the lipid coordinate files, you need to use<br>
<br>
editconf -c -box (x, y, z from .gro file)<br>
<br>
-Justin<br>
<br>
On Thu, Oct 29, 2009 at 4:40 PM, Justin A. Lemkul<br>
<<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>><br></div><div><div></div><div class="h5">
<mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>>> wrote:<br>
<br>
<br>
<br>
sunny mishra wrote:<br>
<br>
okk...here are the commands which I am givng...<br>
<br>
editconf -f 1SU4_cleancg_trans.gro -o 1SU4_newbox.gro -c -box<br>
13.29820 13.29820 6.59160<br>
<br>
I get the 1SU4_newbox.gro which has co-ordinates 13.29820<br>
13.29820 6.59160<br>
<br>
<br>
Have you done the same with lipid.gro? If the protein and<br>
the lipid<br>
are in different coordinate systems (i.e., lipid.gro centered<br>
at the<br>
origin and protein centered within the box, which has its corner<br>
placed at the origin) then you will have a problem.<br>
<br>
-Justin<br>
<br>
cat 1SU4_newbox.gro lipid.gro > system..gro<br>
<br>
Now here system.gro also has co-ordinates 13.29820<br>
13.29820 6.59160<br>
<br>
now inflategro:<br>
<br>
inflategro system.gro 4 DSPC 14 -o inflated_bilayer.gro 5<br>
<br>
Now when i visualize my output file in VMD then protein and<br>
lipid are seperated and even after i scale it to .95 they dnt<br>
meet....they are still apart, I was hoping that lipid will be<br>
scaled and protein shud have remained in the center of<br>
the lipid<br>
but that doesn't happen. I hope you got my question<br>
<br>
On Thu, Oct 29, 2009 at 4:25 PM, Justin A. Lemkul<br>
<<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>><br>
<mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>><br>
<mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>><br>
<mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>>>> wrote:<br>
<br>
<br>
<br>
sunny mishra wrote:<br>
<br>
Yes I did the same and when I see my system.gro<br>
file it<br>
gives me<br>
the dimensions of the protein box half of the lipid<br>
bilayer box<br>
so that means it should now be set in the center<br>
of the lipid<br>
bilayer. But after<br>
<br>
<br>
I don't understand what you mean.<br>
<br>
<br>
that when I run the inflategro script and see my<br>
output<br>
file the<br>
protein and lipid are separated and I dnt know why?<br>
<br>
<br>
Can you post the actual commands you're using?<br>
<br>
-Justin<br>
<br>
<br>
-- ========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<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="http://vt.edu" target="_blank">http://vt.edu</a>> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>><br>
<<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (540)<br>
<br>
231-9080<br>
<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a> <mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>>><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>> <mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>>>><br>
<br>
<br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list.<br>
Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>>>.<br>
<br>
Can't post? Read<br>
<a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br>
<br>
<br>
-- ========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<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="http://vt.edu" target="_blank">http://vt.edu</a>> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (540)<br>
231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a> <mailto:<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/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br>
<br>
<br>
-- ========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<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="http://vt.edu" target="_blank">http://vt.edu</a>> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<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/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br>
<br>
</div></div></blockquote><div><div></div><div class="h5">
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<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/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br></div></div>