Liu, Liang wrote:
This is a "reasonable" answer :)
Thanks for that, and I just tried Gromacs for minimization, and looks
the final structure does not have clashes anymore, and also are very
close to the initial structure.
Another question is if there is a way to add ions automatically, meaning
no need to check the "NOTE" of "System has non-zero total charge" in the
output of grompp command? And also update the topology file
automatically too?
That depends on the system. For a simple solute in water, the running total
charge of the solute is recorded in the qtot column of the .top file. If you
have a complex system with lots of charged things, no, there is no convenient
way to get the charge aside from calculating it by hand or reading the grompp
output.
The topology can be modified automatically by genion after ions are added by
using the -p flag.
-Justin
On Mon, Nov 14, 2011 at 3:49 PM, Justin A. Lemkul <jalem...@vt.edu
<mailto:jalem...@vt.edu>> wrote:
Liu, Liang wrote:
Hi, all,
I am wondering if Gromacs can do the following work?
Assuming I have a pdb file of an RNA molecule. Some atoms may be
too close or even overlap, I am wondering if Gromacs can move
the atoms to reasonable positions and remove the bad contacts?
The final structure is supposed to be the "most" stable
structure with minimal energy. I know AMBER minimization can do
this work, and I am wondering if Gromacs can do the same?
Gromacs is perfectly capable of energy minimization. Whether or not
minimization succeeds and gives a reasonable output is mostly
dependent upon the feasibility of the starting structure. If the
minimization crashes, it's not Gromacs' fault ;)
-Justin
--
==============================__==========
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu <http://vt.edu> | (540) 231-9080
<tel:%28540%29%20231-9080>
http://www.bevanlab.biochem.__vt.edu/Pages/Personal/justin
<http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin>
==============================__==========
--
gmx-users mailing list gmx-users@gromacs.org
<mailto:gmx-users@gromacs.org>
http://lists.gromacs.org/__mailman/listinfo/gmx-users
<http://lists.gromacs.org/mailman/listinfo/gmx-users>
Please search the archive at
http://www.gromacs.org/__Support/Mailing_Lists/Search
<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
<mailto:gmx-users-requ...@gromacs.org>.
Can't post? Read http://www.gromacs.org/__Support/Mailing_Lists
<http://www.gromacs.org/Support/Mailing_Lists>
--
Best,
Liang Liu
--
========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
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 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