<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
Hi,<br><br>I don't understand exactly what is the requested feature here.<br>I am currently reorganizing the force field setup in Gromacs to more cleanly support AMBER and CHARMM<br>and adding some extra functionality. I started a discussion on the gmx-developers list on this topic.<br>On feature I have added is that rtp, hdb, etc files no longer have fixed names and you can have multiple<br>of them. So you can just put, e.g., a file called ligand.itp in your force field or current dir and pdb2gmx<br>will read it.<br><br>Berk<br><br><hr id="stopSpelling">From: alanwilter@gmail.com<br>Date: Fri, 29 Jan 2010 09:20:31 +0000<br>To: js@drugdiscoveryathome.com<br>CC: gmx-users@gromacs.org; av@drugdiscoveryathome.com<br>Subject: [gmx-users] Re: including a custom itp file in topology<br><br>Hi Jack,<div><br></div><div>Look, if I understand well, what you want is to be able to something like tleap does for amber, where you load all the libs you need, including the ones you created for your ligand(s), and then generate the MD input files at once needing just the complexed pdb or mol2 as input.</div>
<div><br></div><div>For pdb2gmx to do the same, you would need to tweak the files ..../gromacs/top/ffamber99sb.rtp, hdb etc., which I once considered doing that for acpypi but then it would add much more complexity for no much gain.</div>
<div><br></div><div>What I hope is for the new *pdb2gmx* (gmx 4.1 maybe? or only for gmx 5.0?) to feature, among other things, this capability as you requested.</div><div><br></div><div>Cheers,</div><div>Alan<br><br><div class="ecxgmail_quote">
On Thu, Jan 28, 2010 at 17:13, Jack Shultz <span dir="ltr"><<a href="mailto:js@drugdiscoveryathome.com">js@drugdiscoveryathome.com</a>></span> wrote:<br><blockquote class="ecxgmail_quote" style="border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div>Hi,</div>
<div> </div>
<div>I was trying to figure out if there is a short-cut for what I'm doing. I have complexes that I'm trying to prep using pdb2gmx. The ligand does not have a standard residue name. The way I know this can work is seperating out the ligand and protein into seperate files and preping the ligand using acpypi and the protein using pdb2gmx. Then incorporating them into a single pdb complex and including a reference to the ligand.itp (generated by acpypi) into a complex topology file. Is there any shortcut to doing this? any way to reference the ligand's itp file when running pdb2gmx?<br clear="all">
<br>-- <br>Jack<br><br><a href="http://drugdiscoveryathome.com">http://drugdiscoveryathome.com</a><br><a href="http://hydrogenathome.org">http://hydrogenathome.org</a><br></div>
</blockquote></div><br><br clear="all"><br>-- <br>Alan Wilter Sousa da Silva, D.Sc.<br>PDBe group, PiMS project <a href="http://www.pims-lims.org/">http://www.pims-lims.org/</a><br>EMBL - EBI, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SD, UK<br>
+44 (0)1223 492 583 (office)<br>
</div> <br /><hr />New Windows 7: Find the right PC for you. <a href='http://windows.microsoft.com/shop' target='_new'>Learn more.</a></body>
</html>