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) 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 to avoid its collapse.<br><br>I want to build biphastic system water-vacum-water where loops would be in water and TM helixes in vacum region.<br>
<br>Something like this I&#39;ve read in old publications about simulation of bacteriorhodopsin <a href="http://ukpmc.ac.uk/abstract/MED/10412722">http://ukpmc.ac.uk/abstract/MED/10412722</a><br><br><br>1- Could you tell me where I could found possible algorithm about builing of such water-vacum-water system?<br>
<br>2- Also I&#39;d like to specify what should I do with the charges residues. I&#39;d like to use AMBER-like or OPLS ff for such simulation As I understood I must neitralize only charges in TM helices and keep residues in LOOP intact. Might this aproache be correct?<br>
<br><br>James<br><br><div class="gmail_quote">2012/1/29 David van der Spoel <span dir="ltr">&lt;<a href="mailto:spoel@xray.bmc.uu.se">spoel@xray.bmc.uu.se</a>&gt;</span><br><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div class="im">On 2012-01-29 17:09, James Starlight wrote:<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Hi, Justin.<br>
<br>
Yes. The GFP chromophore is a part of backbone. It&#39;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&#39;ve made for this CRO residue topology via PRODG server for GROMOS ff.<br>
<br>
Than I&#39;ve imported that new chromophore.top into the topology.top of my<br>
structure in accordance to your tutorial.<br>
<br>
Finally I&#39;ve merged CRO.gro and protein.gro ( I&#39;ve cut CRO from the pdb<br>
for creation of the topoogy for my protein via pdb2gmx)<br>
<br>
Than I&#39;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>
</blockquote>
<br></div>
You need to add the bonds angles etc. The easiest way would be to make a special bond (specbond.dat file). You need an rtp entry for your cro group as well. Then pdb2gmx makes the necessary bonds.<br>
<br>
Of course you can make all the bonds, angles, diehdrals and pairs manually as well, but that is tedious and error prone.<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
James<br>
<br>
2012/1/29 Justin A. Lemkul &lt;<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> &lt;mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>&gt;&gt;<div><div class="h5"><br>
<br>
<br>
<br>
    James Starlight wrote:<br>
<br>
        Hi David!<br>
<br>
        Thanks for reference I&#39;ll study it carefully.<br>
<br>
<br>
        I have some general question about the vacuum simulation<br>
<br>
        1- I&#39;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&#39;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&#39;s GFP protein where<br>
        chromophore group ( like ligand) is covalent bonded to the<br>
        backbone of<br>
        this protein. As I&#39;ve understood in Justins tutorial there are<br>
        no any<br>
        covalent bonds between protein and ligand. But how I could make this<br>
        bond if I operate with TWO topology files ( one for chromophore and<br>
        another for protein itself) ? I&#39;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>
    ==============================<u></u>__==========<div class="im"><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></div>
    jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> &lt;<a href="http://vt.edu" target="_blank">http://vt.edu</a>&gt; | (540) 231-9080<br>
    <a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>__<a href="http://vt.edu/Pages/Personal/justin" target="_blank"><u></u>vt.edu/Pages/Personal/justin</a><br>
    &lt;<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>&gt;<br>
<br>
    ==============================<u></u>__==========<div class="im"><br>
<br>
    --<br>
    gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
    &lt;mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>&gt;<br>
    <a href="http://lists.gromacs.org/__mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/__<u></u>mailman/listinfo/gmx-users</a><div class="im"><br>
    &lt;<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a>&gt;<br>
    Please search the archive at<br></div>
    <a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a><div class="im"><br>
    &lt;<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a>&gt; before posting!<br>
    Please don&#39;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>
    &lt;mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>&gt;.<br>
    Can&#39;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>
    &lt;<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a>&gt;<br>
<br>
<br>
<br>
<br>
</blockquote><div class="im HOEnZb">
<br>
<br>
-- <br>
David van der Spoel, Ph.D., Professor of Biology<br>
Dept. of Cell &amp; 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></div><div class="HOEnZb"><div class="h5">
-- <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&#39;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&#39;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>