Grompy looks like fancy. <br>Is it reliable for conducting GCMC?? <br><br>I also need to perform GCMC along with MD. <br><br>In Grompy, it seems that the Grompy is compatible with gromacs 4.0.5. <br>The version of gromacs I use is 4.5.X.<br>
Is it possible to use Grompy with my Gromacs??<br><br><br><div class="gmail_quote">2012/1/25 Gianluca Interlandi <span dir="ltr"><<a href="mailto:gianluca@u.washington.edu">gianluca@u.washington.edu</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Renè,<br>
<br>
This might be a silly question. In the documentation of GromPy it says that:<br>
<br>
"The system consists of mono-atomic water molceules defined in the MARTINI force field. The intermolecular interactions between such molecules are modelled using the Lennard-Jones potential only."<br>
<br>
My system is all-atom and all-hydrogen and I need to use the GBSA solvation model. Is it still going to work, or is the tutorial just an example of what you can do with GromPy?<br>
<br>
Thanks,<br>
<br>
Gianluca<div><div class="h5"><br>
<br>
On Tue, 24 Jan 2012, René Pool wrote:<br>
<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
By the way, the tutorial can be found on the wiki page at<br>
<a href="https://github.com/GromPy/GromPy/wiki/Tutorial:-Using-GromPy-in-GCMC-mode" target="_blank">https://github.com/GromPy/<u></u>GromPy/wiki/Tutorial:-Using-<u></u>GromPy-in-GCMC-mode</a>.<br>
<br>
==============================<u></u>=======================<br>
René Pool<br>
<br>
Division of Molecular and Computational Toxicology<br>
Department of Chemistry and Pharmaceutical Sciences<br>
Vrije Universiteit Amsterdam<br>
De Boelelaan 1083<br>
1081HV AMSTERDAM, the Netherlands<br>
-----<br>
IBIVU/Bioinformatics<br>
Department of Computer Science<br>
Vrije Universiteit Amsterdam<br>
De Boelelaan 1081a<br>
1081HV AMSTERDAM, the Netherlands<br>
<br>
Room P 2.75<br>
E: <a href="mailto:r.pool@vu.nl" target="_blank">r.pool@vu.nl</a><br>
T: +31 20 598 76 12<br>
F: +31 20 598 76 10<br>
==============================<u></u>=======================<br>
<br>
On 01/24/2012 09:30 AM, Pool, R. wrote:<br>
Hi Gianluca,<br>
<br>
You might want to take a look at GromPy<br>
(<a href="https://github.com/GromPy/GromPy" target="_blank">https://github.com/GromPy/<u></u>GromPy</a>).<br>
This is a python interface to the GROMACS library. Amongst other<br>
applications, there is an option to perform grand-canonical Monte Carlo<br>
using GromPy. The necessary energy evaluations are performed using direct<br>
library calls to GROMACS. In this way it is possible to get around the<br>
possibly prohibitive file I/O in shell implementations that use the<br>
GROMACS executables. A paper on GromPy and GromPy/GCMC has recently been<br>
accepted in JCC, so more details will follow. In the meantime you can<br>
have a look at the code and the tutorial on the above website.<br>
<br>
Good luck!<br>
<br>
Cheers,<br>
Rene<br>
==============================<u></u>=======================<br>
René Pool<br>
<br>
Division of Molecular and Computational Toxicology<br>
Department of Chemistry and Pharmaceutical Sciences<br>
Vrije Universiteit Amsterdam<br>
De Boelelaan 1083<br>
1081HV AMSTERDAM, the Netherlands<br>
-----<br>
IBIVU/Bioinformatics<br>
Department of Computer Science<br>
Vrije Universiteit Amsterdam<br>
De Boelelaan 1081a<br>
1081HV AMSTERDAM, the Netherlands<br>
<br>
Room P 2.75<br>
E: <a href="mailto:r.pool@vu.nl" target="_blank">r.pool@vu.nl</a><br>
T: +31 20 598 76 12<br>
F: +31 20 598 76 10<br>
==============================<u></u>=======================<br>
<br></div></div>
______________________________<u></u>______________________________<u></u>_________________________<div><div class="h5"><br>
From: <a href="mailto:gmx-users-bounces@gromacs.org" target="_blank">gmx-users-bounces@gromacs.org</a> [<a href="mailto:gmx-users-bounces@gromacs.org" target="_blank">gmx-users-bounces@gromacs.org</a><u></u>] on behalf<br>
of Rodrigo Faccioli [<a href="mailto:rodrigo_faccioli@uol.com.br" target="_blank">rodrigo_faccioli@uol.com.br</a>]<br>
Sent: 24 January 2012 05:29<br>
To: Discussion list for GROMACS users<br>
Subject: Re: [gmx-users] Monte Carlo with Gromacs and implicit solvent<br>
<br>
Hi Gianluca,<br>
<br>
I forgot to say that the step of compute the energies is considered from<br>
pdb2gmx until g_energy.<br>
<br>
Sorry about my oblivion.<br>
<br>
Best regards,<br>
<br>
--<br>
Rodrigo Antonio Faccioli<br>
Ph.D Student in Electrical Engineering<br>
University of Sao Paulo - USP<br>
Engineering School of Sao Carlos - EESC<br>
Department of Electrical Engineering - SEL<br>
Intelligent System in Structure Bioinformatics<br>
<a href="http://laips.sel.eesc.usp.br" target="_blank">http://laips.sel.eesc.usp.br</a><br>
Phone: 55 (16) 3373-9366 Ext 229<br>
Curriculum Lattes - <a href="http://lattes.cnpq.br/1025157978990218" target="_blank">http://lattes.cnpq.br/<u></u>1025157978990218</a><br>
Public Profile - <a href="http://br.linkedin.com/pub/rodrigo-faccioli/7/589/a5" target="_blank">http://br.linkedin.com/pub/<u></u>rodrigo-faccioli/7/589/a5</a><br>
<br>
<br>
On Tue, Jan 24, 2012 at 1:07 AM, Rodrigo Faccioli <<a href="mailto:rodrigo_faccioli@uol.com.br" target="_blank">rodrigo_faccioli@uol.com.br</a>><br>
wrote:<br>
Hi Gianluca,<br>
<br>
Thanks your question.<br>
<br>
All steps to obtain a generation with 200 individuals, on average,<br>
is 6 minutes for 1VII, example. In each generation I must <br>
converted my proteins Dihedral to Cartesian representation, compute<br>
the energies, choose the individuals which will go to reproduction,<br>
save previous population. This test was using one core of a<br>
conventional Desktop computer.<br>
<br>
The proteins, that I have tested my algorithm, vary from 20 to 56<br>
amino acids.<br>
<br>
Best regards,<br>
<br>
--<br>
Rodrigo Antonio Faccioli<br>
Ph.D Student in Electrical Engineering<br>
University of Sao Paulo - USP<br>
Engineering School of Sao Carlos - EESC<br>
Department of Electrical Engineering - SEL<br>
Intelligent System in Structure Bioinformatics<br>
<a href="http://laips.sel.eesc.usp.br" target="_blank">http://laips.sel.eesc.usp.br</a><br>
Phone: 55 (16) 3373-9366 Ext 229<br>
Curriculum Lattes - <a href="http://lattes.cnpq.br/1025157978990218" target="_blank">http://lattes.cnpq.br/<u></u>1025157978990218</a><br>
Public Profile -<br>
<a href="http://br.linkedin.com/pub/rodrigo-faccioli/7/589/a5" target="_blank">http://br.linkedin.com/pub/<u></u>rodrigo-faccioli/7/589/a5</a><br>
<br>
<br>
On Tue, Jan 24, 2012 at 12:36 AM, Gianluca Interlandi<br>
<<a href="mailto:gianluca@u.washington.edu" target="_blank">gianluca@u.washington.edu</a>> wrote:<br>
Thanks Enrico!<br>
<br>
Just wondering. How fast is it? Calling gromacs must have a<br>
lot of overhead. Also, do you call mdrun or does g_energy<br>
evaluate the energy?<br>
<br>
Gianluca<br>
<br>
On Tue, 24 Jan 2012, Rodrigo Faccioli wrote:<br>
<br>
Hi,<br>
<br>
Although I don't work with MC simulation, I have<br>
used Gromacs to obtain energies and others<br>
protein features in my Evolutionary Algorithms<br>
(EA).<br>
<br>
In general lines, I created a script that calls<br>
Gromacs programs and the output of these<br>
programs, such as g_energy (I read its xvg file),<br>
I stored it at a specific file which is<br>
read by my algorithm and its values are put into<br>
my data structure. I call Gromacs either<br>
implicit or explicit solvent. <br>
<br>
My project is working with EA until now. However,<br>
the integration with GROMACS is able to<br>
work with other kind of algorithms such as MC.<br>
Actually, we have idea to work with MC. But,<br>
I don't know when it will be possible. If you<br>
want, we can talk about join these work.<br>
<br>
Best regards,<br>
<br>
--<br>
Rodrigo Antonio Faccioli<br>
Ph.D Student in Electrical Engineering<br>
University of Sao Paulo - USP<br>
Engineering School of Sao Carlos - EESC<br>
Department of Electrical Engineering - SEL<br>
Intelligent System in Structure Bioinformatics<br>
<a href="http://laips.sel.eesc.usp.br" target="_blank">http://laips.sel.eesc.usp.br</a><br>
Phone: 55 (16) 3373-9366 Ext 229<br>
Curriculum Lattes -<br>
<a href="http://lattes.cnpq.br/1025157978990218" target="_blank">http://lattes.cnpq.br/<u></u>1025157978990218</a><br>
Public Profile -<br>
<a href="http://br.linkedin.com/pub/rodrigo-faccioli/7/589/a5" target="_blank">http://br.linkedin.com/pub/<u></u>rodrigo-faccioli/7/589/a5</a><br>
<br>
<br>
On Mon, Jan 23, 2012 at 11:39 PM, Gianluca<br>
Interlandi <<a href="mailto:gianluca@u.washington.edu" target="_blank">gianluca@u.washington.edu</a>> wrote:<br>
Hi!<br>
<br>
I would like to use gromacs to perform Monte<br>
Carlo simulations in implicit<br>
solvent of a protein near a surface. The<br>
protein is treated as a rigid body<br>
whereas the surface is fix.<br>
<br>
I see that there are plans to code MC into<br>
gromacs:<br>
<br>
<br>
<a href="http://www.gromacs.org/Developer_Zone/Programming_Guide/Monte_Carlo" target="_blank">http://www.gromacs.org/<u></u>Developer_Zone/Programming_<u></u>Guide/Monte_Carlo</a><br>
<br>
Is there a preliminary version?<br>
<br>
If not, I wonder whether anybody has tried<br>
to do MC with gromacs using IMD and<br>
MDDriver:<br>
<br>
<a href="http://www.baaden.ibpc.fr/projects/mddriver/" target="_blank">http://www.baaden.ibpc.fr/<u></u>projects/mddriver/</a><br>
<br>
Besides this. Is there an easy way to obtain<br>
the force field energy of a system<br>
using gromacs? Would I have to run a 0-steps<br>
MD and read out the energy? I know<br>
that this would have a big overhead in a MC<br>
simulation, but it might be worth<br>
trying.<br>
<br>
Thanks,<br>
<br>
Gianluca<br>
<br>
<br>
-----------------------------<u></u>------------------------<br>
Gianluca Interlandi, PhD<br>
<a href="mailto:gianluca@u.washington.edu" target="_blank">gianluca@u.washington.edu</a><br>
+1 (206) 685 4435<br>
<br>
<a href="http://artemide.bioeng.washington.edu/" target="_blank">http://artemide.bioeng.<u></u>washington.edu/</a><br>
<br>
Postdoc at the Department of Bioengineering<br>
at the University of Washington, Seattle WA<br>
U.S.A.<br>
<br>
-----------------------------<u></u>------------------------<br>
--<br>
gmx-users mailing list <br>
<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<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<br>
<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a><br>
before posting!<br>
Please don't post (un)subscribe requests to<br>
the list. Use the www interface or<br>
send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
<br>
<br>
<br>
<br>
<br>
------------------------------<u></u>-----------------------<br>
Gianluca Interlandi, PhD <a href="mailto:gianluca@u.washington.edu" target="_blank">gianluca@u.washington.edu</a><br>
+1 (206) 685 4435<br>
<a href="http://artemide.bioeng.washington.edu/" target="_blank">http://artemide.bioeng.<u></u>washington.edu/</a><br>
<br>
Postdoc at the Department of Bioengineering<br>
at the University of Washington, Seattle WA U.S.A.<br>
------------------------------<u></u>-----------------------<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/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at<br>
<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'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/<u></u>Support/Mailing_Lists</a><br>
<br>
<br>
<br>
<br>
<br>
</div></div></blockquote><div class="HOEnZb"><div class="h5">
<br>
------------------------------<u></u>-----------------------<br>
Gianluca Interlandi, PhD <a href="mailto:gianluca@u.washington.edu" target="_blank">gianluca@u.washington.edu</a><br>
+1 (206) 685 4435<br>
<a href="http://artemide.bioeng.washington.edu/" target="_blank">http://artemide.bioeng.<u></u>washington.edu/</a><br>
<br>
Postdoc at the Department of Bioengineering<br>
at the University of Washington, Seattle WA U.S.A.<br>
------------------------------<u></u>-----------------------</div></div><br>--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">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">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></blockquote></div><br>