Hello ambers,
I need to model a 4aa peptide crosslinked to DNA (N2-guanine) via
a propano group. I have generated residue definitions for the
various nonstandard components, and I can dock them (BDNA created
with nab). But I can't generate a topology file for the resulting
complex, seemingly because I have substantial atom overlap. Leap
complains of missing angle and dihedral terms, some of which don't
make sense (there are no such combinations in the pdb files).
That's why I think it must be because of the atom clashes. The
prep files for the nonstandard residue definitions pass parmchk.
The peptide does not "fit" in the minor groove of canonical BDNA,
so I don't see an obvious way to create a starting structure that
is free of steric overlap.
Does anyone have a good idea of a (topology file-independent) way
to create a starting structure? If I were clever enough, I might
be able to come up with some appropriate atom expressions to
create an bounds matrix to do some distance geometry runs in nab,
but what those expressions might be is not at all obvious. That's
the only thing I've been able to come up with so far.
Any ideas would be appreciated.
Bud Dodson
--
M. L. Dodson bdodson.scms.utmb.edu
409-772-2178 FAX: 409-747-8608
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Jul 02 2004 - 02:53:00 PDT
