2
$\begingroup$

I have a range of systems consisting of multiple peptide chains, some of which are one residue (i.e. single amino acids), which I'd like to analyse using molecular dynamics and AMBER-flavour forcefields. Each peptide fragment is a connected chain with unique chain IDs, separated by TER records.

Is there a scripted or automated tool out there to add neutral terminal caps - e.g. n-methyl (NME) and methyl-acetate (ACE)?

Some options that work well but don't seem to help with this specifically:

  • tleap: requires manual editing to change the residue name of one atom at the end of each chain into NME and ACE. After this, tleap can add missing atoms, but the key step is manual
  • pdbfixer: can add charged terminal caps, but not NME/ACE
  • pdb2gmx: does not add NME/ACE; requires manual editing of the PDB file
  • CHARMM-GUI: not scripted/automated
  • HTMD / Maestro: apparently these would work but I don't have a license!
$\endgroup$
2
  • $\begingroup$ I’m rusty with this but I had a script programmatically create input for tleap and I believe ACE was right before our first residue. I believe we would have to create a timer between inputs because tleap would start to have issues if input was made all at once. Would something like this be feasible? I’ll see what I can pull off my research computer. $\endgroup$ Commented Aug 12, 2022 at 2:08
  • 1
    $\begingroup$ In Ruby. No docs, sorry. Our docs were taken down years ago when our lab closed. Mine has ACE and NME. gist.github.com/MelanieS/74b578ac2d2a0dd50372 $\endgroup$ Commented Aug 12, 2022 at 2:14

1 Answer 1

2
$\begingroup$

you can try use pymol to do this.

use Ctrl+middle click to pick an atom,e.g.the N atom of first residue,

then >>Build>>adresidue>>Ace. same as NME.

and save the new structure.

$\endgroup$

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .