Re: [AMBER] writing dna-ligand topology using amber t-leap

From: David A Case via AMBER <amber.ambermd.org>
Date: Sat, 25 Jun 2022 08:04:04 -0400

On Sat, Jun 25, 2022, İmren Bayıl via AMBER wrote:

>my problem is that
>occurred complex_clean2_wat.prmtop complex_clean2_wat inpcrd these two
>files include B ligand I couldnt sure because I will use these two files in
>namd md simulation.

It's hard to provide much advice here: (a) we don't know what is in the
complex_clean2.pdb file; (b) we don't know if any errors were reported in
the tleap run.

You can use ParmEd (or a text editor) to see what is inside prmtop files.
Maybe that will answer your questions. You can also use savePdb inside
tleap, or ambpdb afterwards, to create a PDB file that you can examine in
a visualization program.

I recommend avoiding the use of the "combine" command in tleap. Just make
a pdb file that has exactly the components you want, then load that into
tleap.


Below, I've added some comments about your tleap script. Also note that it
is often helpful to do things in steps, instead of trying to create a really
complicated script. Stop after each "saveAmberParm" step, and examine the
results.

>tleap
>source leaprc.constph
>source leaprc.gaff
>gaff = loadamberparams gaff.dat # this is redundant
> lig=loadmol2 B.mol2
>loadamberparams B.frcmod
> saveoff lig B.lib # this and the next line don't do anything
>loadoff B.lib # and could end up overwriting "lig"
>saveamberparm lig B.prmtop B.inpcrd
>source leaprc.water.tip3p
>loadOff atomic_ions.lib # leaprc.water.tip3 already does this
>loadAmberParams frcmod.tip3p # "
>loadAmberParams frcmod.ions1lm_126_tip3p # "
>loadAmberParams frcmod.ionsjc_tip3p # "
> loadAmberParams frcmod.ions234lm_126_tip3p # "
>source leaprc.DNA.bsc1
>loadoff DNA.parmBSC1.lib # leaprc.DNA.bsc1 already does this
> loadoff parmBSC1.lib # "
> loadamberparams frcmod.parmbsc1 # "
> rec=loadpdb complex_clean2.pdb # check your naming conventions: you
                      # are calling this "rec" but the file name has
                      # "complex"
>com=combine {rec lig}
>check com
>saveamberparm com complex_clean2_gas.prmtop complex_clean2_gas.inpcrd
>savepdb com complex_clean2_gas.pdb
>solvatebox com TIP3PBOX 10
>charge com
>addions com Na+ 0
> saveamberparm com complex_clean2_wat.prmtop complex_clean2_wat.inpcrd

...good luck....dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Aug 04 2022 - 13:33:14 PDT
Custom Search