On Sun, Oct 11, 2015, Robert Molt wrote:
>
> source leaprc.ff14SB
> loadamberprep unit_1.prepin
> loadamberparams unit_1.frcmod
> loadamberprep unit_2.prepin
> loadamberparams unit_2.frcmod
> foo=loadpdb file_name.pdb
>
> for unit_1 and 2 are the new nucleic acid residues, file_name.pdb is the
> input geometry to which I would like to assign parameters. When I save a
> pdb to look at it, ....
I'm already lost: what commands did you use to "save a pdb"?
> the above description of overlapping atoms applies.
Also lost: can you be more explicit about what the "above description" means?
What evidence are you using to say that it "applies"?
> Moreover, in leap, when I load the pdb, there are many missing heavy
> atoms in the leap record for each of the new units. Hence why I have
> overlapping atoms.
A big problem with Amber documentation is the paucity of examples, but that
problem also applies here. How about providing an example? In my experience,
if I get a "missing heavy atom" message, it is usually straightforward to
compare the input pdb file with the prepin (or other library) file, and see
what is different. So you are going to have to provide a simple example to
help us see why you can't pinpoint the differences between the prepin and the
pdb files.
> 4.) This would lead me to believe that I must have a mistake in my .mc
> file listed above, that the connections are being made erroneously. I
> assume to just post that, although I can imagine someone asking for any
> of the other many commands involved in this complicated process (but for
> the sake of brevity, I will leave it as this, since short emails=good).
As in all debugging, try to make the "complicated process" simpler: do you
see the problem with just a single prepin file? Can you make a simple test
case for which you provide all the necessary files to reproduce it?
....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Oct 12 2015 - 06:00:06 PDT
