Re: [AMBER] Cyclic peptide

From: Itamar Kass (Med) <"Itamar>
Date: Sun, 24 Mar 2013 21:47:17 +1100

Hi Jason,

Sorry for my late replay, it took me some time to implement you offer, but
it seems to work like charm.

Thanks again,
Itamar


On 12 March 2013 00:21, Jason Swails <jason.swails.gmail.com> wrote:

> On Mon, Mar 11, 2013 at 6:34 AM, Itamar Kass (Med)
> <itamar.kass.monash.edu>wrote:
>
> > Hi,
> >
> > In the last few years I have used GROMAC in my research. Recently I moved
> > to work with amber (amber12).
> >
> > Currently, I am trying to study cyclic peptides, in which the backbone
> > amino group of the first residue has a peptide-bond with the carboxyl
> > backbone group of the last residue. My question is how to do this in
> tleap.
> > I tried to use 'crosslink' and 'bond' but the FF missed parameters. for
> > this new bond, I guess because tleap translate the residues to N-ter and
> > C-ter. Is there anyway to overcome it without tinkering the topology
> file?
> >
>
> Huh, neat. My suggestion: copy the leaprc.ff12SB file
> ($AMBERHOME/dat/leap/cmd/) into your current directory as just "leaprc".
> Then, get rid of the residue mapping that defines terminal residues as the
> N- or C-terminal variants of those residues. It will look something like
> this:
>
> addPdbResMap {
> { 0 "HYP" "NHYP" } { 1 "HYP" "CHYP" }
> { 0 "ALA" "NALA" } { 1 "ALA" "CALA" }
> { 0 "ARG" "NARG" } { 1 "ARG" "CARG" }
> { 0 "ASN" "NASN" } { 1 "ASN" "CASN" }
> { 0 "ASP" "NASP" } { 1 "ASP" "CASP" }
> { 0 "CYS" "NCYS" } { 1 "CYS" "CCYS" }
> { 0 "CYX" "NCYX" } { 1 "CYX" "CCYX" }
> { 0 "GLN" "NGLN" } { 1 "GLN" "CGLN" }
> { 0 "GLU" "NGLU" } { 1 "GLU" "CGLU" }
> { 0 "GLY" "NGLY" } { 1 "GLY" "CGLY" }
> { 0 "HID" "NHID" } { 1 "HID" "CHID" }
> { 0 "HIE" "NHIE" } { 1 "HIE" "CHIE" }
> { 0 "HIP" "NHIP" } { 1 "HIP" "CHIP" }
> { 0 "ILE" "NILE" } { 1 "ILE" "CILE" }
> { 0 "LEU" "NLEU" } { 1 "LEU" "CLEU" }
> { 0 "LYS" "NLYS" } { 1 "LYS" "CLYS" }
> { 0 "MET" "NMET" } { 1 "MET" "CMET" }
> { 0 "PHE" "NPHE" } { 1 "PHE" "CPHE" }
> { 0 "PRO" "NPRO" } { 1 "PRO" "CPRO" }
> ...
>
> Then, manually draw the bond from the C-terminal carbon to the N-terminal
> nitrogen.
>
> While I think the above should work, I've never actually tried it.
>
> Good luck!
> Jason
>
> --
> Jason M. Swails
> Quantum Theory Project,
> University of Florida
> Ph.D. Candidate
> 352-392-4032
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>



-- 
===========================================
| Itamar Kass, Ph.D.
| Research Fellow
|
| Department of Biochemistry and Molecular Biology
| Building 77 Clayton Campus
| Wellington Road
| Monash University,
| Victoria 3800
| Australia
|
| Tel: +61 3 9902 9376
| Fax: +61 3 9902 9500
| E-mail: Itamar.Kass.monash.edu
============================================
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Mar 24 2013 - 04:00:03 PDT
Custom Search