Peter,<br><br>Yes the main reason is the CPU economy for such experiment. My current experiment consist of investigation of the tight paking
forces (primarily vdv effect) in the membrane receptor by inclusion of some point
mutations. <br><br><br>Actyally I think that simple biphastic system ( like in the Justin tutorial but with vacum layer instead of <font face="Arial"><font size="3">cyclohexane</font></font>) is exactly that I need. How I could make such layers system where water dont move into the middle layer?<br>
<br><br>James<br><br><br><br><br><br><div class="gmail_quote">2012/1/30 Peter C. Lai <span dir="ltr"><<a href="mailto:pcl@uab.edu" target="_blank">pcl@uab.edu</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I am not sure of the actual interpretive value of such a methodology.<br>
Are you just trying to save computational time by not having to simulate<br>
an all atom bilayer? The solvent layer is going to contribute to the<br>
majority of the computation cost in the first place. As an example, the<br>
bilayer we are using for GPCR work consisting of 238 POPC molecules adds<br>
30552 atoms out of a total of 99547 atoms, In any case, people have been<br>
conducting all-atom simulations of opsins since at least 2005 because<br>
of the availability of computational resources. See: Lemaitre V., Yeagle, P.,<br>
Watts, A. Biochemistry 2005, 44, 12667-12680 . So it is unlikely that<br>
any reviewers today will accept an opsin system simulated in anything but.<br>
<br>
Perhaps a biphasic implicit solvent model would work better for you (unsure<br>
if gromacs can support an "exterior" dielectric and an interior dielectric).<br>
There may also be some people using a coarse-grained bilayer with an all-atom<br>
protein, but again, unless you have carefully determined the coupling<br>
parameters between a MARTINI bilayer and an all-atom protein, there will<br>
be doubts about the usefulness of such a system.<br>
<br>
If you really must do without the bilayer, you might be able to get away<br>
with using strong i,i+4 distance constraints within your TMs to preserve<br>
helical stability while the entire protein is solvated.<br>
<div><div><br>
<br>
On 2012-01-30 10:50:35AM +0400, James Starlight wrote:<br>
> David, Justin<br>
><br>
> Thanks again for help!<br>
><br>
><br>
> I have just few questions about in vacum sumulation of membrane proteins.<br>
><br>
><br>
> I want to simulate GPCR receptor wich have big transmembrane ( helix)<br>
> domain and some connecting loop regions wich are exposed to the solvent.<br>
><br>
> As I understood in classical vacum simulation all charges must be redused<br>
> to avoid its collapse.<br>
><br>
> I want to build biphastic system water-vacum-water where loops would be in<br>
> water and TM helixes in vacum region.<br>
><br>
> Something like this I've read in old publications about simulation of<br>
> bacteriorhodopsin <a href="http://ukpmc.ac.uk/abstract/MED/10412722" target="_blank">http://ukpmc.ac.uk/abstract/MED/10412722</a><br>
><br>
><br>
> 1- Could you tell me where I could found possible algorithm about builing<br>
> of such water-vacum-water system?<br>
><br>
> 2- Also I'd like to specify what should I do with the charges residues. I'd<br>
> like to use AMBER-like or OPLS ff for such simulation As I understood I<br>
> must neitralize only charges in TM helices and keep residues in LOOP<br>
> intact. Might this aproache be correct?<br>
><br>
><br>
> James<br>
><br>
> 2012/1/29 David van der Spoel <<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a>><br>
><br>
> > On 2012-01-29 17:09, James Starlight wrote:<br>
> ><br>
> >> Hi, Justin.<br>
> >><br>
> >> Yes. The GFP chromophore is a part of backbone. It's formed from Ser Tyr<br>
> >> Gly by cyclisation of the Ser with Gly and subsequent dehydrotation. As<br>
> >> the consequence the mature chromophore has cyclised structure wich named<br>
> >> as the CRO residue in PDB structure.<br>
> >><br>
> >> I've made for this CRO residue topology via PRODG server for GROMOS ff.<br>
> >><br>
> >> Than I've imported that new chromophore.top into the topology.top of my<br>
> >> structure in accordance to your tutorial.<br>
> >><br>
> >> Finally I've merged CRO.gro and protein.gro ( I've cut CRO from the pdb<br>
> >> for creation of the topoogy for my protein via pdb2gmx)<br>
> >><br>
> >> Than I've done minimisation and chromophore have been diffused from my<br>
> >> protein :) It seems that I must add covalent bond between CRO and<br>
> >> protein into the topology. But how I could do it for my multi topology<br>
> >> file ?<br>
> >><br>
> ><br>
> > You need to add the bonds angles etc. The easiest way would be to make a<br>
> > special bond (specbond.dat file). You need an rtp entry for your cro group<br>
> > as well. Then pdb2gmx makes the necessary bonds.<br>
> ><br>
> > Of course you can make all the bonds, angles, diehdrals and pairs manually<br>
> > as well, but that is tedious and error prone.<br>
> ><br>
> >><br>
> >> James<br>
> >><br>
> >> 2012/1/29 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>>><br>
> >><br>
> >><br>
> >><br>
> >><br>
> >> James Starlight wrote:<br>
> >><br>
> >> Hi David!<br>
> >><br>
> >> Thanks for reference I'll study it carefully.<br>
> >><br>
> >><br>
> >> I have some general question about the vacuum simulation<br>
> >><br>
> >> 1- I've found that common GROMOS fields are not suitable for the<br>
> >> vacuum simulation because of its implementation for condensed<br>
> >> phase .<br>
> >> But In some referencces I've found that people use 53.6 ff for<br>
> >> the in<br>
> >> vacuum simulation. In addition Ive done minimisation and<br>
> >> equilibration<br>
> >> in that ff in vacuum and my system have not been collapse :) Is<br>
> >> there<br>
> >> any wy to adopt this ff for the in vacum ?<br>
> >><br>
> >> 2- I have uncommon onject for simulation. It's GFP protein where<br>
> >> chromophore group ( like ligand) is covalent bonded to the<br>
> >> backbone of<br>
> >> this protein. As I've understood in Justins tutorial there are<br>
> >> no any<br>
> >> covalent bonds between protein and ligand. But how I could make<br>
> >> this<br>
> >> bond if I operate with TWO topology files ( one for chromophore and<br>
> >> another for protein itself) ? I've done all steps in accordance to<br>
> >> Justins tutorial but on the minimisation step my chromphore dissuse<br>
> >> out of the protein interior because of absent of backbone group.<br>
> >><br>
> >><br>
> >> The GFP chromophore is part of the backbone of the protein, is it not?<br>
> >><br>
> >> The tutorial I have for a protein-ligand complex should not be taken<br>
> >> to mean that all non-protein entities are physically separate<br>
> >> entities. Plenty of cofactors, chromophores, and the like are<br>
> >> covalently attached to the protein.<br>
> >><br>
> >> -Justin<br>
> >><br>
> >> --<br>
</div></div>> >> ==============================**__==========<br>
<div>> >><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="http://vt.edu" target="_blank">http://vt.edu</a>> | (540) 231-9080<br>
</div>> >> <a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>__**<a href="http://vt.edu/Pages/Personal/justin" target="_blank">vt.edu/Pages/Personal/justin</a><<a href="http://vt.edu/Pages/Personal/justin" target="_blank">http://vt.edu/Pages/Personal/justin</a>><br>
> >> <<a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>**<a href="http://vt.edu/Pages/Personal/justin" target="_blank">vt.edu/Pages/Personal/justin</a><<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>
<div>> >><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>
</div>> >> <a href="http://lists.gromacs.org/__**mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/__**mailman/listinfo/gmx-users</a><<a href="http://lists.gromacs.org/__mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/__mailman/listinfo/gmx-users</a>><br>
> >><br>
> >> <<a href="http://lists.gromacs.org/**mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/**mailman/listinfo/gmx-users</a><<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a>><br>
<div>> >> ><br>
> >> Please search the archive at<br>
</div>> >> <a href="http://www.gromacs.org/__**Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__**Support/Mailing_Lists/Search</a><<a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__Support/Mailing_Lists/Search</a>><br>
> >><br>
> >> <<a href="http://www.gromacs.org/**Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/**Support/Mailing_Lists/Search</a><<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a>>><br>
<div>> >> before 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>
</div>> >> <mailto:<a href="mailto:gmx-users-request@" target="_blank">gmx-users-request@</a>**<a href="http://gromacs.org" target="_blank">gromacs.org</a><<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>><br>
> >> >.<br>
> >> Can't post? Read <a href="http://www.gromacs.org/__**Support/Mailing_Lists" target="_blank">http://www.gromacs.org/__**Support/Mailing_Lists</a><<a href="http://www.gromacs.org/__Support/Mailing_Lists" target="_blank">http://www.gromacs.org/__Support/Mailing_Lists</a>><br>
> >> <<a href="http://www.gromacs.org/**Support/Mailing_Lists" target="_blank">http://www.gromacs.org/**Support/Mailing_Lists</a><<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a>><br>
<div>> >> ><br>
> >><br>
> >><br>
> >><br>
> >><br>
> >><br>
> ><br>
> > --<br>
> > David van der Spoel, Ph.D., Professor of Biology<br>
> > Dept. of Cell & Molec. Biol., Uppsala University.<br>
> > Box 596, 75124 Uppsala, Sweden. Phone: +46184714205.<br>
> > <a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a> <a href="http://folding.bmc.uu.se" target="_blank">http://folding.bmc.uu.se</a><br>
> > --<br>
> > gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
</div>> > <a href="http://lists.gromacs.org/**mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/**mailman/listinfo/gmx-users</a><<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/**" target="_blank">http://www.gromacs.org/**</a><br>
> > Support/Mailing_Lists/Search<<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a>>before posting!<br>
<div>> > 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>
</div>> > Can't post? Read <a href="http://www.gromacs.org/**Support/Mailing_Lists" target="_blank">http://www.gromacs.org/**Support/Mailing_Lists</a><<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a>><br>
<div>> ><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/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
> Please don't post (un)subscribe requests to the list. Use the<br>
> 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/Support/Mailing_Lists</a><br>
<br>
<br>
</div><span><font color="#888888">--<br>
==================================================================<br>
Peter C. Lai | University of Alabama-Birmingham<br>
Programmer/Analyst | KAUL 752A<br>
Genetics, Div. of Research | 705 South 20th Street<br>
<a href="mailto:pcl@uab.edu" target="_blank">pcl@uab.edu</a> | Birmingham AL 35294-4461<br>
(205) 690-0808 |<br>
==================================================================<br>
</font></span><div><div><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/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
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/Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>