Re: [AMBER] error with loading HEM parameters for constph redox calcuations

From: David Case <david.case.rutgers.edu>
Date: Sat, 1 Jun 2019 00:57:39 +0000

On Fri, May 31, 2019, Vaibhav Dixit wrote:

>I created a pdb with connectivity info from prmtop/inpcrd using ambpdb
>command and it shows expected bonds in pymol.

Don't rely (just) on pymol and ambpdb connecctivity. Learn to use
parmed's printBonds command to make sure what bonds are present (or
missing) from the prmtop file.


>It will be still helpful if someone could explain the meaning of the
>message tleap is giving about "*Residues lacking connect*".

An ordinary protein or DNA residue has "head" and "tail" connections,
telling which atoms are use to bond to the previous and following
residues. Residues (like WAT) that aren't part of a polymer are missing
those connections. An N-terminal amino acid will have only a tail
connection, since there is no previous residue. You CYP residue should
have both head and tail, since it is to be a part of a polypeptide
chain.

This information probably should not be in the tleap output, since it is
only *very* occasionally useful, and you can always just look at the
.lib or .off file (or use the desc command in tleap) to find out this
information. But it's been there for three decades, taking up space....

...dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri May 31 2019 - 18:00:02 PDT
Custom Search