<br><br>----- Original Message -----<br>From: Mark Zottola &lt;mzottola@gmail.com&gt;<br>Date: Friday, June 4, 2010 10:40<br>Subject: Re: [gmx-users] Creation of a Non-Standard Residue<br>To: Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt;<br><br><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>Let me clarify.</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>I have no problem determining the bond stretching or angle bending parameters for HCN.&nbsp; This was a very straightforward task.&nbsp; Unit conversion to get those values in a gromacs force field was also a non-issue.</div>  <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>The problems I am having are:</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>1) How do I specify a POLAR hydrogen on a carbon atom?&nbsp; There are polar hydrogens (or so it seems) in the gromacs forcei field, but ostensibly they are on heteroatoms such as nitrogen or oxyen.&nbsp; Is there a way to force gromacs to let a polar hydrogen on a carbon?<br><br>ffgmx has been deprecated for years (as pdb2gmx tells users), so you'll find it hard to interest people with this problem. PRODRG beta works for one of the GROMOS forcefields, however. The general issue is that you're trying to develop an all-atom model, so you need to target an all-atom forcefield.<br></div>  <div>&nbsp;<br></div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>2) I AM looking for advice on the van der Waals parameters.&nbsp; Yes, I can guess and tweak until the cows come home, but obviously there are recommendations which I do not know.&nbsp; So... any suggestions on this.</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>3)&nbsp; The problem is quite simple, I want to develop parameters for HCN so I can use it as both ligand AND as a component n a solvent box.&nbsp; I do have the formula (4 or 5 steps from the web) to create my solvent box, but I need to have HCN properly parameterized.&nbsp; </div>  <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>4)&nbsp; I am asking for this as I cannot use the Dundee server to make the appropriate topology file.</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>The warning blurb was good at saying "Parameter development is done by professionals, so please do not try this at home."&nbsp; I am asking for help in the above issues as I am new to GROMACS but NOT&nbsp;other&nbsp;MD programs.</div><br>Well, then you'll know how to follow up on the concrete advice there - start by reading the literature for the force field of interest to see how VDW parameters were derived. How best to follow up is strongly problem- and forcefield-dependent. You might decide the path of least resistance is to use an AMBER forcefield in GROMACS so that you can use AMBERTools automated parameterization tools.<br><br>Mark