Hi James,
Unfortunately leap cannot build atoms that it doesn't know about. Unlike the automatic capping N and C terminal residues for an entire protein, capping a custom residue will require you to add the ACE and NME groups manually or through another program.
For my ligand and residue building I use the freely available MarvinSketch followed by geometry optimization and RESP charge calculations using the R.E.D. server.
Best,
Parker
-----Original Message-----
From: James Starlight [mailto:jmsstarlight.gmail.com]
Sent: Tuesday, June 03, 2014 8:58 AM
To: AMBER Mailing List
Subject: Re: [AMBER] Parametrization of non-standart residues
Thanks!
Some methodological question about capping of my residue with ACE and NME groups, Does it possible to do it quickly using some tleap (not xleap) script (adding ACE to selected N and NME to C atoms of my non standart residue keeping initial geometry of it?
2014-05-31 16:53 GMT+02:00 Jason Swails <jason.swails.gmail.com>:
>
> On May 31, 2014, at 9:21 AM, James Starlight <jmsstarlight.gmail.com>
> wrote:
>
> > Many thanks for information!
> >
> > I'll try to spend time on RED tutorials before doing my parametrization.
> > However in cases when I need for blank topology of my residue (for
> > qm/calculations for instance) providing explicitly only its bonding
> > to
> the
> > rest of the protein as I understood antechamber parametrization
> > might be enough.
> > Assuming that I'm providing residue with CO in tail and NH in head
> > positions (all atoms corresponds to the complex.pdb) could I provide
> > mol3 file (with head and tail specification) to leap to build my complex?
>
> Indeed, if you plan on treating this residue entirely with QM in a
> QM/MM calculation, the charge derivation scheme does not matter at
> all. There will be no effective difference between antechamber and RED in this case.
> Only the van der Waals parameters are used for atoms in a quantum region.
>
> Check your PDB file -- does it have TER cards between your custom
> residue and the residues on either side? Certain things in the PDB,
> like a TER line for instance, will trigger leap to terminate one chain
> and start a new one, even if HEAD and TAIL atoms are set.
>
> See if you can figure out why tleap is terminating chains around your
> residue. The "desc" command will show you what a unit is (it
> DESCribes the unit). So you can look at the unit in tleap to make
> sure it has the properties you expect it to have.
>
> HTH,
> Jason
>
> --
> Jason M. Swails
> BioMaPS,
> Rutgers University
> Postdoctoral Researcher
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://scanmail.trustwave.com/?c=129&d=y8aN0zeQgumxw1ypTr23qzZaQAAlxiO
> uBNp_Oi7aFQ&u=http%3a%2f%2flists%2eambermd%2eorg%2fmailman%2flistinfo%
> 2famber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://scanmail.trustwave.com/?c=129&d=y8aN0zeQgumxw1ypTr23qzZaQAAlxiOuBNp_Oi7aFQ&u=http%3a%2f%2flists%2eambermd%2eorg%2fmailman%2flistinfo%2famber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Jun 03 2014 - 06:30:02 PDT
