Hi Carlos,
Thank you so much for the lightning-speed response again! I'll try the loadpdb command and continue with the simulation and update here if everything works out well!
Yen Nguyen
Graduate Research Assistant
George Mason University
Department of Chemistry and Biochemistry
________________________________
From: Carlos Simmerling <carlos.simmerling.gmail.com>
Sent: Thursday, November 18, 2021 3:33 PM
To: AMBER Mailing List <amber.ambermd.org>
Subject: Re: [AMBER] How to add an amidated residue at the C-terminal in AMBER tleap ?
when building a system in leap, you can either load the pdb, or specify the
sequence if you don't have a structure (technically you can do both, but
that's advanced).
if you have a pdb file, try replacing the "peptide = sequence" command with
"peptide = loadpdb xx.pdb" where xx.pdb is your pdb file name.
On Thu, Nov 18, 2021 at 3:27 PM Yen Nguyen <ynguyen3.gmu.edu> wrote:
> Hi Carlos,
>
> Thank you so much for the quick response. I just want to clarify if what I
> understand is right.
>
> So to build my modified peptide sequence in tleap, I'd type the following
> command (AA stands for amino acid)
> peptide = sequence {ACE TRP AA2 AA3 ... LYS NHE } ?
>
> The ACE TRP would give me acetylated Trp on the N-terminal and LYS NHE
> would give me amidated Lys on the C-terminal? Is that correct?
>
> I forgot to mention that I'm doing a peptide simulation using the
> guidelines from this tutorial on AMBER
> https://secure-web.cisco.com/1iV0T-fHKHEmJqsLjwxLmkiQDG9bolcox6e78L2Pw29UsoCaEbjT6Qx1hrl0rY7vZxmyOx4qbfnJLrai_el2SaAjmb0OXgmVC1a6UOmhAAqwb1RUerOSmZCGOgX-ES4NfWzqmXhe3qHLR6n6EQnLtK-jRu_uPpYY-GlshJaEVe6SU6cwx8Dnn5AqF7JGWh2sZOhZOsWjL-iXYiO6X0YuhL72kr4ywxG9vHUke2ZFjHoSYgiD1qvRY1Ai2-iKwDCqvHgUVYAm1WP2ZjNDd3yF01ZibBMylLar6qVEz8eCCDTAZ-BasBM3lG_Kj7A9dHt2LwDi27JYYphuXKI3K9A_05A1TQlSZlDDnsO0kL0gXCqyGmkA2Waeo4u6YAxfW6cT94D3pPRFuyfEt-_K3Id0uQhgCgorPbBp3OUfBf_sSnhoY3RLb-KgrtXr5iqvmCh-h/https%3A%2F%2Fambermd.org%2Ftutorials%2Fbasic%2Ftutorial0%2Findex.php%23Solvate_alanine_dipeptide
> Simple Simulation of Alanine Dipeptide - AMBER<
> https://secure-web.cisco.com/1iV0T-fHKHEmJqsLjwxLmkiQDG9bolcox6e78L2Pw29UsoCaEbjT6Qx1hrl0rY7vZxmyOx4qbfnJLrai_el2SaAjmb0OXgmVC1a6UOmhAAqwb1RUerOSmZCGOgX-ES4NfWzqmXhe3qHLR6n6EQnLtK-jRu_uPpYY-GlshJaEVe6SU6cwx8Dnn5AqF7JGWh2sZOhZOsWjL-iXYiO6X0YuhL72kr4ywxG9vHUke2ZFjHoSYgiD1qvRY1Ai2-iKwDCqvHgUVYAm1WP2ZjNDd3yF01ZibBMylLar6qVEz8eCCDTAZ-BasBM3lG_Kj7A9dHt2LwDi27JYYphuXKI3K9A_05A1TQlSZlDDnsO0kL0gXCqyGmkA2Waeo4u6YAxfW6cT94D3pPRFuyfEt-_K3Id0uQhgCgorPbBp3OUfBf_sSnhoY3RLb-KgrtXr5iqvmCh-h/https%3A%2F%2Fambermd.org%2Ftutorials%2Fbasic%2Ftutorial0%2Findex.php%23Solvate_alanine_dipeptide
> >
> In order to run a molecular dynamics simulation, we need to load a force
> field to describe the potential energy of alanine dipeptide. We will use
> the AMBER force field FF19SB1 for proteins. To learn more about force
> fields, please look through section 3.1.1 on page 36 of the Amber 2020
> Manual.. 5.
> ambermd.org
> I don't really see where my pdb file could be used in the entire process.
> But I'm wondering if I could use the command peptide = loadpdb "my peptide
> pdb file" to use it for the next steps in the tutorial, e.g. solvation,
> creating prmtop and rst files, etc.
>
> Thank you.
>
> Yen Nguyen
> Graduate Research Assistant
> George Mason University
> Department of Chemistry and Biochemistry
> ________________________________
> From: Carlos Simmerling <carlos.simmerling.gmail.com>
> Sent: Thursday, November 18, 2021 2:55 PM
> To: AMBER Mailing List <amber.ambermd.org>
> Subject: Re: [AMBER] How to add an amidated residue at the C-terminal in
> AMBER tleap ?
>
> I think you're looking for the "NHE" residue.
> it should work ok if your pdb file has the matching residue and atom names.
> if not, you can rename them following the info below:
>
> NHE INT 1
> NH2 ENDING GROUP
> 4 N N M 3 2 1 1.335 116.600 180.000 -0.46300
> 5 HN1 H E 4 3 2 1.010 119.800 0.000 0.23150
> 6 HN2 H E 4 3 2 1.010 119.800 180.000 0.23150
>
> On Thu, Nov 18, 2021 at 2:48 PM Yen Nguyen <ynguyen3.gmu.edu> wrote:
>
> > Hello,
> > I've been searching for a way to build a peptide sequence with C and N
> > terminal modifications (C -terminal = amidated Lys, N-terminal =
> acetylated
> > Trp in tleap (AMBER), but I couldn't find a way to make/add an amidated
> Lys
> > residue at the C-terminal. AMBER has the code ACE for acetylated
> residues,
> > but I don't know if it has a code for amide terminals?
> > I'm using amber 18 and the protein forcefield leaprc.protein.ff14SBonlysc
> > , which has the following parameter files
> > /amber20/dat/leap/parm/parm10.dat
> > frcmod.ff14SB
> > frcmod.ff99SB14
> > I have a pdb file of the peptide, but it looks like AMBER doesn't utilize
> > pdb files to create prmtop and coordinate files for my peptide.
> > I'd greatly appreciate any help. Thank you!
> >
> > Yen Nguyen
> > Graduate Research Assistant
> > George Mason University
> > Department of Chemistry and Biochemistry
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> >
> http://secure-web.cisco.com/1fh-2UEqh3O0g_ySSr6pZRevbVJm3HxWahVx5vO6YZAbnAYUfRhYs-cTKTG_WOAt7olzl8XL4tirSCJzCX4vtZEvwYbQkkwM-Hk3QidJW9QAsmOS7L8Rv2hScOt1wxlWOIBaXjDiqweXwXNdpFHVSlHK_omImM0Fq6gldzr8PnVc4gwDrr7Hbnukl0sw-cLbVVrpofA7ab3NdbXIe9VhKfdw9YAQ9owiQDCkpsCKN3PH3-law-db_8xvTsTu1qS-mBgk10iDP4lrQxmpv-1dQk2Fr2rfgDv6QaXDaPaSYITMbRdOcNluTMt5-ly7HCGIP9ryDeKK6GedDqAUsUdD_eiKEi4GMn9pzJodzP67Ifq8DZaFKA526vDf5xpx3lZ_q40OH27jmQJWPLVMs_yOrEy77dnypropwjJTH1TB9mSo4GBAlPJgIEIQKu0xfVGu8/http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber
> >
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
>
> http://secure-web.cisco.com/1fh-2UEqh3O0g_ySSr6pZRevbVJm3HxWahVx5vO6YZAbnAYUfRhYs-cTKTG_WOAt7olzl8XL4tirSCJzCX4vtZEvwYbQkkwM-Hk3QidJW9QAsmOS7L8Rv2hScOt1wxlWOIBaXjDiqweXwXNdpFHVSlHK_omImM0Fq6gldzr8PnVc4gwDrr7Hbnukl0sw-cLbVVrpofA7ab3NdbXIe9VhKfdw9YAQ9owiQDCkpsCKN3PH3-law-db_8xvTsTu1qS-mBgk10iDP4lrQxmpv-1dQk2Fr2rfgDv6QaXDaPaSYITMbRdOcNluTMt5-ly7HCGIP9ryDeKK6GedDqAUsUdD_eiKEi4GMn9pzJodzP67Ifq8DZaFKA526vDf5xpx3lZ_q40OH27jmQJWPLVMs_yOrEy77dnypropwjJTH1TB9mSo4GBAlPJgIEIQKu0xfVGu8/http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://secure-web.cisco.com/1rYfmOhb3gHoYbtvURGLGhSrVqhtiPRqF3LqCjVmk59LEBrJnhBhSQ4Pn9STaaNwoDgP4k1luLwxdNKmq9wxGwWsqT_soSA6vqB5x5QzDVVSEdOY4eL605-2v5j_YPDjMizktLx1lWDGsM-ViDWQ7kGkisoYR68lU5lCadJh0kfSaib9BZ2OVcTwn2PLQTAhHKVOw5ZTDgNo-DFbNutFwMgkuCnfIFqihinE8hfNIDNeWa-oNowbPfqp95T948fLPbHR8KPGVI0PUBMSnINSuXhKobyUFdjTaJf8T7sTpxT4l6MPOoRU9wvsyDI0aNfUh7JBYJKNUJHCjhsYhFchq65K8w6n0EtNy6Lgk0I19Qnx8rnH3XPc00JnzlJIkZTW9rvxIEl6VsrMbFfnqGD1KQYRWEsh8C7NutbVF2OeHMmDKkDbXu5BsKb4H6f2IDElX/http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://secure-web.cisco.com/1rYfmOhb3gHoYbtvURGLGhSrVqhtiPRqF3LqCjVmk59LEBrJnhBhSQ4Pn9STaaNwoDgP4k1luLwxdNKmq9wxGwWsqT_soSA6vqB5x5QzDVVSEdOY4eL605-2v5j_YPDjMizktLx1lWDGsM-ViDWQ7kGkisoYR68lU5lCadJh0kfSaib9BZ2OVcTwn2PLQTAhHKVOw5ZTDgNo-DFbNutFwMgkuCnfIFqihinE8hfNIDNeWa-oNowbPfqp95T948fLPbHR8KPGVI0PUBMSnINSuXhKobyUFdjTaJf8T7sTpxT4l6MPOoRU9wvsyDI0aNfUh7JBYJKNUJHCjhsYhFchq65K8w6n0EtNy6Lgk0I19Qnx8rnH3XPc00JnzlJIkZTW9rvxIEl6VsrMbFfnqGD1KQYRWEsh8C7NutbVF2OeHMmDKkDbXu5BsKb4H6f2IDElX/http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Nov 18 2021 - 13:00:03 PST
