Re: [AMBER] cpptraj and nastruct

From: Daniel Roe via AMBER <amber.ambermd.org>
Date: Mon, 10 Jul 2023 10:33:01 -0400

Hi Jing,

Thanks for the files. I have 2 things for you to try. First, your DNA
contains a non-standard protonated cysteine (residue DCP), the name of
which 'nastruct` doesn't recognize by default. You can add 'resmap'
keywords to map non-standard residues back to recognized nucleic acid
templates, e.g.

resmap DCP:C resmap D3C:C

I need to add a warning to cpptraj to make it clear when 'resmap'
might be needed.

Also, the default cutoffs that 'nastruct' currently uses to determine
base pairing haven't been well-tested on non-WC base paired
structures. Try using values that are closer to what 3DNA uses:

origincut 15.0 zcut 1.5 zanglecut 30

So in summary, your 'nastruct' command should look something like:

nastruct NA naout dna resmap DCP:C resmap D3C:C \
  origincut 15.0 zcut 1.5 zanglecut 30

Try that and let me know how it looks.

-Dan

On Fri, Jul 7, 2023 at 2:00 PM Daniel Roe <daniel.r.roe.gmail.com> wrote:
>
> On Fri, Jul 7, 2023 at 11:49 AM David A Case via AMBER
> <amber.ambermd.org> wrote:
> >
> > I think Dan is correct about needing to see a structure. Parallel strand
> > duplexes may have never been tested, and it's not hard to imagine that they
> > could fool an algorithm trying to find "Watson-Crick" base pairs.
>
> The base pair determination algorithm was recently updated (version
> 6.18.0 https://github.com/Amber-MD/cpptraj/pull/1018), so in theory
> non-WC base pairing should also be detected (which can be the case in
> e.g. Z-DNA, G-quadruplexes etc). However, the new "better"
> determination method is only about 4 months old so there is probably
> still room for improvement.
>
> -Dan

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Jul 10 2023 - 08:00:02 PDT
Custom Search