<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body text="#000000" bgcolor="#ffffff">
    On 7/05/2011 7:16 PM, Mr Bernard Ramos wrote:
    <blockquote cite="mid:574624.42921.qm@web161520.mail.bf1.yahoo.com"
      type="cite">
      <table border="0" cellpadding="0" cellspacing="0">
        <tbody>
          <tr>
            <td style="font: inherit;" valign="top">
              <div>Thanks for the response. I will be doing a simple MD
                without&nbsp;the&nbsp;positional restraints. How do I go about
                reviewing the code for posres? What I did was to add new
                ff parameters for a new residue I introduced. May I
                please know what file in the program&nbsp;determines the
                posres so that I can revew them. Thanks.</div>
            </td>
          </tr>
        </tbody>
      </table>
    </blockquote>
    <br>
    Please give full descriptions the first time. I posted a noob
    question on the AMBER mailing list a few months back, and would
    never have gotten the problem solved except that I included some
    details that I thought were irrelevant, and someone spotted my
    problem. We've all got better things to do than play
    question-and-answer. You might have invented some atom names that
    pdb2gmx doesn't know how to cope with - but you haven't said whether
    your apparently erronously generated position restraints are
    specific to your new residue or not. Unfortunately you haven't
    shared your atom names after I suggested the atom names were
    probably the problem, and I'm not going to make guesses. There's
    various databases in $GMXLIB/share/top for various force fields, but
    I don't know your GROMACS version either...<br>
    <br>
    I don't think looking at the code will lead to a solution.<br>
    <br>
    Mark<br>
    <br>
    <blockquote cite="mid:574624.42921.qm@web161520.mail.bf1.yahoo.com"
      type="cite">
      <table border="0" cellpadding="0" cellspacing="0">
        <tbody>
          <tr>
            <td style="font: inherit;" valign="top">
              <div>--- On <b>Sat, 5/7/11, Mark Abraham <i><a class="moz-txt-link-rfc2396E" href="mailto:Mark.Abraham@anu.edu.au">&lt;Mark.Abraham@anu.edu.au&gt;</a></i></b>
                wrote:<br>
              </div>
              <blockquote style="padding-left: 5px; margin-left: 5px;
                border-left: 2px solid rgb(16, 16, 255);"><br>
                From: Mark Abraham <a class="moz-txt-link-rfc2396E" href="mailto:Mark.Abraham@anu.edu.au">&lt;Mark.Abraham@anu.edu.au&gt;</a><br>
                Subject: Re: [gmx-users] On posre.itp<br>
                To: "Discussion list for GROMACS users"
                <a class="moz-txt-link-rfc2396E" href="mailto:gmx-users@gromacs.org">&lt;gmx-users@gromacs.org&gt;</a><br>
                Date: Saturday, May 7, 2011, 5:10 PM<br>
                <br>
                <div id="yiv1037510764">On 7/05/2011 7:01 PM, Mr Bernard
                  Ramos wrote:
                  <blockquote type="cite">
                    <table border="0" cellpadding="0" cellspacing="0">
                      <tbody>
                        <tr>
                          <td valign="top">
                            <div>Hi everyone!</div>
                            <div>&nbsp;</div>
                            <div>I was able to generate a pdb file in
                              which the hydrogens for the molecule are
                              already included (that way I avoided using
                              the hdb file for the pdb2gmx). However,
                              when I checked the posre.itp, the some of
                              the restraints were applied to the
                              hydrogens.&nbsp;Is this correct because I am
                              not sure since&nbsp;I what I know is that&nbsp;the
                              restraints should be applied to heavy
                              atoms.</div>
                          </td>
                        </tr>
                      </tbody>
                    </table>
                  </blockquote>
                  <br>
                  Maybe your atom names are confusing pdb2gmx. Maybe
                  you're comparing apples and oranges with the indices
                  in the coordinate file and the indices in the
                  posres.itp (which are relative to the [atoms], IIRC).<br>
                  &nbsp;
                  <blockquote type="cite">
                    <table border="0" cellpadding="0" cellspacing="0">
                      <tbody>
                        <tr>
                          <td valign="top">
                            <div>Are the positional restraints needed
                              only for the&nbsp;equalibration step but not
                              for the production runs in normal MD
                              simulations. Thanks.&nbsp;</div>
                          </td>
                        </tr>
                      </tbody>
                    </table>
                  </blockquote>
                  <br>
                  Up to you. Do you want to sample under unphysical
                  restraints or not?<br>
                  <br>
                  Mark<br>
                </div>
                <br>
                -----Inline Attachment Follows-----<br>
                <br>
                <div class="plainMail">-- <br>
                  gmx-users mailing list&nbsp; &nbsp; <a moz-do-not-send="true"
href="http://us.mc1615.mail.yahoo.com/mc/compose?to=gmx-users@gromacs.org"
                    ymailto="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
                  <a moz-do-not-send="true"
                    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
                    moz-do-not-send="true"
                    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 moz-do-not-send="true"
href="http://us.mc1615.mail.yahoo.com/mc/compose?to=gmx-users-request@gromacs.org"
                    ymailto="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
                  Can't post? Read <a moz-do-not-send="true"
                    href="http://www.gromacs.org/Support/Mailing_Lists"
                    target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a></div>
              </blockquote>
            </td>
          </tr>
        </tbody>
      </table>
    </blockquote>
    <br>
  </body>
</html>